首页 / 专利库 / 工业自动化和数控机床 / 坐标系 / 扫描仪坐标系与工程测量坐标系的变换方法

扫描仪坐标系与工程测量坐标系的变换方法

阅读:276发布:2020-05-11

专利汇可以提供扫描仪坐标系与工程测量坐标系的变换方法专利检索,专利查询,专利分析的服务。并且本 发明 提供一种 扫描仪 坐标系 与工程测量坐标系的变换方法,步骤是:1)定向标靶中心的工程测量坐标为柔性坐标,其性质是坐标反算的两点之间的长度与对应地表长度成非线性关系,具有柔性性质;地面激光扫描仪获得的点 云 坐标为刚性坐标,其性质是坐标反算的两点之间平面距离与对应地表长度是相等,具有刚性特征;2)将工程测量坐标即柔性坐标系坐标通过消除距离归算和投影的综合 变形 ,转换成 正交 的三维直 角 坐标系坐标,与三维激光扫描坐标系坐标达到一致性;3)建立刚性化点云定向模型,进行定向参数求解;4)扫描站所有点云的工程测量坐标转换计算。本发明能最大限度消除工程测量中归算变形和投影变形等引起的系统误差。,下面是扫描仪坐标系与工程测量坐标系的变换方法专利的具体信息内容。

1.一种扫描仪坐标系与工程测量坐标系的变换方法,其特征在于采用以下步骤:
1)定向标靶中心的工程测量坐标为柔性坐标,其性质是坐标反算的两点之间的长度与对应地表长度成非线性关系,具有柔性性质,具体是:①坐标反算的平距离s,随着到中央子午线的不同、地形起伏变化呈非线性变化,规律是

其中Δs表示综合变形,Δs1是将地面观测的长度归算至参考椭球面产生的归算变形,是高程对长度归算的影响,Hm为定向标靶的平均大地高;Δs2是投影变形,ym表示点到中央子午线的距离,Δy表示定向标靶与扫描站坐标差,Rm表示扫描站地球曲率半径;
②Z坐标受大地水准曲率和大气垂直折光影响,规律是

其中C为球气差系数,Rm表示扫描站地球曲率半径,s为坐标反算的水平距离,K为扫描站处的大气折光系数,K值为0.08-0.14之间;而地面激光扫描仪获得的点坐标为刚性坐标,其性质是坐标反算的两点之间平面距离与对应地表长度是相等,具有刚性特征;
2)柔性坐标系的刚性化:将工程测量坐标即柔性坐标系坐标通过消除距离归算和投影的综合变形,转换成正交的三维直坐标系坐标,与三维激光扫描坐标系坐标达到一致性,减小测量数据所带来的系统误差,具体是:利用扫描站在工程测量坐标系的概略坐标,对定向标靶中心相对于扫描站的工程测量坐标进行归算和投影改正;利用三角高程测量原理,对正常高进行改正,减小似大地水准面曲率影响和大气折光的影响,改正后刚性坐标为


Z'=Z-Cs2

式中 和 为扫描站的概略坐标,X和Y为高斯坐标,Δy为扫描站到定向标靶的Y坐标差,Z为正常高;λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子,Hm为定向标靶的平均大地高,如果没有高程异常,平均大地高可用正常高代替,ym是定向标靶到中央子午线的平均值距离,Rm是扫描站处平均曲率半径,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离,大气垂直折光系数K值在0.08-0.14之间;
3)建立刚性化点云定向模型,简称为柔性模型,进行定向参数求解,柔性模型具体为:

用3个以上定向标靶中心的扫描坐标(x,y,z)和工程测量坐标系坐标(X,Y,Z),平差计算出6个定向参数,公式中参数X'、Y'、Z'、 λ、C和s均已在步骤2)中定义过,具体是:
X'、Y'、Z'为改正后刚性坐标, 为扫描站的概略坐标,λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离,XS、YS、ZS是扫描仪扫描时的位置参数,R是扫描时扫描仪在工程测量坐标系的姿态参数构成的旋转矩阵;
4)扫描站所有点云的工程测量坐标转换计算,首先按柔性模型计算本站所有扫描点云的刚性化坐标 然后按公式 计算点云的工程测量坐标,公式中
参数λ、 C和s均已在步骤2)中定义过,具体是:λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子, 为扫描站的概略坐标,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离。

说明书全文

扫描仪坐标系与工程测量坐标系的变换方法

技术领域

[0001] 本发明涉及一种扫描仪坐标系与工程测量坐标系的变换方法,属于工程测量领域。

背景技术

