crop the area of interest

Added by francesco menniti about 1 year ago

Hi, I have a problem! I'm trying to get time series from ISCE to STAMPS!
The workflow works fine without error but does not cut the area of interest, although I set the bbox. So Isce creates interferograms as big as the whole scene and StaMPS has too many PS to process fast!
could someone kindly help me?

stackSentinel.py -s /home/protepos/Scrivania/immagini_sar_carsule -d /home/protepos/Carsule/dem_carsule/demLat_N41_N44_Lon_E011_E014.dem -b '42.61 42.66 12.52 12.58' -a /home/protepos/sentinel_aux_cal/ -o /home/protepos/Scrivania/orbits/ -W slc

I also wanted to ask if it is possible to cut the area of interest just before entering the command mt_prep_isce 0.6 4 3 50 200


Replies (28)

RE: crop the area of interest - Added by francesco menniti about 1 year ago

Sorry again, can anyone help?
Why this Warning? Maybe the problem is right here..

----2019-07-15 19: 04: 22.207 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_LINES has been set to the default value 1525 which is the number of lines in the slc image.---

Running: resamp_withCarrier
['--slave', '/home/protepos/Scrivania/albero_processi/ASC_Carsule/slaves/20190619', '--master', '/home/protepos/Scrivania/albero_processi/ASC_Carsule/master', '--coregdir', '/home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619', '--azimuth_misreg', '/home/protepos/Scrivania/albero_processi/ASC_Carsule/misreg/azimuth/dates/20190619.txt', '--range_misreg', '/home/protepos/Scrivania/albero_processi/ASC_Carsule/misreg/range/dates/20190619.txt']
Estimated burst offset: 0
Shifts: {0: 0, 1: 0}
Polynomial Order: 0 - by - 0
0
Polynomial Order: 0 - by - 0
0
Chi squared: 0.000096
Polynomial Order: 5 - by - 3
12103.1 -710.779 59.2121 1.50482
-47650.1 2642.01 -171.684 -5.1366
46899.7 -2446.52 108.939 4.02977
2.00459e-08 -2.34221e-08 8.52071e-09 0
-1.45194e-08 6.60111e-09 0 0
4.10614e-09 0 0 0
Misfit radians - Max: 0.00034796811723936116 , Min : -0.0002451525560900336
True
2019-07-15 19:04:22,207 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_LINES has been set to the default value 1525 which is the number of lines in the slc image.
2019-07-15 19:04:22,216 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_WIDTH has been set to the default value 25710 which is the width of the slc image.
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/slaves/20190619/IW3/burst_01.slc.vrt
API open (WR): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_01.slc
No Range Carrier provided.
Assuming zero range carrier.
Polynomial Order: 0 - by - 0
0
Polynomial Order: 0 - by - 0
0
Polynomial Order: 5 - by - 3
12103.1 -710.779 59.2121 1.50482
-47650.1 2642.01 -171.684 -5.1366
46899.7 -2446.52 108.939 4.02977
2.00459e-08 -2.34221e-08 8.52071e-09 0
-1.45194e-08 6.60111e-09 0 0
4.10614e-09 0 0 0
Polynomial Order: 0 - by - 0
0
Polynomial Order: 0 - by - 0
0
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/range_01.off.vrt
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/azimuth_01.off.vrt
Polynomial Order: 0 - by - 2
0.487717 -833.941 354483

