doris提取DEM脚本（ERS）

来源：互联网发布：js中连接符后怎么使用编辑：程序博客网时间：2024/05/16 15:47

########################         第一部分运行脚本        ####################################
1、run -g

2、run -s1 %读取主影像
3、run -s2 %读取辅影像
4、run -s3 %配准
5、run -s4 %重采样
6、run -s5 %生成干涉图
7、run -s6 % 相位解缠
8、run -s7 %地理编码
9、合并height.raw lambda.raw phi.raw

#########################     第二部分 run -s6 相位解缠出错解决方案 #######################

注：run -s6相位解缠时可能出错。方法：修改snaphu_io.c第1474行，将%s改为%d

fprintf(fp,"# Log file generated %d",ctime(t)); /*IvanLJF %s--->%d :fprintf(fp,"# Log file generated %s",ctime(t));*/

修改后重新编译,运行下面的snaphu命令。

snaphu -f snaphu.conf Outdata/master_slave.cint.filtered 2250

#################### 第三部分    合并height.raw lambda.raw phi.raw      #################################

1、lonlathei2ascii Outdata/lambda.raw Outdata/phi.raw Outdata/heights.raw   lonlathei.dat

################### 第四部分    用GMT显示地理编码后的经纬度高程          ################

1、ascii2grdAutoMinMaxHDR lonlathei.dat

2、ascii2ps lonlathei.dat

# lonlathei.dat.nvi lonlathei.dat.hdr可用ENVI打开

################## 第五部分    用GMT显示地理编码后的UTM投影下的xyz坐标    #############################
1、ascii2ascii_UTM lonlathei.dat +utm_zone=15 xyz_utm.dat

2、ascii2grdAutoMinMaxHDR xyz_UTM.dat

3、ascii2ps xyz_UTM.dat

0 0