CN105627948B - 一种大型复杂曲面测量系统进行复杂曲面采样的方法 - Google Patents
一种大型复杂曲面测量系统进行复杂曲面采样的方法 Download PDFInfo
- Publication number
- CN105627948B CN105627948B CN201610070408.5A CN201610070408A CN105627948B CN 105627948 B CN105627948 B CN 105627948B CN 201610070408 A CN201610070408 A CN 201610070408A CN 105627948 B CN105627948 B CN 105627948B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- point
- laser
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000005070 sampling Methods 0.000 title claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 47
- 230000008569 process Effects 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 238000000691 measurement method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
Abstract
本发明涉及一种大型复杂曲面测量系统及其应用，该测量系统包括底座、水平支架、俯仰支架及计算机；水平支架设置在底座上，俯仰支架与水平支架连接，在俯仰支架上设置有相机、激光测距仪和线激光发射器，在激光测距仪的两侧设置两台所述的相机，计算机分别与相机、激光测距仪、线激光发射器电连接。利用该测量系统通过点激光定位、线激光采样、相机获取图像数据的方式进行，通过摄影方式获取图像，再利用摄影测量原理计算像点的实际坐标，最后通过数据拼接在同一坐标系下进行复杂曲面的数据重建，该测量系统及采样方法在保证精度的情况下，简化了测量过程，减少了测量次数，大大提高了测量效率。
Description
技术领域
本发明涉及一种大型复杂曲面测量系统及其应用，尤其是涉及一种摄影测量技术与激光技术相结合获取大型复杂曲面规则三维点阵的测量系统，以及该测量系统的使用方法。
背景技术
随着全球经济的发展，市场竞争日趋激烈，新的技术革命不断取得进展和突破，技术的飞跃发展已经成为推动世界经济增长的重要因素。为了保持和加强产品在市场上的竞争力，产品的开发周期、生产周期越来越短，促使工业产品越来越向多品种、小批量、高质量、低成本的方向发展，具有复杂曲面的产品越来越多，广泛应用于模具、工具、能源、交通、航空航天、航海等领域。复杂曲面的复杂性主要体现在：其一，许多边缘学科、高科技产品领域对产品涉及的曲面造型有很高的精度要求，以达到某些数学特征的高精度为目的；其二，现代社会中，人们在注重产品功能的同时，对产品的外观造型提出了越来越高的要求，以追求美学效果或功能要求为目的。因此对曲面重建的研究也越来越多，逆向工程以先进产品设备的实物等作为研究对象，应用现代设计方法学原理、生产工程学、材料学和有关专业知识进行系统地分析和研究、探索掌握其关键技术，进而开发出同类的更为先进的产品。
逆向工程的关键技术主要包括三个方面，即曲面数字化测量技术、曲面重构技术和曲面加工技术。数字化测量技术就是利用各种数字化测量手段获取原有实物(产品原型或油泥模型)的三维数字模型。逆向工程中复杂曲面实体的测量方法主要有接触式测量法(如三坐标测量法)和非接触测量法(如激光三角形法、距离法、断层图象扫描法、摄影测量法等)两种。
三坐标测量法利用三坐标测量机(CMM)的探测触头(各种不同直径和形状的探针)逐点地捕捉样件的表面数据。该方法对于没有复杂内部型腔、特征几何尺寸多、只有少量特征曲面的零件该测量方法非常有效。但由于该方法是接触式测量，易于损伤探头和划伤被测样件表面，不能对软质材料和超薄型物体进行测量，对细微部分测量精度也受到影响，应用范围受到限制；始终需要人工干预，不可能实现全自动测量；由于测头的半径而存在三维补偿问题；且价格较高，对使用环境有一定要求；测量速度慢，效率低。
激光三角法根据光学三角形原理，将激光投射到被测物体表面，采用点阵布置的光电敏感元件在另一位置接受激光的成像位置，根据光点或光条在物体上成像的偏移，通过被测物体基平面、像点、像距之间的关系计算物体的深度信息。其测量速度快，可直接测量不可接触的高精度工件、软工件等，但是测量精度取决于感光设备的分辨率和被测表面的光学特性。
距离法是利用光线飞行的时间来计算距离，常采用激光和脉冲光束。可以测量大型物体，但是其测速慢、工作量大，测量相当困难。
断层扫描图像法通过对被测物体进行断层界面扫描，以X射线的衰减系数为依据，重构断层截面图像，根据不同位置的断层图像可建立物体的三维信息。该方法可以对被测物体内部的结构和形状进行无损测量，但其造价高，测量系统的分辨率低，获取数据时间长，设备体积大，对运行的环境要求高。
相比之下，摄影测量法根据同一个三维空间点在不同空间位置的两个(或多个)摄像机拍摄的图像中的视差，以及摄像机之间的空间几何关系来获取该点的三维坐标值。摄影测量法可以对处于两个(或多个)摄像机共同视野内的目标特征点进行测量，而无须伺服机构等扫描装置，是一种非破坏性测量，适于测量大型复杂曲面物体。但其存在的最大问题是空间特征点在多幅数字图像中提取与匹配的精度与准确性等问题。此外按照该种方法获得的图像数据，经过处理提取后得到的点云数据达数兆或者数十兆，点之间没有拓扑关系，为了满足曲面重建的需要，这些点阵数据一般需要经过预处理，主要包括：多视拼合、去除噪音点、数据插补、孔洞修补、数据简化、平滑和光顺等。经过预处理的点阵数据，方可用于重构曲面模型。由于点阵数据预处理过程十分繁杂，稳定性较差，难以保证重构曲面的精度与效率要求。
发明内容
针对现有技术的不足，本发明提供一种大型复杂曲面测量系统。
本发明还提供一种利用大型复杂曲面测量系统进行复杂曲面采样的方法。
本发明的技术方案如下：
一种大型复杂曲面测量系统，包括底座、水平支架、俯仰支架及计算机；所述水平支架设置在底座上，俯仰支架通过转轴与水平支架连接，在俯仰支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧设置两台所述的相机，所述计算机分别与相机、激光测距仪、线激光发射器电连接。
优选的，所述水平支架的底部设置有水平分度盘。此设计的好处在于，当旋转水平支架时，通过水平分度盘能够清楚直观地观察到水平支架旋转的角度。
优选的，所述转轴的一端连接有俯仰分度盘。此设计的好处在于，一方面可通过俯仰分度盘来拧动旋转俯仰支架，用于调节俯仰角度，另一方面可通过俯仰分度盘上的刻度值来准确观察调整俯仰角度。
优选的，两台所述的相机对称设置在激光测距仪的两侧。
一种利用上述的大型复杂曲面测量系统进行复杂曲面采样的方法，包括以下步骤，
一、单一测量区域的数据采样及解算：
(1)将被测物体外表面划分为n个测量区域，每一测量区域为一独立测站，测量系统放置在某一测量区域i的正前方，调整测量系统的位置，使测量系统中相机拍摄的边界位置覆盖该测量区域i；
(2)在测量区域i的左上角确定测量的位置起始点及在测量区域i的右下角确定测量的位置终止点，启动测量系统并通过计算机程序进行控制，然后按如下步骤进行测量：
①激光测距仪发射点激光到测量区域i，调节水平支架和俯仰支架使点激光垂直于测量区域i的表面，此时点激光在测量区域i表面的垂足作为该测站的基准位置点，以测量系统在该位置时的系统空间坐标系作为该测站的基准坐标系；
②调节水平支架和俯仰支架使点激光射到位置起始点并获取该位置时测量系统的位姿坐标；再调节水平支架和俯仰支架使点激光射到位置终止点并获取该位置时测量系统的位姿坐标；
③扫描测量：设置单站测量时测量系统的水平步距角和俯仰步距角，从位置起始点开始根据设置的水平步距角和俯仰步距角，按照“正向从左向右，反向从右向左”的顺序逐行扫描被测物外表面，扫描过程中，对于被测物的任意扫描位置点Qj，激光测距仪发射点激光到Qj，记录相应的激光测距距离Lj，同时记录水平支架和俯仰支架分别相对于基准位置点转过的角度φj和ωj；其中一台相机拍照，获得点激光照片Gj，关闭激光测距仪，打开线激光发射器发射线激光，两台相机同时拍照获得两张线激光照片LAj和 LBj；然后关闭线激光发射器，打开激光测距仪，水平支架或俯仰支架转动一个步距角使点激光对准下一扫描位置点Qj+1位置进行扫描测量；依次对测量区域i的每一扫描位置点进行扫描测量，到位置终止点结束，获得相应扫描位置点的点激光照片与线激光照片；
(3)对测量区域i扫描测量完成后，解算测量区域i的三维点阵数据：
①初始化两台相机和激光测距仪点激光中心的元素，其中每台相机都有3个内方位元素和6个外方位元素：内方位元素包括像主点相对于影像中心的位置x0、y0以及镜头中心到影像面的焦距f，外方位元素包括3个用于描述镜头中心相对于系统空间坐标系位置的线元素XS，YS，ZS和3个用于描述影像面在摄影瞬间空中姿态的角元素ω、k，激光测距仪点激光中心在系统空间坐标系中的坐标为(0,0，Zcj)；
②读取点激光照片Gj中位置点Qj对应的图像坐标(x1j，y1j)；根据(x1j，y1j)在线激光照片LAj中确定其对应的图像坐标点①，以点①为首采样点P1，并将其作为基准点，分别沿着y轴正负方向每间隔△y像素进行采样，得到一组采样点Pk(k＝2,3,4,…)，则Pk的y坐标为ykj＝y1j±(k-1)×△y，其x坐标按照最大像素值原则确定；
③取另一相机线激光照片LBj上一组采样点的y坐标与LAj中一组采样点的y坐标相同，按照最大像素值原则确定其x坐标，最终确定该台相机线激光照片LBj的一组采样点Ph(h＝1,2,3,4，…)的图像坐标为(xhj，yhj)；
(4)按如下公式(Ⅰ)，利用两台相机的3个角元素分别计算其各自的转换矩阵R：
再根据各采样点的图像坐标(x，y)，按照公式(Ⅱ)反算其在当前测量位置系统空间坐标系下的三维坐标(X,Y,Z),可分别求得一组采样点Pk(k＝1,2,3,4,…)和Ph (h＝1,2,3,4,…)所对应的空间点列Sk(k＝1,2,3,4,…)和Sh(h＝1,2,3,4,…)；
对空间点列Sk(k＝1,2,3,4,…)和Sh(h＝1,2,3,4,…)上的两个相同位置处的采样点的坐标取平均值，得到该扫描位置点Qj所获取的一组采样点在系统空间坐标系下的三维坐标Ss(s＝1,2,3,4,…)；
(5)利用公式(Ⅲ)将步骤(4)中解算出来的Ss(s＝1,2,3,4,…)的三维坐标转换到该测站的基准坐标系下；
Sj＝C+RjSs (III)
式中，Sj为一组采样点Ss在该测站基准坐标系下的坐标矩阵，Sj＝[Xj，Yj，Zj]T； Rj为变换矩阵R，由扫描位置点Qj对应的角度(φj，ωj，0)代入公式(Ⅰ)求得； Ss为一组采样点Ss的坐标向量，C为常数矩阵，C＝[0,0,Zcj]T；
(6)按步骤(2)-(5)所述方法完成各扫描位置点Qj(j＝1,2,3，…)的图像数据采样及解算，最终得到测量区域i在该测站基准坐标系下的三维点阵；
二、相邻测量区域的数据拼接：
(7)在两个相邻测量区域重叠的公共区域上布置9个拼接控制点；
(8)利用测量系统按照“一、单一测量区域的数据采样及解算”中所述的方法对测量区域i进行测量，得到该测量区域i在该测站基准坐标系下的三维点阵Su (u＝1,2,3,4,…)；保持测量系统在该测站位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点，并利用两台相机获取拼接控制点的图像数据，然后按照“一、单一测量区域的数据采样及解算”中所述的方法计算出9个拼接控制点在该测站基准坐标系下的对应点J1i(XJ1i,YJ1i,ZJ1i)(i＝1,2,3,4,…,9)；
(9)将测量系统移到相邻测量区域i+1的正前方，首先按照“一、单一测量区域的数据采样及解算”中所述的方法获取该测量区域i+1在该测站基准坐标系下的三维点阵 Sd(d＝1,2,3,4,…)，保持测量系统在该测站位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点，并利用两台相机获取拼接控制点的图像数据，然后按照“一、单一测量区域的数据采样及解算”中所述的方法计算出9个拼接控制点在该测站基准坐标系下的对应点J2i(XJ2i,YJ2i,ZJ2i)(i＝1,2,3,4,…,9)；
(10)分别计算拼接控制点J1i(i＝1,2,3,4,…,9)和J2i(i＝1,2,3,4,…,9) 的重心坐标G1(XG1,YG1,ZG1)和G2(XG2,YG2,ZG2)；分别以重心点G1和G2为原点，建立两重心坐标系，两重心坐标系的坐标轴方向与各自的基准坐标系一致，则两重心坐标系之间的转换矩阵R′可由式(Ⅳ)求得：
式中，
(11)利用公式(Ⅴ)将测量区域i的三维点阵Su(u＝1,2,3,4,…)转换到测量区域i+1的基准坐标系下，得到S′u(u＝1,2,3,4,…)，
S′u＝)Su-G1)R′+G2 (V)
式中，S′u为三维点阵Su转换到测量区域i+1的基准坐标系下的坐标矩阵，S′u＝[X′u， Y′u，Z′u]T；Su为三维点阵Su在测量区域i的基准坐标系下的坐标矩阵，Su＝[Xu，Yu，Zu]T； G1和G2分别为重心点G1和G2的坐标矢量，G1＝[XG1,YG1,ZG1]T，G2＝[XG2,YG2,ZG2]T；
(12)按照步骤(7)-(11)所述的方法，依次将各测量区域基准坐标系下的坐标转换到下一测量区域基准坐标系下，即完成各测量区域点阵数据的拼接。
优选的，步骤(3)中，80≤△y≤120。
本发明的有益效果在于：
(1)本发明通过研制用于大型复杂曲面的新型测量系统以及利用该测量系统所实现的测量方法，相比于目前采用的激光扫描技术对被测物整体扫描，本发明的测量方法采用激光和摄影测量相结合的方法，对于具有大型复杂曲面的测量不再需要对被测物持续循环扫描，而是通过点激光定位、线激光采样、相机获取图像数据的方式进行，这样通过摄影方式获取图像，再利用摄影测量原理计算像点的实际坐标，在保证精度的情况下，简化了测量过程，减少了测量次数，大大提高了测量效率。
(2)目前利用扫描或者摄影测量获得的被测物点阵数据是不规则的散乱点，而散乱点阵的曲面重建仍是一大难题，须经过测量数据的压缩、分割、修补和规则化等繁杂的处理，且难以保证重建曲面的精度与效率。本发明的测量方法利用测量系统对被测物外表面按照设定的步距角与采样间距进行逐行扫描采样，获得的采样点是具有确定行列关系的点阵数据，便于后期的数据处理、管理与存储。此外，本发明的测量方法可通过调整步距角与采样间距，有效调整采样密度，保证了测量的灵活性与适应性。
(3)本发明的测量系统配有两台相机获取激光的图像数据，避免因单台相机的系统误差而导致精度降低。相比于目前单相机摄影测量，两台相机的配合工作提高了测量精度，最终获得的数据更可靠。
附图说明
图1为本发明实现获取大型复杂曲面三维点阵的流程图；
图2为本发明测量系统的结构示意图；
图3为本发明测量系统对被测物进行单站测量时的原理图；
图4a为解算其中一台相机线激光图像数据方法示意图；
图4b为解算另一台相机线激光图像数据方法示意图；
图5为相邻两测站解算数据进行拼接的方法示意图。
其中：1、底座；2、水平分度盘；3、水平支架；4、俯仰支架；5、相机；6、激光测距仪；7、线激光发射器；8、相机；9、俯仰分度盘；10、数据线接口；11、数据线； 12、计算机；13、测量系统；14、被测物；15、位置起始点；16、垂足位置点；17、位置终止点；18、点激光；19、测站；20、区域Ⅰ；21、公共区域；22、拼接控制点；23、区域Ⅱ；24、测站。
具体实施方式
下面通过实施例并结合附图对本发明做进一步说明，但不限于此。
实施例1：
如图1和图2所示，一种大型复杂曲面测量系统，该系统是基于双向空间前方交汇原理实现获取大型复杂曲面三维点阵的装置。该测量系统包括底座1、水平分度盘2、水平支架3、俯仰支架4、相机5和相机8、激光测距仪6、线激光发射器7、俯仰分度盘9、数据线接口10、数据线11和计算机12；水平分度盘2与水平支架3固定在一起，二者安装在底座1的上方；相机5和相机8、以及激光测距仪6和线激光发射器7安装在俯仰支架4上；俯仰支架4安装在水平支架3上；俯仰分度盘9固定安装在俯仰支架4上，在水平支架3的右侧；负责程序控制与数据采集的计算机12通过数据线11与位于水平支架3上的数据线接口10连接。
底座1、水平分度盘2和水平支架3在竖直方向共轴线安装，水平分度盘2与水平支架3用螺栓固定连接，水平分度盘2可绕自身轴线与底座1之间实现相对转动，从而实现系统转动位置水平分度。
水平支架3左右两侧各有一个立板，立板上有水平方向的通孔，用来安装俯仰支架4，俯仰支架4左右两端各有一段短转轴，连接到水平支架3立板的通孔上,，来支撑俯仰支架4；激光测距仪6和线激光发射器7布置在俯仰支架4的中心线上，相机5和相机8 对称安装在激光测距仪6的两侧。
俯仰分度盘9与俯仰支架4右端的短转轴固定连接，通过旋转俯仰分度盘9来控制俯仰支架4的俯仰角。
水平支架3内部装有两个步进电机，用来驱动水平分度盘2和俯仰分度盘9，实现系统的水平转动和俯仰。
计算机12上装有系统的驱动控制程序，用来驱动步进电机转动、相机5、相机8、激光测距仪6和线激光发射器7的协调动作，以及采集相机得到的图像数据和系统的位姿坐标。计算机12上还装有后处理程序，用来处理图像数据，和实现解算三维点阵的方法。
本实施例中测量系统的底座1通过连接件与目前测绘领域通用的测量三脚架连接，以固定和支撑测量系统。
该测量系统自带有一个系统空间坐标系，相机5和相机8以及激光测距仪6在该坐标系下的坐标已知，通过初始化相机得到坐标数据。
实施例2：
基于实施例1所述的测量系统，本实施例详细阐述利用实施例1所述的测量系统进行大型复杂曲面的测量方法，方法如下：
一、单一测量区域的数据采样及解算(对目标物体的外表面进行单站测量与解算的过程)：
如附图3所示，将被测物体外表面14分为n个区域，测量系统13正对着某测量区域i放置，激光测距仪发射的点激光18近似垂直于被测物体外表面14，垂足为位置点 16；位置点15和位置点17分别为测量起始点和终止点。单站测量是以设定的步距角，利用激光循环扫描、相机捕捉获取图像数据进行的。激光扫描过程中，是以点激光作为采样基准点，辅以线激光进行扫描测量，在图像处理过程中，在线激光上按照指定采样间距取点，结合两台相机的数据，利用双向空间前方交会原理解算得到规则的三维点阵。具体步骤如下：
(1)将被测物外表面14分为n个区域，每一区域为一独立测站。为了后续便于不同区域采样数据的配准与拼接，相邻区域间存在重叠的公共区域。对于任一被测区域i，调整测量系统13，使其正对区域i放置；
(2)通过驱动控制程序，打开激光测距仪6发射点激光18到被测区域i的位置点16，调节水平支架3水平转动和俯仰支架4俯仰转动，使点激光近似垂直于被测物外表面14，将此时的垂足位置点16作为该站测量的基准点，以测量系统13在该位置时的系统系统空间坐标系作为测量区域i的基准坐标系，这时设置驱动程序复位，水平支架3 和俯仰支架4分别相对于基准垂足位置点16转过的角度φj和ωj设为零；
(3)调节水平支架3和俯仰支架4转动，使点激光打到被测区域i的左上角位置起始点15，获取该位置时测量系统的位姿坐标(位置坐标也即测量系统的线元素、姿态坐标也即测量系统的角元素，简称位姿坐标)；然后调节水平支架和俯仰支架使点激光达到被测区域i的右下角位置终止点17，获取该位置时测量系统的位姿坐标；
(4)设置扫描测量时的水平步距角和俯仰步距角；
(5)扫描测量：启动激光扫描测量，根据附图3中虚线所示的路径，从位置起始点15开始，按照“正向从左向右，反向从右向左”的顺序逐行扫描被测物外表面。扫描过程中，对于被测物的任意扫描位置点Qj，激光测距仪6发射点激光到Qj，记录相应的激光测距距离Lj，同时记录水平支架3和俯仰支架4分别相对于基准垂足位置点16转过的角度φj和ωj；相机5或相机8拍照，获得点激光照片Gj；关闭激光测距仪6，打开线激光发射器7发射线激光，相机5和相机8同时拍照获得两张线激光照片LAj和LBj；然后关闭线激光发射器7，打开激光测距仪6，转动水平支架3和俯仰支架4使点激光对准下一扫描位置点Qj+1位置进行扫描测量，水平方向上的位置点扫描完成后转动俯仰支架到下一行反向进行位置点的扫描。按上述方法对被测区域i的每一扫描位置点进行扫描测量，到位置终止点17结束，获得相应扫描位置点的点激光照片与线激光照片；
(6)测量完成后，解算被测区域i的三维点阵数据。初始化两台相机5、8及激光测距仪6的元素，其中每台相机都有3个内方位元素和6个外方位元素。内方位元素包括像主点相对于影像中心的位置x0、y0以及镜头中心到影像面的焦距f(也称主距)，外方位元素包括3个用于描述镜头中心相对于系统空间坐标系位置的线元素XS，YS，ZS和 3个用于描述影像面在摄影瞬间空中姿态的角元素φ、ω、k。激光测距仪点激光中心系统空间坐标系中的坐标为(0,0，Zcj)；
(7)打开后处理程序，下面的步骤针对某一扫描位置点Qj说明解算点阵的过程；
(8)读取点激光照片Gj中激光点Qj对应的像坐标(x1j，y1j)；
(9)处理线激光图像数据解算点阵。如附图4a、4b所示，根据点激光照片Gj中激光点Qj的图像坐标(x1j，y1j)，在线激光照片LAj中确定其对应的图像坐标点①，以点①为首采样点P1，并将其作为基准点，分别沿着首采样点P1的y轴正负方向每间隔△y(80≤△y≤120)像素进行取样，得到一组采样点Pk(k＝2,3,4,…)，则Pk的y坐标为ykj＝y1j± (k-1)×△y，其x坐标按照最大像素值原则确定；
(10)取另一相机线激光照片LBj上所有采样点的y坐标与上述LAj中一组采样点的y坐标相同，按照最大像素值原则确定其x坐标，确定该台相机线激光照片LBj的采样点 Ph(h＝1,2,3,4，…)的像坐标(xhj，yhj)；
(11)按照如下公式(Ⅰ)，利用两台相机的3个角元素分别计算其转换矩阵R：
(12)根据各采样点的图像坐标(x，y)，按照公式(Ⅱ)反算其在当前测量位置系统空间坐标系下的三维坐标(X,Y,Z),可分别求得一组采样点Pk(k＝1,2,3,4,…)和 Ph(h＝1,2,3,4,…)所对应的空间点列Sk(k＝1,2,3,4,…)和Sh(h＝1,2,3,4,…)；
(13)对空间点列Sk(k＝1,2,3,4,…)和Sh(h＝1,2,3,4,…)各对应点的坐标取平均值，得到该扫描位置点Qj所获取的一组采样点在其系统空间坐标系下的三维坐标Ss(s＝1,2,3,4,…)；
(14)利用公式(Ⅲ)将步骤(13)中解算出来的Ss(s＝1,2,3,4,…)的三维坐标转换到该测站的基准坐标系下；
Sj＝C+RjSs (III)
式中，Sj为一组采样点Ss在该测站基准坐标系下的坐标矩阵，Sj＝[Xj，Yj，Zj]T， Rj为变换矩阵R，由扫描位置点Qj对应的角度(φj，ωj，0)代入公式(Ⅰ)求得， Ss为一组采样点Ss的坐标向量，C为常数矩阵，C＝[0,0,Zcj]T；
(15)按上述方法完成各扫描位置点Qj(j＝1,2,3，…)的图像数据采样及解算，最终得到测量区域i在该测站基准坐标系下的三维点阵；
二、相邻测量区域的数据拼接：
基于在第一大步中得到被测物外表面各区域的三维点阵，接下来详细介绍各区域数据的拼接方法，以将各区域的数据整合到同一坐标系下。如附图5所示，拼接是通过添加公共控制点实现的。区域Ⅰ20和区域Ⅱ23之间有重叠的公共区域21,在公共区域21上布置9个拼接控制点22。将测量系统先后放置于测量区域Ⅰ20的正前方测站19和测量区域Ⅱ23的正前方测站24进行测量。具体的实施方法如下步骤所述：
(16)在公共区域21上布置9个拼接控制点22；
(17)按照第一大步中介绍的方法，利用测量系统在附图5所示的测站19对被测物14的区域Ⅰ20进行测量，得到区域Ⅰ20在该测站19基准坐标系下的三维点阵Su (u＝1,2,3,4,…)；保持测量系统在测站19位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点22，并利用两台相机5、8获取拼接控制点的图像数据，然后按照第一大步中介绍的步骤计算出9个拼接控制点22在该测站19基准坐标系下的对应点J1i (XJ1i,YJ1i,ZJ1i)(i＝1,2,3,4,…,9)；
(18)将测量系统移到附图5所示的测站24的位置，首先按照第一大步中介绍的方法获取区域Ⅱ23在该测站24基准坐标系下的三维点阵Sd(d＝1,2,3,4,…)，保持测量系统在测站24位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点22，并利用两台相机5、8获取拼接控制点22的图像数据，然后按照第一大步中介绍的方法计算出9个拼接控制点22在该测站24基准坐标系下的对应点J2i(XJ2i,YJ2i,ZJ2i) (i＝1,2,3,4,…,9)；
(19)分别计算控制点J1i(i＝1,2,3,4,…,9)和J2i(i＝1,2,3,4,…,9)的重心坐标G1(XG1, YG1,ZG1)和G2(XG2,YG2,ZG2)；分别以重心点G1和G2为原点，建立两重心坐标系，两重心坐标系的坐标轴方向与各自的基准坐标系一致，则两重心坐标系之间的转换矩阵R′可由式(Ⅳ)求得：
式中，
(20)利用公式(Ⅴ)将测量区域i的三维点阵Su(u＝1,2,3,4,…)转换到测量区域i+1的基准坐标系下，得到S′u(u＝1,2,3,4,…)，
S′u＝(Su-G1)R′+G2 (V)
式中，S′u为三维点阵Su转换到测量区域i+1的基准坐标系下的坐标矩阵，S′u＝[X′u， Y′u，Z′u]T，Su为三维点阵Su在测量区域i的基准坐标系下的坐标矩阵，Su＝[Xu，Yu，Zu]T， G1和G2分别为重心点G1和G2的坐标矢量，G1＝[XG1,YG1,ZG1]T，G2＝[XG2,YG2,ZG2]T；
(21)按照步骤(16)-(20)所述的方法，依次将各测量区域基准坐标系下的坐标转换到下一测量区域基准坐标系下，即完成各测量区域点阵数据的拼接。
Claims (5)
1.一种大型复杂曲面测量系统进行复杂曲面采样的方法，所述大型复杂曲面测量系统包括底座、水平支架、俯仰支架及计算机；所述水平支架设置在底座上，俯仰支架通过转轴与水平支架连接，在俯仰支架上设置有相机、激光测距仪和线激光发射器，激光测距仪与线激光发射器位于同一垂线上，在激光测距仪的两侧设置两台所述的相机，所述计算机分别与相机、激光测距仪、线激光发射器电连接，其特征在于，包括以下步骤，
一、单一测量区域的数据采样及解算：
(1)将被测物体外表面划分为n个测量区域，每一测量区域为一独立测站，测量系统放置在某一测量区域i的正前方，调整测量系统的位置，使测量系统中相机拍摄的边界位置覆盖该测量区域i；
(2)在测量区域i的左上角确定测量的位置起始点及在测量区域i的右下角确定测量的位置终止点，启动测量系统并通过计算机程序进行控制，然后按如下步骤进行测量：
①激光测距仪发射点激光到测量区域i，调节水平支架和俯仰支架使点激光垂直于测量区域i的表面，此时点激光在测量区域i表面的垂足作为该测站的基准位置点，以测量系统在该位置时的系统空间坐标系作为该测站的基准坐标系；
②调节水平支架和俯仰支架使点激光射到位置起始点并获取该位置时测量系统的位姿坐标；再调节水平支架和俯仰支架使点激光射到位置终止点并获取该位置时测量系统的位姿坐标；
③扫描测量：设置单站测量时测量系统的水平步距角和俯仰步距角，从位置起始点开始根据设置的水平步距角和俯仰步距角，按照“正向从左向右，反向从右向左”的顺序逐行扫描被测物外表面，扫描过程中，对于被测物的任意扫描位置点Qj，激光测距仪发射点激光到Qj，记录相应的激光测距距离Lj，同时记录水平支架和俯仰支架分别相对于基准位置点转过的角度φj和ωj；其中一台相机拍照，获得点激光照片Gj，关闭激光测距仪，打开线激光发射器发射线激光，两台相机同时拍照获得两张线激光照片LAj和LBj；然后关闭线激光发射器，打开激光测距仪，水平支架或俯仰支架转动一个步距角使点激光对准下一扫描位置点Qj+1位置进行扫描测量；依次对测量区域i的每一扫描位置点进行扫描测量，到位置终止点结束，获得相应扫描位置点的点激光照片与线激光照片；
(3)对测量区域i扫描测量完成后，解算测量区域i的三维点阵数据：
①初始化两台相机和激光测距仪点激光中心的元素，其中每台相机都有3个内方位元素和6个外方位元素：内方位元素包括像主点相对于影像中心的位置x0、y0以及镜头中心到影像面的焦距f，外方位元素包括3个用于描述镜头中心相对于系统空间坐标系位置的线元素XS，YS，ZS和3个用于描述影像面在摄影瞬间空中姿态的角元素ω、k，激光测距仪点激光中心在系统空间坐标系中的坐标为(0,0，Zcj)；
②读取点激光照片Gj中位置点Qj对应的图像坐标(x1j，y1j)；根据(x1j，y1j)在线激光照片LAj中确定其对应的图像坐标点①，以点①为首采样点P1，并将其作为基准点，分别沿着y轴正负方向每间隔△y像素进行采样，得到一组采样点Pk，k＝2,3,4,…，则Pk的y坐标为ykj＝y1j±(k-1)×△y，其x坐标按照最大像素值原则确定；
③取另一相机线激光照片LBj上一组采样点的y坐标与LAj中一组采样点的y坐标相同，按照最大像素值原则确定其x坐标，最终确定该台相机线激光照片LBj的一组采样点Ph，h＝1,2,3,4，…的图像坐标为(xhj，yhj)；
(4)按如下公式(Ⅰ)，利用两台相机的3个角元素分别计算其各自的转换矩阵R：
再根据各采样点的图像坐标(x，y)，按照公式(Ⅱ)反算其在当前测量位置系统空间坐标系下的三维坐标(X,Y,Z),可分别求得一组采样点Pk，k＝1,2,3,4,…和Ph，h＝1,2,3,4,…所对应的空间点列Sk，k＝1,2,3,4,…和Sh，h＝1,2,3,4,…；
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>f</mi>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>f</mi>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mo>=</mo>
<msub>
<mi>Zc</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mi>L</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
对空间点列Sk，k＝1,2,3,4,…和Sh，h＝1,2,3,4,…上的两个相同位置处的采样点的坐标取平均值，得到该扫描位置点Qj所获取的一组采样点在系统空间坐标系下的三维坐标Ss，s＝1,2,3,4,…；
(5)利用公式(Ⅲ)将步骤(4)中解算出来的Ss，s＝1,2,3,4,…的三维坐标转换到该测站的基准坐标系下；
Sj＝C+RjSs (III)
式中，Sj为一组采样点Ss在该测站基准坐标系下的坐标矩阵，Sj＝[Xj，Yj，Zj]T；Rj为变换矩阵R，由扫描位置点Qj对应的角度(φj，ωj，0)代入公式(Ⅰ)求得；Ss为一组采样点Ss的坐标向量，C为常数矩阵，C＝[0,0,Zcj]T；
(6)按步骤(2)-(5)所述方法完成各扫描位置点Qj，j＝1,2,3，…的图像数据采样及解算，最终得到测量区域i在该测站基准坐标系下的三维点阵；
二、相邻测量区域的数据拼接：
(7)在两个相邻测量区域重叠的公共区域上布置9个拼接控制点；
(8)利用测量系统按照“一、单一测量区域的数据采样及解算”中所述的方法对测量区域i进行测量，得到该测量区域i在该测站基准坐标系下的三维点阵Su，u＝1,2,3,4,…；保持测量系统在该测站位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点，并利用两台相机获取拼接控制点的图像数据，然后按照“一、单一测量区域的数据采样及解算”中所述的方法计算出9个拼接控制点在该测站基准坐标系下的对应点J1i(XJ1i,YJ1i,ZJ1i)，i＝1,2,3,4,…,9；
(9)将测量系统移到相邻测量区域i+1的正前方，首先按照“一、单一测量区域的数据采样及解算”中所述的方法获取该测量区域i+1在该测站基准坐标系下的三维点阵Sd，d＝1,2,3,4,…，保持测量系统在该测站位置不变，调整测量系统使其发射的点激光依次对准9个拼接控制点，并利用两台相机获取拼接控制点的图像数据，然后按照“一、单一测量区域的数据采样及解算”中所述的方法计算出9个拼接控制点在该测站基准坐标系下的对应点J2i(XJ2i,YJ2i,ZJ2i)，i＝1,2,3,4,…,9；
(10)分别计算拼接控制点J1i，i＝1,2,3,4,…,9和J2i，i＝1,2,3,4,…,9的重心坐标G1(XG1,YG1,ZG1)和G2(XG2,YG2,ZG2)；分别以重心点G1和G2为原点，建立两重心坐标系，两重心坐标系的坐标轴方向与各自的基准坐标系一致，则两重心坐标系之间的转换矩阵R′可由式(Ⅳ)求得：
<mrow>
<msup>
<mi>R</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<mi>b</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>b</mi>
</mrow>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>V</mi>
<mo>)</mo>
</mrow>
</mrow>
式中，
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>9</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>9</mn>
</munderover>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>9</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>9</mn>
</munderover>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>J</mi>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>G</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>J</mi>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>G</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
(11)利用公式(Ⅴ)将测量区域i的三维点阵Su，u＝1,2,3,4,…转换到测量区域i+1的基准坐标系下，得到S′u，u＝1,2,3,4,…，
S′u＝(Su-G1)R′+G2 (V)
式中，S′u为三维点阵Su转换到测量区域i+1的基准坐标系下的坐标矩阵，S′u＝[X′u，Y′u，Z′u]T；Su为三维点阵Su在测量区域i的基准坐标系下的坐标矩阵，Su＝[Xu，Yu，Zu]T；G1和G2分别为重心点G1和G2的坐标矢量，G1＝[XG1,YG1,ZG1]T，G2＝[XG2,YG2,ZG2]T；
(12)按照步骤(7)-(11)所述的方法，依次将各测量区域基准坐标系下的坐标转换到下一测量区域基准坐标系下，即完成各测量区域点阵数据的拼接。
2.如权利要求1所述的复杂曲面采样的方法，其特征在于，所述水平支架的底部设置有水平分度盘。
3.如权利要求1所述的复杂曲面采样的方法，其特征在于，所述转轴的一端连接有俯仰分度盘。
4.如权利要求1所述的复杂曲面采样的方法，其特征在于，两台所述的相机对称设置在激光测距仪的两侧。
5.如权利要求1所述的复杂曲面采样的方法，其特征在于，步骤(3)中，80≤△y≤120。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610070408.5A CN105627948B (zh) | 2016-01-31 | 2016-01-31 | 一种大型复杂曲面测量系统进行复杂曲面采样的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610070408.5A CN105627948B (zh) | 2016-01-31 | 2016-01-31 | 一种大型复杂曲面测量系统进行复杂曲面采样的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105627948A CN105627948A (zh) | 2016-06-01 |
| CN105627948B true CN105627948B (zh) | 2018-02-06 |
Family
ID=56043101
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610070408.5A Expired - Fee Related CN105627948B (zh) | 2016-01-31 | 2016-01-31 | 一种大型复杂曲面测量系统进行复杂曲面采样的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105627948B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106225671B (zh) * | 2016-08-03 | 2019-01-18 | 湖南航天机电设备与特种材料研究所 | 一种大型浮空器囊体体积测量装置的现场校准方法 |
| CN108050957B (zh) * | 2018-01-13 | 2020-01-07 | 山东科技大学 | 一种波形唇油封唇口曲面形状误差测量装置及测量方法 |
| CN109360267B (zh) * | 2018-09-29 | 2023-06-06 | 杭州蓝芯科技有限公司 | 一种薄物体快速三维重建方法 |
| CN109297433A (zh) * | 2018-11-15 | 2019-02-01 | 青岛星晖昌达智能自动化装备有限公司 | 3d视觉引导拆垛测量系统及其控制方法 |
| CN110458822B (zh) * | 2019-08-07 | 2022-10-11 | 合肥学院 | 一种复杂曲面零件非接触式三维匹配检测方法 |
| CN110533624A (zh) * | 2019-09-11 | 2019-12-03 | 神博(山东)安防科技有限公司 | 一种激光测距辅助拼接系统和图像拼接方法 |
| CN112325795A (zh) * | 2020-10-16 | 2021-02-05 | 华中科技大学鄂州工业技术研究院 | 一种基于机器视觉引导的三维目标飞行时间测量方法、系统和装置 |
| CN112288824B (zh) * | 2020-10-27 | 2024-04-12 | 中国科学院上海微系统与信息技术研究所 | 一种基于现实场景的长焦相机标定装置及标定方法 |
| CN112288810B (zh) * | 2020-10-29 | 2023-04-07 | 铜陵有色金属集团股份有限公司 | 采样定位方法、装置、系统及计算机存储介质 |
| CN113155023B (zh) * | 2021-04-02 | 2023-02-07 | 甘肃旭盛显示科技有限公司 | 液晶基板玻璃翘曲度测量方法及系统 |
| CN113280790B (zh) * | 2021-05-07 | 2023-08-15 | 王亚超 | 一种基于激光定位的建筑测绘装置 |
| CN114894133A (zh) * | 2022-04-01 | 2022-08-12 | 上海钧控机器人有限公司 | 一种工具头非接触姿态测量与控制方法 |
| CN115811653B (zh) * | 2023-01-29 | 2023-05-16 | 苏州苏映视图像软件科技有限公司 | 一种成像系统、位置校准方法及成像方法 |
| CN117686034A (zh) * | 2024-01-31 | 2024-03-12 | 四川飞拓智控科技有限公司 | 一种多功能测量型相机 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000349988A (ja) * | 1999-06-07 | 2000-12-15 | Nkk Corp | 映像信号合成装置 |
| CN201242438Y (zh) * | 2008-08-12 | 2009-05-20 | 杨越三 | 一种激光标线测量仪 |
| CN102914262A (zh) * | 2012-09-29 | 2013-02-06 | 北京控制工程研究所 | 一种基于附加视距的非合作目标贴近测量方法 |
| CN104457569A (zh) * | 2014-11-27 | 2015-03-25 | 大连理工大学 | 一种大型复合板材几何参数视觉测量方法 |
| CN104897060A (zh) * | 2015-06-17 | 2015-09-09 | 大连理工大学 | 采用坐标跟踪控制板的大视场全局测量方法 |
| CN205352322U (zh) * | 2016-01-31 | 2016-06-29 | 山东科技大学 | 一种大型复杂曲面测量系统 |
-
2016
- 2016-01-31 CN CN201610070408.5A patent/CN105627948B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000349988A (ja) * | 1999-06-07 | 2000-12-15 | Nkk Corp | 映像信号合成装置 |
| CN201242438Y (zh) * | 2008-08-12 | 2009-05-20 | 杨越三 | 一种激光标线测量仪 |
| CN102914262A (zh) * | 2012-09-29 | 2013-02-06 | 北京控制工程研究所 | 一种基于附加视距的非合作目标贴近测量方法 |
| CN104457569A (zh) * | 2014-11-27 | 2015-03-25 | 大连理工大学 | 一种大型复合板材几何参数视觉测量方法 |
| CN104897060A (zh) * | 2015-06-17 | 2015-09-09 | 大连理工大学 | 采用坐标跟踪控制板的大视场全局测量方法 |
| CN205352322U (zh) * | 2016-01-31 | 2016-06-29 | 山东科技大学 | 一种大型复杂曲面测量系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105627948A (zh) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105627948B (zh) | 一种大型复杂曲面测量系统进行复杂曲面采样的方法 | |
| CN105157566B (zh) | 三维立体彩色点云扫描的方法 | |
| CN107421465B (zh) | 一种基于激光跟踪仪的双目视觉拼接方法 | |
| CN105698699B (zh) | 一种基于时间转轴约束的双目视觉测量方法 | |
| CN109000582A (zh) | 跟踪式三维扫描装置的扫描方法及系统、存储介质、设备 | |
| CN104034263B (zh) | 一种锻件尺寸的非接触测量方法 | |
| CN105115560B (zh) | 一种船舱舱容的非接触测量方法 | |
| CN107144241B (zh) | 一种基于景深补偿的双目视觉高精度测量方法 | |
| CN109215108A (zh) | 基于激光扫描的全景三维重建系统及方法 | |
| Guidi et al. | Fusion of range camera and photogrammetry: a systematic procedure for improving 3-D models metric accuracy | |
| CN106683068A (zh) | 一种三维数字化图像采集方法及设备 | |
| CN205352322U (zh) | 一种大型复杂曲面测量系统 | |
| CN109559349A (zh) | 一种用于标定的方法和装置 | |
| CN109141226A (zh) | 单相机多角度的空间点坐标测量方法 | |
| CN112927133B (zh) | 一种基于一体化标定参数的图像空间投影拼接方法 | |
| CN107633532B (zh) | 一种基于白光扫描仪的点云融合方法及系统 | |
| CN111189415B (zh) | 一种基于线结构光的多功能三维测量重建系统及方法 | |
| CN102927917A (zh) | 多目铁塔视觉测量方法 | |
| CN107339935A (zh) | 用于全视角扫描测量系统的靶标空间交会测量方法 | |
| CN108614277A (zh) | 双激光单摄像头三维成像扫描台及扫描、成像方法 | |
| CN112254670B (zh) | 一种基于光扫描和智能视觉融合的3d信息采集设备 | |
| CN108413865A (zh) | 基于三维测量和坐标系转换的二次反射镜面型检测方法 | |
| JP2004163271A (ja) | 非接触画像計測装置 | |
| Pulcrano et al. | 3D cameras acquisitions for the documentation of cultural heritage | |
| CN113432558B (zh) | 一种基于激光的不规则物体表面积的测量装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| 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: 20180206 |