CN111811396A - 一种基于平面约束的多边法激光跟踪三维坐标测量方法 - Google Patents

一种基于平面约束的多边法激光跟踪三维坐标测量方法 Download PDF

Info

Publication number
CN111811396A
CN111811396A CN202010530001.2A CN202010530001A CN111811396A CN 111811396 A CN111811396 A CN 111811396A CN 202010530001 A CN202010530001 A CN 202010530001A CN 111811396 A CN111811396 A CN 111811396A
Authority
CN
China
Prior art keywords
measuring
plane
measured
formula
station position
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
Application number
CN202010530001.2A
Other languages
English (en)
Other versions
CN111811396B (zh
Inventor
张福民
张画迪
曲兴华
周伦彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Jicheng Calibration Detection Authentication Co ltd
Tianjin University
Original Assignee
Suzhou Jicheng Calibration Detection Authentication Co ltd
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Jicheng Calibration Detection Authentication Co ltd, Tianjin University filed Critical Suzhou Jicheng Calibration Detection Authentication Co ltd
Priority to CN202010530001.2A priority Critical patent/CN111811396B/zh
Publication of CN111811396A publication Critical patent/CN111811396A/zh
Application granted granted Critical
Publication of CN111811396B publication Critical patent/CN111811396B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/002Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates

Abstract

本发明公开了一种基于平面约束的多边法激光跟踪三维坐标测量方法:(1)选定第一测站位置并布置光学面包板;(2)选定采样点,在第一测站位置对每一个采样点进行测量,获得每一个采样点在球坐标系下的三维坐标数据和每一个采样点与绝对测量式激光跟踪仪的测量中心之间的距离;(3)根据步骤(2)完成其他测站位置的数据测量;(4)根据步骤(2)得到的三维坐标数据,获得每块光学面包板的标准平面;(5)系统参数解算;(6)采用系统参数对待测物体的待测点坐标依照激光跟踪多边法进行解算。本发明通过在传统多边系统中引入平面约束,为系统提供优化方程,提高多边系统的系统参数解算精度,进而提高多边系统对空间点三维坐标测量的精度。

Description

一种基于平面约束的多边法激光跟踪三维坐标测量方法
技术领域
本发明涉及三维坐标测量技术,特别涉及一种基于平面约束的多边法激光跟踪三维坐标测量方法。
背景技术
工业测量是在工业生产和科研的各个环节中,为产品的设计、模拟、测量、放样、仿制、仿真、产品质量控制、产品运动状态等提供测量技术支撑的一门学科。激光以其良好的定向性和能量集中性,在工业测量领域有重要的应用。
激光跟踪仪是一种建立在球坐标上的测量设备,绝对测量式激光跟踪仪依照TOF与激光干涉原理相结合的方式进行距离测量,使用高精度角度编码器进行角度测量。具有精度高、便携性好等特点,在工业测量中使用十分广泛。然而在某些对精度有更高要求的测量场合,特别是在大尺寸测量中,激光跟踪仪角度测量精度相对距离测量精度较低的劣势就会显现出来。多边法系统能够只使用距离测量数据进行空间三维坐标计算从而避免角度测量带来的误差。
发明内容
本发明的目的是克服现有技术中的不足,提供一种基于平面约束的多边法激光跟踪三维坐标测量方法,通过在传统多边系统中引入平面约束,为系统提供优化方程,提高多边系统的系统参数解算精度,进而提高多边系统对空间点三维坐标测量的精度。
本发明所采用的技术方案是:一种基于平面约束的多边法激光跟踪三维坐标测量方法,所采用的测量装置包括绝对测量式激光跟踪仪和跟踪仪靶球,测量方法包括以下步骤:
步骤1,测量场搭建:选定能够测量待测物体的某一位置作为第一测站位置;在所述第一测站位置的设定范围内布置p块光学面包板;
步骤2,数据测量:选定每块所述光学面包板上的q个孔位点共p×q个点作为采样点,采用所述绝对测量式激光跟踪仪在第一测站位置配合所述跟踪仪靶球,对每一个所述采样点进行测量,获得每一个所述采样点在球坐标系下的三维坐标数据和每一个所述采样点与所述绝对测量式激光跟踪仪的测量中心之间的距离,完成第一测站位置的测量;
步骤3,转站测量:选定除第一测站位置以外能够测量待测物体的三处位置,作为第二测站位置、第三测站位置和第四测站位置,每个测站之间保持一距离,按照步骤2分别完成第二测站位置、第三测站位置和第四测站位置的测量;
步骤4,标准平面获取:根据步骤2中所获得的第一测站位置的每块光学面包板上的q个孔位点的三维坐标数据,对于每块所述光学面包板,在每块光学面包板上的q个孔位点中选取6个孔位点进行标准平面拟合,在全部平面结果中,选取平面度误差最小的平面,作为该块所述光学面包板的标准平面,标准平面所对应的6个孔位点为该标准平面上的目标点;
步骤5,系统参数解算:以p块所述光学面包板和4站测站位置组成的系统中,以第一测站位置的测量坐标系为系统坐标系,根据激光跟踪多边法求出第二测站位置、第三测站位置和第四测站位置在系统坐标系下的坐标;
步骤6,待测物体测量:采用步骤5得到的系统参数,对待测物体的待测点坐标依照激光跟踪多边法进行解算。
进一步地,步骤5中,所述的系统参数解算具体包括:
激光跟踪多边法使用平面约束下的距离方程作为解算空间点三维坐标的依据,已知平面约束下的距离方程为:
(xi-xj)2+(yi-yj)2+(zi-zj)2=lij 2 (1)
式中,(xi,yi,zi)为标准平面上第i个目标点的坐标,i=1,2,…,p×6;(xj,yj,zj)为第j个测站位置的坐标,j=1,2,3,4;lij为第i个目标点到第j个测站位置的距离;
将公式(1)写为函数形式如下:
Figure BDA0002534858070000021
式中,fij表示测量距离平方的平差,即理想距离平方和实际距离平方的差。
对公式(2)进行泰勒展开,略去高阶项后结果为:
Figure BDA0002534858070000031
式中,
Figure BDA0002534858070000032
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪测量得到;
Figure BDA0002534858070000033
是第j个测站位置的近似坐标,由各个测站位置对目标点的距离测量值使用后方交汇计算得到;
Figure BDA0002534858070000034
是fij的近似值,表达式如下所示:
Figure BDA0002534858070000035
式中,
Figure BDA0002534858070000036
为第j个测站位置对第i目标点距离的测量值,即第j个测站位置上的绝对测量式激光跟踪仪对第i个目标点测量得到;
对于标准平面约束上的点,存在平面方程:
Axi+Byi+Czi+D=0 (5)
式中,A,B,C,D为平面参数;
将公式(5)写为函数形式如下:
fpi=Axi+Byi+Czi+D (6)
式中,fpi是评定测量目标点坐标代入平面方程时的平差;
将公式(6)写为微分形式为:
Figure BDA0002534858070000037
式中,
Figure BDA0002534858070000038
为fpi的近似值,
Figure BDA0002534858070000039
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪测量得到;其中,m为标准平面序号,m=1,2,…,p;i为目标点序号,每个标准平面上有6个目标点,共计p×6个目标点,对于第m个标准平面,i=6×(m-1)+1,6×(m-1)+2,…,6×(m-1)+6;
将公式(3)和公式(7)写为矩阵形式如下:
Figure BDA00025348580700000310
Figure BDA00025348580700000311
式中,Fq为公式(3)中fij组成的列矩阵;Aq为公式(3)中未知数系数aij、bij、cij所构成的矩阵;dX为系统参数测量误差构成的矩阵,是球坐标系下测量值与实际值之差,dX=[dxi,dyi,dzi,…dxj,dyj,dzi,…]T,其中,(dxi,dyi,dzi)为
Figure BDA0002534858070000041
(dxj,dyj,dzj)为
Figure BDA0002534858070000042
Figure BDA0002534858070000043
为公式(3)中
Figure BDA0002534858070000044
组成的列矩阵;Fp为公式(7)中fpi组成的列矩阵;Ap为公式(7)中未知数系数Am、Bm、Cm所构成的矩阵;
Figure BDA0002534858070000045
为公式(7)中
Figure BDA0002534858070000046
组成的列矩阵;
将公式(8)和(9)合写为:
F=AdX+f0 (10)
式中,F为fpi和fij组成的列矩阵;A为方程组中参数系数矩阵,f0为平差矩阵,
Figure BDA0002534858070000047
采用系统参数测量误差对测站位置的近似坐标进行修正:
Figure BDA0002534858070000048
进一步地,步骤6中,所述的待测物体测量具体包括:
在测得系统参数之后按照式(12)组成方程组对待测物体的待测点坐标进行解算:
(xk-xj)2+(yk-yj)2+(zk-zj)2=rkj 2 (12)
式中,(xk,yk,zk)为待测物体上的待测点坐标;(xj,yj,zj)为步骤5得到的各个测站位置在系统坐标系下的坐标;rkj是通过测量得到的第k个待测点到第j个测站位置的距离。
本发明的有益效果是:本发明采用标准平面约束,有效提高激光跟踪仪对待测点测量精度。标准平面约束为系统在不增加未知数个数的情况下提供优化方程,使多边法测量系统能够在只有三台激光跟踪仪的情况下构建,降低成本。提出一种获取标准平面约束的方法,能够实现在测量现场通过光学面包板获得标准平面约束。
附图说明
图1:本发明一种基于平面约束的多边法激光跟踪三维坐标测量方法流程与图;
图2:本发明所涉及的测量装置布置图;
附图标注:1-绝对测量式激光跟踪仪,2-光学面包板,3-靶座,4-长度标准尺。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配合附图详细说明如下:
如图2所示,本发明所涉及的测量装置包括:绝对测量式激光跟踪仪1,跟踪仪靶球,光学面包板2,用于固定跟踪仪靶球的靶座3以及用于固定光学面包板2和绝对测量式激光跟踪仪1的三脚架。
所述绝对测量式激光跟踪仪1指拥有绝对测距能力进而得到空间测点三维坐标的激光跟踪仪,在本发明提出的测量方法中,需要四站绝对测量式激光跟踪仪1,可以同时使用四台绝对测量式激光跟踪仪1,也可以使用一台绝对测量式激光跟踪仪1在多个位置间进行移站操作。
所述跟踪仪靶球,一般拥有直角三棱锥棱镜结构,该结构能够将激光沿入射方向反射,保证入射光和反射光的平行。具有良好的球度保证跟踪仪靶球在靶座3上朝向不同方向时,球心位置尽可能保持不变。具有良好的制造精度保证棱镜结构的锥顶与球心尽可能重合。
所述靶座3,提供三个支点并具有吸附能力,能够将跟踪仪靶球固定在其上。
光学面包板2,提供孔位阵列,用于辅助跟踪仪靶球在其上定位。具有一定的平面度,提供本发明提出的测量方法中所需的标准平面约束的基底。
所述三角架,用于固定绝对测量式激光跟踪仪1与光学面包板2。
本发明中,绝对测量式激光跟踪仪1型号为LEICA AT960;跟踪仪靶球采用配套的LEICA 1.5”跟踪仪靶球;选取三块普通光学面包板2;标准尺三根,标称长度分别为969.045mm、970.108mm、1000.943mm。另使用三脚架若干,用于在测量场中布置绝对测量式激光跟踪仪1和光学面包板2。
如图1所示,一种基于平面约束的多边法激光跟踪三维坐标测量方法,通过在传统多边系统中引入平面约束,为系统提供优化方程,提高多边系统的系统参数解算精度,进而提高多边系统对空间点三维坐标测量的精度。包括以下步骤:
步骤1,测量场搭建:选定能够测量待测物体的某一位置作为第一测站位置;在所述第一测站位置的设定范围内布置p块光学面包板2。
本实施例中,将绝对测量式激光跟踪仪1放置在测量场中合适的位置,作为第一测站位置。在绝对测量式激光跟踪仪1 5m×5m的范围内布置三块光学面包板2,其中两块使用三脚架进行支撑,一块放置在地面,保证其固定没有晃动,三脚架上能有效避免晃动带来的误差。在绝对测量式激光跟踪仪1 10m×10m的范围内布置9个靶座3,作为系统的9个定向点位置,定向点可以使用在传统多边法系统参数解算中,也可以使用在平面约束多边法系统参数解算中,实际应用可以不设置靶座3。将被测物体放置在一定距离处,本实施例中,为验证本发明测量方法的精度,使用三根长度标准尺4作为被测物体,三根长度标准尺4分别放置在距离第一测站位置约25m、30m、38m的位置上,三根长度标准尺4标称长度分别为969.045mm、970.108mm、1000.943mm。
步骤2,数据测量:选定每块光学面包板2上的q个孔位点共p×q个点、作为采样点,采用所述绝对测量式激光跟踪仪1在第一测站位置配合所述跟踪仪靶球,对每一个所述采样点进行测量,获得每一个所述采样点在球坐标系下的三维坐标数据和每一个所述采样点与所述绝对测量式激光跟踪仪1的测量中心之间的距离,完成第一测站位置的测量。
本实施例中,使用绝对测量式激光跟踪仪1在第一测站位置配合跟踪仪靶球,分别对光学面包板2上的孔位点、定向点以及长度标准尺4两端点进行测量,分别保存所有点三坐标数据以及距离数据。其中每个光学面包板2上测量36个点,三块光学面包板2共108点,9个定向点以及6个长度标准尺4端点,共123点。
步骤3,转站测量:选定除第一测站位置以外能够测量待测物体的三处位置,作为第二测站位置、第三测站位置和第四测站位置,每个测站之间保持一定距离,按照步骤2分别完成第二测站位置、第三测站位置和第四测站位置的测量;
本实施例中,使用一台绝对测量式激光跟踪仪1在多个位置间进行移站操作。在第一站测站位置测量完毕后,移动绝对测量式激光跟踪仪1位置,并改变跟踪仪三脚架的高度,作为第二站测站位置,并重新测量全部点位。共进行三次移站操作,共对全部点位进行4次测量。
步骤4,标准平面获取:根据步骤2中所获得的第一测站位置的每块光学面包板2上的q个孔位点的三维坐标数据,对于每块光学面包板2,在每块光学面包板2上的q个孔位点中选取6个孔位点进行标准平面拟合,在全部平面结果中,选取平面度误差最小的平面,作为该块光学面包板2的标准平面,标准平面所对应的6个孔位点为该标准平面上的目标点。
本发明测量方法中要求提供标准平面作为约束使用,若使用平面度误差较大的平面,会给测量结果造成很大影响,因此,本发明提出一种能够在测量现场获取平面度较好平面的方法。分别利用光学面包板2本身的平面度以及其上提供的大量孔位。光学面包板2本身的平面度为寻找高平面度平面提供了基底;孔位能够帮助跟踪仪靶球在光学面包板2上定位。以121孔位的面包板,选取6点进行平面拟合为例,在其中按一定规律选取36个孔位,尽量覆盖整个面包板范围,计算所有选取6点后组成平面的平面度误差,共有
Figure BDA0002534858070000071
种可能。选取平面度误差最小的那一组孔位作为组成标准平面的点组。
本实施例中,在数据测量中获得了每块光学面包板2上的36点,3块共108点的坐标,在每块光学面包板2上的36点中选取6点进行标准平面拟合,已知
Figure BDA0002534858070000072
在全部1947792个平面结果中,选取平面度误差最小的平面,作为标准平面约束使用。选取6点进行平面拟合的是综合考虑平面度需求和测量点数需求的结果。实验表明,使用6点拟合时得到的三个最小平面度误差分别为7.37×10-5mm、1.91×10-4mm、2.33×10-4mm,均值1.66×10- 4mm,为百nm级。认为可以作为标准平面约束使用。
步骤5,系统参数解算:获得标准平面约束并得到全部采样点坐标以及距离测量值后,可进行系统参数解算。以p块光学面包板2和4站测站位置组成的系统中,以第一测站位置的测量坐标系为系统坐标系,根据激光跟踪多边法求出第二测站位置、第三测站位置和第四测站位置在系统坐标系下的坐标。
以第一测站位置为原点,第一测站位置的测量坐标系为系统坐标系。如果有m个测站,n个目标点。这样能引入m×n个距离方程,n个平面方程。m个测站引入3m-3个未知数,n个目标点引入3n个未知数,若要方程组有解,需满足m×n+n≥3m+3n-3,即
Figure BDA0002534858070000081
在m、n均取正整数的情况下,m最小取3,n最小取6。
激光跟踪多边法使用平面约束下的距离方程作为解算空间点三维坐标的依据,已知平面约束下的距离方程可以写为:
(xi-xj)2+(yi-yj)2+(zi-zj)2=lij 2 (1)
式中,(xi,yi,zi)为标准平面上第i个目标点的坐标,i=1,2,…,p×6;(xj,yj,zj)为第j个测站位置的坐标,j=1,2,3,4;lij为第i个目标点到第j个测站位置的距离。
不难发现,要使用多边法测量空间三维坐标,需要知道各测站之间的位置关系,也就是系统参数。传统多边法采用设置定向点的方法,利用引入定向点时引入的方程数大于未知数个数的特点,对系统参数解算。本发明提出的平面约束方案,使得定向点能够额外引入平面方程,达到优化系统的目的。另外由于方程数进一步大于未知数个数,使得减少测站数量成为可能。绝对测量式激光跟踪仪1在一个测站位置完成对所有测点的测量之后,需要移动位置,作为下一个测站,直至所需要的测站数量都被测量完成。因而减少测站的方案能够大幅度减少多边法测量时的工作量。
将公式(1)写为函数形式如下:
Figure BDA0002534858070000086
式中,fij表示测量距离平方的平差,即理想距离平方和实际距离平方的差。
对公式(2)进行泰勒展开,略去高阶项后结果为:
Figure BDA0002534858070000082
式中,
Figure BDA0002534858070000083
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪1测量得到;
Figure BDA0002534858070000084
是第j个测站位置的近似坐标,由各个测站位置对目标点的距离测量值使用后方交汇计算得到;
Figure BDA0002534858070000085
是fij的近似值,表达式如下所示:
Figure BDA0002534858070000091
式中,
Figure BDA0002534858070000092
为第j个测站位置对第i目标点距离的测量值,即第j个测站位置上的绝对测量式激光跟踪仪对第i个目标点测量得到;
对于标准平面约束上的点,存在平面方程:
Axi+Byi+Czi+D=0 (5)
式中,A,B,C,D为平面参数;
将公式(5)写为函数形式如下:
fpi=Axi+Byi+Czi+D (6)
式中,fpi是评定测量目标点坐标代入平面方程时的平差。
将公式(6)写为微分形式为:
Figure BDA0002534858070000093
式中,
Figure BDA0002534858070000094
为fpi的近似值,
Figure BDA0002534858070000095
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪1测量得到;其中,m为标准平面序号,m的取值为1,2,…,p;i为目标点序号,每个标准平面上有6个目标点,共计p×6个目标点,对于第m个标准平面,i的取值依次为6×(m-1)+1,6×(m-1)+2,…,6×(m-1)+6。
以上表达式中的所有目标点近似坐标
Figure BDA0002534858070000096
均由第一测站位置上的绝对测量式激光跟踪仪1直接测得,所有测站的近似坐标
Figure BDA0002534858070000097
由各测站对目标点的距离测量值使用后方交汇计算得到。
将公式(3)和公式(7)写为矩阵形式如下:
Figure BDA0002534858070000098
Figure BDA0002534858070000099
式中,Fq为公式(3)中fij组成的列矩阵;Aq为公式(3)中未知数系数aij、bij、cij所构成的矩阵;dX为系统参数测量误差构成的矩阵,是球坐标系下测量值与实际值之差,dX=[dxi,dyi,dzi,…dxj,dyj,dzi,…]T,其中,(dxi,dyi,dzi)为
Figure BDA0002534858070000101
(dxj,dyj,dzj)为
Figure BDA0002534858070000102
Figure BDA0002534858070000103
为公式(3)中
Figure BDA0002534858070000104
组成的列矩阵;Fp为公式(7)中fpi组成的列矩阵;Ap为公式(7)中未知数系数Am、Bm、Cm所构成的矩阵;
Figure BDA0002534858070000105
为公式(7)中
Figure BDA0002534858070000106
组成的列矩阵;
公式(3)(7)联立可以得到线性方程组,矩阵形式即为将公式(8)和(9)合写:
F=AdX+f0 (10)
式中,F为fpi和fij组成的列矩阵;A为方程组中参数系数矩阵,f0为平差矩阵,
Figure BDA0002534858070000107
本实施例中,以4测站,3标准平面约束,每标准平面上6个目标点,另设15个目标点(包括9个定向点以及6个长度标准尺4端点)共33目标点为例,结果如下:
Figure BDA0002534858070000108
dX=dX111×1=[dxi,dyi,dzi,…dxj,dyj,dzi,…]T
Figure BDA0002534858070000109
利用上述矩阵形式对超定方程组求最小二乘解可得到各系统参数,即各测站间的位置关系,进而可以使用多边法对待测空间点三维坐标进行计算。
由上述步骤得到的系统参数误差对系统参数的近似值进行修正,得到较为准确的系统参数后即可使用这些系统参数,即准确的测站间位置关系对待测目标的坐标依照多边法基本原理进行解算。
Figure BDA0002534858070000111
步骤6,待测物体测量与进度评定:采用步骤5得到的系统参数,对待测物体的坐标依照激光跟踪多边法进行解算。
在测得系统参数即各测站之间的位置关系之后可按照式(12)组成方程组对待测物体的待测点坐标进行解算:
(xk-xj)2+(yk-yj)2+(zk-zj)2=rkj 2 (12)
式中,(xk,yk,zk)为待测物体上的待测点坐标;(xj,yj,zj)为步骤5得到的各个测站位置在系统坐标系下的坐标;rkj是通过测量得到的第k个待测点到第j个测站位置的距离。
本实施例中,为验证系统的精度,使用三根标称长度分别为969.045mm、970.108mm、1000.943mm的标准尺作为被测物体,并分别放置在约25m、30m、38m的距离上,验证系统在不同测量距离上的精度表现。在得到长度标准尺4两端点坐标后,可对长度标准尺4的长度进行计算,由本发明的案例进行计算得到三根标称长度为969.045mm、970.108mm、1000.943mm的长度标准尺4的测量长度分别为969.035mm,970.081mm,1000.894mm。得到的基于平面约束的多边法激光跟踪测量方法测量误差分别为0.01mm,0.027mm,0.049mm,平均值0.029mm。而使用传统四站多边法进行测量时,得到的平均误差为0.056mm,而使用单站测量时,平均误差更是达到了0.085mm。可知平面约束能够有效的提高系统测量精度。
综上,本发明提出了一种基于平面约束的多边法激光跟踪三维坐标测量方法,同时,提出了一种在测量现场获取标准平面约束的方法,为上述测量方法的实现提供可能。测量方法所采用的测量装置包括绝对测量式激光跟踪仪1,1.5”跟踪仪靶球,光学面包板2,长度标准尺4,三脚架等设备。多边法激光跟踪三维坐标测量系统能够避免绝对测量式激光跟踪仪1角度测量带来的误差,只使用距离测量数据解算空间点的三维坐标。对多边系统参数的准确标定是提高空间三维坐标测量精度的关键。为提高多边法测量系统的精度,本发明提出一种基于平面约束的多边法激光跟踪三维坐标测量方法,采用标准平面约束,为系统在不增加未知数个数的情况下提供优化方程。为能够实现在测量现场获得标准平面约束,又提出一种获取标准平面约束的方法。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。

Claims (3)

1.一种基于平面约束的多边法激光跟踪三维坐标测量方法,所采用的测量装置包括绝对测量式激光跟踪仪(1)和跟踪仪靶球,其特征在于,测量方法包括以下步骤:
步骤1,测量场搭建:选定能够测量待测物体的某一位置作为第一测站位置;在所述第一测站位置的设定范围内布置p块光学面包板(2);
步骤2,数据测量:选定每块所述光学面包板(2)上的q个孔位点共p×q个点作为采样点,采用所述绝对测量式激光跟踪仪(1)在第一测站位置配合所述跟踪仪靶球,对每一个所述采样点进行测量,获得每一个所述采样点在球坐标系下的三维坐标数据和每一个所述采样点与所述绝对测量式激光跟踪仪(1)的测量中心之间的距离,完成第一测站位置的测量;
步骤3,转站测量:选定除第一测站位置以外能够测量待测物体的三处位置,作为第二测站位置、第三测站位置和第四测站位置,每个测站之间保持一距离,按照步骤2分别完成第二测站位置、第三测站位置和第四测站位置的测量;
步骤4,标准平面获取:根据步骤2中所获得的第一测站位置的每块光学面包板(2)上的q个孔位点的三维坐标数据,对于每块所述光学面包板(2),在每块光学面包板(2)上的q个孔位点中选取6个孔位点进行标准平面拟合,在全部平面结果中,选取平面度误差最小的平面,作为该块所述光学面包板(2)的标准平面,标准平面所对应的6个孔位点为该标准平面上的目标点;
步骤5,系统参数解算:以p块所述光学面包板(2)和4站测站位置组成的系统中,以第一测站位置的测量坐标系为系统坐标系,根据激光跟踪多边法求出第二测站位置、第三测站位置和第四测站位置在系统坐标系下的坐标;
步骤6,待测物体测量:采用步骤5得到的系统参数,对待测物体的待测点坐标依照激光跟踪多边法进行解算。
2.根据权利要求1所述的一种基于平面约束的多边法激光跟踪三维坐标测量方法,其特征在于,步骤5中,所述的系统参数解算具体包括:
激光跟踪多边法使用平面约束下的距离方程作为解算空间点三维坐标的依据,已知平面约束下的距离方程为:
(xi-xj)2+(yi-yj)2+(zi-zj)2=lij 2 (1)
式中,(xi,yi,zi)为标准平面上第i个目标点的坐标,i=1,2,…,p×6;(xj,yj,zj)为第j个测站位置的坐标,j=1,2,3,4;lij为第i个目标点到第j个测站位置的距离;
将公式(1)写为函数形式如下:
Figure FDA0002534858060000021
式中,fij表示测量距离平方的平差,即理想距离平方和实际距离平方的差。
对公式(2)进行泰勒展开,略去高阶项后结果为:
Figure FDA0002534858060000022
式中,
Figure FDA0002534858060000023
Figure FDA0002534858060000024
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪(1)测量得到;
Figure FDA0002534858060000025
是第j个测站位置的近似坐标,由各个测站位置对目标点的距离测量值使用后方交汇计算得到;
Figure FDA0002534858060000026
是fij的近似值,表达式如下所示:
Figure FDA0002534858060000027
式中,
Figure FDA0002534858060000028
为第j个测站位置对第i目标点距离的测量值,即第j个测站位置上的绝对测量式激光跟踪仪(1)对第i个目标点测量得到;
对于标准平面约束上的点,存在平面方程:
Axi+Byi+Czi+D=0 (5)
式中,A,B,C,D为平面参数;
将公式(5)写为函数形式如下:
fpi=Axi+Byi+Czi+D (6)
式中,fpi是评定测量目标点坐标代入平面方程时的平差;
将公式(6)写为微分形式为:
Figure FDA0002534858060000029
式中,
Figure FDA00025348580600000210
为fpi的近似值,
Figure FDA00025348580600000211
是标准平面上第i目标点的近似坐标,由第一测站位置上的绝对测量式激光跟踪仪(1)测量得到;其中,m为标准平面序号,m=1,2,…,p;i为目标点序号,每个标准平面上有6个目标点,共计p×6个目标点,对于第m个标准平面,i=6×(m-1)+1,6×(m-1)+2,…,6×(m-1)+6;
将公式(3)和公式(7)写为矩阵形式如下:
Figure FDA0002534858060000031
Figure FDA0002534858060000032
式中,Fq为公式(3)中fij组成的列矩阵;Aq为公式(3)中未知数系数aij、bij、cij所构成的矩阵;dX为系统参数测量误差构成的矩阵,是球坐标系下测量值与实际值之差,dX=[dxi,dyi,dzi,…dxj,dyj,dzi,…]T,其中,(dxi,dyi,dzi)为
Figure FDA0002534858060000033
(dxj,dyj,dzj)为
Figure FDA0002534858060000034
Figure FDA0002534858060000035
为公式(3)中
Figure FDA0002534858060000036
组成的列矩阵;Fp为公式(7)中fpi组成的列矩阵;Ap为公式(7)中未知数系数Am、Bm、Cm所构成的矩阵;
Figure FDA0002534858060000037
为公式(7)中
Figure FDA0002534858060000038
组成的列矩阵;
将公式(8)和(9)合写为:
F=AdX+f0 (10)
式中,F为fpi和fij组成的列矩阵;A为方程组中参数系数矩阵,f0为平差矩阵,
Figure FDA0002534858060000039
采用系统参数测量误差对测站位置的近似坐标进行修正:
Figure FDA00025348580600000310
3.根据权利要求1所述的一种基于平面约束的多边法激光跟踪三维坐标测量方法,其特征在于,步骤6中,所述的待测物体测量具体包括:
在测得系统参数之后按照式(12)组成方程组对待测物体的待测点坐标进行解算:
(xk-xj)2+(yk-yj)2+(zk-zj)2=rkj 2 (12)
式中,(xk,yk,zk)为待测物体上的待测点坐标;(xj,yj,zj)为步骤5得到的各个测站位置在系统坐标系下的坐标;rkj是通过测量得到的第k个待测点到第j个测站位置的距离。
CN202010530001.2A 2020-06-11 2020-06-11 一种基于平面约束的多边法激光跟踪三维坐标测量方法 Active CN111811396B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010530001.2A CN111811396B (zh) 2020-06-11 2020-06-11 一种基于平面约束的多边法激光跟踪三维坐标测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010530001.2A CN111811396B (zh) 2020-06-11 2020-06-11 一种基于平面约束的多边法激光跟踪三维坐标测量方法

Publications (2)

Publication Number Publication Date
CN111811396A true CN111811396A (zh) 2020-10-23
CN111811396B CN111811396B (zh) 2021-05-25

Family

ID=72845746

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010530001.2A Active CN111811396B (zh) 2020-06-11 2020-06-11 一种基于平面约束的多边法激光跟踪三维坐标测量方法

Country Status (1)

Country Link
CN (1) CN111811396B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112665523A (zh) * 2020-11-24 2021-04-16 北京星航机电装备有限公司 一种复杂型面的联合测量方法
CN112683163A (zh) * 2020-12-03 2021-04-20 中国科学院光电技术研究所 一种适用于视觉测量相机的绝对位姿测量精度评估方法
CN113587819A (zh) * 2021-09-01 2021-11-02 哈尔滨理工大学 一种基于三路激光跟踪的大尺度空间位姿动态测量方法与测量精度验证方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090125265A1 (en) * 2008-12-06 2009-05-14 Parker David H Method for Calibrating a Laser-Based Spherical Coordinate Measurement System by a Mechanical Harmonic Oscillator
CN201463834U (zh) * 2009-05-22 2010-05-12 南京航空航天大学 用于大型曲面的激光跟踪仪坐标定位辅助装置
CN201867495U (zh) * 2010-11-19 2011-06-15 二重集团(德阳)重型装备股份有限公司 用于对静态工件进行三维空间定位的装置及专用手持活动式光学逆反射器
CN104240300A (zh) * 2014-08-29 2014-12-24 电子科技大学 基于分布式并行的大规模点云复杂空间曲面重构方法
CN104677280A (zh) * 2015-03-11 2015-06-03 中国科学院光电技术研究所 摆臂式轮廓仪旋转轴空间状态标定方法
CN106052555A (zh) * 2016-06-08 2016-10-26 中国科学院重庆绿色智能技术研究院 一种工业机器人基座坐标测量方法
CN106247932A (zh) * 2016-07-25 2016-12-21 天津大学 一种基于摄影系统的机器人在线误差补偿装置及方法
CN108344361A (zh) * 2018-01-30 2018-07-31 清华大学 一种基于激光跟踪仪的平面法向量测量方法
CN108413988A (zh) * 2018-03-13 2018-08-17 北京卫星环境工程研究所 机器人末端经纬仪坐标系快速标定方法
CN110926337A (zh) * 2019-12-24 2020-03-27 大连理工大学 一种基于多矢量约束的全局测量数据配准方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090125265A1 (en) * 2008-12-06 2009-05-14 Parker David H Method for Calibrating a Laser-Based Spherical Coordinate Measurement System by a Mechanical Harmonic Oscillator
CN201463834U (zh) * 2009-05-22 2010-05-12 南京航空航天大学 用于大型曲面的激光跟踪仪坐标定位辅助装置
CN201867495U (zh) * 2010-11-19 2011-06-15 二重集团(德阳)重型装备股份有限公司 用于对静态工件进行三维空间定位的装置及专用手持活动式光学逆反射器
CN104240300A (zh) * 2014-08-29 2014-12-24 电子科技大学 基于分布式并行的大规模点云复杂空间曲面重构方法
CN104677280A (zh) * 2015-03-11 2015-06-03 中国科学院光电技术研究所 摆臂式轮廓仪旋转轴空间状态标定方法
CN106052555A (zh) * 2016-06-08 2016-10-26 中国科学院重庆绿色智能技术研究院 一种工业机器人基座坐标测量方法
CN106247932A (zh) * 2016-07-25 2016-12-21 天津大学 一种基于摄影系统的机器人在线误差补偿装置及方法
CN108344361A (zh) * 2018-01-30 2018-07-31 清华大学 一种基于激光跟踪仪的平面法向量测量方法
CN108413988A (zh) * 2018-03-13 2018-08-17 北京卫星环境工程研究所 机器人末端经纬仪坐标系快速标定方法
CN110926337A (zh) * 2019-12-24 2020-03-27 大连理工大学 一种基于多矢量约束的全局测量数据配准方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴斌等: "激光跟踪绝对测长多边法三维坐标测量系统", 《红外与激光工程》 *
谢政委等: "基于空间长度约束的坐标控制场精度增强方法", 《中国激光》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112665523A (zh) * 2020-11-24 2021-04-16 北京星航机电装备有限公司 一种复杂型面的联合测量方法
CN112665523B (zh) * 2020-11-24 2022-04-19 北京星航机电装备有限公司 一种复杂型面的联合测量方法
CN112683163A (zh) * 2020-12-03 2021-04-20 中国科学院光电技术研究所 一种适用于视觉测量相机的绝对位姿测量精度评估方法
CN112683163B (zh) * 2020-12-03 2022-06-28 中国科学院光电技术研究所 一种适用于视觉测量相机的绝对位姿测量精度评估方法
CN113587819A (zh) * 2021-09-01 2021-11-02 哈尔滨理工大学 一种基于三路激光跟踪的大尺度空间位姿动态测量方法与测量精度验证方法

Also Published As

Publication number Publication date
CN111811396B (zh) 2021-05-25

Similar Documents

Publication Publication Date Title
CN111811396B (zh) 一种基于平面约束的多边法激光跟踪三维坐标测量方法
CN107703499B (zh) 一种基于自制地基激光雷达对准误差的点云误差校正方法
CN106524905B (zh) 一种基于激光追踪仪多站位测量的四轴机床标定方法
CN103454619B (zh) 一种星载微波跟瞄雷达的电轴光学标定系统及其标定方法
CN112415493B (zh) 三维扫描激光雷达坐标误差修正方法
CN108413988B (zh) 机器人末端经纬仪坐标系快速标定方法
CN108489401A (zh) 分体式校准靶标、具有该靶标的校准装置及其校准方法
CN110487182B (zh) 一种基于多维动态标准器的坐标转换方法
CN101539397B (zh) 物体三维姿态的精密光学测量方法
CN109163658B (zh) 一种可提供位置和角度基准的光学基准件的标定方法
CN108534801A (zh) 三维坐标基准场室内空间测量定位扫描光面校准方法
CN111595238B (zh) 一种基于多站法的激光跟踪仪精度现场评价系统
CN114894094B (zh) 一种空间多靶标布设的激光扫描仪指标评估方法
CN110966934B (zh) 一种调姿器球窝中心点的标定方法
CN114485392B (zh) 基于激光跟踪仪的航天器大尺寸机械基准建立方法和系统
CN109884658A (zh) 基于激光追踪仪多站位测量系统的激光追踪仪站位方法
CN107588929B (zh) 球幕投影/跟踪系统标定方法及标定器
CN115371545A (zh) 激光跟踪仪姿态测量校准装置及方法
CN107727118B (zh) 大型飞行器中的gnc分系统设备姿态测量系统标定方法
CN111189391B (zh) 一种基于测量碳纤维杆轴线中点的坐标统一化方法
CN113587819A (zh) 一种基于三路激光跟踪的大尺度空间位姿动态测量方法与测量精度验证方法
CN110779503B (zh) 一种三维精密控制网测量方法
CN107991684B (zh) 大型飞行器中的gnc分系统设备姿态测量系统
CN113513986B (zh) 几何公差测量设备及其测量方法
CN112050731B (zh) 基于大型虚拟计量器的测量数据配准精度评价方法

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