详细介绍：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