Comparison of Camera Calibration and Measurement Accuracy Techniques for Phase Measuring Deﬂectometry

: Phase measuring deﬂectometry (PMD) is a competitive method for specular surface measurement that offers the advantages of a high dynamic range, non-contact process, and full ﬁeld measurement; furthermore, it can also achieve high accuracy. Camera calibration is a crucial step for PMD. As a result, a method based on the calibration of the entrance pupil center is introduced in this paper. Then, our proposed approach is compared with the most popular photogrammetric method based on Zhang’s technique (PM) and Huang’s modal phase measuring deﬂectometry (MPMD). The calibration procedures of these three methods are described, and the measurement errors introduced by the perturbations of degrees of freedom in the PMD system are analyzed using a ray tracing technique. In the experiment, a planar window glass and an optical planar element are separately measured, and the measurement results of the use of the three methods are compared. The experimental results for the optical planar element (removing the ﬁrst 6 terms of the Zernike polynomial) show that our method’s measurement accuracy reached 13.71 nm RMS and 80.50 nm PV, which is comparable to accuracy values for the interferometer.


Introduction
Phase measuring deflectometry (PMD) with structured light projection and the phaseshifting technique is a highly accurate optical surface measuring method based on the law of light reflection. It has the advantages of high speed, low cost, and a high dynamic range; furthermore, it requires a simple device and is non-destructive [1][2][3][4][5][6][7][8][9][10]. In the PMD system, rays are considered to be emitted from the projection center of a camera lens and reflected off a planar surface under test (SUT) to reach a liquid crystal display (LCD) screen. During the measurement process, the calibration step is carried out to establish the relationships among the screen coordinate system, the coordinate system of the SUT, and the camera coordinate system; then, the slope of the reflection points on the SUT are obtained according to the slope calculation equation [9,10]. Finally, the slope data are integrated to obtain the surface information [3][4][5][6][7][8][9][10]. Among these variables, the most important are the coordinate positions of the screen, the SUT and the projection center of the camera. The calibration error in them will degrade the measurement accuracy of the SUT.
At present, an external pinhole is utilized as the projection center of the camera, whose position can be accurately measured with a laser tracker in PMD or the Software Configurable Optical Test System (SCOTS) [3,11,12]. The projection center of the camera lens can also be calibrated based on the photogrammetry principle [13,14]. Based on Zhang's technique, Xiao et al. proposed a calibration method with a flat mirror to obtain the position of the projection center of the camera [15,16]. Additionally, Hou et al. proposed a method to calculate the position of the projection center in a unified world coordinate system from the external parameters (including the translation vector and rotation matrix) of the camera [16,17]. Zhang's technique requires taking several images from different orientations and different poses, but it is not always accessible in deflectometry, especially for a machine-integrated system.
Optically, the chief rays reflected off the SUT will certainly converge at the entrance pupil center of the camera lens [18][19][20][21]. The role of the projection center in the PMD system, as is mentioned in the references [13][14][15][16][17], is equivalent to this center. For the sake of distinction, we named it the entrance pupil center in our method, but it is still called the projection center in other methods.
For an actual lens, the position of the entrance pupil center can vary as the off-axis angle of an incoming ray increases [18]. However, for most of the lens, this effect is not significant, and thus it can usually be ignored in camera calibration [19]. Therefore, a method based on the calibration of the entrance pupil center is introduced in this paper. The entrance pupil calibration method (EPCM) that we propose is compared with the popular photogrammetric method based on Zhang's technique (PM) [16], and Huang's MPMD method (MPMD) [7]. A planar window glass about 90 mm in diameter and an optical planar element approximately 103 mm in diameter are measured separately by the three methods, and the measurement results are analyzed.
The rest of the paper is organized as follows: Section 2 illustrates the basic principle of the three methods (EPCM, PM, and MPMD). Section 3 describes the measurement procedures of the three methods. Section 4 presents the calculation of the perturbations of the degrees of freedom in the PMD system. Section 5 shows the experimental results, including the comparison of the camera calibration and measurement accuracy. Section 6 discusses and summarizes this paper.

