GAMMA简易教程-第三部分-差分干涉反演形变
差分干涉部分-形变反演
事实上，GAMMA里有两种做差分的思路，一种是直接做差分干涉，还有一种是先做干涉，再做差分。第二种的好处是有一些干涉的其他输出，可以用于高程反演等任务。本流程承接之前的任务（高程反演），在完成干涉的基础上，进行差分。
首先，需要下载研究区域的DEM数据，网站：EarthExplorer 可以免费下载。
本次实验下载了四个DEM：
进一步，需要对DEM进行重采样，生成对应雷达坐标的相位。
-
使用 srtm2dem 命令将每个 .tif 文件转成内部格式（例如 .DEM 文件和对应的 .par 参数文件）。
输入：1.下载的DEM.Tif
输出：1.GAMMA内部格式的.DEM 2.DEM对应的参数文件.par
样例:srtm2dem n35_w112_1arc_v3.tif n35_w112_1arc_v3.dem n35_w112_1arc_v3.dem_par 3
-
DEM的拼接
每个DEM完成导入后，需要拼接成一个大的DEM文件，生成对应参数文件，使用命令mosaic：
输入：1. DEM数量 2. 各个DEM文件与对应参数文件
输出：1.合并的DEM 2.合并DEM的参数文件
样例：mosaic 4 n36_w113_1arc_v3.dem n36_w113_1arc_v3.dem.par n36_w112_1arc_v3.dem n36_w112_1arc_v3.dem.par n35_w113_1arc_v3.dem n35_w113_1arc_v3.dem.par n35_w112_1arc_v3.dem n35_w112_1arc_v3.dem.par mosaic_outputDEM.dem mosaic_outputDEM.dem.par 0 0
-
清洗 DEM 值并进行插值处理（可选）
本部主要针对存在海岸等特殊场景的DEM进行处理，包含以下操作：
1从 demutm.par 中提取 DEM 宽度（行像素数量），为后续处理提供尺寸信息。
set width_dem=`grep width: demutm.par |awk '{print $2}'`
2 将 demutm 中所有值为 -100 的像元替换为 0.01（常用于将异常或空值替换成合理的小值）。
replace_values demutm -100 0.01 demutm_temp1 $width_dem 2
3 对 demutm_temp1 图像进行自适应插值，填补空值区域，输出为 demutm_temp2。
interp_ad demutm_temp1 demutm_temp2 $width_dem
4 再次处理，将值为 0 的区域换成 0.001，生成最终处理后的高程图 optdem。
replace_values demutm_temp2 0 0.001 optdem $width_dem 2
5 ashgt 是显示或计算高程的命令
rashgt optdem - $width_dem
本例没有进行DEM的特殊处理，输出高程图如下（mosaic_outputDEM.dem）：
-
主影像与DEM的几何映射（粗略几何编码）
gc_map 是初始地理编码工具：
参数说明：
输入：
1.主影像（master ima
