doris提取DEM脚本（ERS）

########################         第一部分  运行脚本        ####################################
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