How to call InSARStackOverview in python module

[sky_moon]

How to call InSARStackOverview in python module. where can I find a documentation about snappy

[marpet]

You can find information in the sticky post of this category

About the python category

Put your question regarding the usage of snappy, the snap-python interface or anything else related to Python and SNAP into this category. This could be questions regarding the access of the SNAP API from Python or implementing an operator in Python, all questions like this are best suited here. Question regarding Sen2Cor should go into the sen2cor category. General documentation about snappy can be found in How to use the SNAP API from Python of the Developer Guide. There is also a page which…

and also here

Snappy: Where to start?

Good afternoon, I am a new ESA Trainee and I would like to implement some of the processes of S2-tbx/sen2cor into a online platform. As a programming language I am going to use python. Are there beginner guides, handbooks or well documented scripts available? Until now I found: How to use the SNAP API from Python: https://senbox.atlassian.net/wiki/display/SNAP/How+to+use+the+SNAP+API+from+Python -> I could install snappy properly and run “snappy first contact” Scripts https://github.c…

You can also get to the Developer Guide from the website in the Documentation section.

[sky_moon]

Actually, I only want to use the function of calculate Baseline (it’s InSAR Stack Overview in desktop，I just want use it with python ). I want to find best Interference baseline in many files.
could you tell me how to use the function？ I didn‘t find in java scripts.

[marpet]

you can load the class as follows

import snappy
InSARStackOverview = snappy.jpy.get_type(org.esa.s1tbx.insar.gpf.InSARStackOverview)
stack_overview = InSARStackOverview()
stack_overview.setInput(…)

But the class InSARStackOverview does not have a calculateBaseline method. So I don’t know what you want to call.
InSAR Stack Overview is doing much more than just calling one method.

[sky_moon]

However , I can get perpendicular baseline by InSAR Stack Overview on the desktop. SO, the class InSARStackOverview will be not?

[marpet]

It is not the only one you need to use. But I can’t tell you want you need to do. I’m not working with SAR data.

[sky_moon]

I just want to realize the InSAR Stack Overview function of SNAP DESKTOP by snappy.

[marpet]

Maybe @cwong or @timmoorhouse can shed some light o it.
Or maybe someone else did it already?
But you can have a look at our repository and see what needs to be done.
The entry point is the class org.esa.s1tbx.insar.rcp.dialogs.InSARStackOverviewDialog and its method processStack.

[sky_moon]

Thank you very much. I will try to do it.Thanks

[ge46bex]

Does anyone found out how to call InSARStackOverview from snappy python?

I get the following error while calling the function, see the snippet below:

InSARStackOverview = snappy.jpy.get_type(‘org.esa.s1tbx.insar.gpf.In
SARStackOverview’)
stack_overview = InSARStackOverview()
stack_overview.calculateInSAROverview(insar_stack)

[marpet]

How have you created your insar_stack?
I think it needs to be created like this:

insar_stack=jpy.array('org.esa.snap.core.datamodel.Product', 5)
insar_stack[0] = product1
...
insar_stack[4] = product5
# or loop over it

[ge46bex]

@marpet Thanks for the reply, below is the full code.

slc_paths =  glob.glob(r'C:\Users\Placeholder\Desktop\Pre-Processing_Module\Quering\Download\*.zip')

insar_stack = snappy.jpy.array('org.esa.snap.core.datamodel.Product', len(slc_paths))

for i, slc_path in enumerate(slc_paths):

    slc_product = snappy.ProductIO.readProduct(slc_path)

    insar_stack[i] = slc_product

InSARStackOverview = snappy.jpy.get_type('org.esa.s1tbx.insar.gpf.InSARStackOverview')

stack_overview = InSARStackOverview()

stack_overview.calculateInSAROverview(insar_stack)

I still get the same error. Please let me know if I am doing something wrong.

Thanks in advance.

[marpet]

Now I see. calculateInSAROverview is a static method.
Don’t use the instance stack_overview of InSARStackOverview.

Try directly:
InSARStackOverview.calculateInSAROverview(insar_stack)

[ge46bex]

@marpet Thanks for your quick reply; I was successful in getting the name of the optimal master product. However, I was curious how I could extract the other data like perpendicular baseline, temporal baseline and modelled coherence.

snippet to get the name of optimal_master

InSARStackOverview = snappy.jpy.get_type(‘org.esa.s1tbx.insar.gpf.InSARStackOverview’)

optimal_master = InSARStackOverview.findOptimalMasterProduct(insar_stack).getName()

Thanks in advance

[marpet]

I’m not familiar wir SAR processing.
So, I’m also only guessing by looking at the code.

I think you need to call. This gives you an IfgStack array.
ifgStack = InSARStackOverview.calculateInSAROverview(insar_stack)

Each of the entries has a master and a master_slave. You can get the master_slave by
ifgpairs = ifgStack[0].getMasterSlave()
This is an array of IfgPair instances.
On these IfgPair instances you can call

ifgpairs[0].getPerpendicularBaseline()
ifgpairs[0].getTemporalBaseline()
ifgpairs[0].getCoherence()

and more.
See

senbox-org/s1tbx/blob/53450f70d346bb7cbf4e6ecd80c42d8ddaff1f8f/s1tbx-op-insar/src/main/java/org/esa/s1tbx/insar/gpf/InSARStackOverview.java#L292

[ge46bex]

@marpet Thank you so much, will try them.

[Nureldin]

Hi,
I’m using snap10 and python3.10 in linux, having now this error:
ValueError: Java class ‘eu.esa.microwave.insar.gpf.InSARStackOverview’ not found.
Any ideas.

[Marco_EOM]

@Nureldin

eu.esa.microwave.insar.gpf.InSARStackOverview

Try instead: eu.esa.sar.insar.gpf.InSARStackOverview

The page where this should be mentioned, does not reflect the latest decision for renaming the package.
API and Architectural Changes in SNAP 10 - SNAP - Confluence

@jun_lu @lveci

[Nureldin]

Thank you for the quick reply. It works now.

Regards