PIE-SAR软件利用TanDEM-X干涉对（差分相位补偿法）获取RAW DSM操作流程

1.数据情况说明

TDX干涉对文件名：TDM1_SAR__COS_BIST_SM_S_SRA_20121205T101835_20121205T101843

地理位置：湖北省

外部参考DEM:注意事项：使用全国30m SRTM DEM做出的结果精度不对，采用ALOS DEM。DEM需要做高程异常转换.

2.差分补偿技术流程

1. 数据导入（导入主辅影像）
2. 粗配准生成偏移文件
3. 多视、地理编码（生成模拟影像，精化查找表，RDC DEM用来模拟地形相位）
4. 获取残余地形相位：干涉图计算、基线计算、平地相位去除、模拟地形相位，相位加减（去平后的相位减去模拟的地形相位，得残余地形相位），对残余地形相位进行ADF滤波
5. 相位解缠：对干涉图、相干图进行2倍过采样、利用干涉图生成掩膜文件、最小费用流相位解缠、对解缠结果进行降采样、利用unw_model命令对原始差分干涉图解缠、拟合残余地形相位趋势面、去除趋势面，权重插值。
6. 相高转换：模拟地形相位加上解缠的残余地形相位，选择GCP点，基线精化，相位转高程
7. RAW DSM融合
8. 质量检查，生成精度报告。

3. PIE-SAR差分补偿处理流程

3.1 DEM高程异常转换

下载的SRTM DEM、AW3D DEM、Copernicus DEM（World dem）均是正常高，需要基于EGM96或EGM2008高程异常模型转到大地高。

外部参考DEM用来对主影像进行几何校正处理，高程异常会影响平面定位精度。另外，外部参考DEM也会用来模拟地形相位，未考虑高程异常，会对最终的RAW DSM高程精度有影响。

3.2 数据导入

采用TSX/TDX批量导入功能，同时导入主辅影像，选择输出【原始影像名】。因干涉对是一发双收模式，成像时间一样。在输出文件夹建一个data_import文件夹作为输出目录。

3.3 地理编码

按输出分辨率10米，设置多视视数据为距离向×方位向=4×5。确保后续主影像多视视数，和生成干涉图的多视视数保持一致。

3.3.1 多视处理

多视处理，可选择自定义视数，手动设置距离向视数为4，方位向视数为5。也可选择输入栅格格网大小：10m。

3.3.2 地理编码处理

对多视后的主影像进行地理编码（GTC，地形校正），选择EGM2008大地高外部参考DEM，输出10米分辨率，选择CGCS2000坐标系，其他产品勾选【模拟幅度图】，用于精化查找表。勾选【叠掩与阴影】，用于与外部DEM融合。

3.4 获取雷达坐标DEM

3.4.1 坐标转换（地图坐标->雷达坐标）

将模拟幅度图由地理坐标转到雷达坐标，用于后续配准生成偏移文件，用来精化查找表。因需要从雷达坐标下的文件来获取行列号，因此需要输入雷达坐标多视文件。

【注】：生成的距离多普勒坐标模拟幅度图，行列号与多视主影像行列号相同。

3.4.2 幅度配准

生成的距离多普勒坐标模拟幅度图与多视主影像行列号相同，两者存在方位向和距离向上存在一定的偏差，通过幅度配准算法，求得两者在距离向、方位向的偏移量，用于精化查找表，提升无控定位精度以及地形相位模拟精度。

3.4.3 精化查找表

利用幅度配准得到的偏移文件，对查找表进行精化处理，提升定位精度。

3.4.4 坐标转换（雷达坐标->地图坐标）

利用精化后的查找表重新生成雷达坐标系下的模拟幅度图（只是为了验证一下是否配准好，若配准没问题，可近一步将segment_utm_dem.img转到雷达坐标系下，用于模拟地形相位）

1.利用精化查找表将模拟幅度图转到雷达距离多普勒坐标下，卷帘查看与多视主影像是否贴合，确保幅度配准计算的偏移量是正确的。

2.若雷达距离多普勒坐标的模拟幅度图与多视主影像贴合度没问题，则将SegmentDEM转到距离多普勒坐标系下，用于模拟地形相位。

3.5 获取残余地形相位

3.5.1 粗配准

虽然TDX干涉对已经完成了精配准，但在干涉处理过程中多处都要用到偏移文件，偏移文件存储了干涉图相关的信息，包括干涉图的行列号，多视视数等。单视主影像参数文件配合偏移文件，可作为干涉图的参数文件。

注意事项：

• 工作模式选择了双基成像。对于TH-02、HT-01、TSX/TDX卫星，采用一发双收、多收的成像模式，选择双基成像，对于GF-3重轨卫星，选择单基成像。当然此处仅是为了获取偏移文件，用不到计算的偏移量结果。
• 主辅星不要选反了，不然计算出来的偏移量正负号就会反。

3.5.2 干涉图计算

最终输出的RAW DSM分辨率为10米，设置多视视数为距离向* 方位向==4*5。

3.5.3 提取相位

3.5.4 基线计算

通过主辅影像的轨道矢量，计算干涉对的初始基线，基线文件包含影像中心时刻对应的基线（TCN坐标），以及基线的线性形变速率，用于求解每行的基线矢量。

3.5.5 去除平地相位

去除平地相位，得到地形相位。将模拟的平地相位也同步输出，干涉定标有用到。

3.5.6 地形相位模拟

Notice, that the baseline geometry described in the baseline file is may not be perfect and have residual quadratic phase trends.

3.5.7 相位相减

去平后的干涉图，减去模拟的地形相位，得残余地形相位。

3.6 相位解缠

3.6.1 自适应谱滤波

对残余地形相位进行自适应谱滤波，并获得相干图。

3.6.2 生成掩膜文件

3.6.3 MCF相位解缠

3.6.4 权重插值

3.7 相位补偿

3.7.1 趋势面拟合

3.7.2 相位补偿

模拟地形相位加上解缠残余地形相位，得到用于相高转换的地形相位。

3.8 获取RAW DSM

3.8.1 相高转换

3.8.2 坐标转换

3.8.3 DSM融合

InSAR干涉对获取的RAW DSM，在叠掩、阴影区存在拉花现象，需要用外部参考DEM补充。在融合前，先将地理编码生成的叠掩阴影文件转为掩膜文件，并做形态学运算，对边界进行外扩，融合完做中值滤波。

3.9 精度检查

3.9.1 提取ICESAT控制点

3.9.2 生成精度报告