一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法技术

本发明专利技术公开了一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，利用SAR强度偏移量追踪算法从两景SAR强度影像中生成矿区地表雷达视线向和方位向形变；基于矿区地表水平移动与下沉梯度的比例关系建立东西、南北方向水平移动与下沉之间的函数模型；分别建立矿区地表下沉与视线向和方位向形变的观测方程组；利用最小二乘法估计求解矿区地表下沉；基于求解的地表下沉，使用矿区地表水平移动与下沉梯度的比例关系估计矿区地表在东西、南北方向的水平移动。实现了仅利用两景SAR强度影像估计矿区地表大量级三维形变，有效地克服了传统方法无法从单个InSAR干涉对中获取大量级三维形变的局限。

【技术实现步骤摘要】
一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法
本专利技术涉及一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法。
技术介绍
矿区地表在东西、南北和垂直方向的三维形变监测对于提前评估和控制地下开采导致的潜在地质灾害和建构筑物损坏有着重要作用。合成孔径雷达干涉测量(interferometricsyntheticapertureradar,InSAR)技术是一种新兴的雷达遥感技术。该技术基于两景SAR单视复数(singlelookcomplex,SLC)影像生成的单个InSAR干涉对，实现了地表雷达视线方向(line-of-sight,LOS)形变的大范围、高精度、低成本、非接触监测。然而，由于InSAR技术仅能从单个InSAR干涉对中获取一维LOS向形变，而非真正的地表三维形变。从而大大地制约了InSAR技术在矿区的应用前景。为了获取矿区地表三维形变，传统的方法至少需要三个来自于不同雷达成像几何学SAR数据的InSAR干涉对，即至少需要六景SLC影像。受限于目前较少的可用SAR卫星数量，该方法在实际应用中很难满足。专利技术专利ZL201210440875.4提出了一种利用单个InSAR干涉对获取矿区地表三维形变场的方法，该专利技术有效地克服了传统方法无法从单个InSAR干涉对中获取矿区地表三维形变的局限，大大降低了传统方法对于数据的苛刻要求，为矿区地质灾害和建构筑物风险评估提供了重要的参考依据。由于InSAR技术是基于两景SLC影像的干涉相位的，所以当干涉相位失相关时，InSAR技术则无法准确地获取地表雷达视线向形变。然而，大量级形变常常导致InSAR干涉相位失相关，从而致使专利技术专利ZL201210440875.4无法获取矿区大量级三维形变。
技术实现思路
本专利技术提供了一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，其目的在于克服传统基于单个InSAR干涉对的方法无法获取矿区地表大量级三维形变的局限，利用矿区地表水平移动与下沉梯度的比例关系从SAR强度偏移量追踪算法得到的矿区地表视线向和方位向形变中估计矿区地表沿着垂直、东西和南北方向的三维形变，大大地拓宽了SAR技术在矿区的应用前景。一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，包括以下几个步骤：步骤1：利用SAR强度偏移量追踪算法从两景SAR强度影像中生成待监测矿区地表雷达视线方向dLOS和方位方向形变dAZI；待监测的矿区的形变图的总行数和总列数分别为n和m；步骤2：基于矿区地表水平移动与下沉梯度的比例关系建立东西、南北方向的水平移动UE、UN与下沉W的函数关系；步骤3：分别根据视线向形变dLOS和方位向形变dAZI与地表真实三维形变的投影关系建立矿区地表下沉与利用偏移量追踪算法获取的视线向与方位向形变之间的投影方程；步骤4：按照步骤3构建的投影方程建立待监测矿区地表各像素点的下沉值与视线向形变的观测方程；由于矿区边缘受地下开采影响较小，在待监测矿区最后一行和最后一列的水平移动设置为0；步骤5：求解待监测矿区各像素点的下沉值；步骤6：基于求解的下沉值并利用步骤2中所述的东西、南北方向的水平移动与下沉梯度的关系估计待监测矿区地表在东西、南北方向的水平移动。所述利用偏移量追踪算法获取的视线向与方位向形变中的偏移量追踪算法属于现有技术。进一步地，所述基于矿区地表水平移动与下沉梯度的比例关系建立的东西、南北方向的水平移动UE、UN与下沉W的函数关系如下：其中，i和j为地表点的像素坐标；W(i,j)表示地表点(i,j)处的下沉值，CE和CN分别示地表点(i,j)处在东西、南北方向的下沉梯度比例系数，b、H和β分别表示矿区水平移动系数、采深和主要影响角，b、H和β从监测矿区获取；RE和RN分别表示视线向和方位向形变图在东西和南北方向的空间分辨率，RE和RN从视线向和方位向形变图中获取。进一步地，所述矿区地表下沉与利用偏移量追踪算法获取的视线向与方位向形变之间的投影方程如下：和/或其中，为投影系数矩阵；A11(i,j)＝-A12(i,j)-A13(i,j)，A21(i,j)＝cosθ-A22(i,j)-A23(i,j)；A12(i,j)＝sin(α-3π/2)CE(i,j)，A22(i,j)＝-b·CE(i,j)·sin(α-3π/2)；A13(i,j)＝cos(α-3π/2)CN(i,j)，A23(i,j)＝-b·CN(i,j)·cos(α-3π/2)；α和θ表示雷达飞行方位角和入射角，从SAR强度影像的头文件中获得。用方位向或者视线向监测三维形变，即用任意一个投影方程来估计三维形变时，采用线性方程求解算法直接求解矿区下沉值；两者联合使用时，采用最小二乘求解算法直接求解矿区下沉值，联合使用时增加了多余观测，使得下沉值的求解结果精度更高。进一步地，所述待监测矿区地表各像素点的下沉值与视线向形变的观测方程如下：其中，Wi和Wi+1分别表示待监测的矿区形变图中第i行和第i+1行像素中各像素的下沉值；Wi＝[W(i,1),W(i,2),…,W(i,m)]T，Wi+1＝[W(i+1,1),W(i+1,2),…,W(i+1,m)]T；dLOSi和dAZIi分别表示待监测的矿区形变图中第i行像素中各像素的视线方向和方位方向形变观测值；dLOSi＝[dLOS(i,1),dLOS(i,2),…,dLOS(i,m)]T，dAZIi＝[dAZI(i,1),dAZI(i,2),…,dAZI(i,m)]T；Ci、Ci′、Di及Di′表示观测系数矩阵；Ci′＝diag[A13(i,1)A13(i,2)…A13(i,m-1)0]m×m,Di′＝diag[A23(i,1)A23(i,2)…A23(i,m-1)0]m×m其中，待监测的矿区形变图中最后一行像素中各像素的视线方向和方位方向形变观测值为D″＝diag(cosθcosθ…cosθ)m×m，diag表示对角阵。进一步地，按照下述公式求解待监测矿区各像素点的下沉值：利用最小二乘法求解：其中，Ei和Li表示中间变量矩阵，为Ei的转置矩阵，表示的逆矩阵；其中，最后一行各像素的下沉值按照公式计算获得：Wn＝(D″)-1·dLOSn。有益效果本专利技术公开了一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，利用SAR强度偏移量追踪算法从两景SAR强度影像中生成矿区地表雷达视线向和方位向形变；基于矿区地表水平移动与下沉梯度的比例关系建立东西、南北方向水平移动与下沉之间的函数模型；分别建立矿区地表下沉与视线向和方位向形变的观测方程组；利用最小二乘法估计求解矿区地表下沉；基于求解的地表下沉，使用矿区地表水平移动与下沉梯度的比例关系估计矿区地表在东西、南北方向的水平移动。本专利技术实现了仅利用两景SAR强度影像估计矿区地表大量级三维形变，有效地克服了传统方法无法从单个InSAR干涉对中获取大量级三维形变的局限。此外，该方法数据制约少，监测结果准确有效，大大拓宽了InSAR技术的应用前景。本方法其对于拓宽SAR技术的应用空间，降低矿区三维形变监测成本有着重要意义。此外，其对于指导矿区安全生产、分析矿区沉降机理、预警矿区地表地质灾害以及生态环境保护也起着重要作用。附图说明图1为本专利技术的所述方法的流程示意图。图2本文档来自技高网...

