gmtsar / gmtsar

GMTSAR

INSTRUCTIONS FOR INSTALLING AND RUNNING GMTSAR

INSTALL

1. Go to the GMTSAR WIKI and follow the instructions to install dependencies and download the GMTSAR package. https://github.com/gmtsar/gmtsar/wiki/GMTSAR-Wiki-Page

2. Download orbit files for just ERS and ENVISAT. The other satellites have orbits distributed with the data. https://topex.ucsd.edu/gmtsar/downloads/

3. Go to the gmtsar directory and enter:

   autoconf
   

4. Configure GMTSAR: Run the configure script. To see all options, run

   ./configure --help
   

   Most users will simply run

   ./configure --with-orbits-dir=<orbitsdir> --prefix=<installdir>
   

   For example, you might run

   ./configure --with-orbits-dir=/usr/local/orbits
   

5. To build all executables, type

   make
   make install
   

6. test the commands: gmt, esarp, xcorr, conv, gmtsar.csh, etc. If using C-shell you may have to type rehash first. If this does not work then make sure the is in your system $PATH or $path.

RUN

1. GET DATA. There is an example data set at our website: https://topex.ucsd.edu/gmtsar/downloads/

Uncompress the file and then unpack with tar.

2. ORGANIZE the DISK. The standard GMTSAR run has the following directories.

raw SLC topo intf

raw - contains the original data.

SLC - contains the single look complex images derived from the raw data.

topo - contains a digital elevation model for the area in geographic and ultimately radar coordinates.

intf - contains subdirectories with the possible interferograms (only 1 for this example).

3. MAKE DEM

Go to the following web site and construct a file dem.grd that encloses the SAR frame. https://topex.ucsd.edu/gmtsar/demgen

Place the file dem.grd in the /topo directory.

4. EXAMPLE RUNS

There are complete examples for each data type at the following web site: https://topex.ucsd.edu/gmtsar/downloads/