[0002] 工程对象如大型工程建设设施、一个区域的地形等,是放在工程测量坐标系中进行研究的,为了控制投影变形和归算变形选择,工程坐标系要选择中央子午线和抵偿面,使投影后综合变形小于2.5cm。当用三维激光扫描进行工程测量时,选用高斯投影东坐标为X坐标、北坐标为Y坐标、以正常高为Z坐标,组成3维右手直坐标系,用3个以上控制点的扫描坐标(x,y,z)和工程测量坐标系坐标(X,Y,Z),用(1)式或(2)模型式计算两个坐标系的参数。
[0003]
[0004]
[0005] XS、YS、ZS是扫描仪扫描时的位置参数;R是扫描时扫描仪在工程测量坐标系的姿态参数构成的旋转矩阵。根据工程测量“坐标反算长度与表面实测长度相差小于指定的值”的要求(1)的(2)式有两个不足之处:①公式(1)给三维激光扫描坐标z乘λ是不合适的,缩放因子一般由长度变形引起,而高程没有此变形,所以选用(2)式比较好。②工程测量坐标(X,Y,Z)中的高斯坐标受综合长度变形残余影响,是有系统误差的,在用2个式子所列的误差方程式只对偶然误差进行平差时,存在模型误差,势必影响坐标转换精度,所以在高精度测量中,用(2)式列方程时,需要消除(X,Y,Z)中的系统误差。

发明内容

[0006] 本发明的目的是提供一种能解决上述问题、转换精度高的扫描仪坐标系与工程测量坐标系的变换方法。其技术方案为:
[0007] 一种扫描仪坐标系与工程测量坐标系的变换方法,其特征在于采用以下步骤:
[0008] 1)定向标靶中心的工程测量坐标为柔性坐标,其性质是坐标反算的两点之间的长度与对应地表长度成非线性关系,具有柔性性质,具体是:①坐标反算的平距离s,随着到中央子午线的不同、地形起伏变化呈非线性变化,规律是
[0009]
[0010] 其中Δs表示综合变形,Δs1是将地面观测的长度归算至参考椭球面产生的归算变形,是高程对长度归算的影响,Hm为定向标靶的平均大地高;Δs2是投影变形,ym表示点到中央子午线的距离,Δy表示定向标靶与扫描站坐标差,Rm表示扫描站地球曲率半径;
[0011] ②Z坐标受大地水准曲率和大气垂直折光影响,规律是
[0012]
[0013] 其中C为球气差系数,Rm表示扫描站地球曲率半径,s为坐标反算的水平距离,K为扫描站处的大气折光系数,K值为0.08-0.14之间;而地面激光扫描仪获得的点坐标为刚性坐标,其性质是坐标反算的两点之间平面距离与对应地表长度是相等,具有刚性特征;
[0014] 2)柔性坐标系的刚性化:将工程测量坐标即柔性坐标系坐标通过消除距离归算和投影的综合变形,转换成正交的三维直角坐标系坐标,与三维激光扫描坐标系坐标达到一致性,减小测量数据所带来的系统误差,具体是:利用扫描站在工程测量坐标系的概略坐标,对定向标靶中心相对于扫描站的工程测量坐标进行归算和投影改正;利用三角高程测量原理,对正常高进行改正,减小似大地水准面曲率影响和大气折光的影响,改正后刚性坐标为
[0015]
[0016]
[0017] Z'=Z-Cs2
[0018]
[0019] 式中 和 为扫描站的概略坐标,X和Y为高斯坐标,Δy为扫描站到定向标靶的Y坐标差,Z为正常高;λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子,Hm为定向标靶的平均大地高,如果没有高程异常,平均大地高可用正常高代替,ym是定向标靶到中央子午线的平均值距离,Rm是扫描站处平均曲率半径,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离,大气垂直折光系数K值在0.08-0.14之间;
[0020] 3)建立刚性化点云定向模型,简称为柔性模型,进行定向参数求解,柔性模型具体为:
[0021]
[0022] 用3个以上定向标靶中心的扫描坐标(x,y,z)和工程测量坐标系坐标(X,Y,Z),平差计算出6个定向参数,公式中参数X'、Y'、Z'、 λ、C和s均已在步骤2)中定义过,具体是:X'、Y'、Z'为改正后刚性坐标, 为扫描站的概略坐标,λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离,XS、YS、ZS是扫描仪扫描时的位置参数,R是扫描时扫描仪在工程测量坐标系的姿态参数构成的旋转矩阵;
[0023] 4)扫描站所有点云的工程测量坐标转换计算,首先按柔性模型计算本站所有扫描点云的刚性化坐标 然后按公式 计算点云的工程测量坐标,公式中参数λ、 C和s均已在步骤2)中定义过,具体是:λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子, 为扫描站的概略坐标,C为大气折光误差和似大地水准面弯曲误差的球气差系数,s为定向标靶到扫描站的水平距离。
[0024] 本发明与现有技术相比,其优点在于:1.解决现有模型存在的问题,最大限度地消除工程测量中归算变形和投影变形等引起的系统误差的影响;2.针对工程测量坐标有变形误差影响的问题,提出坐标转换的柔性模型,以用于每公里综合变形大于2.5cm的长距离三维激光扫描高精度测量,为地面三维激光扫描仪应用于高精度工程测量提供技术支持。

