SARscape D-InSAR实战：从哨兵1号数据到地表形变图全流程解析

1. 从数据到形变图：D-InSAR处理的核心脉络

如果你手头有一堆哨兵1号（Sentinel-1）的SAR数据，想知道某个区域的地表是不是在悄悄下沉或者抬升，比如监测矿区沉降、城市地面沉降或者火山活动，那么D-InSAR（差分合成孔径雷达干涉测量）就是你绕不开的技术。而SARscape，作为商业遥感软件中的“瑞士军刀”，把这一整套复杂的流程封装得相对友好。但“友好”不代表没有坑，从数据下载、预处理到最终生成形变图，每一步的选择都直接影响结果的可靠性和精度。我自己在项目里反复折腾过多次，从最初的“跑通就行”到后来的“抠细节求精度”，积累了不少实战心得。这篇文章，我就以哨兵1号数据为例，拆解一遍SARscape里D-InSAR的标准处理流程，重点不是复述菜单怎么点，而是告诉你每一步背后的逻辑、为什么要这么设置参数，以及那些手册里不会写、但能让你少走弯路的经验。

整个流程可以看作一条清晰的流水线：数据准备 -> 生成连接图 -> 干涉处理 -> 轨道精炼与重去平 -> 相位解缠 -> 形变反演与地理编码。听起来步骤不少，但核心思想就一个：利用两景或多景SAR图像之间的相位差，剔除地形和平地效应，最终提取出微小的地表形变信息。哨兵1号数据免费、重访周期短（6天/12天），是做大范围、长时间序列形变监测的绝佳数据源。而SARscape的优势在于它把诸如配准、去平、相位解缠这些高难度动作都做成了相对标准化的模块，让研究者能更专注于结果分析和地质解释。不过，软件再智能，参数理解不到位，出来的结果可能就是“垃圾进，垃圾出”。

2. 流程起点：数据准备与工程环境搭建

在打开SARscape第一个模块之前，充分的准备工作能避免后续90%的莫名报错。这个阶段的核心就两件事：搞对数据，设好环境。

2.1 哨兵1号数据的选择与下载

不是所有的哨兵1号数据都适合做D-InSAR。你需要关注以下几个关键属性：

1. 轨道方向与模式：对于形变监测，优先选择升轨（Ascending）或降轨（Descending）的干涉宽幅（IW）模式数据。IW模式是哨兵1号默认的地面观测模式，提供稳定的、足够宽的幅宽。尽量避免使用波模式（Wave）或超幅宽模式（EW），除非你的研究目标特殊。
2. 极化方式：VV或VH极化？对于地表形变监测，尤其是在城市、裸露岩石区，VV极化信号的后向散射更强，相干性通常更好，是首选。VH极化更多用于植被分类等。
3. 时间基线与空间基线：这是干涉成功的生命线。
   • 时间基线：两景图像拍摄的时间差。理论上越短越好，因为地表特征（尤其是植被、土壤湿度）变化小，相干性高。对于哨兵1（6天重访），时间基线通常是6天、12天、18天等的倍数。尽量选择时间基线短的像对。
   • 空间（垂直）基线：两颗卫星在垂直于飞行方向上的距离。基线越长，对地形越敏感，但相位噪声也越大，容易失相干。SARscape在处理时会计算并建议一个合理的基线阈值（通常几百米以内），超出阈值的像对干涉质量会很差。
4. 数据级别：你需要下载的是SLC（Single Look Complex）数据。这是单视复数数据，保留了完整的振幅和相位信息，是进行干涉处理的唯一选择。千万不要下成GRD（地距检测）数据，那已经丢失了相位。

下载渠道：欧空局的哥白尼数据空间（Copernicus Data Space Ecosystem）或阿拉斯加卫星设施（ASF）都是可靠的免费源。下载时，确保将整个.SAFE格式的文件夹完整保存。一个常见的坑是：下载的数据缺失了部分辅助数据文件（如精确轨道星历文件），这会导致后续轨道信息不准确。通常，SARscape或配套的下载工具能帮你自动补下载精确轨道文件（Precise Orbit Files），务必确保这一步成功。

2.2 SARscape工程与参数预设

打开SARscape，第一步通常是创建一个处理链或直接在各模块中设置工作目录。我强烈建议你为每个项目建立独立的、路径中不含中文和特殊字符的文件夹。目录结构可以这样组织：

