哨兵一号数据处理实战：ENVI+SARscape的InSAR完整流程详解

刚接触哨兵一号数据的朋友，最容易踩的坑就是：兴冲冲把数据下载回来，扔进ENVI里一打开，发现屏幕上是一片灰蒙蒙的东西，跟平时处理的光学影像完全不是一回事。你甚至看不到建筑物和道路，更别提做什么地物分类了。问题不在于数据坏了，而是在于你拿处理光学影像的思路去处理合成孔径雷达数据，从底层逻辑就错了。

哨兵一号（Sentinel-1）是欧洲空间局哥白尼计划里的雷达卫星，全天时全天候成像，能穿透云雨，这些特性让它在滑坡监测、地面沉降、洪涝评估、海洋船舶检测这些场景里非常能打。但想把这些数据真正用起来，单靠ENVI本身那套光学处理工具是不够的，得挂上SARscape这个扩展模块，让两者协同工作。ENVI负责影像基础管理和可视化，SARscape负责专业的SAR处理链路，从数据导入、配准、滤波、干涉到形变提取，一步到位。

这篇文章我把整套流程拆开讲清楚：环境怎么搭、数据怎么导、参数怎么定、报错怎么排、结果怎么用。不是泛泛而谈，是按我实际跑项目的习惯来写的，读者不管是刚入门还是已经跑了几个流程的老手，都应该能从里面拿走点有用的东西。

1. 哨兵一号的数据特性，决定了你不能拿它当光学影像处理

1.1 一张图看清SAR数据与光学数据的差别

哨兵一号是C波段SAR卫星，工作频率约5.405GHz，对应波长5.5厘米左右。它自带主动微波传感器，自己发射电磁波脉冲并接收地面反射回来的回波信号，不需要依赖太阳光照。这个特性带来两个直接结果：白天黑夜都能成像（所以叫全天时），云层和雨雪对它基本透明（所以叫全天候）。

但代价是，SAR影像的物理含义跟光学影像完全不一样。光学影像记录的是地物对太阳光的反射率，人眼看着直观；SAR影像记录的是地物的后向散射强度，灰度值取决于地物的介电常数、粗糙度、几何结构，还有雷达波的入射角。比如平静的水面在SAR影像里是暗的（镜面反射把能量弹走了），而城市建成区因为角反射效应特别亮。这些规律需要专门的知识去解释，不是靠眼睛看就能看明白的。

SAR还有两个光学影像没有的麻烦。第一个是斑点噪声，因为雷达回波是相干叠加，同一块均匀地物在不同分辨单元里也会产生随机的亮暗起伏，看起来像撒了一层盐和胡椒；第二个是几何畸变，侧视成像导致叠掩、透视收缩、阴影等现象，在山区尤其明显，不做几何校正就叠加分析，结果会错得离谱。

1.2 ENVI和SARscape各自扮演什么角色

ENVI本身是个通用遥感平台，影像打开、波段组合、光谱分析、面向对象分类这些它都很强。但涉及SAR的专有处理，比如多视处理、干涉图生成、相位解缠、轨道精炼，ENVI原生功能覆盖不到。

SARscape就是补这个缺口的。它作为ENVI的扩展模块存在，菜单栏里会多出一整套Sarscape工具，本质上是一个第三方开发者针对SAR数据处理开发的插件包，把Synthetic Aperture Radar领域常用的算法都集成进了ENVI的图形界面里。你不需要写代码去调Gamma或者GMTSAR，鼠标点几下就能跑通一套InSAR流程。

实际使用中两者的分工大概是这样的：

关键词是“协同”——SARscape跑计算，ENVI做外围支撑，两者互相配合才能把流程走得顺。

2. 装环境是最容易翻车的一步：版本、顺序与license

2.1 版本匹配关系

我见过太多人卡在第一步：SARscape装上了，菜单也出来了，但一运行就报错，或者干脆找不到模块。十有八九是版本对不上。

