Next Article in Journal
Analytical Modeling of the Temperature Using Uniform Moving Heat Source in Planar Induction Heating Process
Next Article in Special Issue
3D Face Profilometry Based on Galvanometer Scanner with Infrared Fringe Projection in High Speed
Previous Article in Journal
Estimation of Source-Based Aerosol Optical Properties for Polydisperse Aerosols from Receptor Models
Previous Article in Special Issue
Dynamic Phase Measuring Profilometry Based on Tricolor Binary Fringe Encoding Combined Time-Division Multiplexing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Calibration Method for System Parameters in Direct Phase Measuring Deflectometry

1
State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300130, China
2
School of Mechanical Engineering, Heibei University of Technology, Tianjin 300130, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(7), 1444; https://doi.org/10.3390/app9071444
Submission received: 13 March 2019 / Revised: 3 April 2019 / Accepted: 3 April 2019 / Published: 6 April 2019
(This article belongs to the Special Issue High-speed Optical 3D Shape and Deformation Measurement)

Abstract

:

Featured Application

The paper proposes a novel method of calibrating system parameters for a PMD (phase measuring deflectometry) based system. This calibration method can be used in the areas of reverse engineering, car industry and aerospace industry.

Abstract

Phase measuring deflectometry has been widely studied as a way of obtaining the three-dimensional shape of specular objects. Recently, a new direct phase measuring deflectometry technique has been developed to measure the three-dimensional shape of specular objects that have discontinuous and/or isolated surfaces. However, accurate calibration of the system parameters is an important step in direct phase measuring deflectometry. This paper proposes a new calibration method that uses phase information to obtain the system parameters. Phase data are used to accurately calibrate the relative orientation of two liquid crystal display screens in a camera coordinate system, by generating and displaying horizontal and vertical sinusoidal fringe patterns on the two screens. The results of the experiments with an artificial specular step and a concave mirror showed that the proposed calibration method can build a highly accurate relationship between the absolute phase map and the depth data.

1. Introduction

The past several decades have seen revolutionary developments in three-dimensional (3D) optical shape measurement techniques [1,2], which are now widely used across many different fields, including reverse engineering, 3D games, digitization of cultural relics, clothing, biometrics and other biological fields [3,4,5,6,7,8,9]. In the area of structured light techniques, phase-calculation-based fringe pattern projection has made rapid progress so far, because it has the advantages of full-field measurement, non-contact operation, automatic data processing, and fast data acquisition [10,11,12,13]. Nevertheless, there are a large number of specular objects that need to be measured [14]. Current fringe pattern projection techniques cannot deal with this demand, especially in the measurement of specular objects that have discontinuous and/or isolated surfaces.
With the increasing demands of industrial development, much relative research has been done to solve the problem of specular object measurement. Shape reconstruction from gradient data [15,16,17,18] based on phase measuring deflectometry (PMD, also called fringe reflection profilometry) [19,20,21,22,23,24,25,26,27,28] has been widely studied because of its high dynamic range, full-field measurement, fast data acquisition, and automatic data processing. Cross path integration techniques, Fourier transform integration, and area wavefront reconstruction methods have been demonstrated [29]. However, existing PMD methods have the drawback of local error accumulation and are incapable of measuring discontinuous specular objects while reconstructing 3D shapes by integrating gradient data. Many improved algorithms have been proposed to avoid this integration procedure. Huang et al. [30] presented an iterative method that uses discrete cosine transforms to deal with the integration problem, with an incomplete gradient dataset in a Southwell configuration. Yuan et al. [31] presented a method of measuring the absolute height and slope of a specular surface with an incident ray, camera pinhole, and dummy assistant surfaces. Song et al. [32] selected a binary shifting strip as a structured light pattern for 3D reconstruction of the specular surface. Recently, Huang et al. [33] have proposed a model phase measuring deflectometry technique that uses a mathematical method to reconstruct 3D shapes of specular objects.
In all these PMD methods, calibration of the system parameters plays a key role in the measurement of the specular surface because it directly influences the accuracy of the measurement system. Ren et al. [34] presented an iterative optimization algorithm to calibrate stereo deflectometry. It includes two steps: a stepwise method and a self-calibration method; however, errors accumulated via the iterative principle is a potential problem. Xiao et al. [35] introduced a calibration method that uses a markerless flat mirror. This method requires a two-dimensional cosine fringe pattern to be reflected three times. It is easy for the images to become out of focus, which influences the accuracy of extraction of the feature points on the screen. System calibration is thus still a challenging problem in measuring 3D shapes by PMD.
A novel DPMD (direct PMD) method that can obtain the full-field 3D shape of specular objects with discontinuous and/or isolated surfaces has been developed to solve the problem of measuring complex specular surfaces [36,37,38]. This method builds a direct relationship between the absolute phase and the depth information. The DPMD method can calculate depth data from one directional fringe patterns at two display screens, instead of two orthogonal fringe patterns. Moreover, 3D shape data of specular objects could be obtained from two unwrapped phase maps directly, without the procedure of integrating gradient data. Calibrating the system parameters is an important step in measuring 3D shape data accurately, including the distance between the reference plane and the liquid crystal device (LCD) screen, and the distance between the reference plane and the charge coupled device (CCD) camera. In order to calibrate the system parameters, the existing methods used a translating stage and a machine vision-based method. Therefore, the calibration accuracy, on one hand, depends on the accuracy of the stage, but on the other hand, it has a low accuracy due to the limited depth of field (DOF) of the lens.
This paper presents a novel parameter calibration method that uses phase data to accurately determine the two system parameters. The phase information is calculated from sinusoidal fringe patterns, which are unaffected by the imaging lens being out of focus, so the calibrated system parameters have high accuracy. The following section briefly describes the DPMD technique. Section 3 explains the principle of using phase data to calibrate the system parameters. Section 4 describes the experiments used to verify the proposed calibration method and the results. Section 5 presents our conclusions.