The Basic Principle of Phase Measuring Deflectometry
The PMD system is composed of an LCD screen, a camera, and the SUT. Structured light patterns displayed by the screen are reflected by the SUT and captured by the camera. Then, the slope data can be calculated by the system's geometrical model and calibration parameters. Finally, the slope data is integrated to obtain the SUT shape. The coordinates of the screen pixels and the coordinates of the SUT instead of the reference surface are obtained for the triangulation calculation. The slope data of the SUT is calculated as shown in Equation (1) [10,22]: where x m and y m are the coordinates of the reflection points; x s and y s are the coordinates of the screen pixels, which can be calculated by the phase-shifting algorithm [23,24]; x c and y c are the coordinates of the camera; z m2s and z m2c are the distance between the SUT and the screen and distance between the SUT and the camera, respectively; and d m2s = Then, the SUT shape can be obtained by integrating the slope data using the Zernike modal method [25][26][27].

Entrance Pupil Calibration Method (EPCM)
In Figure 1a, because both the entrance pupil center and the external stop lie on the optical axis in the paraxial region of the camera lens, a reference surface and an external stop are used to obtain the entrance pupil center, which is usually a virtual image of the aperture stop of the camera lens and is inside the lens. The origin of the global coordinate In Figure 1a, because both the entrance pupil center and the external stop lie on the optical axis in the paraxial region of the camera lens, a reference surface and an external stop are used to obtain the entrance pupil center, which is usually a virtual image of the aperture stop of the camera lens and is inside the lens. The origin of the global coordinate system − is established as the point , which is the coordinate of the external stop on the camera lens. Additionally, point is the entrance pupil center. Based on the geometry shown in Figure 1a, where the points and are along with the optical axis, from the view of ray tracing in reverse we extend the ray of to the reflection point ( , ) on the reference surface and reach spot on the LCD screen. These four points, , , , and , are along the optical axis. Assuming that the reference surface is an ideal plane, according to the law of reflection, the incident angle is Based on the geometry shown in Figure 1a, where the points E and O are along with the optical axis, from the view of ray tracing in reverse we extend the ray of EO to the reflection point M r0 (x mr0 , y mr0 ) on the reference surface and reach spot S on the LCD screen. These four points, O, E, M r0 , and S, are along the optical axis. Assuming that the reference surface is an ideal plane, according to the law of reflection, the incident angle is equal to the reflection angle. The reflection points of the reference surface can be obtained, as shown in Equation (2): x mr = x sr z m2o +x o z m2s z m2o +z m2s y mr = y sr z m2o +y o z m2s where z m2s and z m2o are the distance between the reference surface and the screen and the distance between the reference surface and point O in the z direction, respectively; x sr and y sr are the coordinates of the screen pixels corresponding to the central area of the image plane (as shown in Figure 1b) Figure 1a. They are perpendicular to the x-axis and the y-axis, respectively. M r1 and M r2 are equal to point O and point E in the Y direction, respectively. Similarly, B 1 and B 2 are equal to point O and point E in the Y direction. As shown in Figure 1a, the angle relationships between the camera and the reference surface can be expressed in triangles M r1 OB 1 and M r0 OM r1 , as shown in Equation (3): Based on triangles M r0 OM r1 and M r0 EM r2 , triangles M r1 OB 1 and M r2 EB 2 are similar to each other. We assume that the distance between point O and point E is d p . According to the above, the coordinates of the entrance pupil center can be determined in the global coordinate system, as shown in Equation (4): where x e , y e , and z e are the coordinates of the entrance pupil center, which are the coordinates of the camera, α 2 = α 1 and β 2 = β 1 . After obtaining the camera coordinates, the reflection points coordinates can be solved using an iterative reconstruction strategy [28].

