3. 矿区SAR数据获取与预处理判读要点

InSAR形变监测的精度上限,往往在数据获取与预处理阶段就已经被决定。矿区地表覆盖着植被、裸露采坑、排土场和运输道路,地形起伏剧烈,且开采活动使地表在短时间内发生快速变化。如果SAR数据参数选择不当、成像几何配置不合理,或者配准与多视处理存在偏差,后续的相位解缠和形变判读就会出现系统性误差,甚至把预处理引入的伪形变误判为真实沉降。因此,本节把"数据获取—几何配置—配准多视—相干性影响—质量检查"作为一条完整链条来讲,重点说明每一步应该看什么、怎么判断。

一、SAR数据参数选择:波段、分辨率与重访周期

矿区监测首先要回答"选什么数据"。波段决定了穿透能力与对植被的敏感程度:C波段(如Sentinel-1,波长约5.6 cm)对植被冠层散射敏感,在植被覆盖区相干性下降明显,但数据免费、重访周期12天,适合大范围普查;L波段(波长约23 cm)穿透植被能力更强,在矿区林地与农作物区能保持较好相干性,但空间分辨率通常较低;X波段(如TerraSAR-X)分辨率高(可达1 m级),适合采坑边坡、工业广场等重点目标的精细监测,但波长最短,对植被和大气扰动最敏感。实际选型时,应把"监测对象尺度"与"地表覆盖类型"放在一起权衡,而不是单纯追求高分辨率。

分辨率与幅宽是一对矛盾参数。条带模式(Stripmap)幅宽窄但分辨率高,适合单个采区;扫描模式(ScanSAR/TOPS)幅宽可达250 km以上,适合矿区群整体普查,但方位向分辨率下降。重访周期则直接决定能否捕捉快速形变:矿区开采沉陷在活跃期月沉降量可达数十厘米,若重访周期过长,相邻两景间的形变梯度过大,会导致干涉相位严重失相干甚至解缠失败。教学示例中,若某矿区活跃采面预计月沉降约30 cm,采用12天重访的Sentinel-1,单干涉对形变约12 cm,在C波段半波长2.8 cm条件下相位梯度已相当大,需要结合短基线策略和相位解缠质量控制来处理——这类量级估算应在选数据阶段就完成。

二、轨道与成像几何:基线、入射角与视线方向

InSAR只能测量沿雷达视线(LOS)方向的一维形变,因此成像几何直接决定形变投影关系。垂直基线决定干涉对的空间相干性:基线过长会使几何去相干加剧,一般应控制垂直基线在临界基线的十分之一以内;时间基线则受地表变化和植被生长影响,矿区应尽量选用时间间隔短的干涉对。入射角影响形变敏感度,LOS形变与垂直形变的关系为 d_LOS = d_v·cosθ(忽略水平分量),入射角越大,对垂直沉降的敏感度越低,但侧视几何对边坡的水平位移更敏感。

轨道方向同样关键。升轨与降轨的LOS方向不同,同一沉降在两种几何下投影符号和量值都不同,联合升降轨可以分解出垂直与东西向形变分量。矿区判读时应明确记录每景数据的轨道方向、飞行方位角与入射角,并在成果中统一说明"正LOS表示朝向卫星运动(抬升),负LOS表示远离卫星(沉降)"这一符号约定,避免后续解译时方向混淆。

矿区SAR成像几何与LOS形变投影关系 SAR卫星 入射角 θ 平地 排土场 采坑边坡 目标点 P 垂直沉降 d_v LOS 投影 d_LOS = d_v·cosθ 升降轨几何差异 升轨 LOS 降轨 LOS 干涉基线约束 B⊥ 短 B⊥ 长 B⊥ 宜 < B_c/10

三、配准与多视处理:把像素对齐到同一地表

配准是干涉处理的基础。主从影像若存在亚像素级错位,干涉相位中会引入与地形和错位量相关的斜坡条纹,严重时整幅干涉图无法使用。工程上通常采用外部DEM辅助的地理编码配准或基于强度互相关的精配准,配准精度应达到0.1~0.2个像元。判读时可直接观察干涉条纹:若条纹呈规则平行斜坡且与地形高度相关,往往提示配准残差或DEM误差,而非真实形变。

多视处理通过距离向和方位向平均降低相位噪声、提高辐射分辨率,但会牺牲空间分辨率。矿区采坑边坡、裂缝等小尺度目标对分辨率敏感,多视比不宜过大;而大范围沉陷盆地可采用较大多视比以压制噪声。教学示例中,Sentinel-1 IW数据距离向像元约2.3 m、方位向约14 m,若采用4×1多视,方位向仍保留较高分辨率,适合线性工程与边坡;若采用1×4多视则方位向被过度平滑,可能漏掉窄小裂缝。多视后应记录最终地面像元尺寸,作为后续形变梯度判读和滤波窗口选择的依据。

四、矿区植被与地形对相干性的影响

相干性是判断干涉质量最直接的指标。矿区相干性损失主要来自三方面:一是植被去相干,农作物、灌木在时间基线内生长或随风摆动,C波段下尤为明显,夏季农田区相干性可能降至0.3以下;二是地形去相干,陡坡区叠掩与阴影导致有效像元缺失,采坑深部边坡常出现阴影区;三是时间与形变去相干,活跃采面在短时间内形变梯度过大,相位梯度超过一个条纹周期即无法可靠解缠。

判读时应把相干图与光学影像、DEM坡度图叠加分析:低相干区若与林地、农田、陡坡阴影高度吻合,属于可解释的系统性失相干,应在成果中圈定为"低置信区";若低相干区出现在平坦裸地且与开采活动时间吻合,则更可能是快速形变引起的失相干,反而提示该处形变剧烈,需要缩短时间基线或改用L波段数据复核。这一区分是矿区判读中容易被忽略却非常关键的一步。

五、预处理质量检查清单

预处理完成后,建议按以下清单逐项核查,任一项不通过都应回溯处理,而不是直接进入解缠与形变解算。

检查项判读依据合格参考
影像对时间/空间基线基线表与轨道文件垂直基线小于临界基线1/10,时间基线尽量短
配准精度干涉条纹形态、配准残差图残差≤0.2像元,无规则斜坡条纹
多视与像元尺寸处理参数记录与监测目标尺度匹配,记录最终地面分辨率
相干性分布相干图叠加光学影像与坡度低相干区可解释,圈定低置信区
相位连续性缠绕相位图条纹连续、无突跳断裂
几何与符号约定轨道方向、入射角记录LOS正负定义统一并写入成果说明
实用提示:把每景数据的轨道方向、入射角、垂直基线、多视比和最终像元尺寸记录在同一张元数据表中,与形变成果一并归档。后续出现"形变方向异常"或"量值偏小"的争议时,这张表往往是快速定位问题的第一手依据。

总体而言,矿区SAR预处理的核心不是把流程跑通,而是让每一步的参数与判读依据都可追溯、可解释。只有在数据参数、成像几何、配准多视和相干性分析都通过质量检查之后,后续的相位解缠与形变解译才具备可靠基础。