CN109696152B - 一种低相干性区域地面沉降量估算方法 - Google Patents
一种低相干性区域地面沉降量估算方法 Download PDFInfo
- Publication number
- CN109696152B CN109696152B CN201910112992.XA CN201910112992A CN109696152B CN 109696152 B CN109696152 B CN 109696152B CN 201910112992 A CN201910112992 A CN 201910112992A CN 109696152 B CN109696152 B CN 109696152B
- Authority
- CN
- China
- Prior art keywords
- ground settlement
- area
- research
- low coherence
- research area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
- Image Processing (AREA)
Abstract
一种低相干性区域地面沉降量估算方法,所述估算方法是研究区遥感影像订购,研究区遥感影像处理,研究区地面沉降量反演,研究区地面沉降量数据检验,研究区低相干性区域地面沉降量插值方法确定,研究区低相干性区域地面沉降量估算,本发明为及时准确估算低相干性区域范围地面沉降量提供了一种新方法,实现了基于地理信息空间技术的低相干性区域地面沉降量的快速准确监测。
Description
技术领域
本发明涉及一种低相干性区域地面沉降量的估算方法,尤其是一种利用地理信息空间插值技术高效准确估算低相干性区域地面沉降量的方法。
背景技术
星载合成孔径雷达干涉测量技术(InSAR)因其高分辨率、高精度、全天候等优点已经成为目前广泛使用的地面沉降测量技术之一,其旨在通过计算两次或多次过境时合成孔径雷达影像的相位差来获取地面沉降数值,主要包括D-InSAR方法、PS-InSAR方法、SBAS-InSAR方法、DS-InSAR方法及MAI方法等。
然而,受成像时刻分辨率单元内散射体变化、传感器姿态变化和雷达波透射比等因素制约,雷达回波信号不同程度受到时空去相干的影响,导致干涉相位随机噪声增加,相干性衰减,地面沉降测量精度降低。低相干性区域目前并没有明确的定义,一般指空间或时间相干性小于某个临界值时的研究对象,如农田或牧草地是典型的低相干性区域。低相干性区域地面沉降量的准确估算成为目前亟需解决且具有挑战性的问题。
从理论上讲,将多源传感器数据进行联合解算以便增加观测数并同时提高时间分辨率是实现低相干性区域地面沉降监测的最好方法,但是,这不仅对不同平台的数据配准提出更高的要求,而且不可避免地衰减空间分辨率,更为麻烦的是,目前可用的传感器数据无法满足同一区域多次观测数的基本要求,因此发明一种准确估算低相干性区域地面沉降量方法具有重要的现实意义。
发明内容
本发明的目的在于提供一种低相干性区域地面沉降量估算方法,以克服传统InSAR技术无法准确提取低相关性区域地面沉降量的不足。
本发明的上述目的是通过以下技术方案来实现的。
一种低相干性区域地面沉降量估算方法,其特征在于:所述估算方法是按下列步骤进行的:
1)研究区遥感影像订购
从专业机构遥感影像订购网站,订购下载研究区域研究期间内多期Sentinel-1A遥感影像,将Sentinel-1A数据转换为SARScape内部格式,并使用研究区域范围矢量边界进行数据裁剪,得到研究区域多时相SLC数据;
2)研究区遥感影像处理
对研究区域多时相SLC数据,进行基线估算及小基线集生成、图像匹配、干涉图生成、去平地效应、干涉图滤波、相干系数计算、相位解缠、轨道精炼和重去平等处理,其中将相位解缠相干系数阈值设定为0.35,相干系数小于0.35的区域,即为低相干性区域;
3)研究区地面沉降量反演
对上述步骤2)处理得到的研究区遥感影像,进行轨道精炼和重去平、形变速率和高程系数估算、地形残余相位和大气相位去除以及地理编码,反演得到研究区域研究期间内不同时间段的地面沉降量;
4)研究区地面沉降量数据检验
对上述步骤3)中反演得到的研究区域研究期间内不同时间段的地面沉降量数据进行正态检验,检验其是否呈正态分布,若呈正态分布,则具有自相关性,可以使用克里金法进行插值,否则对地面沉降量数据进行变换,使其呈正态分布;
5)研究区低相干性区域地面沉降量插值方法确定
对上述步骤4)处理得到的研究区域某一时段地面沉降量,进行普通克里金插值与经验贝叶斯克里金插值,获得低相干性区域的地面沉降量,计算两种插值方法得到的低相干性区域边界地面沉降量与其相邻的高相干区域地面沉降量的均方根误差,均方根误差越小,插值精度越高,进而确定适用于研究区域低相干性区域地面沉降量的克里金插值方法;
6)研究区低相干性区域地面沉降量估算
利用上述步骤5)确定的适用于研究区域的低相干性区域地面沉降量最优插值方法,对上述步骤3)得到研究区域研究期间内不同时间段的地面沉降量进行空间插值,获得低相干区域研究期间内不同时间段的地面沉降量。
本发明上述技术方案通过地球信息空间技术,挖掘了克里金插值方法在低相干性区域地面沉降量估算中的潜力,构建克里金插值方法估算低相干性区域地面沉降量估算中的方法体系,并应用此方法对研究区低相干性区域不同时间段的地面沉降量进行了空间反演,克服了传统的InSAR技术估算低相干性区域地面沉降量的局限性。
本方法填补了估算低相干性区域地面沉降量的空白,为及时准确估算低相干性区域范围地面沉降量提供了一种新方法,实现了基于地理信息空间技术的低相干性区域地面沉降量的快速准确监测。
附图说明
图1是本发明实施案例所述的一种低相干性区域地面沉降量估算方法中计算的2017年5月30日轩岗采煤沉陷区地面沉降量空间分布图。
图2是本发明实施案例所述的一种低相干性区域地面沉降量估算方法中对2017年5月30日轩岗采煤沉陷区地面沉降量进行插值后的地面沉陷量空间分布图。
具体实施方式
下面通过对轩岗采煤沉陷区低相干性区域地面沉降量估算方法的应用案例,对本发明的具体实施方式作进一步说明。
步骤一、轩岗采煤沉陷区遥感影像订购
从专业机构遥感影像订购网站,订购下载覆盖轩岗采煤沉陷区2015年11月—2017年6月期间内的26景Sentinel-1A遥感影像。将Sentinel-1A数据转换为SARScape内部格式,并使用轩岗采煤沉陷区矢量边界进行数据裁剪,得到轩岗采煤沉陷区26景多时相SLC数据。
步骤二、轩岗采煤沉陷区遥感影像处理
对轩岗采煤沉陷区多时相SLC数据,进行基线估算及小基线集生成、图像匹配、干涉图生成、去平地效应、干涉图滤波、相干系数计算、相位解缠、轨道精炼和重去平等处理。其中,将相位解缠相干系数阈值设定为0.35。相干系数小于0.35的区域,即为低相干性区域。
步骤三、轩岗采煤沉陷区地面沉降量反演
对上述步骤二得到的轩岗采煤沉陷区遥感影像,进行轨道精炼和重去平、形变速率和高程系数估算、地形残余相位和大气相位去除以及地理编码等,反演得到轩岗采煤沉陷区研究期间内不同时间段的地面沉降量。图1为2017年5月30日轩岗采煤沉陷区地面沉降量空间分布图(空缺部分为低相干性区域)。
步骤四、轩岗采煤沉陷区地面沉降量数据检验
对上述步骤三得到的轩岗采煤沉陷区研究期间内不同时间段的地面沉降量数据进行正态检验,检验表明,不同时间段的地面沉降量均呈正态分布,具有自相关性,可以使用克里金法进行插值。
步骤五、轩岗采煤沉陷区低相干性区域地面沉降量插值方法确定
对上述步骤四得到的轩岗采煤沉陷区2017年5月30日地面沉降量,进行普通克里金插值与经验贝叶斯克里金插值,估算低相干性区域的地面沉降量。计算两种插值方法得到的低相干性区域边界地面沉降量与其相邻的高相干区域地面沉降量的均方根误差。结果表明,经验贝叶斯克里金插值方法的均方根误差最小,进而将经验贝叶斯克里金法确定适用于轩岗采煤沉陷区低相干性区域地面沉降量的最优插值方法。
步骤六、轩岗采煤沉陷区低相干性区域地面沉降量估算
利用上述步骤五确定的适用于轩岗采煤沉陷区低相干性区域地面沉降量的经验贝叶斯克里金插值方法,对上述步骤三得到的轩岗采煤沉陷区研究期间内不同时间段的地面沉降量进行空间插值,得到轩岗采煤沉陷区低相干区域研究期间内不同时间段的地面沉降量。图2为对2017年5月30日轩岗采煤沉陷区地面沉降量进行插值后的地面沉陷量空间分布图。
Claims (1)
1.一种低相干性区域地面沉降量估算方法,其特征在于:所述估算方法是按下列步骤进行的:
1)研究区遥感影像订购
从专业机构遥感影像订购网站,订购下载研究区域研究期间内多期Sentinel-1A遥感影像,将Sentinel-1A数据转换为SARScape内部格式,并使用研究区域范围矢量边界进行数据裁剪,得到研究区域多时相SLC数据;
2)研究区遥感影像处理
对研究区域多时相SLC数据,进行基线估算及小基线集生成、图像匹配、干涉图生成、去平地效应、干涉图滤波、相干系数计算、相位解缠、轨道精炼和重去平等处理,其中将相位解缠相干系数阈值设定为0.35,相干系数小于0.35的区域,即为低相干性区域;
3)研究区地面沉降量反演
对上述步骤2)处理得到的研究区域遥感影像,进行轨道精炼和重去平、形变速率和高程系数估算、地形残余相位和大气相位去除以及地理编码,反演得到研究区域研究期间内不同时间段的地面沉降量;
4)研究区地面沉降量数据检验
对上述步骤3)中反演得到的研究区域研究期间内不同时间段的地面沉降量数据进行正态检验,检验其是否呈正态分布,若呈正态分布,则具有自相关性,使用克里金法进行插值,否则对地面沉降量数据进行变换,使其呈正态分布;
5)研究区低相干性区域地面沉降量插值方法确定
对上述步骤4)处理得到的研究区域某一时段地面沉降量,进行普通克里金插值与经验贝叶斯克里金插值,得到研究区域低相干性区域的地面沉降量,计算两种插值方法得到的低相干性区域边界地面沉降量与其相邻的高相干区域地面沉降量的均方根误差,均方根误差越小,插值精度越高,进而确定适用于研究区低相干性区域地面沉降量的克里金插值方法;
6)研究区低相干性区域地面沉降量估算
利用上述步骤5)确定的适用于研究区域的低相干性区域地面沉降量最优插值方法,对上述步骤3)得到研究区域研究期间内不同时间段的地面沉降量进行空间插值,得到低相干区域研究期间内不同时间段的地面沉降量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910112992.XA CN109696152B (zh) | 2019-02-13 | 2019-02-13 | 一种低相干性区域地面沉降量估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910112992.XA CN109696152B (zh) | 2019-02-13 | 2019-02-13 | 一种低相干性区域地面沉降量估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109696152A CN109696152A (zh) | 2019-04-30 |
CN109696152B true CN109696152B (zh) | 2021-06-15 |
Family
ID=66234837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910112992.XA Active CN109696152B (zh) | 2019-02-13 | 2019-02-13 | 一种低相干性区域地面沉降量估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109696152B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110058234A (zh) * | 2019-05-20 | 2019-07-26 | 太原理工大学 | 一种解算矿区地表沉降三维形变的方法 |
CN111562557A (zh) * | 2020-05-12 | 2020-08-21 | 太原理工大学 | 一种黄土高原轩岗矿区地面沉降遥感预测方法 |
CN111856463B (zh) * | 2020-07-30 | 2024-03-15 | 首都师范大学 | 一种合成孔径雷达差分干涉测量最大有效干涉基线计算方法 |
CN113008202B (zh) * | 2021-03-31 | 2023-02-21 | 中铁上海设计院集团有限公司 | 一种融合不同合成孔径雷达干涉测量的地面沉降监测方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5847816A (en) * | 1997-01-14 | 1998-12-08 | Mcdonnell Douglas Corporation | Fiber optic micro-doppler ladar system and operating method therefor |
CN102927934A (zh) * | 2012-11-07 | 2013-02-13 | 中南大学 | 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法 |
CN104111456A (zh) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | 一种高速铁路沿线地表形变高分辨率InSAR监测方法 |
CN104990535A (zh) * | 2015-08-04 | 2015-10-21 | 天津市测绘院 | 水准/InSAR一体化监测装置 |
CN205079759U (zh) * | 2015-10-26 | 2016-03-09 | 江苏省电力公司南京供电公司 | 一种电缆户外终端支架防沉降监测装置 |
CN105444730A (zh) * | 2015-11-12 | 2016-03-30 | 中国矿业大学 | 多源数据监测矿区形变的时空特性及越界开采识别方法 |
CN106950556A (zh) * | 2017-05-03 | 2017-07-14 | 三亚中科遥感研究所 | 基于分布式散射体时序干涉sar技术的遗产地形变监测方法 |
CN107918127A (zh) * | 2017-11-20 | 2018-04-17 | 武汉大学 | 一种基于车载InSAR的道路边坡形变检测***及方法 |
CN108007401A (zh) * | 2017-11-20 | 2018-05-08 | 贵州省水利水电勘测设计研究院 | 一种基于船载InSAR平台的河湖库沿岸形变检测装置及方法 |
CN108896009A (zh) * | 2018-04-29 | 2018-11-27 | 天津大学 | 一种基于有效相干性的大型线状结构体沉降监测方法 |
CN109000613A (zh) * | 2018-08-22 | 2018-12-14 | 太原理工大学 | 一种采空区地质沉降三维快速检测*** |
CN109029370A (zh) * | 2018-06-09 | 2018-12-18 | 太原理工大学 | 一种量化区域草地放牧强度的遥感方法 |
CN109238227A (zh) * | 2018-10-31 | 2019-01-18 | 首都师范大学 | 一种表征地面沉降时空演变的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104122553B (zh) * | 2014-07-23 | 2017-01-25 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN106772377A (zh) * | 2017-01-18 | 2017-05-31 | 深圳市路桥建设集团有限公司 | 一种基于InSAR的建筑物变形监测方法 |
CN107515397A (zh) * | 2017-07-17 | 2017-12-26 | 中国南方电网有限责任公司超高压输电公司大理局 | 基于InSAR技术高频率地震地区换流站边坡沉降监测方法 |
CN107389029B (zh) * | 2017-08-24 | 2019-10-29 | 北京市水文地质工程地质大队 | 一种基于多源监测技术融合的地面沉降集成监测方法 |
CN108196252A (zh) * | 2017-12-05 | 2018-06-22 | 中国资源卫星应用中心 | 一种矿沉陷区沉降动态监测方法 |
CN108801211A (zh) * | 2018-04-08 | 2018-11-13 | 山东农业大学 | 一种采煤沉陷区insar影像快速识别方法 |
-
2019
- 2019-02-13 CN CN201910112992.XA patent/CN109696152B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5847816A (en) * | 1997-01-14 | 1998-12-08 | Mcdonnell Douglas Corporation | Fiber optic micro-doppler ladar system and operating method therefor |
CN102927934A (zh) * | 2012-11-07 | 2013-02-13 | 中南大学 | 一种利用单个InSAR干涉对获取矿区地表三维形变场的方法 |
CN104111456A (zh) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | 一种高速铁路沿线地表形变高分辨率InSAR监测方法 |
CN104990535A (zh) * | 2015-08-04 | 2015-10-21 | 天津市测绘院 | 水准/InSAR一体化监测装置 |
CN205079759U (zh) * | 2015-10-26 | 2016-03-09 | 江苏省电力公司南京供电公司 | 一种电缆户外终端支架防沉降监测装置 |
CN105444730A (zh) * | 2015-11-12 | 2016-03-30 | 中国矿业大学 | 多源数据监测矿区形变的时空特性及越界开采识别方法 |
CN106950556A (zh) * | 2017-05-03 | 2017-07-14 | 三亚中科遥感研究所 | 基于分布式散射体时序干涉sar技术的遗产地形变监测方法 |
CN107918127A (zh) * | 2017-11-20 | 2018-04-17 | 武汉大学 | 一种基于车载InSAR的道路边坡形变检测***及方法 |
CN108007401A (zh) * | 2017-11-20 | 2018-05-08 | 贵州省水利水电勘测设计研究院 | 一种基于船载InSAR平台的河湖库沿岸形变检测装置及方法 |
CN108896009A (zh) * | 2018-04-29 | 2018-11-27 | 天津大学 | 一种基于有效相干性的大型线状结构体沉降监测方法 |
CN109029370A (zh) * | 2018-06-09 | 2018-12-18 | 太原理工大学 | 一种量化区域草地放牧强度的遥感方法 |
CN109000613A (zh) * | 2018-08-22 | 2018-12-14 | 太原理工大学 | 一种采空区地质沉降三维快速检测*** |
CN109238227A (zh) * | 2018-10-31 | 2019-01-18 | 首都师范大学 | 一种表征地面沉降时空演变的方法 |
Non-Patent Citations (3)
Title |
---|
PS选择策略及在京津高铁沉降监测中的应用研究;丁尚起 等;《城市勘探》;20131231(第6期);第147-149,155页 * |
基于D-InSAR方法监测羊八井—当雄区域活动断层;李亚平 等;《黑龙江科技信息 科技论坛》;20181231;第61页 * |
高分辨SAR相干系数图像的开采沉陷区提取;马威 等;《测绘科学》;20180531;第43卷(第5期);第115-120页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109696152A (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109696152B (zh) | 一种低相干性区域地面沉降量估算方法 | |
US11269072B2 (en) | Land deformation measurement | |
CN106772342B (zh) | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 | |
Pepe et al. | A review of interferometric synthetic aperture RADAR (InSAR) multi-track approaches for the retrieval of Earth’s surface displacements | |
Trouve et al. | Improving phase unwrapping techniques by the use of local frequency estimates | |
Wang et al. | Improved SAR amplitude image offset measurements for deriving three-dimensional coseismic displacements | |
Kleinherenbrink et al. | A comparison of methods to estimate vertical land motion trends from GNSS and altimetry at tide gauge stations | |
Pawluszek-Filipiak et al. | Monitoring mining-induced subsidence by integrating differential radar interferometry and persistent scatterer techniques | |
CN113340191B (zh) | 时间序列干涉sar的形变量测量方法及sar*** | |
CN112986949B (zh) | 针对角反射器的sar高精度时序形变监测方法和装置 | |
Peltier et al. | Assessing the reliability and consistency of InSAR and GNSS data for retrieving 3D-displacement rapid changes, the example of the 2015 Piton de la Fournaise eruptions | |
CN112051572A (zh) | 一种融合多源sar数据三维地表形变监测方法 | |
WO2009055320A2 (en) | Method and apparatus for determining parameters for a parametric expression characterizing the phase of an acquired signal | |
Yun et al. | Mitigating atmospheric effects in InSAR measurements through high-resolution data assimilation and numerical simulations with a weather prediction model | |
KR20120009186A (ko) | Sar데이터를 이용하여 수치고도모형을 제작하기 위한 방법 | |
Feng et al. | Automatic selection of permanent scatterers-based GCPs for refinement and reflattening in InSAR DEM generation | |
Mao et al. | An InSAR ionospheric correction method based on variance component estimation with integration of MAI and RSS measurements | |
CN116148855B (zh) | 时序InSAR大气相位去除和形变解算的方法及*** | |
Yu et al. | Digital Elevation Model generation using ascending and Descending multi-baseline ALOS/PALSAR radar images | |
Zhang | Temporarily coherent point SAR interferometry | |
Du et al. | Patch-Like Reduction (PLR): A SAR Offset Tracking amplitude filter for deformation monitoring | |
Haque | 3D surface deformation model from sub-pixel correlation of optical imagery and InSAR | |
Abileah et al. | Bathymetry from fusion of multi-temporal Landsat and radar altimetery | |
Chaubey et al. | Improvement in InSAR phase unwrapping using external DEM | |
He et al. | Stereo radargrammetry in south-east Asia using TerraSAR-X stripmap data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |