煤矿沉降 SBAS-InSAR从零看懂陆探时序处理一篇给小白的长文。对照 B 站教程《煤矿沉降 SBAS-InSAR 监测教程基于 LT-1 数据》整理。 感谢B站博主 河马会讲话上一篇讲的是两景差分DInSAR一张图看一场地震。这篇难一档要用二十多景陆探影像组成几十个干涉对反演出半年里煤矿地面一直在沉、沉多快、沉成什么碗。不要求你会遥感。读完你应该能看懂SBAS 和两景差分差在哪、组网为啥要「短空间、长时间」、第一次反演和第二次反演各打什么、以及软件里每一步到底在干什么。目录这篇在讲什么先和上一课对齐DInSAR vs SBAS基础知识比两景差分多出来的那些词时序 InSAR 在回答什么形变模糊度为啥煤矿爱用 L 波段极化HH、VV、交叉极化空间基线再拆一层平行 vs 垂直SBAS 组网小基线子集相干性的两种干涉相干 vs 时序模型相干准相干散射体 / 允许断开块解缠算法MCF vs Delaunay MCF煤矿沉降碗地在往下沉图为啥还是平面的大气的「时间高通、空间低通」完整处理流程① 导入数据② 制作 DEM③ 干涉组网Connection Graph④ 批量干涉处理⑤ 检查干涉图与解缠必要时删对⑥ 第一次反演GCP 精炼 线性速率⑦ 第二次反演去大气⑧ 地理编码⑨ QGIS 制图栅格 矢量点每一步减掉什么、留下什么做到哪算这一步成功和积石山那一课怎么对照术语对照表一句话收束这篇在讲什么教程用的是陆探一号LT-1官方样例山西大同附近煤矿大约2022 年底到 2023 年 5 月半年里二十来景后面报表是21 景L 波段影像。软件是SARscape 6.1和 ENVI 5.6 兼容。数据来自自然资源卫星遥感云服务平台的测试样例。目标不是「震前震后拍一张差」煤矿是持续开采、持续下沉要的是一张年平均或半年沉降速率图毫米 / 年到厘米 / 年量级看值域每个点一条时间序列什么时候开始沉、沉得匀不匀平面上看起来像一口口沉降碗方法是SBAS-InSARSmall Baseline Subset小基线子集。整条链路可以先记成二十多景 SLC → 按「短空间基线」组成几十个干涉对 → 每一对都做差分干涉 滤波 解缠 → 检查、删掉坏对 → 第一次反演去轨道趋势得到粗速率 → 第二次反演按大气的统计特性抠大气 → 地理编码 → 速率图 矢量点和积石山那课相比单对干涉的物理没变还是那五种相位难在对很多、网要连对、还要在时间上把形变和大气拆开。先和上一课对齐DInSAR vs SBAS积石山那课DInSAR这课SBAS影像数量2 景21 景干涉对数量1 对组网后61 对时间跨度4 天同震约半年要的产品一张同震位移图约 7 cm 抬升速率 时间序列煤矿持续下沉组网无核心步骤谁和谁配对轨道趋势一对上做 RPF 控制点每对都可能有斜坡靠 GCP 第一次反演一起扣大气几乎不单独建模4 天太短第二次反演专门去大气形变模型一次事件不拟合速率默认线性速率网密了可用二次项计算量一台电脑也能跑完演示干涉 解缠约占全流程一半时间一句话DInSAR 是「前后两张照片相减」SBAS 是「很多张照片两两精选配对再把所有差分解成一条随时间变化的沉降曲线」。基础知识比两景差分多出来的那些词下面默认你已经知道相位、一个条纹、LOS、失相干、解缠、去平。这些在积石山那篇里写过。这里只把SBAS / 煤矿 / 陆探时序多出来的概念写细。时序 InSAR 在回答什么两景差分回答从 A 日到 B 日地面沿视线一共动了多少时序回答三件更麻烦的事速率平均每天 / 每年沉多少velocity过程是匀速还是某个月突然加快时间序列空间格局哪几个矿在沉、碗有多大、中心在哪要回答它们必须有多余观测。一对干涉只有一个方程分不清「真沉降、轨道斜坡、大气、DEM 误差」。几十对、且每对基线都不大才能把这些东西在最小二乘里拆开。这就是 SBAS 存在的理由。形变模糊度为啥煤矿爱用 L 波段教程一上来就强调半年 20 多景、还是 L 波段对矿区沉降「非常有效」。关键不是「看得更清楚」而是相邻两次成像之间地可以多沉一截而不把条纹缠死。相位每转一圈视线向位移大约是半个波长来回双程哨兵 C 波段陆探 L 波段教程口径公开参数约值波长~5.6 cm21 cm~24 cm半波长一圈~2.8 cm~10 cm~12 cm煤矿沉降可以很快一个月沉十几厘米并不罕见。C 波段下一对干涉就可能转过好几圈山区再一失相干解缠直接崩。L 波段一圈大约 10 cm监测上限大约是哨兵的 45 倍教程原话。所以地震那种「4 天、几厘米」C、L 都能做C 更灵敏煤矿那种「持续、幅度大、地表还可能翻新」L 波段更扛得住教程里把这个叫「成像模糊度 / 形变模糊度」。记住物理就行一圈能代表的位移更大 → 同样的沉降条纹更疏、更好解。极化HH、VV、交叉极化雷达发射和接收可以选水平H或垂直V电场极化符号信号强弱更适合同极化HH、VV强InSAR 常规首选交叉极化HV、VH弱一般不当形变主数据同极化里再分VV对地面微小起伏更敏感海面波浪、草地一类HH对城市、微小构筑物高度变化更合适教程认为更适合沉降监测陆探这套样例只有HH直接用。原则有同极化就用同极化HH / VV 对矿区沉降通常差别不大有什么用什么。空间基线再拆一层平行 vs 垂直积石山那课把「空间基线 741 m」当成一个数。SBAS 组网时要拆开空间基线 ├── 平行基线沿视线 → 主要制造 轨道 / 平地相位残余 └── 垂直基线垂直于视线→ 主要制造 地形残差相位垂直基线大 山陡→ 差分减完 DEM 仍花干涉条纹跟着山走平行基线不准 / 轨道文件差→ 整景斜坡陆探特别常见所以组网不是「基线越小越好」这么简单而是陆探轨道管粗垂直基线容易大 →优先短空间基线L 波段模糊度高、时间失相干慢 →时间基线可以放长。平原矿区教程举安徽反过来地形残差不是主矛盾大气和地表变化更烦更该缩短时间基线。SBAS 组网小基线子集名字里的 Small Baseline 就是不要让每一景都去跟超级主影像配那是老式 PS 思路而是只连接「空间基线够短、时间基线可接受」的像对连成一张网。每个点的沉降从这张网里的所有差分解出来。一张时空基线图怎么读纵轴空间基线米 点 一景影像 横轴时间 线 一对干涉线越短、越平空间基线小山区友好线在时间上拉得很开时间基线长L 波段往往仍能相干某景孤立、没线连出去它不参加解算网从左到右断成两截两段时间各沉各的中间对不上除非开「允许断开块」去插教程这套参数下的报表务必留下印象量数值影像21 景平均空间基线200 m最大空间基线878 m平均时间基线42.5 天每景大约连接6 对向前 3 向后 3总干涉对61经验约等于 景数 × 3对照积石山单对741 m组网之后平均只有 200 m这就是 SBAS「小基线」四个字的意义——不是陆探突然飞准了是你拒绝了那些又肥又远的配对。官方提示每景连接数要大于 5反演才稳。教程习惯设6。冗余度选项含义High凡是落在「时间 ≤ 100 天、空间 ≤ 临界基线的 5%」里的对能连就连数据量大Low按规则瘦身。可按最短时间连前后各 3 景或在 100 天内挑空间基线最短的若干对陆探最大的麻烦是空间基线 → 教程选 Low并在 100 天窗口里优先空间最短。为何不能「每景只与下一景配一对」组网时最容易产生的误解是既然影像已经按时间排好每一景只与下一景做一对干涉1–2、2–3、3–4…再把相邻位移累加不就得到整条沉降曲线了吗这样连成的是一条时间链不是一张网。对 SBAS 反演来说链条在方程数目、误差传播和参数可区分性上都不成立。未知量与观测方程SBAS 最终要估计的是各成像时刻相对于某一参考日的视线向位移或等价的平均速率。(N) 景影像选定一景为参考后待估形变未知量约为(N-1)个。每一对干涉提供一个观测方程含义仅为这两个时刻之间的相位差位移差是多少。若只连接相邻两景干涉对数目 (M \approx N-1)与未知量个数相当。方程组接近恰定能写出一个解但几乎没有多余观测无法做最小二乘意义上的「多数表决」。噪声、解缠跳周、单日大气都会直接写进结果没有第二条路径把它们识别出来。类比测量一条折线的总长若每两段相邻点之间只量一次中间任一把尺子偏了首尾对不上也无从判断是哪一段出错。相干性的两种干涉相干 vs 时序模型相干这是整课最容易混的点。软件里好几处都叫 coherence不是同一个数。干涉相干.cc时序模型相干第一次反演产品何时出现每一对干涉图把所有对反演成速率之后在问什么这两景回波有多像01观测到的形变和你选的模型线性 / 二次贴不贴低意味着植被、叠掩、大基线、噪声非线性沉降、解缠错、大气没去净、模型选错掩膜用途解缠时去掉烂像素地理编码时去掉「速率不可信」的点教程后半段把时序相干阈值调到0.15是怕 0.2 把矿心挖空——矿心沉降快、地表破碎时序相干往往偏低低不代表那里没沉代表那里难拟合。准相干散射体 / 允许断开块传统做法某个像素只要在任何一景里是空值后面所有产品这个位置全废。SARscape 的优势之一Allow disconnected blocks true准相干散射体一类思路。意思是10 景里有 1 景这块失相干不必整条时间序列扔掉。用其余 90% 有值的观测拟合趋势把空档插上。组网时打开它第一次反演里还有一个配套参数用 90% 的有效值去内插 10% 的无效值min valid percent 90。太贪比如 50% 也插会编造沉降。解缠算法MCF vs Delaunay MCFSBAS 要对几十张干涉图解缠算法选错会又慢又错。MCF最小费用流Delaunay MCFDMCF做法整幅都解连低相干、水体也走先按相干性掩膜再在好像素的三角网上解适合相干性好、噪声少相干性差、山区、矿区破碎速度更快更慢、更稳教程建议电脑一般、相干还行电脑好、想稳 → 用 DMCF和积石山那课「调相干阈值再解」是同一件事的工业化版本DMCF 把「先掩再解」写进算法。并行SARscape 5.6.2 以后可同时解多对。经验约1000×1000 像元占 1 GB。30 m 产品大约一景解缠3 GB。32 线程、112 GB 内存的机器教程用16 线程约 48 GB避免把内存打满。叠掩 阴影掩膜务必打开true。陆探基线大叠掩区相位完全混杂阴影区没有回波不解比乱解强。煤矿沉降碗地在往下沉图为啥还是平面的和积石山「向四周扩散」是同一类问题地质机制不同。开采前 开采后 地表 —————— 地表 ——\______/—— ← 沉降碗 岩层 岩层 塌 陷 煤层 ████████ 煤层 采空中心下沉最多往外逐渐回到 0 → 平面图上是一口碗 / 一圈圈等值线以垂直沉降为主边缘会有指向碗心的水平位移地面往中间挤雷达斜视LOS 里会混进这一点水平分量平原矿区仍可近似勾Vertical Displacement山区矿区 陡坡水平分量不能假装没有解释要谨慎所以最终图上「一块红/蓝往外变淡」不是土在平面上往四周流是下沉量从中心向外变小。大气的「时间高通、空间低通」第二次反演不靠气象模型靠一条统计假设教程原话的意思维度大气长什么样滤波时间不是天天下雨今天湿、后天干像随机跳高通去掉时间上跳来跳去的成分空间同一天、同一景附近往往同晴同阴低通去掉单景里空间上很平滑的成分形变煤矿通常是时间上持续、空间上局部成碗。大气则是时间上不连续、空间上成片。两者频谱不同才能滤。风险滤波窗口不合适会把缓慢、成片的真沉降当成大气抠掉或把噪声留着。教程认为山西这块大气不严重而且数据只有半年所以空间窗口不要默认那么狠时间窗口从 365 天改成 180 天数据集就这么长。完整处理流程软件ENVI 5.6 SARscape 6.1。路径只能英文、不能空格、尽量不要小数点空格用下划线代替。源代码不认中文路径。① 导入数据干什么把二十多景陆探文件夹批量读成 SLC。哨兵可以对着 zip 读陆探必须先解压成文件夹。菜单Import Data → Spaceborne Single Sensor → 陆探一号课里技巧不要一景景点进去。用Select file list / 选文件夹指到上一级目录软件递归进子文件夹读meta.xml。极化选HH。一景导入后大约 8 个文件文件去留无后缀或_SLC.sml.hdr必须留二进制 SLC 两套头ENVI SARscape XMLshp kml覆盖范围示意拷贝数据时可不带本质是 ENVI 的.dat栅格外加 SARscape 的 XML。结束时你有什么input_data一类英文目录里21 套 SLC 三件套。② 制作 DEM干什么差分时减地形、地理编码时投坐标都要同一份 SARscape 认识的椭球 DEM。两种来源你手里的 DEM怎么转已是 ENVI 格式ENVI → SARscape选Ellipsoid DEMGeoTIFFOther format → TIFF同样选 Ellipsoid DEM硬性规定输出文件名必须以_dem结尾不要再加.tif这类英文点后缀用 Change output name 改坐标系尽量WGS84不要 CGCS2000、不要 UTM开源 DEM 最好大约30 m教程用 Copernicus DEM干涉多视分辨率应和 DEM 一致。DEM 30 m、影像硬输出 5 m等于用 30 m 高程去内插 5 m地形是假细的结束时你有什么一份xxxx_dem后面组网、干涉、精炼、地理编码反复用。③ 干涉组网Connection Graph干什么决定谁和谁做干涉。这是 SBAS 相对 DInSAR 多出来的第一步也是陆探最要动脑子的一步。菜单SBAS 工作流里选全部 SLC → Connection Graph不要让它瞎自动课里参数山区煤矿 陆探参数课里取值为什么最大时间基线100 天L 波段时间可以放长空间基线阈值临界基线的5%陆探临界大约 2 万米轨道间距又常能收到 1000 m 内先设大一点再靠 Low 瘦身冗余度Low不需要 High 那么多对每景连接数6官方建议 5向前 3 向后 3或 100 天内空间最短的 6 条连接偏好短空间山区地形残差是主矛盾Allow disconnected blockstrue准相干90% 有值也能出点平原矿区请改策略短时间基线优先不要死磕最短空间。看什么时空基线图纵轴跨度尽量短横轴可以拉长报表平均 / 最大空间基线、平均时间基线、总对数主影像日期这套样例落在2023-04-10附近报表会写坐标系陷阱connection graph文件夹里的矢量范围是WGS84 地理坐标SLC 本身是像元坐标左上角为 0方位向一条轴、距离向一条轴没有投影。两者叠不上是正常的不是数据坏了。结束时你有什么一张网 61 对清单。还没有彩虹图。④ 批量干涉处理干什么对网上每一对做配准 → 干涉 → 去平 → 减 DEM → 多视 → 滤波 → 解缠。积石山那课的 ③④⑥ 在这里对 61 对批量发生所以能占掉全流程一半时间。菜单Interferogram generationSBAS 流程里的干涉一步头文件用组网生成的那个DEM 用 ②。分辨率设30 m与 Copernicus DEM 一致。软件自动算多视比。矿碗尺度远大于 30 m要看碗沿更锐的条纹可试 15 m计算量和噪声都会上去。教程为演示用 30 m。Goldstein 滤波陆探相干已经很好不要滤太狠参数课里理由滤波窗口32偏小大窗口会糊积石山单对曾用 64那是山区噪声和解缠在打架Smooth window偏大细节多留一点Alpha自适应强度2.0默认常 2.5滤弱一点相干性估计窗口7默认 5窗口太小会高估相干、把噪声当好像素解缠相干好 → MCF相干差、机器强 →DMCF并行线程按内存估不要打满Apply layover and shadow mask true两个输出目录别看混位置内容工程根目录interferogram stacking快视图 TIFF方便不靠 SARscape 预览work / interferogram stacking真正拿来算的数据.cc.dint.fint.upha等快视图上的大黑边陆探轨道漂移各景覆盖不完全重合软件只留公共区域其余是空值。单对文件含义和积石山同一套文件含义original_dint已去平、已减地形未滤波.fintGoldstein 之后mask_cc叠掩 / 阴影已涂黑.upha解缠结果未滤波就能看清条纹79 天、L 波段相干很好。矿区形变包也还在。花掉的地方多是叠掩、透视收缩 这对空间基线偏大。沿距离向、与方位向平行的宽条纹 轨道相位残差。陆探初始轨道质量较差精密轨道会好一些但时序里仍常要靠后面 GCP 精炼不能指望自动 refinment那是哨兵更常用的。结束时你有什么61 套干涉产品。工程文件夹已经很大这一步完成后再挪盘是允许的。⑤ 检查干涉图与解缠必要时删对干什么SBAS 不怕少一对怕一对解缠错把整网基准带偏。所以不要闭眼进反演。解缠错误长什么样解缠假设整张图从一个零点平滑长出来。错了常见两种局部割裂某条缝左右亮度 / 基准突然不一致跳了整整一截整幅对切一半亮一半暗像被刀切开 —— 参考点在两个区域里对不上矿区在碗里、错误在远处的山可以留。错误落在你要报的矿必须处理。两种修法做法何时改解缠参数重跑这一对多数对都有问题或还能救SBAS Tools → Edit Connection Graph→ 选中这对Remove只有个别对坏教程举例主 1 与辅 3即 2023-01-24 与 2023-02-01 一类删对后不要忘了后面第一次反演若已跑过相关步骤要重跑。结束时你有什么一张「坏对已剔除 / 错误不在矿区」的网。⑥ 第一次反演GCP 精炼 线性速率干什么三件事捆在一起用控制点去掉每对上还在的轨道斜坡陆探几乎必须手动 GCP不要用哨兵那套 auto refinement在时序上按你选的模型默认线性解出粗速率可选估残余高程这套矿区教程关掉这是积石山那课「轨道精炼 重去平」的时序版钉子仍是稳定点但方程变成「很多对一起解」。选 GCP和地震那课同一套纪律更要避开矿碗尽量全局散开不要点在形变区沉降碗里绝对不点不要点在地形剧变处平坦、高相干、相位干净、不要空值三个真值得讲的参数参数课里含义产品相干阈值0.1时序模型相干不是解缠.cc。观测和线性模型的拟合程度形变模型Linear景不多、质量一般时更稳。组网密度到每景 1012 对才考虑Quadratic挖非线性估计残差高false矿区不是城市楼高Copernicus DEM约 2015且 30 m 对齐够用有效比例90%用 90% 有效观测插 10% 空缺解缠相干阈值与干涉一步对齐课里说干涉若为 0.1这里用 0.2 一类以菜单实际为准仍用.cc挡烂像素先看精炼后的干涉图轨道条纹应明显变弱。若某对教程举例 814精炼后斜坡还在回到 ⑤ 删掉再重跑这一步。演示课没抠到极致论文级要更狠。再看三个产品文件看什么velocity粗速率。直方图应接近零均值正态矿区应负值一侧略多以下沉为主precision_velocity越大越不准。用来掩膜coherence时序模型相干直方图若出现双峰 / 多个鼓包多半是解缠基准不一致不是真的有两种地质。用阈值「试砍」决定后面地理编码怎么掩教程现场调的结果请当经验不要当真理阈值现象他怎么改精度 40矿心被盖住太多放到60时序相干 0.2矿心一块掉成空洞放到0.15记住这两个数⑧ 要填进去。对称拉伸看图矿碗应干净右上角若还有缓坡是轨道没去净可再细抠 GCP 或删对。山上和地形相关的花纹若不在关注区可以忍。结束时你有什么粗速率、精度、时序相干轨道大斜坡基本压住。大气还在。⑦ 第二次反演去大气干什么按「时间上跳、空间上成片」的统计假设做时空滤波把大气从粗速率 / 时间序列里抠掉。参数课里怎么动参数默认倾向课里有效比例90不动空间低通窗口米偏小教程说默认约 1200滤得太狠不要滤太强山西大气不重时间高通窗口常 365 天改成180 天数据只有半年执行完看second inversion里的velocity也可叠 LOS 速度。直方图仍应零均值、负值略多正值一侧许多解缠噪声会被压下去。山上仍可能残留大基线造成的解缠错——不在矿区可忽略。结束时你有什么去大气后的速率这是准备地理编码的主产品。⑧ 地理编码干什么斜距产品投到经纬度输出栅格 点矢量。把 ⑥ 试砍得到的阈值填回来参数课里说明时序相干0.15保矿心速率精度60太大则烂点进图太小则矿心被挖高程精度很大SBAS 顺带估的地表高程和沉降不是同一回事放宽以免误伤形变点水体掩膜0这片没有湖有湖再开垂直向平原煤矿、以沉降碗为主勾Vertical Displacement。LOS 仍会混一点水平分量解释时写清楚假设。矢量每个 shp 最多约100 万点600 万点就会拆成 6 个文件。other format选 QGIS 能读的类型。文件夹名称内容未加 VDLOS 结果加了 VD垂直向结果结束时你有什么地图坐标下的速率栅格 一堆点 shp。⑨ QGIS 制图栅格 矢量点栅格速率图若已有单独教程这课主讲矢量点。推荐流程把所有VD的 shp 拖进 QGIS符号化字段选velocity对称、等间隔分类类别和色带按值域对齐色带可用教程分享的干涉彩虹rainbow一类在样式管理器里导入 XML设好一个图层后复制样式 → 其它图层粘贴样式保证所有分文件同一套色点大小用「地图单位 / 米」不要用毫米ArcGIS 常用毫米缩放时点会乱变。QGIS 用米和地面比例尺绑死。目的点大小课里全局图、让矿碗跳出来40 m干涉分辨率 30 m略大于像元局部叠影像、看点和地物套合20 m常见圆点效果15 m 圆形符号底图可叠 Google Earth 一类影像再加标注。每一步减掉什么、留下什么步骤主要减掉 / 压掉主要留下③ 组网空间基线过大、时间过长的「肥对」平均约 200 m 的小基线网④ 干涉地形DEM、椭球平地、噪声滤波、叠掩阴影每对上的形变 轨道残余 大气⑤ 删对解缠崩掉的个别对基准还能对齐的网⑥ 第一次反演轨道斜坡不符合线性模型的乱跳部分粗速率 矿碗⑦ 第二次反演时间上跳、空间上成片的大气更干净的速率 / 序列⑧ 地理编码低时序相干、低精度、可选水体地图上的沉降碗五种相位在时序里的分工平地/轨道 → ④ 去平 ⑥ GCP 精炼 地形 → ④ 减 DEM组网时已经用短垂直基线预防残差 形变 → ⑥⑦ 解出来的速率和序列 大气 → ⑦ 时空滤波 噪声 → ④ 滤波 ⑤ 删对 ⑧ 掩膜做到哪算这一步成功导入21 套 SLC 三件套路径纯英文极化 HHDEM文件名_demWGS8430 m组网平均空间基线明显小于单对乱配每景约 6 连山区纵轴短、横轴可长报表对数约 景数×3干涉未滤波短时间对条纹已清晰矿碗位置有局部条纹平行于方位向的宽条纹允许有后面精炼扣解缠矿区内部连续远处山的割裂可忍整幅对切的对要删第一次反演精炼后斜坡大减velocity直方图单峰、近零均值、负值略多矿碗干净记下相干 0.15、精度 60 这类试砍结果第二次反演正值噪声被压矿碗还在没有把碗心滤没地理编码LOS / VD 分清掩膜后矿心仍在制图分文件同一色带点大小用米碗和影像能套上和积石山那课怎么对照同一套陆探、同一类「轨道差、基线肥、L 波段电离层 / 斜坡」对策从「一对上死磕」变成「网上优胜劣汰」。问题积石山 DInSAR煤矿 SBAS空间基线 741 m 这种肥对只有一对只能硬做组网时尽量不连平均压到 ~200 m整景轨道斜坡RPF 控制点多项式每对 GCP 精炼坏对直接删除滤波窗口曾用 64为了救解缠相干已经很好窗口 32、alpha 调低怕削矿碗解缠阈值0.05因为相干从原始图估、山区要多留批量 DMCF / 叠掩阴影掩膜大气4 天几乎不管第二次反演专门滤垂直向假设逆冲抬升沉降碗平原可勾 VD产品一张位移图速率 点 可导出时间序列组网口诀可以钉在墙上山区陆探短空间、时间可长。平原煤矿短时间更重要。每景至少连 56 对。宁删一对勿留一张整幅切开来的解缠。术语对照表口播 语音识别会错。对照字幕 / 口语正确说法space inside / SPARS / SPARA / sparksSBAS-InSARSASK / SaaS skip / SAAS蒂普SARscapeN 维 / N维他ENVI路探 / 鹿探 / 录探 / 路泰陆探一号 LT-1同级化 / 同计划同极化 HH / VV交叉计划交叉极化成像模糊度 / 形变模糊度半波长对应的形变上限一圈空间极限 / 时间期限空间基线 / 时间基线临界极限临界基线组网冗余度Connection graph redundancy准相干散射体 / allow disconnected block允许断开块部分失相干仍拟合时空基线图时间–空间基线网络图ghost st / ghost down / Go downGoldstein 滤波解馋 / 检查 / snapping相位解缠snaphu / MCF德劳内是 MCFDelaunay MCF蝶眼 / 点点 / 叠眼叠掩 layover透视收缩foreshorteningauto reflementauto refinement 自动轨道精炼first in versionfirst inversion 第一次反演possession / v position / repositionvelocity precision 速率精度残差高residual height 残余高程QU / 二次方模型quadratic 非线性项second step / second inversion第二次反演去大气low pass / high PS空间低通 / 时间高通jo codinggeocoding 地理编码沉降晚沉降碗QG4QGIS阿克西斯ArcGIS比例米地图单位米coo p d e mCopernicus DEMELLIPSODMEllipsoid DEM建议的学习顺序先有积石山那课的底座五种相位、条纹疏密、LOS、解缠、控制点不点形变再搞懂为啥要组网一对方程拆不开速率 / 轨道 / 大气小基线网才拆得开陆探组网策略山区短空间L 波段允许时间拉长每景 ≥ 5 连两种相干.cc管解缠时序相干管「线性模型可信不可信」两次反演第一次钉轨道和粗速率第二次按大气频谱抠天气最后才是 QGIS 点大小那是表达不改计算结果数量级手算和矿碗对得上才踏实陆探一圈约10 cmLOS教程口径若某对只隔 8 天、碗心已经转过两圈说明这 8 天沉了约 20 cm 量级C 波段会更难解半年速率图上若碗心只有每年几毫米而现场矿一直在采先怀疑掩膜、模型、大气滤波把真信号吃了一句话收束这课不是在把积石山流程重复 21 遍而是在讲用半年、二十多景陆探 HH 条带在山区煤矿做 SBAS。因为 L 波段一圈大约 10 cm、时间失相干慢组网可以时间放到 100 天、空间尽量短因为轨道差必须手动 GCP坏对直接从网上剪掉因为要的是持续沉降而不是一次地震必须先线性反演、再按「时间跳、空间成片」去大气最后在沉降碗以垂直为主的假设下投到地图上。SBAS 很多小基线干涉对组成的网 最小二乘时序。组网 主动丢掉肥基线换山区能用的相干。第一次反演 轨道精炼的时序版 粗速率。第二次反演 用统计滤波打大气。沉降碗 采空区上覆岩层塌陷中心沉得多、向外变少——平面图上的圈仍是等位移线不是土在往四周流。对照教程整理。参数以 SARscape 6.1 菜单和这套大同样例为准换平原矿区、换更长时序时组网和时间窗口要改不要把 100 天 / 180 天 / 0.15 / 60 当成全世界通用的魔法数。 
长河君 · 2026/11/6 5:22:59 · 阅读 1.2k 小白看过来，陆探数据嗷——煤矿沉降 SBAS-InSAR：从零看懂陆探时序处理
煤矿沉降 SBAS-InSAR从零看懂陆探时序处理一篇给小白的长文。对照 B 站教程《煤矿沉降 SBAS-InSAR 监测教程基于 LT-1 数据》整理。 感谢B站博主 河马会讲话上一篇讲的是两景差分DInSAR一张图看一场地震。这篇难一档要用二十多景陆探影像组成几十个干涉对反演出半年里煤矿地面一直在沉、沉多快、沉成什么碗。不要求你会遥感。读完你应该能看懂SBAS 和两景差分差在哪、组网为啥要「短空间、长时间」、第一次反演和第二次反演各打什么、以及软件里每一步到底在干什么。目录这篇在讲什么先和上一课对齐DInSAR vs SBAS基础知识比两景差分多出来的那些词时序 InSAR 在回答什么形变模糊度为啥煤矿爱用 L 波段极化HH、VV、交叉极化空间基线再拆一层平行 vs 垂直SBAS 组网小基线子集相干性的两种干涉相干 vs 时序模型相干准相干散射体 / 允许断开块解缠算法MCF vs Delaunay MCF煤矿沉降碗地在往下沉图为啥还是平面的大气的「时间高通、空间低通」完整处理流程① 导入数据② 制作 DEM③ 干涉组网Connection Graph④ 批量干涉处理⑤ 检查干涉图与解缠必要时删对⑥ 第一次反演GCP 精炼 线性速率⑦ 第二次反演去大气⑧ 地理编码⑨ QGIS 制图栅格 矢量点每一步减掉什么、留下什么做到哪算这一步成功和积石山那一课怎么对照术语对照表一句话收束这篇在讲什么教程用的是陆探一号LT-1官方样例山西大同附近煤矿大约2022 年底到 2023 年 5 月半年里二十来景后面报表是21 景L 波段影像。软件是SARscape 6.1和 ENVI 5.6 兼容。数据来自自然资源卫星遥感云服务平台的测试样例。目标不是「震前震后拍一张差」煤矿是持续开采、持续下沉要的是一张年平均或半年沉降速率图毫米 / 年到厘米 / 年量级看值域每个点一条时间序列什么时候开始沉、沉得匀不匀平面上看起来像一口口沉降碗方法是SBAS-InSARSmall Baseline Subset小基线子集。整条链路可以先记成二十多景 SLC → 按「短空间基线」组成几十个干涉对 → 每一对都做差分干涉 滤波 解缠 → 检查、删掉坏对 → 第一次反演去轨道趋势得到粗速率 → 第二次反演按大气的统计特性抠大气 → 地理编码 → 速率图 矢量点和积石山那课相比单对干涉的物理没变还是那五种相位难在对很多、网要连对、还要在时间上把形变和大气拆开。先和上一课对齐DInSAR vs SBAS积石山那课DInSAR这课SBAS影像数量2 景21 景干涉对数量1 对组网后61 对时间跨度4 天同震约半年要的产品一张同震位移图约 7 cm 抬升速率 时间序列煤矿持续下沉组网无核心步骤谁和谁配对轨道趋势一对上做 RPF 控制点每对都可能有斜坡靠 GCP 第一次反演一起扣大气几乎不单独建模4 天太短第二次反演专门去大气形变模型一次事件不拟合速率默认线性速率网密了可用二次项计算量一台电脑也能跑完演示干涉 解缠约占全流程一半时间一句话DInSAR 是「前后两张照片相减」SBAS 是「很多张照片两两精选配对再把所有差分解成一条随时间变化的沉降曲线」。基础知识比两景差分多出来的那些词下面默认你已经知道相位、一个条纹、LOS、失相干、解缠、去平。这些在积石山那篇里写过。这里只把SBAS / 煤矿 / 陆探时序多出来的概念写细。时序 InSAR 在回答什么两景差分回答从 A 日到 B 日地面沿视线一共动了多少时序回答三件更麻烦的事速率平均每天 / 每年沉多少velocity过程是匀速还是某个月突然加快时间序列空间格局哪几个矿在沉、碗有多大、中心在哪要回答它们必须有多余观测。一对干涉只有一个方程分不清「真沉降、轨道斜坡、大气、DEM 误差」。几十对、且每对基线都不大才能把这些东西在最小二乘里拆开。这就是 SBAS 存在的理由。形变模糊度为啥煤矿爱用 L 波段教程一上来就强调半年 20 多景、还是 L 波段对矿区沉降「非常有效」。关键不是「看得更清楚」而是相邻两次成像之间地可以多沉一截而不把条纹缠死。相位每转一圈视线向位移大约是半个波长来回双程哨兵 C 波段陆探 L 波段教程口径公开参数约值波长~5.6 cm21 cm~24 cm半波长一圈~2.8 cm~10 cm~12 cm煤矿沉降可以很快一个月沉十几厘米并不罕见。C 波段下一对干涉就可能转过好几圈山区再一失相干解缠直接崩。L 波段一圈大约 10 cm监测上限大约是哨兵的 45 倍教程原话。所以地震那种「4 天、几厘米」C、L 都能做C 更灵敏煤矿那种「持续、幅度大、地表还可能翻新」L 波段更扛得住教程里把这个叫「成像模糊度 / 形变模糊度」。记住物理就行一圈能代表的位移更大 → 同样的沉降条纹更疏、更好解。极化HH、VV、交叉极化雷达发射和接收可以选水平H或垂直V电场极化符号信号强弱更适合同极化HH、VV强InSAR 常规首选交叉极化HV、VH弱一般不当形变主数据同极化里再分VV对地面微小起伏更敏感海面波浪、草地一类HH对城市、微小构筑物高度变化更合适教程认为更适合沉降监测陆探这套样例只有HH直接用。原则有同极化就用同极化HH / VV 对矿区沉降通常差别不大有什么用什么。空间基线再拆一层平行 vs 垂直积石山那课把「空间基线 741 m」当成一个数。SBAS 组网时要拆开空间基线 ├── 平行基线沿视线 → 主要制造 轨道 / 平地相位残余 └── 垂直基线垂直于视线→ 主要制造 地形残差相位垂直基线大 山陡→ 差分减完 DEM 仍花干涉条纹跟着山走平行基线不准 / 轨道文件差→ 整景斜坡陆探特别常见所以组网不是「基线越小越好」这么简单而是陆探轨道管粗垂直基线容易大 →优先短空间基线L 波段模糊度高、时间失相干慢 →时间基线可以放长。平原矿区教程举安徽反过来地形残差不是主矛盾大气和地表变化更烦更该缩短时间基线。SBAS 组网小基线子集名字里的 Small Baseline 就是不要让每一景都去跟超级主影像配那是老式 PS 思路而是只连接「空间基线够短、时间基线可接受」的像对连成一张网。每个点的沉降从这张网里的所有差分解出来。一张时空基线图怎么读纵轴空间基线米 点 一景影像 横轴时间 线 一对干涉线越短、越平空间基线小山区友好线在时间上拉得很开时间基线长L 波段往往仍能相干某景孤立、没线连出去它不参加解算网从左到右断成两截两段时间各沉各的中间对不上除非开「允许断开块」去插教程这套参数下的报表务必留下印象量数值影像21 景平均空间基线200 m最大空间基线878 m平均时间基线42.5 天每景大约连接6 对向前 3 向后 3总干涉对61经验约等于 景数 × 3对照积石山单对741 m组网之后平均只有 200 m这就是 SBAS「小基线」四个字的意义——不是陆探突然飞准了是你拒绝了那些又肥又远的配对。官方提示每景连接数要大于 5反演才稳。教程习惯设6。冗余度选项含义High凡是落在「时间 ≤ 100 天、空间 ≤ 临界基线的 5%」里的对能连就连数据量大Low按规则瘦身。可按最短时间连前后各 3 景或在 100 天内挑空间基线最短的若干对陆探最大的麻烦是空间基线 → 教程选 Low并在 100 天窗口里优先空间最短。为何不能「每景只与下一景配一对」组网时最容易产生的误解是既然影像已经按时间排好每一景只与下一景做一对干涉1–2、2–3、3–4…再把相邻位移累加不就得到整条沉降曲线了吗这样连成的是一条时间链不是一张网。对 SBAS 反演来说链条在方程数目、误差传播和参数可区分性上都不成立。未知量与观测方程SBAS 最终要估计的是各成像时刻相对于某一参考日的视线向位移或等价的平均速率。(N) 景影像选定一景为参考后待估形变未知量约为(N-1)个。每一对干涉提供一个观测方程含义仅为这两个时刻之间的相位差位移差是多少。若只连接相邻两景干涉对数目 (M \approx N-1)与未知量个数相当。方程组接近恰定能写出一个解但几乎没有多余观测无法做最小二乘意义上的「多数表决」。噪声、解缠跳周、单日大气都会直接写进结果没有第二条路径把它们识别出来。类比测量一条折线的总长若每两段相邻点之间只量一次中间任一把尺子偏了首尾对不上也无从判断是哪一段出错。相干性的两种干涉相干 vs 时序模型相干这是整课最容易混的点。软件里好几处都叫 coherence不是同一个数。干涉相干.cc时序模型相干第一次反演产品何时出现每一对干涉图把所有对反演成速率之后在问什么这两景回波有多像01观测到的形变和你选的模型线性 / 二次贴不贴低意味着植被、叠掩、大基线、噪声非线性沉降、解缠错、大气没去净、模型选错掩膜用途解缠时去掉烂像素地理编码时去掉「速率不可信」的点教程后半段把时序相干阈值调到0.15是怕 0.2 把矿心挖空——矿心沉降快、地表破碎时序相干往往偏低低不代表那里没沉代表那里难拟合。准相干散射体 / 允许断开块传统做法某个像素只要在任何一景里是空值后面所有产品这个位置全废。SARscape 的优势之一Allow disconnected blocks true准相干散射体一类思路。意思是10 景里有 1 景这块失相干不必整条时间序列扔掉。用其余 90% 有值的观测拟合趋势把空档插上。组网时打开它第一次反演里还有一个配套参数用 90% 的有效值去内插 10% 的无效值min valid percent 90。太贪比如 50% 也插会编造沉降。解缠算法MCF vs Delaunay MCFSBAS 要对几十张干涉图解缠算法选错会又慢又错。MCF最小费用流Delaunay MCFDMCF做法整幅都解连低相干、水体也走先按相干性掩膜再在好像素的三角网上解适合相干性好、噪声少相干性差、山区、矿区破碎速度更快更慢、更稳教程建议电脑一般、相干还行电脑好、想稳 → 用 DMCF和积石山那课「调相干阈值再解」是同一件事的工业化版本DMCF 把「先掩再解」写进算法。并行SARscape 5.6.2 以后可同时解多对。经验约1000×1000 像元占 1 GB。30 m 产品大约一景解缠3 GB。32 线程、112 GB 内存的机器教程用16 线程约 48 GB避免把内存打满。叠掩 阴影掩膜务必打开true。陆探基线大叠掩区相位完全混杂阴影区没有回波不解比乱解强。煤矿沉降碗地在往下沉图为啥还是平面的和积石山「向四周扩散」是同一类问题地质机制不同。开采前 开采后 地表 —————— 地表 ——\______/—— ← 沉降碗 岩层 岩层 塌 陷 煤层 ████████ 煤层 采空中心下沉最多往外逐渐回到 0 → 平面图上是一口碗 / 一圈圈等值线以垂直沉降为主边缘会有指向碗心的水平位移地面往中间挤雷达斜视LOS 里会混进这一点水平分量平原矿区仍可近似勾Vertical Displacement山区矿区 陡坡水平分量不能假装没有解释要谨慎所以最终图上「一块红/蓝往外变淡」不是土在平面上往四周流是下沉量从中心向外变小。大气的「时间高通、空间低通」第二次反演不靠气象模型靠一条统计假设教程原话的意思维度大气长什么样滤波时间不是天天下雨今天湿、后天干像随机跳高通去掉时间上跳来跳去的成分空间同一天、同一景附近往往同晴同阴低通去掉单景里空间上很平滑的成分形变煤矿通常是时间上持续、空间上局部成碗。大气则是时间上不连续、空间上成片。两者频谱不同才能滤。风险滤波窗口不合适会把缓慢、成片的真沉降当成大气抠掉或把噪声留着。教程认为山西这块大气不严重而且数据只有半年所以空间窗口不要默认那么狠时间窗口从 365 天改成 180 天数据集就这么长。完整处理流程软件ENVI 5.6 SARscape 6.1。路径只能英文、不能空格、尽量不要小数点空格用下划线代替。源代码不认中文路径。① 导入数据干什么把二十多景陆探文件夹批量读成 SLC。哨兵可以对着 zip 读陆探必须先解压成文件夹。菜单Import Data → Spaceborne Single Sensor → 陆探一号课里技巧不要一景景点进去。用Select file list / 选文件夹指到上一级目录软件递归进子文件夹读meta.xml。极化选HH。一景导入后大约 8 个文件文件去留无后缀或_SLC.sml.hdr必须留二进制 SLC 两套头ENVI SARscape XMLshp kml覆盖范围示意拷贝数据时可不带本质是 ENVI 的.dat栅格外加 SARscape 的 XML。结束时你有什么input_data一类英文目录里21 套 SLC 三件套。② 制作 DEM干什么差分时减地形、地理编码时投坐标都要同一份 SARscape 认识的椭球 DEM。两种来源你手里的 DEM怎么转已是 ENVI 格式ENVI → SARscape选Ellipsoid DEMGeoTIFFOther format → TIFF同样选 Ellipsoid DEM硬性规定输出文件名必须以_dem结尾不要再加.tif这类英文点后缀用 Change output name 改坐标系尽量WGS84不要 CGCS2000、不要 UTM开源 DEM 最好大约30 m教程用 Copernicus DEM干涉多视分辨率应和 DEM 一致。DEM 30 m、影像硬输出 5 m等于用 30 m 高程去内插 5 m地形是假细的结束时你有什么一份xxxx_dem后面组网、干涉、精炼、地理编码反复用。③ 干涉组网Connection Graph干什么决定谁和谁做干涉。这是 SBAS 相对 DInSAR 多出来的第一步也是陆探最要动脑子的一步。菜单SBAS 工作流里选全部 SLC → Connection Graph不要让它瞎自动课里参数山区煤矿 陆探参数课里取值为什么最大时间基线100 天L 波段时间可以放长空间基线阈值临界基线的5%陆探临界大约 2 万米轨道间距又常能收到 1000 m 内先设大一点再靠 Low 瘦身冗余度Low不需要 High 那么多对每景连接数6官方建议 5向前 3 向后 3或 100 天内空间最短的 6 条连接偏好短空间山区地形残差是主矛盾Allow disconnected blockstrue准相干90% 有值也能出点平原矿区请改策略短时间基线优先不要死磕最短空间。看什么时空基线图纵轴跨度尽量短横轴可以拉长报表平均 / 最大空间基线、平均时间基线、总对数主影像日期这套样例落在2023-04-10附近报表会写坐标系陷阱connection graph文件夹里的矢量范围是WGS84 地理坐标SLC 本身是像元坐标左上角为 0方位向一条轴、距离向一条轴没有投影。两者叠不上是正常的不是数据坏了。结束时你有什么一张网 61 对清单。还没有彩虹图。④ 批量干涉处理干什么对网上每一对做配准 → 干涉 → 去平 → 减 DEM → 多视 → 滤波 → 解缠。积石山那课的 ③④⑥ 在这里对 61 对批量发生所以能占掉全流程一半时间。菜单Interferogram generationSBAS 流程里的干涉一步头文件用组网生成的那个DEM 用 ②。分辨率设30 m与 Copernicus DEM 一致。软件自动算多视比。矿碗尺度远大于 30 m要看碗沿更锐的条纹可试 15 m计算量和噪声都会上去。教程为演示用 30 m。Goldstein 滤波陆探相干已经很好不要滤太狠参数课里理由滤波窗口32偏小大窗口会糊积石山单对曾用 64那是山区噪声和解缠在打架Smooth window偏大细节多留一点Alpha自适应强度2.0默认常 2.5滤弱一点相干性估计窗口7默认 5窗口太小会高估相干、把噪声当好像素解缠相干好 → MCF相干差、机器强 →DMCF并行线程按内存估不要打满Apply layover and shadow mask true两个输出目录别看混位置内容工程根目录interferogram stacking快视图 TIFF方便不靠 SARscape 预览work / interferogram stacking真正拿来算的数据.cc.dint.fint.upha等快视图上的大黑边陆探轨道漂移各景覆盖不完全重合软件只留公共区域其余是空值。单对文件含义和积石山同一套文件含义original_dint已去平、已减地形未滤波.fintGoldstein 之后mask_cc叠掩 / 阴影已涂黑.upha解缠结果未滤波就能看清条纹79 天、L 波段相干很好。矿区形变包也还在。花掉的地方多是叠掩、透视收缩 这对空间基线偏大。沿距离向、与方位向平行的宽条纹 轨道相位残差。陆探初始轨道质量较差精密轨道会好一些但时序里仍常要靠后面 GCP 精炼不能指望自动 refinment那是哨兵更常用的。结束时你有什么61 套干涉产品。工程文件夹已经很大这一步完成后再挪盘是允许的。⑤ 检查干涉图与解缠必要时删对干什么SBAS 不怕少一对怕一对解缠错把整网基准带偏。所以不要闭眼进反演。解缠错误长什么样解缠假设整张图从一个零点平滑长出来。错了常见两种局部割裂某条缝左右亮度 / 基准突然不一致跳了整整一截整幅对切一半亮一半暗像被刀切开 —— 参考点在两个区域里对不上矿区在碗里、错误在远处的山可以留。错误落在你要报的矿必须处理。两种修法做法何时改解缠参数重跑这一对多数对都有问题或还能救SBAS Tools → Edit Connection Graph→ 选中这对Remove只有个别对坏教程举例主 1 与辅 3即 2023-01-24 与 2023-02-01 一类删对后不要忘了后面第一次反演若已跑过相关步骤要重跑。结束时你有什么一张「坏对已剔除 / 错误不在矿区」的网。⑥ 第一次反演GCP 精炼 线性速率干什么三件事捆在一起用控制点去掉每对上还在的轨道斜坡陆探几乎必须手动 GCP不要用哨兵那套 auto refinement在时序上按你选的模型默认线性解出粗速率可选估残余高程这套矿区教程关掉这是积石山那课「轨道精炼 重去平」的时序版钉子仍是稳定点但方程变成「很多对一起解」。选 GCP和地震那课同一套纪律更要避开矿碗尽量全局散开不要点在形变区沉降碗里绝对不点不要点在地形剧变处平坦、高相干、相位干净、不要空值三个真值得讲的参数参数课里含义产品相干阈值0.1时序模型相干不是解缠.cc。观测和线性模型的拟合程度形变模型Linear景不多、质量一般时更稳。组网密度到每景 1012 对才考虑Quadratic挖非线性估计残差高false矿区不是城市楼高Copernicus DEM约 2015且 30 m 对齐够用有效比例90%用 90% 有效观测插 10% 空缺解缠相干阈值与干涉一步对齐课里说干涉若为 0.1这里用 0.2 一类以菜单实际为准仍用.cc挡烂像素先看精炼后的干涉图轨道条纹应明显变弱。若某对教程举例 814精炼后斜坡还在回到 ⑤ 删掉再重跑这一步。演示课没抠到极致论文级要更狠。再看三个产品文件看什么velocity粗速率。直方图应接近零均值正态矿区应负值一侧略多以下沉为主precision_velocity越大越不准。用来掩膜coherence时序模型相干直方图若出现双峰 / 多个鼓包多半是解缠基准不一致不是真的有两种地质。用阈值「试砍」决定后面地理编码怎么掩教程现场调的结果请当经验不要当真理阈值现象他怎么改精度 40矿心被盖住太多放到60时序相干 0.2矿心一块掉成空洞放到0.15记住这两个数⑧ 要填进去。对称拉伸看图矿碗应干净右上角若还有缓坡是轨道没去净可再细抠 GCP 或删对。山上和地形相关的花纹若不在关注区可以忍。结束时你有什么粗速率、精度、时序相干轨道大斜坡基本压住。大气还在。⑦ 第二次反演去大气干什么按「时间上跳、空间上成片」的统计假设做时空滤波把大气从粗速率 / 时间序列里抠掉。参数课里怎么动参数默认倾向课里有效比例90不动空间低通窗口米偏小教程说默认约 1200滤得太狠不要滤太强山西大气不重时间高通窗口常 365 天改成180 天数据只有半年执行完看second inversion里的velocity也可叠 LOS 速度。直方图仍应零均值、负值略多正值一侧许多解缠噪声会被压下去。山上仍可能残留大基线造成的解缠错——不在矿区可忽略。结束时你有什么去大气后的速率这是准备地理编码的主产品。⑧ 地理编码干什么斜距产品投到经纬度输出栅格 点矢量。把 ⑥ 试砍得到的阈值填回来参数课里说明时序相干0.15保矿心速率精度60太大则烂点进图太小则矿心被挖高程精度很大SBAS 顺带估的地表高程和沉降不是同一回事放宽以免误伤形变点水体掩膜0这片没有湖有湖再开垂直向平原煤矿、以沉降碗为主勾Vertical Displacement。LOS 仍会混一点水平分量解释时写清楚假设。矢量每个 shp 最多约100 万点600 万点就会拆成 6 个文件。other format选 QGIS 能读的类型。文件夹名称内容未加 VDLOS 结果加了 VD垂直向结果结束时你有什么地图坐标下的速率栅格 一堆点 shp。⑨ QGIS 制图栅格 矢量点栅格速率图若已有单独教程这课主讲矢量点。推荐流程把所有VD的 shp 拖进 QGIS符号化字段选velocity对称、等间隔分类类别和色带按值域对齐色带可用教程分享的干涉彩虹rainbow一类在样式管理器里导入 XML设好一个图层后复制样式 → 其它图层粘贴样式保证所有分文件同一套色点大小用「地图单位 / 米」不要用毫米ArcGIS 常用毫米缩放时点会乱变。QGIS 用米和地面比例尺绑死。目的点大小课里全局图、让矿碗跳出来40 m干涉分辨率 30 m略大于像元局部叠影像、看点和地物套合20 m常见圆点效果15 m 圆形符号底图可叠 Google Earth 一类影像再加标注。每一步减掉什么、留下什么步骤主要减掉 / 压掉主要留下③ 组网空间基线过大、时间过长的「肥对」平均约 200 m 的小基线网④ 干涉地形DEM、椭球平地、噪声滤波、叠掩阴影每对上的形变 轨道残余 大气⑤ 删对解缠崩掉的个别对基准还能对齐的网⑥ 第一次反演轨道斜坡不符合线性模型的乱跳部分粗速率 矿碗⑦ 第二次反演时间上跳、空间上成片的大气更干净的速率 / 序列⑧ 地理编码低时序相干、低精度、可选水体地图上的沉降碗五种相位在时序里的分工平地/轨道 → ④ 去平 ⑥ GCP 精炼 地形 → ④ 减 DEM组网时已经用短垂直基线预防残差 形变 → ⑥⑦ 解出来的速率和序列 大气 → ⑦ 时空滤波 噪声 → ④ 滤波 ⑤ 删对 ⑧ 掩膜做到哪算这一步成功导入21 套 SLC 三件套路径纯英文极化 HHDEM文件名_demWGS8430 m组网平均空间基线明显小于单对乱配每景约 6 连山区纵轴短、横轴可长报表对数约 景数×3干涉未滤波短时间对条纹已清晰矿碗位置有局部条纹平行于方位向的宽条纹允许有后面精炼扣解缠矿区内部连续远处山的割裂可忍整幅对切的对要删第一次反演精炼后斜坡大减velocity直方图单峰、近零均值、负值略多矿碗干净记下相干 0.15、精度 60 这类试砍结果第二次反演正值噪声被压矿碗还在没有把碗心滤没地理编码LOS / VD 分清掩膜后矿心仍在制图分文件同一色带点大小用米碗和影像能套上和积石山那课怎么对照同一套陆探、同一类「轨道差、基线肥、L 波段电离层 / 斜坡」对策从「一对上死磕」变成「网上优胜劣汰」。问题积石山 DInSAR煤矿 SBAS空间基线 741 m 这种肥对只有一对只能硬做组网时尽量不连平均压到 ~200 m整景轨道斜坡RPF 控制点多项式每对 GCP 精炼坏对直接删除滤波窗口曾用 64为了救解缠相干已经很好窗口 32、alpha 调低怕削矿碗解缠阈值0.05因为相干从原始图估、山区要多留批量 DMCF / 叠掩阴影掩膜大气4 天几乎不管第二次反演专门滤垂直向假设逆冲抬升沉降碗平原可勾 VD产品一张位移图速率 点 可导出时间序列组网口诀可以钉在墙上山区陆探短空间、时间可长。平原煤矿短时间更重要。每景至少连 56 对。宁删一对勿留一张整幅切开来的解缠。术语对照表口播 语音识别会错。对照字幕 / 口语正确说法space inside / SPARS / SPARA / sparksSBAS-InSARSASK / SaaS skip / SAAS蒂普SARscapeN 维 / N维他ENVI路探 / 鹿探 / 录探 / 路泰陆探一号 LT-1同级化 / 同计划同极化 HH / VV交叉计划交叉极化成像模糊度 / 形变模糊度半波长对应的形变上限一圈空间极限 / 时间期限空间基线 / 时间基线临界极限临界基线组网冗余度Connection graph redundancy准相干散射体 / allow disconnected block允许断开块部分失相干仍拟合时空基线图时间–空间基线网络图ghost st / ghost down / Go downGoldstein 滤波解馋 / 检查 / snapping相位解缠snaphu / MCF德劳内是 MCFDelaunay MCF蝶眼 / 点点 / 叠眼叠掩 layover透视收缩foreshorteningauto reflementauto refinement 自动轨道精炼first in versionfirst inversion 第一次反演possession / v position / repositionvelocity precision 速率精度残差高residual height 残余高程QU / 二次方模型quadratic 非线性项second step / second inversion第二次反演去大气low pass / high PS空间低通 / 时间高通jo codinggeocoding 地理编码沉降晚沉降碗QG4QGIS阿克西斯ArcGIS比例米地图单位米coo p d e mCopernicus DEMELLIPSODMEllipsoid DEM建议的学习顺序先有积石山那课的底座五种相位、条纹疏密、LOS、解缠、控制点不点形变再搞懂为啥要组网一对方程拆不开速率 / 轨道 / 大气小基线网才拆得开陆探组网策略山区短空间L 波段允许时间拉长每景 ≥ 5 连两种相干.cc管解缠时序相干管「线性模型可信不可信」两次反演第一次钉轨道和粗速率第二次按大气频谱抠天气最后才是 QGIS 点大小那是表达不改计算结果数量级手算和矿碗对得上才踏实陆探一圈约10 cmLOS教程口径若某对只隔 8 天、碗心已经转过两圈说明这 8 天沉了约 20 cm 量级C 波段会更难解半年速率图上若碗心只有每年几毫米而现场矿一直在采先怀疑掩膜、模型、大气滤波把真信号吃了一句话收束这课不是在把积石山流程重复 21 遍而是在讲用半年、二十多景陆探 HH 条带在山区煤矿做 SBAS。因为 L 波段一圈大约 10 cm、时间失相干慢组网可以时间放到 100 天、空间尽量短因为轨道差必须手动 GCP坏对直接从网上剪掉因为要的是持续沉降而不是一次地震必须先线性反演、再按「时间跳、空间成片」去大气最后在沉降碗以垂直为主的假设下投到地图上。SBAS 很多小基线干涉对组成的网 最小二乘时序。组网 主动丢掉肥基线换山区能用的相干。第一次反演 轨道精炼的时序版 粗速率。第二次反演 用统计滤波打大气。沉降碗 采空区上覆岩层塌陷中心沉得多、向外变少——平面图上的圈仍是等位移线不是土在往四周流。对照教程整理。参数以 SARscape 6.1 菜单和这套大同样例为准换平原矿区、换更长时序时组网和时间窗口要改不要把 100 天 / 180 天 / 0.15 / 60 当成全世界通用的魔法数。