CN113176544B - 一种边坡雷达影像与地形点云的误匹配校正方法 - Google Patents
一种边坡雷达影像与地形点云的误匹配校正方法 Download PDFInfo
- Publication number
- CN113176544B CN113176544B CN202110244384.1A CN202110244384A CN113176544B CN 113176544 B CN113176544 B CN 113176544B CN 202110244384 A CN202110244384 A CN 202110244384A CN 113176544 B CN113176544 B CN 113176544B
- Authority
- CN
- China
- Prior art keywords
- radar
- point
- sub
- terrain
- coordinate
- 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
- 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/40—Means for monitoring or calibrating
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Abstract
本发明公开了一种边坡形变监测雷达影像与地形点云基于少量人工控制点的误匹配校正方法,根据少量人工控制点先初步地估计偏差模式,作为最小二乘迭代初值,利用几何映射方法完成粗匹配。进而利用雷达图像中的明显地物和周围的相关性分析提取更多的地面控制点,利用大量控制点优化最小二乘变换模型参数完成精匹配实现误匹配校正。本发明仅需至少3个人工控制点即可完成边坡形变监测雷达影像与地形点云的误匹配纠正,具有计算量小,达像元级匹配精度,易于编程实现等优点,所述方法不受雷达视角和设站位置的影响,可有效地定位监测数据显示的异常形变区域,对边坡工程安全监测,遏制重大安全生产事故发生具有重要的意义。
Description
技术领域
本发明涉及一种边坡形变监测雷达影像与地形点云基于少量人工控制点的误匹配校正方法,属于露天矿边坡形变监测技术领域。
背景技术
近年来,全世界地质灾害多发频发,各种地质原因和人类活动等因素均可能会诱发区域形变并衍生灾害。地基边坡形变监测雷达(GB-InSAR)技术越来越多地应用于露天矿边坡、水利水电工程、自然山体滑坡、崩塌、泥石流灾害监测场景。直线重复轨道雷达更是以其高时空分辨力和高形变测量精度,可连续大面积监测,受雨雪雾恶劣天气影响较小等优势,在各类边坡形变监测中得到日益广泛的应用。雷达在二维极坐标空间构像,对于不熟悉雷达成像原理的监测人员存在不直观的问题,实际使用需变换至熟悉的几何以直观定位形变区。
现有的技术方案用来实现雷达图像和地形点云匹配主要有平面相似变换法和几何映射方法。平面相似变换法为星载雷达方法在地基雷达领域应用的延伸,需要雷达图像精确的正射纠正,但常常难以获得纠正参数。几何映射通过测量轨道信息和成像参数,数据获取简单得到广泛应用。然而由于现场条件的特殊性,单独依靠几何映射方法匹配结果会发生较大偏差。目标边坡多在持续变形极度危险,使得人员难以深入布设大量地面控制点校正偏差,因此亟需依靠少量人工控制点校正几何映射匹配偏差。现已提出的雷达图像和地形点云匹配方法多是测量直线轨道端点然后依据成像参数进行几何映射的粗匹配方法,无法应对边坡监测准确定位异常形变区的目标。
发明内容
为克服现有技术的不足,本发明提供一种基于少量人工目标控制、显著目标识别增选控制点和最小二乘优化的边坡形变监测雷达影像与地形点云误匹配校正方法。
本发明为解决上述技术问题采用以下技术方案:
一种边坡雷达影像与地形点云的误匹配校正方法,包括以下具体步骤:
步骤2,使用几何映射方法完成雷达图像和三维地形点云数据粗匹配,得到粗匹配映射表Trough;
步骤3,选取雷达图像中若干人工目标,测量得到人工目标的二维坐标集合{d1(r,θ),…,dn(r,θ)},并定位人工目标在三维地形点云数据中的三维坐标集合{d1(X,Y,Z)T,…,dn(X,Y,Z)T}以作为人工目标的真实坐标;
步骤4,查找粗匹配映射表Trough获得{d1(r,θ),…,dn(r,θ)}在粗匹配映射表中匹配的三维坐标集合{d1(X,Y,Z)S,…,dn(X,Y,Z)S}作为人工目标的偏差坐标;
步骤5,将人工目标的真实坐标和偏差坐标输入空间坐标变换方程,得到空间坐标变换方程的变换参数初值;
步骤6,进行显著地物和周边区域的相关性分析提取控制点对;
步骤7,利用步骤5中的变换参数初值、步骤6中的控制点对及最小二乘迭代估计空间坐标变换方程的变换参数的最优值,对粗匹配映射表Trough进行变换完成误匹配校正。
进一步,步骤1中天线行迹矢量为:其中,为雷达轨道中轴向量,Δh为俯仰调节旋钮中心距离轨道轴心高度,Lantenna为天线中心至俯仰调节主轴距离,θsight为雷达视主轴与水平面夹角,和分别为以雷达合成孔径中心建立的坐标系Os-XsYsZs中Ys和Zs轴方向的空间单位向量。
进一步,坐标系Os-XsYsZs中XsOsYs为水平面,Zs为铅垂线反方向,Xs轴方向为合成孔径波束的指向在水平面的投影,Ys轴与Xs轴垂直并按右手坐标系定出。
进一步,步骤2具体为:
通过插值搜索方法获得取雷达图像I(r,θ)中插值点的像素索引,得到粗匹配映射表Trough,完成粗匹配。
进一步,Pterrain中第i个顶点Ai相对于孔径中心Os的斜距:
以相对竖直中心平面左侧为方位负角度,Pterrain中第i个顶点Ai相对于孔径中心Os的方位角为:
进一步,步骤5中空间坐标变换方程为:其中为真实坐标,[TX TY TZ]T为平移因子,λ为缩放因子,Rx(ωX)、Ry(ωY)、Rz(ωZ)为各坐标轴的旋转矩阵,ωX、ωY、ωZ为各坐标轴的偏转角,为偏差坐标;空间坐标变换方程的变换参数为平移因子、缩放因子和各坐标轴的偏转角。
进一步,步骤6具体为:
提取雷达图像内带有监测场景内的显著地物的子图像Isub以及三维地形点云数据中包含这些显著地物的子点云Psub,查找Trough内Isub图像中心点匹配的三维坐标dsub(X,Y,Z)T;
计算Psub各点相对于雷达的斜距和方位角,斜距和方位角二维网格化生成网格,使用雷达入射角的经验散射模型计算网格点反射率,生成Psub的反射率图Iσ,获取Psub各点与Iσ中各点的映射表Tpcl;
使用相关性分析方法分析Isub和Iσ,提取Isub图像中心坐标与相关性峰值点坐标Isub(r,θ);
查找映射表Tpcl内Isub(r,θ)对应的三维坐标dsub(X,Y,Z)S,dsub(X,Y,Z)T与dsub(X,Y,Z)S形成一组非人工设置的控制点对;
重复上述步骤,获取包括若干控制点对的集合S。
进一步,步骤7具体为:
利用步骤5中的变换参数初值、步骤6中的控制点,使用最小二乘迭代估计空间坐标变换方程的变换参数的最优值,将变换参数的最优值代入空间坐标变换方程形成变换方程f2;
将三维地形点云数据Pterrain代入变换方程f2得到误匹配校正参考点集Ptransform,两者之间的映射关系记为gtransform;
利用gtransform遍历Pterrain与Ptransform中各点的最近欧式距离,形成最近欧式距离vdist向量及最近欧式距离索引vidx向量;
剔除Ptransform内最近欧式距离超过雷达图像像元空间分辨单元大小的点,保证误匹配纠正结果为像元级,将粗匹配结果Prough内每点的形变值赋给Pterrain内索引为vidx的点,形成Pcorrected完成误匹配纠正。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)利用辅助量测设备结合雷达设计参数计算天线行迹矢量,避免了常规方法仅仅测量轨道两端坐标造成的几何映射较大误差;
(2)通过场景内少量人工目标的匹配偏差确定空间变换方程参数初值,解决了空间变换方程初值选取困难,常规方法参数置零造成迭代难以收敛的问题;
(3)通过少量人工目标控制,利用场景内显著地物与周边区域的相关性分析提取大量控制点,算法简单,易于实现,解决了目标区域持续变形缺少足够的人工目标难题;
(4)通过大量控制点最小二乘迭代优化变换方程参数,并加以匹配距离约束,可达像元级匹配校正精度。
附图说明
图1是以雷达理想孔径中心为原点建立的用于匹配的测量观测坐标系示意图;
图2是边坡形变监测雷达系统与斜坡地形数据俯视视角几何关系模型示意图;
图3是本发明实施例提供的一种边坡形变监测雷达影像与地形点云基于少量人工控制点的误匹配校正方法的流程图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明一种边坡形变监测雷达影像与地形点云基于少量人工控制点的误匹配校正方法,如图3所示,包括以下步骤:
步骤1,利用辅助量测设备测量雷达轨道两端点坐标P1、P2(即雷达行进起始点和终止点限位坐标),结合雷达的设计参数计算天线行迹矢量,使用几何映射方法完成雷达图像和地形数据粗匹配。所述雷达的设计参数包括合成孔径长度Ls、步进间隔ΔL、方位向采样点数m、俯仰调节旋钮中心距离轨道轴心高度Δh、天线中心至俯仰调节主轴距离Lantenna。
步骤2,使用几何映射方法完成雷达图像和地形数据粗匹配,得到粗匹配映射表Trough。
天线行迹矢量作为几何映射算法输入参数,利用常规距离-多普勒算法计算地形点云数据Pterrain中各点与雷达的相对斜距R3D和相对方位角θ3D,通过常规最邻近搜索获得取雷达形变图像I(r,θ)中最近邻近像素索引(R2D,θ2D),得到粗匹配映射表Trough,完成粗匹配。得到粗匹配结果Prough,Prough包括坐标点和每点对应形变值。
步骤3,依据场景内的雷达和地形点云均可精确量测的少量人工目标的匹配偏差确定空间变换方程参数迭代初值,并进行显著地物和周边区域的相关性分析提取更多控制点。
通过识别雷达图像中为高亮度像元定位场景内具有少量人工的雷达波的强反射目标,获得二维坐标d1(r,θ),…,dn(r,θ)。
通过辅助量测设备定位人工目标在三维地形点云数据中的三维坐标d1(X,Y,Z)T,…,dn(X,Y,Z)T作为真实坐标。
查找粗匹配映射表Trough,获得二维坐标d1(r,θ),…,dn(r,θ)在粗匹配映射表中匹配的d1(X,Y,Z)S,…,dn(X,Y,Z)S作为人工目标的偏差坐标。
将人工目标真实坐标和偏差坐标输入空间坐标变换方程f:获得七个空间坐标变换方程参数的初值λ0,方程中为真实坐标,λ为缩放因子,[TX TY TZ]T为平移因子,Rx(ωX)、Ry(ωY)、Rz(ωZ)为坐标系间各坐标轴的旋转矩阵,ωX、ωY、ωZ为偏转角,为偏差坐标。
提取雷达图像内带有监测场景内的显著地物的子图像Isub以及地形点云内包含这些显著地物的子点云Psub,查找Trough内Isub图像中心点匹配的三维坐标dsub(X,Y,Z)T。
计算Psub各点相对于雷达的斜距和方位角,网格化生成二维网格,使用雷达入射角的经验散射模型计算网格点反射率生成各子点云区的反射率图Iσ,获取Psub各点与Iσ图像网格点的映射表Tpcl。
使用常规相关性分析方法分析Isub和Iσ,提取Isub图像中心坐标与相关性峰值点坐标Isub(r,θ)。
查找映射表Tpcl内Isub(r,θ)对应的三维坐标dsub(X,Y,Z)S,集合dsub(X,Y,Z)T与dsub(X,Y,Z)S形成一组识别的非人工设置的控制点对。
重复上述步骤,获取识别大量控制点对集合S。
步骤4,利用空间坐标变换方程参数初值、控制点对及最小二乘迭代估计空间变换方程最优参数,对几何映射方法的粗匹配结果进行变换完成误匹配校正。
将Pterrain带入变换方程f2得到误匹配校正参考点集Ptransform,映射关系记为gtransform。
利用映射gtransform遍历Pterrain与Ptransform各点的最近欧式距离,形成最近欧式距离vdist向量及最近欧式距离索引vidx向量。
查找vidx,剔除vdist超过雷达图像像元空间分辨单元大小的点,保证误匹配纠正结果为像元级,将Prough每点形变量赋给Pterrain内索引为vidx点,形成Pcorrected完成误匹配纠正。
实施例
图1为本发明实施例提供的以雷达理想孔径中心为原点建立的用于匹配的测量观测坐标系示意图
利用全站仪、三维激光扫描、静态全球卫星导航定位系统(GPS/GNSS)等常用辅助量测设备测量雷达轨道两端坐标P1、P2,P1、P2分别为雷达行进起始点和终止点限位坐标,共观测np次,np≥3,第一次测量结果为第np次测量结果为平均测量结果
所述雷达的设计参数,包括合成孔径长度Ls、步进间隔ΔL、方位向采样点数m、俯仰调节旋钮中心距离轨道轴心高度Δh、天线中心至俯仰调节主轴距离Lantenna。
天线行迹矢量计算:其中,为雷达轨道中轴向量。具体实施时通过对雷达方位向采样点坐标xm线性拟合,然后在向量投影计算,xm=mΔL-Ls/2,m=0,1,2,...,M-1,M为雷达方位向采样总点数。θsight为雷达视主轴与水平面夹角,分别为以雷达合成孔径中心建立的坐标系各轴空间单位向量,所述坐标系Os-XsYsZs为地基雷达架设施工使用的测量坐标系,如图1所示,孔径中心Os,XsOsYs为水平面,Zs为铅垂线反方向,Xs方向为合成孔径波束的指向在水平面的投影,Ys轴与Xs垂直并按右手坐标系定出。
如图2,为本发明实施例提供的边坡形变监测雷达系统与斜坡地形数据俯视视角几何关系模型示意图。借助斜坡地形点云Pterrain的各点量化几何映射方法。遍历Pterrain中每个顶点相对于孔径中心Os的斜距:
遍历地形点云Pterrain中每个顶点相对于孔径中心Os的方位角,此处以相对竖直中心平面左侧为方位负角度,则方位角为:
然后通过常规双线性插值方法即可获得三维空间点在二维雷达图像中的像素位置,将对应像素的形变测量结果赋值给三维空间点即实现了本实施例的几何映射粗匹配方法,得到粗匹配映射表Trough,完成粗匹配。
通过识别雷达图像中为高亮度像元定位场景内具有少量人工目标如三面角反射器、金属靶球、金属箱、混凝土体等雷达波的强反射目标,获得二维坐标集合{d1(r,θ),…,dn(r,θ)},人工目标数量n≥3。
通过全站仪、三维激光扫描、静态GPS/GNSS等辅助量测设备定位人工目标在三维地形点云数据中的三维坐标集合{d1(X,Y,Z)T,…,dn(X,Y,Z)T},作为真实坐标。
进一步,查找粗匹配映射表Trough,获得二维坐标{d1(r,θ),…,dn(r,θ)}在粗匹配映射表中匹配的{d1(X,Y,Z)S,…,dn(X,Y,Z)S}作为人工目标的偏差坐标。
将人工目标真实坐标和偏差坐标输入空间坐标变换方程f,本发明实施例以空间右手直角坐标系为例,首先将坐标绕Xs轴旋转ωX,得到旋转矩阵Rx(ωX);再将坐标轴绕Ys轴旋转ωY,得旋转矩阵Ry(ωY);最后将坐标轴绕Zs轴旋转ωZ,得到旋转矩阵Rz(ωZ)。三维坐标转换模型方程为:
三个人工目标的真实坐标集合和偏差坐标集合带入方程f,可求得七个变换参数的初值 λ0,方程f中为真实坐标,λ为缩放因子,[TX TY TZ]T为平移因子,Rx(ωX)、Ry(ωY)、Rz(ωZ)为坐标系间各坐标轴的旋转矩阵,ωX、ωY、ωZ为偏转角,为偏差坐标。
方程f参数优化需要更多控制点,本发明实施例提取雷达图像内带有监测场景内的显著地物如建筑、金属房、巨型岩石、挡墙等的子图像Isub以及地形点云内包含这些显著地物的子点云Psub,查找Trough内Isub图像中心点匹配的三维坐标dsub(X,Y,Z)T。
计算Psub各点相对于雷达的斜距和方位角,网格化生成二维网格,使用雷达入射角的经验散射模型,常用的经验散射模型有余弦模型、余弦平方模型,本实施例选取经验模型:σiδ为二维网格点的值,ηiδ代表雷达波入射角。本实施例可通过对Psub各点与邻近的k个点拟合小面元平面D,计算Os与Psub各点所成向量与小面元D的法向量夹角求取雷达波入射角,k≥6。
计算网格点反射率生成各子点云区的反射率图Iσ,获取Psub各点与Iσ图像网格点的映射表Tpcl。
使用归一化互相关性分析方法分析Isub和Iσ,:
式中,C为互相关系数,在0.1~1之间;Iσ(x,y)为图像Iσ网格点(x,y)的灰度值;Isub(x,y)为Isub图像像点(x,y)的灰度值;μT为Iσ灰度均值;μI为图像Isub的灰度均值。提取Isub图像中心坐标与相关性峰值点坐标Isub(r,θ)。
查找映射表Tpcl内Isub(r,θ)对应的三维坐标dsub(X,Y,Z)S,集合dsub(X,Y,Z)T与dsub(X,Y,Z)S形成一组识别的非人工设置的控制点对。重复上述步骤,获取识别大量控制点对集合S。
线性化变换公式(5)得到误差方程:
最小二乘迭代令公式(6)误差最小。通过迭代计算控制舍入误差,舍入误差选取单位权中误差自由度为误差方程数量与实际求解方程所需变量数的差,各控制点视为独立观测量权阵P取单位矩阵,求出七个变换参数的最佳估值。实施例迭代计算过程:
将参数初值带入XT,R′,组成式6误差方程。n个控制点将组成3n个误差方程。
利用最小二乘方法求取变换参数改正数x(∈+1)(∈代表迭代次数)④检核参数的改正数是否小于给定的限差,本实施例取平移量mm,偏转角10-3″和缩放系数10-7为限差,若未达到限差则将X(∈)=X(∈-1)+x(∈)作为新的初值,重复迭代;符合限差,则计算结束,将X∈作为参数最佳估值,获得优化的最终变换参数带回f形成变换方程f2。
将Pterrain带入变换方程f2得到误匹配校正参考点集Ptransform,映射关系记为gtransform。
利用映射gtransform遍历Pterrain与Ptransform各点最近欧式距离。形成距离vdist向量及最近邻索引vidx向量。本实施例引入常用的点云结果化处理方法k-d树(k dimensionaltree)和k近邻(k-nearest neighbor,KNN)聚类方法。在处理过程中,利用KNN查找k个最接近的点,令k=1,即求得Ptransform与Pterrain中最邻近点。
查找vidx,剔除Ptransform内最近欧式距离超过雷达图像像元空间分辨单元大小的点,本发明实施例空间分辨单元大小指雷达二维图像像元匹配至三维空间距离向和方位向覆盖范围的平均长度,小于该长度的匹配点将保留,保证误匹配纠正结果为像元级,将Prough每点形变量赋给Pterrain内索引为vidx点,形成Pcorrected完成误匹配纠正。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (10)
1.一种边坡雷达影像与地形点云的误匹配校正方法,其特征在于,包括以下具体步骤:
步骤2,使用几何映射方法完成雷达图像和三维地形点云数据粗匹配,得到粗匹配映射表Trough;
步骤3,选取雷达图像中若干人工目标,测量得到人工目标的二维坐标集合{d1(r,θ),…,dn(r,θ)},并定位人工目标在三维地形点云数据中的三维坐标集合{d1(X,Y,Z)T,…,dn(X,Y,Z)T}以作为人工目标的真实坐标;其中,n为人工目标数量;
步骤4,查找粗匹配映射表Trough获得{d1(r,θ),…,dn(r,θ)}在粗匹配映射表中匹配的三维坐标集合{d1(X,Y,Z)S,…,dn(X,Y,Z)S}作为人工目标的偏差坐标;
步骤5,将人工目标的真实坐标和偏差坐标输入空间坐标变换方程,得到空间坐标变换方程的变换参数初值;
步骤6,进行显著地物和周边区域的相关性分析提取控制点对;
步骤7,利用步骤5中的变换参数初值、步骤6中的控制点对及最小二乘迭代估计空间坐标变换方程的变换参数的最优值,对粗匹配映射表Trough进行变换完成误匹配校正。
3.如权利要求2所述的一种边坡雷达影像与地形点云的误匹配校正方法,其特征在于,坐标系Os-XsYsZs中XsOsYs为水平面,Zs为铅垂线反方向,Xs轴方向为合成孔径波束的指向在水平面的投影,Ys轴与Xs轴垂直并按右手坐标系定出。
8.如权利要求1所述的一种边坡雷达影像与地形点云的误匹配校正方法,其特征在于,步骤6具体为:
提取雷达图像内带有监测场景内的显著地物的子图像Isub以及三维地形点云数据中包含这些显著地物的子点云Psub,查找Trough内Isub图像中心点匹配的三维坐标dsub(X,Y,Z)T;
计算Psub各点相对于雷达的斜距和方位角,斜距和方位角二维网格化生成网格,使用雷达入射角的经验散射模型计算网格点反射率,生成Psub的反射率图Iσ,获取Psub各点与Iσ中各点的映射表Tpcl;
使用相关性分析方法分析Isub和Iσ,提取Isub图像中心坐标与相关性峰值点坐标Isub(r,θ);
查找映射表Tpcl内Isub(r,θ)对应的三维坐标dsub(X,Y,Z)S,dsub(X,Y,Z)T与dsub(X,Y,Z)S形成一组非人工设置的控制点对;
重复上述步骤,获取包括若干控制点对的集合S。
10.如权利要求1所述的一种边坡雷达影像与地形点云的误匹配校正方法,其特征在于,步骤7具体为:
利用步骤5中的变换参数初值、步骤6中的控制点,使用最小二乘迭代估计空间坐标变换方程的变换参数的最优值,将变换参数的最优值代入空间坐标变换方程形成变换方程f2;
将三维地形点云数据Pterrain代入变换方程f2得到误匹配校正参考点集Ptransform,两者之间的映射关系记为gtransform;
利用gtransform遍历Pterrain与Ptransform中各点的最近欧式距离,形成最近欧式距离Vdist向量及最近欧式距离索引Vidx向量;
剔除Ptransform内最近欧式距离超过雷达图像像元空间分辨单元大小的点,保证误匹配纠正结果为像元级,将粗匹配结果Prough内每点的形变值赋给Pterrain内索引为Vidx的点,形成Pcorrected完成误匹配纠正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110244384.1A CN113176544B (zh) | 2021-03-05 | 2021-03-05 | 一种边坡雷达影像与地形点云的误匹配校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110244384.1A CN113176544B (zh) | 2021-03-05 | 2021-03-05 | 一种边坡雷达影像与地形点云的误匹配校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113176544A CN113176544A (zh) | 2021-07-27 |
CN113176544B true CN113176544B (zh) | 2022-11-11 |
Family
ID=76921856
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110244384.1A Active CN113176544B (zh) | 2021-03-05 | 2021-03-05 | 一种边坡雷达影像与地形点云的误匹配校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113176544B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113740844B (zh) * | 2021-09-09 | 2024-04-02 | 雷添杰 | 面向坝体三维形变监测的两台地基雷达联合观测方法 |
CN115146745B (zh) * | 2022-09-01 | 2022-12-02 | 深圳市城市公共安全技术研究院有限公司 | 点云数据坐标点位置校正的方法、装置、设备及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770027B (zh) * | 2010-02-05 | 2012-05-16 | 河海大学 | 基于InSAR与GPS数据融合的地表三维形变监测方法 |
CN105678689B (zh) * | 2015-12-31 | 2020-01-31 | 百度在线网络技术(北京)有限公司 | 高精地图数据配准关系确定方法及装置 |
CN105717527B (zh) * | 2016-02-18 | 2017-11-17 | 河海大学 | 一种利用边坡变形数据快速确定监测点移动轨迹的方法 |
CN108828606B (zh) * | 2018-03-22 | 2019-04-30 | 中国科学院西安光学精密机械研究所 | 一种基于激光雷达和双目可见光相机联合测量方法 |
CN110415342B (zh) * | 2019-08-02 | 2023-04-18 | 深圳市唯特视科技有限公司 | 一种基于多融合传感器的三维点云重建装置与方法 |
CN112102458A (zh) * | 2020-08-31 | 2020-12-18 | 湖南盛鼎科技发展有限责任公司 | 基于激光雷达点云数据辅助的单镜头三维图像重构方法 |
-
2021
- 2021-03-05 CN CN202110244384.1A patent/CN113176544B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113176544A (zh) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li | Potential of high-resolution satellite imagery for national mapping products | |
Dall’Asta et al. | Unmanned Aerial Systems and DSM matching for rock glacier monitoring | |
CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
CN113176544B (zh) | 一种边坡雷达影像与地形点云的误匹配校正方法 | |
CN104952107A (zh) | 基于车载LiDAR点云数据的桥梁三维重建方法 | |
Hu et al. | DEM extraction from worldview-3 stereo-images and accuracy evaluation | |
Glira et al. | Hybrid orientation of airborne lidar point clouds and aerial images | |
Toutin | Comparison of 3D physical and empirical models for generating DSMs from stereo HR images | |
Kumari et al. | Adjustment of systematic errors in ALS data through surface matching | |
Eltner et al. | Integrated processing of high resolution topographic data for soil erosion assessment considering data acquisition schemes and surface properties | |
Lee et al. | A study on scan data matching for reverse engineering of pipes in plant construction | |
CN113358091A (zh) | 一种利用三线阵立体卫星影像生产数字高程模型方法 | |
CN114689015A (zh) | 一种提高光学卫星立体影像dsm高程精度的方法 | |
Toutin | Radarsat-2 DSM generation with new hybrid, deterministic, and empirical geometric modeling without GCP | |
Zhang et al. | Generation of digital surface model from high resolution satellite imagery | |
RU2406071C1 (ru) | Способ навигации движущихся объектов | |
Zhang et al. | Matching of Ikonos stereo and multitemporal GEO images for DSM generation | |
Lee et al. | Point-to-surface matching for DEM correction using icesat data | |
Tian et al. | Automatic calibration method for airborne lidar systems based on approximate corresponding points model | |
Shin et al. | A study on airborne LiDAR calibration and operation techniques for bathymetric survey | |
Zhao et al. | Comparison and analysis of accuracy of elevation extraction based on the ZY-3 01 and 02 satellites stereoscopic images | |
CN102735225B (zh) | 月球控制网的建立方法 | |
Pan et al. | A general method of generating satellite epipolar images based on RPC model | |
Fang et al. | Beach topography surveying using unmanned aerial vehicle photogrammetry technology | |
Trisakti et al. | DEM generation from stereo ALOS PRISM and its quality improvement |
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 |