<< Resample one image to another image coordinates >>
Input Image Dimensions: 
25710 pixels
1525 lines
Output Image Dimensions: 
25710 pixels
1525 lines
Number of threads:            8
Complex data interpolation
Initializing Sinc interpolator
Azimuth Carrier Poly
Polynomial Order: 5 - by - 3
12103.1 -710.779 59.2121 1.50482
-47650.1 2642.01 -171.684 -5.1366
46899.7 -2446.52 108.939 4.02977
2.00459e-08 -2.34221e-08 8.52071e-09 0
-1.45194e-08 6.60111e-09 0 0
4.10614e-09 0 0 0
Range Carrier Poly
Polynomial Order: 0 - by - 0
0
Range offsets poly
Polynomial Order: 0 - by - 0
0
Azimuth offsets poly
Polynomial Order: 0 - by - 0
0
Reading in the image
At line 1000
Interpolating image
At line 1000
Elapsed time: 7.89062500
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/range_01.off.vrt
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/azimuth_01.off.vrt
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/slaves/20190619/IW3/burst_01.slc.vrt
API close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_01.slc
API open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_01.slc
API close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_01.slc
Adjust valid samples
Before: 564 23638
Offsets : 9.73828411102 11.0706443787
After: 557 23630
Adjust valid lines
Before: 29 1468
Offsets : 5.69881296158 5.72307872772
After: 28 1461
Chi squared: 0.000776
Polynomial Order: 5 - by - 3
12103.6 -697.923 109.886 2.35371
-47652 2616.74 -271.837 -6.56573
46901.6 -2446.44 108.941 4.02967
1.79518e-08 -2.13603e-08 7.95563e-09 0
-1.29839e-08 5.93388e-09 0 0
3.68163e-09 0 0 0
Misfit radians - Max: 0.002587864568340592 , Min : -0.0013358887154026888
True
2019-07-15 19:04:42,257 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_LINES has been set to the default value 1525 which is the number of lines in the slc image.
2019-07-15 19:04:42,257 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_WIDTH has been set to the default value 25710 which is the width of the slc image.
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/slaves/20190619/IW3/burst_02.slc.vrt
API open (WR): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_02.slc
No Range Carrier provided.
Assuming zero range carrier.
Polynomial Order: 0 - by - 0
0
Polynomial Order: 0 - by - 0
0
Polynomial Order: 5 - by - 3
12103.6 -697.923 109.886 2.35371
-47652 2616.74 -271.837 -6.56573
46901.6 -2446.44 108.941 4.02967
1.79518e-08 -2.13603e-08 7.95563e-09 0
-1.29839e-08 5.93388e-09 0 0
3.68163e-09 0 0 0
Polynomial Order: 0 - by - 0
0
Polynomial Order: 0 - by - 0
0
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/range_02.off.vrt
GDAL open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/azimuth_02.off.vrt
Polynomial Order: 0 - by - 2
0.616637 -1504.33 905632
<< Resample one image to another image coordinates >>
Input Image Dimensions: 
25710 pixels
1525 lines
Output Image Dimensions: 
25710 pixels
1525 lines
Number of threads:            8
Complex data interpolation
Initializing Sinc interpolator
Azimuth Carrier Poly
Polynomial Order: 5 - by - 3
12103.6 -697.923 109.886 2.35371
-47652 2616.74 -271.837 -6.56573
46901.6 -2446.44 108.941 4.02967
1.79518e-08 -2.13603e-08 7.95563e-09 0
-1.29839e-08 5.93388e-09 0 0
3.68163e-09 0 0 0
Range Carrier Poly
Polynomial Order: 0 - by - 0
0
Range offsets poly
Polynomial Order: 0 - by - 0
0
Azimuth offsets poly
Polynomial Order: 0 - by - 0
0
Reading in the image
At line 1000
Interpolating image
At line 1000
Elapsed time: 7.98437500
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/range_02.off.vrt
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/azimuth_02.off.vrt
GDAL close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/slaves/20190619/IW3/burst_02.slc.vrt
API close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_02.slc
API open (R): /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_02.slc
API close: /home/protepos/Scrivania/albero_processi/ASC_Carsule/coreg_slaves/20190619/IW3/burst_02.slc
Adjust valid samples
Before: 564 23638
Offsets : 9.72492218018 11.0448322296
After: 557 23630
Adjust valid lines
Before: 28 1470
Offsets : 5.69615364075 5.72022771835
After: 28 1461
creating /home/protepos/Scrivania/albero_processi/ASC_Carsule/stack

RE: crop the area of interest - Added by Jerry Kuang 12 months ago

Hi francesco,

I have the same problem about cropping the AOI even though I already set the bounding box in the running of stackSentinel.py.

May I ask, have you found the solution for this problem? How can I crop the area with the bounding box?

Thanks a lot.

Best regards,

Jerry

RE: crop the area of interest - Added by francesco menniti 12 months ago

Hi Jerry
Unfortunately I could not find a solution to this problem!
In September I will try to solve the problem! if I succeed, I will write the solution here!

Best regards

RE: crop the area of interest - Added by Jerry Kuang 12 months ago

Hi francesco,

Thanks for your reply.

I found an option to limit the swath number in the stackSentinel.py running, but not sure if it can make the bounding box option work in the running. Will give it a try.

When I looked into the code, it should read the bbox as input. However, it would process much larger areas instead of the bbox areas.

Thanks again.

Best regards,

Jerry

RE: crop the area of interest - Added by francesco menniti 11 months ago

Hi Jerry
May I ask, have you found the solution for this problem?
Best regards

Francesco

RE: crop the area of interest - Added by Jerry Kuang 11 months ago

Hi Francesco,

Seems like even the updated github code doesn't fix this problem. However, I found an indirect solution to cut the data.

After running the mt_prep_isce with slected number of patches in range and azimuth, you can modify the path_list file to choose the pathes which jointly cover the aoi you want.

That solution could significantly reduce the time of processing and coverage.

Best regards,

Jerry

RE: crop the area of interest - Added by francesco menniti 11 months ago

Hi Jerry
Thanks for your reply.
I have also tried this solution, but I fear there might be some problems in unwrapping step.
what do you think about it?
Best regards,

Francesco

RE: crop the area of interest - Added by francesco menniti 11 months ago

HI!
the problem comes from this warning! But I still haven't figured out how to solve it!
Has anyone else encountered this problem?

2019-09-16 12:37:35,784 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_LINES has been set to the default value 1510 which is the number of lines in the slc image.
2019-09-16 12:37:35,794 - isce.stdproc.resamp_slc - WARNING - The variable INPUT_WIDTH has been set to the default value 25957 which is the width of the slc image.

RE: crop the area of interest - Added by Collin Cong 4 months ago

Hi Francesco
Have you found the solution for this problem? I have the same problem.

Collin