/Your_Project/ ├── /Data/ # 存放原始的.SAFE数据 ├── /Orbit/ # 存放下载的精密轨道文件 ├── /Processing/ # SARscape各模块输出目录 │ ├── /Coregistration/ │ ├── /Interferogram/ │ └── ... └── /Output/ # 最终成果图

在**/Interferogram Formation**（干涉图生成）或**/Interferometric Stacking**（干涉堆栈）模块中，你需要首先设置一些全局参数：

• 多视数（Looks）：为了抑制斑点噪声，需要对图像进行多视处理。Range Looks（距离向视数）和Azimuth Looks（方位向视数）的设置会影响最终图像的分辨率和像素大小。一个常见的起始设置是5:1（距离向:方位向）。这个比例大致对应于使输出像素接近正方形（例如，约20m x 20m）。增加视数可以提高信噪比和相干性，但会损失空间细节。你需要根据研究区域的尺度和细节要求做权衡。
• DEM选择：D-InSAR需要高精度的数字高程模型（DEM）来去除地形相位。SARscape内置支持SRTM（90米/30米）、ASTER GDEM等。对于一般区域，SRTM 1Sec（约30米）通常足够。对于山区或需要高精度的情况，可以考虑使用更高分辨率的DEM（如ALOS World 3D）。务必确保DEM覆盖你的整个研究区，且没有空洞。
• 输出像素大小：这通常由多视数决定，软件会自动计算。你可以接受默认值，也可以根据需求微调。

注意：在处理大面积区域或长时间序列时，数据量会非常大。确保你的计算机有足够的硬盘空间（至少是原始数据量的10-20倍）和内存（32GB以上会更顺畅）。处理中间文件尤其占空间。

3. 干涉处理的核心三步：配准、生成干涉图与滤波

数据准备好了，就进入核心的干涉处理阶段。这一步的目标是生成高质量的干涉相位图。

3.1 精密配准：让像素“对齐”

两景SAR图像即使覆盖同一区域，由于卫星轨道、姿态的微小差异，同一个地面目标在图像上的像素位置也会有偏移。配准的目的就是将副图像（Slave）的几何精确对齐到主图像（Master）上，精度要达到亚像素级（通常优于0.1个像素）。

在SARscape中，这通常在/InSAR工具箱下的/Coregistration模块中完成。流程一般是：

1. 粗配准：基于卫星轨道参数，计算一个初步的偏移多项式。
2. 精配准：在整幅图像上选取大量的、分布均匀的连接点（Tie Points），通过窗口匹配算法（如幅度互相关）精确计算每个连接点的偏移量。SARscape会自动完成这一步。
3. 重采样：根据计算出的精确偏移模型，对副图像进行重采样，使其网格与主图像完全一致。

这里的关键经验和坑点：

• 连接点数量与质量：软件默认的连接点数量有时不够，尤其在山区、水体或低相干区。如果配准精度报告显示残差很大（RMS误差 > 0.2像素），你需要手动增加连接点数量或调整连接点搜索策略。配准精度差是后续所有问题的根源，会导致干涉条纹混乱甚至无法解缠。
• 主图像选择：通常选择时间序列中间时刻的影像作为主图像，可以使所有像对相对于主图像的空间基线最小化。SARscape的堆栈分析工具可以帮助你自动选择最优主图像。

3.2 干涉图生成与去平地效应

配准完成后，就可以生成干涉图了。在/Interferogram Formation模块中，软件会：

1. 将配准好的主副图像复数相乘：Interferogram = Master * conj(Slave)。结果是一个复数图像，其相位值就是两景图像之间的相位差。
2. 去除平地相位：这是非常关键的一步。即使是在完全平坦的地面上，由于雷达的侧视几何，也会产生随距离向线性变化的相位条纹，这被称为“平地效应”。如果不去除，它会掩盖真正的地形和形变相位。SARscape会利用卫星轨道和DEM信息模拟并减去这个相位。

生成的干涉图，其相位值（-π 到 +π）包含了三种主要成分：

干涉相位 = 地形相位 + 形变相位 + 大气相位 + 噪声

我们的目标是通过后续步骤，一步步分离出“形变相位”。

3.3 自适应滤波与相干性计算

原始的干涉图噪声很大，主要是由时间去相干、空间去相干和热噪声引起的。直接在这样的图上进行相位解缠几乎不可能成功。因此，必须进行滤波。