Photogrammetric Method Based on Zhang's Technique (PM)
In PM, the camera's internal and external parameters are calculated from feature points with known world coordinates and corresponding image coordinates. As a result, external parameters (rotation matrix R and translation vector T) are obtained using calibration technology. Assuming that the world coordinates of the point M = (X W , Y W , Z W , 1) T and its image coordinates m = (u, v, 1) T are known, the relationship between M and m can be described as where f x and f y are the equivalent focal lengths in pixels along the x and y directions, α is a skew factor, and (u 0 , v 0 ) are the coordinates of the principal point. These parameters are called the camera intrinsic parameters A. The rotation matrix R and translation vector T are the external parameters. In Hou's method [17], assuming the coordinates of a point are M c in the camera coordinate system and M w in the world coordinate system, the relationship between M c and M w is expressed by the following equation: The coordinates of the projection center in the camera coordinate system are (0, 0, 1), and its world coordinates can be obtained by the following equation, In E's method [10], which can achieve a high measurement accuracy, a checkboard is displayed on another LCD screen with several different orientations, the camera takes pictures of the checkboard in different poses, and matrix A and [R T] are obtained with the calibration toolbox [29]. Hence, Equation (8) can be obtained by substituting M = (X W , Y W , Z W , 1) T in Equation (5) It is worth noting that the last pose of the checkerboard plane is located at the identical position as the SUT for getting rid of misalignment, as it makes the local coordinate system on the last checkerboard parallel to each axis of the global coordinate system. Since H is a homography matrix in Equation (8), the mirror pixel coordinates can be determined using the following equation: In addition, the system errors can be eliminating by measuring a high-quality reference surface that has a known shape and is aligned at an identical position to the SUT [10]. This measurement can be mathematically stated as follows: where P t and P r are the measurement of the SUT and the reference surface, respectively, and are directly obtained from the PMD measurement method, which contains the system errors error_t and error_r, respectively. I t is the measured surface shape of the test flat after the system errors are eliminated. I r is the surface shape of the reference element and can be measured by Fizeau interferometer.

Modal Phase Measuring Deflectometry (MPMD)
In MPMD, the height and slopes of the SUT are represented by the well-established mathematical model and updated by optimizing the model coefficient c to minimize the discrepancy between the re-projection in ray tracing and the actual measurement [7]. As shown in Figure 2, the incident ray → p is reflected off the SUT to the screen, and the reflected ray → r can be obtained with Equation (11): where → n , → p , and → r are the unit vector of the normal, the incident ray and the reflected ray, respectively, and <> represents the dot product.
The position of the screen relative to the camera (rotation matrix R m , translation vector T m ) can be obtained through a pre-calibration procedure, which includes camera calibration, screen calibration, and geometric calibration. Once the pose of the screen is initialized via pre-calibration, the corresponding intersection S τ of the screen is calculated. The coefficient vector c is initialized using the mathematical model of the SUT.
The coefficient vector c and the pose of the screen (R m , T m ) are optimized by minimizing the differences between the re-projection S τ and the measurement S on the screen among the N pairs of correspondence in the least squares sense. The nonlinear least squares problem in Equation (12) can be solved using the Levenberg-Marquardt algorithm. ) among the N pairs of correspondence in the least squares sense. The nonlinear least squares problem in Equation (12) can be solved using the Levenberg-Marquardt algorithm.

EPCM
According to the principle of the entrance pupil calibration, the position of an external stop on the camera is measured. The EPCM is introduced in detail as follows:

EPCM
According to the principle of the entrance pupil calibration, the position of an external stop on the camera is measured. The EPCM is introduced in detail as follows: The coordinates of point E can be obtained by Equations (3) and (4) in the global coordinates; (d) The points of the reference surface are used as the initial values of the reflection points; (e) The slope is calculated using Equation (1) and the reconstruction height w from the slope. It is output after the height is less than the threshold; (f) The surface shape is reconstructed.
The calibration procedure of EPCM is illustrated in Figure 3.  (3) and (4) in the global coordinates; (d) The points of the reference surface are used as the initial values of the reflection points; (e) The slope is calculated using Equation (1) and the reconstruction height w from the slope. It is output after the height is less than the threshold; (f) The surface shape is reconstructed.
The calibration procedure of EPCM is illustrated in Figure 3.

PM
The procedure of PM is introduced in detail as follows: (a) The test system is established. It is worth noting that the last pose of the checkerboard plane is located at an identical position to the SUT, and it makes the local coordinate

PM
The procedure of PM is introduced in detail as follows: (a) The test system is established. It is worth noting that the last pose of the checkerboard plane is located at an identical position to the SUT, and it makes the local coordinate system on the last checkerboard parallel to each axis of the global coordinate system; (b) The rotation matrix R and translation vector T are obtained by camera calibration; (c) The projection center is obtained with Equation (7), and the coordinates of the SUT are calculated with Equation (9); (d) The slope is calculated with Equation (1), and the surface shape of the reference surface and the SUT are reconstructed; (e) The systematic error is eliminated with Equation (10).
The procedure of PM is illustrated in Figure 4a.

