一种利用TanDEM-X双站InSAR提取地震滑坡体积的方法
【技术领域】
[0001] 本发明属于卫星遥感应用领域，具体涉及一种利用TanDEM-X双站InSAR提取地震 滑坡体积的方法，利用TanDEM-X双站InSAR获取地震前后地表高程变化，进而定量估算区 域地震滑坡体积，可在全球区域拓展应用。
【背景技术】
[0002] 近年来，地震容易引发大量滑坡、泥石流等次生地质灾害，对人民居住环境和生命 财产安全构成了巨大的威胁，同时也是铁路、公路等交通设施的重要病害，是震区面临最严 重的两类地质灾害，进行区域地震滑坡体积提取和泥石流危险性评估是地震灾区恢复重建 和经济发展的紧迫工作。其中滑坡是指斜坡上的土体或者岩体，受河流冲刷、地下水活动、 雨水浸泡、地震及人工切坡等因素影响，在重力作用下，沿着一定的软弱面或者软弱带，整 体地或者分散地顺坡向下滑动的自然现象。其中泥石流是指在山区或者其他沟谷深壑，地 形险峻的地区，因为暴雨、暴雪或其他自然灾害引发的山体滑坡并携带有大量泥沙以及石 块的特殊洪流。地震滑坡体积作为泥石流物源定量计算的重要参数，开展地震滑坡体积定 量研究对估算泥石流物源具有重要的实际应用价值和科学意义，为震区泥石流灾害的风险 评价和预测预报等提供科学依据。地震发生后，地震滑坡、泥石流等次生灾害频发，迫切需 要一种能够大范围、定量、精确地提取地震滑坡体积的方法，进而能够地震滑坡等地质灾害 进行预警，减少百姓财产损失和保障人员安全。
[0003] 目前，地震滑坡体积提取的技术手段主要有GPS实地测量、三维激光LiDAR、光学 遥感等。GPS实地测量和三维激光LiDAR具有观测精度高等优势，尤其三维激光LiDAR能够 获取整个滑坡面的高分辨率点云地形数据，可精确获取地震滑坡体积，但地震滑坡个数通 常上万个，甚至几十万个，该方法提取区域范围的地震滑坡体积效率较低；光学遥感可以大 范围获取地震滑坡中心点位置、边界、面积等分布参数，通过经验公式或模型等提取地震滑 坡体积，但经验公式或模型对最终地震滑坡体积带来不可避免的误差，难以定量、精确的获 取地震滑坡体积。且地震发生后，较长一段时间会存在持续降雨天气，光学遥感的地震滑坡 体积提取工作难以开展。
[0004] TanDEM-X双站InSAR可获取全球尺度的高分辨率、高精度地形信息，结合震前地 形资料，能够精确计算地震地震滑坡体积，为地震滑坡风险性评价以及地震滑坡灾害防治 提供重要的基础数据。
【发明内容】
[0005] 本发明要解决的技术问题是针对目前地震滑坡体积提取方法存在的上述不足，提 供一种利用TanDEM-X双站InSAR提取地震滑坡体积的方法，能够大范围、定量提取地震滑 坡体积。
[0006] 本发明为解决上述技术问题所采用的具体步骤是：
[0007] -种利用TanDEM-X双站InSAR提取地震滑坡体积的方法，根据TanDEM-X双站 InSAR获取的地震前后地表高程变化结果，结合地震滑坡边界编目数据和地震前后光学遥 感数据，求取地震滑坡体积，具体包括以下步骤：
[0008] 1)双站SAR数据影像配准：使两幅复图像中同一位置的像素对应地面上的同一回 波点；
[0009] 2)双站InSAR干涉相位计算：利用TanDEM-X和TerraSAR-X两颗卫星获取的高分 辨率SAR影像，通过干涉图生成数据处理，计算双站InSAR干涉相位；
[0010] 3)双站InSAR差分相位估算：双站InSAR干涉相位与外部DEM模拟相位进行差分 处理；
[0011] 4)双站InSAR差分相位解缠：相位解缠采用兼顾运算速度和准确性的最小费流量 法（MinimumCostFlow，MCF，双站InSAR差分相位解缠主要解决干涉条纹过密引起的解缠 相位跳跃、不连续，以及轨道基线误差引起的长波趋势面等问题）；
[0012] 5)双站InSAR差分相位校正：采用高精度GPS或ICESat点进行绝对相位偏移校 正；
[0013] 6)地震前后地表高程变化结果估算：双站InSAR差分相位经过相位高程转换，估 算地震前后地表高程变化结果；
[0014] 7)地震滑坡体积求取：地震前后地表高程变化结果经过地震滑坡边界编目数据 筛选，以及地震前后光学遥感数据进一步识别获取，求取地震滑坡体积。
[0015] 按上述方案，所述步骤1)中影像配准是为获得干涉相位具有较高的信噪比，使两 幅复图像中同一位置的像素对应地面上的同一回波点的过程。
[0016] 按上述方案，所述步骤2)中干涉图生成数据处理是影像经过精确配准后，进行复 共辄相乘生成复干涉纹图，复干涉相位生成公式如下：
[0017]
[0018] 式中，uint是干涉结果，5；(#)和戶)分别为两颗卫星获取的雷达信号，AR为两 颗卫星到地面P点的距离之差，A为雷达波长。
[0019] 按上述方案，所述步骤3)中外部DEM模拟相位通过外部DEM经过相位高程转换获 取。
[0020] 按上述方案，所述步骤3)中双站InSAR干涉相位与外部DEM模拟相位进行差分处 理公式如下：
[0021] ^^Differences- ^ ^TanDEM X_^ ^SRTMDEM
[0022] 式中，AW Differences为差分相位，AW TanDEM X是双站InSAR干涉相位，AW SRTMDEM为外 部DEM模拟相位。
[0023] 按上述方案，所述步骤5)中高精度GPS或ICESat点是覆盖研究区域的高精度测 量点。
[0024] 按上述方案，所述步骤6)中相位高程转换为根据相位P和高程z(y)之间的关系 获取相应高程变化结果，满足下式：
[0025]
[0026] 式中，Bi为两天线之间的空间矢量在平行主影像上的分量，0为视角，/?(/》为两 颗卫星到地面P点的平均距离。
[0027] 按上述方案，所述步骤7)中地震滑坡边界编目数据通过光学遥感、实地观测技术 手段获取。
[0028] 按上述方案，所述步骤7)中地震前后光学遥感数据是覆盖地震滑坡区域的地震 发生前后的星载、基站的光学遥感数据。
[0029] 按上述方案，所述外部DEM包括30m和90m分辨率的SRTM-CDEM、30m分辨率的 SRTM-XDEM以及其它测量技术手段获取的DEM结果。
[0030] 本发明与现有技术相比具有以下优点：
[0031] 1、TerraSAR-X和TanDEM-X为两颗高分辨率SAR卫星，地震前后地表高程变化结 果是校正后的双站InSAR差分相位经过相位高程方法生成，相位级的高程转换使地震前后 地表高程变化结果具有高精度特性。地震前后地表高程变化结果又经过地震滑坡边界编目 数据筛选和地震前后光学遥感数据识别，生成的地震滑坡体积具有高分辨率、高精度特征， 故本发明利用TanDEM-X双站InSAR提取地震滑坡体积的方法具有高分辨率、高精度优势；
[0032] 2、InSAR具有全天时全天候成像能力，不受云层及昼夜影响。尤其在地震发生之 后，震区多数发生长时间、大规模的降雨，导致光学遥感无法成像，而SAR在阴雨天气等恶 劣环境仍能进行影像获取，进而对地震滑坡持续监测；
[0033] 3、本发明具有面状、大区域的、滑坡地震滑坡体积提取能力，而GPS、三维激光 LiDAR等地面测量手难以进行；
[0034] 4、本发明成本低，不需要布设地面测量设备。
【附图说明】
[0035] 图1是本发明的方法流程图；
[0036] 图2是本发明方法应用的一实例图；
【具体实施方式】
[0037] 下面结合附图和实施例对本发明作详细具体的说明。
[0038] 参照图1所示，本发明所述的利用TanDEM-X双站InSAR提取地震滑坡体积的方 法，通过TanDEM-X双站InSAR干涉相位、外部DEM模拟相位、高精度GPS或ICESat校正点 获取校正后的TanDEM-X双站InSAR差分相位，进而获取地震前后地表高程变化结果，再结 合地震滑坡边界编目数据、地震前后光学遥感数据求取地震滑坡体积，结合滑坡面积参数， 求取高精度高分辨率地震滑坡物源厚度。
[0039] 双站InSAR干涉相位由TanDEM-X和TerraSAR-X两颗卫星SAR数据经过影像配 准、干涉图计算等步骤获取，外部DEM模拟相位通过SRTMDEM等外部DEM经过相位高程转 换获取，高精度GPS或ICESat校正点是覆盖研究区域的高精度测量点，校正后的双站InS