SARscape不是独立安装就完事的，它对ENVI的版本有严格依赖。不光主版本要匹配，64位和32位也不能错，现在基本都是64位了。常见的匹配关系大概是这样的（以我实际用过的为准，不同时期可能微调）：

不是说你拿SARscape 5.6往ENVI 5.3里装就一定跑不了，但我建议别挑战这种兼容性问题。SARscape安装程序在安装过程中要识别ENVI的安装目录，版本不匹配时经常会出现识别不到或者初始化失败的状况，浪费时间不说，还容易让人误判是数据出了问题。

2.2 安装顺序与配置要点

正确顺序是：先装ENVI，再装SARscape，最后装license文件。严格按这个顺序来。

安装SARscape时它会自动探测ENVI的安装路径，然后往ENVI的扩展目录里写入模块文件。装完之后启动ENVI，正常情况下菜单栏里会出现“Sarscape”这一项。如果没有，去ENVI的Extensions目录下看一眼，SARscape相关文件夹在不在，绝大多数情况都是没装全或者被安全软件拦了。

License这块要单独说。SARscape需要专门的许可文件，许可文件里注明了模块权限，比如基本的“SARscape Basic”和进阶的“SARscape Interferometry”是分开授权的。你要是只装了Basic授权，菜单里的干涉测量相关工具就是灰色不可点击的。所以安装完成后，第一件事就是确认授权模块覆盖了你需要的功能范围。

提示：装完SARscape之后先不要急着导数据，建议重启一次ENVI，让软件重新扫描扩展模块。这一步能省掉很多莫名其妙的初始化问题。

2.3 安装完成后必须做的三件事

第一，检查菜单栏是否出现Sarscape，并随意点开一个工具，确认能正常弹出参数面板。

第二，点开SARscape的Preferences（偏好设置），确认临时目录和缓存路径设置合理。默认路径可能在C盘系统盘，如果数据量大，C盘空间不够就会在跑干涉时报错。

第三，测试一下GPU加速。SARscape部分模块支持GPU并行计算，如果你的机器有独立显卡，在偏好设置里把GPU选项打开，能明显加快滤波和干涉的计算速度。但注意，不是所有模块都吃GPU，有时候开了GPU反而报错，那就先关掉，用CPU跑，稳定性优先。

3. 数据源与导入：拿到SAFE包之后该怎么喂给SARscape

3.1 哨兵一号有哪几种模式和数据级别

哨兵一号的工作模式有四种：条带模式（SM）、干涉宽幅模式（IW）、超宽幅模式（EW）和波浪模式（WV）。其中IW是陆地区域默认观测模式，也是做地面形变监测最常用的，幅宽250公里，分辨率5米乘20米。EW主要用在极地冰盖和海洋区域，WV用于海洋动力参数反演。

数据级别上，哨兵一号官方提供两大类产品：SLC（单视复数）和GRD（地距多视）。

做干涉测量、形变反演这类要做相位分析的，必须用SLC数据，因为相位信息就在SLC里。GRD是把相位信息抹掉了的多视强度产品，做后向散射分析、洪水制图、船舶检测可以用，但做不了干涉。

所以下载数据之前先想清楚目的。大多数人做沉降监测、滑坡排查，别犹豫，直接下IW模式的SLC数据，双极化（VV+VH）通常都下，反正存储也不贵，用不上就扔着。

3.2 导入环节最容易被忽略的参数

SARscape对哨兵一号的数据导入路径是：Sarscape → Import → Import Data → SENTINEL-1 → Sentinel-1 SLC。

选择好SAFE文件夹后，系统会自动识别里面的清单文件。这里有几个参数需要仔细确认：

极化方式要选对。IW模式一般包含VV和VH两个极化通道。做干涉分析通常用同极化VV，因为相干性更好，尤其在植被覆盖区，VV的干涉质量比VH高。导入时你可以只导VV，也可以两个都导，但后面的干涉流程里要用同一个极化的数据。

输出路径不要放在中文或带空格的目录下。这在后面单独说，但导入阶段就要养成习惯。