2. Direct Phase Measuring Deflectometry

Figure 1 is a schematic diagram of the proposed full-field 3D shape measurement system for specular surfaces using the DPMD technique. The system consists of two LCD screens, one mirror, one CCD camera, and one beam splitter plate.
The beam splitter was adjusted to a suitable position that makes a virtual image LCD1 of LCD1 parallel to LCD2, as if the two LCD screens were located in two parallel positions. The two screens were parallel to a reference mirror located in the reference plane. The parallelism was established by repeated adjustments of the relative orientations of the three components (the two LCD screens and reference mirror). A checkerboard could be generated and displayed on the LCD1 and LCD2 screens, which were both located on the angular displacement platforms. As a result, the screens could change position and direction randomly and accurately. A CCD camera captured the checkerboard images on the two screens. The extrinsic parameters of LCD1 and LCD2 were calculated from the captured checkerboard image. Then, based on their differences, the relative orientation between the virtual image of LCD1 and LCD2 was adjusted to be parallel, using the angular displacement platforms. This procedure was repeated several times until the three components were nearly parallel. And then there is also Zhang’s method [39] for compensation by generating deformed fringes by software to be equivalent to the exact parallel alignment of two display screens.
In the diagram, h is the depth of the measured point with respect to the reference mirror, d1 is the distance between LCD1 and the reference mirror, and d2 is the distance between LCD2 and the reference mirror.
As illustrated in Figure 1, two rays of light were displayed and reflected into the CCD camera via the measured surface and the mirror at the reference position. The two incident rays corresponded to the same reflection light. The absolute phases of the two incident rays are denoted as ϕ 1 and ϕ 2 on the reference plane and ϕ 1 and ϕ 2 on the measured specular surface. θ and 2 θ + β are the angle between the incident ray and the normal vector of the reference and the angle between the incident ray and the normal vector of the measured specular surface, respectively. The period of the displayed fringe pattern on the LCD screen is denoted as q. l is the distance on LCD1 between the two incident rays because of the height and gradient of the measured surface. Parameter h stands for the height of the measured specular surface with respect to the reference plane. d1 is the distance between LCD1 and the reference mirror, and d2 is the distance between LCD2 and the reference mirror.
According to geometric triangle relationships, which are illustrated in Figure 1, Equations (1)–(4) can be acquired respectively.
( d 1 d 2 ) tan β = ( ϕ 1 ϕ 2 ) × q 2 π
( d 1 d 2 ) tan ( 2 θ + β ) = ( ϕ 1 ϕ 2 ) × q 2 π
( d 1 + h ) tan θ + l = ( d 1 h ) tan ( 2 θ + β )
l = ( ϕ 1 ϕ 1 ) × q 2 π .
Equation (5) can be derived based on the above Equations (1)–(4).
h = d 1 ( ϕ 2 ϕ 2 ) d 2 ( ϕ 1 ϕ 1 ) ( ϕ 1 ϕ 2 ) + ( ϕ 1 ϕ 2 ) .
This equation shows clearly that the depth value can be calculated directly from the captured fringe patterns as long as the two parameters d1 and d2 and phase information on the reference mirror are known beforehand. The calibration method used to obtain the system parameters is explained in the following section.

3. Calibration of System Parameters

An important step in obtaining 3D shape data is to build the relationship between the absolute phase and the depth, which is known as calibration. In the mathematical model derived in Equation (5), the two parameters d1 and d2, need to be determined beforehand. When the extrinsic parameters (orientations) of the mirror, LCD1, and LCD2 have been acquired in the same camera coordinate system, d1 can be calibrated from the relative distance between the mirror and LCD1, which can be obtained through matrix transformation; d2 can be calibrated in the same way. A novel parameter calibration method for the DPMD system has been proposed by using a phase target to solve the existing problems of dependency on a translating stage and the limited depth of field (DOF) of the lens.

3.1. Calibration of Internal Parameters

