Jak wypełnić luki w obrazach SLC-Off Landsat 7 ETM + za pomocą QGIS?

Wyjaśnienie problemu:

W dniu 31 maja 2003 r. Czujnik Landsat 7 Enhanced Thematic Mapper (ETM +) miał awarię korektora linii skanowania (SLC). Od tego czasu wszystkie obrazy Landsat ETM + miały kliny po obu stronach każdej sceny, co spowodowało utratę danych o około 22%. Te obrazy są dostępne do bezpłatnego pobrania ze strony USGS GloVis i znajdują się w kolekcji L7 SLC-off.

źródło: Uniwersytet Yale

Czy mogę użyć Qgis zamiast oprogramowania zastrzeżonego (ENVI, Erdas ....), aby wypełnić luki?

Możesz wypróbować narzędzie gdal_fillnodata , które jest również dostępne w QGIS poprzez menu Raster-> Analiza-> Wypełnij nodata. Wykorzystuje interpolację odwrotnego ważenia odległości (IDW). Właśnie wypróbowałem zarówno tę metodę, jak i interpolację triangulacyjną pojedynczej daty (w ENVI) i gdal_fillnodata wyglądały znacznie lepiej. Jeśli chcesz scalić wiele dat, być może będziesz musiał zaimplementować zlokalizowaną metodologię dopasowania liniowego histogramu w pythonie (GDAL + Numpy), aby wyglądała dobrze.

Jeden ze współautorów tego artykułu zapewnia implementację open source w IDL dla swojego algorytmu do wypełniania luk w danych SLC poza danymi, nazwany algorytmem NSPI (Neighborhood Similar Pixel Interpolator).