还有一个容易被忽略的是“数据输出格式”。SARscape导入后会生成自己的.si文件和配套的.sml参数文件，这些是后续处理的基础，文件夹不要手动去改名字或移动位置，SARscape是按照路径管理的。

3.3 导入后的快速检查

导入完成后，用ENVI打开生成的.sml文件，肉眼看一下影像范围对不对、后向散射强度分布是否合理。正常陆地场景应该是城市亮、水体暗、山体有阴影起伏。如果整幅图都是暗的或者出现大面积条纹，可能是导入时选择的极化方式不对，或者原始数据有问题。

同时检查一下每个像元的大小。SLC数据的像元尺寸是2.3米乘14米左右（多视处理后不一样），SARscape在导入时会自动写入这些元数据，你在文件头（Header）信息里能看到完整的描述。如果像元大小显示异常，后面的地理编码肯定偏。

导入这个环节是整套流程里最简单但最关键的一步。数据没导对，后面所有处理都是白费力气。我见过有人把GRD当SLC导，程序跑一半报出相位信息缺失的错误，回头检查才发现数据级别选错了。

4. 干涉测量全流程：连接图、滤波、解缠与地理编码

4.1 主影像选择与连接图生成

干涉测量的原理，是把同一地区两段时间不同、几何位置略有差异的雷达影像组合起来，计算它们之间的相位差。这个相位差里包含地形相位、形变相位、大气延迟相位和噪声，后续处理的本质就是一步步把这些成分拆开。

首先要用Sarscape → Interferometry → Interferogram Generation → Connection Graph工具做干涉像对的连接图。核心工作有两块：选主影像和确定像对组合。

主影像的选择标准主要是时间基线、空间基线和多普勒质心频率差。SARscape的连接的图工具会自动评估这些参数，并给出一个建议的最佳主影像（它会自动选一个能让所有像对相干性整体最好的时间点）。在没有特殊需求的情况下，直接用自动推荐的主影像就行。

但如果你有特定目的，比如做某一灾区的应急形变分析，那就要围绕灾害发生前后这个时间窗口来选像对，不要机械地接受默认推荐。

空间基线是个硬指标。正交基线越大，干涉相位对地表微小形变的敏感性越高，但相干性会下降。SARscape里用临界基线来过滤像对，超过临界基线的组合会被标记为不可用。对哨兵一号来说，它在设计上就是一个小基线星座，重访周期6天（两颗星轮流），基线普遍很小，所以绝大多数像对都能用，这也是哨兵一号天生适合做InSAR的原因。

4.2 干涉图生成和滤波参数怎么定

连接图确认后，下一步是跑干涉图生成。这一步内部做的事情很多：根据精密轨道数据对像对做配准、重采样，计算复干涉图，做去平地效应，还可以顺手生成相干性图。

配准时SARscape默认的算法精度很高，普通场景直接认默认参数就行。但在地形起伏大的区域，如果配准质量不好，后面干涉条纹会非常破碎。这时候要把配准里的DEM选项打开，用外部DEM辅助配准，能显著提升山区像对的配准质量。

干涉图生成后会有一堆输出：复干涉图、强度图、相干性图、去平后的干涉图等。先看相干性图，它反映了两个时相地物散射特征的相似程度。相干性高（0.4以上）的区域，形变信号可信度高；相干性低（低于0.2）的区域，基本就是噪声主导，解缠都很困难。

滤波是为了抑制斑点噪声对干涉图的干扰。SARscape里常用的是Goldstein滤波，参数Alpha值控制滤波强度，默认值是0.4到0.5之间。取值越大，平滑越强，条纹越连续，但也容易把细节抹掉。城区这种相干性高的地方，Alpha取0.3左右能保留更多细节；植被区相干性差，取0.5或更高才能让条纹看起来像样。

4.3 相位解缠的选参策略

相位解缠是整个InSAR流程里最烧脑也最容易出问题的环节。原始干涉相位是被包裹在-π到π之间的，要恢复真实的相位变化，就得进行解缠。SARscape里的解缠算法主要是基于最小费用流（MCF）的，它会在低相干性区域自动设置“障碍”，让解缠路径绕开不可靠的区域。