To apply a machine vision-based method to determine distance d1 and d2 by a CCD camera, the internal parameters need to be calibrated by using the following equation:
λ [ u v 1 ] = A [ R T ] [ X W Y W Z W 1 ] ,
where R is a matrix representing the three rotation angles and T = [tx, ty, tz] is a vector representing the three linear translations. [Xw, Yw, Zw] is the coordinate vector of a point in the world coordinate system, while [u, v] is a coordinate vector of the corresponding point in the pixel coordinate system. λ is an arbitrary scaling factor. A is a matrix of the CCD camera internal parameters, including two focal lengths (Fu and Fv), two principal point coordinates (Pu and Pv), and four image radial and tangential distortion coefficients (K1, K2, K3, and K4). The eight internal parameters need to be calibrated beforehand and Zhang’s camera calibration method is used to obtain their values by using a checkerboard at different positions [40].

3.2. Orientation of the Mirror

A mirror with a matrix of discrete hollow ring markers on its surface was used to determine the distance between the reference plane and the two LCD screens in orientating the reference plane in the camera coordinate system. After the mirror plane was adjusted to be parallel to the screen surface, the CCD camera captured an image of the hollow ring marker matrix, as illustrated in Figure 2a.
After extracting the inner and outer edges of each ring marker, the center of each of the markers was determined using the principles of cross ratio and the pole–polar relationship [41]. Figure 2b shows the center positions determined for each marker in the matrix, illustrated by red dots. Because the distance between neighboring ring markers is known beforehand, the external parameters of the mirror, rotation matrix R, and translating vector T could be obtained in the camera coordinate system to define the orientation of the mirror [42].
Obviously, the location accuracy of each center position was the most influential factor in the above process, and had an effect on determining the orientation of the mirror. This paper used Chen’s [30] method to calculate the real projection center to reduce the error as much as possible, this method is more accurate and stable than using the ellipse fitting method, as the radius of inner circle increases.

3.3. Orientations of the Two Screens

The luminance nonlinearity caused by the gamma effect of the digital screen and camera is the major error source of digital phase shift technology [43]. As a result, gamma correction is necessary to reduce the influence for the following calibration process. This paper adopted a direct, simple, and effective method by encoding initial fringe pattern based on Wang’s paper [44]. Then, horizontal and vertical encoded sinusoidal fringe patterns with the optimum fringe numbers were subsequently generated and displayed on LCD1 and LCD2 to determine the orientations of the two LCD screens in the same camera coordinate system. However, the fringe patterns were reflected by the reference mirror and captured by the CCD camera from another position for post-processing. Because phase data were calculated from the captured sinusoidal fringe patterns, they were insensitive to the imaging lens being out of focus. The phase of each pixel on the two LCD screens could be known using certain fringe numbers.
Absolute phase data were calculated using the optimum three-fringe number selection method and a four-step phase shifting algorithm. However, the phase data in the location of the circles could not be obtained correctly due to the low fringe modulation, as illustrated in Figure 3a for the horizontal fringe pattern. Because all phase data were located on the plane mirror, a plane could be fitted by using the obtained correct phase data without circles. This process can eliminate the influence of the inaccurate phase data in the location of the circles, as shown in Figure 3b.
A two-dimensional feature point matrix can be designed for the LCD screen to match the calibration requirements, as illustrated in Figure 4. The world coordinate and corresponding phase data of each feature point can be determined by following equations:
x w ϕ x = n x × p 2 π ,
y w ϕ y = n y × p 2 π ,
where p is the size of LCD screen pixel pitch, nx is the number of LCD screen pixels per horizontal fringe, and ny is the number of LCD screen pixels per vertical fringe. According to Equations (7) and (8), one camera pixel can uniquely locate its corresponding physical position (xw, yw) in the world coordinate system based on its phase value (φx, φy).
As a result, the position of each feature point in the unwrapped phase map was determined using a phase matching algorithm. The absolute phase of every pixel was invariable in the image captured from different camera viewpoints. The phase data can thus be used to build the relationship of the feature points between the world coordinate system and camera coordinate system.
The orientations of the two LCD screens in the camera coordinate system were determined from the parameter matrix of the LCD screen and the calibrated internal parameters of the CCD camera.

3.4. Calculation of the two system parameters

In order to calculate the distance, the orientations of the reference mirror and the two screens LCD1 and LCD2 needed to be transformed into the same camera coordinate system. After obtaining the orientations of the reference mirror and screens LCD1 and LCD2 in the same camera coordinate system, the coordinate transformations for the LCD screens and the mirror were given by the following equations.
R M C × P M W + T M C = P M C ,
R L 1 C × P L 1 W + T L 1 C = P L 1 C ,
R L 2 C × P L 2 W + T L 2 C = P L 2 C ,
where, R M C , R L 1 C , and R L 2 C are the rotation matrices for the mirror, LCD1, and LCD2 from the world coordinate system to the camera coordinate system; P M W , P L 1 W , and P L 2 W are the world coordinates of the mirror, LCD1, and LCD2; T M C , T L 1 C , and T L 2 C are the translation vectors of the mirror, LCD1, and LCD2 from the world coordinate system to the camera coordinate system; and P M C , P L 1 C , and P L 2 C are the position coordinates of the mirror, LCD1, and LCD2 in the camera coordinate system.
The distances between the two LCD screens and the mirror were calculated to give the value of the two system parameters d1 and d2, and a plane was fitted based on the data for P M C . In Figure 5, d1 and d2 are the distances between P L 1 C and the fitted plane and between P L 2 C and the fitted plane, respectively, allowing the two system parameters to be calibrated in the same camera coordinate system.
In order to describe the proposed calibration method of this paper, a flowchart with the whole calibration procedure is shown in Figure 6.

