I’m trying to use SNAP with Python. I tried the snappy examples and everything worked fine. But I would want to perform radiometric calibration (i.e., calculating sigma0) and terrain correction with a Sentinel-1 product. Where can I find the functions for that? I already looked into http://step.esa.int/docs/v2.0/apidoc/engine/ and found the “Orthorectifier” function, if that’s the right one. But I can’t find anything about “Calibration”.
Is there some kind of user-guide or handbook for using the Python functions? I already found this one: https://senbox.atlassian.net/wiki/display/SNAP/How+to+use+the+SNAP+API+from+Python
but that’s not very supportive for further functions.
Maybe someone can help.
I was in the same situation as you. You can find snappy Operators when you typ “gpt -h” in your command window. To get more informations about possible input parameters you can search for “gpt -h < operators >” (e.g. “gpt -h Calibration”)
I already implemented some of them in a python script (bandMath, kMeans, PCA, Subset, Reprojection,…). If you can use them…
I would like to calibrate SAR (out of TC operation) by using the python.
I do appreciate if you could advise me how I can reach python or snappy scripts for “Calibration"(from SLC to sigma0) and "filtering”, if possible.
I did not create a script for Calibration or Filtering. But you assume it can create the script easily yourself:
Use the GLCM script as basic
Replace for example func_code=GLCM with func_code=Calibration (name of the Functions, when you call gpt -h)
Use the correct parameters as described in gpt -h Calibration.
I hope this helps. If it does not work, then there are maybe other functions behind it.
You can also call these functions via a command prompt. In this case you need to write the parameters into an XML-File (also described in gpt -h Calibration) and call it. But I did not use this approach until now. So I can not tell you how exactly it works
I have tried Calibration code and it can be opened correctly in terms of sigma0 value.
But the problem is that I can not use it as the input for the Terrain Correction due to the fact that Tie Point Grids (e.g., “Slant range time” and “Incidence Angle”) are missing. Here is the code, does anyone know by any chance what goes wrong with?
Your advice will be a great help.
Thanks
I was wondering if you finally solved your issue ? I also need to do a radiometric calibration, terrain correction and speckle filtering on Radarsat-2 product and I would like to write a snappy script but I can’t find examples of scripts for these processing.