CN115494500A - 基于遥感干涉测量的采空区快速探测方法、系统及应用 - Google Patents
基于遥感干涉测量的采空区快速探测方法、系统及应用 Download PDFInfo
- Publication number
- CN115494500A CN115494500A CN202211177774.2A CN202211177774A CN115494500A CN 115494500 A CN115494500 A CN 115494500A CN 202211177774 A CN202211177774 A CN 202211177774A CN 115494500 A CN115494500 A CN 115494500A
- Authority
- CN
- China
- Prior art keywords
- deformation
- phase
- image
- filtering
- interference
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 37
- 238000005305 interferometry Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 69
- 238000007670 refining Methods 0.000 claims abstract description 37
- 230000000694 effects Effects 0.000 claims abstract description 28
- 238000000605 extraction Methods 0.000 claims abstract description 12
- 238000007781 pre-processing Methods 0.000 claims abstract description 12
- 230000009466 transformation Effects 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003384 imaging method Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012732 spatial analysis Methods 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 238000012952 Resampling Methods 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000013138 pruning Methods 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 20
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000002265 prevention Effects 0.000 abstract description 2
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011045 prefiltration Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01C—CONSTRUCTION OF, OR SURFACES FOR, ROADS, SPORTS GROUNDS, OR THE LIKE; MACHINES OR AUXILIARY TOOLS FOR CONSTRUCTION OR REPAIR
- E01C23/00—Auxiliary devices or arrangements for constructing, repairing, reconditioning, or taking-up road or like surfaces
- E01C23/01—Devices or auxiliary means for setting-out or checking the configuration of new surfacing, e.g. templates, screed or reference line supports; Applications of apparatus for measuring, indicating, or recording the surface configuration of existing surfacing, e.g. profilographs
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Image Analysis (AREA)
Abstract
本发明属于遥感与公路防灾减灾数据识别技术领域,公开了基于遥感干涉测量的采空区快速探测方法、系统及应用。所述基于遥感干涉测量的采空区快速探测方法包括:通过影像预处理、基线估算、主辅影像配准、干涉图生成、去平地效应、干涉滤波、相位解缠、轨道精炼、相位转形变及地理编码、公路沿线形变提取,提取公路沿线形变图,根据所述公路沿线形变图对公路沿线进行沉降探测,判别公路发生沉陷位置是否位于采空区。与现有技术相比本发明利用遥感技术对公路沿线进行沉降监测具有实现简单、快捷、费用低、监测精度高、监测范围大、自动化程度高等优点。本发明借助遥感卫星,通过重复对目标区域的测量,能够快速获取到公路沿线微小的形变。
Description
技术领域
本发明属于遥感与公路防灾减灾数据识别技术领域,尤其涉及基于遥感干涉测量的采空区快速探测方法、系统及应用。
背景技术
通过上述分析,现有技术存在的问题及缺陷为:
(1)传统的监测方法的监测其覆盖范围十分有限。如果需要实时的监测数据,则需要在道路周围布设监测设备,这样就需要大量的设备成本和维护成本。
(2)现有技术对公路沿线进行沉降监测程序复杂、运行时间长、监测精度低、自动化程度低,而且不能快速获取到公路沿线微小的形变。
发明内容
为克服相关技术中存在的问题,本发明公开实施例提供了基于遥感干涉测量的采空区快速探测方法、系统及应用。
所述技术方案如下:基于遥感干涉测量的采空区快速探测方法包括:
通过影像预处理、基线估算、主辅影像配准、干涉图生成、去平地效应、干涉滤波、相位解缠、轨道精炼、相位转形变及地理编码、公路沿线形变提取,提取公路沿线形变图;
根据所述公路沿线形变图对公路沿线进行沉降探测,判别公路发生沉陷位置是否位于采空区。
所述基线估算包括:根据遥感影像轨道参数,计算两幅影像的空间基线,判断其是否小于临界基线。
在一个实施例中,所述主辅影像配准包括:
(1)以主影像为参考,根据主辅影像每个像元之间的坐标关系,通过相干系数法对辅影像的像元进行方位向和距离向重采样,并对辅影像的坐标进行变换,使得辅影像与主影像的像元对齐;
(2)进行粗配准与精配准,粗配准利用卫星的轨道数据选取大于等于3个特征点计算辅影像像元相对于主影像像元的偏移量;精配准在粗配准的基础上按照相干系数法提取主辅影像的同名点,计算上述同名点的偏差,并表示为主影像坐标的函数关系式,利用函数关系式对辅影像进行重采样;该主影像坐标的函数关系式表示为:
在一个实施例中,所述干涉图生成包括:形变前后的两幅影像经过精配准,配准精度达到1/10个像元;然后进行两幅影像的共轭相乘,得到条纹相间的干涉图纹理;对获得的干涉图进行相干性计算,得到相干性系数图,该相干性系数图用于表示干涉像对的干涉质量。
在一个实施例中,所述去平地效应包括:采用参考椭球和基线模型模拟平地相位,去除平地效应;
在一个实施例中,所述轨道精炼包括:对解缠后的干涉图进行轨道精炼处理,采用多项式精细化方法通过选取地面控制点进行轨道精炼;二次多项式为:
所述相位转形变及地理编码包括:将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下;
干涉相位是两次雷达成像时天线分别到监测的地面目标点的斜距的函数,在这两次成像期间如果地面目标点发生了形变,相位就会有形变相位的作用,电磁波在大气中传播会有大气延迟的作用,此时的干涉相位由地形相、平地相位、地表形变相位、大气延迟相位和噪声构成,表达式为:
表示第一次成像时雷达天线到地面观测点的斜距;表示卫星第一次成像时的高程;表示雷达两次成像时多普勒中心连线的距离,称为空间基线;是空间基线与水平方向线之间的倾角;表示第一次观测目标点的侧视角;为地面点的高程;
本发明的另一目的在于提供一种实现所述基于遥感干涉测量的采空区快速探测方法的系统,该基于遥感干涉测量的采空区快速探测系统包括:
基线估算模块,用于根据遥感影像轨道参数,计算两幅影像的空间基线,判断其是否小于临界基线;
主辅影像配准模块,用于先进行粗配准再进行精配准;
干涉图生成模块,用于形变前后的两幅影像经过精配准,达到配准的精度要求进行两幅影像的共轭相乘,得到条纹相间的干涉图纹理;对获得的干涉图进行相干性计算,得到相干性系数图,该相干性系数图用于表示干涉像对的干涉质量;
去平地效应模块,用于采用参考椭球和基线模型模拟平地相位,去除平地效应;
轨道精炼模块,用于对解缠后的干涉图进行轨道精炼处理,采用多项式精细化方法通过选取地面控制点进行轨道精炼;
相位转形变及地理编码模块,用于将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下;
本发明的另一目的在于提供一种计算机设备,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行所述的基于遥感干涉测量的采空区快速探测方法。
本发明的另一目的在于提供一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,使得所述处理器执行所述的基于遥感干涉测量的采空区快速探测方法。
结合上述的所有技术方案,本发明所具备的优点及积极效果为:
第一、针对上述现有技术存在的技术问题以及解决该问题的难度,紧密结合本发明的所要保护的技术方案以及研发过程中结果和数据等,详细、深刻地分析本发明技术方案如何解决的技术问题,解决问题之后带来的一些具备创造性的技术效果。通过裁剪原始影像,提高了数据处理的整体速率而且避免了大面积处理引入的误差;主辅影像的配准达到了十分之一个像元,大大提高了干涉图的质量;通过干涉滤波,进一步减少数据处理中的误差。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本公开的实施例,并与说明书一起用于解释本公开的原理。
图1是本发明实施例提供的基于遥感干涉测量的采空区快速探测方法流程图;
图2是本发明实施例提供的基于遥感干涉测量的采空区快速探测系统示意图;
图3是本发明实施例提供的遥感干涉测量模型示意图;
图4是本发明实施例提供的基于遥感干涉测量的采空区快速探测方法原理图;
图5是本发明实施例提供的裁剪后的遥感影像图;
图6是本发明实施例提供的干涉图;
图7是本发明实施例提供的去平后的干涉图;
图8是本发明实施例提供的滤波后的干涉图;
图9是本发明实施例提供的相位解缠图;
图10是本发明实施例提供的沉降监测形变图;
图11是本发明实施例提供的公路沿线沉降监测形变图;
图中:1、影像预处理模块;2、基线估算模块;3、主辅影像配准模块;4、干涉图生成模块;5、去平地效应模块;6、干涉滤波模块;7、相位解缠模块;8、轨道精炼模块;9、相位转形变及地理编码模块;10、公路沿线形变提取模块。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其他方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施的限制。
一、解释说明实施例:
如何快速对公路沿线进行沉降观测,判别公路发生沉陷灾害的可能性是本发明所要解决的技术问题,为此,本发明的目的是通过卫星SAR影像预处理、基线估算、主辅影像配准、干涉图生成、去平地效应、干涉滤波、相位解缠、轨道精炼、相位转形变及地理编码、公路沿线形变提取,提供了一种基于遥感干涉测量的采空区快速探测方法,具体涉及一种基于遥感干涉测量技术的公路沿线采空区快速探测技术。
实施例1
如图1所示,本发明实施例提供的基于遥感干涉测量的采空区快速探测方法包括以下步骤:
干涉像对的基线,特别是空间垂直基线对平地相位的去除起到至关重要的作用,空间基线过大容易造成干涉图失相干。根据遥感影像轨道参数计算两幅影像的空间基线,判断其是否小于临界基线,只有小于临界基线才可以进行下一步工作。
先进行粗配准再进行精配准;要想利用遥感影像进行干涉处理获取地面形变,首要的步骤就是实现主辅影像的精配准。主辅影像的配准精度直接决定能否生成干涉图纹,配准有误差就会产生噪声,影响干涉图的相干性,当配准误差达到一个像素时,就会造成干涉图完全失相干。主辅影像的配准就是以主影像为参考,根据主辅影像每个像元之间的坐标关系通过相干系数法,对辅影像的像元进行方位向和距离向进行重采样,并对辅影像的坐标进行变换,使得辅影像与主影像的像元对齐。一般先进行粗配准再进行精配准。粗配准利用卫星的轨道数据来选取少许(大于等于3)的特征点来计算辅影像像元相对于主影像像元的偏移量。精配准在粗配准的基础之上采用相干系数法提取主辅影像的同名点,计算这些同名点的偏差,并表示为主影像坐标的函数关系式,利用这一关系对辅影像进行重采样。一般要求精配准的误差不能超过1/10个像元。
在参考椭球面上高程不变的平地引起的干涉图纹在距离向和方位向呈现重复性变化的现象称为平地效应。由于平地效应的影响,使得干涉图纹变得密集而且覆盖了地形变化引起的相位。必须去除平地相位才能得到精确的地形信息,一般采用参考椭球和基线模型模拟平地相位,得到去除的目的。
由于图像本生的斑点噪声以及受到失相关的影响,会使得干涉图不清晰。噪声的来源和形成机理不同一般消除噪声的方法也不一样。总体分为前置滤波和后置滤波两部分,前置滤波指对干涉前的影像进行滤波处理,而后置滤波指的是对干涉图进行滤波处理。前置滤波又分为距离向滤波和方位向滤波。后置滤波有均值滤波、中值滤波和方法等。
相位解缠是整个干涉处理中重要的一环,解缠结果的好坏直接影像最后形变结果的精度。相位值是通过主辅影像共轭相乘得到的,其干涉相位丢失了的整数倍,造成了相位的解缠,因此只有还原这些解缠的相位,才能更好的反演地表的形变。目前比较经典的解缠算法有枝切树法、最小二乘法、网络流法及最小费用流法等。
为了消除干涉图中可能存在的轨道误差,需要对解缠后的干涉图进行轨道精炼处理。采用多项式精细化方法通过选取地面控制点进行轨道精炼。
用于将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下。
实施例2
所述主辅影像配准包括:
(1)主辅影像的配准:以主影像为参考,根据主辅影像每个像元之间的坐标关系通过相干系数法对辅影像的像元进行方位向和距离向进行重采样,并对辅影像的坐标进行变换,使得辅影像与主影像的像元对齐;
(2)进行粗配准与精配准,粗配准利用卫星的轨道数据选取少许的特征点计算辅影像像元相对于主影像像元的偏移量;精配准在粗配准的基础上采用相干系数法提取主辅影像的同名点,计算上述同名点的偏差,并表示为主影像坐标的函数关系式,如下所示,利用所述函数关系式对辅影像进行从采样。
实施例3
所述相位转形变及地理编码包括:将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下。
实施例4
基于本发明实施例提供的基于遥感干涉测量的采空区快速探测方法,进一步地,在步骤中,干涉相位是两次雷达成像时天线分别到监测的地面目标点的斜距的函数,在这两次成像期间如果地面目标点发生了形变干涉相位就会有形变相位的作用,电磁波在大气中传播会有大气延迟的作用。此时的干涉相位是由地形相、平地相位、地表形变相位、大气延迟相位和噪声构成:
表示第一次成像时雷达天线到地面观测点的斜距;表示卫星第一次成像时的高程;表示雷达两次成像时多普勒中心连线的距离,称为空间基线;是空间基线与水平方向线之间的倾角;表示第一次观测目标点的侧视角;为地面点的高程;
实施例5
如图2所示,本发明实施例提供的基于遥感干涉测量的采空区快速探测系统包括:
影像预处理模块1:用于根据公路范围裁剪遥感影像到合适的大小。
基线估算模块2:用于根据遥感影像轨道参数计算两幅影像的空间基线,判断其是否小于临界基线,只有小于临界基线才可以进行下一步工作。
主辅影像配准模块3,用于先进行粗配准再进行精配准。精配准在粗配准的基础之上采用一定的算法提取主辅影像的同名点,计算这些同名点的偏差,并表示为主影像坐标的函数关系式,利用这一关系对辅影像进行从采样。
干涉图生成模块4,用于形变前后的两幅影像经过精配准,达到配准的精度要求就进行两幅影像的共轭相乘,得到条纹相间的干涉图纹理。对获得的干涉图进行相干性计算,得到相干性系数图,该相干性系数图用于表示干涉相对的干涉质量。
去平地效应模块5,用于采用参考椭球和基线模型模拟平地相位去除平地效应。
轨道精炼模块8,用于消除干涉图中可能存在的轨道误差,对解缠后的干涉图进行轨道精炼处理。采用多项式精细化方法通过选取地面控制点进行轨道精炼。
相位转形变及地理编码模块9,用于将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
上述装置/单元之间的信息交互、执行过程等内容,由于与本发明方法实施例基于同一构思,其具体功能及带来的技术效果,具体可参见方法实施例部分,此处不再赘述。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
二、应用实施例:
应用例1
本发明应用实施例还提供了一种计算机设备,该计算机设备包括:至少一个处理器、存储器以及存储在所述存储器中并可在所述至少一个处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任意各个方法实施例中的步骤。
应用例2
本发明应用实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时可实现上述各个方法实施例中的步骤。
应用例3
本发明应用实施例还提供了一种信息数据处理终端,所述信息数据处理终端用于实现于电子装置上执行时,提供用户输入接口以实施如上述各方法实施例中的步骤,所述信息数据处理终端不限于手机、电脑、交换机。
应用例4
本发明应用实施例还提供了一种服务器,所述服务器用于实现于电子装置上执行时,提供用户输入接口以实施如上述各方法实施例中的步骤。
应用例5
本发明应用实施例提供了一种计算机程序产品,当计算机程序产品在电子设备上运行时,使得电子设备执行时可实现上述各个方法实施例中的步骤。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质至少可以包括:能够将计算机程序代码携带到拍照装置/终端设备的任何实体或装置、记录介质、计算机存储器、只读存储器、随机存取存储器、电载波信号、电信信号以及软件分发介质。例如U盘、移动硬盘、磁碟或者光盘等。
三、实施例相关效果的实验证据:
如图4所示,本发明实施例提供的基于遥感干涉测量的采空区快速探测方法包括以下步骤:
步骤2:干涉像对的基线,特别是空间垂直基线对平地相位的去除起到至关重要的作用。本实验基于轨道参数方法对干涉像对进行基线估算,估算结果如下:
表1基线估算结果
步骤3:主辅影像的精度配准直接决定能否生成干涉图纹,配准有误差就会产生噪声,影响干涉图的相干性,当配准误差达到一个像素时,就会造成干涉图完全失相干。主辅影像的配准就是以主影像为参考,根据主辅影像每个像元之间的坐标关系通过相干系数法,对辅影像的像元进行方位向和距离向进行重采样,并对辅影像的坐标进行变换,使得辅影像与主影像的像元对齐。本实验精配准采用交叉相关性方法对数据进行处理配准结果误差在1/10个像元之内。
步骤4:主辅影像经过配准后进行共轭相乘生成差分干涉图,干涉图如图6。
步骤6:去平地相位后的差分干涉图存在有相位噪声,而滤波的目的就是减少相位噪声,本实验采用的滤波方法是—“枝切”—算法。这种滤波方法运用可变的滤波器,可以提高干涉条纹的清晰度,减少干涉相对失相干造成的噪声,滤波后干涉图如图8。
步骤7:相位解缠是整个干涉处理中重要的一环,解缠结果的好坏直接影像最后形变结果的精度。相位值是通过主辅影像共轭相乘得到的,其干涉相位丢失了的整数倍,造成了相位的缠绕,因此只有还原这些缠绕的相位,才能更好的反演地表的形变。本实验利用最小费用流法,相干性阈值设置为0.35,对滤波后的干涉相位进行解缠,这种方法掩膜处理相干性小于阈值的像元,对于有大面积的低相干区域时,可以取得更好的解缠效果。相位解缠图如图9。
步骤8:为了消除干涉图中可能存在的轨道误差,需要对解缠后的干涉图进行轨道精炼处理。本实验采用多项式精细化方法通过选取地面控制点进行轨道精炼。本发明沿某市主城区道路选取干涉图中相干性高的点进行轨道精炼,选择轨道精炼模型为二次多项式,表达式为:
轨道精炼结果见表2所示。
表2轨道精炼结果
步骤9:据干涉相位与形变相位之间的关系估算形变相位,此时的坐标系还是斜距坐标系。要利用轨道数据以及控制点数据等对这一坐标系进行地理编码,即实现斜距坐标系到地理坐标系的转换。经过地理编码后得到各期形变结果,见图10所示。
以上所述,仅为本发明较优的具体的实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种基于遥感干涉测量的采空区快速探测方法,其特征在于,该方法包括:
通过影像预处理、基线估算、主辅影像配准、干涉图生成、去平地效应、干涉滤波、相位解缠、轨道精炼、相位转形变及地理编码、公路沿线形变提取,提取公路沿线形变图;
根据所述公路沿线形变图对公路沿线进行沉降探测,判别公路发生沉陷位置是否位于采空区。
3.根据权利要求1所述的基于遥感干涉测量的采空区快速探测方法,其特征在于,所述主辅影像配准包括:
(1)以主影像为参考,根据主辅影像每个像元之间的坐标关系,通过相干系数法对辅影像的像元进行方位向和距离向重采样,并对辅影像的坐标进行变换,使得辅影像与主影像的像元对齐;
(2)进行粗配准与精配准,粗配准利用卫星的轨道数据选取大于等于3个特征点计算辅影像像元相对于主影像像元的偏移量;精配准在粗配准的基础上按照相干系数法提取主辅影像的同名点,计算上述同名点的偏差,并表示为主影像坐标的函数关系式,利用函数关系式对辅影像进行重采样;该主影像坐标的函数关系式表示为:
7.根据权利要求1所述的基于遥感干涉测量的采空区快速探测方法,其特征在于,所述轨道精炼包括:对解缠后的干涉图进行轨道精炼处理,采用多项式精细化方法通过选取地面控制点进行轨道精炼;二次多项式为:
所述相位转形变及地理编码包括:将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下;
干涉相位是两次雷达成像时天线分别到监测的地面目标点的斜距的函数,在这两次成像期间如果地面目标点发生了形变,相位就会有形变相位的作用,电磁波在大气中传播会有大气延迟的作用,此时的干涉相位由地形相、平地相位、地表形变相位、大气延迟相位和噪声构成,表达式为:
表示第一次成像时雷达天线到地面观测点的斜距;表示卫星第一次成像时的高程;表示雷达两次成像时多普勒中心连线的距离,称为空间基线;是空间基线与水平方向线之间的倾角;表示第一次观测目标点的侧视角;为地面点的高程;
8.一种实现如权利要求1-7任意一项所述基于遥感干涉测量的采空区快速探测方法的系统,其特征在于,该基于遥感干涉测量的采空区快速探测系统包括:
基线估算模块(2),用于根据遥感影像轨道参数,计算两幅影像的空间基线,判断其是否小于临界基线;
主辅影像配准模块(3),用于先进行粗配准再进行精配准;
干涉图生成模块(4),用于形变前后的两幅影像经过精配准,达到配准的精度要求进行两幅影像的共轭相乘,得到条纹相间的干涉图纹理;对获得的干涉图进行相干性计算,得到相干性系数图,该相干性系数图用于表示干涉像对的干涉质量;
去平地效应模块(5),用于采用参考椭球和基线模型模拟平地相位,去除平地效应;
轨道精炼模块(8),用于对解缠后的干涉图进行轨道精炼处理,采用多项式精细化方法通过选取地面控制点进行轨道精炼;
相位转形变及地理编码模块(9),用于将经过滤波以及解缠后的形变相位转换为形变值;然后通过地理编码将SAR斜距坐标系下的形变数据转换到地理坐标系下;
9.一种计算机设备,其特征在于,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行权利要求1-7任意一项所述的基于遥感干涉测量的采空区快速探测方法。
10.一种计算机可读存储介质,存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,使得所述处理器执行权利要求1-7任意一项所述的基于遥感干涉测量的采空区快速探测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211177774.2A CN115494500A (zh) | 2022-09-27 | 2022-09-27 | 基于遥感干涉测量的采空区快速探测方法、系统及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211177774.2A CN115494500A (zh) | 2022-09-27 | 2022-09-27 | 基于遥感干涉测量的采空区快速探测方法、系统及应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115494500A true CN115494500A (zh) | 2022-12-20 |
Family
ID=84472866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211177774.2A Pending CN115494500A (zh) | 2022-09-27 | 2022-09-27 | 基于遥感干涉测量的采空区快速探测方法、系统及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115494500A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115993601A (zh) * | 2023-03-22 | 2023-04-21 | 四川省公路规划勘察设计研究院有限公司 | 一种强盐渍土区域公路变形的时序InSAR监测方法 |
-
2022
- 2022-09-27 CN CN202211177774.2A patent/CN115494500A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115993601A (zh) * | 2023-03-22 | 2023-04-21 | 四川省公路规划勘察设计研究院有限公司 | 一种强盐渍土区域公路变形的时序InSAR监测方法 |
CN115993601B (zh) * | 2023-03-22 | 2023-06-09 | 四川省公路规划勘察设计研究院有限公司 | 一种强盐渍土区域公路变形的时序InSAR监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107102333B (zh) | 一种星载InSAR长短基线融合解缠方法 | |
KR100529401B1 (ko) | 합성 개구 레이더 영상을 이용한 수치표고자료 제작 장치및 그 방법 | |
White et al. | Lidar-derived national shoreline: empirical and stochastic uncertainty analyses | |
CN109212522B (zh) | 一种基于双基星载sar的高精度dem反演方法及装置 | |
CN113340191A (zh) | 时间序列干涉sar的形变量测量方法及sar系统 | |
CN112711021B (zh) | 一种多分辨率InSAR交互干涉时序分析方法 | |
CN113311433B (zh) | 一种质量图和最小费用流结合的InSAR干涉相位两步解缠方法 | |
CN115201825B (zh) | 一种InSAR震间形变监测中的大气延迟校正方法 | |
CN109239710B (zh) | 雷达高程信息的获取方法及装置、计算机可读存储介质 | |
CN112882030B (zh) | InSAR成像干涉一体化处理方法 | |
CN111239736A (zh) | 基于单基线的地表高程校正方法、装置、设备及存储介质 | |
CN115494500A (zh) | 基于遥感干涉测量的采空区快速探测方法、系统及应用 | |
CN117523404B (zh) | 一种堆土场动态变化监测方法、装置、终端及存储介质 | |
CN110927765B (zh) | 激光雷达与卫星导航融合的目标在线定位方法 | |
CN113238228A (zh) | 一种基于水准约束的InSAR三维地表形变获取方法、系统及装置 | |
CN116740151A (zh) | 一种InSAR点云的配准方法及终端设备 | |
CN114280608B (zh) | 一种DInSAR高程相关大气效应去除方法及系统 | |
CN113311432B (zh) | 一种基于相位导数方差的InSAR长短基线融合相位估计方法 | |
CN115546264A (zh) | 一种星载InSAR图像精细配准与立体测量方法 | |
CN115113202A (zh) | 一种基于二维高斯模型的干涉相位迭代解缠方法 | |
CN110021041B (zh) | 基于双目相机的无人驾驶场景增量式网格化结构重建方法 | |
CN111580101A (zh) | 一种基于外部DEM的InSAR基线误差无控改正方法及装置 | |
Li et al. | Atmospheric phase delay correction of D-InSAR based on Sentinel-1A | |
CN117970332B (zh) | 一种联合光学影像的InSAR城区地表形变危险性精细识别方法 | |
KR102275168B1 (ko) | 비전 센서에 기반한 차량 내비게이션 방법 |
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 |