GAMMA2stamps(全分辨率)

#滤波处理

mk_adf_2d rslc_tab itab rmli/MLI_ave diff0 7 0.4 64 4

#################################### 生成经纬度文件 ##############################

建立经纬度文件

(1)生成pointlist

mkgrid plist 5800 2000 1 1 1 1

(2)将栅格转为点格式

data2pt cop30.rdc.hgt 20231227.rmli.par plist 20231227.rmli.par pdem 1 2
# cop30.rdc.hgt 雷达坐标系下的dem

(3)生成经纬度数据

pt2geo plist - 20231227.rmli.par - pdem cop30.utm.par diff_par.in 1 1 plist_map pmap_coord plat_lon phgt 0 pmask_dem
# cop30.utm.par utm 坐标系下dem的par文件

(4)swap

swap_bytes plat_lon plat_lon4 4

(5)python拆分

python3 lonlat.py

‘’’

import numpy as np

cors = np.fromfile('plat_lon4', np.float32)

lon = cors[0:len(cors):2]
lat = cors[1:len(cors):2]

lon.tofile('lon.raw')
lat.tofile('lat.raw')

‘’‘’

(6)swap

swap_bytes lon.raw lon4.raw 4
swap_bytes lat.raw lat4.raw 4

#################################### 转到stamps ##############################
#准备格式
INSAR_20231227
./SMALL_BASELINES/20231227.rlsc
./SMALL_BASELINES/20210120.rslc
./SMALL_BASELINES/20231227_20210120.diff
./SMALL_BASELINES/20231227_20210120.base

./geo/*.den_rdc # 编码查找表
./geo/20231227.lon
./geo/20231227.lat

./dem/*_seg.par # utm坐标系下的dem的par文件

mkdir SMALL_BASELINES

创建文件夹

run_all.pl ../bperp_file 'mkdir SMALL_BASELINES/$2_$3'  

#复制文件
run_all.pl ../bperp_file 'ln ../diff0/$2_$3.adf.diff SMALL_BASELINES/$2_$3/$2_$3.diff' 
run_all.pl ../bperp_file 'ln ../diff0/$2_$3.base SMALL_BASELINES/$2_$3/$2_$3.base'


run_all.pl ../bperp_file 'ln ../rslc/$2.rslc SMALL_BASELINES/$2_$3/$2.rslc'
run_all.pl ../bperp_file 'ln ../rslc/$2.rslc.par SMALL_BASELINES/$2_$3/$2.rslc.par'
run_all.pl ../bperp_file 'ln ../rslc/$3.rslc SMALL_BASELINES/$2_$3/$3.rslc'
run_all.pl ../bperp_file 'ln ../rslc/$3.rslc.par SMALL_BASELINES/$2_$3/$3.rslc.par'

pwd # 查看绝对路径

mt_prep_gamma 20231227 /media/estar/SSGC/INSAR/GAMMA2STAMPS_singlelook/INSAR_20231227 0.6 6 2 50 50 

对于前面dem文件有不理解的同学，可以B站视频了解（与上方同名）

经纬度的文件整理相对繁琐，这里有博主整理的便捷bsh脚本，下载脚本后在终端直接prep_lonlat即可一键实现文件的转换（有疑问可通过主页的联系方式进行咨询）。

https://pan.baidu.com/s/1ym9e2u1W3cVwt_NZEJ2OMQ?pwd=xhrp