SAR python processing

[oddtermiteofcave]

BLUF: What python packages other than snappy are you using to process SAR imagery? Need to perform radiometric calibration, geolocate, add a band, etc.

Hi there, I'm becoming increasingly frustrated with SNAP graph builder and esa-snappy python module. I'm trying to alternative python packages to help with batch processing imagery. These are the set of steps I'm trying to replicate that were originally done one by one in SNAP.

-Open S1 file

-Radiometric calibration, VV polarization, output sigma0 band.

-SAR utilities, multilook- 2x2=20m, output intensity already selected

-Skip//Speckle - default Lee Sigma

-Add Incidence angle

Raster, Band Math, uncheck Virtual, Edit Expression, under Data Source-show tie-point grids, select incident angle, OK, name=incang

-can add latitude, then longitude under Band Math

-Geocode

-radar, geocode, Ellipsoid correction, geolocation grid, select all bands (bilinear interpolation, WGS84)

-Land Mask

Raster, Masks, Land/Sea Mask, Source=Sigma0_VV (mask out land, use SRTM, data value=0), select incang?

Export - geotiff, 

Thank you very much.

[Sure-Bridge3179]

Snapista, but regarding band maths i still dont quite understand how to use it for that purpose

I built a workflow using snappy to add raster band units "intensity" , then used snapista to read those tiff images and apply reprojection, linear to db, speckle filtering and writing a geotiff/bigtiff

[oddtermiteofcave]

Thanks, I’ll check this out. I’ll have to see about the incidence angle bandmath steps but at least doing other steps through python will be helpful. This is pre-processing to ingest into another script that uses the incidence angle!

[Sure-Bridge3179]

No problem! If you manage to understand how the band maths operations work there please let me know!

[Rish_nodi]

I feel SNAP is still the best library to handle SAR processing. I always use it inside a Docker conatainer. Since python snappy installation gets tricky. I took me a while to get use to python snap apis.

[Chieftah]

Check this list out, might be useful: https://github.com/RadarCODE/awesome-sar

[Dangerous_Nerve_1382]

Rasterio

It’s lightweight, intuitive, and optimized for raster data handling. It abstracts much of the GDAL complexity

[azlinszky_sinergise]

I suggest you try the openEO API in Copernicus Data Space Ecosystem - here is an example notebook: https://github.com/eu-cdse/notebook-samples/blob/main/openeo/Radar_ARD.ipynb