1. 为什么我最终选择了ENVISARscape这套组合1.1 从一次地面沉降监测项目说起去年接手了一个城区地面沉降监测的活儿面积大概覆盖了主城区加周边工业园区约200平方公里。甲方要求提供近三年的形变速率图精度要能识别出每年5毫米以上的沉降趋势。一开始我尝试用开源方案SNAP做预处理、StaMPS做时序分析折腾了将近两周光是InSAR数据的配准和相位解缠就反复出问题尤其是StaMPS的编译环境和依赖库版本冲突换了两台机器才勉强跑通结果出来的点密度还不太理想。后来一个做地质灾害监测的老同行跟我说你试试ENVI加SARscape虽然要装的东西不少但流程化程度高PS和SBAS两条线都有现成的向导式操作参数调优的空间也够用。我抱着试一试的心态装了ENVI 5.6配上SARscape 5.6从数据导入到出形变结果第一个实验区跑通只用了三天。这个效率差距让我决定把整套流程固定下来后面几个项目都沿用这套方案。这篇文章就是把我从环境搭建到PS-InSAR和SBAS-InSAR两条技术路线跑通的完整过程整理出来包括中间踩过的坑、参数怎么调、结果怎么判断靠不靠谱。如果你也在做地表形变监测不管是矿区沉降、城市基础设施监测还是滑坡隐患识别这套流程应该都能直接拿来用。1.2 PS和SBAS到底该选哪个很多人一开始会纠结PS-InSAR和SBAS-InSAR到底用哪个我的经验是这两个不是二选一的关系而是互补的。PS-InSARPersistent Scatterer InSAR永久散射体干涉测量的核心思路是找那些在时间序列上相位和幅度都保持稳定的点——比如建筑物的角反射器、裸露的岩石、人工设施等。这些点在多期影像中回波特性变化很小所以可以拿来精确反演形变。它的优势是精度高单个点的形变测量能到毫米级缺点是在植被覆盖区、农田、水面这些地方几乎找不到稳定的PS点点密度会很低。SBAS-InSARSmall Baseline Subset InSAR小基线集干涉测量走的是另一条路。它不要求单点稳定而是通过组合短基线、短时间间隔的干涉对降低去相干的影响然后通过奇异值分解把多个干涉对联合起来求解形变时间序列。它的优势是在自然地表、植被区也能得到不错的点密度缺点是精度相对PS略低而且对干涉对的选择比较敏感。我通常的做法是城市建成区以PS为主城乡结合部和植被区用SBAS补充。如果项目经费和时间允许两条线都跑一遍交叉验证结果会更可靠。对比维度PS-InSARSBAS-InSAR适用地表城市、岩石、人工设施植被区、农田、自然地表点密度建成区高植被区极低整体较均匀形变精度毫米级亚厘米级到毫米级所需影像数一般≥20景一般≥15景计算量较大中等对DEM精度要求高中等1.3 哨兵一号数据为什么适合这个组合哨兵一号Sentinel-1是C波段SAR卫星波长约5.6厘米重访周期12天单星IW模式下的幅宽250公里分辨率5米×20米。这个参数组合对于区域尺度的形变监测来说非常合适——幅宽够大一景就能覆盖一个中等城市分辨率够用城市建筑和基础设施的PS点能识别出来重访周期短时间采样密度足够捕捉季节性形变。更重要的是哨兵一号的数据是免费开放的通过欧空局的Copernicus开放访问中心就能下载。我一般用SLC格式的IW数据配合精密轨道文件POD Precise Orbit和AUX_POEORB辅助数据。下载的时候注意选对极化方式VV极化对地表形变监测最常用如果是城区也可以考虑VH辅助判断地表类型。提示下载哨兵一号数据时精密轨道文件通常比影像晚几天发布建议影像下载后等一周左右再下载对应的轨道文件避免用快速轨道Restituted Orbit导致基线误差。2. 环境搭建ENVI 5.6与SARscape的安装细节2.1 安装前的准备工作ENVI 5.6和SARscape的安装不算复杂但有几个地方特别容易出问题我一个个说。首先是系统环境。ENVI 5.6官方支持Windows 10 64位专业版及以上我实测Windows 11也能跑但要注意.NET Framework版本。SARscape 5.6需要.NET Framework 4.6.2以上如果系统里装的是更老的版本安装程序会直接报错退出。建议安装前先在启用或关闭Windows功能里确认.NET版本或者直接去微软官网下载最新的4.8版本装上。其次是安装路径。ENVI和SARscape的安装路径里绝对不要有中文和空格。我见过有人装在Program Files下面路径里有空格结果SARscape的某些命令行工具调用时路径解析出错。建议直接装在D盘或E盘根目录下比如D:\ENVI56和D:\SARscape56简单干净。第三是许可证。ENVI 5.6用的是FlexNet授权SARscape是单独的许可证文件。安装顺序是先装ENVI再装SARscape最后配置许可证。SARscape的许可证需要绑定网卡MAC地址如果你换了网卡或者用了虚拟网卡许可证会失效。我一般建议在物理机上装虚拟机里跑虽然也能用但计算性能会打折扣尤其是SBAS的奇异值分解那一步虚拟机里跑大区域数据会非常慢。2.2 安装步骤与常见报错处理安装过程本身是向导式的一路下一步就行但有几个关键节点需要注意。ENVI 5.6安装运行安装程序后选择安装目录等待文件复制完成。安装完成后不要急着启动先安装许可证管理器。许可证文件一般由供应商提供放到指定目录后在许可证管理器中指向该文件即可。SARscape 5.6安装SARscape的安装程序会检测ENVI的安装路径如果检测不到会提示手动指定。这里要注意SARscape的版本必须和ENVI版本匹配ENVI 5.6对应SARscape 5.6不能混用。安装完成后SARscape会自动在ENVI的扩展模块目录下注册。常见的报错及处理方式Cannot find ENVI installationENVI没有正确安装或者注册表信息缺失。解决办法是重新运行ENVI安装程序选择修复或者手动在SARscape安装向导中指定ENVI路径。License checkout failed许可证问题。检查许可证文件是否过期、MAC地址是否匹配、系统时间是否正确。系统时间不对也会导致许可证校验失败这个坑我踩过。IDL runtime errorSARscape依赖IDL运行时环境如果IDL版本不匹配会报这个错。ENVI 5.6自带IDL 8.8SARscape 5.6需要IDL 8.8确认版本一致即可。Cannot open temporary file临时目录权限问题。SARscape在运行时会往系统临时目录写大量中间文件如果临时目录空间不足或权限不够会报错。建议把临时目录设置到一个空间充足的盘至少留100GB以上。注意SARscape处理哨兵一号数据时中间文件体积可能是原始数据的5到10倍。一个200平方公里的项目如果用了30景影像中间文件可能占到200GB以上。提前规划好磁盘空间别跑到一半发现盘满了。2.3 哨兵一号数据的导入与预处理数据导入是SARscape流程的第一步。SARscape支持直接导入哨兵一号的SLC数据操作路径是SARscape Import SAR Spaceborne Sentinel-1。导入时需要提供以下文件哨兵一号的SLC数据文件.zip或解压后的.SAFE文件夹精密轨道文件AUX_POEORB可选的AUX_CAL校准文件导入过程中SARscape会自动完成几件事读取元数据、提取轨道信息、计算基线、生成SARscape格式的影像文件。导入完成后你会得到一个包含幅度图和相位图的文件对以及一个参数文件。这里有个细节值得注意哨兵一号的SLC数据是复数格式包含幅度和相位两部分。幅度图用来做PS点的初选相位图用来做干涉。导入时SARscape会自动把这两部分分开存储后续处理时根据需要调用。导入完成后建议先做一次多视处理Multi-looking。多视的目的是降低斑点噪声、提高信噪比代价是空间分辨率降低。对于形变监测一般方位向和距离向各做4:1或2:2的多视具体取决于你的监测目标尺度。如果是城市建筑级别的监测多视比可以小一些如果是区域性的地面沉降多视比可以大一些。3. PS-InSAR流程的完整拆解3.1 主影像选择与干涉对组合PS-InSAR的第一步是选主影像。主影像的选择直接影响后续所有干涉对的质量所以这一步不能随便。SARscape提供了自动选主影像的功能依据是时间基线、空间基线和多普勒质心频率的综合评分。评分越高的影像作为主影像与其他影像的相干性越好。我一般会让SARscape自动选然后人工检查一下选出来的主影像是不是在时间序列的中间位置。如果自动选出来的主影像偏前或偏后可以手动调整。选完主影像后SARscape会自动生成所有从影像与主影像的干涉对。对于PS-InSAR所有从影像都和主影像做干涉形成一个星形拓扑结构。这也是PS-InSAR和SBAS最大的区别之一——PS是单主影像SBAS是多主影像。干涉对生成后需要检查一下基线分布。理想情况下空间基线应该控制在临界基线的1/3以内时间基线尽量短。哨兵一号的临界基线大约在5公里左右实际处理中空间基线控制在150米以内比较稳妥。如果某些影像的基线太大可以考虑剔除避免引入过多的去相干。3.2 差分干涉与相位解缠干涉对生成后下一步是差分干涉。差分干涉的目的是去掉地形相位得到形变相位。SARscape的差分干涉流程包括生成干涉图把主从影像的相位做差得到干涉相位。去平地效应去掉参考椭球面引起的相位趋势。去地形相位用外部DEM模拟地形相位并减去。滤波降低干涉图的噪声提高相干性。这里DEM的选择很关键。SARscape默认使用SRTM 1弧秒DEM对于大多数区域够用。但如果你的研究区地形起伏很大建议用更高精度的DEM比如ALOS World 3D 30米DEM或者本地实测DEM。DEM的精度直接影响差分干涉的质量DEM误差会直接混入形变相位。相位解缠是PS-InSAR里最考验经验的一步。SARscape提供了多种解缠方法我常用的是最小费用流法Minimum Cost Flow和枝切法Branch Cut。最小费用流法对噪声的鲁棒性更好适合相干性一般的区域枝切法速度快适合相干性好的区域。解缠参数里有一个解缠相干性阈值默认是0.2。这个值设得太低噪声区域也会被解缠导致解缠错误传播设得太高低相干区域会被掩掉点密度下降。我的经验是城区设0.3到0.4植被区设0.15到0.25根据实际相干性图调整。3.3 PS点候选与形变反演解缠完成后SARscape会进行PS点的候选提取。PS点的判断依据主要是幅度离散指数Amplitude Dispersion IndexADI和相干性。ADI越小说明该点在时间序列上幅度越稳定越可能是PS点。SARscape里PS点候选的阈值一般设ADI小于0.25相干性大于0.7。这两个阈值可以根据实际情况调整。如果点密度太低可以适当放宽如果噪声点太多可以收紧。PS点候选出来后进行形变反演。反演的核心是建立相位与形变之间的函数模型然后通过最小二乘或奇异值分解求解。SARscape会自动完成这一步输出每个PS点的形变速率和累计形变量。反演完成后一定要做残差检查。残差是观测相位与模型预测相位之间的差异残差越小说明模型拟合越好。如果残差普遍偏大可能是解缠出了问题或者形变模型选得不对比如实际形变是非线性的但用了线性模型。3.4 结果导出与精度验证PS-InSAR的结果可以导出为多种格式我一般导出为CSV和Shapefile两种。CSV方便做统计分析Shapefile方便在GIS里叠加其他图层。精度验证是最后一步也是很多人容易忽略的一步。验证方法主要有两种与GNSS数据对比如果有研究区的GNSS连续观测站可以把PS点的形变速率和GNSS解算的形变速率做对比。两者差异在2毫米/年以内说明结果可靠。交叉验证用SBAS的结果和PS的结果做交叉验证如果两条线在重叠区域的形变趋势一致说明结果可信。我一般两种方法都用。GNSS对比看绝对精度SBAS交叉验证看点密度和空间一致性。4. SBAS-InSAR流程的关键差异点4.1 干涉对组合策略的不同SBAS和PS最大的区别在于干涉对的组合方式。PS是单主影像星形拓扑SBAS是多主影像短基线集。SBAS的干涉对选择遵循两个原则空间基线尽量短、时间基线尽量短。SARscape里可以设置空间基线和时间基线的阈值自动生成满足条件的干涉对。我一般设空间基线小于150米时间基线小于60天。这样生成的干涉对数量适中既能保证相干性又不会因为干涉对太少导致解算不稳定。干涉对的数量对SBAS的解算影响很大。干涉对太少奇异值分解的秩不够形变时间序列解不出来干涉对太多计算量爆炸而且冗余的干涉对可能引入误差。经验值是干涉对数量在影像数量的2到3倍之间比较合适。比如30景影像干涉对控制在60到90对。4.2 相位解缠与轨道精炼SBAS的相位解缠和PS类似但因为干涉对更多解缠的工作量更大。SARscape支持批量解缠但建议先做一次轨道精炼Orbital Refinement把轨道误差去掉再解缠。轨道精炼的原理是利用地面控制点GCP来估计和校正轨道误差。GCP的选择很关键要选那些形变已知且稳定的点比如远离形变区的基岩露头。SARscape里可以手动选GCP也可以用自动选点功能。我一般手动选10到20个GCP分布在研究区四周确保轨道误差估计的稳定性。轨道精炼完成后再做一次去残余相位把大气延迟和轨道残差进一步去掉。这一步对SBAS的结果影响很大尤其是大区域项目大气延迟的空间变化可能达到几厘米不去掉的话形变结果会被大气信号淹没。4.3 奇异值分解与时间序列解算SBAS的核心是奇异值分解SVD。因为SBAS的干涉对是冗余的观测方程是超定的需要用SVD来求解形变时间序列。SARscape里SVD的参数主要有两个奇异值截断阈值和时间相干性阈值。奇异值截断阈值控制保留多少个奇异值阈值太小会保留噪声太大可能丢掉有用的信号。我一般设0.01到0.05之间根据奇异值谱的拐点来定。时间相干性阈值控制哪些点参与解算一般设0.6到0.7。SVD解算完成后SARscape会输出每个点的形变时间序列。这个时间序列是SBAS最有价值的输出之一因为它能反映形变的时间演化过程而不仅仅是平均速率。比如你可以看到某个区域是持续沉降还是季节性波动这对判断形变机制很有帮助。4.4 PS与SBAS结果的融合与对比PS和SBAS的结果各有优劣融合两者的结果能得到更全面的形变场。融合的方法有几种简单叠加、加权平均、贝叶斯融合。我一般用加权平均权重根据点的时间相干性和幅度离散指数来定。PS点的时间相干性高、ADI小权重给大一些SBAS点的时间相干性稍低权重给小一些。融合后的结果可以在GIS里可视化用颜色表示形变速率用箭头表示形变方向。如果研究区有已知的形变特征比如矿区的沉降漏斗、滑坡体的滑动方向可以叠加对比验证结果的合理性。对比项PS-InSAR结果SBAS-InSAR结果融合结果点密度城区高中高点密度植被区低中高中高形变精度高中高高时间序列完整性好好好计算时间长中长5. 实操中踩过的坑与参数调优经验5.1 基线估计的常见错误基线估计是InSAR处理的基础基线估错了后面全错。我踩过的坑主要有两个第一个坑是轨道文件用错。哨兵一号的轨道文件有快速轨道Restituted Orbit和精密轨道Precise Orbit两种。快速轨道在影像获取后几小时就能下载但精度差基线误差可能达到几厘米精密轨道要等两周左右但精度高基线误差在毫米级。我一开始图快用了快速轨道结果干涉图上的条纹明显不对后来换成精密轨道才正常。第二个坑是基线阈值设得太宽。SARscape默认的空间基线阈值是临界基线的50%对于哨兵一号大约是2500米。这个阈值太宽了实际处理中空间基线超过150米相干性就明显下降。我后来把阈值收紧到150米干涉图质量明显提升。5.2 解缠错误的识别与修复解缠错误是InSAR处理中最常见的问题表现为形变结果中出现不合理的跳变或大面积异常值。识别解缠错误的方法看解缠后的相位图如果出现明显的条纹断裂或相位跳变说明解缠出了问题。另外看残差图如果残差在某个区域突然增大也可能是解缠错误导致的。修复解缠错误的方法调整解缠相干性阈值把阈值调高把低相干区域掩掉避免错误传播。增加滤波强度在解缠前加强滤波提高干涉图信噪比。手动编辑SARscape支持手动编辑解缠结果对于局部错误可以手动修正。换解缠方法最小费用流法不行就换枝切法或者反过来。我的经验是解缠错误80%以上可以通过调整相干性阈值和滤波参数解决剩下的20%需要手动编辑或者换方法。5.3 大气延迟的去除技巧大气延迟是InSAR形变监测的主要误差源之一尤其是C波段的哨兵一号大气延迟可能达到几厘米足以淹没毫米级的形变信号。去除大气延迟的方法主要有三种经验模型法用气象数据温度、气压、湿度建立大气延迟模型然后从干涉相位中减去。SARscape支持导入外部气象数据但精度有限。相位堆叠法假设大气延迟在时间上是随机的通过多期干涉图的堆叠平均来削弱大气影响。这个方法简单但对长时间序列的形变监测效果有限。GACOS数据法GACOSGeneric Atmospheric Correction Online Service提供全球的大气延迟校正数据精度比经验模型高。SARscape支持导入GACOS数据我一般都会用。提示GACOS数据是免费下载的但需要提前申请账号。下载后是栅格格式需要转换成SARscape能识别的格式。转换时注意坐标系和分辨率要匹配。5.4 计算资源与效率优化SARscape处理大区域、多影像的InSAR数据时计算资源是瓶颈。我总结了几条优化经验分块处理把大区域切成几个小块分别处理后再拼接。SARscape支持分块处理可以显著降低内存占用。并行计算SARscape支持多线程在参数设置里把线程数调到CPU核心数。但注意不要超过物理核心数超线程反而会降低效率。SSD硬盘中间文件的读写非常频繁用SSD比机械硬盘快很多。我实测同样的数据SSD上跑比机械硬盘快40%以上。内存分配SARscape默认的内存分配可能不够可以在配置文件里调大。一般建议分配系统内存的70%到80%。6. 形变结果的解读与报告输出6.1 形变速率图的判读要点拿到形变速率图后怎么判断结果靠不靠谱我一般看几个方面第一看整体趋势。形变速率图应该呈现出合理的空间分布特征。比如城区沉降一般集中在几个漏斗状区域滑坡体形变一般沿滑动方向呈条带状分布。如果形变速率图看起来像随机噪声那结果肯定有问题。第二看异常值。形变速率图里如果有孤立的、极端的异常值比如±100毫米/年大概率是解缠错误或相位跳变导致的需要检查原始干涉图。第三看与已知信息的吻合度。如果研究区有已知的地面沉降区、采矿活动区、地下水超采区形变结果应该和这些已知信息吻合。如果不吻合要么是数据处理有问题要么是形变机制比预想的复杂。6.2 时间序列曲线的分析时间序列曲线是SBAS的强项也是判断形变机制的重要依据。我一般会挑几个典型点把它们的形变时间序列画出来看看形变是线性的还是非线性的。线性形变时间序列近似一条直线说明形变速率稳定可能是地下水持续超采或矿区持续开采导致的。季节性波动时间序列呈现周期性起伏说明形变受季节性因素影响比如冻融循环、季节性灌溉。突变时间序列在某个月份突然跳变说明可能有突发事件比如地震、滑坡启动。时间序列分析需要结合研究区的实际情况不能只看曲线形状。比如同样是季节性波动在北方可能是冻融在南方可能是降雨入渗。6.3 报告输出的规范与建议形变监测报告一般包括以下几个部分项目概况研究区范围、数据源、时间跨度、处理方法。技术路线PS和SBAS的流程简述关键参数设置。结果展示形变速率图、时间序列曲线、典型区域放大图。精度验证与GNSS对比、交叉验证结果。结论与建议形变趋势判断、风险区域识别、后续监测建议。报告里的图件要注意坐标系和比例尺确保读者能准确对应到实际地理位置。形变速率图的颜色分级要合理一般用红-黄-蓝渐变红色表示沉降蓝色表示抬升中间用白色或浅色过渡。提示报告里最好附上误差分析说明形变结果的精度水平和不确定性范围。这不仅是专业性的体现也是让甲方信任结果的关键。7. 从项目实战中积累的几点个人体会7.1 数据质量决定结果上限InSAR形变监测的结果好坏七分靠数据三分靠处理。如果原始影像的相干性差、轨道误差大、大气延迟严重再好的算法也救不回来。所以我在项目开始前一定会先做一次数据质量评估检查影像的相干性、基线分布、轨道精度如果数据质量不达标宁可换数据源或者调整时间窗口也不硬着头皮处理。7.2 参数没有万能值SARscape里的参数很多但没有一套参数能适用于所有场景。同样是PS-InSAR城区的参数和植被区的参数完全不同同样是SBAS小区域和大区域的参数也不一样。我的经验是先跑一遍默认参数看结果然后根据结果调整参数迭代两到三次基本就能找到合适的参数组合。7.3 交叉验证是信心的来源单一方法的结果总让人心里没底PS和SBAS交叉验证再加上GNSS对比三方面的结果一致才能放心把报告交出去。我做过一个项目PS和SBAS的结果在大部分区域一致但在一个局部区域差异很大后来查出来是那个区域有施工活动地表变化剧烈PS点失相干SBAS的点也受到了影响。这种细节不做交叉验证根本发现不了。7.4 持续学习是必须的InSAR技术发展很快新的算法、新的数据源、新的处理工具不断涌现。SARscape的版本也在更新每个版本都会增加新功能、优化旧流程。我一般会关注SARscape的官方更新日志看看有没有值得尝试的新功能。另外多和同行交流也很重要很多实操中的小技巧文档里不会写只有做过的人才知道。最后分享一个我常用的检查清单每次跑完InSAR流程后按这个清单过一遍能避免大部分低级错误轨道文件是否用了精密轨道空间基线是否控制在150米以内解缠相干性阈值是否根据区域调整残差是否在合理范围内形变速率图是否有不合理异常值时间序列曲线是否符合预期是否做了GNSS对比或交叉验证报告里的坐标系和比例尺是否正确这套流程我跑了十几个项目从城市沉降到矿区监测到滑坡识别基本都能覆盖。ENVI加SARscape的组合虽然不是免费的但流程化程度高、参数可控、结果可靠对于需要出正式报告的项目来说效率和质量都比纯开源方案省心不少。如果你刚开始接触InSAR建议先从一个小区域、少量影像跑通流程熟悉了之后再上大项目循序渐进踩的坑会少很多。
