ISCE (stripmapstack) to StaMPS (PS) - Error in llh2local (line 39). z=llh(2,:)~=0;

342 views
Skip to first unread message

Tonian Robinson

unread,
Jun 11, 2021, 10:48:30 AM6/11/21
to MAINSAR
Hello,

I am using ISCE ( stackStripMap ) to create StaMPS - PS Processing with TerraSAR-X datasets. With the following steps:

1.       In ISCE prep_TSX.py then stackStripMap.py processing with workflow  –W slc:

               i.          !stackStripMap.py -s /home/toni/WCF_TerraSARX_processing/scenes/HS_021/SLC/ -d /home/toni/WCF_TerraSARX_processing/scenes/HS_021/DEM/demLat*.dem.wgs84 -t 44 -a 2 -r 2 -u snaphu -W slc --nofocus --zero

               ii.         Ran all files run_files 1-7 using: !run_files/run_0*

2.      In terminal StaMPS processing dir:

               i.          $ make_single_reference_stack_isce  input_file

Input_file:
source_data slc_stack
slc_stack_path /home/toni/WCF_TerraSARX_processing/scenes/HS_021/merged/SLC/
slc_stack_reference 20161118
slc_stack_geom_path /home/toni/WCF_TerraSARX_processing/scenes/HS_021/merged/geom_reference/
slc_stack_baseline_path /home/toni/WCF_TerraSARX_processing/scenes/HS_021/merged/baselines/
maskfile /home/toni//WCF_TerraSARX_processing/scenes/HS_021/merged/geom_reference/shadowMask.rdr
range_looks 2
azimuth_looks 2
aspect_ratio 4
lambda 0.056
slc_suffix
geom_suffix

 

              ii.          $ cd INSAR_20161118/

              iii.          $ mt_prep_isce 0.4

              iv.          $ matlab

3.      In Matlab command window: >> stamps

              i.          After running stamps it goes to the PATCH_1 dir /home/toni/WCF_TerraSARX_processing/scenes/HS_021/StaMPS/INSAR_20161118/PATCH_1

              ii.          Then error message:

STAMPS: Will process patch subdirectories

 

STAMPS: ########################################

STAMPS: ################ Step 1 ################

STAMPS: ########################################

STAMPS: Directory is PATCH_1

 

Found inc angle file, will load the data

Loading data into Matlab...

SETPARM: heading = -169.7074

SETPARM: lambda = 0.056

Index in position 1 exceeds array bounds.

 

Error in llh2local (line 39)

   z=llh(2,:)~=0;

Error in ps_load_initial_isce (line 171)

xy=llh2local(lonlat',ll0)*1000;

 

Error in stamps (line 281)

                    ps_load_initial_isce(data_inc) 


Has anyone ever received this error message? Any idea how to address it?
I have attached the directory layout.
ISCEtoSTAMPS_error_in_llh2local_dir_structure.png

Suri Babu

unread,
Jun 12, 2021, 12:05:08 AM6/12/21
to mai...@googlegroups.com
try this in matlab

setparm("merge_resample_size", 0)

and run the stamps 

--
You received this message because you are subscribed to the Google Groups "MAINSAR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mainsar+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/mainsar/7312e6d8-cf18-4ebd-9467-0dc736f7907cn%40googlegroups.com.

Tonian Robinson

unread,
Jun 12, 2021, 3:46:39 PM6/12/21
to MAINSAR
Thank you but I am still getting the same error message. I also noticed that my pscands.1.ll file in PATCH_1/ is empty which I assume is the root cause of this error message. Any idea why?

Tonian Robinson

unread,
Aug 30, 2021, 2:09:54 PM8/30/21
to MAINSAR
Still haven't solved this.
I've realized that the selpsc_patch.c script is producing empty pscands.1.ij and pscands.1.da files which in return makes the pscands.1.hgtpscands.1.ll, and pscands.1.ph files empty. Any help will be appreciated!

Tonian Robinson

unread,
Oct 12, 2021, 2:00:04 PM10/12/21
to MAINSAR
I noticed that outside of the empty pscands.* files I was also seeing a Segmentation Fault (core dumped) error message while running mt_prep_isce
I referenced this snap forum to address this.

Solved by 
  1. Installing g++-7 and gcc-7 (just incase, to make sure both are present)
  2. Changing the compiler from CC = g++ to CC = g++-7 in the StaMPS_master/src/Makefile 
  3. cd StaMPS_master/src
  4. make
  5. make install

Suri Babu

unread,
Oct 12, 2021, 9:52:50 PM10/12/21
to mai...@googlegroups.com
Thank you

--
You received this message because you are subscribed to the Google Groups "MAINSAR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mainsar+u...@googlegroups.com.
Reply all
Reply to author
Forward
0 new messages