GMTSAR的D-InSAR基本处理过程

使用总的p2p，即可得到D-InSAR 干涉处理结果，如下图：处理步骤及文件夹放置位置

1.准备raw文件夹

在ASF官网下载哨兵数据和哥白尼官网下载轨道数据

ASF下载哨兵数据官网

esa下载数据官网

2.准备topo文件

在gmtsar官网下载对应的dem数据（4度以内）

 gmtsar官网dem

或者gmt下载dem

make_dem.csh

3.准备配置文件

如下图

pop_config.csh S1_TOPS > config.s1a.txt

 得到的配置文件内容如下

#
# This is an example configuration file for p2p_S1_TOPS.csh
#
# all the comments or explanations are marked by "#"
# The parameters in this configuration file is distinguished by their first word so 
# user should follow the naming of each parameter.
# the parameter name, "=" sign, parameter value should be separated by space " ". 
# leave the parameter value blank if using default value. 
#
  
#####################
# processing stage  #
#####################
# 1 - start from preprocess
# 2 - start from align SLC images
# 3 - start from make topo_ra 
# 4 - start from make and filter interferograms 
# 5 - start from unwrap phase
# 6 - start from geocode  
proc_stage = 1
skip_stage = 

##################################
#   parameters for preprocess    #
#   - pre_proc.csh               #
##################################
# num of patches
num_patches = 

# earth radius 
earth_radius =

# near_range
near_range = 

# Doppler centroid 
fd1 = 

################################################
#   parameters for focus and align SLC images  #
#   - align.csh                                #
################################################
# region to cut in radar coordinates (leave it blank if process the whole image)
# example 300/5900/0/25000
region_cut =

#
#####################################
#   parameters for make topo_ra     #
#   - dem2topo_ra.csh               #
#####################################
# subtract topo_ra from the phase
#  (1 -- yes; 0 -- no)
topo_phase = 1
# if above parameter = 1 then one should have put dem.grd in topo/

# topo_ra shift (1 -- yes; 0 -- no)
shift_topo = 0

####################################################
#   parameters for make and filter interferograms  #
#   - intf.csh                                     #
#   - filter.csh                                   #
####################################################
# switch the master and slave when doing intf. 
# put "1" if assume master as repeat and slave as reference 
# put "0" if assume master as reference and slave as repeat [Default]
# phase = repeat phase - reference phase
switch_master = 0

# filters 
# look at the filter/ folder to choose other filters
# for tops processing, to force the decimation factor
# recommended range decimation to be 8, azimuth decimation to be 2
filter_wavelength = 200

# decimation of images 
# decimation control the size of the amplitude and phase images. It is either 1 or 2.
# Set the decimation to be 1 if you want higher resolution images.
# Set the decimation to be 2 if you want images with smaller file size.
# 
dec_factor = 2 
range_dec = 8
azimuth_dec = 2
#
# make ionospheric phase corrections using split spectrum method
correct_iono = 0
iono_filt_rng = 1.0
iono_filt_azi = 1.0
iono_dsamp = 1
# 
# set the following parameter to skip ionospheric phase estimation
iono_skip_est = 1 
#
#####################################
#   parameters for unwrap phase     #
#   - snaphu.csh                    #
#####################################
# correlation threshold for snaphu.csh (0~1)
# set it to be 0 to skip unwrapping.
threshold_snaphu = 0.1

# interpolate masked or low coherence pixels with their nearest neighbors, 1 means interpolate, 
# others or blank means using original phase, see snaphu.csh and snaphu_interp.csh for details
# this could be very slow in case a large blank area exist
near_interp = 1

# mask the wet region (Lakes/Oceans) before unwrapping (1 -- yes; else -- no)
mask_water = 1

#
# Allow phase discontinuity in unrapped phase. This is needed for interferograms having sharp phase jumps.
# defo_max = 0 - used for smooth unwrapped phase such as interseismic deformation
# defo_max = 65 - will allow a phase jump of 65 cycles or 1.82 m of deformation at C-band
#
defomax = 0

#####################################
#   parameters for geocode          #
#   - geocode.csh                   #
#####################################
# correlation threshold for geocode.csh (0< threshold <=1), set 0 to skip
threshold_geocode = 0.1

4.运行p2p_S1_TOPS_Frame.csh命令

运行这个命令即可

p2p_S1_TOPS_Frame.csh S1A_IW_SLC__1SDV_20210520T112620_20210520T112704_037971_047B4F_0A0A.SAFE S1A_OPER_AUX_POEORB_OPOD_20210609T121911_V20210519T225942_20210521T005942.EOF S1B_IW_SLC__1SDV_20210526T112538_20210526T112628_027075_033C0F_3031.SAFE S1B_OPER_AUX_POEORB_OPOD_20210615T111837_V20210525T225942_20210527T005942.EOF config.s1a.txt vv 1

 准备数据图片如下，readme.sh即上面终端的命令，可以直接在终端运行，也可以sh readme.sh运行

5.输出结果主要在merge文件夹

merge主要输出如下图的文件

主要看以下几个文件 ，后缀带ll的为地理编码以后的经纬度坐标，不带ll的是在雷达坐标系下的坐标

corr_ll.grd 相干性文件

phasefilt_mask_ll.grd 缠绕图

unwrap_ll.grd 解缠图

los_ll.grd  los向位移图 

unwrap_ll.grd转换到los_ll.grd公式为

grdmath unwrap_mask_ll.grd 0.0554658 MUL -79.58 MUL = los_ll.grd

-4.413968364

gmt grdfilter los_ll.grd -Glos_ll_filter.grd -D1 -Fg10+h -V

参考网站：GMTSAR