融合无人机DOM和星载SAR影像的地表三维形变提取方法技术

本发明专利技术公开了一种融合无人机DOM和星载SAR影像的地表三维形变提取方法，适用于地表形变及地质灾害监测领域。首先利用SAR或InSAR技术获取目标区地表LOS向形变；利用无人机获取目标区地表影像数据，生成具有相同分辨率的数字正射影像图DOM；利用精配准方法计算两期DOM上同名像素点在东西向、南北向的坐标偏移量，结合DOM分辨率，求得各像素点对应地表点的东西和南北方向的水平移动；将利用DOM获取的东西、南北向水平移动及LOS向形变代入SAR三维形变模型，解算地表竖向下沉值W，从而得到地表三维形变。本发明专利技术结合无人机DOM和星载SAR影像获取地表三维形变，其覆盖范围广、非接触地表、效果好，为地表三维形变监测提供了一种新方法。法。法。

【技术实现步骤摘要】
融合无人机DOM和星载SAR影像的地表三维形变提取方法


[0001]本专利技术涉及一种融合无人机DOM和星载SAR影像获取地表三维形变的方法，属于地表形变及灾害监测领域。
技术背景
[0002]我国幅员辽阔，自然环境多样，每年因地下资源开发导致的地表沉降、地质灾害量大、面广。比如：城市地区地表及地下工程建设、地下水资源开发导致地表沉降，进而影响地表建构筑物安全运营；煤炭、石油、金属等矿产资源开发后，会对矿区环境产生严重损害，形成地表塌陷、裂缝、滑坡等灾害；我国西部地区，尤其是云南、贵州、四川、青海等省份，地形起伏大，各类地质灾害频发。这些地表形变、地质灾害本质上是各地表点移动轨迹的综合反映，可以在三维空间中进行投影，分解为竖向移动和水平移动。竖向移动为下沉或隆起，水平移动可以按照垂直或平行于某一断面进行设置，比如：南北向及东西向水平移动。
[0003]传统地表形变监测方法如：GNSS、水准测量、全站仪等存在工作量大、点位密度不足、点位易破坏、成本高、不便于进行连续测量和自动化测量等缺点。合成孔径雷达(SAR)测量技术自上世纪90年代发展以来，有效弥补了传统监测技术的不足。目前SAR已广泛应用于区域灾害探测、监测领域。然而，由于该技术只能获取沿雷达视线方向地表形变，单轨道SAR影像在不借助外部数据或数学模型的情况下难以获取地表三维形变，极大地限制了其在建构筑物形变及地质灾害监测方面的应用。随着无人机技术的成熟与普及，无人机摄影测量已广泛应用于各行各业，其优势在于机动灵活、非接触被摄物体、分辨率高、速度快、精度高，但该技术难以获取地表竖直沉降，也缺乏直接将其用于地表水平移动获取的广泛应用。
[0004]为此，本专利技术结合无人机影像及SAR影像的优势，提出一种融合无人机DOM和星载SAR影像的地表三维形变提取方法，可以快速、高精度获取地表及建构筑物三维形变，应用前景广阔。

技术实现思路

