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.
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’
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
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?
#####################################################################
#####################################################################
***** 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.
I did run it and I found no issues on that AOI. Which issue was given in step7?
Hello
in the elaboration I noticed results that are different from those expected in the various articles. Especially in the areas circled in green in the last image. it seems to eliminate some of the deformation.
I have attached the plots.
Do you have any suggestions?
many Thanks
Francesco Menniti
Hi
I run the python sbas_topsar_ifg.py project.txt command but it seems to have not created the IFGs.
I continued with the python sbas_export.py project.txt command and the command seems to work fine.
When I run the command mt_prep_snap … it gives me this mean amplitude 0.
pening /home/protepos/ELABORAZIONE_BAIA_SNAP_SBAS/ASC/INSAR_20180508/SMALL_BASELINES/20180108_20200813/rslc/20180108.rslc…
WARNING : SLC /home/protepos/ELABORAZIONE_BAIA_SNAP_SBAS/ASC/INSAR_20180508/SMALL_BASELINES/20180108_20200813/rslc/20180108.rslchas ZERO mean amplitude
Mean amplitude = 0
And when I run stamps(1,5) i get this error
Index in position 2 exceeds array bounds.
Error in sb_load_initial_gamma (line 97)
rg=rgn+ij(:,3)*rps;
Error in stamps (line 222)
sb_load_initial_gamma;
can you help me @dongyusen @ABraun
thank you very much for your availability