SARscape提供了多种滤波器，最常用的是Goldstein自适应滤波器。它的原理是根据局部相干系数的大小动态调整滤波强度：在相干性高的区域（如城市）弱滤波以保留细节；在相干性低的区域（如植被）强滤波以抑制噪声。

滤波强度的选择（Alpha值）是个经验活。默认值（如0.5）是个不错的起点。如果滤波后条纹变得模糊不清，可以适当降低Alpha值（如0.3）；如果噪声仍然很多，可以增大（如0.7）。一个重要的检查方法是：对比滤波前后的干涉图和相干图。好的滤波应该在平滑噪声的同时，尽可能保持清晰的干涉条纹轮廓。

与此同时，软件会计算相干系数图。这个值范围在0到1之间，表示两景图像对应像素的相似程度。1代表完全相干（如稳定的建筑物），0代表完全失相干（如流动的水体）。相干图是评估干涉质量、指导后续解缠和掩膜生成的黄金标准。

4. 相位优化：轨道精炼、重去平与相位解缠

经过滤波，我们得到了一个“干净”一些的干涉图，但里面仍然混叠着残余的轨道误差相位、地形相位和非线性形变。这一步就是要进一步提纯。

4.1 轨道精炼与重去平

即使使用了精密轨道，卫星轨道的微小误差仍会引入一个在空间上呈线性或低阶多项式变化的相位趋势，这被称为“轨道残差相位”。它会表现为干涉图中大范围的、弯曲的条纹背景。

/Orbit Refinement and Re-flattening模块就是用来解决这个问题的。它的原理是：

1. 在干涉图中，选择一些高相干、且假定没有形变的区域作为“控制点”。通常，这些点会选择在远离形变区的稳定地面（如基岩出露区、老城区）。
2. 假设这些控制点上的相位残余主要由轨道误差和可能的大气效应构成，拟合一个多项式曲面（如一阶或二阶）来模拟这个残余相位场。
3. 从整个干涉图中减去这个拟合的曲面。

这一步的成败关键在于控制点的选取：

• 自动选择：SARscape可以基于相干系数自动选择点。你需要设置一个较高的相干阈值（如 > 0.7）。
• 手动检查与编辑：永远不要完全相信自动选择！务必打开控制点叠加在干涉图上的图层，逐一检查。删除那些落在明显形变区（如沉降漏斗）、水体或相位跳变边缘的点。误选的控制点会扭曲拟合曲面，反而引入误差。
• 多项式阶数：通常从一阶开始尝试。如果减去后，干涉图背景仍有明显的弯曲，可以尝试二阶。但阶数越高，过拟合的风险越大，可能会把真实的缓慢形变信号也当误差去掉。

4.2 相位解缠：从缠绕相位到连续形变

这是D-InSAR流程中技术难度最高、也最容易出问题的一步。干涉相位被包裹在（-π, π]区间内，呈现周期性的条纹。而真实的地形和形变是连续的。相位解缠就是要把这些“缠绕”的相位恢复成连续的绝对相位。

SARscape主要提供Minimum Cost Flow (MCF)和Region Growing等算法。MCF是目前最主流和稳健的算法之一。

解缠的核心参数与实战技巧：

1. 解缠掩膜：这是最重要的输入。你需要基于相干图生成一个掩膜，告诉解缠算法只在哪些像素上进行解缠。通常，设置一个相干阈值（如0.3），低于此值的区域（水体、茂密植被）被掩膜掉。在掩膜内部，解缠更容易成功。
2. 初始解缠点（Seed Point）：你需要手动或自动指定一个起点，这个点的解缠相位被设为零（或一个已知值）。这个点必须位于一个高相干、相位质量极好的位置。通常选在图像中心附近一个稳定的、无条纹的亮目标上（如角反射器或永久散射体）。
3. 残差点处理：干涉图中存在一些相位循环不闭合的点，称为残差点。它们可能由噪声、低相干或相位突变引起。MCF算法通过设置“枝切线”来阻止解缠路径穿过残差点。你需要设置一个残差点密度阈值，密度过高的区域（通常是低相干区）可能被自动排除。
4. 多尺度解缠：对于大区域或复杂条纹，可以启用多尺度解缠。它先在低分辨率图像上完成全局解缠，再将结果作为约束引导高分辨率解缠，提高成功率和速度。

