CN109631863A - 一种空地结合的潮间带一体化测绘方法 - Google Patents
一种空地结合的潮间带一体化测绘方法 Download PDFInfo
- Publication number
- CN109631863A CN109631863A CN201910056969.3A CN201910056969A CN109631863A CN 109631863 A CN109631863 A CN 109631863A CN 201910056969 A CN201910056969 A CN 201910056969A CN 109631863 A CN109631863 A CN 109631863A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- axis
- coordinate
- data
- plane
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
Abstract
本发明提供了一种空地结合的潮间带一体化测绘方法,通过原始数据采集下载后,进行点云时空融合,完成从仪器坐标系到地心地固坐标系的转换;通过获取单波束点云与激光扫描仪点云、激光扫描仪点云与机载激光探测与测量系统点云数据之间的点、线、面等几何特征来建立对应关系,求解旋转矩阵和平移量进行点云匹配。本方法集成了机载激光探测与测量系统和拖曳测量系统,可以获取潮间带区域水下地形和水上地形数据,解决了潮间带地形测量问题。两者能完美覆盖测量“盲区”,实现潮间带区域一体化无缝测量。
Description
技术领域
本发明属于海洋测绘技术领域,具体涉及以一种空地结合的方式进行潮间带区域水上和水下地形快速获取的测绘方法。
背景技术
潮间带属于海岸带的一部分,易受到潮汐的影响,测量条件比较困难。当前,海岸带测绘技术主要有人工实地测量、船载测量和航测遥感等方式。自然海岸带通常存在礁石、滩涂、湿地等特殊地理环境,以传统人工实地测量和船载测量方式难以高效地进行作业,甚至存在危及人身安全的作业风险,无法保障海岸带基础地理信息的快速更新。航测遥感技术手段具有范围广、频度高和实时的优点,被广泛应用于海洋测绘,但是由于受到卫星重复周期、海岸带云雨天气、分辨率低等影响,难以满足动态监测和水下地形的高精度测量需求,对潮水覆盖区域尤其是大面积浅滩和宽阔滩涂区域无能为力。
船载多传感器水岸一体化综合测量技术的发展,为海岸带一体化测量提供了依据。但当海岸带地形以滩涂为主且潮差较小时,由于测量船和测深仪存在吃水问题使得基于船载激光测量无法对潮间带进行扫测,出现潮间带的测量“盲区”。
机载LiDAR技术具有高便捷、无控制、高精度、高密集、高效率、高分辨率、多回波等特点,作为精确、快速地获取地面三维数据的工具已得到广泛的认同,但由于潮间带地形特殊性,激光信号受地物遮挡、水面反射等影响,也会出现测量“盲区”。
发明内容
针对现有方法中存在的问题,本发明提出一种空地结合的潮间带一体化测绘方法,机载激光探测与测量系统结合拖曳一体化测量系统,解决了潮间带常见的测量问题,具有良好的效果。
本发明采用以下的技术方案:
一种空地结合的潮间带一体化测绘方法,包括测量装置,测量装置包括无人机和载体车,无人机上搭载激光探测与测量系统,载体车上设置有激光扫描仪、单波束测深仪、惯性导航系统和卫星定位系统;
测绘方法包括以下步骤:
步骤1:原始数据采集下载后,检查各测量仪器观测数据的完整性、正确性,进行观测数据文件的格式转换,进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,通过点云数据预处理,形成标准LAS数据文件;
步骤2:进行点云时空融合,完成从仪器坐标系到地心地固坐标系的转换;
步骤3:通过获取单波束点云与激光扫描仪点云、激光扫描仪点云与机载LIDAR点云数据之间的点、线、面等几何特征来建立对应关系,求解旋转矩阵和平移量进行点云匹配。
优选地,所述步骤1包括以下子步骤:
步骤1.1:利用硬件同步控制器进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,硬件同步控制器的时间保持与卫星定位系统的时间同步;硬件同步控制器控制各测量仪器的触发及监测其反馈信号,信号的传输都采用差分方式;
步骤1.2:利用IE进行基准站和流动站数据结算,输出POS数据;
步骤1.3:利用VSursPROCESS软件进行点云数据预处理,形成标准LAS数据文件。
优选地,所述步骤2包括以下子步骤:
步骤2.1:从激光扫描仪坐标系到载体坐标系转换;
根据标定的公共点计算坐标系转换需要的的6参数:lx、ly、lz、ω、κ;
lx、ly、lz分别为坐标转换过程中点坐标沿载体坐标系x轴、y轴、z轴的平移量;
ω、κ分别为坐标转换过程中点坐标绕载体坐标系x轴、y轴、z轴的旋转角度;
设激光扫描仪坐标系下点坐标为在载体坐标系下的坐标为则有,
其中,代表从激光扫描仪坐标系到载体坐标系转换的旋转矩阵,Rx(ω)、Rz(κ)分别为绕载体坐标系x轴、y轴、z轴旋转的旋转矩阵;
步骤2.2:从载体坐标系到站心坐标系转换;
惯导坐标系前进方向为y轴、向右为x轴、向上为z轴;
惯导记录姿态角:侧滚角Roll、俯仰角Pitch、偏航角Heading;
侧滚角Roll:惯导x轴与水平方向之间的夹角,载体右侧向下为正;
俯仰角Pitch:惯导y轴与水平方向之间的夹角,载体向上为正;
偏航角Heading:惯导前进方向,即xy平面与正北方向之间的夹角,顺时针为正;
设Roll、Pitch、Heading分别为r、p、y;设点在站心坐标系下的坐标为载体坐标系下坐标向站心坐标系转换;
①先绕z轴旋转y;
②再绕x轴旋转p;
③最后绕y轴旋转r;
则有,
其中,代表从载体坐标系到站心坐标系转换的旋转矩阵,Rx(p)、Ry(r)、Rz(y)分别为绕站心坐标系x轴、y轴、z轴旋转的旋转矩阵;
步骤2.3:从站心坐标系到地心地固坐标系转换;
站心坐标系原点在WGS84下的经纬度分别为L和B;
设扫描点在地心地固坐标系下的坐标为站心坐标系下的坐标转换为WGS84坐标系下的坐标:
①先绕x轴旋转
②再绕z轴旋转
③最后将站心坐标系原点平移到WGS84坐标系原点;
则有,
其中,代表从站心坐标系到地心地固坐标系转换的旋转矩阵;
其中,N代表卯酉圈曲率半径,h代表大地高,a、b代表椭球长、短半径;
得
为站心坐标系原点在地心地固坐标系下的空间直角坐标;
最后得,
优选地,所述步骤3包括:
首先介绍一下所需要的理论基础:
平面拟合原理
已知平面方程:a′x+b′y+c′z=d (6)且a′2+b′2+c′2=1,d≥0,(x y z)是平面上的任意点;
其中,(a′ b′ c′)代表垂直于平面且远离坐标原点方向的方向余弦;
同样的,它也表示垂直于平面远离坐标原点的单位法向量的分量;
d表示了平面和坐标原点间的垂直距离;
最小二乘算法基本原理
最小二乘法的思想就是估算一个最佳数学模型,求得最佳模型参数使估算模型点的值和实际观察值的差值的平方和最小;
在二维数据中,假设给出模型函数g(xi),g(xi)与真实值(xi,yi)的误差可用三种形式表述:
①误差最大值,即其中n为点总个数;
②误差绝对值之和,即
③误差平方和,即
最小二乘算法的原理便是寻找一个最佳参数模型,对于给定的数据(xi,yi)(i=0,1…n),误差平方和最小;
最小二乘平面拟合算法
基于最小二乘法的平面拟合算法是在三维点云数据的基础上进行的拟合算法,以公式(6)为三维模型,确定a′,b′,c′,d四个参数使三维模型的估算值和测量值之间的差值平方和E2最小;由于E2是非负数,因此它存在极小值,E2对每个参数的偏导数为零,利用偏导数为零构建的方程组便可求解出a′,b′,c′,d,并使估算值和测量值之间的差值平方和E2最小;
假设点云数据的点集为{pi|pi∈oxyz,i=1,2,3...,N},其中xi和yi设为无误差变化的自变量,zi是包含误差的因变量,则z和x、y的函数相关性可以假设为函数表达式:
z=f(x,y;p,q,r)=p+qx+ry (7)
且其中的p,q,r可以通过最小二乘算法得到;
根据公式(6)、(7),令d=c′p,a′=-qc′,b′=-rc′且则公式(7)可以描述为平面方程;
根据上述定义f函数的公式,假如自变量的值为xi和yi,被测量值z应该为p+qxi+ryi,然而因为误差的存在测量值z为zi,则第i个点的测量误差为:
ei=f(xi,yi;p,q,r)-zi (8)
根据最小二乘的原理,p,q,r应该满足公式(8)的平方和最小,即E2最小:
若使得E2最小,应满足:
即:
有:
根据公式(11)可以求解出p,q,r,同样的,公式(6)中a′,b′,c′,d四个参数也就可通过d=c′p,a′=-qc′,b′=-rc′,求出;
求解旋转矩阵
在上述算法的计算过程中,提取3对以上对应的平面特征量后,需要求取其各自的平面法向量;
假设求解的一对平面为a1x+b1y+z=d1和a2x+b2y+z=d2,则其法向量对可表示为P=(a1,b1,-1)、Q=(a2,b2,-1),法向量的旋转矩阵R满足P=RQ,采用罗德里格旋转法求解;上述P,Q矩阵,其夹角为:
P,Q矩阵的叉乘为:
若已知单位向量利用法向量将其旋转θ角度后,根据罗德里格旋转法求解其三维旋转矩阵为:
即可求出的三维变换矩阵C1;
求解平移量
利用三维变换矩阵C1可以对点云数据进行旋转,相应的平面特征量对相互平行,在两片点云数据中寻找一对以上的特征点,或利用对应面的的质心坐标的差值,计算点云数据的平移量即:
本发明具有的有益效果是:
本发明提供的一种空地结合的潮间带一体化测绘方法,集成了机载激光探测与测量系统和拖曳测量系统,可以获取潮间带区域水下地形和水上地形数据,解决了潮间带地形测量问题。两者能完美覆盖测量“盲区”,实现潮间带区域一体化无缝测量。
附图说明
图1为本发明提供的潮间带载体车外业采集流程图。
图2为本发明提供的机载激光探测与测量系统外业采集流程图。
图3为本发明提供的单波束数据格式转换流程图。
图4为本发明提供的数据解算流程图。
图5为本发明提供的单波束测深仪数据、激光扫描仪数据、机载激光探测与测量系统数据空间整合流程图。
具体实施方式
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
结合图1至图5,一种空地结合的潮间带一体化测绘方法,包括测量装置,测量装置包括无人机和载体车,无人机上搭载激光探测与测量系统,载体车上设置有激光扫描仪、单波束测深仪、惯性导航系统和卫星定位系统。
测绘方法包括以下步骤:
步骤1:原始数据采集下载后,检查各测量仪器观测数据的完整性、正确性,进行观测数据文件的格式转换,进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,通过点云数据预处理,形成标准LAS数据文件。
各测量仪器观测数据包括激光扫描仪数据、激光探测与测量系统数据、单波束测深仪数据、惯性导航系统数据和卫星定位系统数据等。
具体的,步骤1包括以下子步骤:
步骤1.1:如图1和图2所示,分别为载体车和机载激光探测与测量系统数据采集流程图。图3为单波束数据格式转换流程图。
利用硬件同步控制器进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,它作用是为系统中的激光扫描仪、单波束测深仪等仪器提供相应的同步信号,并记录相应的时间,为后期的数据处理提供一个时间同步基准,硬件同步控制器的时间保持与卫星定位系统的时间同步;硬件同步控制器控制各测量仪器的触发及监测其反馈信号,考虑到环境因素复杂,为保证系统的高可靠性,信号的传输都采用差分方式;
步骤1.2:如图4所示,利用IE进行基准站和流动站数据结算,输出POS数据;
步骤1.3:利用VSursPROCESS软件进行点云数据预处理,形成标准LAS数据文件。
步骤2:进行点云时空融合,完成从仪器坐标系到地心地固坐标系的转换。
步骤2包括以下子步骤:
步骤2.1:从激光扫描仪坐标系到载体坐标系转换;
根据标定的公共点计算坐标系转换需要的的6参数:lx、ly、lz、ω、κ;
lx、ly、lz为坐标转换过程中点坐标沿载体坐标系x轴、y轴、z轴的平移量;
ω、κ为坐标转换过程中点坐标绕载体坐标系x轴、y轴、z轴的旋转角度;
设激光扫描仪坐标系下点坐标为在载体坐标系下的坐标为则有,
其中,代表从激光扫描仪坐标系到载体坐标系转换的旋转矩阵,Rx(ω)、Rz(κ)分别为绕载体坐标系x轴、y轴、z轴旋转的旋转矩阵。
步骤2.2:从载体坐标系到站心坐标系转换;
惯导坐标系前进方向为y轴、向右为x轴、向上为z轴;
惯导记录姿态角:侧滚角Roll、俯仰角Pitch、偏航角Heading。
侧滚角Roll:惯导x轴与水平方向之间的夹角,载体右侧向下为正;
俯仰角Pitch:惯导y轴与水平方向之间的夹角,载体向上为正;
偏航角Heading:惯导前进方向,即xy平面与正北方向之间的夹角,顺时针为正;
设Roll、Pitch、Heading分别为r、p、y;设点在站心坐标系下的坐标为载体坐标系下坐标向站心坐标系转换;
①先绕z轴旋转y;
②再绕x轴旋转p;
③最后绕y轴旋转r;
则有,
其中,代表从载体坐标系到站心坐标系转换的旋转矩阵,Rx(p)、Ry(r)、Rz(y)分别为绕站心坐标系x轴、y轴、z轴旋转的旋转矩阵。
步骤2.3:从站心坐标系到地心地固坐标系转换;
站心坐标系原点在WGS84下的经纬度分别为L和B;
设扫描点在地心地固坐标系下的坐标为站心坐标系下的坐标转换为WGS84坐标系下的坐标:
①先绕x轴旋转
②再绕z轴旋转
③最后将站心坐标系原点平移到WGS84坐标系原点;
则有,
代表从站心坐标系到地心地固坐标系转换的旋转矩阵。
其中,N代表卯酉圈曲率半径,h代表大地高,a、b代表椭球长、短半径。
得
为站心坐标系原点在地心地固坐标系下的空间直角坐标;
最后得,
步骤3:通过获取单波束点云与激光扫描仪点云、激光扫描仪点云与机载LIDAR点云数据之间的点、线、面等几何特征来建立对应关系,求解旋转矩阵和平移量进行点云匹配。
步骤3包括:
首先介绍一下所需要的理论基础:
平面拟合原理
已知平面方程:
a′x+b′y+c′z=d (6)
且a′2+b′2+c′2=1,d≥0,(x y z)是平面上的任意点;
其中(a′ b′ c′)代表垂直于平面且远离坐标原点方向的方向余弦;
同样的,它也表示垂直于平面远离坐标原点的单位法向量的分量;
d表示了平面和坐标原点间的垂直距离。
最小二乘算法基本原理
最小二乘法的思想就是估算一个最佳数学模型,求得最佳模型参数使估算模型点的值和实际观察值的差值的平方和最小;
在二维数据中,假设给出模型函数g(xi),g(xi)与真实值(xi,yi)的误差可用三种形式表述:
①误差最大值,即其中n为点总个数;
②误差绝对值之和,即
③误差平方和,即
最小二乘算法的原理便是寻找一个最佳参数模型,对于给定的数据(xi,yi)(i=0,1…n),误差平方和最小;
最小二乘平面拟合算法
基于最小二乘法的平面拟合算法是在三维点云数据的基础上进行的拟合算法,以公式(6)为三维模型,确定a′,b′,c′,d四个参数使三维模型的估算值和测量值之间的差值平方和E2最小;由于E2是非负数,因此它存在极小值,E2对每个参数的偏导数为零,利用偏导数为零构建的方程组便可求解出a′,b′,c′,d,并使估算值和测量值之间的差值平方和E2最小;
假设点云数据的点集为{pi|pi∈oxyz,i=1,2,3...,N},其中xi和yi设为无误差变化的自变量,zi是包含误差的因变量,则z和x、y的函数相关性可以假设为函数表达式:
z=f(x,y;p,q,r)=p+qx+ry (7)
且其中的p,q,r可以通过最小二乘算法得到;
根据公式(6)、(7),令d=c′p,a′=-qc′,b′=-rc′且则公式(7)可以描述为平面方程;
根据上述定义f函数的公式,假如自变量的值为xi和yi,被测量值z应该为p+qxi+ryi,然而因为误差的存在测量值z为zi,则第i个点的测量误差为:
ei=f(xi,yi;p,q,r)-zi (8)
根据最小二乘的原理,p,q,r应该满足公式(8)的平方和最小,即E2最小:
若使得E2最小,应满足:
即:
有:
根据公式(11)可以求解出p,q,r,同样的,公式(6)中a′,b′,c′,d四个参数也就可通过d=c′p,a′=-qc′,b′=-rc′,求出。
求解旋转矩阵
在上述算法的计算过程中,提取3对以上对应的平面特征量后,需要求取其各自的平面法向量;
假设求解的一对平面为a1x+b1y+z=d1和a2x+b2y+z=d2,则其法向量对可表示为P=(a1,b1,-1)、Q=(a2,b2,-1),法向量的旋转矩阵R满足P=RQ,采用罗德里格旋转法求解;上述P,Q矩阵,其夹角为:
P,Q矩阵的叉乘为:
若已知单位向量利用法向量将其旋转θ角度后,根据罗德里格旋转法求解其三维旋转矩阵为:
即可求出的三维变换矩阵C1;
求解平移量
利用三维变换矩阵C1可以对点云数据进行旋转,相应的平面特征量对相互平行,在两片点云数据中寻找一对以上的特征点,或利用对应面的的质心坐标的差值,计算点云数据的平移量即:
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (4)
1.一种空地结合的潮间带一体化测绘方法,其特征在于,包括测量装置,测量装置包括无人机和载体车,无人机上搭载激光探测与测量系统,载体车上设置有激光扫描仪、单波束测深仪、惯性导航系统和卫星定位系统;
测绘方法包括以下步骤:
步骤1:原始数据采集下载后,检查各测量仪器观测数据的完整性、正确性,进行观测数据文件的格式转换,进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,通过点云数据预处理,形成标准LAS数据文件;
步骤2:进行点云时空融合,完成从仪器坐标系到地心地固坐标系的转换;
步骤3:通过获取单波束点云与激光扫描仪点云、激光扫描仪点云与机载LIDAR点云数据之间的点、线、面等几何特征来建立对应关系,求解旋转矩阵和平移量进行点云匹配。
2.根据权利要求1所述的一种空地结合的潮间带一体化测绘方法,其特征在于,所述步骤1包括以下子步骤:
步骤1.1:利用硬件同步控制器进行各测量仪器高精度时间同步,各测量仪器数据的时、空融合,硬件同步控制器的时间保持与卫星定位系统的时间同步;硬件同步控制器控制各测量仪器的触发及监测其反馈信号,信号的传输都采用差分方式;
步骤1.2:利用IE进行基准站和流动站数据结算,输出POS数据;
步骤1.3:利用VSursPROCESS软件进行点云数据预处理,形成标准LAS数据文件。
3.根据权利要求1所述的一种空地结合的潮间带一体化测绘方法,其特征在于,所述步骤2包括以下子步骤:
步骤2.1:从激光扫描仪坐标系到载体坐标系转换;
根据标定的公共点计算坐标系转换需要的的6参数:lx、ly、lz、ω、κ;
lx、ly、lz分别为坐标转换过程中点坐标沿载体坐标系x轴、y轴、z轴的平移量;
ω、κ分别为坐标转换过程中点坐标绕载体坐标系x轴、y轴、z轴的旋转角度;
设激光扫描仪坐标系下点坐标为在载体坐标系下的坐标为则有,
其中,代表从激光扫描仪坐标系到载体坐标系转换的旋转矩阵,Rx(ω)、Rz(κ)分别为绕载体坐标系x轴、y轴、z轴旋转的旋转矩阵;
步骤2.2:从载体坐标系到站心坐标系转换;
惯导坐标系前进方向为y轴、向右为x轴、向上为z轴;
惯导记录姿态角:侧滚角Roll、俯仰角Pitch、偏航角Heading;
侧滚角Roll:惯导x轴与水平方向之间的夹角,载体右侧向下为正;
俯仰角Pitch:惯导y轴与水平方向之间的夹角,载体向上为正;
偏航角Heading:惯导前进方向,即xy平面与正北方向之间的夹角,顺时针为正;
设Roll、Pitch、Heading分别为r、p、y;设点在站心坐标系下的坐标为载体坐标系下坐标向站心坐标系转换;
①先绕z轴旋转y;
②再绕x轴旋转p;
③最后绕y轴旋转r;
则有,
其中,代表从载体坐标系到站心坐标系转换的旋转矩阵,Rx(p)、Ry(r)、Rz(y)分别为绕站心坐标系x轴、y轴、z轴旋转的旋转矩阵;
步骤2.3:从站心坐标系到地心地固坐标系转换;
站心坐标系原点在WGS84下的经纬度分别为L和B;
设扫描点在地心地固坐标系下的坐标为站心坐标系下的坐标转换为WGS84坐标系下的坐标:
①先绕x轴旋转
②再绕z轴旋转
③最后将站心坐标系原点平移到WGS84坐标系原点;
则有,
其中,代表从站心坐标系到地心地固坐标系转换的旋转矩阵;
其中,N代表卯酉圈曲率半径,h代表大地高,a、b代表椭球长、短半径;
得
为站心坐标系原点在地心地固坐标系下的空间直角坐标;
最后得,
4.根据权利要求1所述的一种空地结合的潮间带一体化测绘方法,其特征在于,所述步骤3包括:
首先介绍一下所需要的理论基础:
平面拟合原理
已知平面方程:a′x+b′y+c′z=d (6)
且a′2+b′2+c′2=1,d≥0,(x y z)是平面上的任意点;
其中,(a′ b′ c′)代表垂直于平面且远离坐标原点方向的方向余弦;
同样的,它也表示垂直于平面远离坐标原点的单位法向量的分量;
d表示了平面和坐标原点间的垂直距离;
最小二乘算法基本原理
最小二乘法的思想就是估算一个最佳数学模型,求得最佳模型参数使估算模型点的值和实际观察值的差值的平方和最小;
在二维数据中,假设给出模型函数g(xi),g(xi)与真实值(xi,yi)的误差可用三种形式表述:
①误差最大值,即其中n为点总个数;
②误差绝对值之和,即
③误差平方和,即
最小二乘算法的原理便是寻找一个最佳参数模型,对于给定的数据(xi,yi)(i=0,1…n),误差平方和最小;
最小二乘平面拟合算法
基于最小二乘法的平面拟合算法是在三维点云数据的基础上进行的拟合算法,以公式(6)为三维模型,确定a′,b′,c′,d四个参数使三维模型的估算值和测量值之间的差值平方和E2最小;由于E2是非负数,因此它存在极小值,E2对每个参数的偏导数为零,利用偏导数为零构建的方程组便可求解出a′,b′,c′,d,并使估算值和测量值之间的差值平方和E2最小;
假设点云数据的点集为{pi|pi∈oxyz,i=1,2,3...,N},其中xi和yi设为无误差变化的自变量,zi是包含误差的因变量,则z和x、y的函数相关性可以假设为函数表达式:
z=f(x,y;p,q,r)=p+qx+ry (7)
且其中的p,q,r可以通过最小二乘算法得到;
根据公式(6)、(7),令d=c′p,a′=-qc′,b′=-rc′且则公式(7)可以描述为平面方程;
根据上述定义f函数的公式,假如自变量的值为xi和yi,被测量值z应该为p+qxi+ryi,然而因为误差的存在测量值z为zi,则第i个点的测量误差为:
ei=f(xi,yi;p,q,r)-zi (8)
根据最小二乘的原理,p,q,r应该满足公式(8)的平方和最小,即E2最小:
若使得E2最小,应满足:
即:
有:
根据公式(11)可以求解出p,q,r,同样的,公式(6)中a′,b′,c′,d四个参数也就可通过d=c′p,a′=-qc′,b′=-rc′,求出;
求解旋转矩阵
在上述算法的计算过程中,提取3对以上对应的平面特征量后,需要求取其各自的平面法向量;
假设求解的一对平面为a1x+b1y+z=d1和a2x+b2y+z=d2,则其法向量对可表示为P=(a1,b1,-1)、Q=(a2,b2,-1),法向量的旋转矩阵R满足P=RQ,采用罗德里格旋转法求解;上述P,Q矩阵,其夹角为:
P,Q矩阵的叉乘为:
若已知单位向量利用法向量将其旋转θ角度后,根据罗德里格旋转法求解其三维旋转矩阵为:
即可求出的三维变换矩阵C1;
求解平移量
利用三维变换矩阵C1可以对点云数据进行旋转,相应的平面特征量对相互平行,在两片点云数据中寻找一对以上的特征点,或利用对应面的的质心坐标的差值,计算点云数据的平移量即:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056969.3A CN109631863A (zh) | 2019-01-22 | 2019-01-22 | 一种空地结合的潮间带一体化测绘方法 |
PCT/CN2019/098770 WO2020151213A1 (zh) | 2019-01-22 | 2019-08-01 | 一种空地结合的潮间带一体化测绘方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056969.3A CN109631863A (zh) | 2019-01-22 | 2019-01-22 | 一种空地结合的潮间带一体化测绘方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109631863A true CN109631863A (zh) | 2019-04-16 |
Family
ID=66062911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910056969.3A Pending CN109631863A (zh) | 2019-01-22 | 2019-01-22 | 一种空地结合的潮间带一体化测绘方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109631863A (zh) |
WO (1) | WO2020151213A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046211A (zh) * | 2019-04-16 | 2019-07-23 | 重庆市地理信息中心 | 基于语义完整性的测绘成果目录自动检查与入库发布方法 |
CN110516304A (zh) * | 2019-07-26 | 2019-11-29 | 同济大学 | 一种室内空间建模方法 |
CN114018293A (zh) * | 2021-11-19 | 2022-02-08 | 中交第一航务工程局有限公司 | 多波束测深系统的精度检测方法 |
CN114018230A (zh) * | 2021-12-15 | 2022-02-08 | 山东柏远复合材料科技股份有限公司 | 一种空地一体化三维测绘装置 |
CN114739369A (zh) * | 2022-03-09 | 2022-07-12 | 中铁大桥勘测设计院集团有限公司 | 基于无人机航测与测深仪水下测量的滩涂测图方法及设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012032273A (ja) * | 2010-07-30 | 2012-02-16 | Ministry Of Land Infrastructure & Transport Hokkaido Regional Development Bureau | 港湾構造物計測装置 |
CN104050474A (zh) * | 2014-06-10 | 2014-09-17 | 上海海洋大学 | 一种基于LiDAR数据的海岛岸线自动提取方法 |
US20150009060A1 (en) * | 2006-05-16 | 2015-01-08 | Underground Imaging Technology | Sensor Cart Positioning System and Method |
CN105352476A (zh) * | 2015-11-23 | 2016-02-24 | 青岛秀山移动测量有限公司 | 船载水岸线水上水下一体化测量系统集成方法 |
CN106643671A (zh) * | 2016-12-01 | 2017-05-10 | 江苏省测绘工程院 | 一种基于机载LiDAR测深系统的水下点云去噪方法 |
CN108120987A (zh) * | 2017-12-21 | 2018-06-05 | 云南大学 | 一种大江大河的水下河道地形测量装置及测量方法 |
CN108919274A (zh) * | 2018-04-11 | 2018-11-30 | 华南理工大学 | 一种基于单波束的浅水随波扫描探测系统及其工作方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5467122A (en) * | 1991-10-21 | 1995-11-14 | Arete Associates | Underwater imaging in real time, using substantially direct depth-to-display-height lidar streak mapping |
US8903134B2 (en) * | 2010-07-21 | 2014-12-02 | Ron Abileah | Methods for mapping depth and surface current |
CN103711050B (zh) * | 2013-12-31 | 2016-04-06 | 中交第二公路勘察设计研究院有限公司 | 一种激光雷达道路改扩建勘测设计方法 |
-
2019
- 2019-01-22 CN CN201910056969.3A patent/CN109631863A/zh active Pending
- 2019-08-01 WO PCT/CN2019/098770 patent/WO2020151213A1/zh active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150009060A1 (en) * | 2006-05-16 | 2015-01-08 | Underground Imaging Technology | Sensor Cart Positioning System and Method |
JP2012032273A (ja) * | 2010-07-30 | 2012-02-16 | Ministry Of Land Infrastructure & Transport Hokkaido Regional Development Bureau | 港湾構造物計測装置 |
CN104050474A (zh) * | 2014-06-10 | 2014-09-17 | 上海海洋大学 | 一种基于LiDAR数据的海岛岸线自动提取方法 |
CN105352476A (zh) * | 2015-11-23 | 2016-02-24 | 青岛秀山移动测量有限公司 | 船载水岸线水上水下一体化测量系统集成方法 |
CN106643671A (zh) * | 2016-12-01 | 2017-05-10 | 江苏省测绘工程院 | 一种基于机载LiDAR测深系统的水下点云去噪方法 |
CN108120987A (zh) * | 2017-12-21 | 2018-06-05 | 云南大学 | 一种大江大河的水下河道地形测量装置及测量方法 |
CN108919274A (zh) * | 2018-04-11 | 2018-11-30 | 华南理工大学 | 一种基于单波束的浅水随波扫描探测系统及其工作方法 |
Non-Patent Citations (4)
Title |
---|
梁爽等: "工业测量中标准二次曲面的一种拟合方法", 《组合机床与自动化加工技术》 * |
王丽荣: "《机载LiDAR数据误差处理理论与方法》", 31 December 2013 * |
赵夫群: "《点云模型的优化配准方法研究》", 31 July 2018 * |
顾燕: "潮间带地形遥感动态监测体系研究", 《中国博士学位论文全文数据库基础科学辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110046211A (zh) * | 2019-04-16 | 2019-07-23 | 重庆市地理信息中心 | 基于语义完整性的测绘成果目录自动检查与入库发布方法 |
CN110516304A (zh) * | 2019-07-26 | 2019-11-29 | 同济大学 | 一种室内空间建模方法 |
CN110516304B (zh) * | 2019-07-26 | 2023-05-02 | 同济大学 | 一种室内空间建模方法 |
CN114018293A (zh) * | 2021-11-19 | 2022-02-08 | 中交第一航务工程局有限公司 | 多波束测深系统的精度检测方法 |
CN114018293B (zh) * | 2021-11-19 | 2024-03-15 | 中交第一航务工程局有限公司 | 多波束测深系统的精度检测方法 |
CN114018230A (zh) * | 2021-12-15 | 2022-02-08 | 山东柏远复合材料科技股份有限公司 | 一种空地一体化三维测绘装置 |
CN114018230B (zh) * | 2021-12-15 | 2022-07-05 | 山东柏远复合材料科技股份有限公司 | 一种空地一体化三维测绘装置 |
CN114739369A (zh) * | 2022-03-09 | 2022-07-12 | 中铁大桥勘测设计院集团有限公司 | 基于无人机航测与测深仪水下测量的滩涂测图方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
WO2020151213A1 (zh) | 2020-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110779498B (zh) | 基于无人机多视点摄影的浅水河流水深测绘方法及系统 | |
CN109631863A (zh) | 一种空地结合的潮间带一体化测绘方法 | |
CN101551455B (zh) | 干涉合成孔径雷达三维地形成像系统及其高程测绘方法 | |
Li et al. | A comparative study of shoreline mapping techniques | |
CN105352476A (zh) | 船载水岸线水上水下一体化测量系统集成方法 | |
Singh et al. | Microbathymetric mapping from underwater vehicles in the deep ocean | |
CN107918127A (zh) | 一种基于车载InSAR的道路边坡形变检测系统及方法 | |
CN105488852B (zh) | 一种基于地理编码和多维校准的三维图像拼接方法 | |
Dix et al. | Accuracy evaluation of terrestrial LiDAR and multibeam sonar systems mounted on a survey vessel | |
CN109100698B (zh) | 一种用于海上编队的雷达目标球面投影方法 | |
Gruno et al. | Determining sea surface heights using small footprint airborne laser scanning | |
CN111398980A (zh) | 一种机载LiDAR数据处理的方法及装置 | |
Horta et al. | Can recreational echosounder-chartplotter systems be used to perform accurate nearshore bathymetric surveys? | |
Mamajonova et al. | Main Methods of Creating Topographic Maps: A Comprehensive Review | |
CN111487621B (zh) | 一种基于雷达图像的海表流场反演方法及电子设备 | |
Julzarika | Utilization of LAPAN Satellite (TUBSAT, A2, and A3) in supporting Indonesia’s potential as maritime center of the world | |
Spore et al. | Collection, processing, and accuracy of mobile terrestrial lidar survey data in the coastal environment | |
CN116106925A (zh) | 严密的光子计数机制激光雷达水下测深点坐标计算方法 | |
CN116105685A (zh) | 基于声光遥感和滚轮的潮间带地形无缝一体化测量方法 | |
Boeder et al. | Initial experience with the integration of a terrestrial laser scanner into the mobile hydrographic multi sensor system on a ship | |
CN113483739B (zh) | 海上目标位置测量方法 | |
Prempraneerach et al. | Hydrographical survey using point cloud data from laser scanner and echo sounder | |
CN110136487A (zh) | 地图导航防触礁安全助航技术应用模型 | |
CN114234932A (zh) | 一种获取海底控制点数据的水下导线测量方法及装置 | |
Woolard et al. | Shoreline mapping from airborne lidar in Shilshole Bay, Washington |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190416 |