4. Experiments and Results

4.1. System Hardware

Calibration experiments were carried out to verify the proposed method using the developed DPMD measurement system, as illustrated in Figure 7. The setup consisted of a CCD camera, two LCD screens, a beam splitter plate, and a mirror with a matrix of hollow ring markers. The camera was the model ECO655CVGE from SVS (Bremen, Germany), and had a resolution of 2448 × 2048 pixels. Two small size LCD screens were used in the DPMD system. The two LCD screens were the model LP097QX2 from LG (Seoul, Korea) and had a resolution of 2048 × 1536 pixels. However, they were low quality and the flatness was not so good. The three angular displacement platforms had a minimum readout of 5′ from Micro-nano Optics Automation Equipment Cooperation (Beijing, China). The mirror with the matrix of hollow ring markers and the beam splitter plate were manufactured by Jiaite Photoelectric Corporation (Shenzhen, China) to match the experimental requirements.

4.2. Calibration experiments

In order to determine the internal parameters of the CCD camera, the checkerboard was placed at eighteen random positions and orientations with a large angle between the imaging axis and the board’s normal. The checkerboard size along the rows and columns had the same value of 6 mm, as illustrated in Figure 8.
At each position, the CCD camera captured the checkerboard image. Using the known checkerboard size in the eighteen captured images, the internal parameters of the CCD camera were obtained by the Camera Calibration Toolbox for Matlab [45], the results are shown in Table 1. The radial and tangential distortions of the following captured images were corrected by the obtained internal parameters. The average reprojection error in direction x and y was 0.0856 pixel and 0.0668 pixel, respectively, as illustrated in Figure 9.
After obtaining the internal parameters of the CCD camera, the calibration procedure included the following steps: (1) The camera captured an image of the mirror with the hollow ring marker matrix to calibrate its extrinsic parameters and determine its orientation in the camera coordinate system; (2) Two groups of orthogonal phase-shifted fringe patterns with the optimum three fringe numbers were displayed on screens LCD1 and LCD2. The camera captured these fringe pattern images in sequence. The absolute maps were obtained from the captured fringe patterns using a four-step phase shifting algorithm and the optimum three-fringe number selection method [46,47]. A matrix with 12 × 9 feature points was designed on the LCD screens to enable the corresponding points on the absolute map to be searched for. The orientations of the two screens were determined in the camera coordinate system based on the phase matching and fringe order constraint principle; (3) The two system parameters d1 and d2 were calculated based on the orientations of the two LCD screens and the reference mirror; (4) Using the known distance between neighboring markers, calibration was achieved by determining the relationship between each pixel position and its XY coordinates.

4.3. Evaluation of Results

A concave mirror, from Micro-nano Optical Corporation (Beijing, China), was measured using the calibrated DPMD system. The center of the concave surface was obtained according to the measured 3D shape data by fitting all the points. The radius could be calculated as a true value. The difference between every point and the obtained true radius was calculated. The average value of all the differences was the calculated systematic error of 26 µm. To evaluate the calibrated system parameters quantitatively, the calibrated DPMD system was also tested using a standard artificial step with specular surfaces, the photos of the objects are illustrated in Figure 10. The specular step was designed and manufactured with known distances between neighboring steps, as shown in the first column in Table 2. Twelve sinusoidal fringe patterns with fringe numbers of 81, 80, and 72 were generated and displayed sequentially on the two LCD screens. The fringe patterns reflected by the specular objects being tested were deformed by the slope and shape of the measured surface. The reflected fringe patterns captured by the CCD camera are shown in Figure 11a for the step and Figure 12a for the concave mirror. Figure 11b and Figure 12b illustrate single absolute phase maps for the specular step and the mirror, respectively. The absolute phase map was then converted into depth data, and the calibrated parameters d1 and d2 were used to obtain the 3D shape data for the specular step and the concave mirror, as shown in Figure 11c and Figure 12c. Figure 11d and Figure 12d show the profile maps of the specular step and concave mirror being tested.
All the points measured on one step surface were fitted onto a plane to calculate the distance between neighboring steps. The measured distance between neighboring steps was calculated using the average of the distance from all the points obtained on the other step surface to the fitted plane. The measured distance, the absolute error (absolute difference between the average distance measured and the calibration distance), and the standard deviation are listed in the second to fourth columns of Table 2.
The experimental results demonstrate that the proposed calibration method can reconstruct the 3D shape of specular objects with high precision and reliability. The absolute error in the distance between neighboring steps was below 22 µm, while the value is 52 µm by using the existing calibration method [37].

