CN112986988B - 一种特殊地形条件下sar干涉测量形变量解析方法 - Google Patents
一种特殊地形条件下sar干涉测量形变量解析方法 Download PDFInfo
- Publication number
- CN112986988B CN112986988B CN202110095336.0A CN202110095336A CN112986988B CN 112986988 B CN112986988 B CN 112986988B CN 202110095336 A CN202110095336 A CN 202110095336A CN 112986988 B CN112986988 B CN 112986988B
- Authority
- CN
- China
- Prior art keywords
- phase
- deformation
- angle
- sight
- processing
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 10
- 238000005305 interferometry Methods 0.000 claims description 10
- 230000001154 acute effect Effects 0.000 claims description 2
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及遥感影像处理领域,公开了一种特殊地形条件下SAR干涉测量形变量解析方法,该方法步骤如下:步骤1:对前后时相的SAR影像进行配准处理;步骤2:对配准后的影像对进行差分干涉处理;步骤3:对差分干涉处理后得到的相位差进行相位处理;将差分相位转换为视线方向的形变量;步骤4:对坡向视线角α进行计算。坡向视线角定义为点位梯度方向与雷达视线方向的夹角;步骤5:根据坡向视线角及形变量进行判断解析;步骤6:得到形变量解析结果。本发明提供的解析方法能实现特殊地形条件下,SAR干涉测量形变量的正确解析,对城市地面沉降监测有重要意义。
Description
技术领域
本发明涉及遥感影像处理领域,具体说是一种用于特殊地形条件下SAR干涉测量形变量的解析方法。
背景技术
现有的SAR干涉测量形变量解析方法通常不考虑地面的地形条件。其通常采用前后时相的SAR影像进行配准、干涉处理、相位处理等步骤获取沿视线方向的点位形变量,然后根据其符号来判断点位的升降情况。但在特殊地形条件下,当地形梯度方向与雷达视线方向夹角为锐角时,传统的判断方法与真实情况呈现相悖的状态,使得传统方法应用局限极大。
发明内容
本发明所要解决的技术问题是:针对上述存在的问题,提供了一种用于特殊地形条件下的SAR干涉测量形变量解析方法。该方法在利用SAR干涉测量获得形变量的基础上,通过对梯度方向与雷达视线方向的位置关系的分析,实现点位沉降或抬升情况的判断。该发明克服了传统方法在特殊地形条件下不能正确判断点位升降情况的缺陷,具有更广泛的适用性。
本申请提供一种用于特殊地形条件下SAR干涉测量形变量的解析方法,采用的技术方案如下:
步骤1:对前后时相的SAR影像进行配准处理;
步骤2:对配准后得到的影像对进行差分干涉处理;差分干涉处理包括两个部分,首先对影像复数数据做干涉运算,然后将干涉相位与地形相位做差值运算;
步骤3:对差分干涉处理后得到的相位差进行相位处理,将差分相位转换为视线方向的形变量;
步骤4:对坡向视线角进行计算;其中,所述坡向视线角为点位梯度方向和雷达视线方向的夹角;
坡向视线角α的计算方式如下:
步骤5:根据坡向视线角及形变量进行判断解析;
当α≥90°,d<0或α<90°,d>0时,形变为沉降其中,d为形变量;当α≥90°,d>0或α<90°,d<0时,形变为抬升;当d=0时,不存在相对形变;
步骤6:得到形变量的分析结果。
进一步的,配准处理精度要达到亚像素级别,配准处理是进行干涉测量的基础,所以需确保其精度。
进一步的,所述步骤3中的相位处理包括:相位滤波、相位截仓和相高转换,其中,相位滤波采用自适应滤波方法,相位解缠使用最小费用流方法,通过相高转换后得到形变量。
与现有技术相比,采用上述技术方案的有益效果为:本发明能实现特殊地形条件下,SAR干涉测量形变量的正确解析,对城市地面沉降监测有重要意义。
附图说明
图1是本发明实施例提供的一种用于特殊地形条件下SAR干涉测量形变量的解析方法的流程示意图。
具体实施方式
下面结合附图对本发明做进一步描述。
实施例1
如图1所示,本发明实施例提供一种用于特殊地形条件下SAR干涉测量形变量的解析方法,该方法步骤包括以下步骤:
步骤1:对前后时相的SAR影像进行配准处理。
步骤2:对配准后的影像对进行差分干涉处理。
差分干涉处理包括两个部分,首先对影像复数数据做干涉运算,然后将干涉相位与地形相位做差值运算。
步骤3:对得到的相位差进行相位处理。包括相位滤波、相位解缠、相高转换等处理。目的是将差分相位转换为视线方向的形变量d。
步骤4:对坡向视线角α进行计算。坡向视线角定义为点位梯度方向与雷达视线方向的夹角。
步骤5:根据坡向视线角α及形变量d进行判断解析。
步骤6:得到形变量解析结果。
具体的,在本实施例中,坡向视线角α的计算方式如下:
具体的,在本实施例中,根据坡向视角和形变量判断解析具体如下:
当α≥90°,d<0或α<90°,d>0时,形变为沉降;当α≥90°,d>0或α<90°,d<0时,形变为抬升;当d=0时,不存在相对形变。据此可实现特殊地形条件下的形变量解析。
实施例2
本实施例与实施例1基本一致,优选地,在本实施例中,配准处理精度要达到亚像素级别,配准处理是进行干涉测量的基础,必须保证其精度。
实施例3
本实施例与实施例1或实例2基本一致,优选地,在本实施例中,其中干涉运算是幅度相乘和相位相减。此时相位相减的结果中不仅包含形变相位还包含地形相位。其中,地形相位由DEM数据反演得到。通过干涉相位与地形相位做差值运算后得到的差分相位保留的基本上是两时相影像地表形变导致的相位变化。
实施例4
本实施例与实施例1或实施例2或实施例3基本一致,优选地,在本实施例中,相位滤波采用自适应滤波方法,相位解缠使用最小费用流方法,然后通过相高转换后得到形变量d。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员,在不脱离本发明的精神所做的非实质性改变或改进,都应该属于本发明权利要求保护的范围。
Claims (5)
1.一种特殊地形条件下SAR干涉测量形变量解析方法,其特征在于,所述特殊地形条件下为在当地形梯度方向与雷达视线方向夹角为锐角时,包括:
步骤1:对前后时相的SAR影像进行配准处理;
步骤2:对配准后得到的影像对进行差分干涉处理;所述步骤2中的差分干涉处理包括:
先对影像的复数数据做干涉运算,然后将得到的干涉相位与地形相位做差值运算;所述干涉运算包括对影像的复数数据的幅度相乘和相位相减;
步骤3:对差分干涉处理后得到的相位差进行相位处理,将差分相位转换为视线方向的形变量;
步骤4:对坡向视线角进行计算;其中,所述坡向视线角为点位梯度方向和雷达视线方向的夹角;所述步骤4中对坡向视线角的计算方式如下:
步骤5:根据坡向视线角及形变量进行判断解析;所述步骤5中判断解析具体包括:
当α≥90°,d<0或α<90°,d>0时,形变为沉降,其中,d为形变量;
当α≥90°,d>0或α<90°,d<0时,形变为抬升;当d=0时,不存在相对形变;
步骤6:得到形变量的分析结果。
2.根据权利要求1所述的一种特殊地形条件下SAR干涉测量形变量解析方法,其特征在于,所述步骤1中配准处理的精度为亚像素级别。
3.根据权利要求1所述的一种特殊地形条件下SAR干涉测量形变量解析方法,其特征在于,所述步骤3中的相位处理包括:相位滤波、相位解缠和相高转换。
4.根据权利要求3所述的一种特殊地形条件下SAR干涉测量形变量解析方法,其特征在于,所述相位滤波采用自适应滤波方法。
5.根据权利要求3所述的一种特殊地形条件下SAR干涉测量形变量解析方法,其特征在于,所述相位解缠使用最小费用流方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110095336.0A CN112986988B (zh) | 2021-01-25 | 2021-01-25 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110095336.0A CN112986988B (zh) | 2021-01-25 | 2021-01-25 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112986988A CN112986988A (zh) | 2021-06-18 |
CN112986988B true CN112986988B (zh) | 2022-09-20 |
Family
ID=76344759
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110095336.0A Active CN112986988B (zh) | 2021-01-25 | 2021-01-25 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112986988B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113534154B (zh) * | 2021-09-16 | 2021-11-30 | 成都理工大学 | 一种sar视线向变形与坡度坡向敏感度计算方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7589662B1 (en) * | 2007-06-20 | 2009-09-15 | The United States Of America As Represented By The Secretary Of The Navy | Synthetic aperture radar systems and methods |
CN102508245A (zh) * | 2011-11-18 | 2012-06-20 | 北京航空航天大学 | 一种星载多频率与多基线InSAR高程估计精度等效性确定方法 |
CN102608576A (zh) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 |
CN103675790A (zh) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
CN106226764A (zh) * | 2016-07-29 | 2016-12-14 | 安徽理工大学 | 一种基于D‑InSAR的煤矿开采地沉陷区域的测定方法 |
CN107132539A (zh) * | 2017-05-03 | 2017-09-05 | 中国地质科学院探矿工艺研究所 | 一种基于小基线集的时间序列InSAR的滑坡早期识别方法 |
JP2017219483A (ja) * | 2016-06-09 | 2017-12-14 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
CN109001731A (zh) * | 2018-06-06 | 2018-12-14 | 中国电子科技集团公司第二十九研究所 | 一种sar干涉相位解缠参考点选取方法、设备及存储介质 |
CN112014841A (zh) * | 2020-08-31 | 2020-12-01 | 中国矿业大学 | 一种基于DS-InSAR技术监测油田区地表形变分析方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102645650A (zh) * | 2012-03-06 | 2012-08-22 | 北京北科安地科技发展有限公司 | 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术 |
KR101490981B1 (ko) * | 2012-12-28 | 2015-02-09 | 서울시립대학교 산학협력단 | 위성레이더 간섭도의 이온왜곡 보정방법 및 그 장치 |
JP6716435B2 (ja) * | 2016-11-21 | 2020-07-01 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
US11782141B2 (en) * | 2018-02-05 | 2023-10-10 | Centre Interdisciplinaire De Developpement En Cartographie Des Oceans (Cidco) | Method and apparatus for automatic calibration of mobile LiDAR systems |
CN112051571A (zh) * | 2020-08-14 | 2020-12-08 | 中冶(北京)交通科技发展有限公司 | 一种新型差分InSAR的LOS向形变量估计方法 |
-
2021
- 2021-01-25 CN CN202110095336.0A patent/CN112986988B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7589662B1 (en) * | 2007-06-20 | 2009-09-15 | The United States Of America As Represented By The Secretary Of The Navy | Synthetic aperture radar systems and methods |
CN102508245A (zh) * | 2011-11-18 | 2012-06-20 | 北京航空航天大学 | 一种星载多频率与多基线InSAR高程估计精度等效性确定方法 |
CN102608576A (zh) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 |
CN103675790A (zh) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
JP2017219483A (ja) * | 2016-06-09 | 2017-12-14 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
CN106226764A (zh) * | 2016-07-29 | 2016-12-14 | 安徽理工大学 | 一种基于D‑InSAR的煤矿开采地沉陷区域的测定方法 |
CN107132539A (zh) * | 2017-05-03 | 2017-09-05 | 中国地质科学院探矿工艺研究所 | 一种基于小基线集的时间序列InSAR的滑坡早期识别方法 |
CN109001731A (zh) * | 2018-06-06 | 2018-12-14 | 中国电子科技集团公司第二十九研究所 | 一种sar干涉相位解缠参考点选取方法、设备及存储介质 |
CN112014841A (zh) * | 2020-08-31 | 2020-12-01 | 中国矿业大学 | 一种基于DS-InSAR技术监测油田区地表形变分析方法 |
Non-Patent Citations (2)
Title |
---|
一种基于GPU的机载SAR数据高效处理方法;阎冬 等;《第六届高分辨率对地观测学术年会论文集(上)》;20190920;全文 * |
高山峡谷区滑坡灾害隐患InSAR早期识别——以雅砻江中段为例;戴可人 等;《雷达学报》;20200622;第9卷(第3期);第554-568页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112986988A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113340191B (zh) | 时间序列干涉sar的形变量测量方法及sar系统 | |
CN106651828B (zh) | 一种工业小尺度运动模糊成像条件下的产品尺寸亚像素测量方法 | |
CN103091676A (zh) | 矿区地表开采沉陷合成孔径雷达干涉测量的监测及解算方法 | |
CN101916373B (zh) | 基于小波检测和脊线跟踪的道路半自动提取方法 | |
CN112986988B (zh) | 一种特殊地形条件下sar干涉测量形变量解析方法 | |
CN108053416B (zh) | 一种基于单幅卫星图像的最大储油量提取系统 | |
CN102279973A (zh) | 基于高梯度关键点的海天线检测方法 | |
CN103400388A (zh) | 一种利用RANSAC消除Brisk关键点错误匹配点对的方法 | |
CN113484860B (zh) | 基于多普勒中心异常的sar图像海洋锋检测方法及系统 | |
CN112051572A (zh) | 一种融合多源sar数据三维地表形变监测方法 | |
CN108562900B (zh) | 一种基于高程校正的sar图像几何配准方法 | |
CN115755043B (zh) | 一种基于x波段非相参雷达的波浪场重构及预测方法 | |
CN116047519A (zh) | 一种基于合成孔径雷达干涉测量技术的选点方法 | |
CN108230365B (zh) | 基于多源差异图像内容融合的sar图像变化检测方法 | |
Zhang et al. | A study on coastline extraction and its trend based on remote sensing image data mining | |
CN111562575A (zh) | 一种用于地面沉降的监测方法 | |
CN111896953A (zh) | 基于自适应KF的时序InSAR流域库岸形变预测方法 | |
CN108876829B (zh) | 基于非线性尺度空间及径向基函数的sar高精度配准方法 | |
CN114114257A (zh) | 一种坝区形变与水位相关性检测方法和装置 | |
Shah | Automated lake shoreline mapping at subpixel accuracy | |
CN113505764A (zh) | 顾及局部入射角的高山峡谷区sar几何畸变识别方法 | |
CN107369163B (zh) | 一种基于最佳熵双阈值分割的快速sar图像目标检测方法 | |
CN115457022B (zh) | 基于实景三维模型正视影像的三维形变检测方法 | |
CN116908854A (zh) | 海浪谱分析法和基于Canny算子的雷达图像几何阴影统计法相结合的海浪参数反演方法 | |
CN103065297A (zh) | 一种基于傅里叶变换的图像边缘检测方法 |
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 |