搞InSAR的人应该都体会过市面上的干涉合成孔径雷达处理软件要么是付费商业套件要么是学术免费但安装配置极其折磨人的开源工具。而ENVI配合SARscape这套组合几乎是行业里默认的“正规军”配置从哨兵一号Sentinel-1数据导入、配准、干涉图生成、滤波、解缠到形变图输出全流程在一个图形界面里就能完成不需要写一行代码。如果你正要上手哨兵一号数据的干涉处理或者已经在处理过程中被各种报错折磨到怀疑人生这篇文章就是按我自己的实操经验整理的一份完整流程记录从软件环境、数据准备到每一步参数设置再到高频报错排查一次性讲透。先说清楚这篇文章适合谁。刚接触InSAR、手里有哨兵一号数据但不知道怎么下手的新手可以直接按流程抄作业已经跑通过一遍但总在某个步骤卡住的老手重点看第五部分的报错排查和参数心得应该能帮你省不少时间。无论是做地表形变监测、滑坡隐患识别还是矿区沉降分析这套流程都适用差异只在于研究区的DEM选择和相干性阈值调整。1. ENVI与SARscape这对组合为什么能成为行业默认1.1 这对工具组合解决了什么问题先理清一个概念ENVI是图像处理平台SARscape是运行在ENVI上的一个扩展模块。没有SARscape的ENVI也能打开栅格数据、做分类、做光谱分析但处理SAR数据时尤其是InSAR干涉测量这一整套流程就力不从心了。SARscape是Sarmap公司开发的商业化SAR处理工具以ENVI插件的形式存在启动后在ENVI的工具栏里会多出SARscape这个专属菜单里面从数据导入、滤波、配准、干涉、解缠到地理编码全流程都有现成模块。选择这套组合而不是别的方案最大的优势是省心。比如GAMMA或ISCE这类开源工具处理结果确实更可控但环境配置、数据格式转换、命令行操作每一项都在消耗精力。SARscape把所有步骤封装成了图形界面参数设置得相对合理默认值就能跑出不错的结果。对大多数做工程应用的人来说能用图形界面解决的问题没必要去折腾命令行。再一个优势是生态联动。ENVI本身的光学影像处理能力可以和SARscape形成互补SAR提取的形变结果可以直接叠加在光学底图上做分析ENVI的纹理特征提取、面向对象分类、土地利用转移矩阵等功能可以和SAR的后向散射强度图结合做更深度的应用。这也是我后面单独花篇幅写延展应用的原因。1.2 哨兵一号数据本身有什么特殊性哨兵一号是欧洲航天局哥白尼计划的雷达卫星星座由A星和B星组成工作在C波段波长约5.6厘米重访周期在双星运行时可以做到6天。对InSAR来说C波段的优势是相干性比X波段好对植被区域的穿透能力也更强劣势是波长较长对微小形变的敏感度不如X波段但在大多数地表形变监测场景下完全够用。哨兵一号的成像模式主要有三种条带模式SM、干涉宽幅模式IW和超宽幅模式EW。做InSAR最常用的是IW模式它采用TOPSTerrain Observation with Progressive Scans成像技术幅宽约250公里地面分辨率5米乘20米。TOPS模式有个特点影像在方位向和距离向都存在显著的几何畸变数据导入和配准的精度要求比传统条带模式更高这也是为什么SARscape专门为哨兵一号做了TOPSAR专用的导入和处理流程。很多新手跳过这个专用流程直接按普通SAR数据导出来的干涉图质量就会很差。2. 环境准备ENVI版本、SARscape安装与License问题2.1 版本匹配关系不容忽视ENVI和SARscape之间存在严格的版本对应关系这是很多人安装时踩的第一个坑。SARscape的版本号通常与对应的ENVI版本号一致比如SARscape 5.6对应ENVI 5.6SARscape 5.7对应ENVI 5.7。版本不匹配的后果很直接SARscape菜单加载不出来或者处理过程中直接崩溃。具体的匹配关系建议以SARscape官方发布说明为准一般来说最新的SARscape都支持对应大版本的ENVI同时向下兼容某个旧版本。比如5.6版本的SARscape可能支持ENVI 5.3到5.6但安装了5.6的ENVI去配5.5的SARscape就可能出问题。我自己的建议是装ENVI的时候直接选当前SARscape依赖的那个版本不要用太新的ENVI去配旧版SARscape。下载和授权方面ENVI和SARscape都是商业软件正规渠道是通过官方代理商获取安装包和试用授权。SARscape的License文件是单独放置的在ENVI里通过“帮助-许可管理”或者SARscape菜单下的License工具来加载。注意安装路径和工程路径不要出现中文也不要用带空格的文件夹名称否则后续处理很容易出现“SARscape process unexpectedly terminated”这类莫名其妙的中断问题。这不是玄学是SARscape底层对路径解析比较脆弱我在第五部分会细说。2.2 安装SARscape后的环境检查安装完成后先别急着导数据做两个检查第一个是检查菜单完整性。打开ENVI确认主菜单栏出现“SARscape”菜单项点击展开后能看到Import、Filtering、Focus、InSAR等子模块。如果菜单没有出现多半是版本不匹配或者安装时没有以管理员权限运行。第二个是检查系统资源。SARscape处理哨兵一号数据时非常吃内存和磁盘建议内存不低于16GB磁盘预留至少100GB空间。一景IW模式的哨兵一号SLC数据大约4到5GB干涉处理会产生大量中间文件磁盘不够会在处理到一半的时候直接中断。这个坑很多人等到数据处理到一半才遇到那时候再腾硬盘已经晚了。3. 数据准备哨兵一号影像的获取、轨道文件与DEM选择3.1 哨兵一号数据下载渠道与产品选择哨兵一号数据是免费开放的主流下载渠道有欧空局的Copernicus Data Space Ecosystem、美国阿拉斯加卫星站ASF DAAC等。国内用户访问ASF DAAC通常更稳定而且ASF还提供了一个很实用的功能可视化搜索直接在地图上框选研究区设置时间范围、轨道方向、产品类型就能筛选出合适的影像对。下载时要注意产品类型。做InSAR需要的是SLC单视复数数据不是GRD地距影像。GRD是已经做过幅度检波和多视处理的强度影像适合做目视解译和分类不适合做干涉。下载时产品类型选SLCIW模式极化方式根据研究目标选VV或VH都可以但保持一景和参考景的极化一致。选择像对时注意时间基线和空间基线。时间间隔越短地表去相关的影响越小空间基线过长干涉条纹会非常密集解缠难度增大。SARscape的基线估算模块会给出具体数值一般建议空间基线控制在临界基线之下时间基线看研究区的地表类型城区可以放宽植被覆盖区尽量选冬季或时间间隔短的像对。3.2 精密轨道文件的下载逻辑哨兵一号的轨道数据分几种实时轨道POD预测轨道精度较低适合快速查看事后精密轨道POD精密轨道在数据获取后大约21天可以拿到精度在5厘米以内。做正式InSAR处理时强烈建议用精密轨道文件它会把基线误差对干涉相位的影响降到最低。轨道文件的命名和下载很简单去欧空局的数据平台输入影像的绝对轨道号就能找到对应的文件。下载后无需解压直接把压缩包路径告诉SARscape就行。在导入数据或者生成干涉图时软件会自动读取与影像时间匹配的轨道文件。如果处理时刻还没到精密轨道的发布时间可以用存档轨道凑合但精度会受到影响能换精密轨道时建议从头再跑一遍。3.3 DEM的选择SRTM还是外部DEMInSAR流程里需要DEM来做地形相位模拟和地理编码SARscape默认支持自动下载或指定本地DEM。常用的全球DEM有SRTM、ASTER GDEM、Copernicus DEM等。对哨兵一号这类C波段数据来说SRTM 1弧秒约30米分辨率在大部分区域都能满足需求SARscape内置的SRTM 3弧秒版本在数据覆盖面更广但在山区或地形起伏大的区域精度有限。如果研究区的地形复杂建议准备一个30米或更高分辨率的局部DEM通过“External DEM”选项导入。一个小技巧SARscape在导入数据时可以同时完成行列坐标到地理坐标的映射DEM必须覆盖整个研究的影像范围否则后续的干涉图生成会因为缺少地形相位而报错。下载DEM时最好比研究区范围向外多扩一点。3.4 用TOPSAR模块导入哨兵一号数据的操作细节启动SARscape的导入功能路径是SARscape Import SAR Spaceborne Sentinel-1 TOPSAR。这个选项是专门为哨兵一号IW模式的TOPS成像设计的不要选成别的模式。操作步骤并不复杂在“Input Files”里选择哨兵一号的SAFE格式数据文件夹或压缩包。在“Output Files”里设置输出目录建议新建一个项目文件夹把SLC数据、轨道文件、DEM分目录存放。设置极化方式VV或VH单独导入不能一次把多个极化混在一个输出里。“Parameter”选项里可以设置结果输出的类型保持默认SLC即可。点击执行后软件会生成SLC数据文件和对应的元数据文件.sml这个.sml文件是后续处理的基础索引不要随便改名字或移动位置。导入的时间取决于数据量一景数据通常需要几分钟到十几分钟属于正常范围。实操心得导入多个像对时建议把同一轨道的影像按时间先后排列并在文件名里加上日期标记。后期处理时需要频繁对应参考影像和辅影像井然有序的命名能极大降低出错概率。我习惯用“S1A_IW_日期”这种格式之后在SARscape里选影像时一目了然。4. 干涉处理主流程一步一步从像对到形变图哨兵一号的干涉处理在SARscape里通常走这样一条流程基线估算、干涉图生成、干涉图滤波、相位解缠、轨道精炼与重去平、形变图生成、地理编码。每一步的输出都会成为下一步的输入一个环节的参数没设置好后面的结果质量就跟着崩。下面按顺序拆解。4.1 基线估算第一关也是最重要的一关打开SARscape InSAR Baseline Estimation模块分别选择参考影像Reference和辅影像Secondary执行后软件会在图形窗口中显示干涉像对的基线参数包括空间基线长度、时间基线天数、多普勒质心差异等。重点看两个指标空间基线是否在临界基线范围内多普勒质心差异是否显著。对于TOPS模式临界基线大约在500米左右实际处理时建议空间基线不超过临界基线的20%到30%特别是植被覆盖区基线太长会让干涉条纹密集到无法解缠。时间基线的选择没有绝对标准做形变监测时一般控制在一到两个重访周期内也就是6到12天做长时间序列分析时可以根据数据情况放宽到几个月。基线估算通过之后把参考影像和辅影像对应的SLC文件在后续流程中固定使用。SARscape里的基本操作原则是参考影像选时间较早的那景或者选质量更好的那景后续的干涉图、滤波、解缠、地理编码都围绕这个组合展开。4.2 干涉图生成与去平核心逻辑和参数选择干涉图生成的入口是SARscape InSAR Interferometry Interferogram Generation。这一步做三件事影像配准、干涉相位计算、地形相位去除。配准是多步骤流程SARscape会做粗配准和精配准对于TOPS模式还会做burst间的配准。配准精度直接决定干涉条纹质量。在“Parameter”选项里配准窗口大小和搜索半径保持默认即可但要注意如果参考影像和辅影像的轨道差别大可以适当增大搜索半径。干涉图上叠加了地形相位大部分情况下地形相位的贡献远大于形变相位必须通过外部DEM模拟并减掉这就是“去平”或“去地形”过程。在生成干涉图时“DEM”选项选择之前下载的SRTM或外部DEMSARscape会读取DEM进行相位模拟。效果好的干涉图条纹主要出现在形变区域和残余地形误差区域总体看起来干净平滑。参数方面多视比Multi-looking值得说一下。多视处理本质是对方位向和距离向做平均以减少斑点噪声代价是分辨率下降。对哨兵一号IW模式常用的多视比设置在“Range Looks: 4Azimuth Looks: 1”左右实际的空间分辨率大约在20米量级对于区域形变监测来说已经是很好的平衡。如果想保留更多细节可以少多视但滤波时要加一点强度。4.3 干涉图滤波与相干性评估干涉图生成后一般不会直接用原始的干涉相位包含大量斑点噪声解缠前必须滤波。SARscape里常用的滤波方法有两种Boxcar滤波和Goldstein滤波。Goldstein是自适应滤波对高相干区平滑强度大对低相干区保持细节应用更广。滤波的窗口大小选择很关键窗口过大虽然噪声压得厉害但会损失细节甚至把残余条纹给抹掉窗口过小噪声仍然明显。面对一般数据我用默认的5乘5窗口如果干涉条纹看起来特别破碎可以适当增大到7乘7但不要一路加到底。在InSAR流程中滤波的同时通常会生成相干性图它反映了参考影像和辅影像在局部区域的相似程度数值在0到1之间。相干性是衡量干涉质量的核心指标大于0.4的区域通常是可靠的0.2到0.4之间解缠时可能出风险低于0.2的区域基本没法解缠。这部分信息在后面做轨道精炼和解缠时都要用到。4.4 相位解缠从缠绕相位到绝对相位的跳跃解缠是InSAR流程中最容易出问题的环节。干涉图里的相位值实际是缠绕在“-π到π”之间的相对值形变反演需要知道绝对相位解缠就是利用相邻像素之间的相位差来恢复整周模糊度。SARscape的解缠入口是SARscape InSAR Phase Unwrapping。在“Parameter”面板里选择解缠方法。常用的是区域增长法Region Growing和最小费用流法MCF。区域增长速度快适合相干性较高的区域最小费用流对低相干区域更稳健但计算量大一些。解缠时可以把相干性阈值设到0.2到0.25把这个阈值作为掩膜低于阈值的像素不参与解缠。这样可以避免在噪声严重区域生成不可靠的解缠结果。解缠完成后颜色渐变连续的相位图基本算是成功如果出现明显的“跳变斑块”或颜色锯齿可能是解缠失败需要调整相干性阈值或换解缠方法重跑。4.5 轨道精炼与重去平消除残余轨道误差和相位偏置解缠后先别急着生成形变图轨道误差和残余平地相位会以线性趋势的形式混在相位里需要做轨道精炼Refinement。这一步的逻辑是在影像上选一些地面控制点GCP这些点应该位于高相干区域没有残余条纹没有地形条纹不在形变区内。通过拟合这些点的相位估算并扣除轨道误差和相位偏置。SARscape里入口是SARscape InSAR Orbital Refinement软件会在解缠图上叠加相干性图你手动在图像上选点。选GCP是经验活建议遵循以下几点优先选相干性大于0.4的区域。不要在残余干涉条纹密集的地方选点条纹代表相位变化不适合作为“无形变”的控制点。避开明显的地表形变区域比如塌陷区、滑坡体。全局均匀分布不要集中在一个角落。选好GCP后执行软件会输出轨道精炼后的参数文件并对之前的解缠结果做重去平得到更准确的绝对形变相位。这一步不认真做你后面得到的形变图会莫名多出一个沿轨道的线性梯度看起来像整体抬升或沉降实则是轨道误差。4.6 形变图生成与地理编码输出重去平后的相位就可以转换成形变图了入口是SARscape InSAR Phase to Displacement。选择“Displacement”输出SARscape会基于雷达波长哨兵一号为C波段波长约5.6厘米把相位换算成视线向位移。注意这是视线向形变不是垂直向形变如果研究区地形起伏明显需要用入射角做投影转换才能得到垂直向分量。最后一步是地理编码入口在SARscape Geocoding。输入形变图、DEM和轨道参数设置输出的投影坐标系建议和你的底图一致比如UTM或CGCS2000投影执行后得到带地理坐标的形变栅格图。这个结果可以直接在ArcGIS或ENVI里叠加底图做分析了。5. 高频报错排查那些年我踩过的坑5.1 “SARscape process unexpectedly terminated”是最常见的拦路虎上网搜“SARscape process unexpectedly terminated”能跳出一堆求助帖这个报错几乎伴随每个新手必经的阶段。根据我的经验这个报错背后通常对应以下几种原因内存不足。处理大区域、多景哨兵一号数据时SARscape的进程需要的内存很容易超过系统可用内存。比如处理一整景IW模式数据时干涉图生成和解缠阶段占用内存常在8GB以上。如果是用32位环境的旧版本ENVI尤其容易触发。磁盘空间不足。中间文件写入失败时也会以这个报错表达。路径包含中文、特殊字符或空格。这个在前面已经提过。License连接中断或授权过期。输入数据格式不正确比如误把GRD数据当SLC导入或者SAFE数据不完整。排查思路很直接先看日志输出窗口SARscape在报错前会打印最后一步的操作记录比如哪个文件读取失败或者是内存分配失败再看系统资源监控确认内存和磁盘是否够用。如果都看不出问题把任务放到后台运行试试很多前台运行的进程因为在GUI线程和计算线程之间切换出现竞争而崩溃后台运行反而能顺利跑完。注意处理过程中不要随意点击ENVI主界面的其他工具SARscape处理时是单线程占用当前工作区的频繁操作鼠标和切换面板会加大进程被终止的概率。这个习惯改掉之后我的报错率低了很多。5.2 前后台运行模式大任务的一个保命技巧SARscape在参数设置面板里提供了一个选项处理是在前台Foreground还是后台Background执行。后台运行模式下任务提交给独立的进程ENVI界面还能继续做别的操作任务进度在“Job Status”窗口里查看。我的习惯是一景左右的轻量任务前台跑方便实时看到进度两景以上的大数据量任务尤其是做解缠和地理编码时一律后台跑。后台跑的好处不只是抗中断还能避免某些内存管理上的冲突。比如在后台运行时即使你切换到其他应用处理进程的优先级也不会被抢占。5.3 干涉条纹“炸毛”是怎么造成的很多新手在生成干涉图后发现条纹极其混乱噪声像撒了一把盐。这类问题一般出在以下几个方面一是没做配准就直接干涉TOPS数据对配准要求极高不做精配准的干涉图完全没法看二是多视比设置太小噪声没有被平均掉三是滤波窗口太小或没滤波四是DEM与影像范围不完全匹配导致残余地形相位大量混入。排查顺序建议是确认导入时选的是TOPSAR专用流程确认配准质量报告看配准误差是否在合理范围然后看DEM覆盖范围最后再调滤波参数。绝大多数情况下前两个环节稍加修正条纹质量就会有质的提升。5.4 解缠后出现奇怪的“斑块”解缠结果里如果出现大面积的块状跳变通常与相干性低或解缠参数设置有关。可以分几种情况处理低相干区太多解缠时阈值内有效像素太少把相干性阈值从0.2降到0.15让更多像素参与解缠。解缠方法不适合当前数据换成MCF方法试试。轨道精炼时选的GCP落在了形变区导致重去平后相位重新混入误差。如果是最后一个原因需要回到轨道精炼步骤重新选点确保GCP分布避开形变区域。很多人忽略了GCP的重要性实际上选点好坏直接决定了解缠后相位能否正确恢复。6. 处理成果还能怎么用从形变图到更多应用场景6.1 用ENVI提取SAR影像的纹理特征哨兵一号数据做完干涉流程后SLC数据会保留多视强度信息。即使不做InSAR这套数据本身也有很高价值。ENVI里提供了基于灰度共生矩阵GLCM的纹理提取功能路径是ENVI Filter Texture Co-occurrence Measures。可以选择均值、方差、同质性、对比度、二阶矩等纹理特征输出为单独的栅格图层。SAR图像上的纹理特征可以用于地物分类和地表参数反演。比如在城区的建筑区后向散射强度高且纹理对比度大水体区域后向散射弱且纹理平坦植被区域的纹理呈现粗糙特征。把纹理特征影像叠加到光学影像上参与分类能显著提升分类精度尤其适合那些光学影像容易混淆的地类比如裸地和建成区。我试过用哨兵一号的VV极化强度图提取纹理在山地丘陵地区能有效区分梯田、林地与居民地。纹理特征提取时建议在ENVI里设置合理的窗口大小和灰度量化级别窗口3乘3到5乘5比较适合高分辨率数据灰度级别可以设到64级兼顾计算速度和纹理细节。6.2 与土地利用分类和转移矩阵结合SAR的后向散射强度对地表粗糙度和介电特性敏感在土地利用分类中可以作为光学影像的重要补充特别是对光学数据覆盖差的云雨地区。你可以把哨兵一号的多视强度影像进行滤波和地理编码后导入ENVI的监督分类工具比如支持向量机或随机森林中与光学波段一起训练。因为雷达对植被结构、建筑几何、水面平滑度的响应不同分类结果的边界往往比单纯光学分类更清晰。更有意思的是将不同时期的土地利用分类结果叠到一起做转移矩阵分析。比如分别对两期影像做分类用ENVI的Change Detection工具或直接做分类后比较得到土地利用类型转移矩阵可以看到哪块林地变成了建设用地哪块农田变成了水域。把SAR形变图叠加在土地利用转移结果上还能分析出地面沉降是否发生在新建城区这种组合分析对城市规划、环境监测都有很强的现实价值。7. 一个真正能帮到你的小建议把参数记录当成项目标配我想在最后分享一个常年积累的习惯。SARscape处理流程涉及的参数太多从数据导入到地理编码一次全流程跑下来可能设置了几十个参数如果把参数记录下来后面换数据、换区域甚至换同事接手时都能快速复现。我会建立一个文本文件记录每个步骤的关键参数影像对名称、轨道方向、轨道文件路径、DEM名称、滤波窗口大小、相干性阈值、解缠方法、GCP数量、坐标系信息。看起来繁琐但当你处理到第20个像对时会发现这份记录比任何官方文档都有用。另外有个经验是关于GCP选点的我最初总想着用最少的时间尽快把GCP选完后来发现解缠结果不好、形变图上有明显的轨道趋势误差回头检查全是GCP选点不当惹的祸。认真选出30个左右分布均匀、相干性好的点比随便点几百个点有效得多。尤其在地形起伏大的区域选点更要仔细尽量避开山谷阴影和叠掩区域那些地方的相位本身就是无效的。做InSAR本来就是一个“七分数据、三分处理”的活数据质量决定了上限处理细节决定了你能不能接近这个上限。ENVI和SARscape降低了很多上手门槛但每个参数背后的原理还是值得花时间去理解的。希望这份流程记录能让你少走一些弯路把更多精力放在对结果的解读和应用上。