MPMD
The procedure of MPMD is introduced in detail as follows: (a) The incident ray ⃗ is calibrated using the system calibration. The pose of the screen relative to the camera ( , ) can be obtained through the pre-calibration procedure; (b) The initial coefficient matrix c and Zernike polynomial are used to represent the surface shape of the SUT; (c) The incident ray ⃗ is reflected from the SUT to the screen. The reflected ray ⃗ can be calculated with Equation (5), and the intersection on the screen is obtained; (d) The coefficient vector c and pose of the screen ( , ) are optimized by minimizing the differences between the re-projection and the measurement S of the screen among the N pairs of correspondence in the least squares sense; (e) Output the shape coefficient and reconstruct the surface shape.
The procedure of MPMD is illustrated in Figure 4b.

Accuracy Analysis
The following section describes the analysis of the measurement errors due to the

MPMD
The procedure of MPMD is introduced in detail as follows: (a) The incident ray → p is calibrated using the system calibration. The pose of the screen relative to the camera (R m , T m ) can be obtained through the pre-calibration procedure; (b) The initial coefficient matrix c and Zernike polynomial are used to represent the surface shape of the SUT; (c) The incident ray → p is reflected from the SUT to the screen. The reflected ray → r can be calculated with Equation (5), and the intersection S τ on the screen is obtained; (d) The coefficient vector c and pose of the screen (R m , T m ) are optimized by minimizing the differences between the re-projection S τ and the measurement S of the screen among the N pairs of correspondence in the least squares sense; (e) Output the shape coefficient and reconstruct the surface shape.
The procedure of MPMD is illustrated in Figure 4b.

Accuracy Analysis
The following section describes the analysis of the measurement errors due to the perturbations of the degrees of freedom in the PMD system. Generally, the measurement result depends on the accuracy of the measurement tool and device. Therefore, it is necessary to analyze measurement errors caused by each coordinate value, which is called the degree of freedom (DOF). The error analysis aids in the identification of which DOF need to be tightened and which might be loosened. Thus, there are three DOF (x c , y c , z c ) for the camera coordinates. For the coordinate measurement of the screen pixels, there are six DOFs; these are translations (x s , y s , z s ) and rotations (screen tilt x, screen tilt y, screen tilt z) in three directions. For the coordinate measurement of the SUT, the rotation of the SUT around the z axis does not cause measurement error; only the measurement results need to be rotated. Therefore, there are five DOFs with three translations (x m , y m , z m ) and two rotations (flat tilt x, flat tilt y) for the SUT. Regardless of the errors by the LCD screen, such as thermal distortion, roughness, and noise, the 14 DOFs need to be considered [9].
As mentioned above, the DOFs are analyzed in the simulation. The schematic diagram of the PMD system is shown in Figure 5. Assuming that the pinhole O coordinates are the origin of the world coordinate system, the distance between the SUT and the camera O is 1525 mm, and the distance between the LCD screen and the SUT is 1500 mm in the z direction, while the angle between the camera and the z direction is 5 • . Korsch ray tracing [30] is used to calculate the slope of an ideal planar element in the x and y directions. For DOF, initial values are considered to be the accurate values. For the coordinate values of the camera (x c , y c , z c ), the screen pixel (x s , y s , z s ), and the SUT (x m , y m , z m ), a perturbation of 1 mm is artificially added to the accurate value of each DOF, respectively. Similarly, a perturbation of 1 mrad is artificially added to the accurate value for the five DOFs of the rotations, respectively. When error is added artificially to each DOF, the others remain unchanged.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 16 two rotations (flat tilt x, flat tilt y) for the SUT. Regardless of the errors by the LCD screen, such as thermal distortion, roughness, and noise, the 14 DOFs need to be considered [9]. As mentioned above, the DOFs are analyzed in the simulation. The schematic diagram of the PMD system is shown in Figure 5. Assuming that the pinhole coordinates are the origin of the world coordinate system, the distance between the SUT and the camera is 1525 mm, and the distance between the LCD screen and the SUT is 1500 mm in the z direction, while the angle between the camera and the z direction is 5°. Korsch ray tracing [30] is used to calculate the slope of an ideal planar element in the x and y directions. For DOF, initial values are considered to be the accurate values. For the coordinate values of the camera( , , ), the screen pixel( , , ), and the SUT( , , ), a perturbation of 1 mm is artificially added to the accurate value of each DOF, respectively. Similarly, a perturbation of 1 mrad is artificially added to the accurate value for the five DOFs of the rotations, respectively. When error is added artificially to each DOF, the others remain unchanged. Based on the above analysis, the Zernike polynomial coefficients are analyzed after the surface shape is reconstructed with the Zernike polynomial [25,26] under different DOFs with errors. The differences from the 4th to 11th terms of the Zernike polynomial coefficient between the results of the accurate surface figure and the results gained using different perturbations are shown in Table 1. The terms from the 4th to 11th Zernike polynomial coefficients are defocus, oblique astigmatism, vertical astigmatism, vertical coma, horizontal coma, vertical trefoil, oblique trefoil, and primary spherical, respectively. The root sum square (RSS) is the value for evaluating the effects on the measurement results due to the perturbation of each DOF or in every term of the Zernike polynomial coefficient.  Based on the above analysis, the Zernike polynomial coefficients are analyzed after the surface shape is reconstructed with the Zernike polynomial [25,26] under different DOFs with errors. The differences from the 4th to 11th terms of the Zernike polynomial coefficient between the results of the accurate surface figure and the results gained using different perturbations are shown in Table 1. The terms from the 4th to 11th Zernike polynomial coefficients are defocus, oblique astigmatism, vertical astigmatism, vertical coma, horizontal coma, vertical trefoil, oblique trefoil, and primary spherical, respectively. The root sum square (RSS) is the value for evaluating the effects on the measurement results due to the perturbation of each DOF or in every term of the Zernike polynomial coefficient. In general, the first eight terms of the Zernike polynomial coefficients are greatly affected when the DOF can be not obtained accurately. Therefore, according to the pre-set perturbation for all DOFs shown in the first column in Table 1, the accuracy of the alignment can be estimated in this measuring system. When the positions of the camera, the screen, and the SUT are separately changed under a given perturbation, these alignment errors for the 4th to the 11th Zernike polynomial coefficients can be calculated, as shown in Table 1. This table shows that the coordinate error of z e and z s will induce some relatively large errors in terms of defocus, while screen tilt error about the z axis will produce relatively significant errors in terms of astigmatism; however, the other terms are so small that they can be ignored.

