#滤波处理
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