To process a pair interferogram, try p2p_processing.csh, or p2p_S1_TOPS_Frame.csh for Sentinel-1 TOPS data and p2p_ALOS2_SCAN_Frame.csh for ALOS-2 ScanSAR data. For batch processing, see instructions in the GMTSAR documentation (http://topex.ucsd.edu/gmtsar/tar/GMTSAR_2ND_TEX.pdf) or the Sentinel-1 time-series recipe (http://topex.ucsd.edu/gmtsar/tar/sentinel_time_series_5.pdf).

gmtsar's Issues

gmtsar_sharedir.csh is not generated when building GMTSAR

This is somewhat related to #5. When building GMTSAR with CMake, the script gmtsar/csh/gmtsar_sharedir.csh is not created from its template gmtsar/csh/gmtsar_sharedir.csh.in.

preproc_batch_tops_esd.csh doesn't work

Dear Xiaohua,
Due to clear burst boundaries when using preproc_batch_tops.csh, I used preproc_batch_tops_esd.csh. Unfortunately the preproc_batch_tops_esd.csh doesn't work. I used the same input files as preproc_batch_tops.csh.
How to fix this problem.
Below the processing history:
preproc_batch_tops_esd.csh data.in dem.grd 2

Usage: preproc_tops_batch.csh data.in dem.grd mode esd_mode
preprocess and align a set of tops images in data.in, precise orbits required

format of data.in:
image_name:orbit_name

example of data.in
s1a-iw1-slc-vv-20150626...001:s1a-iw1-slc-vv-20150626...001:s1a-iw1-slc-vv-20150626...001:S1A_OPER_AUX_POEORB_V20150601_20150603.EOF
s1a-iw1-slc-vv-20150715...001:s1a-iw1-slc-vv-20150715...001:s1a-iw1-slc-vv-20150715...001:S1A_OPER_AUX_POEORB_V20150625_20150627.EOF

outputs:
baseline.ps align_table.ra (contains info for precise geomatric alignment)
*.PRM *.LED *.SLC(mode 2)

esd_mode:
0 for average, 1 for median, 2 for interpolation

Note:
The names must be in time order in each line to be stitched together

Workflow for ENVI SLC stacking

Hello on the website https://topex.ucsd.edu/gmtsar/downloads/ there are examples for ENVISAT, ENVISAT SLCs and stack ENVI, but no ENVI_SLC stack. It seems the workflow shown in the ENVI stack does not work for the ENVISAT_SLC. Is it possible to stack or do SBAS time series on SLC with GMTSAR? If so, does someone have an example work flow?

• Nathan

Unwrapping Errors

Hello,

After merging the subswaths and sacking the coherence and masking water surfaces, I tried unwrapped using the script unwrap_intf.csh.

I checked the forums for similar issues, but those solutions did not work for me.
For reference, my unwrap_intf.csh script is as follows:

#!/bin/csh -f

intflist contains a list of all date1_date2 directories.

foreach line (awk '{print $1}' intflist)
cd $line
ln -s ../mask_def.grd
snaphu_interp.csh 0.001 0/8618/0/1839
cd ..
end

What am I doing wrong?
Thanks,
Mieszko

Gaps in merged swaths

Hello,

After running the unwrap_int.csh command, I checked the unwrap.pdf files of a few merges and they all have gaps between the swaths. Example:
unwrap.pdf

Looking at the phasefilt.grd files, most of them appear fine but some have gaps.
Example:

What could be causing this issue?

Thanks,
Mieszko

GPS observations to compare displacement

Hello.

This is more of a question than an issue. Sorry to post this here, but after reading the GMTSAR documentation and the associated papers, I couldn't figure out how to proceed.

I want to use GPS observations to compare the displacement obtained from a SBAS workflow of Santo Domingo city (velocity map below; I have the displacement maps too). I also attached the workflow followed, which might help someone else.

I generated the velocity map processing 13 S1 scenes from the period 2019 (August) to 2020 (June). I only have 3 GPS observations of the time period analyzed, and I'm aware that they aren't enough for validation, but those observations are the only ground data I have. My question is: How can I use my GPS data to perform a minimum comparison? I read in one of the SBAS reference papers that the points need to be converted to LOS, so I guess I'll need to use proj_ll2ra_ascii.csh, but I'm not fully aware on how to proceed.

Any insights would be greatly appreciated.

executed.txt

How to extract point values of deformation and plot them as a time series in gmt?

Dear Xiaohua,

I want to create a graph of time-series deformation from output of sbas processing. How do we extract point values of the deformation and plot them as time-series graph?

I appreciate your help.

merge doesn't work

Hello, I have a problem of merging data. I use p2p_S1_TOPS_Frame.csh to process Sentinel-1 data of Idaho earthquake.
The data are:

S1A_IW_SLC__1SDV_20200407T013632_20200407T013656_032015_03B2C9_2AB4.SAFE S1B_IW_SLC__1SDV_20200320T013552_20200320T013619_020769_027618_9910.SAFE
S1A_OPER_AUX_RESORB_OPOD_20200407T035944_V20200406T234559_20200407T030329.EOF
S1B_OPER_AUX_POEORB_OPOD_20200409T110831_V20200319T225942_20200321T005942.EOF

When p2p_S1_TOPS_Frame.csh done, I check the merge directory, I found some space between two swath. like this:

Problem with merge of subwath

Hi, use merge_batch.csh to merge two subwaths, but there is a separation between them, which becomes evident when viewing the correlation and unwrap, this only happens with some interferograms of the stack, also in the terminal no processing error is registered.

cheers.

Error when trying to install GMTSAR (issues with ./configure)

Hello,

I am using currently using Ubuntu 20.04 and have tried solving this problem all day with no luck.

after typing this in the command line:

1. git clone --branch 5.8 https://github.com/gmtsar/gmtsar GMTSAR

2. autoconf

3. ./configure --with-orbits-dir=/usr/local/orbits
   I get the following text...
   issue.txt

If anyone could help me it would be much appreciated!

José

Rank deficiency of matrix in SBAS

Dear GMTSAR user,

I am working on SBAS with GMTSAR, but I feel puzzled because I am informed that the matrix is rank-deficient.

Pre-processing and alignment of SAR images and making interferograms seem to have worked fine, but when I run SBAS like
% sbas intf.tab scene.tab 10 5 2710 2520
then I get the screen output as below, informing that the matrix is rank-deficient.

I believe that the matrix is define regardless of the data (interferogram in this case), so this matrix deficiency might indicate that my scene.tab and intf.tab (see attached) might be wrong. However, I cannot see how I am wrong.

Could you inform me of how I went wrong and how I can fix it?

Best Regards,
Yosuke Aoki
Earthquake Research Institute, UTokyo

---

setting smoothing to 0.000
setting radar wavelength to 0.236 m
setting radar incidence angle to 37.000 degree
setting range to 866000.000 m
Memory Allocation Successful...
read table file ...
number of SAR scenes is 5
number of interferograms is 10
read phase and correlation grids ...
0.000000 0.000102 0.000000 -22.023777
fill the G matrix ...
run least-squares problem over 2710 by 2520 pixel (0) ...
matrix is rank-deficient: 4

write output ...

intf.tab.zip
scene.tab.zip

GMTSAR, gmtlib

Hi
You told me to chech
inside your GMTSAR/config.mk, the GMT_LIB has a path for you to find the gmt library.

I found the attached, but I dont know what shuld i do next

—

SBAS aps

Dear all
I have a problem with my results with ALOS data. I processed 10 interferograms with 1yr temporal baseline in a desert with ocean to the west. i made the batch processing for Level 0 data Alos (2008-2011), then i ran sbas with atm 5 iterations, but my results are wrong in some cm deformation( interseismic periods have a small deformación) , i think the atmospheric delay Is greater than deformation, but my aps results Is 0 in z. I dont understsnd, Could be also some ionospheric delay, or the large temporal baseline in sbas not working?. Thank
Simon

Orbits files for ENVISAT

I am not sure if this is an appropriate place to ask for this kind of help as it is not an issue with the code. If it is not, I will delete this post.

I was hoping someone knew where I could find the orbit files for ENVISAT. I can find the data through SSARA pretty easily, but I have not found a download for the orbits.

On another note, if someone has any pointers on processing with ENVISAT data I would gladly take any notes.

Nathan

mai

is there any way to process MAI with GMTSAR?

SBAS: Dimension don't match!

Hello,

When running sbas, I get the error "dimension don't match! ../merge/2018220_2018268/corr.grd"
When I run gmt grdinfo on that file, I notice that the xdim and ydim are significantly smaller on the unwrap.grd file than on phasefilt.grd or corr.grd. corr.grd and phasefilt.grd have dimensions of 8618x1839 while unwrap.grd has dimensions of 1077x920.

Did something go wrong with my unwrapping step?

My unwrap_intf (converted to txt file) is attached.
unwrap_intf.txt

Thank you,
Mieszko

intf_tops.csh: no landmask is applied with region cut to snaphu

Greetings Dr. Xu, Dr. Sandwell, or anyone else who might know how to help me!

I'm trying to interfere and unwrap multiple Sentinel 1 interferograms with intf_tops.csh. I've been using GMT6.1.1 with GMTSAR 6; we installed the most up to date version of GMTSAR 6 this weekend. This area we are surveying are in the Aleutian Islands. Running this program with a landmask seems important for minimal unwrapping errors.

I can't seem to apply the landmask prior to calling to snaphu_interp.csh. I am seeing lines 193-205 in intf_tops.csh have some commented out sections; as well, landmask.csh doesn't seem to have a $region_cut being applied to it.

I've tried to edit this script by removing the #'s in line 195-197 out;
I've also tried to apply the $region_cut to landmask.csh in line 200
I found an older forum post that seems to highlight the same issue. When I tried this manual method, my landmask_ra.grd was removed from the topo directory during processing, and again no landmask was applied prior to the unwrap.

Is there a known solution, even a manual one, that might help me unwrap with the landmask_ra.grd applied?

Thanks in advance.

Logan, in Fairbanks AK

CMake fails due to missing source files

Hi, I'm trying to compile GMTSAR on Debian 10 (Debian Buster). I've installed all dependencies necessary for compiling GMTSAR. When running cmake -DCMAKE_INSTALL_PREFIX=/usr/local -DCMAKE_BUILD_TYPE=RelWithDebInfo .., cmake fails in the end with:

~~~ Output shortened ~~~

*  Options:
*  GMT library                : /usr/lib/x86_64-linux-gnu/libgmt.so
*  GMT include dir            : /usr/include/gmt
*  TIFF library               : /usr/lib/x86_64-linux-gnu/libtiff.so
*  LAPACK library             : /usr/lib/x86_64-linux-gnu/liblapack.so;/usr/lib/x86_64-linux-gnu/libblas.so /usr/lib/x86_64-linux-gnu/liblapack.so
*
*  Locations:
*  Installing GMTSAR in       : /usr/local

-- Configuring done
CMake Error at preproc/ALOS_preproc/CMakeLists.txt:42 (add_executable):
  Cannot find source file:

    ALOS_filter_SLC/ALOS_filter_SLC.c

  Tried extensions .c .C .c++ .cc .cpp .cxx .cu .m .M .mm .h .hh .h++ .hm
  .hpp .hxx .in .txx


CMake Error at preproc/ALOS_preproc/CMakeLists.txt:42 (add_executable):
  No SOURCES given to target: ALOS_filter_SLC


-- Build files have been written to: /home/user/projects/gmtsar/build

Error in organize_files_tops_linux.csh

Hello.

I tried to run the recipe "Sentinel Time Series Recipe", and found and error in the organize_files_tops_linux.csh script. I used the following command:

tcsh -x /usr/local/GMTSAR/bin/organize_files_tops_linux.csh SAFE_filelist pins.ll 1

The resulting output log is attached. The content of SAFE_filelist is included below (ascending track frames):

/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20190806T224456_20190806T224523_028455_033731_BD53.SAFE
/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20191017T224459_20191017T224526_029505_035B79_CEBA.SAFE
/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20191204T224458_20191204T224525_030205_0373D1_7BCE.SAFE
/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20200214T224455_20200214T224522_031255_039842_87A2.SAFE
/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20200402T224456_20200402T224523_031955_03B0A7_A5B5.SAFE
/media/jose/HD-PNTU3/radar/gmtsar_course/processed/ts_sd/asc/data/S1A_IW_SLC__1SDV_20200613T224459_20200613T224526_033005_03D2AA_1436.SAFE

And this is the content of pins.ll:

-70.033945 18.285100
-70.122071 18.569698

As you may notice in the log file, the script informed that only the last scence in the SAFE_filelist file ... will be re-organized, while for the remaining scenes it returned the message SKIP ######, as it does not have enough scenes. So I decided to run the same command using a SAFE_filelist with only the first scene in it, and the script returned "Frames on date 20190806 will be re-organized...". I guess that it is a bug related with how the SAFE_filelist is built, which impacts the foreach loop.

I'm running GMTSAR on Ubuntu 18.04.

Thanks in advance.

logfile_organize_issue.txt

stack smash error when organize_file_tops_linux.csh runs on ubuntu (SAT_llt2rat)

Found on ubuntu 16.0, 18.0, etc. Unable to reproduce on Mac

From Quanwen Li:

Saving to:
‘S1A_OPER_AUX_POEORB_OPOD_20180311T120551_V20180218T225942_20180220T005942.E
OF’
S1A_OPER_AUX_POEORB_OPOD_2018031 100%
[=========================================================>] 4.21M 931KB/s
in 4.6s
2020-01-24 15:01:16 (931 KB/s) -
‘S1A_OPER_AUX_POEORB_OPOD_20180311T120551_V20180218T225942_20180220T005942.E
OF’ saved [4409491/4409491]
Reading in 5053 lines of info from XML...
Writing 17 lines for orbit...
Writing 58 lines of precise orbit for the LED file...
*** stack smashing detected ***: SAT_llt2rat terminated
*** stack smashing detected ***: SAT_llt2rat terminated
Reading in 5053 lines of info from XML...
Writing 17 lines for orbit...
Writing 58 lines of precise orbit for the LED file...
*** stack smashing detected ***: SAT_llt2rat terminated
*** stack smashing detected ***: SAT_llt2rat terminated

Coregistration Stacking

Hi dear all,
I have a question about the coregistration in the alignment process to Envisat SLC stacking. I follow the alignment with modified align_ALOS_SLC.csh in xcorr parameters, and its working with a secondary alignment slave and supermaster, then I followed intf_tops.csh, with a set the interferograms in intf.in. Then i get the interferograms in folder intf_all, but the results i don't know if they are correct. I compare with the interferogram in the same dates but only one process (p2p_processing ENVI_SLC), with same parameters but that is different wrapped and unwrapped phase. Why i get different results? i think the alignment secondary is different but it's the correct step for the SBAS method or i can follow the SBAS with interferograms made in individual process (p2p_processing). Thank you

How to change to highest resolution in proj_ra2ll.csh instead of 60m default?

Hello @Xiaohua-Eric-Xu

I need to project vel.grd but the default spatial resolution is 60m, Is there any option to change to the highest resolution for Sentinel-1?

Thank you.
Daniel Hermosillo

Problems about SBAS: output vel.grd

Hello Eric,
When I process the SBAS with command:

sbas intf.tab scene.tab 8 5 5788 2309 -range 898971 -incidence 37 -wavelength 0.0554658 -smooth 0.1 -rms -dem

The output files (dem.grd, disp_0000000.grd, rms.grd, vel.grd) contain only zero values.
I noticed your answer in [http://gmt.soest.hawaii.edu/boards/6/topics/7951?r=8288#message-8288].But the problem is still not solved after I do unwrap with the mask_def.grd. And I think many pixel have values in all images after I did this.
I do not know where the problem lies on. I wonder if you could help me with this.
Thank you.

Best
Jiang

mask_def:0.075
snaphu_interp.csh 0.001

sbas log:
compute RMS misfit
compute DEM error

setting smoothing to 0.100
setting radar wavelength to 0.055 m
setting radar incidence angle to 37.000 degree
setting range to 898971.000 m
Memory Allocation Successful...
read table file ...
number of SAR scenes is 5
number of interferograms is 8
read phase and correlation grids ...
0.100000 0.000419 0.000000 -33.408131
fill the G matrix ...
run least-squares problem over 5788 by 2309 pixel (0) ...
matrix is rank-deficient:
write output ...

intf (复件).tab.txt
scene (复件).tab.txt

Error while installing GMTSAR

I got an error message as attached herewith while doing make after configuring the software GMTSAR. Kindly suggest the way to resolve it if any.

INSTALLATION IN MACBOOK AIR

Hey
Please guide me through to install GMTSAR in my MacBook Air

Problems with generation files

Dear developers,
I have a strange problem with GMTSAR. I changed my SO to Ubuntu 20.04 LTS and installed the GMT 6.0 and GMTSAR branch 6, both building the code with cmake. So everything is working, but when I process my data in intf.csh make a error hermite: interpolation not in center interval , then the ifgs are bad. The same process with the same SAR images in a cluster with Ubuntu 18.04 is correct. I also noticed that the baseline.ps is wrong in batch process:

That's still happened with select_pairs.csh, so i see the text created for select_pairs and noticed that date is write with "," and not with "." like a select_pairs in cluster

So when i write manually the first line awk '{print 2014+$3/365.25, $5, $1}' < baseline_table.dat the results is with ","

And for example in cluster with other time series is this

So I think that some mathematics operation with the system is changed from "." to "," then is the wrong interpolation with ifg? or you think that is other issue?. Thanks.

Problem merge frames Sentinel-1

Hi, I joined images using "create_frame_tops.csh", but when doing the interferogram (intf_tops.csh) the phase looks like this. I made sure the dem and orbit were fine. What could be causing the problem?

Bugs: Ionospheric correction on sample InSAR data

Hi Xiaohua,

I ran ionospheric correction on sample InSAR data downloaded from here https://topex.ucsd.edu/gmtsar/downloads/

The p2p_processing.csh failed to generate SLCH, SLCL, params1, and params2 (lines 335-350).

Best,
Rino

folder yyyyddd_yyyyddd format in intf/ or intf_all/ folder being not consistent among satellites

e.g. Sentinel-1/TSX/CMD/RS2 counts ddd from 0 while others from 1.

libgmt

Hello
I installed GMT 6.1 and GMTSAR. The test for gmt worked well, but when i try the following commands to test GMTSAR, it doesnot work.

% esarp

% SAT_baseline

% p2p_ALOS.csh

I receive the following error:

mahsaafra@mahsaafra-Latitude-5400 DING!:~> esarp
esarp: error while loading shared libraries: libgmt.so.6: cannot open shared object file: No such file or directory

It seems there is something wrong with libgmt ...

(also

see attached file)

Intf.in Wrong file names listed

When running "select_pairs.csh baseline_table.dat 50 100" in my F1/F2/F3 directories, I can successfully generate an intf.in file, but rather than having the expected structure, I see the following structure:

s1b-iw1-slc-vv-20180926t172407-20180926t172413-012889-017cd5-004:s1a-iw1-slc-vv-20181107t172449-20181107t172455-024485-02af4c-004
s1a-iw1-slc-vv-20181002t172449-20181002t172455-023960-029dde-004:s1b-iw1-slc-vv-20181113t172408-20181113t172413-013589-019281-004
s1b-iw1-slc-vv-20180809t172405-20180809t172411-012189-01673f-004:s1b-iw1-slc-vv-20180926t172407-20180926t172413-012889-017cd5-004

Where did I go wrong?

Regards,
Mieszko

detrend_before_unwrap.csh

Dear all,

I would like to use "detrend_before_unwrap.csh". I used the command line as
detrend_before_unwrap.csh S1_20170222_ALL_F2 1-5 0/2375/0/1800.

I got a message "F1-5/intf/*/: No match."

I am wondering about input parameters below
detrend_before_unwrap.csh master # bounds

I really appreciate our help.

Thanks

SBAS problem

I got some problem when I run SBAS in GMTSAR. Here is a command I use sbas intf.tab scene.tab 33 16 24452 16212 -wavelength 0.0554658 -rms -dem and the result I got is

compute RMS misfit
compute DEM error

setting smoothing to 0.000
setting radar wavelength to 0.055 m
setting radar incidence angle to 37.000 degree
setting range to 866000.000 m
memory allocation! phi

prep_data_linux.csh working in F1 and F2 but not F3

Hello all,

I have been working on a small time series for GMTSAR practice and I am getting a strange error. I am up to the point in which I need to merge the subswaths for F1 and F2. When I try to get things rolling in the F3 directory with prep_data_linux.csh the script ends early. Here is the ouput of the script

(base) nathan@nathan-MS-7A93:/media/nathan/massstorage/GMTSAR_DATA/Salton_sea_proc/F3/raw$ prep_data_linux.csh 
rm: cannot remove 'data.in': No such file or directory
--2020-07-28 15:28:46--  https://s1qc.asf.alaska.edu/aux_poeorb/
Resolving s1qc.asf.alaska.edu (s1qc.asf.alaska.edu)... 137.229.87.170, 137.229.87.163
Connecting to s1qc.asf.alaska.edu (s1qc.asf.alaska.edu)|137.229.87.170|:443... connected.
HTTP request sent, awaiting response... 503 Service Unavailable
2020-07-28 15:28:47 ERROR 503: Service Unavailable.

Looping over all the lines
Setting a new name for rec
Writing record 20191003
if: Missing file name.

F3_ls.txt

Here is the contents of my directory. Notice that the orbits.list and text.dat is created, but never removed. Now, I do not have wgetasf set up, but this was not a problem in the previous directories that have the exact same setup. I also added the hacky script I have to set up the directories in case that reveals anything.
create_dirs.txt

Vertical displacement vector with LOS of ascending and descending orbit

Hi. I have investigated ways to calculate the vertical displacement vector using the LOS data of both orbits, but I have not seen any solution to the inversion of the matrices that are proposed in the articles.
Wen S., et al. 2016 talks about a simpler solution but the document is in Chinese. Do you have any suggestion?

prep_batch_tops.csh

Dear all,

I have a problem with prep_batch_tops.csh, getting errors as below

preproc_batch_tops.csh data.in dem.grd 2
grdfilter (api_alloc_grid): Error: Could not reallocate memory [33527.74 Gb, 9000034800016 items of 4 bytes]
grd2xyz: File flt.grd not found
[Session gmt (0)]: Error returned from GMT API: GMT_FILE_NOT_FOUND (16)
rm: No match.

How do I fix the errors.

I appreciate for helps.

preproc_batch_tops.csh produces wrong baseline_table.dat

Error when defining region_cut in intf_tops.csh...and possible solution

I needed to batch generate interferogram and unwrapped phase grids for a small region of a single subswath, so I tried to run intf_tops.csh using the region_cut parameter. For this, I calculated the radar coordinates of an AOI defined by lon/lat/height using SAT_llt2rat. The script performed well until it reached the snaphu part, which returned an error regarding an invalid region. I guess the error was due to the fact that the same region_cut parameter defined in the config file, applied also to the snaphu processing, so when snaphu tried to unwrap the phase of the already cut image, the region was no longer valid.

In order to continue with the workflow, I tried a solution but I am not sure if it's consistent. I forked the gmtsar repo and there I added a new parameter called region_cut_snaphu to both intf_tops.csh and batch_tops.config files (in the latter, I left it blank). I'm not sure if this is a consistent solution, or if it will "play well" with the rest of the scripts, but at least in my case I could unwrap the phase of all my interferograms.

Any feedback would be greatly appreaciated.

ERS_SLC or ENVI_SLC for SBAS

Dear Xiaohua(Eric):
I want to process ENVI_SLC and ERS_SLC for interseismic periods in Mejillones Peninsula, Chile, with time-series SBAS method. In GMTSAR there is only the option to make Sentinel 1, ALOS_SLC, ENVI or ERS, but for SLCs not working. How do I make this possible?

error installing GMTSAR: no file tiffio.h

I am trying to install GMTSAR but I am running into a problem at the "make" step. It seems that there is no file or directory "tiffio.h".

I tried to use apt to install libtiff and it says that it cannot be found. I have version 5.4.3 of gmt installed with apt, and I am on Ubuntu 18.04.

Any help is appreciated!

Cheers,

Mitch

S1A problm error - could not open file S1_20200616_150043_F1.PRM

Generate DEM error

I have an area which cross 0 Latitude,
When I try to generate DEM from https://topex.ucsd.edu/gmtsar/demgen/ as the instruction says,
It gives me an error which says on the status "Error, Area cannot cross 0 Latitude"

How to solve this? Is there any workaround for this?
Thank You

Suggested improvement

On multiple occasions, while running some of the scripts of GMTSAR (e.g. intf_tops.csh), I receive the error "If expression syntax". I recently realized that this is because I use to edit one parameter while preserving a copy of it commented in the config file. For example, something like this:

#region_cut = 20/30/100/110
region_cut = 10/20/90/100

I guess that the easy solution is avoid adding copies of the parameters. However, I guess that adding a "^" symbol before the current pattern of grep expressions, would prevent many errors like the one commented (hopefully) without affecting the rest of the code. For example, in intf_tops.csh, the grep expression that captures the region_cut parameter could be rewritten to something like this:

set region_cut = `grep '^region_cut' $2 | awk '{print $3}'`

error with gmt movie: FFmpeg not found

I installed gmt version 6.1 on Friday and the first test case worked fine, however, when I ran the second one I got an error with gmt movie telling me that ffmpeg is not installed. I uninstalled and purged it from my machine and then tried to reinstall it freshly but that did not fix the problem. You can also see that if I try to just run ffmpeg from the command line I get an error with a shared library but I am not sure if that is the problem or not since gmt seems to think it is not installed.

I installed gmt using anaconda and I am currently running on Ubuntu 18.04. Any help is appreciated!

Cheers,

Mitch

gmt lib

When I do ls-ltG, there are many files. Could you please tell which files should be executable? Only the highlighted one in the attached?
Thanks

Merge frames ERS/ENVISAT

HI Xiaohua,
I have a question about the merge frames of one track. I attached a screenshot of two frame phases, but this don't match totally. I processed each frame in the same parameters. Do you know anything for a better match? or any script to make one interferogram of two frames in one process with these satellites.
Thank you,
Simon

Ionosphere correction for TOPS data may require a flattened spectrum along range.

One way is to flatten the spectrum. The other is to re-compute the corresponding central frequency for band-passed SLCs.

error in resamp while processing the SAR data

I am using GMTSAR version 5.8 in ubuntu 14.04.
I have earlier process many SAR data.
But recently i am getting error in resamp.
resamp: error while loading shared libraries: libgfortran.so.1: cannot open shared object file: No such file or directory.
i am also attaching the log file.
please let me know where i am wrong.
Sanjay
log.txt

preproc_batch_tops.csh data.in dem.grd 1 problem

HI Xiaohua,

I have a problem processing Sentinel-1B images in "preproc_batch_tops.csh data.in dem.grd 1". During the this step, I got like this ：

"Reading in 5464 lines of info from XML...
Writing 16 lines for orbit...
Writing 6 lines of precise orbit for the LED file...
Successfully opened junk1
hermite: interpolation not in center interval
hermite: interpolation not in center interval
hermite: interpolation not in center interval".

I checked my EOF, they do cover the images.

data.in:
s1b-iw1-slc-vv-20190808t220938-20190808t220941-017500-020e97-004:S1B_OPER_AUX_POEORB_OPOD_20190828T110923_V20190807T225942_20190809T005942.EOF
s1b-iw1-slc-vv-20191019t220941-20191019t220943-018550-022f3e-004:S1B_OPER_AUX_POEORB_OPOD_20191108T110545_V20191018T225942_20191020T005942.EOF

"organize_files_tops_linux.csh SAFE_filelist pins.ll 1" failed after I changed the oreder of lon and lat in the pins.ll file, so I guess the original one is correct.

original pins.ll:
120.58 49.24
120.45 49.18

Could you let me know how to deal with this problem? Look forward to your reply.

Best
Jiang

Some file name length is limited, e.g. make_s1a_tops.c tmp_str

This usually fails some codes when you have a "path+name" being 200+ characters long