Extracting the central strip from LANDSAT 7 imagery

Extracting the central strip from LANDSAT 7 imagery

February 8, 2016

Here is a simple Python code to extract the central strip from Landsat 7 imagery (SLC-off),  that is not affected by the SLC failure. The algorithm shrinks the striping zones through a morphological filter (erosion) and creates a new shapefile AOI that extracts the desired raster extent without striping effects. The code is based on Python for ArcGIS (arcpy) – so you require a ArcGIS license.

General steps:

  1. Loop through all Landsat 7 data folders
  2. Stack bands for each image
  3. Create a mask
  4. Erode the mask by 20 pixels
  5. Convert the mask to polygon
  6. Create a minimum bounding box
  7. Clip the original raster through the bbox

 

import arcpy
from arcpy.sa import *

import sys,os

#  Environment settings (Activate Spatial Analyst, Overwrite Outputs allowed and TIFF compression is LZW)
arcpy.CheckOutExtension("spatial")
arcpy.env.overwriteOutput = True
arcpy.env.compression = 'LZW'

# this is your main directory with all unzipped Landsat datasets
 rootdir = "D:\\DATA\\Landsat7\\"

# create scratch folder "temp" 
temp = "D:\\DATA\\temp\\"

# loop through directory with all unzipped Landsat 7 folders
 for files in os.listdir(rootdir):   
    files = os.path.join(rootdir, files)   
    
    # for each loop the subdir "files" is now the current workspace 
    # (e.g. LE71520322015157-SC20160224113319) that contains the Landsat bands
    arcpy.env.workspace = files  
    rasters = arcpy.ListRasters("*", "TIF")  
    
    # create empty array
    stack_liste = []  
    # loop through all rasters in subdir
    for raster in rasters:   

        image = arcpy.Raster(raster) 
        name  = image.name 
        index = name.split("_")[0]  

        # fill up the array only with the actual spectral bands        
        sr = "_sr_band"  
        if sr in raster:   
            stack_liste.append(raster)             

    # now stack all bands within the array
    stack_name = files + "\\" + index + "_stack.tif"    
    arcpy.CompositeBands_management(stack_liste, stack_name)  

    # convert the image stack to a mask by logical operation with an absurd value that will result in an output "0"
    con = EqualTo(stack_name, 123456789)  

    # now shrink the raster mask with value "0" by 20 pixels
    shrink = temp + "shrink"  
    shrinking = Shrink(con, 20, 0) 
    shrinking.save(shrink)  

    zone = temp + "zone.shp" 
    bbox = temp + "bbox.shp"  

    # conver the shrunk mask to polygon and create a minimum bounding box
    arcpy.RasterToPolygon_conversion(shrink, zone, "NO_SIMPLIFY", "VALUE") 
    arcpy.MinimumBoundingGeometry_management(zone, bbox, "RECTANGLE_BY_WIDTH", "NONE")  

    # now use that bounding box as a mask to cut out the central nadir strip from the original stack
    # Final result 
    extract = files + "\\" + index + "_aoi.tif"  
    ExtractByMask = arcpy.sa.ExtractByMask(stack_name, bbox) 
    ExtractByMask.save(extract)

 

you may also like:

Bridging Scales: How Radar Satellites supports Crop Monitoring

Bridging Scales: How Radar Satellites supports Crop Monitoring

In an era of climate uncertainty and increasing pressure on agricultural systems, understanding how crops grow and respond to environmental stress is more important than ever. A new study led by researchers from Martin-Luther-University Halle-Wittenberg, in close...

Upcoming PhD Defense by Ariane Droin

Upcoming PhD Defense by Ariane Droin

Ariane Droin will defend her PhD thesis "Permeabilität und Erreichbarkeit lokaler Nachbarschaften im urbanen Kontext. Eine geographische Analyse auf Basis räumlicher Netzwerke." on September 16th at 4 p.m. at the John-Skilton Straße 4a, Seminar Room 2/00.B.03.  ...

PhD Defense by Dorothee Stiller

PhD Defense by Dorothee Stiller

Dorothee Stiller will defend her PhD thesis "Potential of Remote Sensing Data and Methods for Urban Transport Research" on 15th of September at 4 p.m. at the John-Skilton Straße 4a, seminar Room 2/00.B.03. Everyone who is interested is cordially invited to join her...

Strengthening Collaboration with SANParks for Conservation Research

Strengthening Collaboration with SANParks for Conservation Research

Our long-standing collaboration with Dr. Corli Coetsee and Dr. Ben Wigley from SANParks is moving forward with promising new research activities. The joint work is focusing on mapping savanna features more accurately such as trees, paths, or animals through innovative...

Exploring Future Collaborations on Fire Research in African Savannas

Exploring Future Collaborations on Fire Research in African Savannas

During recent discussions, new opportunities for collaboration emerged between Navashni Govender, Senior Conservation Manager at SANParks in Kruger National Park, Prof. Katharina Breininger, head of the Pattern Recognition Lab in Informatics, and Dr. Mirjana Bevanda...

Understanding Urban Heat in Germany

Understanding Urban Heat in Germany

In a world where summers grow ever hotter, understanding and combating urban heat islands is becoming more urgent than ever. A recent study by our Prof. Hannes Taubenboeck sheds new light on this challenge—and at its helm is Dr. Tobias Leichtle, Dr. Thilo Erbertseder...