CN116908789B - 地基合成孔径雷达干涉测量建筑立面变形信息提取方法 - Google Patents

地基合成孔径雷达干涉测量建筑立面变形信息提取方法 Download PDF

Info

Publication number
CN116908789B
CN116908789B CN202311179142.4A CN202311179142A CN116908789B CN 116908789 B CN116908789 B CN 116908789B CN 202311179142 A CN202311179142 A CN 202311179142A CN 116908789 B CN116908789 B CN 116908789B
Authority
CN
China
Prior art keywords
building
deformation
image pixels
radar
sar image
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
Application number
CN202311179142.4A
Other languages
English (en)
Other versions
CN116908789A (zh
Inventor
李双平
陈文夫
王鹏
张斌
段杭
柯传芳
陈一鸣
郑敏
史波
于胜利
徐福斌
刘德军
许强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changjiang Spatial Information Technology Engineering Co ltd
China Three Gorges Construction Engineering Co Ltd
Original Assignee
Changjiang Spatial Information Technology Engineering Co ltd
China Three Gorges Construction Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changjiang Spatial Information Technology Engineering Co ltd, China Three Gorges Construction Engineering Co Ltd filed Critical Changjiang Spatial Information Technology Engineering Co ltd
Priority to CN202311179142.4A priority Critical patent/CN116908789B/zh
Publication of CN116908789A publication Critical patent/CN116908789A/zh
Application granted granted Critical
Publication of CN116908789B publication Critical patent/CN116908789B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
    • 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/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
    • 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
    • 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
    • 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
    • G01S7/414—Discriminating targets with respect to background clutter
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models

Abstract

一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法，首先建立地理编码查找表，该地理编码查找表为GB‑SAR成像投影面上GB‑SAR影像像元与三维点云的匹配关系，然后获得建筑立面区域内的多个GB‑SAR影像像元，测量各GB‑SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB‑SAR影像像元相匹配的各三维点云，再根据各GB‑SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB‑SAR影像像元的立面变形值，最后基于各GB‑SAR影像像元的立面变形值得到建筑立面的变形分布。本设计能够实现全天时、全天候的建筑立面变形实时监测。

Description

地基合成孔径雷达干涉测量建筑立面变形信息提取方法
技术领域
本发明属于工程建筑变形监测技术领域，具体涉及一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法。
背景技术
建筑物受到外部环境（周边施工、恶劣气象环境等）的影响会发生变形，这种变形超过一定限度会直接影响建筑结构安全，如建筑周边的基坑开挖施工，邻近建筑物往往会朝向基坑开挖区域进行倾斜变化，由于在不同的高度、不同的位置其变形大小是不一样的，传统工程建筑变形监测技术的测点一般极为稀疏，很难同时获取建筑物整个立面上的高精度变形信息，且通常需要去往变形体范围内布设棱镜等协作目标进行观测，这在受地质灾害影响严重的高危建筑监测应用中是非常危险的，地基合成孔径雷达干涉测量（Ground-based Synthetic Aperture Radar Interferometry，GB-InSAR）技术的出现为非接触遥感变形监测提供了新的途径。因此，亟需一种将GB-SAR技术应用到建筑立面变形分布信息监测中的方法。
发明内容
本发明的目的是针对现有技术存在的上述问题，提供一种能够实现建筑立面变形分布监测的地基合成孔径雷达干涉测量建筑立面变形信息提取方法。
为实现以上目的，本发明的技术方案如下：
地基合成孔径雷达干涉测量建筑立面变形信息提取方法，所述提取方法包括：
S1、建立地理编码查找表，所述地理编码查找表为GB-SAR成像投影面上GB-SAR影像像元与三维点云的匹配关系；
S2、获得建筑立面区域内的多个GB-SAR影像像元，然后测量各GB-SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB-SAR影像像元相匹配的各三维点云；
S3、根据各GB-SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB-SAR影像像元的立面变形值；
S4、基于各GB-SAR影像像元的立面变形值得到建筑立面的变形分布。
所述步骤S1包括：
先计算三维点云在GB-SAR成像平面上的斜距投影坐标、GB-SAR影像像元在GB-SAR成像平面上的直角坐标，然后在GB-SAR成像投影面上搜索与各GB-SAR影像像元距离最小的三维点云，完成各GB-SAR影像像元与各三维点云的一一匹配。
所述GB-SAR影像像元在GB-SAR成像平面上的直角坐标根据以下公式计算得到：
；
；
上式中，为影像像元在GB-SAR成像平面上的直角坐标，/>为影像像元相对于雷达中心线的偏角，/>为方位向分辨率，/>为方位向分辨单元序号，/>为首列像元偏角。
所述三维点云在GB-SAR成像平面上的斜距投影坐标根据以下公式计算得到：
；
；
；
；
上式中，为三维点云（x，y，z）在GB-SAR成像平面上的斜距投影坐标，即三维点云（x，y，z）依据GB-SAR成像几何进行斜距投影至雷达中心线所在平面获得的坐标，/>为三维点云（x，y，z）相对于雷达中心的斜距，/>为斜距投影后的坐标缩放系数，为三维点云（x，y，z）相对于GB-SAR成像平面的三维坐标，所述三维点云（x，y，z）是原始三维点云经过原点平移和y轴旋转后形成的雷达坐标系下的三维坐标，该雷达坐标系下以雷达中心为坐标原点，以雷达中心线为y轴并指向成像目标，以雷达运行导轨的中心线为x轴；ω、κ分别为雷达姿态的俯仰角、侧倾角，所述俯仰角为雷达中心视线与水平面之间的夹角，所述侧倾角为雷达运行轨道中轴线与水平面之间的夹角。
步骤S2中，所述建筑立面区域内的多个GB-SAR影像像元根据以下方法获得：
A1、先计算获得的建筑GB-SAR影像像元的热信噪比，预选出热信噪比高于对应设定阈值的建筑GB-SAR影像像元，然后分别计算预选出的建筑GB-SAR影像像元的振幅离散指数、相关系数，保留振幅离散指数、相关系数均高于对应设定阈值的建筑GB-SAR影像像元；
A2、对得到的多个GB-SAR影像像元进行建筑立面区域剪裁，以剔除噪点和建筑立面区域外的测点，获得建筑立面区域内的多个GB-SAR影像像元。
步骤A1中，所述热信噪比根据以下公式计算得到：
；
上式中，为热信噪比，/>为接收机热噪声功率，/>为影像像元复数，/>为取平均值的计算；
所述振幅离散指数根据以下公式计算得到：
；
上式中，为振幅离散指数，/>、/>分别为影像像元振幅值A的标准差、均值；
所述相关系数根据以下公式计算得到：
；
上式中，为相关系数，i、j分别为影像像元的行、列号，m、n分别为方形窗口的长和宽，对m行、n列的方形窗口内所有影像像元均进行相关系数计算；*表示取共轭复数，M、S分别代表干涉对的主影像和从影像。
步骤S3中，所述GB-SAR影像像元的立面变形值根据以下公式计算得到：
；
；
；
上式中，d为影像像元的立面变形值，为由雷达技术测定的影像像元在雷达视线向的变形值，/>为影像像元至雷达中心的斜距，z为与该影像像元相匹配的三维点云（x，y，z）中的高程，所述三维点云（x，y，z）是原始三维点云（X，Y，Z）经过原点平移至P点和y轴旋转后形成的雷达坐标系下的三维坐标，/>为建筑立面的倾斜角度。
步骤S4中，获得各GB-SAR影像像元的立面变形值后，利用Kriging插值方法得到建筑立面的变形分布。
所述步骤S4还包括：
基于各GB-SAR影像像元的立面变形值得到建筑立面上任一纵线的相对倾斜度。
步骤S4中，指定建筑立面某一位置的纵线，该位置的相对倾斜度根据以下公式计算得到：
；
上式中，为相对倾斜度，/>为变形点g相对于雷达中心的高差，/>为变形点g处的立面变形值，G为当前纵线上所选取的变形点总数。
与现有技术相比，本发明的有益效果为：
1、本发明一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法，首先建立地理编码查找表，该地理编码查找表为GB-SAR成像投影面上GB-SAR影像像元与三维点云的匹配关系，然后获得建筑立面区域内的多个GB-SAR影像像元，测量各GB-SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB-SAR影像像元相匹配的各三维点云，再根据各GB-SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB-SAR影像像元的立面变形值，最后基于各GB-SAR影像像元的立面变形值得到建筑立面的变形分布，最终实现全天时、全天候的建筑立面变形实时监测。因此，本发明能够实现全天时、全天候的建筑立面变形实时监测。
2、本发明一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法中，先计算三维点云在GB-SAR成像平面上的斜距投影坐标、GB-SAR影像像元在GB-SAR成像平面上的直角坐标，然后在GB-SAR成像投影面上搜索与各GB-SAR影像像元距离最小的三维点云，完成各GB-SAR影像像元与各三维点云的一一匹配，构建地理编码查找表，在计算三维点云在GB-SAR成像平面上的斜距投影坐标时充分考虑了雷达姿态倾角的影响，即雷达的俯仰角与雷达运行导轨的侧倾角，从而准确地实现GB-SAR影像像元的坐标变换，提高最终建筑变形监测结果的精度和可靠性。因此，本发明在应用三维点云辅助GB-SAR影像像元坐标变换时充分考虑了雷达姿态倾角的影响，提高了最终建筑变形监测结果的精度和可靠性。
3、本发明一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法中，先计算获得的建筑GB-SAR影像像元的热信噪比，预选出热信噪比高于对应设定阈值的建筑GB-SAR影像像元，滤除大部分虚假信号和弱信号的影像像元后，分别计算预选出的建筑GB-SAR影像像元的振幅离散指数、相关系数，保留振幅离散指数、相关系数均高于对应设定阈值的建筑GB-SAR影像像元，再对保留多个GB-SAR影像像元进行建筑立面区域剪裁，以剔除噪点和建筑立面区域外的测点，获得建筑立面区域内的多个GB-SAR影像像元；经过以上筛选后，不仅最终获得建筑立面区域内GB-SAR影像像元的质量高，而且能够减少后续计算量，显著提高计算效率。因此，本发明不仅能够获得高质量的GB-SAR影像像元，而且能够显著提高计算效率。
附图说明
图1为本发明的流程图。
图2为实施例中高质量GB-SAR影像像元的获取流程图。
具体实施方式
下面结合具体实施方式以及附图对本发明作进一步详细的说明。
参见图1，地基合成孔径雷达干涉测量建筑立面变形信息提取方法，所述提取方法包括：
S1、建立地理编码查找表，所述地理编码查找表为GB-SAR成像投影面上GB-SAR影像像元与三维点云的匹配关系；
S2、获得建筑立面区域内的多个GB-SAR影像像元，然后测量各GB-SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB-SAR影像像元相匹配的各三维点云；
S3、根据各GB-SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB-SAR影像像元的立面变形值；
S4、基于各GB-SAR影像像元的立面变形值得到建筑立面的变形分布。
所述步骤S1包括：
先计算三维点云在GB-SAR成像平面上的斜距投影坐标、GB-SAR影像像元在GB-SAR成像平面上的直角坐标，然后在GB-SAR成像投影面上搜索与各GB-SAR影像像元距离最小的三维点云，完成各GB-SAR影像像元与各三维点云的一一匹配。
所述GB-SAR影像像元在GB-SAR成像平面上的直角坐标根据以下公式计算得到：
；
；
上式中，为影像像元在GB-SAR成像平面上的直角坐标，/>为影像像元相对于雷达中心线的偏角，/>为方位向分辨率，/>为方位向分辨单元序号，/>为首列像元偏角。
所述三维点云在GB-SAR成像平面上的斜距投影坐标根据以下公式计算得到：
；
；
；
；
上式中，为三维点云（x，y，z）在GB-SAR成像平面上的斜距投影坐标，即三维点云（x，y，z）依据GB-SAR成像几何进行斜距投影至雷达中心线所在平面获得的坐标，/>为三维点云（x，y，z）相对于雷达中心的斜距，/>为斜距投影后的坐标缩放系数，为三维点云（x，y，z）相对于GB-SAR成像平面的三维坐标，所述三维点云（x，y，z）是原始三维点云经过原点平移和y轴旋转后形成的雷达坐标系下的三维坐标，该雷达坐标系下以雷达中心为坐标原点，以雷达中心线为y轴并指向成像目标，以雷达运行导轨的中心线为x轴；ω、κ分别为雷达姿态的俯仰角、侧倾角,所述俯仰角为雷达中心视线与水平面之间的夹角，所述侧倾角为雷达运行轨道中轴线与水平面之间的夹角。
步骤S2中，所述建筑立面区域内的多个GB-SAR影像像元根据以下方法获得：
A1、先计算获得的建筑GB-SAR影像像元的热信噪比，预选出热信噪比高于对应设定阈值的建筑GB-SAR影像像元，然后分别计算预选出的建筑GB-SAR影像像元的振幅离散指数、相关系数，保留振幅离散指数、相关系数均高于对应设定阈值的建筑GB-SAR影像像元，将其作为GB-SAR影像像元；
A2、对得到的多个GB-SAR影像像元进行建筑立面区域剪裁，以剔除噪点和建筑立面区域外的测点，获得建筑立面区域内的多个GB-SAR影像像元。
步骤A1中，所述热信噪比根据以下公式计算得到：
；
上式中，为热信噪比，/>为接收机热噪声功率，/>为影像像元复数，/>为取平均值的计算；
所述振幅离散指数根据以下公式计算得到：
；
上式中，为振幅离散指数，/>、/>分别为影像像元振幅值A的标准差、均值；
所述相关系数根据以下公式计算得到：
；
上式中，为相关系数，i、j分别为影像像元的行、列号，m、n分别为方形窗口的长和宽，对m行、n列的方形窗口内所有影像像元均进行相关系数计算；*表示取共轭复数，M、S分别代表干涉对的主影像和从影像。
步骤S3中，所述GB-SAR影像像元的立面变形值根据以下公式计算得到：
；
；
上式中，d为影像像元的立面变形值，为由雷达技术测定的影像像元在雷达视线向的变形值，/>为影像像元至雷达中心的斜距，z为与该影像像元相匹配的三维点云（x，y，z）中的高程，所述三维点云（x，y，z）是原始三维点云（X，Y，Z）经过原点平移至P点和y轴旋转后形成的雷达坐标系下的三维坐标，/>为建筑立面的倾斜角度。
步骤S4中，获得各GB-SAR影像像元的立面变形值后，利用Kriging插值方法得到建筑立面的变形分布。
所述步骤S4还包括：
基于各GB-SAR影像像元的立面变形值得到建筑立面上任一纵线的相对倾斜度。
步骤S4中，指定建筑立面某一位置的纵线，该位置的相对倾斜度根据以下公式计算得到：
；
上式中，为相对倾斜度，/>为变形点g相对于雷达中心的高差，/>为变形点g处的立面变形值，G为当前纵线上所选取的变形点总数。
本发明的原理说明如下：
本发明一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法将GB-SAR技术与建筑立面变形分布监测相结合，由于地基合成孔径雷达测站点布设灵活，成像时空分辨率高，干涉测量精度高，十分适用于对工程结构和地质灾害的全天时、全天候的建筑立面变形实时监测。
GB-SAR技术特殊的成像方式给目标位置的准确识别以及变形的解译分析带来了较大困难，一般需借助测区三维数据（如三维激光扫描点云、数字高程模型或倾斜摄影测量模型等）恢复GB-SAR影像的三维空间关系，但传统的GB-SAR影像像元坐标与三维测量数据坐标的变换方法中往往忽略了雷达传感器姿态带来的影响，成像投影面一般均认定为水平面，而在实际GB-SAR的监测应用中，常常会遇到监测区域同设站点高差较大使得水平视角无法对准核心变形区域，也可能会由于设站条件限制、施工因素或人为错误导致雷达传感器姿态和运行导轨不水平，此外，地基真实孔径雷达的转动平面倾斜、线性GB-SAR导轨倾斜、ArcSAR连杆旋转平面和馈源天线面板的倾斜均会导致成像投影面的倾斜；上述诸多因素使得雷达传感器无法处于理想水平姿态，直接影响中短距离建筑结构变形区位的准确识别；而本发明在充分考虑雷达传感器姿态的影响下对点云斜距投影坐标进行计算，再在此基础上，根据最小距离准则，对GB-SAR影像像元和三维点云数据进行匹配计算，建立了更为可靠的地理编码查找表，从而准确地实现GB-SAR影像像元坐标与三维点云坐标的变换，进而准确提取建筑立面变形和建筑纵线相对倾斜等变形信息。
实施例1：
参见图1，一种地基合成孔径雷达干涉测量建筑立面变形信息提取方法，依次按照以下步骤进行：
S1、根据以下步骤建立地理编码查找表，所述地理编码查找表为在GB-SAR成像投影面上GB-SAR影像像元与三维点云的匹配关系：
首先，计算三维点云在GB-SAR成像平面上的斜距投影坐标、GB-SAR影像像元在GB-SAR成像平面上的直角坐标，其中，所述GB-SAR影像像元在GB-SAR成像平面上的直角坐标根据以下公式计算得到：
；
；
上式中，为影像像元在GB-SAR成像平面上的直角坐标，/>为影像像元相对于雷达中心线的偏角，/>为方位向分辨率，/>为方位向分辨单元序号，/>为首列像元偏角；
依据GB-SAR成像几何，将各三维点云坐标进行斜距投影至雷达波束中心线所在平面，当雷达传感器姿态同时存在俯仰角ω和侧倾角κ时，三维点云在GB-SAR成像平面上的斜距投影坐标为：
；
；
；
；
上式中，为三维点云（x，y，z）在GB-SAR成像平面上的斜距投影坐标，即三维点云（x，y，z）依据GB-SAR成像几何进行斜距投影至雷达中心线所在平面获得的坐标，/>为三维点云（x，y，z）相对于雷达中心的斜距，/>为斜距投影后的坐标缩放系数，为三维点云（x，y，z）相对于GB-SAR成像平面的三维坐标，所述三维点云（x，y，z）是原始三维点云（X，Y，Z）经过原点平移至P点/>和y轴旋转后形成的雷达坐标系下的三维坐标，该雷达坐标系下以雷达中心为坐标原点，以雷达中心线为y轴并指向成像目标，以雷达运行导轨的中心线为x轴；ω、κ分别为雷达姿态的俯仰角、侧倾角，所述俯仰角为雷达中心视线与水平面之间的夹角，所述侧倾角为雷达运行轨道中轴线与水平面之间的夹角，令俯仰角ω等于零即可进行仅存在侧倾角κ的计算，令侧倾角κ等于零即可进行仅存在俯仰角ω的计算，更常见的情况是雷达姿态同时存在俯仰角ω和侧倾角κ；
然后，在GB-SAR成像投影面上搜索与GB-SAR影像像元距离最小的三维点云，完成各GB-SAR影像像元与各三维点云的一一匹配，在选定待计算的GB-SAR影像像元后，假设该GB-SAR影像像元的极坐标为（r，θ1），给定一个搜索范围（r±Δr，θ1±Δθ1），计算该搜索范围内的所有三维点云的斜距投影坐标到GB-SAR影像像元的直角坐标的距离，距离最小的那个三维点云即为该GB-SAR影像像元的匹配点，最终得到地理编码查找表；
S2、获得建筑立面区域内的多个GB-SAR影像像元，然后通过干涉计算和时序分析计算得到各GB-SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB-SAR影像像元相匹配的各三维点云；
其中，所述建筑立面区域内的多个GB-SAR影像像元根据以下方法获得：
A1、先计算获得的建筑GB-SAR影像像元的热信噪比，预选出热信噪比高于对应设定阈值的建筑GB-SAR影像像元，然后分别计算预选出的建筑GB-SAR影像像元的振幅离散指数、相关系数，保留振幅离散指数、相关系数均高于对应设定阈值的建筑GB-SAR影像像元，此时获得的GB-SAR影像像元为高质量的GB-SAR影像像元；
所述热信噪比根据以下公式计算得到：
；
上式中，为热信噪比，/>为接收机热噪声功率，/>为影像像元复数，/>为取平均值的计算；
所述振幅离散指数根据以下公式计算得到：
；
上式中，为振幅离散指数，/>、/>分别为影像像元振幅值A的标准差、均值；
所述相关系数根据以下公式计算得到：
；
上式中，为相关系数，i、j分别为影像像元的行、列号，m、n分别为方形窗口的长和宽，在相关系数计算时会以当前选定的影像像元为中心，对m行、n列的方形窗口内所有影像像元均进行计算，一般m、n均取大于2的奇数；*表示取共轭复数，M、S分别代表干涉对的主影像和从影像；
A2、选定建筑立面范围，对得到的多个GB-SAR影像像元进行建筑立面区域剪裁，以剔除噪点和建筑立面区域外的测点，获得建筑立面区域内的多个GB-SAR影像像元；
S3、根据各GB-SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB-SAR影像像元的立面变形值，所述GB-SAR影像像元的立面变形值根据以下公式计算得到：
；
；
；
上式中，d为影像像元的立面变形值，为由雷达技术测定的影像像元在雷达视线向的变形值，/>为影像像元至雷达中心的斜距，z为与该影像像元相匹配的三维点云（x，y，z）中的高程，/>为建筑立面的倾斜角度，该建筑立面的倾斜角度是指建筑的初始外形相对于铅垂线的倾斜角度，一般设置为零；
S4、基于各GB-SAR影像像元的立面变形值得到建筑立面的变形分布以及建筑立面上任一纵线的相对倾斜度；
其中，所述建筑立面的变形分布为利用Kriging插值方法得到，建筑立面的变形分布指的是立面变形值在建筑的外立面上是如何分布的，可以通过建筑立面变形值大小的分布情况来分析建筑的安全状态；
所述建筑立面上任一纵线的相对倾斜度根据以下公式计算得到：
；
上式中，为指定建筑立面某一位置的纵线的相对倾斜度，该建筑纵线相对倾斜度是用于以GB-SAR首次测量值为基准，评价建筑相对于GB-SAR首次测量值的倾斜变化情况，/>为变形点g相对于雷达中心的高差，/>为变形点g处的立面变形值，G为当前纵线上所选取的变形点总数。/>

Claims (10)

1.地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
所述提取方法包括以下步骤：
S1、建立地理编码查找表，所述地理编码查找表为GB-SAR成像投影面上GB-SAR影像像元与三维点云的匹配关系；
S2、获得建筑立面区域内的多个GB-SAR影像像元，然后测量各GB-SAR影像像元的LOS变形值，以及在地理编码查找表中找到与各GB-SAR影像像元相匹配的各三维点云；
S3、根据各GB-SAR影像像元的LOS变形值以及与其相匹配的三维点云，计算各GB-SAR影像像元的立面变形值；
S4、基于各GB-SAR影像像元的立面变形值得到建筑立面的变形分布。
2.根据权利要求1所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
所述步骤S1包括：
先计算三维点云在GB-SAR成像平面上的斜距投影坐标、GB-SAR影像像元在GB-SAR成像平面上的直角坐标，然后在GB-SAR成像投影面上搜索与各GB-SAR影像像元距离最小的三维点云，完成各GB-SAR影像像元与各三维点云的一一匹配。
3.根据权利要求2所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
所述GB-SAR影像像元在GB-SAR成像平面上的直角坐标根据以下公式计算得到：
θ＝Δθ·n+θ0；
上式中，(xradar，yradar)为影像像元在GB-SAR成像平面上的直角坐标，θ为影像像元相对于雷达中心线的偏角，Δθ为方位向分辨率，n为方位向分辨单元序号，θ0为首列像元偏角。
4.根据权利要求2或3所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
所述三维点云在GB-SAR成像平面上的斜距投影坐标根据以下公式计算得到：
上式中，(xradar′，yradar′)为三维点云(x，y，z)在GB-SAR成像平面上的斜距投影坐标，即三维点云(x，y，z)依据GB-SAR成像几何进行斜距投影至雷达中心线所在平面获得的坐标，S为三维点云(x，y，z)相对于雷达中心的斜距，k′为斜距投影后的坐标缩放系数，(x′，y′，z′)为三维点云(x，y，z)相对于GB-SAR成像平面的三维坐标，所述三维点云(x，y，z)是原始三维点云经过原点平移和y轴旋转后形成的雷达坐标系下的三维坐标，该雷达坐标系下以雷达中心为坐标原点，以雷达中心线为y轴并指向成像目标，以雷达运行导轨的中心线为x轴；ω、κ分别为雷达姿态的俯仰角、侧倾角，所述俯仰角为雷达中心视线与水平面之间的夹角，所述侧倾角为雷达运行轨道中轴线与水平面之间的夹角。
5.根据权利要求1-3中任一项所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
步骤S2中，所述建筑立面区域内的多个GB-SAR影像像元根据以下方法获得：
A1、先计算获得的建筑GB-SAR影像像元的热信噪比，预选出热信噪比高于对应设定阈值的建筑GB-SAR影像像元，然后分别计算预选出的建筑GB-SAR影像像元的振幅离散指数、相关系数，保留振幅离散指数、相关系数均高于对应设定阈值的建筑GB-SAR影像像元；
A2、对得到的多个GB-SAR影像像元进行建筑立面区域剪裁，以剔除噪点和建筑立面区域外的测点，获得建筑立面区域内的多个GB-SAR影像像元。
6.根据权利要求5所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
步骤A1中，所述热信噪比根据以下公式计算得到：
上式中，TSNR为热信噪比，为接收机热噪声功率，Sm为影像像元复数，E为取平均值的计算；
所述振幅离散指数根据以下公式计算得到：
上式中，DA为振幅离散指数，σA、mA分别为影像像元振幅值A的标准差、均值；
所述相关系数根据以下公式计算得到：
上式中，γ为相关系数，i、j分别为影像像元的行、列号，m、n分别为方形窗口的长和宽，对m行、n列的方形窗口内所有影像像元均进行相关系数计算，*表示取共轭复数，M、S分别代表干涉对的主影像和从影像。
7.根据权利要求1-3中任一项所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
步骤S3中，所述GB-SAR影像像元的立面变形值根据以下公式计算得到：
z＝Z-Zc；
上式中，d为影像像元的立面变形值，dLOS为由雷达技术测定的影像像元在雷达视线向的变形值，S′为影像像元至雷达中心的斜距，z为与该影像像元相匹配的三维点云(x，y，z)中的高程，所述三维点云(x，y，z)是原始三维点云(X，Y，Z)经过原点平移至P点(Xc，Yc，Zc)和y轴旋转后形成的雷达坐标系下的三维坐标，ε为建筑立面的倾斜角度。
8.根据权利要求1-3中任一项所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
步骤S4中，获得各GB-SAR影像像元的立面变形值后，利用Kriging插值方法得到建筑立面的变形分布。
9.根据权利要求1-3中任一项所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
所述步骤S4还包括：
基于各GB-SAR影像像元的立面变形值得到建筑立面上任一纵线的相对倾斜度。
10.根据权利要求9所述的地基合成孔径雷达干涉测量建筑立面变形信息提取方法，其特征在于：
步骤S4中，指定建筑立面某一位置的纵线，该位置的相对倾斜度根据以下公式计算得到：
上式中，k为相对倾斜度，hg为变形点g相对于雷达中心的高差，dg为变形点g处的立面变形值，G为当前纵线上所选取的变形点总数。
CN202311179142.4A 2023-09-13 2023-09-13 地基合成孔径雷达干涉测量建筑立面变形信息提取方法 Active CN116908789B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311179142.4A CN116908789B (zh) 2023-09-13 2023-09-13 地基合成孔径雷达干涉测量建筑立面变形信息提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311179142.4A CN116908789B (zh) 2023-09-13 2023-09-13 地基合成孔径雷达干涉测量建筑立面变形信息提取方法

Publications (2)

Publication Number Publication Date
CN116908789A CN116908789A (zh) 2023-10-20
CN116908789B true CN116908789B (zh) 2023-12-05

Family

ID=88353465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311179142.4A Active CN116908789B (zh) 2023-09-13 2023-09-13 地基合成孔径雷达干涉测量建筑立面变形信息提取方法

Country Status (1)

Country Link
CN (1) CN116908789B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772377A (zh) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 一种基于InSAR的建筑物变形监测方法
CN110244302A (zh) * 2019-07-05 2019-09-17 苏州科技大学 地基合成孔径雷达影像像元坐标三维变换方法
WO2019242174A1 (zh) * 2018-06-21 2019-12-26 华南理工大学 基于激光雷达的建筑结构自动测量及3d模型生成方法
WO2020068295A1 (en) * 2018-09-28 2020-04-02 Futurewei Technologies, Inc. System and method for determining line of sight (los)
WO2021227423A1 (zh) * 2020-05-13 2021-11-18 深圳大学 基于动态基线的InSAR数字高程模型构建方法及系统
WO2022213673A1 (zh) * 2021-04-06 2022-10-13 中国矿业大学 融合无人机dom和星载sar影像的地表三维形变提取方法
CN115951354A (zh) * 2023-02-14 2023-04-11 中国铁道科学研究院集团有限公司铁道建筑研究所 一种融合升降轨的D-InSAR形变监测方法
CN116168176A (zh) * 2023-04-23 2023-05-26 深圳大学 联合InSAR与激光点云的建筑物几何与变形提取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3457357B1 (en) * 2017-09-13 2021-07-07 Tata Consultancy Services Limited Methods and systems for surface fitting based change detection in 3d point-cloud

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772377A (zh) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 一种基于InSAR的建筑物变形监测方法
WO2019242174A1 (zh) * 2018-06-21 2019-12-26 华南理工大学 基于激光雷达的建筑结构自动测量及3d模型生成方法
WO2020068295A1 (en) * 2018-09-28 2020-04-02 Futurewei Technologies, Inc. System and method for determining line of sight (los)
CN110244302A (zh) * 2019-07-05 2019-09-17 苏州科技大学 地基合成孔径雷达影像像元坐标三维变换方法
WO2021227423A1 (zh) * 2020-05-13 2021-11-18 深圳大学 基于动态基线的InSAR数字高程模型构建方法及系统
WO2022213673A1 (zh) * 2021-04-06 2022-10-13 中国矿业大学 融合无人机dom和星载sar影像的地表三维形变提取方法
CN115951354A (zh) * 2023-02-14 2023-04-11 中国铁道科学研究院集团有限公司铁道建筑研究所 一种融合升降轨的D-InSAR形变监测方法
CN116168176A (zh) * 2023-04-23 2023-05-26 深圳大学 联合InSAR与激光点云的建筑物几何与变形提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
IBIS-L变形监测中低质量影像的分析去除探讨;骆林喜 等;矿山测量;第45卷(第02期);第33-36页 *
Microdeformation monitoring by permanent scatterer GB-SAR interferometry based on image subset series with short temporal baselines: The Geheyan Dam case study;Peng Wang 等;Measurement(第184期);第1-16页 *
基于GB-SAR的建筑物微变形测量研究;周校 等;测绘地理信息;第37卷(第05期);第40-43页 *

Also Published As

Publication number Publication date
CN116908789A (zh) 2023-10-20

Similar Documents

Publication Publication Date Title
US20200103530A1 (en) Method for extracting elevation control point with assistance of satellite laser altimetry data
CN113091599B (zh) 融合无人机dom和星载sar影像的地表三维形变提取方法
CN104820217B (zh) 一种多法向平面的多元线阵探测成像激光雷达的检校方法
CN102034238B (zh) 基于光学成像测头和视觉图结构的多摄像机系统标定方法
CN108983232B (zh) 一种基于邻轨数据的InSAR二维地表形变监测方法
CN112904337B (zh) 一种基于Offset Tracking技术的边坡形变时序监测方法
CN103093459B (zh) 利用机载LiDAR点云数据辅助影像匹配的方法
CN109633724B (zh) 基于单星与多地面站联合测量的无源目标定位方法
CN105929393B (zh) 大型人工建筑物三维形变的雷达测量方法及测量装置
CN109613583A (zh) 基于单星与地面站测向及联合测时差的无源目标定位方法
KR101804522B1 (ko) 다중 변위 추정 커널을 이용한 위성레이더 오프셋트래킹 장치 및 그 방법
CN106525054B (zh) 一种采用星上推扫遥感图像信息的单星自主测定轨方法
EP3156818B1 (en) Automatic three-dimensional geolocation of sar targets and simultaneous estimation of tropospheric propagation delays using two long-aperture sar images
CN105866777B (zh) 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法
US20230350067A1 (en) Robust laser scanning for generating a 3d model
CN103744081A (zh) 一种机载圆迹合成孔径雷达高精度三维成像补偿方法
CN114608465B (zh) 基于计算机视觉的土木工程结构多测点空间位移实时测量方法
CN105783754A (zh) 基于三维激光扫描的GBInSAR三维位移场提取方法
Yu et al. Automatic extrinsic self-calibration of mobile LiDAR systems based on planar and spherical features
CN113176544B (zh) 一种边坡雷达影像与地形点云的误匹配校正方法
CN116908789B (zh) 地基合成孔径雷达干涉测量建筑立面变形信息提取方法
CN105403886A (zh) 一种机载sar定标器图像位置自动提取方法
CN113124821B (zh) 一种基于曲面镜和平面镜的结构测量方法
RU2406071C1 (ru) Способ навигации движущихся объектов
CN108957500B (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