1.参考网址:https://github.com/scottyhq/insar_scripts/tree/master/TSX
2.修改.bashrc,将sar-0.6的make_slc_tsx可执行路径/usr/local/bin提前,不然会调用GMT5SAR的make_slc_tsx同名可执行程序.
3.将下载的get_height.pl放到ROC_PAC的scr目录下:/home/lu/StaMPS/StaMPS_v3.3b1/ROI_PAC_3_0_1/ROI_PAC/INT_SCR
4.工作目录下放两景TSX数据.执行,用prepare_tsx.py执行以下命令
$./prepare_tsx.py ./dims_op_oc_dfd2_203830440_1
$./prepare_tsx.py ./dims_op_oc_dfd2_204492801_1
5.准备DEM
$mkdir DEM
整理DEM文件夹下的文件,编辑dem.hdr dem.rsc文件内容.
dem.hdr内容如下
set_name (&) = DTE
set_rows (pixels) = 6000
set_cols (pixels) = 6000
set_hddr (bytes) = 0
set_tail (bytes) = 0
val_size (bytes) = 2
val_frmt (&) = INTEGER*2
val_ordr (&) = LITTLE ENDIAN
dte.row_hddr (bytes) = 0
dte.row_tail (bytes) = 0
dte.col_hddr (bytes) = 0
dte.col_tail (bytes) = 0
dem.rsc内容如下:
FILE_DIR /media/lu/94FE61A5FE617FFA/04_ROI_PAC_Demo/DEM
WIDTH 6000
FILE_LENGTH 6000
XMIN 0
XMAX 6000
YMIN 0
YMAX 6000
X_FIRST 115
Y_FIRST 40
X_STEP 0.000833333
Y_STEP -0.000833333
X_UNIT degres
Y_UNIT degres
Z_OFFSET 0
Z_SCALE 1
PROJECTION LATLON
5.编辑int.proc文件内容
SarDir1=090327
SarDir2=091205
IntDir=int_090327_091205
#DEM=/media/lu/94FE61A5FE617FFA/04_ROI_PAC_Demo/DEM/tianjin_srtm.dem
#flattening=topo
OrbitType=HDR
pixel_ratio=1
很重要的一个参数pixel_ratio=1,对TSX来说,该值为1才能保证距离向,方位向分辨率大致相同.
对ERS来说,pixel_ratio=5,因为ERS方位向像素大小为4m,距离向地距分辨率为20m.(20/4=5)
其他参数设置参考ROI_PAC手册P33.
注:从baseline.pl知,$orbit_type =~ /(ODR|PRC|HDR)/ or die “Orbit type must be ODR, PRC or HDR\n”;
注:注释掉了DEM flattening两个参数,因为模拟影像不对,考虑到影像范围也无地形
注:若要重新执行,需要将GEO SIM int_090327_091205均删除
$process_2pass.pl int.proc roi_prep orbbase
$process_2pass.pl int.proc slcs done_sim_removal
mdx.pl 显示