Comparison of Camera Calibration
The coordinates of the screen are obtained by the phase-shifting algorithm in three methods. For EPCM, the distance d p between the external stop and the front surface is obtained as d p = 23.96 mm (see Appendix A). The experimental setup of the PMD is shown in Figure 6, where the LCD screen on the optical experiment platform displays a sixteenstep sinusoidal fringe pattern. The resolution of the CCD camera is 1296 × 966 pixels, and one pixel is 3.75 µm. The light emitted from the LCD screen is reflected by a high-quality flat element and captured by a camera with an external stop, which can be removed. The coordinates of the external stop can be accurately positioned with CMM and the camera focus on the element regardless of the screen being out of focus.
Three methods, including EPCM, PM, and MPMD, are selected for comparison. In EPCM, the camera coordinates are obtained with Equation (4) according to the geometric system in Figure 1a. In PM, the 16 poses of the checkerboard are captured (the last one is restored at an identical position to the SUT); then, the internal and external parameters are obtained by the calibration toolbox [29]. Finally, the projection center of the camera is calculated with Equation (7). Likewise, in MPMD, Xiao's algorithm [15] is used to obtain the relationship between the camera and the screen. Additionally, we then calculate the camera coordinates with Equation (7), which is equivalent to the PM in principle. According to the calibration processes mentioned above, the camera coordinates are obtained in EPCM, MPMD, and PM, respectively. Additionally, they are transformed into the same coordinate system, as shown in Table 2. It can be seen from Table 2 that the camera coordinates obtained by using the above three methods are different in the x, y, and z directions. According to Table 1, if the camera coordinate in the x direction is inaccurate, the error will mainly be induced in the 4th and 6th terms of the Zernike polynomial; inaccuracy in the y direction will be largely embodied in the 5th term; and error in the z direction will be expressed in the 4th term and will have a great influence on the reconstruction accuracy. Therefore, if the three sets of coordinates in Table 1 are separately used to measure the SUT, the inaccuracy of the camera coordinates will be mostly reflected in the defocus and astigmatism error. ment can be estimated in this measuring system. When the positions of the camera, the screen, and the SUT are separately changed under a given perturbation, these alignment errors for the 4th to the 11th Zernike polynomial coefficients can be calculated, as shown in Table 1. This table shows that the coordinate error of and will induce some relatively large errors in terms of defocus, while screen tilt error about the z axis will produce relatively significant errors in terms of astigmatism; however, the other terms are so small that they can be ignored.