5. Conclusions

This paper proposes a novel method of calibrating system parameters for the PMD technique. The method uses a reference mirror with a matrix of hollow ring markers to determine the orientations of two LCD screens in the same camera coordinate system, by displaying orthogonal sinusoidal fringe patterns on the screens using phase data. This allows the two system parameters d1 and d2 to be calibrated in the camera coordinate system. The calibrated DPMD system was tested with an artificial specular mirror with known distances between steps, and a standard concave mirror of known radius. The experimental results validated the effectiveness and reliability of the proposed calibration method for the system parameters.
The proposed calibration method has the following advantages: (1) The parameters of the DPMD system are calibrated using phase information calculated from sinusoidal fringe patterns. The digital phase calibration board has been designed to calibrate the position of LCD screens instead of geometric feature points. Because the plane and screen are different distances from the imaging device, they cannot be clearly captured, given the limited DOF of the lens. Using the phase target can solve this problem efficiently; (2) The calibrated DPMD system can measure specular objects that have discontinuous and/or isolated surfaces; (3) The mirror surface is selected as the reference plane, making it unnecessary to replace components, and ensuring an accurate reference position. Also, the proposed calibration method can obtain the two parameters in the same way, so that it does not need to remove the board in the measuring system after calibration, which is more convenient.
Of course, this paper also considered several factors which would affect the accuracy of whole calibration process. Firstly, Chen’s method was used instead of the traditional ellipse fitting method to calculate the real projection center to reduce the error as much as possible, as it has been proven to have a higher accuracy. Secondly, gamma correction by encoding the initial fringes was applied to reduce the influence of nonlinear distortion. However, the influence of a finite width of the beam splitter, the accuracy evaluation of phase matching principle and other error sources for calibration also need to be studied in further experiments.

Author Contributions

Conceptualization, Z.Z., N.G. and X.D.; methodology, Z.Z.; software, X.D.; validation, X.D.; formal analysis, Z.Z., N.G. and X.D.; investigation, X.D.; resources, Z.Z. and N.G.; data curation, Z.Z., N.G. and X.D.; writing—original draft preparation, X.D.; writing—review and editing, Z.Z., N.G. and X.D.; visualization, X.D.; supervision, Z.Z.and N.G.; project administration, Z.Z.; funding acquisition, Z.Z.

Funding

