Please test the SBAS_snap2stamps

Hi, all
I write some code to generate the SBAS to StaMPS, it seems working now.

Please test it, thx

I had finished the data processing part in SNAP, but now has trouble in export the SBAS dataset to StaMPS.
Does anybody know the data structure of SBAS ?

Thx

it seems working now:

add 1: (1) import with mt_prep_snap; (2) modify the parms to small_baseline_flag to y, (3)modify line 141 in sb_load_initial_gamma.m, (4) run stamps in matlab

3 Likes

Such exciting news for the last long period, Congratulations to you for your breaking through the barrier between snap and SBAS!I’ll have a try.

@dongyusen ,

Thanks for sharing.
I have been looking at the graphs and it seems to me that coregistration is done in a single pair (master-slave) basis.

But I guess you know that for SBAS it is needed that the images had been coregistered first with a supermarket.

Can you please bring some light here and explain it to us?
Best

Fist, all the slaves are coregistrated to the super-master, and all the others. Then followed this post to re-organize the pairs:

I am not sure it’s right or not, because I have no isce or other softwares to double check finial results.

If u have any results, pls let me know or post it.

The new sbas_export.py had uploaded, it is more easy to run stamps now. please update the sbas_exprot.py from:

Hi dongyusen,

Thank you for sharing your script and informing this great news! Could you please show some initial results of your processing so that we have a look?

The data structure of SBAS for StaMPS is very similar to PS, except that you should create a new folder named ‘SMALL_BASELINES’ to store all the interferograms.

Regarding the double check, I have an idea that you can check whether the phase closure is consistent without any multilooking or spatial filtering.

Many thanks,

Fei.

sbas_snap2stamps 1.0, now tested!!

Hi and thank you very much for sharing sbas_snap2stamps.
I tried to process using snap2stamps by mdelgado (thanks for sharing) an area of Campi Flegrei, but on STAMPS I was unable to obtain the expected result due to problems in step7.
So I would like to test the SBAS snap technique.
But when I use sbas_topsar_coreg.py I get this error

