Preprocessing of Sentinel-1 SAR data via Snappy Python module

Preprocessing of Sentinel-1 SAR data via Snappy Python module

m

August 1, 2016

This chapter demonstrates the Snappy Python module for the automatization of the ESA SNAP tool.

Code examples will be shown for an automated processing chain for the preprocessing of Sentinel-1 SAR data including Calibration, Subsetting and Terrain Correction of GRD (Ground Range Detected data).

A detailed installation tutorial for snappy can be found here: https://senbox.atlassian.net/wiki/display/SNAP/How+to+use+the+SNAP+API+from+Python

First, import the needed Python modules:

 

import snappy

from snappy import ProductIO
from snappy import HashMap

import os, gc   
from snappy import GPF

GPF.getDefaultInstance().getOperatorSpiRegistry().loadOperatorSpis()
HashMap = snappy.jpy.get_type('java.util.HashMap')

Now loop through all Sentinel-1 data sub folders that are located within a super folder (of course, make sure, that the data is already unzipped):

path = "D:\\SENTINEL\\"
 for folder in os.listdir(path):

   gc.enable()
   
   output = path + folder + "\\"  
   timestamp = folder.split("_")[4] 
   date = timestamp[:8]

Then, read in the Sentinel-1 data product:

   sentinel_1 = ProductIO.readProduct(output + "\\manifest.safe")    
   print sentinel_1

If polarization bands are available, spolit up your code to process VH and VV intensity data separately. The first step is the calibration procedure by transforming the DN values to Sigma Naught respectively. You can specify the parameters to output the Image in Decibels as well.

   pols = ['VH','VV'] 
   for p in pols:  
      polarization = p    
    
      ### CALIBRATION
  
      parameters = HashMap() 
      parameters.put('outputSigmaBand', True) 
      parameters.put('sourceBands', 'Intensity_' + polarization) 
      parameters.put('selectedPolarisations', polarization) 
      parameters.put('outputImageScaleInDb', False)  

      calib = output + date + "_calibrate_" + polarization 
      target_0 = GPF.createProduct("Calibration", parameters, sentinel_1) 
      ProductIO.writeProduct(target_0, calib, 'BEAM-DIMAP')

Next, specify a subset AOI to reduce the data amount and processing time. The AOI specified by its outer polygon corners and is formatted through a Well Known Text (WKT).

      ### SUBSET

      calibration = ProductIO.readProduct(calib + ".dim")    
      WKTReader = snappy.jpy.get_type('com.vividsolutions.jts.io.WKTReader')

      wkt = "POLYGON((12.76221 53.70951, 12.72085 54.07433, 13.58674 54.07981, 
                      13.59605 53.70875, 12.76221 53.70951))"

      geom = WKTReader().read(wkt)

      parameters = HashMap()
      parameters.put('geoRegion', geom)
      parameters.put('outputImageScaleInDb', False)

      subset = output + date + "_subset_" + polarization
      target_1 = GPF.createProduct("Subset", parameters, calibration)
      ProductIO.writeProduct(target_1, subset, 'BEAM-DIMAP')

Apply a Range Doppler Terrain Correction to correct for layover and foreshortening effects, by using the SRTM 3 arcsecond product (90m) that is downloaded automatically. You could also specify an own DEM product with a higher spatial resolution from a local path:

      ### TERRAIN CORRECTION
 
      parameters = HashMap()     
      parameters.put('demResamplingMethod', 'NEAREST_NEIGHBOUR') 
      parameters.put('imgResamplingMethod', 'NEAREST_NEIGHBOUR') 
      parameters.put('demName', 'SRTM 3Sec') 
      parameters.put('pixelSpacingInMeter', 10.0) 
      parameters.put('sourceBands', 'Sigma0_' + polarization)
 
      terrain = output + date + "_corrected_" + polarization 
      target_2 = GPF.createProduct("Terrain-Correction", parameters, subset) 
      ProductIO.writeProduct(target_2, terrain, 'GeoTIFF')

Fergana_Sentinel

you may also like:

New EAGLEs take off into the Winter Term 2025/26

New EAGLEs take off into the Winter Term 2025/26

As in previous years, the next generation of EAGLE Master's students from around the world gathered at the Earth Observation Research Center (EORC) on the first day of the winter term to begin their studies at the University of Würzburg. Prof. Dr. Tobias Ullmann...

Recording the Sounds of a River

Recording the Sounds of a River

Over the weekend, EORC PI Florian Betz met with Martina Cecchetto and Riccardo Fumigalli from the University of Padua to conduct ambient sound recordings and collect photographs of the Lech River, one of the major tributaries of the upper Danube. The photographs and...

Our PhD Wall is Growing — and So Is Our Research Family!

Our PhD Wall is Growing — and So Is Our Research Family!

It’s been a remarkable year for our research team! The PhD Wall of Fame, showcasing all past and current doctoral researchers, has officially reached its limits — and we’ve had to expand it to make room for even more success stories. So far six PhD defenses have taken...

🎉 A Sweet Surprise for a Special Birthday!

🎉 A Sweet Surprise for a Special Birthday!

At our department, we not only work hard together — we also celebrate the milestones that make our team so special. This week, we had the joy of surprising our wonderful secretary Tine Linge on her 60th birthday! Early in the morning, colleagues gathered to prepare a...

Contribution at SilviLaser Conference in Quebec

Contribution at SilviLaser Conference in Quebec

At SilviLaser 2025 in Québec City, PhD candidate Julia Rieder (EORC, University of Würzburg and staff member of EO4CAM) presented her work on "European Beech under Drought: Effects of Topography, Competition and Soil Water Availability." Her study uses LiDAR to reveal...

EORC at Remote Sensing Symposium in Darmstadt

EORC at Remote Sensing Symposium in Darmstadt

On 2 October 2025, Dr. John Friesen and Dr. Julian Fäth from the Earth Observation Research Cluster (EORC) at the University of Würzburg and staff members of EO4CAM took part in the symposium "Vom Orbit zur Entscheidung: Satellitenfernerkundung in der...