具体实施方式

[0025] 实验所用扫描仪器为Riegl VZ-1000,使用球形定向标靶,共扫描4站,每站设5个球形标靶,其中相邻扫描站最少2个重合球形标靶,每站约4000万个扫描点。球形标靶中心扫描坐标由三维激光数据处理软件自动测算出。标靶对中地面控制点的工程测量坐标用RTK技术测量标靶地面点平面坐标和高程,通过量取标靶高,将工程测量的高程引入到标靶中心。工程测量坐标系为80西安坐标系,高斯投影的中央子午线为117度,扫描区离中央子午线1°40',高程系统选用85高程基准。用刚性模型公式和柔性模型公式计算的标靶转换后各项误差见表1。实现刚性坐标与柔性坐标变换的具体步骤是:
[0026] 步骤1)定向标靶中心的工程测量坐标为柔性坐标,而地面激光扫描仪获得的点云坐标为刚性坐标,柔性坐标系相对于刚性坐标系的性质有:
[0027] ①坐标反算的水平距离s,随着到中央子午线的不同、地形起伏变化呈非线性变化,规律是
[0028]
[0029] 其中Δs表示综合变形,Δs1是将地面观测的长度归算至参考椭球面产生的归算变形,是高程对长度归算的影响,Hm为控制标靶的平均大地高;Δs2是投影变形,ym表示点到中央子午线的距离,Δy表示标靶与扫描站坐标差,Rm表示扫描站地球曲率半径;
[0030] ②Z坐标受大地水准曲率和大气垂直折光影响,规律是
[0031]
[0032] 其中C为球气差系数,K为扫描站处的大气折光系数,K值为0.08-0.14之间。
[0033] 步骤2)柔性坐标系的刚性化:利用扫描站在工程测量坐标系的概略坐标,对定向标靶中心相对于扫描站的工程测量坐标,进行投影改正和投影改正;利用三角高程测量原理,对正常高进行改正,减小似大地水准面曲率影响和大气折光的影响,坐标改正公式为[0034]
[0035]
[0036] Z'=Z-Cs2
[0037]
[0038] 式中 和 为扫描站的概略坐标,X和Y为高斯坐标,Δy为扫描站到定向标靶的Y坐标差,Z为正常高;λ扫描站处的长度归算变形和投影变形和的系数即综合变形因子,Hm为控制标靶的平均大地高,如果没有高程异常,平均大地高可用正常高代替,ym是标靶到中央子午线的平均值距离,Rm是扫描站处平均曲率半径,C为大气折光误差和似大地水准面弯曲误差的球气差系数,S为定向标靶到扫描站的水平距离;大气垂直折光系数K值在0.08-0.14之间;
[0039] 步骤3)中坐标改正刚性化定向参数求解的数学模型即柔性模型为[0040]
[0041] 用3个以上标靶中心的扫描坐标(x,y,z)和工程测量坐标系坐标(X,Y,Z),平差计算出6个定向参数。
[0042] 步骤4)依据步骤3)确定的定向参数解,按柔性模型计算本站所有扫描点云的刚性化坐标 然后按公式 计算点云的工程测量坐标。
[0043] 对定向后扫描点位置精度进行评定。精度指标有X方向中误差、Y方向中误差、高程中误差、平面点位中误差和空间点位中误差。使用现有刚性模型和本发明的柔性模型计算的5个精度指标如下:
[0044]中误差名称 刚性模型(mm) 柔性模型(mm)
X方向中误差 43.2 10.3
Y方向中误差 35.0 8.2
高程中误差 9.9 9.6
平面点位中误差 55.6 13.2
空间点位中误差 56.5 16.5
[0045] 实验数据证明:1.解决现有模型存在的问题,最大限度地消除工程测量坐标中归算变形和投影变形等引起的系统误差的影响;2.柔性模型有利于提高点云坐标转换后的精度,为地面三维激光扫描仪应用于精密工程测量提供技术支持。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