【技术保护点】
一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，其特征在于，包括以下几个步骤：步骤1：利用SAR强度偏移量追踪算法从两景SAR强度影像中生成待监测矿区地表雷达视线方向d

【技术特征摘要】
1.一种基于两景SAR强度影像的矿区地表大量级三维形变估计方法，其特征在于，包括以下几个步骤：步骤1：利用SAR强度偏移量追踪算法从两景SAR强度影像中生成待监测矿区地表雷达视线方向dLOS和方位方向形变dAZI；待监测的矿区的形变图的总行数和总列数分别为n和m；步骤2：基于矿区地表水平移动与下沉梯度的比例关系建立东西、南北方向的水平移动UE、UN与下沉W的函数关系；步骤3：分别根据视线向形变dLOS和方位向形变dAZI与地表真实三维形变的投影关系建立矿区地表下沉与利用偏移量追踪算法获取的视线向与方位向形变之间的投影方程；步骤4：按照步骤3构建的投影方程分别建立待监测矿区地表各像素点的下沉值与视线向形变和方位向形变的观测方程组；步骤5：求解待监测矿区各像素点的下沉值；步骤6：基于求解的下沉值并利用步骤2中所述的东西、南北方向的水平移动与下沉梯度的关系估计待监测矿区地表在东西、南北方向的水平移动。2.根据权利要求1所述的方法，其特征在于，所述基于矿区地表水平移动与下沉梯度的比例关系建立的东西、南北方向的水平移动UE、UN与下沉W的函数关系如下：其中，i和j为地表点的像素坐标；W(i,j)表示地表点(i,j)处的下沉值，CE和CN分别示地表点(i,j)处在东西、南北方向的下沉梯度比例系数，H和β分别表示矿区水平移动系数、采深和主要影响角，b、H和β从监测矿区获取；RE和RN分别表示视线向和方位向形变图在东西和南北方向的空间分辨率，RE和RN从视线向和方位向形变图中获取。3.根据权利要求2所述的方法，其特征在于，所述矿区地表下沉与利用偏移量追踪算法获取的视线向与方位向形变之间的投影方程如下：和/或

【专利技术属性】
技术研发人员：杨泽发，朱建军，李志伟，易辉伟，胡俊，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43