Comparison of Camera Calibration
The coordinates of the screen are obtained by the phase-shifting algorithm in three methods. For EPCM, the distance between the external stop and the front surface is obtained as = 23.96 mm (see Appendix A). The experimental setup of the PMD is shown in Figure 6, where the LCD screen on the optical experiment platform displays a sixteen-step sinusoidal fringe pattern. The resolution of the CCD camera is 1296 × 966 pixels, and one pixel is 3.75 μm. The light emitted from the LCD screen is reflected by a high-quality flat element and captured by a camera with an external stop, which can be removed. The coordinates of the external stop can be accurately positioned with CMM and the camera focus on the element regardless of the screen being out of focus.

Comparison of Measurement Accuracy
In the experiments, the three methods (EPCM, PM, and MPMD) are used to measure the planar window glass with a diameter of about 90 mm. According to the above three methods, the surface shape results of the planar window glass are obtained after removing the first three terms of Zernike polynomials, as shown in Figure 7. the planar window glass with a diameter of about 90 mm. According to the above three methods, the surface shape results of the planar window glass are obtained after removing the first three terms of Zernike polynomials, as shown in Figure 7. Figure 7a-e shows the reconstruction results obtained by the EPCM, PM with systematic errors, the PM after eliminating systematic errors, and the MPMD. The measurement result obtained by the interferometer are shown in Figure 7e, and its root-mean-square (RMS) and peak-to-valley (PV) values are 519.77 nm and 2328.59 nm.  The surface shape shown in Figure 7a is closest to the interferometer, and its RMS value and PV value are 558.04 nm and 2505.05 nm, respectively. Compared with the interferometer, the RMS and PV values of the surface shape in Figure 7b are 317.58 nm and 1479.31 nm, and the inaccuracy of the defocus term is directly presented. However, the surface shape of Figure 7c is obtained by eliminating the systematic error of Figure 7b, and its RMS value and PV value are 559.63 nm and 2516.01 nm, respectively. This is closer to the interferometer than that in Figure 7b. In Figure 7d, compared with interference, it can be seen that there are serious errors in both defocus and astigmatism.
According to the analysis in Table 1, not only the inaccuracy of the camera coordinates but also the rotation of the screen around the z-axis has a great impact on the reconstruction results, which introduce errors in astigmatism. As a result, the reconstruction results of the planar window glass are obtained after removing the first 6 terms of Zernike polynomials, as shown in Figure 8. I It can be seen from Figure 8 that Figure 8a-c are relatively closer to Figure 8e, while the low older error in Figure 8d is not yet removed due to the initial value choice of the shape coefficient and accuracy limit of the method. In a word, for the PM with systematic errors and MPMD, error of the low-order term is generated due to inaccurate coordinates, but the EPCM still has a high measurement accuracy.
To further compare their accuracy performance, an optical planar element with a diameter of about 103 mm is measured. The surface shapes are obtained after removing the first three terms of the Zernike polynomials, as shown in Figure 9. Compared with the interferometer, the errors in Figure 9a,c are mainly reflected in the astigmatism term. It can be seen from Table 1 that the error is mainly caused by the rotation of the screen around the z-axis. Both Figure 9b,d have serious errors in terms of the defocus. These are mainly caused by the inaccuracy of z c and z s , according to Table 1. In order to remove the errors mentioned above, the surface shapes are obtained after removing the first six terms of the Zernike polynomials, as shown in Figure 10. It can be directly shown from Figure 10 that Figure 10a,c are comparable to the interferometers in shapes and magnitudes. Due to the limitation of calibration accuracy, both Figure 10b,d have significant errors compared with the interferometer. Additionally, in our previous work (see Appendix A), for the same element, the results of removing the first four terms were obtained, and the RMS and PV were 27.70 nm and 187.40 nm (f /8), respectively, which are very close to the interferometer results (24.86 nm RMS and 154.73 nm PV). In order to compare the accuracy of the three methods more clearly in these two experiments, the measurement results are summarized in Table 3. struction results, which introduce errors in astigmatism. As a result, the reconstruction results of the planar window glass are obtained after removing the first 6 terms of Zernik polynomials, as shown in Figure 8. I It can be seen from Figure 8 that Figure 8a-c ar relatively closer to Figure 8e, while the low older error in Figure 8d is not yet removed due to the initial value choice of the shape coefficient and accuracy limit of the method. In a word, for the PM with systematic errors and MPMD, error of the low-order term is gen erated due to inaccurate coordinates, but the EPCM still has a high measurement accu racy. To further compare their accuracy performance, an optical planar element with a di ameter of about 103 mm is measured. The surface shapes are obtained after removing th first three terms of the Zernike polynomials, as shown in Figure 9. Compared with th interferometer, the errors in Figure 9a,c are mainly reflected in the astigmatism term. I can be seen from Table 1 that the error is mainly caused by the rotation of the screen around the z-axis. Both Figure 9b,d have serious errors in terms of the defocus. These ar mainly caused by the inaccuracy of and , according to Table 1. In order to remove th errors mentioned above, the surface shapes are obtained after removing the first six term of the Zernike polynomials, as shown in Figure 10  . In order to compare the accuracy o the three methods more clearly in these two experiments, the measurement results ar summarized in Table 3.