解缠的关键参数有两个：解缠方法和相干性阈值。

解缠方法一般选Region Growing（区域生长）或Minimum Cost Flow。区域生长速度快，适合相干性整体较高的区域；最小费用流稳定性好，适合场景复杂、相干性分布不均衡的地区。我个人的习惯是，数据质量好就用区域生长，植被多、山区多就用MCF。

相干性阈值默认是0.2，意思是相干性低于这个值的像元不参与解缠，相当于被掩膜掉。阈值设得太低，噪声会被纳入解缠，结果出现明显的“跳变”现象——相位图上一块一块的不连续；阈值设得太高，有效解缠区域太少，很多地方是空洞。这个参数没有绝对的黄金值，建议跑完看一眼结果，如果跳变严重就把阈值往上调，0.25到0.35之间都可以试。

4.4 轨道精炼、重去平与地理编码

解缠做完不等于结束，还有最关键的一步叫轨道精炼（Refinement），目的是去除因为轨道误差引起的残余相位趋势。SARscape的做法是让你在残余相位图或形变图上选地面控制点（GCP），用这些点的已知信息来拟合轨道误差并扣除。

选GCP有几个原则：不能选在有残余条纹的地方（说明相位还没解干净）、不能选在形变区（否则会把真实形变当误差去掉）、不能选在低相干区域。最理想的控制点是稳定、相干性高、无明显形变的裸露地表或人工建筑。实际操作中我会在控制点表里加至少十几个点，尽量均匀分布在全图范围，而不是挤在一块。

轨道精炼完成后，如果用的是单像对流程，SARscape会生成一个“Phase to Displacement”的步骤，把残余相位转换成雷达视线向（LOS）的形变量。最后用地理编码工具把结果投影到标准地理坐标系（通常选UTM投影，和你的DEM投影一致），输出成GeoTIFF，就能在ENVI里直接叠加到光学影像或地形图上做分析了。

5. 高频报错排查：process unexpectedly terminated和它的亲戚们

5.1 "SARscape process unexpectedly terminated"究竟是哪一环出了问题

这句话应该是运行SARscape的人最常见的噩梦。SARscape是个在ENVI进程里跑的插件，一旦底层算法计算出问题或者资源不足，整个进程就会崩掉，弹窗里只有这句干巴巴的提示，没有任何具体错误码。

我总结下来，最常引发这个报错的原因有三个。

一是数据路径问题。路径里有中文、有空格、有特殊符号，SARscape底层算法对这些支持极差，一旦读不到文件就崩溃。解决方案很粗暴：所有数据放到纯英文路径下，目录层级尽量浅，比如D:\Data\S1\，下面建个子文件夹存数据就完了。

二是磁盘空间不够。干涉处理中间文件非常大，尤其是一组时序影像同时处理时，中间产物动辄几十GB。SARscape在写临时文件时发现磁盘满了，不会优雅地提示，而是直接让整个进程死掉。建议预留至少原始数据三倍以上的可用磁盘空间。

三是内存不足。SARscape的很多步骤不支持分块处理，需要在内存里加载整个景或大块数据。32GB内存的机器跑标准景勉强够，但做时序分析一次性载入几十景就很容易爆内存。这时候没什么好办法，加内存，或者减少同时处理的影像数量。

5.2 内存和磁盘的账，得提前算

这里我给出一个粗略的估算方法：一景哨兵一号IW模式的SLC数据，导入后大约3到5GB，干涉处理产生的中间文件（配准后的SLC、滤波后的干涉图、解缠结果）加起来大约是原始数据的5到8倍。如果你做12景的时序处理，光中间文件就可能超过200GB。

所以开始处理一个项目前，先做三件事：确认磁盘剩余空间是预估数据量的2倍以上，确认内存足够大（16GB是下限，32GB起步舒服，64GB可以流畅处理时序）、确认系统临时目录所在盘符空间充足——SARscape的临时文件默认放在系统Temp目录，不在你的数据盘上。

