CN103808277A - 一种多传感器点云拼接误差的修正方法 - Google Patents
一种多传感器点云拼接误差的修正方法 Download PDFInfo
- Publication number
- CN103808277A CN103808277A CN201310724678.XA CN201310724678A CN103808277A CN 103808277 A CN103808277 A CN 103808277A CN 201310724678 A CN201310724678 A CN 201310724678A CN 103808277 A CN103808277 A CN 103808277A
- Authority
- CN
- China
- Prior art keywords
- scanning
- actual value
- data
- scan
- layer
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title abstract description 36
- 238000012937 correction Methods 0.000 claims abstract description 38
- 238000003754 machining Methods 0.000 claims abstract description 15
- 238000013519 translation Methods 0.000 claims description 35
- 230000009466 transformation Effects 0.000 claims description 32
- 229910000838 Al alloy Inorganic materials 0.000 claims description 19
- 238000002715 modification method Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 239000000956 alloy Substances 0.000 abstract 1
- 229910045601 alloy Inorganic materials 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 11
- 238000005259 measurement Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 229910052602 gypsum Inorganic materials 0.000 description 3
- 239000010440 gypsum Substances 0.000 description 3
- 239000004411 aluminium Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012113 quantitative test Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
一种多传感器点云拼接误差的修正方法,是采用激光三维人体扫描仪作为扫描工具,采用空心铝合金圆柱体作为标准物体,基于Matlab软件平台实现,具体包括如下步骤:真实值的确立;扫描数据与真实值的圆函数间关系的确立;建立误差修正模型;应用误差修正模型进行补偿。本发明的作为标准物体的空心圆柱体形状简单且能够确保加工精度,圆柱体与扫描高度等高,可对整个扫描区域的多传感器点云拼接误差进行有效校正;修正参数根据每层的全局误差规律得到,与被测物位置无关,可以对整个扫描区域内任何位置的物体的点云拼接误差进行有效修正,所以非常适合用于激光三维人体扫描仪现场。本发明求出的修正值能够有效降低截面拟合尺寸的误差。
Description
技术领域
本发明涉及一种误差的修正方法。特别是涉及一种应用于激光三维扫描中的多传感器点云拼接误差的修正方法。
背景技术
三维数字化技术把现实物体的信息转化为能够直接被计算机处理的数字信息,通过建立各种3D测量理论和算法对这些信息进行处理和分析,最终重构出空间三维物体的数字化模型。随着数字化技术和计算机技术的飞速发展,三维数字化模型在产品设计与制造、虚拟现实、逆向工程、人体工程学等领域具有重要的应用价值和广阔的发展前景。
以三维数字化为核心的激光三维扫描***将被测物转换为点数据的集合,称为点云。随着激光三维扫描技术的快速发展及应用,得到的被测物点云数据的精度问题以及点云的误差分析问题变得尤为重要。激光三维扫描***中,点云误差可能来源于硬件、软件或被测物等多方面,但是很难确定所有误差来源和每种误差所占比例及具体量值,所以应该综合考虑全局误差。估计扫描***中全局误差的适当方法是在充分可控的环境中对不同的标准物体进行测量。如果整个***中的测量误差能够确定,就可以用校正公式对测量结果进行修正。
文献“Establishing calibration models and error envelopes of a3D whole body scanning system”(Lin Y C,Wang M J J,Chen T.Proceedings of the5th WSEAS international conference on Instrumentation,measurement,circuits and systems.World Scientific and Engineering Academy and Society(WSEAS),2006:129-134)开发的标准物体包括一个钢板平台、一个垂直于钢板的标准支架和5个垂直于支架安装的不同尺寸的块规。用德国Vitronic公司的三维人体扫描仪进行扫描,得到不同支架位置的不同块规安装位置的测量数据。最终根据多组扫描数据建立宽度和深度两个方向的畸变校正模型,并给出了误差包络线,对扫描结果进行修正。文献“A Performance evaluation test for laser line scanners on CMMs”(Van Gestel N,Cuypers S,Bleys P,et al.Optics and Lasers in Engineering,2009,47(3):336-342)利用一个平行平板作为标准物体,用夹持在CMMS上的激光扫描器从不同的扫描位置和不同的方向,对标准物体进行扫描。可以从扫描深度、内面角和外面角三个方面分别得到激光扫描器的***误差和随机误差。文献“Error compensation for three-dimensional line laser scanning data”(Xi F,Liu Y,Feng H Y.The International Journal of Advanced Manufacturing Technology,2001,18(3):211-216)把标准球面和标准平板以固定间距组合在一起作为参考物体,固定在平台上,然后用线激光扫描器沿平行于平台的方向进行移动扫描。用多个块规提高平台的一端形成投影角。变换投影角度和扫描深度进行多次采集,然后用最小二乘拟合的方法求出表征误差与扫描深度、投影角度之间关系的经验公式,最后用迭代接近的方法对误差进行补偿。相比而言,国内相关的研究非常少。
上述的第一种方法虽然可以对人体扫描仪的误差从宽度和深度两方面进行修正,但是使用的标准物体复杂,建立的修正算法不能方便的、快速的进行误差修正。而后两种方法虽然标准物体简单,但是只能建立单独一个线激光扫描器的扫描深度、角度与误差的关系,不能 对多个扫描器或传感器之间的误差进行修正。
发明内容
本发明所要解决的技术问题是,提供一种能够有效减小激光三维扫描仪多传感器点云拼接误差影响的多传感器点云拼接误差的修正方法。
本发明所采用的技术方案是:一种多传感器点云拼接误差的修正方法,是采用激光三维人体扫描仪作为扫描工具,采用空心铝合金圆柱体作为标准物体,基于Matlab软件平台实现,具体包括如下步骤:
1)真实值的确立:将空心铝合金圆柱体作为标准物体,放置于激光三维人体扫描仪扫描范围中心区域,进行5次以上的扫描,每一次扫描都得到分层显示的点云数据,根据每一次扫描结果确定作为该次扫描的真实值的圆函数的圆心坐标,并根据空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标求得作为该次扫描的真实值的圆函数,从而得到5次以上的扫描的真实值的圆函数;
2)扫描数据与真实值的圆函数间关系的确立:将每一次扫描得到的扫描数据的每一层向该次扫描的真实值的圆函数进行平移刚性变换,得到该次扫描的每一层x方向和y方向的平移变换参数;
3)建立误差修正模型:将每一次扫描的每一层x方向和y方向的平移变换参数与该层竖直坐标Zw一同保存,构成每一次扫描的整个扫描范围内的误差修正模型,对5次以上扫描的修正模型的相对应层的平移变换参数求取平均值,作为最终误差修正模型中的该层的修正值;
4)应用误差修正模型进行补偿:对每层扫描数据进行如下操作:读取该层扫描数据的竖直坐标Zw,在误差修正模型中寻找最接近的竖直坐标Zw值,读取该竖直坐标Zw对应的修正值,利用修正值对扫描数据进行平移变换,得到修正后的坐标,与该层扫描数据的竖直坐标Zw一并保存。
步骤1)中所述的确定作为真实值的圆函数的圆心坐标,是对每一次扫描得到的点云数据的每一层分别经过拟合求得圆心坐标,再对该次扫描数据求得的所有层的圆心坐标求取平均值,该平均值就是该次扫描的真实值的圆函数的圆心坐标。
步骤1)中所述的根据空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标求得作为该次扫描的真实值的圆函数,是将空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标代入如下公式,求得该次扫描的真实值的圆函数,
x2+y2-2ax-2by+a2+b2-r2=0
式中,x为该层点云的横坐标,y为该层点云的纵坐标,a为圆心坐标的横坐标,b为圆心坐标的纵坐标,r为空心铝合金圆柱体的加工半径。
步骤2)中所述的x方向和y方向的平移变换参数,是利用Matlab的nlinfit函数求出,所求出的结果即为x方向和y方向的平移变换参数。
步骤3)中所述的该层竖直坐标Zw,是在三维人体扫描仪中自上而下逐层扫描且每层数据互相无干扰,在5次以上的扫描中相同层的数据有相同的竖直坐标Zw。
本发明的一种多传感器点云拼接误差的修正方法,设计出的作为标准物体的空心圆柱体 形状简单且能够确保加工精度,圆柱体与扫描高度等高,可对整个扫描区域的多传感器点云拼接误差进行有效校正;修正参数根据每层的全局误差规律得到,与被测物位置无关,可以对整个扫描区域内任何位置的物体的点云拼接误差进行有效修正,所以非常适合用于激光三维人体扫描仪现场;提出平移刚性变换的方法来求扫描数据和真实值之间的修正值,平移变换不会改变像素之间的关系,只是改变点云的位置,所以能够快速、有效的求得扫描数据和真实值之间的变换关系,即修正值。定量分析,本发明求出的修正值能够有效降低截面拟合尺寸的误差;本发明的修正方法可以使任意形状的多传感器点云间的拼接更加光顺、平滑。
具体实施方式
下面结合实施例对本发明的一种多传感器点云拼接误差的修正方法做出详细说明。
本发明的一种多传感器点云拼接误差的修正方法,采用激光三维人体扫描仪作为扫描工具,采用空心铝合金圆柱体作为标准物体,基于Matlab软件平台实现。其中激光三维人体扫描仪可以采用采用专利号为200510013085.8所给出的激光三维人体扫描仪或德国的VitusSmart XXL三维人体扫描仪。
本发明的一种多传感器点云拼接误差的修正方法包括如下步骤:
1)真实值的确立:
估计扫描***中全局误差的一种适当方法是在充分可控的环境中对标准物体进行扫描,建立标准物体与扫描数据之间的误差关系。所以需要通过确立标准物体来确立真实值。标准物体需满足以下要求:由于误差补偿必须在扫描仪调试现场使用,因此标准物体既要容易加工且精度高,又要方便运输和安装调整;根据标准物体建立的补偿方法需要能够对整个扫描范围内位于任何位置物体的扫描结果都可补偿;由于***误差不会随着被测物的形状而发生变化,所以需要标准物体的形状便于拟合和计算。
具体是:为了减小随机误差影响,将空心铝合金圆柱体作为标准物体,放置于激光三维人体扫描仪扫描范围中心区域,进行5次以上的扫描,每一次扫描都得到分层显示的点云数据,根据每一次扫描结果确定作为该次扫描的真实值的圆函数的圆心坐标,所述的确定作为真实值的圆函数的圆心坐标,是对每一次扫描得到的点云数据的每一层分别经过拟合求得圆心坐标,再对该次扫描数据求得的所有层的圆心坐标求取平均值,该平均值就是该次扫描的真实值的圆函数的圆心坐标(a,b)。再并根据空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标求得作为该次扫描的真实值的圆函数,具体是将空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标代入如下公式,求得该次扫描的真实值的圆函数,
x2+y2-2ax-2bu+a2+b2-r2=0
式中,x为该层点云的横坐标,y为该层点云的纵坐标,a为圆心坐标的横坐标,b为圆心坐标的纵坐标,r为空心铝合金圆柱体的加工半径。从而得到5次以上的扫描的真实值的圆函数;
2)扫描数据与真实值的圆函数间关系的确立:
将每一次扫描得到的扫描数据的每一层向该次扫描的真实值的圆函数进行平移刚性变 换,得到该次扫描的每一层x方向和y方向的平移变换参数,所述的x方向和y方向的平移变换参数,是利用Matlab的nlinfit函数求出,所求出的结果即为x方向和y方向的平移变换参数。
常用空间变换中的刚体变换是一种典型的线性变换,即只有物体的位置(平移变换)和朝向(旋转变换)发生改变,而不对像素间的相对空间关系作任何改变。将单个CCD的单层测量结果作为整体向作为真实值的圆函数做刚性变换,可以得到平移和旋转参数。假设旋转的角度为M,x方向平移量为D1,y方向平移量为D2,这样变换后能够最接近作为真实值的圆函数,存在误差的扫描数据用表示,修正后的数据用表示,则
即x2=(cosM×x1-sinM×y1)+D1
y1=(sinM×x1+cosM×y1)+D2
修正后的x2和y2和真值圆函数上的相应的点之间的差值为残差v,其平方和叫做残差平方和,它的值越小,就说明拟合的程度越好。残差平方和目标函数为
其中,n是修正后点的数量,vi为第i个点修正后的残差方程,通过求导可以得到线性方程组,对方程组求解,可以得到系数值M,D1,D2。由此可见刚性变换可以分为只有旋转变换无平移变换、只有平移变换无旋转变换和平移旋转都有的变换。为了具体确定采取哪种刚性变换方法,选取圆柱体的一层截面测量数据,分别用这三种方法相对于作为真实值的圆函数做变换,求出系数M,D1,D2以及变换后的残差平方和如表1所示。表中可见,既有旋转又有平移的情况和只有平移的情况的残差平方和完全相同,然而只有旋转的残差平方和非常大,约等于另外两种的32.952倍,所以排除只有旋转的方法。在圆柱体和长方棱柱体的扫描数据中选取多层数据分别用另外两种方法修正,对修正后的效果进行对比和选择。只平移变换的方法对于两种形状的截面数据都有很好的修正效果,且根据实验可知修正时间比旋转加平移的方法短。本发明选择平移的方法进行每一次扫描得到的扫描数据的每一层向该次扫描的真实值的圆函数的变换。
表1不同刚性变换类型求出的修正值及残差平方和
3)建立误差修正模型:
将每一次扫描的每一层x方向和y方向的平移变换参数与该层竖直坐标Zw一同保存,构成每一次扫描的整个扫描范围内的误差修正模型,对5次以上扫描的修正模型的相对应层的平移变换参数求取平均值,作为最终误差修正模型中的该层的修正值。所述的该层竖直坐标Zw,是在三维人体扫描仪中自上而下逐层扫描且每层数据互相无干扰,在5次以上的扫描中相同层的数据有相同的竖直坐标Zw。
即,求出每一次扫描数据的每一层与该次扫描的真实值的圆函数的x方向的平移变换参数D1,y方向的平移变换参数D2。然后,与该层竖直坐标Zw一同保存,构成每一次扫描的整个扫描范围内的误差修正模型,对5次以上扫描的修正模型的相对应层的平移变换参数求取平均值,作为最终误差修正模型中的该层的修正值。
4)应用误差修正模型进行补偿:
对每层扫描数据进行如下操作:读取该层扫描数据的竖直坐标Zw,在误差修正模型中寻找最接近的竖直坐标Zw值,读取该竖直坐标Zw对应的修正值,利用修正值对扫描数据进行平移变换,得到修正后的坐标,与该层扫描数据的竖直坐标Zw一并保存。逐层对扫描数据进行修正,得到修正后的整体点云。
本发明求出的修正值能够有效降低截面拟合尺寸的误差,如表2、4所示,圆柱体和长方棱柱体截面修正前后拟合直径的相对误差对比列表;直观观察,本发明的修正方法可以使任意形状的多传感器点云间的拼接更加光顺、平滑。
下面给出具体实例:
本实施例所采用的激光三维人体扫描仪是采用专利号为200510013085.8中所给出的激光三维人体扫描仪。其包括4根立柱,每根立柱上安装一组传感器,每组传感器由一个线激光光源和上下两个对称分布的CCD构成。四个激光光源位于同一水平面,在人体表面形成一个环带。在同步扫描***作用下,4组传感器从上到下按一定扫描间隔分层扫描。扫描范围为直径1000mm、高度2000mm,水平分辨率为2mm,深度分辨率为2mm,垂直分辨率为4mm,扫描时间为16.7秒。采用本发明的修正算法,首先建立误差修正模型,然后分别对铝合金圆柱体、长方棱柱体以及石膏人体模特进行修正实验,并给出实验结果。
(1)误差修正模型的建立。本发明所设计的标准物体-铝合金空心圆柱体,其高度为2000mm,横截面圆直径设计值为176.5mm,实际测量值为176.29mm至176.38mm,在整个高度范围内误差为-0.21mm至-0.12mm。把圆柱体放置在扫描区域的中间位置s1处,进行10次重复扫描。把每次扫描得到的扫描数据的每一层进行圆拟合求得圆心坐标,将该次扫描的所有层圆心坐标求平均得到平均圆心坐标,加之圆柱体加工半径为截面真实值的圆半径,求得该次扫描的作为真实值的圆函数,然后把每层扫描数据向作为真实值的圆函数进行平移刚性变换,求得每层的平移变换参数,与该层的竖直坐标一起保存,建立起该次扫描的整个扫描范围内的误差修正模型。对10次重复扫描的修正模型的相对应层的平移变换参数求取平均值,作为最终误差修正模型中的该层的修正值。
(2)不同位置圆柱体扫描数据修正实验及结果。为了验证误差修正模型是否对不同位置的圆柱体扫描数据具有好的修正效果,把标准圆柱体放置在扫描区域的s2和s3位置进行扫描,并用求得的修正值对扫描数据进行补偿修正。为了对修正前后点云拼接效果进行比较, 对s1、s2、s3三个位置的扫描结果,从第一层开始,Zw每隔200mm进行一次抽样,得到9个截面的数据,用matlab软件对每个截面数据进行外凸壳法非线性最小二乘圆拟合求得直径,结果如表2所示,其中D为直径,RE为相对误差。
表2三个位置修正前后圆拟合直径对比列表
表中数据看出,s1位置补偿前后的相对误差平均值由3.4%降低为1.3%;s2位置补偿前后的相对误差平均值由2.9%降为1.3%;s3位置补偿前后的相对误差平均值由2.9%降为1.1%。虽然补偿后截面直径的误差有所减小,但与实际值相比仍有2mm左右的误差,且误差总是正误差,即测量值总是大于实际值。分析原因,一是本发明的方法是利用平移刚性变换法将每层扫描数据向真值圆函数进行逼近,会有残留误差;二是截面直径的测量是先采用外凸壳法测得截面周长,再由周长求得直径,外凸壳法会导致截面直径相对于真实值偏大。可以采用内凸壳法和外凸壳法分别测量截面点云的外部轮廓和内部轮廓直径,再求取两个测量值的平均值作为最后测量值,会更接近实际值,表3给出表2中s1位置处修正后截面直径分别用内凸壳法和外凸壳法求得的测量值及其平均值。可以看出,平均后的结果误差较小,更加接近真实值。
表3s1位置修正后内、外凸壳法拟合直径及二者平均值
以上实验结果说明本发明建立的误差修正模型对不同位置的圆柱体扫描误差的补偿是有 效的,其补偿效果与圆柱摆放的位置无关。
(3)不同位置长方棱柱体扫描数据修正实验及结果。为了验证平均修正值对长方棱柱体的扫描数据的补偿效果,把截面边长理论值为133mm×85mm,周长436mm,高2000mm的长方棱柱铝型材放置在扫描区域中的s1、s2和s3位置进行扫描,并对扫描结果进行修正。经实际测量,整个高度内铝型材截面长度的误差为0.22mm至0.30mm,宽度的误差为0.15mm至0.40mm,截面周长的误差为0.74mm至1.40mm。如对长方棱柱铝型材扫描结果选取三个截面修正前后外凸壳法求得的截面周长结果如表4所示,其中P为周长,RE为相对误差。从表中数据看出,s1位置补偿前后的相对误差平均值由3.1%降低为1.4%;s2位置补偿前后的相对误差平均值由2.9%降为1.4%;s3位置补偿前后的相对误差平均值由2.8%降为1.5%。结果与圆柱体测量结果类似,也存在一定的正误差,原因也与上面分析一样,对s1位置修正后的截面周长,同时采用外凸壳法和内凸壳法测量的截面尺寸及平均尺寸结果如表5所示,可以看出,平均尺寸更接近真实值。
表4三个位置棱柱修正前后外凸壳拟合周长对比列表
表5s1位置修正后内、外凸壳法拟合周长及二者平均值
说明本发明建立的误差修正模型对不同位置的长方棱柱体扫描误差的补偿是有效的,且补偿效果与长方棱柱体的摆放位置无关。
(4)石膏人体模特扫描数据修正实验及结果。为了验证本发明建立的误差修正模型是否对不规则形状点云的多传感器拼接误差仍然有好的修正效果,对放置在扫描区域中s1位置的石膏人体模特进行扫描,并对扫描结果进行修正。可看出,修正前,由于误差的影响,上下传感器的点云没有重合在一起,相邻传感器点云没有很好的拼接在一起,导致点云轮廓线产生一定的宽度。修正后,点云轮廓线的宽度便细,说明上下传感器的重合程度、相邻传感器的拼接效果都变好。
综上所述,本发明提出的基于圆柱标准物体的多传感器点云拼接误差修正算法简单、有效,非常适合用于激光三维人体扫描仪现场使用。
Claims (5)
1.一种多传感器点云拼接误差的修正方法,其特征在于,采用激光三维人体扫描仪作为扫描工具,采用空心铝合金圆柱体作为标准物体,基于Matlab软件平台实现,具体包括如下步骤:
1)真实值的确立:将空心铝合金圆柱体作为标准物体,放置于激光三维人体扫描仪扫描范围中心区域,进行5次以上的扫描,每一次扫描都得到分层显示的点云数据,根据每一次扫描结果确定作为该次扫描的真实值的圆函数的圆心坐标,并根据空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标求得作为该次扫描的真实值的圆函数,从而得到5次以上的扫描的真实值的圆函数;
2)扫描数据与真实值的圆函数间关系的确立:将每一次扫描得到的扫描数据的每一层向该次扫描的真实值的圆函数进行平移刚性变换,得到该次扫描的每一层x方向和y方向的平移变换参数;
3)建立误差修正模型:将每一次扫描的每一层x方向和y方向的平移变换参数与该层竖直坐标Zw一同保存,构成每一次扫描的整个扫描范围内的误差修正模型,对5次以上扫描的修正模型的相对应层的平移变换参数求取平均值,作为最终误差修正模型中的该层的修正值;
4)应用误差修正模型进行补偿:对每层扫描数据进行如下操作:读取该层扫描数据的竖直坐标Zw,在误差修正模型中寻找最接近的竖直坐标Zw值,读取该竖直坐标Zw对应的修正值,利用修正值对扫描数据进行平移变换,得到修正后的坐标,与该层扫描数据的竖直坐标Zw一并保存。
2.根据权利要求1所述的一种多传感器点云拼接误差的修正方法,其特征在于,步骤1)中所述的确定作为真实值的圆函数的圆心坐标,是对每一次扫描得到的点云数据的每一层分别经过拟合求得圆心坐标,再对该次扫描数据求得的所有层的圆心坐标求取平均值,该平均值就是该次扫描的真实值的圆函数的圆心坐标。
3.根据权利要求1所述的一种多传感器点云拼接误差的修正方法,其特征在于,步骤1)中所述的根据空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标求得作为该次扫描的真实值的圆函数,是将空心铝合金圆柱体的加工半径和该次扫描的真实值的圆函数的圆心坐标代入如下公式,求得该次扫描的真实值的圆函数,
x2+y2-2ax-2by+a2+b2-r2=0
式中,x为该层点云的横坐标,y为该层点云的纵坐标,a为圆心坐标的横坐标,b为圆心坐标的纵坐标,r为空心铝合金圆柱体的加工半径。
4.根据权利要求1所述的一种多传感器点云拼接误差的修正方法,其特征在于,步骤2)中所述的x方向和y方向的平移变换参数,是利用Matlab的nlinfit函数求出,所求出的结果即为x方向和y方向的平移变换参数。
5.根据权利要求1所述的一种多传感器点云拼接误差的修正方法,其特征在于,步骤3)中所述的该层竖直坐标Zw,是在三维人体扫描仪中自上而下逐层扫描且每层数据互相无干扰,在5次以上的扫描中相同层的数据有相同的竖直坐标Zw。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310724678.XA CN103808277B (zh) | 2013-12-23 | 2013-12-23 | 一种多传感器点云拼接误差的修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310724678.XA CN103808277B (zh) | 2013-12-23 | 2013-12-23 | 一种多传感器点云拼接误差的修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103808277A true CN103808277A (zh) | 2014-05-21 |
CN103808277B CN103808277B (zh) | 2016-07-06 |
Family
ID=50705392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310724678.XA Expired - Fee Related CN103808277B (zh) | 2013-12-23 | 2013-12-23 | 一种多传感器点云拼接误差的修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103808277B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104111054A (zh) * | 2014-07-08 | 2014-10-22 | 哈尔滨工业大学 | 一种基于回归分析的反射面精度分析方法 |
CN104296659A (zh) * | 2014-10-27 | 2015-01-21 | 苏州天准精密技术有限公司 | 一种利用倾斜标准面修正拼接误差方法 |
CN105115426A (zh) * | 2015-09-06 | 2015-12-02 | 合肥芯硕半导体有限公司 | 一种激光直接成像设备图形拼接误差的检测方法 |
CN107644121A (zh) * | 2017-08-18 | 2018-01-30 | 昆明理工大学 | 一种路面材料骨架结构的逆向三维重构与实体建模方法 |
CN107784138A (zh) * | 2016-08-26 | 2018-03-09 | 上海宝冶集团有限公司 | 基于结构力学分析的点云重力变形修正方法 |
CN108180834A (zh) * | 2018-02-05 | 2018-06-19 | 中铁二十二局集团有限公司 | 一种工业机器人同三维成像仪位姿关系现场实时标定方法 |
CN108663029A (zh) * | 2018-02-01 | 2018-10-16 | 深圳市建设综合勘察设计院有限公司 | 一种获取水下圆柱形基桩信息的方法、存储介质及终端 |
CN108871268A (zh) * | 2018-07-13 | 2018-11-23 | 湖南联智桥隧技术有限公司 | 一种基于激光点云的隧道超欠挖数值计算方法 |
CN110111381A (zh) * | 2019-03-13 | 2019-08-09 | 中山易裁剪网络科技有限公司 | 一种远程确定衣服尺寸***及其确定方法 |
CN114729849A (zh) * | 2019-11-19 | 2022-07-08 | 惠普发展公司,有限责任合伙企业 | 确定扫描仪的优选区域 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4801207A (en) * | 1985-06-14 | 1989-01-31 | The Broken Hill Proprietary Company Limited | Method and apparatus for the optical determination of surface profiles |
JPH1074272A (ja) * | 1996-08-30 | 1998-03-17 | Minolta Co Ltd | 3次元形状データ処理装置 |
JP2007121126A (ja) * | 2005-10-28 | 2007-05-17 | Yokohama Rubber Co Ltd:The | 3次元形状測定装置及び3次元形状測定方法 |
JP3945279B2 (ja) * | 2002-03-15 | 2007-07-18 | ソニー株式会社 | 障害物認識装置、障害物認識方法、及び障害物認識プログラム並びに移動型ロボット装置 |
JP2012018076A (ja) * | 2010-07-08 | 2012-01-26 | Jeol Ltd | 走査プローブ顕微鏡を用いた試料の磁気力測定方法及び装置 |
WO2012018076A1 (ja) * | 2010-08-04 | 2012-02-09 | 株式会社ブリヂストン | タイヤ外形計測データの補正方法、及びタイヤ外観検査装置 |
CN102494657A (zh) * | 2011-12-12 | 2012-06-13 | 北京建筑工程学院 | 一种曲面轮廓测量及检测的测头半径补偿方法 |
CN103106632A (zh) * | 2012-11-29 | 2013-05-15 | 华中科技大学 | 一种基于均值漂移的不同精度三维点云数据的融合方法 |
CN103453849A (zh) * | 2013-07-18 | 2013-12-18 | 黑龙江科技大学 | 多光学传感器协同的复杂曲面零件三维测量方法与*** |
-
2013
- 2013-12-23 CN CN201310724678.XA patent/CN103808277B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4801207A (en) * | 1985-06-14 | 1989-01-31 | The Broken Hill Proprietary Company Limited | Method and apparatus for the optical determination of surface profiles |
JPH1074272A (ja) * | 1996-08-30 | 1998-03-17 | Minolta Co Ltd | 3次元形状データ処理装置 |
JP3945279B2 (ja) * | 2002-03-15 | 2007-07-18 | ソニー株式会社 | 障害物認識装置、障害物認識方法、及び障害物認識プログラム並びに移動型ロボット装置 |
JP2007121126A (ja) * | 2005-10-28 | 2007-05-17 | Yokohama Rubber Co Ltd:The | 3次元形状測定装置及び3次元形状測定方法 |
JP2012018076A (ja) * | 2010-07-08 | 2012-01-26 | Jeol Ltd | 走査プローブ顕微鏡を用いた試料の磁気力測定方法及び装置 |
WO2012018076A1 (ja) * | 2010-08-04 | 2012-02-09 | 株式会社ブリヂストン | タイヤ外形計測データの補正方法、及びタイヤ外観検査装置 |
CN102494657A (zh) * | 2011-12-12 | 2012-06-13 | 北京建筑工程学院 | 一种曲面轮廓测量及检测的测头半径补偿方法 |
CN103106632A (zh) * | 2012-11-29 | 2013-05-15 | 华中科技大学 | 一种基于均值漂移的不同精度三维点云数据的融合方法 |
CN103453849A (zh) * | 2013-07-18 | 2013-12-18 | 黑龙江科技大学 | 多光学传感器协同的复杂曲面零件三维测量方法与*** |
Non-Patent Citations (3)
Title |
---|
F. XI ET AL: "Error Compensation for Three-Dimensional Line Laser Scanning", 《TNT J ADV MAMIF TECHNOL》 * |
左超等: "一种新的点云拼接算法", 《中国激光》 * |
葛宝臻等: "激光三维彩色人体扫描的数据重建技术", 《激光与光电子学进展》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104111054A (zh) * | 2014-07-08 | 2014-10-22 | 哈尔滨工业大学 | 一种基于回归分析的反射面精度分析方法 |
CN104296659A (zh) * | 2014-10-27 | 2015-01-21 | 苏州天准精密技术有限公司 | 一种利用倾斜标准面修正拼接误差方法 |
CN104296659B (zh) * | 2014-10-27 | 2017-01-18 | 苏州天准科技股份有限公司 | 一种利用倾斜标准面修正拼接误差方法 |
CN105115426A (zh) * | 2015-09-06 | 2015-12-02 | 合肥芯硕半导体有限公司 | 一种激光直接成像设备图形拼接误差的检测方法 |
CN107784138A (zh) * | 2016-08-26 | 2018-03-09 | 上海宝冶集团有限公司 | 基于结构力学分析的点云重力变形修正方法 |
CN107784138B (zh) * | 2016-08-26 | 2021-03-26 | 上海宝冶集团有限公司 | 基于结构力学分析的点云重力变形修正方法 |
CN107644121B (zh) * | 2017-08-18 | 2020-07-31 | 昆明理工大学 | 一种路面材料骨架结构的逆向三维重构与实体建模方法 |
CN107644121A (zh) * | 2017-08-18 | 2018-01-30 | 昆明理工大学 | 一种路面材料骨架结构的逆向三维重构与实体建模方法 |
CN108663029A (zh) * | 2018-02-01 | 2018-10-16 | 深圳市建设综合勘察设计院有限公司 | 一种获取水下圆柱形基桩信息的方法、存储介质及终端 |
CN108180834A (zh) * | 2018-02-05 | 2018-06-19 | 中铁二十二局集团有限公司 | 一种工业机器人同三维成像仪位姿关系现场实时标定方法 |
CN108871268B (zh) * | 2018-07-13 | 2021-02-02 | 湖南联智科技股份有限公司 | 一种基于激光点云的隧道超欠挖数值计算方法 |
CN108871268A (zh) * | 2018-07-13 | 2018-11-23 | 湖南联智桥隧技术有限公司 | 一种基于激光点云的隧道超欠挖数值计算方法 |
CN110111381A (zh) * | 2019-03-13 | 2019-08-09 | 中山易裁剪网络科技有限公司 | 一种远程确定衣服尺寸***及其确定方法 |
CN114729849A (zh) * | 2019-11-19 | 2022-07-08 | 惠普发展公司,有限责任合伙企业 | 确定扫描仪的优选区域 |
Also Published As
Publication number | Publication date |
---|---|
CN103808277B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103808277A (zh) | 一种多传感器点云拼接误差的修正方法 | |
CN102003938B (zh) | 大型高温锻件热态在位检测方法 | |
Teza et al. | Geometric characterization of a cylinder-shaped structure from laser scanner data: Development of an analysis tool and its use on a leaning bell tower | |
CN103292748A (zh) | 一种基于激光测量的多基板拼合检测方法 | |
CN103026310B (zh) | 实现从机床加工点到安装测量基准点的空间转换方法 | |
CN108198224A (zh) | 一种用于立体视觉测量的线阵相机标定装置及标定方法 | |
CN106017319A (zh) | 一种基于高精度点位测量的三维激光扫描数据坐标转换工具及方法 | |
CN109990703A (zh) | 一种预制构件的尺寸检测方法及*** | |
CN109631754B (zh) | 一种测量装置坐标系标定的方法及相关装置 | |
CN103486984A (zh) | 一种风洞内型面同轴度的检测方法 | |
CN104729529B (zh) | 地形图测量***误差判断的方法和*** | |
Wiedemann et al. | Reconstruction of historical buildings based on images from the Meydenbauer archives | |
CN112857315A (zh) | 基于三维激光扫描的钢管柱垂直度测量方法 | |
KR101373139B1 (ko) | 레이저 간섭계를 이용한 직각도 측정 방법 | |
JP2010151757A (ja) | トータルステーションを用いた測量方法 | |
CN208061260U (zh) | 一种用于立体视觉测量的线阵相机标定装置 | |
CN108801187A (zh) | 基于坐标变换的导轨滑台运动的几何误差辨识方法 | |
Welkenhuyzen et al. | Investigation of the kinematic system of a 450 kV CT scanner and its influence on dimensional CT metrology applications | |
Zhang et al. | Large-scale shape measurement by a combined method based on three instruments | |
CN114485438A (zh) | 一种大型模块钢结构圆立柱间距测量方法 | |
UEHAN et al. | Development of an aerial survey system and numerical analysis modeling software for unstable rock blocks | |
Mill | Simulation of terrestrial laser scanning errors occurring during deformation monitoring | |
Kang et al. | Computing method for estimating strain and stress of steel beams using terrestrial laser scanning and FEM | |
KR101255901B1 (ko) | 시공 중인 구조물의 수직 변위량 측정 장치 및 시공 중인 구조물의 기둥 축소량에 대한 보정값 산출 장치 | |
Zhang et al. | Novel three-dimensional data conversion technique and profile measurement system for engine cylinder head blank |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 |