如何判断解缠是否成功？

• 视觉检查：解缠后的相位图应该是相对平滑、连续的表面，不应该有突然的、不合理的剧烈跳变（除了真正的形变边界，如断层）。
• 查看解缠报告：SARscape会生成报告，显示解缠成功率、残差点数量等。成功率一般应高于85%-90%。
• 一致性检查：如果处理了多个像对，可以观察不同像对在同一稳定区域的解缠相位是否一致（应接近零或一个常数）。

如果解缠失败（出现“条纹”状的错误），你需要回溯：检查配准精度、提高滤波强度、调整解缠掩膜（缩小范围或提高相干阈值）、或更换初始解缠点。

5. 形变反演、地理编码与成果解读

经过千辛万苦的解缠，我们得到了连续的相位差。最后一步，就是把这相位差转化为有物理意义的形变量，并放到地图坐标系里。

5.1 从相位到形变：反演计算

在/Phase to Displacement或/Inversion模块中，SARscape会进行如下计算：

1. 去除地形相位残余：即使使用了DEM去平，由于DEM误差或配准误差，解缠相位中可能仍有与地形相关的残余相位。这一步可以利用外部高精度DEM或通过模型进行进一步改正（可选，但对山区很重要）。
2. 相位到形变的转换：对于D-InSAR（两景图），形变（沿雷达视线方向，LOS）的计算公式非常简单：LOS_Displacement = (解缠相位 / (4π)) * 波长哨兵1号的波长（C波段）约为5.6厘米。因此，一个完整的2π相位循环对应着约2.8厘米的形变（因为波往返一次，相位变化是双程的）。正相位值通常表示目标远离卫星（沉降），负值表示靠近卫星（抬升），但这个约定取决于软件处理流程，务必查看软件说明。

5.2 地理编码：从斜距坐标系到地图坐标系

SAR图像原本处于斜距-地距坐标系。为了与其它GIS数据叠加分析，必须进行地理编码，将其转换到标准的地图投影坐标系（如WGS84 UTM）。

在/Geocoding模块中，你需要：

• 指定输出投影和像素大小。
• 输入用于地理编码的DEM（通常与去平时使用的DEM一致）。

地理编码后，你得到的就是最终的产品：一张地理坐标的形变栅格图，单位通常是米或厘米。

5.3 成果验证、可视化与常见问题排查

拿到形变图，工作只完成了一半。解读和验证同样重要。

• 量级合理性检查：根据你研究区域的地质背景和已知信息，判断形变量级是否合理。例如，年沉降速率通常从几毫米到几十厘米不等，如果出现米级的异常值，很可能是解缠错误或大气影响。
• 与已知点对比：如果有水准测量、GPS等地面实测数据，一定要进行对比验证。这是评估D-InSAR精度的最直接方法。
• 大气相位校正：对于单对D-InSAR，大气延迟（尤其是水汽）是主要的误差源，其空间尺度可能很大，容易与大规模缓慢形变混淆。SARscape提供了基于气象模型或滤波的方法来估计和去除大气相位，但对于单一时相，校正效果有限。更可靠的方法是使用时间序列InSAR（如PSI, SBAS）技术。
• 可视化技巧：在GIS软件中，使用合适的色带（如彩虹色或发散色带）来显示形变，并设置科学的拉伸范围以突出细节。叠加道路、建筑、行政区划等矢量图层，可以帮助定位形变中心。

典型问题与排查思路：

• 整个图像出现条带状周期性误差：很可能是配准不精，返回检查配准报告的RMS误差。
• 形变图边缘出现巨大异常值：通常是解缠在低相干区域失败导致的相位“爆炸”。用相干图生成掩膜，在解缠和反演时排除这些区域。
• 形变图案与地理特征完全不相关：检查地理编码是否正确，DEM是否有问题。也可能是轨道精炼时选择了错误的控制点。
• 结果看起来“一片模糊”，没有细节：可能是滤波过度，或者多视数设置过大，损失了空间分辨率。

最后，记住D-InSAR测量的是沿雷达视线方向的一维形变。要获取真实的三维形变，至少需要结合升轨和降轨两个方向的数据进行分解，或者与其它测量手段融合。对于哨兵1号这样的C波段数据，它对垂直方向的形变最敏感，但对东西向的水平形变几乎不敏感。