CN105825548B - 使用近心镜头的双棱镜单相机三维数字图像相关重构方法 - Google Patents

使用近心镜头的双棱镜单相机三维数字图像相关重构方法 Download PDF

Info

Publication number
CN105825548B
CN105825548B CN201610150858.5A CN201610150858A CN105825548B CN 105825548 B CN105825548 B CN 105825548B CN 201610150858 A CN201610150858 A CN 201610150858A CN 105825548 B CN105825548 B CN 105825548B
Authority
CN
China
Prior art keywords
biprism
camera
image
coordinate
point
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.)
Active
Application number
CN201610150858.5A
Other languages
English (en)
Other versions
CN105825548A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201610150858.5A priority Critical patent/CN105825548B/zh
Publication of CN105825548A publication Critical patent/CN105825548A/zh
Application granted granted Critical
Publication of CN105825548B publication Critical patent/CN105825548B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10004—Still image; Photographic image

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及使用近心镜头的双棱镜单相机三维数字图像相关重构方法，属于光测实验力学、三维数字图像相关技术领域，该方法包括：采集试样变形前、后的图像、棱边和别图像和标定板图像，标定获得相机内部、外部参数；建立平行于双棱镜后表面的参考世界坐标系,根据外部参数确定其与相机坐标系之间的转换关系；确定双棱镜棱边和原点标记的空间位置，建立空间坐标系；通过对标定板角点空间坐标的重构运算与试算确定空间坐标系的位置；根据试样图像中配对的像素坐标，分别对试样变形前、后的形貌进行重构运算，确定被测试样变形产生的位移。本方法大大降低了重构方法的模型误差。同时，本方法运算量小、编程方便、容易实现大批量数据的自动化处理。

Description