RE: crop the area of interest - Added by francesco menniti 4 months ago

Hi Collin

Unfortunately I was unable to understand where the problem came from and how to fix it.
Can anyone help us out?

Francesco

RE: crop the area of interest - Added by Collin Cong 4 months ago

Hi Francesco

You can try to use the -x parameter to bound the ROI.
eg.

stackStripMap.py -s ../SLC/ -x ‘33.05 34.10 108.68 109.74’ -d ../DEM/demLat_N33_N35_Lon_E108_E110.dem.wgs84 -m 20171127 -t 360 -b 500 -u snaphu --nofocus

Best regard with you.

Collin

RE: crop the area of interest - Added by francesco menniti 3 months ago

Hi Collin,
thanks for the reply.

I have to use the stackSentinel.py for elaborate stack of coregistred interferogram for use STAMPS.

I am using the

stackSentinel.py s ..........  -d..........  -o..............  -b-------- -W slc

I cant use the -x in to stack sentinel.py

Can you help me?

Best regard with you.

Francesco

RE: crop the area of interest - Added by francesco menniti 3 months ago

Hi,
I think these warnings derive from errors in the installation of the software and in the setting of the PATHs.

RE: crop the area of interest - Added by Collin Cong 3 months ago

Hi Francesco,
When you process the Sentinel1 data, you'd better make use of the topsStack.py script. You must onfigure environment variables before using the new script.
The topsStack parameter '-b' can help to crop the AOI.

Best regard with you.
Collin

RE: crop the area of interest - Added by francesco menniti 3 months ago

Hi Collin,
thanks for the answer.
I used anaconda for setting ISCE.
however I have used stackSentinel.py

stackSentinel.py -s /root/............ -d /root/demLat_N42_N46_Lon_E008_E014.dem.wgs84 -a /......./sentinel_aux_cal -o /orbits -b '44.1735 44.1941 9.9159 9.9591' -C NESD -W slc

I also entered the -b parameter to perform the processing.
I also discussed it in another discussion trying to interpret the warnings to try to trace the problem of crop area of interest, but but I still can't understand where I'm wrong and why ISCE gives me these warnings.

do you think i might have set the environment variables wrong?

Best regard with you.
Francesco

RE: crop the area of interest - Added by Collin Cong 3 months ago

Sorry Francesco, I haven't used the stackSentinel.py script. It may be an old version.
I think you should ask some others who familiar with the script for help.

Best regard to you.
Collin

RE: crop the area of interest - Added by Heresh Fattahi 3 months ago

Hi Francesco,

Just to be sure, your "-b" option should be "SOUTH NORTH WEST EAST". Is that set properly?

In your working directory, how does your configs/config_master file looks like?

Heresh

RE: crop the area of interest - Added by francesco menniti 3 months ago

hi Heresh
thanks for the quick response.

yes, I used the -b '44.1735 44.1941 9.9159 9.9591'

South=ymin= wgs84= 44.1735
North=ymax=wgs84= 44.1941
West=xmin=wgs84=9.9159
East=xmax=wgs84=9.9591

Attached there is the config_master and the picture of area of interest.

When I enter BBOX, should IISCE cut out my area of interest? right?

Best regard to you.

Francesco

TRYINSAR1.jpg (75 kB)

config_master_forum (737 Bytes)

RE: crop the area of interest - Added by Eric Fielding 3 months ago

Hi Francesco,

I think you may need to process a larger area with stackSentinel.py. Your area is only a few kilometers across.

++Eric

RE: crop the area of interest - Added by Heresh Fattahi 3 months ago

The smallest processing unit is a busrt in ISCE. We don't allow to subset within a burst. Based on your boudning box, the busrsts that cover it get extracted. In your case, your bbox is too samll. Choose a larger one. If your problem is with your time-series analysis, then you should be able to crop the coregistered SLCs yourself after processing with ISCE is done.

Heresh

RE: crop the area of interest - Added by francesco menniti 3 months ago

thanks to both of you. Now I get it! to cut the area you don't have to use the stackSentinel.py bbox.
StackSentinel is only used to extract the burst, not to crop the area of interest.
P.s.I tried to use a bigger bbox but the warning is released anyway

isce.stdproc.resamp_slc - WARNING - The variable INPUT_LINES has been set to the default value 1510 which is the number of lines in the slc image.

thank you so much for your quick reply

RE: crop the area of interest - Added by Eric Fielding 3 months ago

The Sentinel-1 TOPS burst is about 20 km along-track, so you need to use a bbox bigger than 20 km in the north-south direction.

RE: crop the area of interest - Added by Heresh Fattahi 3 months ago

The warning is fine. That's the total number of lines in a burst.

RE: crop the area of interest - Added by francesco menniti 2 months ago

thank you very much again! now I understand better how the software works

RE: crop the area of interest - Added by jun sun about 1 month ago

Hello,everyone, I encountered a similar problem.

I use stackSentinel.py processing to get some results (see attachment)

My crop area: the north-south span is 0.40°, and the east-west span is 0.45°

What caused this? Is it just because the crop area is too small?

1 2 Next » (1-25/28)