CN104765026A - 合成孔径雷达干涉测量数据中地面属性数据的提取方法 - Google Patents
合成孔径雷达干涉测量数据中地面属性数据的提取方法 Download PDFInfo
- Publication number
- CN104765026A CN104765026A CN201510214592.1A CN201510214592A CN104765026A CN 104765026 A CN104765026 A CN 104765026A CN 201510214592 A CN201510214592 A CN 201510214592A CN 104765026 A CN104765026 A CN 104765026A
- Authority
- CN
- China
- Prior art keywords
- data
- ground
- synthetic aperture
- aperture radar
- point
- 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
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
- 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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2927—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value
-
- 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
-
- 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/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
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)
- Radar Systems Or Details Thereof (AREA)
- Image Processing (AREA)
Abstract
本发明涉及合成孔径雷达干涉测量数据中地面属性数据的提取方法。目的在于提出了PS数据中地面属性数据的提取方法。步骤如下:导入地面地理数据;将地面边界的数据沿东西方向内收Tx,将地面边界的数据沿南北方向内收Ty,重新确定地面新边界;判断每个PS点是否落入地面新边界,如果是,提取该数据;如果否,剔除该数据;获得提取后的第一数据集;判断第一数据集内PS点图像的灰度值是否在Vmin与Vmax之间,如果是,进一步提取该数据;如果否,剔除该数据。经过本发明方法对地面属性数据提取,PS点的正确率为95.1%。
Description
技术领域:
本发明涉及合成孔径雷达技术,进一步涉及合成孔径雷达干涉测量数据中地面属性数据的提取方法。
背景技术:
作为一种重要的遥感对地监测手段,合成孔径雷达干涉测量(Interferometry Synthetic Aperture Radar,InSAR)可以监测的目标包含有地表各种复杂的地物,如建筑物、地面、桥梁、地面等。不同的地物覆盖基础结构差异极大,地面沉降表现出来的特征受多种因素的影响。对大区域尺度的地面形变的分析,可近似认为像元点的形变率代表了地面的真实形变,然而,对建筑物和工程设施密集的城市,这种假设的合理性就不存在。以建筑物为例,建筑物的桩基结构千差万别,特别是高层建筑物自身地基处理时有桩基加固措施,桩基的着力点在地层深处。建筑自身的重力所表现出的下沉现象,反映的是建筑物通过桩基与深层地层的相互作用的结果,建筑物与地表的支撑力没有直接关系,与无附载的地表所表现的地面下沉有完全的区别,无附载地表的沉降与建筑物的沉降无论在时间与速率上都有不同的表现。因此从沉降目标分析,有必要将建筑物、地面等不同类型的永久散射体(Permanent Scatter,PS)点区分开来。
从沉降机理角度分析,地面沉降类型按成因可划分为自然因素和人为因素类型。自然因素类型主要为地质构造运动和土次固结引发的沉降等;人为因素类型主要有超量抽取地下水、工程建设引发的地面沉降等。且两者表现出来的特征也不相同。前者具有整体性,其对应的监测信息的特点在于一定范围内具有变化缓慢、梯度小等特点,具有可内插性;后者具有局部性,其对应的监测信息的特点在于相对于周围地物而言具有变化梯度大、影响范围有限等特点。因此通过从沉降机理的理论角度进行分析可知,不同类型地物目标的沉降机理存在差异,沉降信息特征也因此表现出不同的特征,因此应当开展PS点的属性分类工作。
发明内容:
本发明的目的在于提出基于地理信息数据库,提出了PS数据中地面属性数据的提取方法。
合成孔径雷达干涉测量数据中地面属性数据的提取方法,包括如下步骤:
导入地面地理数据,至少包括地面边界的经度和纬度;
将地面边界的数据沿东西方向内收Tx,将地面边界的数据沿南北方向内收Ty,重新确定地面新边界;所述Tx为合成孔径雷达干涉测量数据东西方向的定位误差,Ty为合成孔径雷达干涉测量数据东西方向的定位误差;
判断每个PS点是否落入地面新边界,如果是,提取该数据;如果否,剔除该数据;获得提取后的第一数据集;
判断第一数据集内PS点图像的灰度值是否在Vmin与Vmax之间,如果是,进一步提取该数据;如果否,剔除该数据;Vmin为地面属性PS点的最小亮度值,Vmax为地面PS点的最大亮度值;
获得进一步提取后的第二数据集,该数据集内的PS点属性为地面。
作为优选方案,Vmin=60;Vmax=122。
本发明相对于现有技术的优点在于:
(一)经过时序InSAR数据采用本发明方法对地面属性数据提取,PS点的正确率为95.6%。
(二)属性数据提取后进行的地面沉降监测结果与水准测量相比的符合精度优于5mm,能满足PS-InSAR开展高精度地面沉降监测的需求。
附图说明:
图1代表本发明在东西方向、南北方向内收的示意图;图中,外部的矩形代表地理数据所确定的地面边界,内部矩形代表内收后所确定的地面边界。
具体实施方式:
实施例:
(1)地理编码误差分析
在实验区范围内共选择27个地面点进行其地理编码误差分析,误差统计如表1所示。X向(东西向)中误差为1.19米,Y向(南北向)中误差为0.94米,平面总中误差为1.52米。
(2)二维空间分析
以实验区内地面为例,说明数据提取的过程:
覆盖面积1.3平方公里,经过时序InSAR数据处理后共有PS点2931个。
首先,导入地面地理数据,至少包括地面边界的经度和纬度;
将地面边界的数据沿东西方向内收1.19米,将地面边界的数据沿南北方向内收0.94米,重新确定地面新边界;
判断2931个PS点中每个PS点是否落入地面新边界,如果是,提取该数据;如果否,剔除该数据;
获得提取后的第一数据集,该数据集内包含1265个PS点;经过分析异常点为206个,正确率为83.7%。
(3)影像特征分析
选择实验区内选择14个区域进行地面的亮度信息分析,亮度统计信息如表2所示。采用最值分析法可知所有地面的最小值为53,最大值为159。
判断第一数据集内PS点图像的灰度值是否在53与159之间,如果是,进一步提取该数据;如果否,剔除该数据;
获得进一步提取后的第二数据集,该数据集内的PS点属性为地面。最终获得属于该地面的高可信度的PS点1075个,经过分析,其中异常点为53个,正确率为95.1%。
最终获得属于该地面的高可信度的PS点1075个,监测密度为826个/平方公里。最大沉降点为-14.76mm/yr,平均沉降量为-5.77mm/yr。
表1 地面点地理编码误差分析
No | X90 | Y90 | GCP_X | GCP_Y | dx | dy |
1 | 141825.66 | 285283.88 | 141824.56 | 285284.8638 | 1.10 | -0.98 |
2 | 141680.20 | 284635.45 | 141680.17 | 284634.2978 | 0.03 | 1.16 |
3 | 142326.89 | 284810.62 | 142326.53 | 284811.3343 | 0.36 | -0.71 |
4 | 141306.14 | 285282.01 | 141305.26 | 285281.0458 | 0.89 | 0.96 |
5 | 142948.60 | 285089.38 | 142948.05 | 285088.6009 | 0.55 | 0.78 |
6 | 143579.62 | 285769.78 | 143578.99 | 285770.4259 | 0.63 | -0.65 |
7 | 143883.68 | 285426.99 | 143882.79 | 285427.1965 | 0.89 | -0.20 |
8 | 142260.31 | 284121.51 | 142261.19 | 284120.9835 | -0.88 | 0.52 |
9 | 141959.80 | 284112.65 | 141961.39 | 284112.2042 | -1.59 | 0.45 |
10 | 141759.70 | 282279.89 | 141760.61 | 282278.213 | -0.90 | 1.68 |
11 | 141918.86 | 281340.69 | 141917.59 | 281341.6432 | 1.28 | -0.95 |
12 | 143453.57 | 282094.66 | 143453.01 | 282095.2933 | 0.56 | -0.63 |
13 | 144492.34 | 284237.73 | 144493.57 | 284238.26 | -1.22 | -0.54 |
14 | 144915.82 | 284589.45 | 144914.09 | 284590.05 | 1.73 | -0.61 |
15 | 145030.94 | 284421.13 | 145029.22 | 284420.89 | 1.72 | 0.24 |
16 | 144967.24 | 283936.09 | 144967.86 | 283935.37 | -0.62 | 0.72 |
17 | 145867.18 | 283228.88 | 145864.56 | 283226.75 | 2.62 | 2.13 |
18 | 146540.98 | 283741.61 | 146540.98 | 283742.09 | 0.00 | -0.47 |
19 | 146424.90 | 283028.60 | 146426.53 | 283028.32 | -1.63 | 0.28 |
20 | 148425.46 | 282882.72 | 148424.18 | 282882.97 | 1.29 | -0.25 |
21 | 148475.11 | 283086.69 | 148474.82 | 283087.91 | 0.29 | -1.22 |
22 | 144435.61 | 286397.67 | 144436.90 | 286396.80 | -1.29 | 0.86 |
23 | 148725.11 | 285271.97 | 148726.11 | 285273.76 | -1.00 | -1.79 |
24 | 147020.34 | 285098.09 | 147019.18 | 285099.78 | 1.16 | -1.69 |
25 | 147187.60 | 285620.43 | 147188.07 | 285620.09 | -0.48 | 0.33 |
26 | 146288.99 | 286535.71 | 146287.28 | 286535.61 | 1.70 | 0.10 |
27 | 146959.38 | 286427.22 | 146957.96 | 286426.77 | 1.42 | 0.45 |
X向中误差 | 1.19 | Y向中误差 | 0.94 | 总中误差 | 1.52 |
表中,X90代表,Y90代表,GCP_X代表,GCP_Y代表,dx代表,dy代表。
表2 地面点灰度值分析表
No | 最小值 | 最大值 | 均值 |
1 | 65 | 144 | 103 |
2 | 80 | 100 | 90 |
3 | 57 | 159 | 104 |
4 | 53 | 159 | 105 |
5 | 70 | 159 | 124 |
6 | 104 | 159 | 136 |
7 | 86 | 159 | 126 |
8 | 90 | 157 | 118 |
9 | 68 | 159 | 115 |
10 | 56 | 158 | 101 |
11 | 75 | 150 | 115 |
12 | 62 | 159 | 108 |
13 | 77 | 156 | 116 |
14 | 75 | 159 | 121 |
总计 | 53 | 159 | 113 |
Claims (2)
1.合成孔径雷达干涉测量数据中地面属性数据的提取方法,其特征在于,包括如下步骤:
导入地面地理数据,至少包括地面边界的经度和纬度;
将地面边界的数据沿东西方向内收Tx,将地面边界的数据沿南北方向内收Ty,重新确定地面新边界;所述Tx为合成孔径雷达干涉测量数据东西方向的定位误差,Ty为合成孔径雷达干涉测量数据东西方向的定位误差;
判断每个PS点是否落入地面新边界,如果是,提取该数据;如果否,剔除该数据;获得提取后的第一数据集;
判断第一数据集内PS点图像的灰度值是否在Vmin与Vmax之间,如果是,进一步提取该数据;如果否,剔除该数据;Vmin为地面属性PS点的最小亮度值,Vmax为地面PS点的最大亮度值;
获得进一步提取后的第二数据集,该数据集内的PS点属性为地面。
2.根据权利要求1所述合成孔径雷达干涉测量数据中地面属性数据的提取方法,其特征在于,Vmin=53;Vmax=159。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510214592.1A CN104765026A (zh) | 2015-04-29 | 2015-04-29 | 合成孔径雷达干涉测量数据中地面属性数据的提取方法 |
US15/051,674 US20160320479A1 (en) | 2015-04-29 | 2016-02-24 | Method for extracting ground attribute permanent scatter in interferometry synthetic aperture radar data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510214592.1A CN104765026A (zh) | 2015-04-29 | 2015-04-29 | 合成孔径雷达干涉测量数据中地面属性数据的提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104765026A true CN104765026A (zh) | 2015-07-08 |
Family
ID=53646976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510214592.1A Pending CN104765026A (zh) | 2015-04-29 | 2015-04-29 | 合成孔径雷达干涉测量数据中地面属性数据的提取方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20160320479A1 (zh) |
CN (1) | CN104765026A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108537885A (zh) * | 2018-04-19 | 2018-09-14 | 天津市测绘院 | 山体创面三维地形数据的获取方法 |
CN111308469A (zh) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | 一种基于PSInSAR技术的建筑物高程测量方法 |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT201600127152A1 (it) * | 2016-12-15 | 2018-06-15 | Ids Georadar S R L | Metodo e apparecchiatura per il monitoraggio di deformazioni superficiali di uno scenario |
WO2019126972A1 (zh) * | 2017-12-26 | 2019-07-04 | 深圳市城市公共安全技术研究院有限公司 | 基于InSAR的形变信息提取方法、终端及存储介质 |
US11835619B2 (en) | 2018-05-16 | 2023-12-05 | Nec Corporation | Synthetic aperture radar signal analysis device, synthetic aperture radar signal analysis method, and synthetic aperture radar signal analysis program |
CN109039422B (zh) * | 2018-06-28 | 2021-04-13 | 上海卫星工程研究所 | 深空探测高增益天线在轨定标***及方法 |
CN109840249B (zh) * | 2019-01-17 | 2021-05-04 | 武汉大学 | 一种流域库岸形变监控数据处理与集成方法 |
CN109752715B (zh) * | 2019-01-24 | 2023-05-12 | 深圳市数字城市工程研究中心 | 一种sar数据全散射体探测方法及装置 |
CN110888132B (zh) * | 2019-11-22 | 2022-05-20 | 深圳市城市公共安全技术研究院有限公司 | 基于InSAR监测的桥梁形变分析方法及*** |
CN112526515A (zh) * | 2020-11-05 | 2021-03-19 | 山西省交通科技研发有限公司 | 一种基于合成孔径雷达干涉测量技术的地表形变检测方法 |
CN112698328B (zh) * | 2020-11-30 | 2021-08-10 | 四川大学 | 一种用于大坝及滑坡变形gb-sar监测的相位解缠方法及*** |
CN112986994B (zh) * | 2021-02-06 | 2023-11-17 | 中国人民解放军国防科技大学 | 一种sar层析参考网快速生成方法 |
CN113446989B (zh) * | 2021-06-10 | 2023-04-18 | 中铁隧道局集团有限公司 | 基于dic与合成孔径雷达的周围土体变形变位的空间检测方法 |
CN114265062B (zh) * | 2021-11-11 | 2023-11-10 | 电子科技大学 | 一种基于相位梯度估计网络的InSAR相位解缠方法 |
CN115201822B (zh) * | 2022-07-07 | 2023-03-14 | 长沙理工大学 | 一种钻井水溶岩盐矿区采卤量估计方法 |
CN115982968B (zh) * | 2022-12-14 | 2023-07-14 | 首都师范大学 | 一种揭示地面时空间演化差异的算法 |
CN116908853B (zh) * | 2023-09-13 | 2023-11-17 | 北京观微科技有限公司 | 高相干点选取方法、装置和设备 |
CN117572420B (zh) * | 2023-11-14 | 2024-04-26 | 中国矿业大学 | 一种基于深度学习的InSAR相位解缠优化方法 |
CN117310706B (zh) * | 2023-11-28 | 2024-02-02 | 中山大学 | 一种地基雷达间断形变监测方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6583751B1 (en) * | 1999-05-25 | 2003-06-24 | Politecnico Di Milano | Process for radar measurements of the movement of city areas and landsliding zones |
CN103970932A (zh) * | 2014-02-28 | 2014-08-06 | 杭州师范大学 | 一种高分辨率的建筑物和背景分离的永久散射体建模方法 |
-
2015
- 2015-04-29 CN CN201510214592.1A patent/CN104765026A/zh active Pending
-
2016
- 2016-02-24 US US15/051,674 patent/US20160320479A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6583751B1 (en) * | 1999-05-25 | 2003-06-24 | Politecnico Di Milano | Process for radar measurements of the movement of city areas and landsliding zones |
CN103970932A (zh) * | 2014-02-28 | 2014-08-06 | 杭州师范大学 | 一种高分辨率的建筑物和背景分离的永久散射体建模方法 |
Non-Patent Citations (5)
Title |
---|
丁尚起 等: "PS选择策略及在京津高铁沉降监测中的应用研究", 《城市勘测》 * |
兰恒星 等: "城市复杂地面沉降永久干涉雷达监测属性分类研究", 《工程地质学报》 * |
兰恒星 等: "城市建筑物沉降永久散射体干涉雷达监测", 《天津大学学报》 * |
周立凡: "城市重大工程区高分辨率永久散射体雷达干涉地表形变监测", 《中国博士学位论文全文数据库 基础科学辑》 * |
陈基伟: "PS-InSAR技术地面沉降研究与展望", 《测绘科学》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108537885A (zh) * | 2018-04-19 | 2018-09-14 | 天津市测绘院 | 山体创面三维地形数据的获取方法 |
CN111308469A (zh) * | 2019-11-27 | 2020-06-19 | 北京东方至远科技股份有限公司 | 一种基于PSInSAR技术的建筑物高程测量方法 |
CN111308469B (zh) * | 2019-11-27 | 2021-07-16 | 北京东方至远科技股份有限公司 | 一种基于PSInSAR技术的建筑物高程测量方法 |
Also Published As
Publication number | Publication date |
---|---|
US20160320479A1 (en) | 2016-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104765026A (zh) | 合成孔径雷达干涉测量数据中地面属性数据的提取方法 | |
CN1918451B (zh) | 用于从单一图像对3d对象进行测量、建模和地图制作的***和方法 | |
CN101661497B (zh) | 遥感土地利用变化检测方法及*** | |
CN102426395B (zh) | 一种基于三维地质遥感判释技术的滑坡体体积参数量测方法 | |
CN103049655B (zh) | 基于多数据源的矿区土地生态损伤信息获取方法 | |
CN103047970B (zh) | 一种确定塌陷耕地损害边界的方法 | |
CN101714211B (zh) | 高分辨率遥感影像道路中心线的检测的方法 | |
Triantafilis et al. | Digital soil-class mapping across the Edgeroi district using numerical clustering and gamma-ray spectrometry data | |
Khalid et al. | Application of Kriging for development of SPT N value contour maps and USCS-based soil type qualitative contour maps for Islamabad, Pakistan | |
Yan et al. | Detection and numerical simulation of potential hazard in oil pipeline areas based on UAV surveys | |
CN109593958A (zh) | 三维遥感技术辅助离子型稀土矿原地浸取工艺的方法 | |
CN110111376B (zh) | 一种采煤沉陷水域面积计算方法 | |
CN104765887A (zh) | 合成孔径雷达干涉测量数据中道路属性数据的提取方法 | |
An et al. | Ground subsidence monitoring in based on UAV-LiDAR technology: a case study of a mine in the Ordos, China | |
CN114386137B (zh) | 基于倾斜摄影技术的道路地形曲面优化设计方法 | |
CN104793207A (zh) | 合成孔径雷达干涉测量数据中建筑物属性数据的提取方法 | |
Kadhim et al. | Digital geotechnical maps of Basrah city using geographical information systems technique | |
Rajabi | Sinkhole Detection and Quantification Using LiDAR Data | |
CN104597518A (zh) | 确定海域油气勘探有利区的方法 | |
CN113653085A (zh) | 工程桩桩长预估方法及其预估*** | |
CN110412631A (zh) | 一种线性工程gps快速放样施工方法 | |
CN110008302A (zh) | 多影响因子降雨入渗系数定量评价方法 | |
Gumilar et al. | 3D modelling of Mt. Talaga Bodas Crater (Indonesia) by using terrestrial laser scanner for volcano hazard mitigation | |
Xu et al. | Simulation study of rockfall deposition based on UAV-PCAS-PFC: a case study on the rockfall deposition of Xiaomaopo in Kaiyang County | |
Cao et al. | Analysis on rock probability of embedment pile into rock mass for pile group at building sites in Karst Terrain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150708 |
|
RJ01 | Rejection of invention patent application after publication |