使用近心镜头的双棱镜单相机三维数字图像相关重构方法
技术领域
本发明属于光测实验力学、三维数字图像相关技术领域，特别涉及一种使用近心镜头的双棱镜单相机三维数字图像相关(BSL 3D DIC)重构方法。
背景技术
三维数字图像相关(3D DIC)方法是光测实验力学领域中常用的全场非接触变形测量方法，被广泛地应用于航天航空材料的性能测试和结构表征中。该方法使用多个相机同时观察被测区域，通过立体视觉原理以对物体的三维形貌进行重构，并同时获得其三维位移分量，具有适用范围广泛、测试结果可靠的优势。然而，传统的3D DIC方法仍然存在一些局限，例如：由多个相机构成测试系统需要占据较大的空间，而且各个相机的性能通常存在差异；实验时系统进行标定较为复杂；难以保证各个相机同步采集图像等等。这些因素都可能对测量结果产生影响。
基于双棱镜的单相机三维数字图像相关(BSL 3D DIC)方法是一种特殊的3D DIC方法。该方法通过在相机前放置双棱镜进行分光，仅通过单个相机即可完成上述的三维测量，有效地克服了传统多相机3D DIC方法的局限，是3D DIC技术的重要发展。作为BSL 3DDIC方法的核心技术，其重构方法的精确程度直接影响着BSL 3D DIC测量结果的准确性。
一个典型的BSL 3D DIC测试系统由图像采集设备、镜头和一个双棱镜构成，当使用近心镜头时，图像采集设备和镜头的组合可以等效为一个针孔成像模型。在一些研究中(如L.F.Wu,et al.,MeasSciTechnol,2014；L.F.Wu,et al.,Appl Opt,2015)，为了便于坐标描述，在图像采集设备的感光元件上定义像素坐标系PCS、在双棱镜上定义空间坐标系SCS；此外，考虑到双棱镜的两侧分别可看作楔形棱镜，Wu等人还在楔形棱镜上定义了局部坐标系LCS，以便于光线折射的分析。在测试之前需要在试样表面制备散斑。在测试过程中，将试样放置在BSL 3D DIC测试系统的视场之内，并依次采集试样变形前的图像和试样变形后的图像。通过BSL 3D DIC测试系统采集的图像，其左、由两侧分别为试样通过双棱镜的两侧后形成的图像，可以看作两个子图像。通常，将双棱镜中形成左侧子图像的一侧规定为正半侧，将形成右侧子图像的一侧规定为负半侧。通过数字图像相关(DIC)运算，可在左侧子图像与右侧子图像之间进行像素匹配。在测试结束后，将双棱镜从测试系统中移除，并将一块标定板放置于相机视场中，通过采集多幅标定板图像以实现对上述针孔成像模型内部参数和畸变参数的标定；同时，利用上述标定结果可以建立以相机光心为原点的相机坐标系CCS。对于每一幅标定图像，都可在上述在相机坐标系CCS中计算描述图像中标定板空间位置的外部参数。
目前，在计算机视觉领域关于棱镜单相机立体视觉技术的研究中，已经提出了多种可用于物体景深检测和三维形貌重构的模型，例如：虚像点模型(D.H.Lee,et al.,IEEETrans Rob Autom,2000)、虚相机模型(K.B.Lim,et al.,J Electron Imaging,2005)、几何光学模型(X.Y.Cui,et al.,J Opt Soc Am A，2012)、透视投影模型(X.Y.Cui,et al.,OptExpress,2015)等。近年来，实验力学领域的研究者提出并发展了若干新型的重构方法，如：Genovese等使用空间误差函数描述BSL 3D DIC系统中存在的系统误差，从而在不需要特定重构模型的前提下，实现了构件三维形貌和变形测量(K.Genovese,et al.,Opt LasersEng,2013)。但是，这种重构方法需要进行误差函数标定，实验操作相对复杂。Wu等人分析了双棱镜形成虚像点的规律，对传统的虚像点模型进行了改进，提高了其重构精度(L.F.Wu,et al.,MeasSciTechnol,2014)；同时，他们利用双远心镜头，发展了一种精度更高，且不需要标定的双远心BSL 3D DIC方法(L.F.Wu,et al.,Appl Opt,2015)。但是，由于双远心测试系统的视场范围、焦距、工作距离等不易调节，且成本较高，目前得到广泛研究和应用的多为使用近心镜头的BSL 3D DIC方法。
当BSL 3D DIC被用于精确测量材料变形产生的小位移(1000微应变以内)时，其重构方法要满足特定的条件。首先，重构方法应充分考虑测试系统中各个原件的对准误差及其各自产生的图像畸变，并尽可能予以消除，以获得更为精确的测量结果。此外，鉴于全场测量的数据点较多，重构方法应具有较高的效率，且容易实现自动化数据处理。然而，上述重构方法都存在一定的局限，例如：Wu等人提出的改进虚点模型中假设相机和双棱镜是对准的，而且忽略了光线在双棱镜中的三维折射效应，因而存在模型误差。为了满足材料变形产生的小位移测量的需求，BSL 3D DIC重构方法还有待进一步发展和完善。
发明内容
本发明为了进一步提高近心BSL 3D DIC方法的测量精度、增强其在材料变形产生的小位移测量和力学性能表征中的适用性，提出了一种使用近心镜头的双棱镜单相机三维数字图像相关重构方法。本方法大大降低了重构方法的模型误差。同时，本方法运算量小、编程方便、容易实现大批量数据的自动化处理。在本方法的基础上，使用BSL 3D DIC测试系统可以对构件发生小变形时产生的三维位移场进行精确测量。
本发明提出的一种使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其主要步骤包括：
S1：使用BSL 3D DIC测量系统采集试样变形前的图像IRef、试样变形后的图像IDef、棱边识别图像IB、标定板错位图像I0、标定板不同空间位置和姿态的直接观测图像I1-In，其中，图像I1中的标定板与双棱镜的后表面平行；使用图像I1-In对相机针孔成像模型参数进行标定，获得其焦距、主点像素坐标和畸变参数；根据与双棱镜后表面平行的标定板图像I1，计算相机外部参数，并建立参考世界坐标系RWCS；
S2：基于步骤S1计算得到的外部参数，建立相机坐标系CCS，确定相机坐标系CCS和上述参考世界坐标系RWCS的坐标转换关系；
S3：使用图像IB，确定双棱镜棱边的空间位置，计算双棱镜棱边在双棱镜后表面上的投影线的方程，并建立空间坐标系SCS；提取图像I0中角点对的像素坐标，使用步骤S1中获得的畸变参数，对角点对像素坐标进行镜头畸变误差修正后，通过重构运算计算对应角点的空间坐标；通过拟合与试算确定空间坐标系SCS与参考世界坐标系RWCS的间距Z0，及其转换关系；
S4：通过DIC运算，分别在试样变形前的图像IRef、试样变形后的图像IDef的左侧子图像和右侧子图像之间进行像素匹配，获得若干像点对像素坐标；使用步骤S1中的畸变参数，对像点对像素坐标进行镜头畸变误差修正，得到修正后的像点对像素坐标；
S5：基于修正后的像点对像素坐标，通过重构运算计算试样表面若干物点的空间坐标；分别对试样变形前、后的形貌进行重构运算后，根据对应物点的三维坐标增量，确定被测试样发生变形时产生的位移。
其中，步骤S3中所述的重构运算，其具体步骤包括：
RS1：对于镜头畸变误差修正后的角点对的像素坐标，令其坐标值为(x+，y+)和(x-，y-)，在所述相机坐标系CCS中分别计算其成像光线的单位方向矢量和以及位于双棱镜后表面上的出射点和的坐标；
RS2：在空间坐标系SCS中，计算所述出射点和处对应的棱镜厚度t+和+-；对每个出射点，分别建立局部坐标系LCS，确定其与在空间坐标系SCS之间的坐标变换关系；通过坐标变换，将所述和转换到各自的局部坐标系LCS中，分别记为和
RS3：根据t+和t-以及和分别在所述局部坐标系LCS中计算入射光线的单位方向矢量和以及入射点和的坐标；
RS4：通过坐标变换，将所述和以及所述入射点和的坐标转换到空间坐标系SCS中，分别记为和以及A+和A-；分别构建入射光线的空间方程，并联立求解获得对应物点或对应角点的空间坐标。
本发明的特点及有益效果：
本方法在现有BSL 3D DIC测量系统及其坐标系定义的基础上，在测量系统前方增设了一个可以平移的标定板，并在平行于双棱镜后表面的标定板上定义了参考世界坐标系RWCS；通过相机标定，确定了参考世界坐标系RWCS与相机坐标系CCS、空间坐标系SCS和局部坐标系LSC的空间位置关系，从而在重构过程中考虑了相机与双棱镜之间存在的对准误差。另外，本方法利用BSL 3D DIC测量系统成像光路可逆的特性，通过三维空间中的反向光线追迹，考虑了光线在双棱镜中的三维折射效应，对被测构件表面物点的空间坐标进行准确求解，大大降低了重构方法的模型误差。同时，本方法运算量小、编程方便、容易实现大批量数据的自动化处理。在本方法的基础上，使用BSL 3D DIC测试系统可以对构件发生小变形时产生的三维位移场进行精确测量。
附图说明
图1为典型的BSL 3D DIC测量系统的构成示意图；
图2为本发明使用近心镜头的双棱镜单相机三维数字图像相关重构方法流程图；
图3为本发明的实施例中使用BSL 3D DIC测量系统采集的棱边识别图像IB图片；
图4为计算双棱镜棱边上一点在其后表面上投影点位置的方法示意图；
图5为在局部坐标系LCS中确定入射光线的点向式空间方程的方法示意图；
具体实施方式
以下将通过附图及实施例详细描述本发明提出的使用近心镜头的双棱镜单相机三维数字图像相关重构方法。本实施例为一个承载拉伸试样轴向变形的测量。下面通过附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
本发明方法的实施例中基于典型的BSL 3D DIC测量系统，如图1所示。该系统包括：标定板1、双棱镜2、图像采集设备和近心镜头组合称为相机。其中：标定板1的空间位置和姿态可以任意调整；本实施例中使用的标定板为一块方格边长1.6mm的棋盘格标定板。相机用针孔成像模型3等效描述。相机光轴4与双棱镜2的后表面π之间存在对准误差。双棱镜2的最大厚度为t0，双棱镜两侧楔角分别表示为α+和α-；本实施例中变量上角标“+”表示该变量对应于测试系统正半侧光路，上角标“-”表示该参量对应于测试系统负半侧光路；本实施例中t0为11.52mm，楔角α+和α-均为19°50′。在双棱镜2的前表面棱边L上设置了棱边标记线和原点标记投影线MN为双棱镜2的前表面棱边L在其后表面π上的投影。OW-XWYWZW、O-XYZ、OC-XCYCZC和分别代表在BSL 3D DIC测量系统中定义的参考世界坐标系RWCS、空间坐标系SCS、相机坐标系CCS和局部坐标系LCS，其中OC为相机的光心。Γ为针孔成像模型3中的像平面；λ和(cx，cy)分别表示针孔成像模型3的焦距和主点像素坐标。折线表示了某一标定板角点(或空间物点)P通过BSL 3D DIC测量系统中双棱镜的正半侧成像时的光路，(x+，y+)为上述标定板角点(或空间物点)P在像平面Γ上形成的像点的像素坐标。
基于上述BSL 3D DIC测量系统，本实施例的使用近心镜头的双棱镜单相机三维数字图像相关重构方法流程如图2所示。包括以下步骤：
S1：使用BSL 3D DIC测量系统采集试样变形前的图像IRef、试样变形后的图像IDef、棱边识别图像It、标定板1错位图像I0、标定板1不同空间位置和姿态的直接观测图像I1-In(，n为正整数)；其中，图像I1中的标定板1与双棱镜2的后表面π平行；使用图像I1-In对相机针孔成像模型3的参数进行标定，获得其焦距λ、主点像素坐标(cx，cy)和畸变参数；根据与双棱镜2的后表面π平行时的标定板1的图像I1，计算相机外部参数，并建立参考世界坐标系RWCS；
结合图1，具体实现为：在采集图像IB时，将一个浅色物体置于所述BSL 3D DIC测量系统视场区域内，以便在图像中清晰地显示出棱边标记线和原点标记如图3所示，图3中深色线条为棱边标记线的图像，深色圆点为原点标记的图像，以棱边标记线为边界，采集的整体图像划分为左侧子图像和右侧子图像两个部分。采集图像I0时，双棱镜未从系统中被移除，且标定板1与双棱镜2的后表面π平行；图像I1-In是将双棱镜2移除后，通过相机直接采集的标定板1的图像；其中采集图像I1时，标定板1与双棱镜2的后表面π仍保持平行；而采集图像I2-In前，需要分别改变标定板1的位置和姿态。本实施例中共采集了12幅标定板直接观测图像，故n为12；图像IRef、图像IDef、图像IB以及图像I0的左侧子图像和右侧子图像分别由所述BSL 3D DIC测试系统的正半侧光路负半侧光路形成；本实施例中获得的针孔成像模型3的焦距λ为12375.03，(cx，cy)为(612.32,344.08)，畸变参数包括相机径向、切向畸变参数；参考世界坐标系RWCS由与双棱镜2后表面π平行时标定板1上的角点阵列确定，其原点OW为标定板上的一个角点，OWXW、OWYW坐标轴分别沿着角点阵列行、列的方向，OWZW垂直于平行标定板平面向外；
S2：基于步骤S1计算得到的外部参数，建立相机坐标系CCS，确定相机坐标系CCS和参考世界坐标系RWCS的坐标转换关系；
结合图1所示，具体实现为：外部参数包括旋转矩阵RCW和平移向量TC；所述相机坐标系CCS的原点位于相机光心OC，OCXC、OCYC坐标轴分别沿着像平面Γ上像素阵列行、列的方向，OCZC坐标轴沿着背离相机光心OC的方向；相机坐标系CCS和参考世界坐标系RWCS的坐标转换关系为：
vC＝RCWvW, (2.a)
其中，X代表一个空间坐标，v代表一个空间向量，各变量的下角标“C”表示该变量是在相机坐标系CCS中描述的；各变量的下角标“W”表示该变量是在RWCS中描述的；本实施例中：
TC＝[-7.7357 112665 384.6421]T。
S3：使用图像IB，确定棱边L的空间位置，计算棱边L在双棱镜2后表面π上的投影线MN的方程，并建立空间坐标系SCS；提取图像I0中角点对的像素坐标，使用步骤S1中获得的畸变参数，对角点对像素坐标进行镜头畸变误差修正后，通过重构运算计算对应角点的空间坐标；通过拟合与试算确定空间坐标系SCS与参考世界坐标系RWCS的间距Z0，及其转换关系；
以下结合图4对步骤S3的具体实现方法进行说明。图4中双棱镜2及其棱边L和后表面π、参考世界坐标系RWCS、空间坐标系SCS、棱边标记线原点标记相机光心OC、双棱镜2的最大厚度为t0以及投影线MN的定义均与图1相同。折线EBOC表示棱边L上的一点E透过双棱镜2后被相机接收并成像时的光路。
S31：对图像IB(图3)中棱边标记线中心线上的点以及原点标记的中心点进行识别，并进行镜头畸变修正；将得到图像IB各像素坐标记为(xB，yB)。；
S32：给定间距Z0的估测值，对于像素坐标(xB，yB)，在相机坐标系CCS中计算其成像光线的单位方向矢量nBC＝{nBCx，nBCy，nBCz}T以及该成像光线在双棱镜2后表面π上的出射点坐标B(XBx，YBC，ZBC)；其中nBCx、nBCy、nBCz通过下列方程组解得：
XBC、YBC、ZBC通过下列方程组解得：
r13(XBC-TCx)+r23(XBC-TCy)+r33(XBC-TCz)＝Z0； (4.b)
式中，r13、r23、r33、TCx、TCy和TCy是外部参数RCW和TC中的元素。通过(1)、(2)式可以将所述单位方向矢量nBC和坐标(XBC，YBC，ZBC)转换至参考世界坐标系RWCS中，分别记为nBW＝{nBWx nBWy nBWz)T和(XBW，YBW，Z0)；
S33：通过式(5.a)、(5.b)计算与出射点B对应的棱边L上的一点E在双棱镜2后表面π上的投影点世界坐标D(XDW，YDW，Z0)
式中，修正量δ可由双棱镜2的最大厚度为t0、双棱镜折射率N以及对应光线在B处发生折射时的入射角i和折射角r计算得到：δ＝t0tan i、本实施例中N取值为1.4662；
S34：根据S31中得到的各像素的坐标，通过步骤S32～S33中的计算，获得棱边L上的若干点在双棱镜2后表面π上的投影点坐标；利用所述投影点坐标进行线性拟合，获得投影线MN方程：
结合原点标记投影点的坐标以及步骤S32中给定的间距Z0的估测值，建立一个估测的空间坐标系SCS。如图4所示，估测的空间坐标系SCS的原点O为原点标记在双棱镜2后表面π上的投影点，OY坐标轴与棱边L在表面π上的投影线相重合，OZ坐标轴垂直于表面π，OX坐标轴与所述OY、OZ坐标轴构成右手坐直角标系。此外，所述估测的空间坐标系SCS的OXY平面与所述参考世界坐标系RWCS的OWXWYW平面相互平行；进而，所述间距Z0，是上述平行平面OXY与OWXWYW的间距。估测的空间坐标系SCS与参考世界坐标系RWCS的转换关系为：
vS＝RSWvW, (8.a)
式中，某一变量的下角标“S”表示该变量是在空间坐标系SCS中描述的，下同；sinθ和cosθ的值可由(6)式中的tan6确定，进而得到旋转矩阵RSW的值；平移向量TW(O)的值可由所述原点标记在RWCS中的坐标确定；
S35：提取I0中角点对的像素坐标，使用步骤S1中计算得到的畸变参数进行镜头畸变误差修正后，通过重构运算计算对应角点的空间坐标。
所述重构运算具体包括以下步骤：
RS1：对于镜头畸变误差修正后的角点对的像素坐标，令其坐标值为(x+,y+)和(x-,y-)，在所述相机坐标系CCS中分别计算其成像光线的单位方向矢量和以及位于双棱镜2后表面π上的出射点和的坐标；
以图1中所示的像素坐标(x+,y+)为例，和的计算通过(3)、(4)式完成，将式中的变量xB、yB分别替换为x+、y+即可。
同理可由(x-,y-)计算得到和
RS2：在空间坐标系SCS中，计算所述出射点和处对应的棱镜厚度t+和t-；对于每个出射点，分别建立局部坐标系LCS，确定其与在空间坐标系SCS之间的坐标变换关系；通过坐标变换，将所述和转换到各自的局部坐标系LCS中，分别记为和
以图1中出射点和步骤RS1中获得的为例，步骤RS2具体实现为：
RS21：通过(1)、(7)式，将RS1中计算得到的坐标由相机坐标系CCS转换至空间坐标系SCS，记为相应地，t+的值可由双棱镜2的最大厚度t0与其正半侧楔角α+计算得到：
RS22：以为原点建立局部坐标系LCS，坐标轴的方向分别与空间坐标系SCS的OY、OZ坐标轴方向相同，坐标轴平行于空间坐标系SCS的OX坐标轴，且指向点所对应的一侧双棱镜厚度减小的方向。则从局部坐标系LCS到世界坐标系SCS的坐标变换关系如下：
式中，某一变量的上标“^”表示该变量是在局部坐标系LCS中描述的，下同；和矩阵Sign由在空间坐标系SCS中的坐标确定；
RS23：利用(2)、(8)、(11)式，将由相机坐标系CCS转换到局部坐标系LCS中，并计算得到
RS24：同理由和计算得到t-和
RS3：根据t+和t-和和分别在所述局部坐标系LCS中计算入射光线的单位方向矢量和以及入射点和的坐标；
以下结合图5，以步骤RS2中得到的t+和为例，对步骤RS3的具体实现方法进行说明。图5中坐标系是以为原点建立的局部坐标系LCS。在上述坐标系中，图1中的双棱镜2可局部地视为一个楔形棱镜，t+表示该楔形棱镜的最大厚度。为成像光线5的单位方向向量，为折射光线6的单位方向向量，为入射光线7的单位方向向量。楔角α+和双棱镜2后表面π的定义均与图1相同。
RS31：在局部坐标系LCS中，根据成像光线5的单位方向向量为：
以及所述楔形棱镜的最大厚度t+、楔角α+和双棱镜的折射率N，通过求解下列方程组，计算入射点的坐标
具体而言，由方程(12.a)解得，其正负性与相同；将代入方程(12.b)、(12.c)可分别解得
RS32：利用可计算得到折射光线6的单位方向矢量
RS33：假设入射光线7的单位方向矢量为其三个分量和可通过下列方程组求解得到：
具体而言，的值可以通过方程(14.a)直接获得；方程(14.c)是关于的一元二次方程，考虑到需为正数，故将其取为方程(14.c)的正根；将代入方程(14.b)，可以确定未知量的值；进而，可以确定假设量的值。
RS34：同理可由t-和计算得到和入射点的坐标。
RS4：通过坐标变换，将所述和以及所述入射点和的坐标转换到空间坐标系SCS中，分别记为和以及A+和分别构建入射光线的空间方程，并联立求解获得对应角点的空间坐标；
以步骤RS3中得到的入射点的坐标和为例，步骤RS4具体实现为：
RS41：利用(10)、(11)式将入射点的坐标和由局部坐标系LCS转换到空间坐标系SCS，记为和
RS42：利用和在空间坐标系SCS下构建图5中所示的入射光线7的点向式空间方程：
其中：(X,Y,Z)为入射光线7上任意一点的空间坐标。
RS43：同理可由入射点的坐标和计算得到和并构建对应于测试系统负半侧光路的入射光线点向式空间方程：
其中：(X,Y,Z)为对应于测试系统负半侧光路的入射光线上任意一点的空间坐标；
RS44：由于步骤RS42和步骤RS43中确定的入射光线是由同一角点物点发射出的，联立(15)、(16)式，即可对角点的空间坐标进行求解；
S36：将所有的角点空间坐标拟合为一个正交排列的、离面坐标均为Z0的空间点阵，并计算拟合残差；
S37：改变Z0的取值，并重复步骤S32～S36，通过多次试算后，选取使得拟合残差达到最小值时的Z0，作为空间坐标系SCS与参考世界坐标系RWCS的间距Z0的最优值。本实施例中间距Z0的最优值为72.065mm；
S38：利用步骤S37中确定的间距Z0的最优值，重复步骤S31～S34，最终确定空间坐标系SCS，并建立其与参考世界坐标系RWCS的转换关系。本实施例中，间距Z0取值为72.065mm时，得到投影线MN在参考世界坐标系RWCS中的方程为：
原点标记在参考世界坐标系RWCS中的坐标为：(8.529,7.017,72.065)。进而，
TW(O)＝[8.529 7.017 72.065]T。；
S4：通过DIC运算，分别在试样变形前的图像IRef、试样变形后的图像IDef的左侧子图像和右侧子图像之间进行像素匹配，获得若干像点对像素坐标；使用步骤S1中的畸变参数，对像点对像素坐标进行镜头畸变误差修正，得到修正后的像点对像素坐标；
S5：基于修正后的像点对像素坐标，通过重构运算计算对应物点的空间坐标；分别对试样变形前、后的形貌进行重构后，根据对应物点的三维坐标增量，确定被测试样发生变形时产生的位移。
具体而言：步骤S5中的“重构运算”与步骤S35中的“重构运算”基本相同。在进行物点的空间坐标重构时，只需要将上述“重构运算”步骤RS1中的坐标(x+，y+)和(x-,y-)替换为步骤S4中得到的修正后的像点对像素坐标即可。相应地，步骤S5中的“重构运算”可获得与像点对相对应的试样表面若干物点的空间坐标。
本实施例中试样拉伸产生的轴向变形为125微应变，其产生的最大轴向相对位移量不足5μm，尽管如此，采用本发明提出的方法，试样拉伸位移沿轴向线性分布的规律仍被BSL 3D DIC系统准确检测出。从而说明本发明提出的使用近心镜头的双棱镜单相机三维数字图像相关重构方法的可行性与可靠性。