Error: [NodeId: Enhanced-Spectral-Diversity] Operator ‘SpectralDiversityOp’: Unknown element ‘useSuppliedRangeShift’

  • org.jblas INFO Deleting /tmp/jblas7422041840348640256/libjblas_arch_flavor.so
  • org.jblas INFO Deleting /tmp/jblas7422041840348640256/libjblas.so
  • org.jblas INFO Deleting / tmp / jblas7422041840348640256
    033 [1; 35m Error in processing … coreg_ifg2run

Could you give me a hand to correct this Error? (i am using SNAP 6.0)baia.txt (1.1 KB)

Many thanks,
Francesco Menniti

have you seen this solution?

ok I solved it by updating because the sbas_topsar_coreg.xml file was already updated! So it was enough to update SNAP.

But now I have an error. I try to solve it … otherwise I will write on the forum

Executing processing graph
INFO: org.hsqldb.persist.Logger: dataFileCache open start
…INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
INFO: org.esa.snap.engine_utilities.download.DownloadableContentImpl: http retrieving http://step.esa.int/auxdata/dem/egm96/ww15mgh_b.zip
WARNING: org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp: RangeShiftOp: Cross-correlation failed for all bursts, set range shift to 0
0
0
90% done.

Many thanks,
Francesco Menniti

has the process finished? In most cases it is just a warning.

You can search for error messages in the search bar of the forum - many points have already been discussed and solved.

no he didn’t finish the process

Processing … 1 of 227 dataset…
Processing …20180508_20180108…

*****IW1 of 1IWs :20180508_20180108 coregistration

INFO: org.esa.snap.core.gpf.operators.tooladapter.ToolAdapterIO: Initializing external tool adapters
SEVERE: org.esa.s2tbx.dataio.gdal.activator.GDALDistributionInstaller: The environment variable LD_LIBRARY_PATH is not set. It must contain the current folder ‘.’.
INFO: org.esa.snap.core.util.EngineVersionCheckActivator: Please check regularly for new updates for the best SNAP experience.
Executing processing graph
INFO: org.hsqldb.persist.Logger: dataFileCache open start
…INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
INFO: org.esa.snap.engine_utilities.download.DownloadableContentImpl: http retrieving http://step.esa.int/auxdata/dem/egm96/ww15mgh_b.zip
WARNING: org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp: RangeShiftOp: Cross-correlation failed for all bursts, set range shift to 0
0
0
90% done.
org.esa.snap.core.gpf.OperatorException: 0
at org.esa.snap.core.gpf.graph.GraphProcessor$GPFImagingListener.errorOccurred(GraphProcessor.java:363)
at com.sun.media.jai.util.SunTileScheduler.sendExceptionToListener(SunTileScheduler.java:1646)
at com.sun.media.jai.util.SunTileScheduler.scheduleTile(SunTileScheduler.java:921)
at javax.media.jai.OpImage.getTile(OpImage.java:1129)
at com.sun.media.jai.util.RequestJob.compute(SunTileScheduler.java:247)
at com.sun.media.jai.util.WorkerThread.run(SunTileScheduler.java:468)
Caused by: org.esa.snap.core.gpf.OperatorException: 0
at org.esa.snap.engine_utilities.gpf.OperatorUtils.catchOperatorException(OperatorUtils.java:415)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.computeTileStack(SpectralDiversityOp.java:483)
at org.esa.snap.core.gpf.internal.OperatorImageTileStack.computeRect(OperatorImageTileStack.java:116)
at org.esa.snap.core.gpf.internal.OperatorImageTileStack.computeTile(OperatorImageTileStack.java:85)
at com.sun.media.jai.util.SunTileScheduler.scheduleTile(SunTileScheduler.java:904)
… 3 more
0
0
Caused by: org.esa.snap.core.gpf.OperatorException: 0
at org.esa.snap.engine_utilities.gpf.OperatorUtils.catchOperatorException(OperatorUtils.java:415)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.estimateAzimuthOffset(SpectralDiversityOp.java:783)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.computeTileStack(SpectralDiversityOp.java:457)
… 6 more
Caused by: java.lang.ArrayIndexOutOfBoundsException: 0
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.outputESDEstimationToFile(SpectralDiversityOp.java:1328)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.estimateAzimuthOffset(SpectralDiversityOp.java:778)
… 7 more

Error: 0
SEVERE: org.esa.snap.core.util.SystemUtils$SnapImagingListener: JAI error occurred: ‘Problem occurs when computing a tile by the owner.’ at com.sun.media.jai.util.SunTileScheduler@1e80cafd
org.esa.snap.core.gpf.OperatorException: 0
at org.esa.snap.engine_utilities.gpf.OperatorUtils.catchOperatorException(OperatorUtils.java:415)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.computeTileStack(SpectralDiversityOp.java:483)
at org.esa.snap.core.gpf.internal.OperatorImageTileStack.computeRect(OperatorImageTileStack.java:116)
at org.esa.snap.core.gpf.internal.OperatorImageTileStack.computeTile(OperatorImageTileStack.java:85)
at com.sun.media.jai.util.SunTileScheduler.scheduleTile(SunTileScheduler.java:904)
at javax.media.jai.OpImage.getTile(OpImage.java:1129)
at com.sun.media.jai.util.RequestJob.compute(SunTileScheduler.java:247)
at com.sun.media.jai.util.WorkerThread.run(SunTileScheduler.java:468)
Caused by: org.esa.snap.core.gpf.OperatorException: 0
at org.esa.snap.engine_utilities.gpf.OperatorUtils.catchOperatorException(OperatorUtils.java:415)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.estimateAzimuthOffset(SpectralDiversityOp.java:783)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.computeTileStack(SpectralDiversityOp.java:457)
… 6 more
Caused by: java.lang.ArrayIndexOutOfBoundsException: 0
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.outputESDEstimationToFile(SpectralDiversityOp.java:1328)
at org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp.estimateAzimuthOffset(SpectralDiversityOp.java:778)
… 7 more

– org.jblas INFO Deleting /tmp/jblas2107755608166537103/libjblas_arch_flavor.so
– org.jblas INFO Deleting /tmp/jblas2107755608166537103/libjblas.so
– org.jblas INFO Deleting /tmp/jblas2107755608166537103
033[1;35m Error in processing … coreg_ifg2run

hmmm, " Cross-correlation failed for all bursts" sounds bad…

Did you check if the split products have the same extent?

No I didn’t check. Sorry but I’m a beginner, but how can I check if the split products have the same extent?

open them in SNAP and look if they contain data and if all cover the same burst

1 Like

they all cover the same burst.

Do I have to have at least 2 burst cover? ((Resolved: ok with two Burst it works))

In black my area

the WARNING: org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp: RangeShiftOp: Cross-correlation failed for all bursts, set range shift to 0
is it a error or is it continuing to work?

#####################################################################

TOPSAR SBAS coregistration @@

#####################################################################

***** Totally 227 pairs !!! ****

Processing … 1 of 227 dataset…
Processing …20180508_20180108…

*****IW1 of 1IWs :20180508_20180108 coregistration

INFO: org.esa.snap.core.gpf.operators.tooladapter.ToolAdapterIO: Initializing external tool adapters
SEVERE: org.esa.s2tbx.dataio.gdal.activator.GDALDistributionInstaller: The environment variable LD_LIBRARY_PATH is not set. It must contain the current folder ‘.’.
INFO: org.esa.snap.core.util.EngineVersionCheckActivator: Please check regularly for new updates for the best SNAP experience.
Executing processing graph
INFO: org.hsqldb.persist.Logger: dataFileCache open start
.INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E013.SRTMGL1.hgt.zip
INFO: org.esa.snap.core.dataop.dem.ElevationFile: http retrieving http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
WARNING: org.esa.snap.core.dataop.dem.ElevationFile: http error:http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip on http://step.esa.int/auxdata/dem/SRTMGL1/N39E014.SRTMGL1.hgt.zip
WARNING: org.esa.s1tbx.sentinel1.gpf.SpectralDiversityOp: RangeShiftOp: Cross-correlation failed for all bursts, set range shift to 0
…10%…20%…30%…40%…50%…60%…70%…80%…90% done.
– org.jblas INFO Deleting /tmp/jblas3698547542201402362/libjblas_arch_flavor.so
– org.jblas INFO Deleting /tmp/jblas3698547542201402362/libjblas.so
– org.jblas INFO Deleting /tmp/jblas3698547542201402362
Tue Oct 6 17:07:41 2020

Finished coregistration processing of 20180508_20180108 in 2.12205718358 mins.

sounds good!

Now please check some of the interferograms for correct coregistration. You see an example here: Sentinel-1 TOPS interferometry (Figures 10 and 11)

I think the coregistration was not successful.

what could have gone wrong?

another question
is this the correct procedure? @dongyusen
1)splitting_master_free
2)splitting_slave_free
3)sbas_topsar_coreg
4)sbas_topsar_ifg
5)sbas_export.py

Many thanks,
Francesco Menniti

actually, this looks good to me. Water is never homogenous, so the important think is that the white areas have sharp boundaries. Both images are well aligned.

1 Like