Problem pikseli NoData w rastrach

Problem

Mamy kilka plików ASCII grid NMT (w formacie ARC/INFO ASCII GRID). Wszystkie z tym samym 32-bitowym typem piksela z wartością zmiennoprzecinkową i rozdzielczością 5 m. Chcemy połączyć te pliki w jeden plik o rozdzielczości 5 m. Tworzymy plik mozaiki. Niestety na krawędziach niektórych elementów mozaiki znajdują się luki z pikselami o wartościach NoData.

Raster NoData Gaps
Fig. 1. Połączone rastry z widocznymi lukami z pikselami o wartościach NoData

Rozwiązanie

Rozwiązanie polega na zamianie wartości pikseli NoData wartościami średniej arytmetycznej obliczonymi na podstawie najbliższego sąsiedztwa pikseli luk. W tym przypadku wartości pikseli zostaną obliczone na podstawie średniej wyciągnietej z matrycy 5 × 5 pikseli. W tym celu wybierzmy:

ArcToolbox > Spatial Analyst Tools > Map Algebra > Raster Calculator

W kalkulatorze rastrów wprowadź wyrażenie:

Con(IsNull("raster"), FocalStatistics("raster", NbrRectangle(5,5, "CELL"), "MEAN"), "raster")

Uwaga!: słowo raster zamień nazwą zbioru rastrowego z tabeli zawartości

Źródła

How To: Remove and replace no data values within a raster using statistical information from the surrounding data values.