Claims (5)

1.一种使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其特征在于，该方法主要步骤包括：
S1：使用BSL 3D DIC测量系统采集试样变形前的图像IRef、试样变形后的图像IDef、棱边识别图像IB、标定板错位图像I0、标定板不同空间位置和姿态的直接观测图像I1-In，其中，图像I1中的标定板与双棱镜的后表面平行；使用图像I1-In对相机针孔成像模型参数进行标定，获得其焦距、主点像素坐标和畸变参数；根据与双棱镜后表面平行的标定板图像I1，计算相机外部参数，并建立参考世界坐标系RWCS；
S2：基于步骤S1计算得到的外部参数，建立相机坐标系CCS，确定相机坐标系CCS和上述参考世界坐标系RWCS的坐标转换关系；
S3：使用图像IB，确定双棱镜棱边的空间位置，计算双棱镜棱边在双棱镜后表面上的投影线的方程，并建立空间坐标系SCS；提取图像I0中角点对的像素坐标，使用步骤S1中获得的畸变参数，对角点对像素坐标进行镜头畸变误差修正后，通过重构运算计算对应角点的空间坐标；通过拟合与试算确定空间坐标系SCS与参考世界坐标系RWCS的间距Z0，及空间坐标系SCS与参考世界坐标系RWCS的转换关系；
S4：通过DIC运算，分别在试样变形前的图像IRef、试样变形后的图像IDef的左侧子图像和右侧子图像之间进行像素匹配，获得若干像点对像素坐标；使用步骤S1中的畸变参数，对像点对像素坐标进行镜头畸变误差修正，得到修正后的像点对像素坐标；
S5：基于修正后的像点对像素坐标，通过重构运算计算试样表面若干物点的空间坐标；分别对试样变形前、后的形貌进行重构运算后，根据对应物点的三维坐标增量，确定被测试样发生变形时产生的位移；
所述重构运算，具体步骤包括：
RS1：对于镜头畸变误差修正后的角点对或像素点对的像素坐标，令其坐标值为(x+,y+)和(x-,y-)，在所述相机坐标系CCS中分别计算其成像光线的单位方向矢量和以及位于双棱镜后表面上的出射点和的坐标；
RS2：在空间坐标系SCS中，计算所述出射点和处对应的棱镜厚度t+和t-；对每个出射点，分别建立局部坐标系LCS，确定各局部坐标系LCS与在空间坐标系SCS之间的坐标变换关系；通过坐标变换，将所述和转换到各自的局部坐标系LCS中，分别记为和
RS3：根据t+和t-以及和分别在所述局部坐标系LCS中计算入射光线的单位方向矢量和以及入射点和的坐标；
RS4：通过坐标变换，将所述和以及所述入射点和的坐标转换到空间坐标系SCS中，分别记为和以及A+和A-；分别构建入射光线的空间方程，并联立求解获得对应角点或对应物点的空间坐标。
2.如权利要求1所述使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其特征在于，所述步骤S1中还包括将所述BSL 3D DIC测量系统中的双棱镜后表面面向相机，在采集棱边识别图像IB时，将一个浅色物体置于所述BSL 3D DIC测量系统视场区域内，便于在图像中清晰地显示出棱边标记线和原点标记。
3.如权利要求1所述使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其特征在于，所述步骤S1中在采集标定板错位图像I0时，双棱镜置于系统中，且所述标定板与所述双棱镜后表面平行；在采集标定板直接观测图像I1-In是时，将所述双棱镜移除后，通过相机直接采集的所述标定板的图像；其中采集图像I1时，所述标定板与双棱镜后表面保持平行。
4.如权利要求1所述使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其特征在于，所述步骤S2中所述相机坐标系CCS的原点OC位于相机光心，OCXC、OCYC坐标轴分别沿着像素阵列行、列的方向，OCZC坐标轴沿着背离相机光心的方向；所述参考世界坐标系RWCS是由位于所述平行于双棱镜后表面的标定板上的角点阵列确定的坐标系，其原点OW为标定板上的一个角点，OWXW、OWYW坐标轴分别沿着角点阵列行、列的方向，OWZW垂直于平行标定板平面向外；所述参考世界坐标系RWCS的作用在于建立相机坐标系CCS与空间坐标系SCS的转换关系，以减小相机与双棱镜之间存在的对准误差。
5.如权利要求2所述使用近心镜头的双棱镜单相机三维数字图像相关重构方法，其特征在于，所述步骤S3中所述空间坐标系SCS的原点O为所述原点标记在双棱镜后表面上的投影点，OY坐标轴与所述棱边标记线在双棱镜后表面上的投影线相重合，OZ坐标轴垂直于双棱镜后表面，OX坐标轴与所述OY、OZ坐标轴构成右手坐直角标系；所述空间坐标系SCS的OXY平面与所述参考世界坐标系RWCS的OWXWYW平面相互平行；进而所述间距Z0是上述平行平面OXY与OWXWYW的间距。
CN201610150858.5A 2016-03-16 2016-03-16 使用近心镜头的双棱镜单相机三维数字图像相关重构方法 Active CN105825548B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610150858.5A CN105825548B (zh) 2016-03-16 2016-03-16 使用近心镜头的双棱镜单相机三维数字图像相关重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610150858.5A CN105825548B (zh) 2016-03-16 2016-03-16 使用近心镜头的双棱镜单相机三维数字图像相关重构方法