Discussion and Conclusions
We compared three methods in this study. Inaccurate system calibration is one of the major error sources in PMD, which limits its practical applications. In EPCM, the camera position is no longer calibrated by Zhang's technology; the entrance pupil center is accurately calculated as a point source by EPCM. The measurement accuracy is high for measuring the low order terms of the surface if the experimental misalignment can be controlled very well. In PM, the accuracy of Zhang's technique is easily affected by the quality of the calibration target and the accuracy of feature point extraction. There are also too many parameters to be addressed in the optimization of the non-linear camera model, which makes it difficult to accurately obtain the projection center of the camera lens. The projection center obtained by this method is unstable, but the accuracy is high after eliminating the systematic error [10]. In MPMD, the pre-calibration result is only used as an initial value in the post-optimization, and the screen pose can be optimized with the shape estimation for a mono-PMD system. If a good pre-calibration is obtained, the MPMD method will obtain a relatively low shape error. However, the reconstruction result will be worse if the pre-calibration is performed poorly, including an initial guess far from the reality.
It is generally believed that good results of camera calibration imply a higher accuracy of the reconstructed surface shape. Hence, some methods [7,8,11,15,21] have been developed to calibrate parameters by reducing the calibration error as much as possible. This paper shows that EPCM has a higher measurement accuracy in comparison with the other two methods.

Discussion and Conclusions
We compared three methods in this study. Inaccurate system calibration is one of the major error sources in PMD, which limits its practical applications. In EPCM, the camera position is no longer calibrated by Zhang's technology; the entrance pupil center is accurately calculated as a point source by EPCM. The measurement accuracy is high for measuring the low order terms of the surface if the experimental misalignment can be controlled very well. In PM, the accuracy of Zhang's technique is easily affected by the quality of the calibration target and the accuracy of feature point extraction. There are also too many parameters to be addressed in the optimization of the non-linear camera model, which makes it difficult to accurately obtain the projection center of the camera lens. The projection center obtained by this method is unstable, but the accuracy is high after eliminating the systematic error [10]. In MPMD, the pre-calibration result is only used as an initial value in the post-optimization, and the screen pose can be optimized with the shape estimation for a mono-PMD system. If a good pre-calibration is obtained, the MPMD method will obtain a relatively low shape error. However, the reconstruction result will be worse if the pre-calibration is performed poorly, including an initial guess far from the reality.
It is generally believed that good results of camera calibration imply a higher accuracy of the reconstructed surface shape. Hence, some methods [7,8,11,15,21] have been developed to calibrate parameters by reducing the calibration error as much as possible. This paper shows that EPCM has a higher measurement accuracy in comparison with the other two methods.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
In Figure A1, the experimental setup for measuring the entrance pupil position is composed of a light source, the camera lens, and a microscope, which is used to measure the distance between the entrance pupil and the front surface of the camera lens. In the measurement, when the microscope is adjusted to focus on the entrance pupil and the first surface of the camera lens, respectively, two positions are recorded. Considering the distance between the external stop and the front surface of camera lens, the distance (d p ) between the external stop and the entrance pupil of the camera lens can be obtained, as d p = 23.96 mm.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
In Figure A1, the experimental setup for measuring the entrance pu composed of a light source, the camera lens, and a microscope, which is u the distance between the entrance pupil and the front surface of the cam measurement, when the microscope is adjusted to focus on the entrance pu surface of the camera lens, respectively, two positions are recorded. Consi tance between the external stop and the front surface of camera lens, th between the external stop and the entrance pupil of the camera lens can = 23.96 mm.

Microscope
The lens light source Guide rail