专利名称：一种基于准直光路的角秒级三维变形测量装置与方法
技术领域：
本发明属于光学变形测量技术领域，特别涉及一种基于准直光路的角秒级三维变形测量装置与方法。
背景技术：
光学变形测量技术可广泛用于大型运动平台变形测量和高精度姿态基准传递等方面，具有测量精度高、直接测量、稳定性好等优点。目前的光学变形测量大多采用准直光路方法实现。中国专利CN97251378. 7、美国专利US7227627、《中国惯性技术学报》杂志2006年第3期中“一种舰船变形测量方法的研究与应用”和《长春理工大学学报》杂志2006年第 3期中“车载平台变形的激光自准直测量方法研究”论文公开了一类基于自准直光路、平面反射镜和面阵CXD (Charge Coupled Device，电荷耦合器件)的二维变形测量装置，该装置如图I所示，主要由光源I、光阑2、分光镜3、面阵CCD探测器4、光学系统5和平面反射镜6构成。该装置可以实时测量两个载体之间的相对变形，为了方便起见将其中一个载体称为待测物，另外一个载体称为参考基准。可将该测量装置的平面反射镜6与待测物固连，其余部分与参考基准固连，反之亦然。这种测量装置的测量光学变形的方法是利用光源I照射光阑2，透过光阑2的光部分透过分光镜3，然后通过光学系统5对光阑2进行成像，所成的光阑像经过平面反射镜6、分光镜3和光学系统5后成像在面阵CXD探测器4上，根据面阵CCD探测器的光阑图像位置解算变形。该装置可以测量与光轴垂直方向的二维变形，不能测量绕光轴方向的旋转角变形。该装置的二维变形测量精度可以达到角秒量级。美国专利US7697127公开了一种测量绕光轴旋转角的装置，该装置如图2所示，主要由光源I、光阑2和面阵探测器8构成。该装置的光阑2采用周期结构的光阑，由光阑2和面阵探测器8实现绕光轴7旋转角测量，具体测量方法在光源I照射下光阑2直接投影到面阵探测器8上，利用图像处理方法处理面阵探测器8探测光阑图像得到旋转角变化。该装置可以测量绕光轴7的旋转角，但不能测量与光轴垂直方向的角度变化。由于没有准直光路和直接投影的图像质量限制，该装置的旋转角测量精度难以达到角秒量级。中国专利CN200410010856. 3、《计量技术》杂志2005年第4期中“干涉法在船体角度变形测量中的应用”、《光学技术》杂志2005年第5期中“船体三维角度变形的自准直干涉测量方法”和《飞行器测控学报》2009年第3期中“基于莫尔条纹测量法设计的船体变形测量系统”论文公开了一类干涉类型的三维变形测量装置，该装置组成与图I所示装置类似，如图3所不，由第一光源I、光阑2、第一分光棱镜3、第二分光棱镜12、第三分光棱镜7、光学系统5、反射镜6、第一面阵(XD探测器41、第二光源9、自准直光栅10、猫准光栅11和第二面阵CCD探测器42构成。该装置的测量方法是利用第一光源I、光阑2、第一分光棱镜3、第二分光棱镜12、第三分光棱镜7、光学系统5、反射镜6、第一面阵(XD探测器41构成一套自准直光路，采用与图I类似的方法测量与光轴垂直方向的二维变形角；利用第二光源9、自准直光栅10、第二分光棱镜12、第三分光棱镜7、光学系统5、反射镜6、猫准光栅11和第二面阵CCD探测器42构成一套光栅干涉光路，采用莫尔条纹方法测量绕光轴方向的旋转角变形。该方法测量的与光轴垂直方向的二维变形(方位和俯仰方向变形)精度较高，能达到角秒量级；绕光轴方向的变形(横滚方向变形)测量精度相对较低，难以达到角秒量级。该装置采用部件较多且结构较复杂，具有制作难度大、稳定性相对较低和成本高等特点。《光学精密工程》杂志2008年第11期中“船体变形测量的新技术”和《长春理工大学学报》杂志2006年第3期中“基于测量船角变形光电测量系统的研究”论文介绍了一种采用大钢管的三维变形测量装置，该装置的主要结构如图4所示，主要由第一发射光管13、第二发射光管14、第一接收光管19、第二接收光管18、大钢管16、钢管安装座24、第一测扭光管20、第一反射镜15、第二测扭光管21和第二反射镜17构成，用于测量待测物22与参考基准23之间的三维变形。其中，第一发射光管13、第二发射光管14、第一测扭光管20、第二测扭光管21和第二反射镜17都与待测物22固连；第一接收光管19、第二接收光管18、第二测扭光管21都与参考基准23固连；第一反射镜15与大钢管16的一端固连，第二反射镜17与大钢管16的另一端固连；大钢管16采用特殊材料设计加工，具有温度稳定性好和 横滚方向变形小等特点；大钢管16由钢管安装座24安装，钢管安装座24能够支撑大钢管16并且不会使大钢管16在横滚方向产生变形。在第一接收光管19、第二接收光管18、第一测扭光管20和第二测扭光管21内部都分别安装了线阵CCD探测器。该装置通过第一发射光管13和第一接收光管19构成一套自准直光路，测量待测物22与参考基准之间的俯仰方向变形；通过第二发射光管14和第二接收光管18构成一套自准直光路，测量待测物22与参考基准之间的方位方向变形；通过第一测扭光管20和第一反射镜15测量待测物22与大钢管16 —端之间在横滚方向的变形，通过第二测扭光管21和第二反射镜17测量参考基准23与大钢管16另一端在横滚方向的变形，由此计算待测物22与参考基准23在横滚方向的变形。该装置能够实时测量待测物22与参考基准23之间的三维变形，方位、俯仰和横滚方向变形的测量精度都可以达到角秒量级，但该装置具有结构复杂、成本高、加工难度大和占用空间大等不足。如何既能实时测量待测物和参考基准之间的三维变形(方位、俯仰和横滚方向变形)，又具有测量精度高、结构简单、占用空间小、体积小和成本低等特点，是本领域技术人员极为关注的技术问题。
发明内容
本发明要解决的技术问题是提出一种基于准直光路的角秒级三维变形测量装置与方法，可以同时测量三个方向(方位、俯仰和横滚方向)变形，并且具有测量精度高、结构简单、占用空间小和成本低等优点。本发明三维变形测量装置由发射�？椤⒔邮漳？楹痛砟？槿蟛糠止钩�。发射�？橛氩慰蓟脊塘邮漳？橛氪馕锕塘７⑸淠？橛虢邮漳？橹淠芄煌ㄊ樱⑸淠？榉⑸涞墓馐苯咏虢邮漳？�。该装置能实时测量待测物与参考基准的三维变形。发射�？橛晒庠�、光阑和发射光学系统构成。采用光学领域常用的光学固定座和调整架将光源、光阑和发射光学系统固定在参考基准上。光阑位于光源和发射光学系统之间，且位于发射光学系统的焦平面上，光源和光阑的中心都在发射光学系统的光轴上，由此构成发射准直光路。光源照射光阑，发射光学系统将透过光阑的光束准直发射到接收�？橹�。光源可以是半导体激光器或LED或灯泡等。发射光学系统可以采用双胶合透镜或双分离透镜或多片透镜组等类型光学系统。光阑是在玻璃等透明材料上利用光刻等工艺加工而成的，加工后的光阑2有些部分透光，有些部分不透光。光阑2的透光部分采用LJL。为大于I的常数)条直线结构，各直线呈放射状分布，直线结构为透光部分，其余部分不透光。两直线结构光阑(如图6(a)所示)的透光部分由两条直线构成，两直线可互相垂直或不垂直。每条直线由两个线段构成，其中线段Lla、Llb构成直线LI，线段L2a、L2b构成直线L2。采用多直线结构光阑可以提高变形测量精度，三条直线结构光阑如6(b)所示，四条直线结构光阑如6(c)所示，四条以上的直线结构光阑依次类推，直线越多，测量精度越高，但横滚方向的形变测量范围会减小。接收�？橛山邮展庋低澈兔嬲筇讲馄鞴钩�。采用光学领域常用的光学固定座和调整架将接收光学系统和面阵探测器固定在待测物上。面阵探测器位于接收光学系统的焦平面上，且其中心位于接收光学系统的光轴上，由此构成接收准直光路。接收光学系统可采 用双胶合透镜或双分离透镜或多片透镜组等类型光学系统，面阵探测器可采用灰度或彩色的面阵CCD或CMOS等类型探测器。发射�？榉⑸涞淖贾惫馐邮展庋低郴峋酆蟪上裨诿嬲筇讲馄魃希嬲筇讲馄魈讲獾氖枪饫坏闹毕呓峁构饫煌枷�。处理�？榭刹捎霉た鼗�、嵌入式计算机、台式计算机或笔记本电脑等设备，在处理�？橹邪沧傲吮湫谓馑闳砑�。处理模块通过千兆以太网、高速USB或1394等接口与面阵探测器连接，实时采集面阵探测器探测的光阑图像，变形解算软件计算三维变形。按图5所不方式定义坐标系，即发射�？楣庵岱较蛭猋轴，与Y轴垂直且朝上的方向为Z轴，X轴与Y轴、Z轴垂直。当参考基准和待测物之间发生变形时，发射模块与接收�？橹浠岱⑸灞湫巍Ｏ喽杂诜⑸淠？椋邮漳？檎迦芞轴旋转产生方位角变形，绕Y轴旋转产生横滚角变形，绕X轴旋转产生俯仰角变形。当发射模块的光轴与接收�？榈墓庵崞叫星曳⑸淠？橹泄饫坏乃街毕哂虢邮漳？榈拿嬲筇讲馄鞯乃街行南咂叫惺保⑸淠？橛虢邮漳？橹涞姆轿�、俯仰和横滚方向变形都为零。此时，面阵探测器探测的光阑图像中心与面阵探测器中心重合，光阑图像中直线结构的某一直线(称为水平直线)与面阵探测器的水平中心线平行，如图7所示(图中以两直线结构光阑为例，与面阵探测器的水平中心线平行的直线为水平直线)。当发射�？橛虢邮漳？橹浯嬖诜轿缓透┭龇较虮湫问保嬲筇讲馄魈讲獾墓饫煌枷穹⑸揭疲饫煌枷裰行钠朊嬲筇讲馄髦行模缤�8所示(图中以两直线结构光阑为例)。将光阑图像中心与面阵探测器中心在X轴方向的偏移量记为ΛΧ，Z轴方向的偏移量记为Λ Z ;将接收�？楣庋低车慕咕嗉俏猣。根据ΛΧ、ΛΖ和f可以计算方位方向变形V、俯仰方向变形Θ。当发射�？橛虢邮漳？橹浯嬖诤峁龇较虮湫问保嬲筇讲馄鞯墓饫煌枷穹⑸缤�9所示(图中两直线结构光阑为例)。光阑图像中直线方向与面阵探测器中心线存在夹角P，该角度为横滚方向变形。采用本发明所述装置进行变形测量的具体方法如下第一步，发射�？榈墓庠凑丈涔饫唬⑸涔庋低辰腹饫坏墓馐贾狈⑸涞浇邮漳？椤５诙剑邮漳？榈慕邮展庋低辰邮辗⑸淠？榉⑸涞墓馐⒔浠峋鄣矫嬲筇讲馄魃�。
第三步，面阵探测器将接收的光束转换成光阑图像，此光阑图像为数字化的彩色图像或灰度图像。第四步，处理�？榇用嬲筇讲馄鹘邮展饫煌枷�。第五步，处理�？橹斜湫谓馑闳砑莨饫煌枷癜匆韵虏街杓扑闳湫�5. I若面阵探测器是彩色类型探测器，则将彩色图像转换成灰度图像，可参照RafaelC. Gonzalez等编著、阮秋琦等翻译、电子工业出版社2005年出版的《数字图像处理(MATLAB版)》中151 160页介绍的方法进行转换；若面阵探测器是灰度类型探测器，则直接进入步骤5. 2。5. 2对光阑图像进行图像校正，可采用R. Y. Tsai在1986年IEEE Proceedings ofConference on Computer Vision and Pattern Recognition 发表的 An Efficient andAccurate Camera Calibration Technique For 3D Machine Vision 论文介绍的传统标 定或张正友在 2000 年第 11 期 IEEE Transactions on Pattern Analysis and MachineIntelligence发表的A Flexible New Technique for Camera Calibration论文介绍的自标定方法对光阑图像进行校正。发射光学系统畸变、接收光学系统畸变和面阵探测器畸变等会导致光阑图像产生畸变，对光阑图像的畸变进行校正可提高变形测量精度。5. 3对校正后的光阑图像进行亮度分析，并滤除背景噪声。对光阑图像的所有像素根据亮度进行排序，根据光阑的直线结构的直线数量L。选取L3个(L3通常为L。的10 100倍)最亮像素，根据选择的最亮像素计算平均亮度，将该平均亮度作为光阑图像的最大亮度I1AP根据最大亮度Ism和面阵探测器噪声大小设置亮度阈值ITH，亮度阈值I 大于面阵探测器噪声，且小于最大亮度Imax的20%。根据Ith滤除背景噪声，将小于Ith的所有像素的亮度置为零。5. 4计算光阑图像的各像素行的各直线中心点坐标。滤除背景噪声后，光阑图像的每条直线在包含该直线的图像像素行上对应一段非零的亮度数据，根据该亮度数据可以计算该直线在图像像素行上的中心点坐标。将面阵探测器的光阑图像行序号记为m，I 为光阑图像的像素总行数，从下向上依次由小到大编号；将光阑图像列序号记为η，1 < η < N，N为光阑图像的像素总列数，从左向右依次由小到大编号。计算第I至第M行光阑图像的各直线中心点坐标，具体方法为5. 4. I 令 m=l ；5. 4. 2读取第m行光阑图像数据，根据各像素的亮度大小和位置计算所有非零区段；根据各非零区段长度判断区段的有效性，若非零区段的像素长度大于5且小于图像行总长的1/4则该非零区段有效，否则该非零区段无效；5. 4. 3计算第m行光阑图像的有效非零区段总数Km。建立面阵探测器坐标系，坐标系原点O为面阵探测器中心，X轴和Z轴分别沿面阵探测器的水平中心线和垂直中心线。将第m行、η列图像像素的灰度记为I (m，η)。计算第m行光阑图像的第I至第Km个有效非零区段的中心坐标，具体方法为5. 4. 3. I 令 km=l ；5. 4. 3. 2采用行亮度重心方法或行曲线拟合方法计算第Ii111个有效非零区段的中心坐标；行亮度重心方法为按式(一)和式(二)分别计算光阑图像第m行图像第 <个有效非零区段的中心水平坐标Xmkm和垂直坐标Zmkm，其中nkms为第km个有效非零区段的起始图像列序号，η-为第 <个有效非零区段的结束图像列序号。由于km作为Xm的下标后km的下标中的“m”特别小，导致显示不清楚，故在用km作为下标时均写成“km”，如“Xmkm”所示，以下所有以km作为下标的标识都采用这种办法。
Σ ",(爪，〃)
权利要求
1.一种基于准直光路的角秒级三维变形测量装置，其特征在于三维变形测量装置由发射�？�(33)、接收模块(34)和处理�？�(35)三大部分构成，发射�？�(33)与参考基准(23)固连，接收�？�(34)与待测物(22)固连，发射�？�(33)与接收模块(34)之间能够通视，发射�？�(33)发射的光束直接进入接收�？�(34)； 发射模块(33)由光源(I)、光阑⑵和发射光学系统(31)构成,光源(I)、光阑⑵和发射光学系统(31)固定在参考基准(23)上；光阑(2)位于光源(I)和发射光学系统(31)之间，且位于发射光学系统(31)的焦平面上,光源(I)和光阑(2)的中心都在发射光学系统(31)的光轴上，由此构成发射准直光路；光源(I)照射光阑(2)，发射光学系统(31)将透过光阑(2)的光束准直发射到接收�？�(34)中；光源(I)是半导体激光器或LED或灯泡；发射光学系统(31)采用双胶合透镜或双分离透镜或多片透镜组；光阑(2)是在玻璃等透明材料上利用光刻工艺加工而成的，加工后的光阑(2)有些部分透光，有些部分不透光，光阑(2)的透光部分采用L。条直线结构，各直线呈放射状分布，直线结构为透光部分，其余部分不透光，Lc为大于I的常数；每条直线由两个线段构成，直线越多，测量精度越高，但横滚方向的形变测量范围会减�。� 接收�？�(34)由接收光学系统(32)和面阵探测器(8)构成，接收光学系统(32)和面阵探测器(8)固定在待测物(22)上；面阵探测器(8)位于接收光学系统(32)的焦平面上，且其中心位于接收光学系统(32)的光轴上，由此构成接收准直光路；接收光学系统(32)采用双胶合透镜或双分离透镜或多片透镜组，面阵探测器(8)采用灰度或彩色的面阵CCD或CMOS探测器；发射�？�(33)发射的准直光束经过接收光学系统(32)会聚后成像在面阵探测器⑶上，面阵探测器⑶探测光阑(2)的直线结构光阑图像； 处理�？�(35)采用工控机、嵌入式计算机、台式计算机或笔记本电脑，在处理�？�(35)中安装了变形解算软件，处理�？�(35)通过千兆以太网、高速USB或1394接口与面阵探测器(8)连接，实时采集面阵探测器(8)探测的光阑图像，变形解算软件计算参考基准(23)与待测物(22)之间的三维变形，包括方位角、俯仰角和横滚角变形。
2.一种采用如权利要求I所述基于准直光路的角秒级三维变形测量装置进行三维变形测量的方法，其特征在于包括以下步骤 第一步，发射�？�(33)的光源(I)照射光阑(2),发射光学系统(31)将透过光阑(2)的光束准直发射到接收�？�(34)； 第二步，接收�？�(34)的接收光学系统(32)接收发射�？榉⑸涞墓馐⒔浠峋鄣矫嬲筇讲馄�(8)上； 第三步，面阵探测器(8)将接收的光束转换成光阑图像，此光阑图像为数字化的彩色图像或灰度图像； 第四步，处理�？�(35)从面阵探测器(8)接收光阑图像； 第五步，处理�？�(35)中变形解算软件根据光阑图像按以下步骤计算三维变形 .5. I若面阵探测器(8)是彩色类型探测器，则将彩色图像转换成灰度图像；若面阵探测器(8)是灰度类型探测器，则直接进入步骤5. 2 ； . 5. 2采用传统标定或自标定方法对发射光学系统畸变、接收光学系统畸变和面阵探测器畸变引起的光阑图像畸变进行校正； . 5. 3对校正后的光阑图像进行亮度分析，并滤除背景噪声；` 5. 4将面阵探测器(8)的光阑图像行序号记为m，I < m < M，M为光阑图像的像素总行数，从下向上依次由小到大编号；将光阑图像列序号记为n，I < η < Ν，Ν为光阑图像的像素总列数，从左向右依次由小到大编号，计算第I至第M行光阑图像的各直线中心点坐标，具体方法为` 5. 4. I 令 m=l ； ` 5. 4. 2读取第m行光阑图像数据，根据各像素的亮度大小和位置计算所有非零区段；根据各非零区段长度判断区段的有效性，若非零区段的像素长度大于5且小于图像行总长的1/4则该非零区段有效，否则该非零区段无效； ` 5. 4. 3统计第m行光阑图像的有效非零区段总数Km ;建立面阵探测器坐标系，坐标系原点O为面阵探测器(8)中心，X轴和Z轴分别沿面阵探测器的水平中心线和垂直中心线，将第m行、η列图像像素的灰度记为I (m，η)，计算第m行光阑图像的第I至第Km个有效非零区段的中心坐标，具体方法为 ` 5. 4. 3. I 令 Iini= I ； `5. 4. 3. 2采用行亮度重心方法计算第Iini个有效非零区段的中心坐标，行亮度重心方法为按式(一)和式(二)分别计算光阑图像第m行图像第km个有效非零区段的中心水平坐标Xmkm和垂直坐标Znikni,其中ntas为第km个有效非零区段的起始图像列序号，nkme为第Ii111个有效非零区段的结束图像列序号
3.如权利要求2所述的采用如权利要求I所述基于准直光路的角秒级三维变形测量装置进行三维变形测量的方法，其特征在于5. 3步对校正后的光阑图像进行亮度分析，并滤除背景噪声的方法是对光阑图像的所有像素根据亮度进行排序，根据光阑(2)的直线结构的直线数量L。选取L3个最亮像素，根据选择的最亮像素计算平均亮度，将该平均亮度作为光阑图像的最大亮度Imax ;根据Imax和面阵探测器噪声大小设置亮度阈值Ith，根据Ith滤除背景噪声，将小于Ith的所有像素的亮度置为零，Ith大于面阵探测器噪声且小于Imax的20%, L3 为 Lc 的 10 100 倍。
4.如权利要求2所述的采用如权利要求I所述基于准直光路的角秒级三维变形测量装置进行三维变形测量的方法，其特征在于5. 4. 3. 2步采用行曲线拟合方法计算第km个有效非零区段的中心坐标，行曲线拟合方法为按式(三)对光阑图像第m行图像第< 个有效非零区段的亮度进行P阶多项式拟合，其中Jkmp为P阶多项式的第P个系数，I < P < P，阶数P满足使拟合曲线计算的亮度与第km个有效非零区段的亮度的相对误差小于5% ；
5.如权利要求2所述的采用如权利要求I所述基于准直光路的角秒级三维变形测量装置进行三维变形测量的方法，其特征在于5. 8. 3. 2步采用列曲线拟合方法计算第kn个有效非零区段的中心坐标，列曲线拟合方法为按式(十二)对光阑图像第η列图像第< 个有效非零区段的亮度进行P阶多项式拟合，其中Jknp为P阶多项式的第P个系数，I阶数P满足使拟合曲线计算的亮度与第kn个有效非零区段的亮度的相对误差小于5% ；
6.如权利要求4或5所述的采用如权利要求I所述基于准直光路的角秒级三维变形测量装置进行三维变形测量的方法，其特征在于若第km次使用行亮度重心方法，则第km+l次也使用行亮度重心方法，若第km次使用行曲线拟合方法，则第km+l次也使用行曲线拟合方法;若第kn次使用列亮度重心方法，则第kn+l次也使用列亮度重心方法，若第kn次使用列曲线拟合方法，则第kn+l次也使用 列曲线拟合方法。
全文摘要
本发明公开了一种基于准直光路的角秒级三维变形测量装置及方法。装置由发射模块、接收�？楹痛砟？楣钩桑⑸淠？橛晒庠�、光阑和发射光学系统构成，光源和光阑的中心都在发射光学系统的光轴上；光阑的透光部分采用直线结构；接收�？橛山邮展庋低澈兔嬲筇讲馄鞴钩桑嬲筇讲馄魑挥诮邮展庋低车慕蛊矫嫔希渲行奈挥诮邮展庋低车墓庵嵘希淮砟？橹邪沧傲吮湫谓馑闳砑�。测量方法是发射光学系统将透过光阑的光束准直发射到接收模块，接收光学系统将光束会聚到面阵探测器上，面阵探测器将光束转换成光阑图像，变形解算软件根据光阑图像计算三维变形。采用本发明能同时测量方位角、俯仰角和横滚角变形，精度高，能实现角秒级三维变形测量。
文档编号G01B11/16GK102901459SQ20121037768
公开日2013年1月30日 申请日期2012年10月8日 优先权日2012年10月8日
发明者胡春生, 王省书, 秦石乔, 周金鹏, 高旸 申请人:中国人民解放军国防科学技术大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