Publications (2)

Publication Number Publication Date
CN105825548A CN105825548A (zh) 2016-08-03
CN105825548B true CN105825548B (zh) 2018-08-10

Family

ID=56523866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610150858.5A Active CN105825548B (zh) 2016-03-16 2016-03-16 使用近心镜头的双棱镜单相机三维数字图像相关重构方法

Country Status (1)

Country Link
CN (1) CN105825548B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109819235A (zh) * 2018-12-18 2019-05-28 同济大学 一种具备跟踪功能的轴向分布感知集成成像方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106441135A (zh) * 2016-08-29 2017-02-22 清华大学 高温环境下单相机三维变形与温度同步测量装置及方法
CN107121333B (zh) * 2017-04-24 2019-04-30 苏州汇才土水工程科技有限公司 一种整体-局部相结合的试样变形测量装置及方法
CN109781141B (zh) * 2018-12-29 2020-10-27 深圳航星光网空间技术有限公司 立方棱镜组的标定方法及装置
CN112797923B (zh) * 2021-01-05 2022-12-30 上海交通大学 校正粒子衍射图像图案中心和欧拉角的方法、系统及介质
CN112904766B (zh) * 2021-01-18 2022-06-21 福州大学 基于像素误差反馈的旋转双棱镜指向控制系统与方法
CN113703490B (zh) * 2021-08-30 2023-06-23 福州大学 一种基于粒子群算法的旋转双棱镜指向校正方法
CN114543667B (zh) * 2022-01-21 2023-11-07 同济大学 一种基于神经网络的单相机-双棱镜三维测量系统及测量方法
CN115063567B (zh) * 2022-08-19 2022-10-28 中国石油大学(华东) 一种双棱镜单目立体视觉系统三维光路分析方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072877A (zh) * 2010-12-13 2011-05-25 清华大学 一种实时残余应力测量系统及方法
TW201140494A (en) * 2010-05-03 2011-11-16 Nat Univ Tsing Hua Calibration method of three dimensional digital image correlation (3D-DIC)
CN102589476A (zh) * 2012-02-13 2012-07-18 天津大学 高速扫描整体成像三维测量方法
US8339603B1 (en) * 2008-10-03 2012-12-25 J.A. Woollam Co., Inc. Mapping ellipsometers and polarimeters comprising polarization state compensating beam directing means, and method of use
CN103394973A (zh) * 2013-08-06 2013-11-20 山东大学 数字图像相关法测量高速正交切削过程应变装置及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8339603B1 (en) * 2008-10-03 2012-12-25 J.A. Woollam Co., Inc. Mapping ellipsometers and polarimeters comprising polarization state compensating beam directing means, and method of use
TW201140494A (en) * 2010-05-03 2011-11-16 Nat Univ Tsing Hua Calibration method of three dimensional digital image correlation (3D-DIC)
CN102072877A (zh) * 2010-12-13 2011-05-25 清华大学 一种实时残余应力测量系统及方法
CN102589476A (zh) * 2012-02-13 2012-07-18 天津大学 高速扫描整体成像三维测量方法
CN103394973A (zh) * 2013-08-06 2013-11-20 山东大学 数字图像相关法测量高速正交切削过程应变装置及方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A modified virtual point model of the 3D DIC technique using a single camera and a bi-prism;Lifu Wu et al.;《Measurement Science and Technology》;20141013;第1-14页 *
Single-lens 3D digital image correlation system based on a bilateral telecentric lens and a bi-prism: Systematic error analysis and correction;Lifu Wu et al.;《Optics and Lasers in Engineering》;20160302;第87卷;第129-138页 *
Single-lens 3D digital image correlation system based on a bilateral telecentric lens and a bi-prism: validation and application;LIFU WU et al.;《Applied optics》;20150910;第54卷(第26期);第7842-7850页 *
高速三维数字图像相关系统及其动载三维变形测量;曾祥福 等;《北京理工大学学报》;20120430;第32卷(第4期);第364-369页 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109819235A (zh) * 2018-12-18 2019-05-28 同济大学 一种具备跟踪功能的轴向分布感知集成成像方法

