Monday 16 November 2015

remote sensing - How to fill gaps in SLC-Off Landsat 7 ETM+ Images with QGIS?


Explanation of the problem:



On 31 May 2003 the Landsat 7 Enhanced Thematic Mapper (ETM+) sensor had a failure of the Scan Line Corrector (SLC). Since that time all Landsat ETM+ images have had wedge-shaped gaps on both sides of each scene, resulting in approximately 22% data loss. These images are available for free download from the USGS GloVis website and are found in the L7 SLC-off collection.



source: Yale University


Can I use Qgis instead of proprietary software (ENVI, Erdas....) to fill the gaps?



Answer



You could try the gdal_fillnodata tool which is also available in QGIS via the Raster->Analysis->Fill nodata menu. It uses an inverse distance weighting (IDW) interpolation. I just tried both that method and the single date Triangulation interpolation (in ENVI) and gdal_fillnodata looked much better. If you want to merge multiple dates, you might have to implement the localised linear histogram matching methodology in python (GDAL+Numpy) to get it to look good.



No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...