5.3 轨道文件下载失败的退路

干涉处理需要用到精确轨道信息。哨兵一号之前提供了两套轨道产品：辅助轨道（AUX_POEORB，精度高但晚了几天）和预报轨道（AUX_RESORB，出炉快但精度差一截）。SARscape会自动尝试下载精确轨道文件，但经常因为网络连接不稳定而失败。

处理办法很简单。预测到要跑干涉前，先手动把对应时段的轨道文件下载好，放到SARscape可识别的目录里。SARscape的偏好设置里可以指定轨道文件的存放目录，下载好的文件放进去，程序运行时会优先读本地文件。

如果手动下载的轨道文件也报错，多数是格式不匹配或者文件名不规范。命名要按照官方规则来，不要随便改文件名。

6. 典型应用：城区地表沉降的时序形变分析

6.1 从一组影像到沉降速率图的完整路线

哨兵一号时序形变分析，最典型的方向是城市地面沉降监测。过去做沉降监测要布设水准点或GPS连续站，成本高、点密度稀。用InSAR技术，一发卫星影像就能提供上百万个测点，密度和覆盖范围完全不在一个量级。

我自己跑过的案例是华东某城市区，用了12景IW模式SLC数据，时间跨度8个月，跑的是小基线集（SBAS）路线，SARscape里有独立的SBAS模块，要分几步走：先做连接图和干涉像对生成，再做轨道精炼和重去平，然后做SBAS反演，最后做地理编码。

这里要特别提醒的是，SARscape的SBAS流程对像对数量有要求。不是说两景影像就能跑SBAS，至少需要一组影像（通常8景以上），像对组合越多，反演的形变序列越平滑可靠。如果你只有三五景数据，老老实实做单像对干涉就行，别强行上时序。

SBAS反演完成后，输出的沉降速率图通常叠加到光学影像或高精度地形底图上，在ENVI里再做个简单的阈值分类，把沉降速率超过-10毫米/年的区域标红，就是一张非常直观的灾害风险图。

6.2 精度验证的实操做法

用InSAR测出来的形变速率，不能直接当成真值。常用的验证方法有两种。

第一种是跟外部测量数据比。如果研究区里有水准点、GPS站或者IGS站，就把它们坐标处的InSAR视线向形变值提取出来，和目标水准数据进行回归分析。理论上，如果在理想条件下（极小的大气误差、完美的解缠），两者的R²应该在0.8以上。实际能做到0.7左右，就算比较可信了。

第二种是看时序形变的内部一致性。同一地区的形变速率不管用哪几个像对组合来算，结果都应该大致吻合。如果在同一时间窗口下，不同像对组合算出来的形变速率差超过3毫米/年，说明处理质量有问题，需要回到滤波或解缠环节排查。

6.3 进一步拓展方向

哨兵一号时序数据不只是吃沉降一个方向。SARscape的成果可以拿来接后续的很多分析：

比如在滑坡隐患早期识别里，用InSAR测出来的微小形变速率识别潜在滑动体，再结合光学影像判断边界，比单纯靠肉眼解释卫星影像高效得多。

再比如城市基础设施监测，地铁沿线、高架桥、大坝这类人工建筑的形变热点，都能从时序形变成果里观测到。配合PS-InSAR技术（SARscape也有专门模块），还能提取到单个建筑物级别的稳定散射体形变信息。

海冰监测和洪水快速制图，用GRD数据的后向散射变化就行，处理流程比干涉简单得多，但也能和干涉结果形成互补。

从我个人的使用体会来说，ENVI加SARscape这套组合最大的优点是门槛低、路径清晰。缺点也很明显——流程中间过程的参数选择确实需要经验积累，很多参数没有绝对的对错，只能根据场景反复试。

最后分享一个实用小技巧：处理大批量时序数据的时候，把每一步的中间结果都保存好，不要急着覆盖。因为某个参数调了之后效果反而变差的情况太常见了，能随时回到上一步重新来，比从头再跑一遍省太多时间。