This research was funded by the National Key R&D Program of China (under Grant No. 2017YFF0106404), the National Natural Science Foundation of China (under Grant No. 51675160), and the Talents Project Training Funds in Hebei Province (under Grant No. A201500503).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ding, Y.; Xi, J.T.; Yu, Y.G.; Chicharo, J. Recovering the absolute phase maps of two fringe patterns with selected frequencies. Opt. Lett. 2011, 36, 2518–2520. [Google Scholar] [CrossRef]
  2. Zuo, C.; Chen, Q.; Gu, G.H.; Feng, S.J.; Feng, F.X.Y. High-speed three-dimensional profilometry for multiple objects with complex shapes. Opt. Express 2012, 20, 19493–19510. [Google Scholar] [CrossRef]
  3. Li, D.; Kofman, J. Adaptive fringe-pattern projection for image saturation avoidance in 3D surface-shape measurement. Opt. Express 2014, 22, 9887–9901. [Google Scholar] [CrossRef]
  4. Xu, J.; Xi, N.; Zhang, C.; Zhao, J.G.; Gao, B.T.; Shi, Q. Rapid 3D surface profile measurement of industrial parts using two-level structured light patterns. Opt. Lasers Eng. 2011, 49, 907–914. [Google Scholar] [CrossRef]
  5. Chatterjee, A.; Bhatia, V.; Pracash, S. Anti-spoof touchless 3D fingerprint recognition using deflectometry and biospeckle analysis. Opt. Lasers Eng. 2017, 95, 1–7. [Google Scholar] [CrossRef]
  6. Huang, S.J.; Zhang, Z.H.; Zhao, Y.; Dai, J.; Chen, C.; Xu, Y.J.; Zhang, E.; Xie, L.L. 3D fingerprint imaging system based on full-field fringe projection profilometry. Opt. Lasers Eng. 2014, 52, 123–130. [Google Scholar] [CrossRef]
  7. Agarwal, S.; Kumal, V.; Shakher, C. Analysis of red blood cell parameters by Talbot projected fringes. J. Biomed. Opt. 2017, 10, 106009. [Google Scholar] [CrossRef]
  8. Wang, Z.Y.; Zhang, Z.H.; Gao, N.; Xiao, Y.J.; Gao, F.; Jiang, X.Q. Single-shot 3D shape measurement of discontinuous objects based on coaxial fringe projection system. Appl. Opt. 2019, 58, A169–A178. [Google Scholar] [CrossRef]
  9. Singh, P.; Chatterjee, A.; Bhatia, V.; Prakash, S. Fringe projection profilometry based secured fingerprint sensor. In Proceedings of the 2018 3rd International Conference on Microwave and Photonics (ICMAP), Dhanbad, India, 9–11 February 2018; pp. 8–11. [Google Scholar]
  10. Jeffrey Kuo, C.F.; Chang, A.; Joseph Kuo, P.C.; Lee, C.L.; Wu, H.C. Applying innovative stripes adaptive detection to three-dimensional measurement of color fringe profilometry. Opt. Commun. 2016, 381, 116–126. [Google Scholar] [CrossRef]
  11. Liu, X.L.; Xiang, P.; Chen, H.L.; He, D.; Gao, B.Z. Strategy for automatic and complete three-dimensional optical digitization. Opt. Lett. 2012, 37, 3126–3128. [Google Scholar]
  12. Chen, C.; Gao, N.; Wang, X.J.; Zhang, Z.H.; Gao, F.; Jiang, X.Q. Generic exponential fringe model for alleviating phase error in phase measuring profilometry. Opt. Laser Eng. 2018, 110, 179–185. [Google Scholar] [CrossRef]
  13. Xiao, Y.L.; Su, X.Y.; You, Z.S. Pose transfer geometrical calibration for fringe-reflection optical three-dimensional measurement. Opt. Commun. 2013, 305, 143–146. [Google Scholar] [CrossRef]
  14. Sun, X.M.; Liu, Y.; Yu, X.Y.; Wu, H.B.; Zhang, N. Three-dimentional measurement for specular reflection surface based on reflection component separation and priority region filling theory. Sensors 2017, 17, 215. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, Y.K.; Lehtonen, P.; Su, X.Y. High-accuracy measurement for small scale specular objects based on PMD with illuminated film. Opt. Laser Technol. 2012, 44, 459–462. [Google Scholar] [CrossRef]
  16. Tang, Y.; Su, X.Y.; Liu, Y.K.; Jing, H.L. 3D shape measurement of the aspheric mirror by advanced phase measuring deflectometry. Opt. Express 2008, 16, 15090–15096. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, L.; Xue, J.; Gao, B.; Zuo, C.; Idir, M. Zonal wavefront reconstruction in quadrilateral geometry for phase measuring deflectometry. Appl. Opt. 2017, 56, 5139–5144. [Google Scholar] [CrossRef]
  18. Niu, Z.Q.; Gao, N.; Zhang, Z.H.; Gao, F.; Jiang, X.Q. 3D shape measurement of discontinuous specular objects based on advanced PMD with bi-telecentric lens. Opt. Express 2018, 26, 1615–1632. [Google Scholar] [CrossRef]
  19. Knauer, M.C.; Kaminski, J.; Hausler, G. Phase measuring deflectometry: A new approach to measure specular free-form surfaces. Proc. SPIE 2004, 5457, 366–376. [Google Scholar]
  20. Häusler, G.; Faber, C.; Olesch, E.; Ettl, S. Deflectometry vs. interferometry. Proc. SPIE 2013, 8788. [Google Scholar] [CrossRef]
  21. Huang, L.; Idir, M.; Zuo, C.; Kaznatcheev, K.; Zhou, L.; Asundi, A. Comparison of two-dimensional integration methods for shape reconstruction from gradient data. Opt. Lasers Eng. 2015, 64, 1–11. [Google Scholar] [CrossRef]
  22. Huang, L.; Asundi, A. Framework for gradient integration by combining radial basis functions method and least-squares method. Appl. Opt. 2013, 52, 6016–6021. [Google Scholar] [CrossRef]
  23. Huang, L.; Xue, J.; Gao, B.; Zuo, C.; Idir, M. Spline based least squares integration for two-dimensional shape or wavefront reconstruction. Opt. Laser Eng. 2017, 91, 221–226. [Google Scholar] [CrossRef]
  24. Huang, L.; Xue, J.; Gao, B.; Mcpherson, C.; Beverage, J.; Idir, M. Model mismatch analysis and compensation for modal phase measuring deflectometry. Opt. Express 2017, 25, 881–887. [Google Scholar] [CrossRef]
  25. Wu, Y.; Yue, H.; Yi, J.; Li, M.; Liu, Y. Dynamic specular surface measurement based on color-encoded fringe reflection technique. Opt. Eng. 2016, 55, 024104. [Google Scholar] [CrossRef]
  26. Zhang, X.; Jiang, L.; Zhang, G. Novel method of positioning optical freeform surfaces based on fringe deflectometry. CIRP Ann.-Manuf. Technol. 2017, 66, 507–510. [Google Scholar] [CrossRef]
  27. Liu, Z.; Yin, Y.; Wu, Q.; Li, X.J.; Zhang, G.J. On-site calibration method for outdoor binocular stereo vision sensors. Opt. Lasers Eng. 2016, 86, 75–82. [Google Scholar] [CrossRef]
  28. Huang, L.; Ng, C.S.; Asundi, A.K. Dynamic three-dimensional sensing for specular surface with monoscopic fringe reflectometry. Opt. Express 2011, 19, 12809–12814. [Google Scholar] [CrossRef]
  29. Jing, H.L.; Su, X.Y.; Liu, Y.K.; Wu, F. Specular surface measurement based on fringe reflection and analysis of 3D shape reconstruction technique. Opto-Electron. Eng. 2008, 35, 37–41. [Google Scholar]
  30. Huang, L.; Idir, M.; Zuo, C.; Kaznatcheev, K.; Zhou, L.; Asundi, A. Shape reconstruction from gradient data in an arbitrarily-shaped aperture by iterative discrete cosine transforms in Southwell configuration. Opt. Lasers Eng. 2015, 67, 176–181. [Google Scholar] [CrossRef]
  31. Yuan, T.; Zhang, F.; Tao, X.P.; Fu, J.J. Three-dimensional shape measuring for specular surface based on phase measuring deflectometry. Acta Opt. Sin. 2016, 36, 0212004. [Google Scholar] [CrossRef]
  32. Song, Z.; Jiang, H.L.; Lin, H.B.; Tang, S.M. A high dynamic range structured light means for the 3D measurements of specular surface. Opt. Lasers Eng. 2017, 95, 8–16. [Google Scholar] [CrossRef]
  33. Huang, L.; Xue, J.; Gao, B.; McPherson, C.; Beverage, J.; Idir, M. Modal phase measuring deflectometry. Opt. Express 2016, 24, 24649–24664. [Google Scholar] [CrossRef]
  34. Ren, H.Y.; Gao, F.; Jiang, X.Q. Iterative optimization calibration method for stereo deflectometry. Opt. Express 2015, 23, 22060–22068. [Google Scholar] [CrossRef]
  35. Xiao, Y.L.; Su, X.Y.; Chen, W.J. Flexible geometrical calibration for fringe-reflection 3D measurement. Opt. Lett. 2012, 37, 620–622. [Google Scholar] [CrossRef]
  36. Zhang, Z.H.; Liu, Y.; Huang, S.J.; Niu, Z.Q.; Guo, J.; Gao, N.; Gao, F.; Jiang, X.Q. Full-field 3D shape measurement of specular surfaces by direct phase to depth relationship. In Proceedings of the SPIE of Optical Metrology and Inspection for Industrial Applications IV, Photonics Asia 2016, Beijing, China, 12–14 October 2016; pp. 12–14. [Google Scholar]
  37. Liu, Y.; Huang, S.J.; Zhang, Z.H.; Gao, N.; Gao, F.; Jiang, X.Q. Full-field 3D shape measurement of discontinuous specular objects by direct phase measuring deflectometry. Sci. Rep. 2017, 7, 10293. [Google Scholar] [CrossRef]
  38. Zhao, P.; Gao, N.; Zhang, Z.H.; Gao, F.; Jiang, X.Q. Performance analysis and evaluation of direct phase measuring deflectometry. Opt. Laser Eng. 2018, 103, 24–33. [Google Scholar] [CrossRef]
  39. Zhang, Z.H.; Guo, J.; Wang, Y.M.; Huang, S.J. Parallel-alignment and correction of two displays in three-dimensional measuring system of specular surfaces. Opt. Precis. Eng. 2017, 2, 289–296. [Google Scholar] [CrossRef]
  40. Zhang, Z.Y. A flexible new technique for camera calibration. IEEE Trans. Pattern Anal. 2000, 11, 1330–1334. [Google Scholar] [CrossRef]
  41. Chen, X.Y.; Ma, Z.; Hu, Y.; Chen, Y.Q.; Bi, F.L. A new method for accurate location of concentric circles in visual measurement. J. Optoelectron.·Laser 2013, 24, 1524–1528. [Google Scholar]
  42. Zhang, Z.H.; Huang, S.J.; Meng, S.S.; Gao, F.; Jiang, X.Q. A simple, flexible and automatic 3D calibration method for a phase calculation-based fringe projection imaging system. Opt. Express 2013, 21, 12218–12227. [Google Scholar] [CrossRef]
  43. Hoang, T.; Pan, B.; Nguyen, D.; Wang, Z.Y. Generic gamma correction for accuracy enhancement in fringe projection profilometry. Opt. Lett. 2010, 35, 1992–1994. [Google Scholar] [CrossRef]
  44. Wang, Z.Y.; Nguyen, D.A.; Barnes, J.C. Some practical considerations in fringe projection profilometry. Opt. Lasers Eng. 2010, 48, 218–225. [Google Scholar] [CrossRef]
  45. Bouguet, J.V. Camera Calibration Toolbox for Matlab. Available online: http://www.vision.caltech.edu/bouguetj/calib_doc/ (accessed on 8 February 2019).
  46. Creath, K. Phase-measurement interferometry techniques. Prog. Opt. 1998, 26, 349–393. [Google Scholar]
  47. Zhang, Z.H.; Towers, C.E.; Towers, D.P. Time efficient color fringe projection system for 3-D shape and colour using optimum 3-frequency selection. Opt. Express 2006, 14, 6444–6455. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the measuring system using direct phase measuring deflectometry (DPMD).