[0005]本专利技术所要解决的技术问题是：提供一种融合无人机DOM和星载SAR影像的地表三维形变提取方法，解决了单轨SAR影像难以获取地表三维形变、无人机影像难以得到竖向下沉的问题，具有精度高、成本低、不接触被测物体、范围广、易操作等优点。
[0006]为实现上述技术目的，本专利技术的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于具体步骤如下：
[0007]S1，利用卫星SAR/InSAR技术解算目标地区的视线向形变场，记为：LOS；
[0008]S2，利用一架无人机按照一样的航路获取两个不同时期的目标地区地表影像数据，处理无人机影像生成数字正射影像图DOM，且两期DOM空间分辨率相同；
[0009]S3，将第一期DOM1作为主影像，第二期DOM2作为从影像，利用精配准方法计算两期DOM同名像点南北、东西方向的像素偏移量，从该偏移量中去除因为两次无人机获得的影像整体偏移量得到因地表移动而引起的像素偏移量；影像整体偏移量为两次航拍产生的误差
偏移量；
[0010]S4，利用S3得到的南北、东西方向的像素偏移量和影像地面分辨率计算每个同名像点相应地表实际水平移动，每个同名像点实际水平移动包括：南北方向水平移动量U
N
、东西水平移动量U
E
；
[0011]S5，根据SAR三维形变分解模型，结合卫星获取的视线向形变场LOS以及南北方向水平移动量U
N
、东西水平移动量U
E
，解算目标区域地表竖向下沉值W，从而得到地表实际三维形变。
[0012]2.根据权利要求1所述的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于，步骤S1中SAR/InSAR技术采用的是经典的偏移量跟踪算法、子带干涉方法、DInSAR、时序InSAR，可利用上述方法获取地表视线向形变量，记为LOS。
[0013]3.根据权利要求1所述的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于，步骤S2中生成的两期DOM地面分辨率与SAR影像分辨率应相同，否则需要对DOM进行重采样。
[0014]步骤S3中的精配准方法包括：归一化互相关匹配方法、最小二乘匹配方法、特征匹配方法；影像整体偏移量由非变形区域偏移量拟合的二次曲面获取，其主要是配准算法、噪声影响产生的系统误差。
[0015]步骤S4中计算各同名像素点的南北方向水平移动量U
N
、东西水平移动量U
E
，移动量的单位为像素点个数，具体为：
[0016]设第一期DOM1与第二期DOM2匹配到的同名点对为p1(x1,y1)和p2(x2,y2)，点p1(x1,y1)位于第一期DOM1上，点p2(x2,y2)位于第二期DOM2上，(x1,y1)、(x2,y2)分别是点p1(x1,y1)、点p2(x2,y2)在各自影像坐标系中的坐标，影像坐标系的原点为DOM的左上角，原点向右的方向为影像坐标系X轴方向，原点向下的方向为影像坐标系Y轴方向；利用公式U
N
(x,y)＝GSD*(y2‑
y1)和U
E
(x,y)＝GSD*(x2‑
x1)分别计算第一期DOM1与第二期DOM2中记录的地表水平移动；U
N
(x,y)是同名点对p1(x1,y1)和p2(x2,y2)在南北方向上的水平移动，U
E
(x,y)是同名点对p1(x1,y1)和p2(x2,y2)在东西方向上的水平移动，GSD为DOM的地面分辨率。
[0017]步骤S5中利用SAR三维形变分解模型解算地表竖直沉降W的公式为：
[0018][0019]其中，θ为雷达卫星入射角；α
h
为卫星航向角；U
N
、U
E
为无人机DOM计算得到的南北、东西方向水平移动；LOS为SAR/InSAR技术得到的雷达视线向地表形变。
[0020]有益效果
[0021]单轨道SAR技术只能得到高精度雷达视线向(LOS)形变，不能分解到竖直、东西及南北方向的三维形变，无人机影像形成两期DEM做差只能得到低精度竖直向变形，水平移动也缺乏研究和应用。本专利技术融合了无人机影像及SAR影像的优势，利用无人机影像精确配准方法得到水平移动，带入SAR视线向变形分解方程，便可将LOS形变分解获取高精度的竖直向变形，弥补了两者各自难以获取高精度地表三维形变的不足，有效获取了地表及建构筑物的三维形变，破解了单轨道SAR影像只能获取视线向形变的难题，也拓展了无人机摄影测量的应用领域，具有精度高、成本低、不接触被测物体、范围广、易操作等优点，为地表及建构筑物三维形变信息提取、地质灾害监测与预警提供了一种新的技术手段。
附图说明
[0022]图1为本专利技术融合无人机DOM和星载SAR影像的地表三维形变提取方法的实施流程图。
[0023]图2为本专利技术使用的模拟地表三维形变图。
[0024]图3为本专利技术解算出的地表三维本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于具体步骤如下：S1，利用卫星SAR/InSAR技术解算目标地区的视线向形变场，记为：LOS；S2，利用一架无人机按照一样的航路获取两个不同时期的目标地区地表影像数据，处理无人机影像生成数字正射影像图DOM，且两期DOM空间分辨率相同；S3，将第一期DOM1作为主影像，第二期DOM2作为从影像，利用精配准方法计算两期DOM同名像点南北、东西方向的像素偏移量，从该偏移量中去除影像整体偏移量得到因地表移动而引起的像素偏移量；影像整体偏移量为两次航拍产生的误差偏移量；S4，利用S3得到的南北、东西方向的像素偏移量和影像地面分辨率计算每个同名像点相应地表实际水平移动，每个同名像点实际水平移动包括：南北方向水平移动量U
N
、东西水平移动量U
E
；S5，根据SAR三维形变分解模型，结合卫星获取的视线向形变场LOS以及南北方向水平移动量U
N
、东西水平移动量U
E
，解算目标区域地表竖向下沉值W，从而得到地表实际三维形变。2.根据权利要求1所述的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于：步骤S1中SAR/InSAR技术采用的是经典的偏移量跟踪算法、子带干涉方法、DInSAR、时序InSAR，可利用上述方法获取地表视线向形变量，记为LOS。3.根据权利要求1所述的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于：步骤S2中生成的两期DOM地面分辨率与SAR影像分辨率应相同，否则需要对DOM进行重采样。4.根据权利要求1所述的融合无人机DOM和星载SAR影像的地表三维形变提取方法，其特征在于：步骤S3中的精配准方法包括：归一化互相关匹配方法、最小二乘匹配方法、特征匹配方法；影像整体偏移量由非变形区...

【专利技术属性】
技术研发人员：范洪冬，邹英杰，李运江，庄会富，谭志祥，邓喀中，
申请(专利权)人：内蒙古昊盛煤业有限公司，
类型：发明
国别省市：