Also Published As

Publication number Publication date
CN105825548A (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
CN105825548B (zh) 使用近心镜头的双棱镜单相机三维数字图像相关重构方法
Pan et al. Review of single-camera stereo-digital image correlation techniques for full-field 3D shape and deformation measurement
CN103559735B (zh) 一种三维重建方法及系统
CN110514143B (zh) 一种基于反射镜的条纹投影系统标定方法
CN108151671B (zh) 一种三维数字成像传感器、三维扫描系统及其扫描方法
CN106127745B (zh) 结构光3d视觉系统与线阵相机的联合标定方法及装置
JP6319329B2 (ja) プレノプティック・カメラを使った表面属性の推定
US9123116B2 (en) Multiview face capture using polarized spherical gradient illumination
CN102944191B (zh) 一种基于平面圆靶标的三维视觉测量数据拼接方法及装置
JP6580761B1 (ja) 偏光ステレオカメラによる深度取得装置及びその方法
CN106500629B (zh) 一种显微三维测量装置及系统
Wu et al. An accurate method for shape retrieval and displacement measurement using bi-prism-based single lens 3D digital image correlation
CN108629841A (zh) 一种基于激光散斑多视点三维数据测量方法及系统
CN109919856A (zh) 基于双目视觉的沥青路面构造深度检测方法
Luo et al. Design and validation of a demand-oriented single-camera stereo-DIC system with a four-mirror adapter
US10151581B2 (en) Method and device for determining the position and orientation of a specular surface forming a diopter
CN110728745B (zh) 一种基于多层折射图像模型的水下双目立体视觉三维重建方法
CN108175535A (zh) 一种基于微透镜阵列的牙科三维扫描仪
CN110044266A (zh) 基于散斑投影的摄影测量系统
CN109506562A (zh) 一种用于太阳翼展开锁定深度检测的双目视觉测量装置
Chai et al. Epipolar constraint of single-camera mirror binocular stereo vision systems
Wieneke et al. Lagrangian particle tracking in the presence of obstructing objects
CN112562007B (zh) 一种基于三目约束的包裹相位非展开快速立体匹配方法
CN110146032A (zh) 基于光场分布的合成孔径相机标定方法
Kadambi et al. Frequency domain TOF: encoding object depth in modulation frequency

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