Figure 1. Schematic diagram of the measuring system using direct phase measuring deflectometry (DPMD).
Applsci 09 01444 g001
Figure 2. Image of hollow ring marker matrix captured by the camera: (a) texture image and (b) red dots representing the center of each hollow ring marker.
Figure 2. Image of hollow ring marker matrix captured by the camera: (a) texture image and (b) red dots representing the center of each hollow ring marker.
Applsci 09 01444 g002
Figure 3. Absolute phase on the mirror surface for a horizontal fringe pattern: (a) absolute map with inaccurate phase data for the hollow ring markers and (b) absolute map processed with a 3D surface fitting algorithm.
Figure 3. Absolute phase on the mirror surface for a horizontal fringe pattern: (a) absolute map with inaccurate phase data for the hollow ring markers and (b) absolute map processed with a 3D surface fitting algorithm.
Applsci 09 01444 g003
Figure 4. Feature point matching using absolute phase data.
Figure 4. Feature point matching using absolute phase data.
Applsci 09 01444 g004
Figure 5. Calculation of d1 and d2.
Figure 5. Calculation of d1 and d2.
Applsci 09 01444 g005
Figure 6. The flowchart of the process for the whole calibration.
Figure 6. The flowchart of the process for the whole calibration.
Applsci 09 01444 g006
Figure 7. Setup for the experiments using the DPMD technique.
Figure 7. Setup for the experiments using the DPMD technique.
Applsci 09 01444 g007
Figure 8. The image of the checkerboard.
Figure 8. The image of the checkerboard.
Applsci 09 01444 g008
Figure 9. Reprojection error of the calibrated camera.
Figure 9. Reprojection error of the calibrated camera.
Applsci 09 01444 g009
Figure 10. Photos of the specular objects being measured: (a) manufactured artificial step and (b) concave mirror.
Figure 10. Photos of the specular objects being measured: (a) manufactured artificial step and (b) concave mirror.
Applsci 09 01444 g010
Figure 11. Measurements on the concave mirror: (a) one captured fringe pattern reflected by the concave mirror; (b) one absolute phase map from the reflected fringe patterns; (c) measured 3D shape of the mirror and (d) profile map of the mirror.
Figure 11. Measurements on the concave mirror: (a) one captured fringe pattern reflected by the concave mirror; (b) one absolute phase map from the reflected fringe patterns; (c) measured 3D shape of the mirror and (d) profile map of the mirror.
Applsci 09 01444 g011
Figure 12. Measurements on the artificial specular step: (a) one captured fringe pattern reflected by the step; (b) one absolute phase map from the reflected fringe patterns; (c) measured 3D shape of the step and (d) profile map of the step.
Figure 12. Measurements on the artificial specular step: (a) one captured fringe pattern reflected by the step; (b) one absolute phase map from the reflected fringe patterns; (c) measured 3D shape of the step and (d) profile map of the step.
Applsci 09 01444 g012
Table 1. The calibrated internal parameters of the camera.
Table 1. The calibrated internal parameters of the camera.
CoefficientsFocal Length (pixel)Principal Point (pixel)Radial DistortionTangential Distortion
FuFvPuPvK1K2K3K4
Data10137.45210137.8111204.7831011.174−0.2712−1.8758−0.00020.0003
Table 2. Evaluation results for the artificial specular step (units: mm).
Table 2. Evaluation results for the artificial specular step (units: mm).
Calibration DistanceMeasured DistanceAbsolute ErrorRoot Mean Square Error
3.9873.9650.0180.022
7.0257.0440.0150.019
5.0065.0300.0220.024
6.0996.0790.0210.020

Share and Cite

MDPI and ACS Style

Deng, X.; Gao, N.; Zhang, Z. A Calibration Method for System Parameters in Direct Phase Measuring Deflectometry. Appl. Sci. 2019, 9, 1444. https://doi.org/10.3390/app9071444

AMA Style

Deng X, Gao N, Zhang Z. A Calibration Method for System Parameters in Direct Phase Measuring Deflectometry. Applied Sciences. 2019; 9(7):1444. https://doi.org/10.3390/app9071444

Chicago/Turabian Style

Deng, Xiaoting, Nan Gao, and Zonghua Zhang. 2019. "A Calibration Method for System Parameters in Direct Phase Measuring Deflectometry" Applied Sciences 9, no. 7: 1444. https://doi.org/10.3390/app9071444

APA Style

Deng, X., Gao, N., & Zhang, Z. (2019). A Calibration Method for System Parameters in Direct Phase Measuring Deflectometry. Applied Sciences, 9(7), 1444. https://doi.org/10.3390/app9071444

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop