Abstract
A catadioptric camera with a double-mirror system, composed of a pinhole camera and two planar mirrors, can capture multiple catadioptric views of an object. The catadioptric point sets (CPSs) formed by the contour points on the object lie on circles, all of which are coaxial parallel. Based on the property of the polar line of the infinity point with respect to a circle, the infinity points in orthogonal directions can be obtained using any two CPSs, and a pole and polar pair with respect to the image of the absolute conic (IAC) can be obtained through inference of the Laguerre theorem. Thus, the camera intrinsic parameters can be solved. Furthermore, as the five points needed to fit the image of a circle are not easy to obtain accurately, only sets in which five points can be located can be obtained, whereas the points on the line of intersection between the two plane mirrors and the ground plane can easily be obtained accurately. An optimization method based on the analysis of neighboring point sets to compare the intersection points with an image of the center of multiple circle images fitted using the point sets is proposed. Bundle adjustment is then applied to further optimize the camera intrinsic parameters. The feasibility and validity of the proposed calibration methods and their optimization were confirmed through simulation and experiments. Two primary innovations were obtained from the results of this study: (1) by applying coaxial parallel circles to the double-plane-mirror catadioptric camera model, a variety of calibration methods were derived, and (2) we found that the overall model could be optimized by analyzing the features of the neighboring point set and bundle adjustment.
1. Introduction
Vision technology, a modern research hotspot, has been applied in various engineering fields [1–3]. The acquisition of camera intrinsic parameter data and the solving of these parameters through the process of “camera calibration” has become one of the most important issues in the field of vision technology. Camera calibration is an essential step in the process of recovering 3D spatial information from 2D image plane information [4]. To expand the small field-of-view (FOV) of traditional cameras, Baker and Nayar [5] proposed the catadioptric camera system, which can expand the FOV by adding a mirror in front of the camera lens. The catadioptric camera is extensively used in several fields, such as 3D modeling [6–8], virtual monitoring [9, 10], and robot navigation [11].
A plane catadioptric system [12, 13] consists of one or more plane mirrors with a pinhole camera. Camera calibration and 3D reconstruction can be performed using a catadioptric system with a single plane mirror and a pinhole camera [14]; however, the system can only provide two views of an object, limiting the acquisition of image information. Therefore, double-plane mirrors are used to obtain multiple catadioptric views of an object. Several methods have been proposed for determining the internal parameters of a double-mirror camera system. Gluckman and Nayar [15] only solved the focal length of the camera by determining the relationship between the corresponding points between views. As it is not easy to extract the corresponding points experimentally, Forbes et al. [16] extracted the silhouettes of objects from images with multiple views and used the epipolar tangents of these silhouettes to estimate the internal camera parameters, avoiding the extraction of the corresponding points between views; however, the physical coordinates of the silhouettes need to be precisely located. Zhang et al. [17] presented the geometric relationship between the double-plane-mirror catadioptric system and circular motion. Hernández et al. [18] further improved on the theoretical system developed by Zhang et al. [17] by proposing the use of circular motion to determine the camera’s internal parameters, but this method requires the direction of the rotation axis or observation. In addition, Hu et al. [19] used the relationship between the corresponding points in a catadioptric camera with a double-mirror system and the silhouettes of the views to solve the camera intrinsic parameters. Ying et al. [20] presented the pole-polar relationship of the silhouettes between five views and solved the camera intrinsic parameters by first determining the coordinates of the vanishing point corresponding to the three coordinate axes and then using the approaches developed in the studies described above. Circles are common in daily life, and the usage of circles in camera calibration has been extensively researched [21–24]. Zhao and Li [25] used the imaged circular points (ICPs) and the orthogonal vanishing points to solve the camera intrinsic parameters. Although some suggestions have been made in the recent literature, the problem of solving for the vanishing point in the direction of the double mirrors remains to be solved.
Referring to the above literature, this study proposes camera calibration methods based on coaxial parallel circles. In a double-plane-mirror catadioptric camera system, the catadioptric point set (CPS) formed by any point on an object is located on the same circle [21], and these circles are coaxial parallel. The camera intrinsic parameters can be solved according to the corresponding linear constraints. This paper proposes an approach that improves upon the methods described in the literature review above by incorporating the use of coaxial parallel circles into the double-plane-mirror catadioptric camera model; furthermore, a variety of calibration methods are applied to obtain more accurate camera intrinsic parameters through the use of nonlinear optimization.
2. Imaging Model
This section briefly reviews the imaging model of a double-plane-mirror catadioptric camera [20], the camera projection process, and the property of the polar of the infinity point with respect to a circle.
2.1. Imaging Model of a Double-Plane-Mirror Catadioptric Camera
The imaging model comprises a real pinhole camera and two real plane mirrors and , as shown in Figure 1 [20]. When the included angle of the double mirrors ranges from to , single internal reflection occurs. Here, two virtual planar mirrors are donated as and , the image of in is virtual mirror , and the image of in is virtual mirror . represents reflection of in the real and virtual mirrors. The images of in and are virtual cameras and , respectively. The image of in is virtual camera , and the image of in is virtual camera . As an object will form four catadioptric virtual images, an image is captured with camera to obtain five views of the object. View is the image of the real object with , views and are the images reflected by the object in and , respectively, and views and are the images of and in and , respectively. According to the principle of planar mirror imaging, the four groups of views (), (), (), and () are symmetrical to , , , and , respectively. Similarly, cameras and views have the same symmetry. Let the intersection lines between , , , and and the ground plane be , , , and , respectively.

2.2. Camera Projection Process
In Figure 1, the world coordinate system from Section 2.1, , and is on the same plane and takes the intersection point between this plane and the line of intersection between and as the origin, . The direction of the line connecting and is the -axis, the -axis is perpendicular to the plane of , and the -axis is determined by the right-hand coordinate system. Defining in as the homogeneous coordinates of a space point and as the homogeneous coordinates of the image point on the corresponding 2D image, the projection process can be expressed as where is the scale factor, and is a projection matrix, which can be expressed as
In a double-mirror catadioptric camera, the real camera also has four virtual reflections. It realizes a stereo vision system with fixed internal parameters and varying external parameters. Therefore, where is the projection matrix of , is the 3D rotation matrix, and is the translation vector. and are the external parameters of camera ; is the camera intrinsic parameter matrix of camera , which is a upper triangular matrix: where is the homogeneous coordinate of the principal point, and are the focal lengths in the and directions, respectively, and is the skew factor between the and axes of the image plane.
2.3. Property of the Polar Line of the Infinity Point with respect to a Circle
Theorem 1. [26]. If and are a group of conjugate diameters of a circle, infinity point in the direction of and infinity point in the direction of are a pair of orthogonal infinity points.
Proposition 2. The direction of an infinity point is orthogonal to the direction of the polar line of the infinity point with respect to a circle.
Proof. In Figure 2, the infinity point is denoted as , and the polar line of with respect to circle is the diameter . The infinity point in the direction of diameter is denoted as , and the polar line of with respect to circle is denoted as diameter . From polar transformation in the projective space [26], it is known that is on diameter ; and are a group of conjugate diameters of circle , which intersect at the center . According to Theorem 1, points and are a pair of infinity points in orthogonal directions.

3. Calibration of the Double-Mirror Catadioptric Camera
This section details the calibration methods proposed in this study. We discuss the relationship between the image of the CPS (ICPS) and the orthogonal vanishing point, the ICPs, and the pole and polar of the image of the absolute conic (IAC), respectively.
As shown in Figure 3, in the double-mirror catadioptric camera system, let be any real point in space and the images of point in mirrors be points , respectively, where and are the reflected points of ; and are the catadioptric points of .
Definition 3. If is any real point in space and is the corresponding catadioptric points in mirrors , respectively, is called a CPS.

In Figure 3, and are the circles fitted by and , respectively. Based on the symmetry of plane-mirror imaging, and are coaxial parallel circles [21]. It is known that the line passing through and is parallel to the line passing through and and intersects at infinity point . Similarly, the line passing through and is parallel to the line passing through and and intersects at infinity point . The infinity line through points and is . The polar lines of points and with respect to circle intersect infinity line at infinity points and , respectively. According to Proposition 2 mentioned in Section 2.3, and are two pairs of infinity points in orthogonal directions.
The description above and in Figure 3 represents the spatial configuration of an image with corresponding spatial constraints. See the Appendix for the corresponding definitions, propositions, and proofs.
4. Nonlinear Optimization
This section includes two parts: the first involves the optimization of the fitting of and and application to the calibration algorithm. The next involves the consideration of the camera intrinsic parameters obtained through the optimization algorithm as the initial value and iteration for optimizing the camera intrinsic parameter calculation.
4.1. Optimization of the Conics Fitted in the Image
Analyzing the particularity of biplane-mirror imaging, only five fitting points can be selected from the image of a circle connecting the object point and its views; therefore, it is crucial to extract these five points accurately.
Definition 4. On the image plane, the set of pixels that represent a space point in a view is called the effective pixels of the point.
In Figure 1, the four intersecting lines , , , and under camera can be extracted using the least-squares fitting method [27] to calculate the coordinates of the pairwise combinatorial intersection points, and their average value can be obtained as . Let and be any two points on an object in contact with the ground plane. The image of and in and is and and and , respectively. The image of and in is and , respectively; the image of and in is and , respectively.
Select , , as the effective pixels of images , of space points , on view , and fit , conics, respectively (i.e., the two families of conics). Among them, randomly select two conics and and calculate their generalized eigenvectors, three of which represent three points in the image [28]. The points in the two conics are denoted as . The optimization objective is as follows: where dist represents the Euclidean distance between two points in {}, and min represents taking the minimum value in {}. corresponding to , and , corresponding to , i.e., the optimal conic, can be obtained from the above two families of conics. Then, the pixel points required for fitting and can be obtained.
Similarly, using the two conics and and applying the same method, the edge space points along the outskirts of the object can be optimized, including certain points in contact without the ground, to fit circle images of their projections.
The vanishing line and the circular-point images are recovered using the images of the two concentric circles, and the images of the dual circular points denoted as are calculated. Then, in view , according to the effective pixels corresponding to the space points, multiple conics can be fitted, i.e., a family of conics in which the value of is related to the number of selected effective pixels. The intersection points among a family of conics and the vanishing lines are solved, respectively, to obtain the dual of their intersection points, denoted as . Then, the optimization objective is as follows: determine corresponding to and corresponding to , which is the optimal conic selected from the family of conics. Thus, the five pixels required for fitting can be obtained.
4.2. Optimization of the Camera Internal Parameters
Through bundle adjustment [29], the camera intrinsic parameters can be optimized using the reprojection points on a circle. In general, the coordinates of the center of a circle and its radius can determine a circle.
Therefore, images and fitting conics are selected for each image. According to the settings of the world coordinate system depicted in Figure 1, the coordinates of the center of each circle can be represented as , with . The coordinates of points on the circle corresponding to views can be set as the angle parameter ; therefore, points can be expressed as follows:
The above optimization problem can be changed as where is the scale factor; is the projection point of , which can be obtained as described in Section 4.1, and and can be obtained [20]. The next step is to solve the minimum value of the cost function to optimize the internal parameter matrix.
The LM method is used to solve this optimization problem by selecting the initial value of the function (8) and its expected optimization target. Using the optimization method described in Section 4.1, some more accurate pixel points can be obtained to fit the circle images. As the calculated initial values of the camera intrinsic parameters do not deviate considerably from the actual values, the obtained camera intrinsic parameters can be used as the initial values.
5. Proposed Algorithm for Recovering the Euclidean Structure
In accordance with the above discussion, the main steps of the calibration algorithm are as follows:
Input: adjust the camera position to obtain images
Output: matrix of the intrinsic parameters of the camera,
Step 1. Extract the effective pixel coordinates of the two ICPSs and and fit the two conics and , according to the least squares method [30] and the optimization method presented in Section 4.1.
Step 2. Use points and on and to solve the vanishing line , according to (A.1), (A.2), and (A.3).
Step 3. Solve the orthogonal vanishing points and using (A.4) and (A.5).
Step 4. Solve ICPs and using (A.1) and (A.2); obtain vanishing points and vanishing line using (17–21).
Step 5. If and are obtained, solve using (A.14) (method 1); if and are obtained, solve using (A.16) (method 2); if and are obtained, solve using (A.13) (method 3).
Step 6. Obtain the camera intrinsic parameters through the Cholesky factorization of and the solution of the inverse matrix.
Step 7. Further, optimize the internal parameters as described in Section 4.2.
6. Experimental
To verify the feasibility and effectiveness of the proposed algorithm, we performed simulation and experiments. Ying et al. [21] proposed a self-calibration theory based on the analysis of the structural characteristics obtained by imaging using a double-plane-mirror camera; Zhao and Li [25] subsequently modified this approach to meet the requirements of specific applications and described some newer methods developed by other researchers. We therefore compared results obtained using Ying et al.’s [21] and Zhao and Li’s [25] methods.
6.1. Simulated Data
In the simulation, the camera intrinsic parameters were set as follows:
The angle between the two plane mirrors was set within . Two points in the included corner were selected randomly, and their reflections were generated and projected on the image plane. Gaussian noise with zero mean and standard deviation (noise level) was added to each image point, and was varied from pixels. For each , 100 independent experiments were performed, and the absolute error of the camera intrinsic parameters was calculated. The obtained results are shown in Figures 4(a)–4(e). It can be observed that the absolute errors of increase with the increase in the . At the same , the absolute errors of our algorithms are similar to those of Ying et al.’s [21] and Zhao and Li’s [25] algorithms.

(a)

(b)

(c)

(d)

(e)
Figure 5 depicts the absolute error of the camera intrinsic parameters under different , where optimization method 1, method 2, and method 3 are denoted as Method1 OP, method 2 OP, and method 3 OP, respectively. Under the same , as the absolute errors are nearly equal, Figure 5 shows only one, similar to . After optimization, it can be seen that the absolute error of the calculated camera parameters is lower.

(a)

(b)

(c)
6.2. Actual Experiment
We first took pictures of the checkerboard grid from three different poses, as shown in Figures 6(a)–6(c), using rectangular mirrors, which could be stood up more easily. Camera parameters calculated using reference [31] were used as the ground truth. To carry out calibration, a doll placed was placed between the two plane mirrors with the intrinsic parameters left unchanged. A double mirror with an included angle of (within ) was placed on a horizontal desktop such that there was no occlusion between the five views. Multiple pictures were captured from different positions, and three images with a better effect were selected to solve the camera intrinsic parameters. The three selected images are displayed in Figures 7(a)–7(c) with an image resolution of (pixels). The results of Canny edge detection [32] are shown in Figures 7(d)–7(f). The silhouette and the silhouette coordinates of the five views in each image are shown in Figures 7(g)–7(i) and 7(j)–7(l), respectively.

(a)

(b)

(c)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)
The silhouette outlines and epipolar tangency lines of the five views in Figures 7(a) and 7(b) are shown in Figures 8(a) and 8(b), respectively. The image coordinates of the silhouette outlines and the intersections of the epipolar tangency lines in the image are shown in Figures 8(c) and 8(d). Therefore, the coordinates of the ICPSs can be obtained in each image [25].

(a)

(b)

(c)

(d)
According to the abovementioned procedure, the three algorithms proposed in this study as well as the calibration algorithms in references [21, 25] were used to perform 100 independent experiments. In addition, the method described in Chapter 4 was used to optimize and average the three algorithms proposed in this study. The average values of the results are listed in Table 1. It is seen that the results are similar, with the results of method 1 OP, method 2 OP, and method 3 OP closer to the ground truth, proving the effectiveness of the respective calibration algorithms.
According to the characteristics of plane-mirror imaging and the reconstruction method in reference [16], we used the visible silhouettes of the five views to obtain the epipolar geometric relationship between the views, determine the relationship between the rotation and translation between views, and reconstruct the approximate 3D structure of the object. In this study, we selected the images in Figures 8(a) and 8(b) to reconstruct the 3D model of the doll. Because the calculation results obtained using the different methods were similar and those obtained experimentally using method 1 were better, only the camera intrinsic parameters data obtained using method 1 and method 1 OP (Table 1) were for used for the 3D reconstruction. The reconstruction results are shown in Figures 9. Table 2 lists the norm of deviation values between the real and reconstructed dolls for several corresponding feature points in the same reference frame. From Figures 9 and Table 2, it is easy to see that the reconstruction results are very similar to the doll depicted in Figures 7(a) and 7(b). The right and back views of the reconstruction results obtained using method 1 (Figures 9(a) and 9(b), respectively) are highly similar to those obtained using method 1 OP (Figures 9(c) and 9(d), respectively).

(a)

(b)

(c)

(d)
7. Conclusions
In this study, we proposed methods for obtaining the camera intrinsic parameters of a double-mirror catadioptric camera system. When the included angle of the two mirrors is between , five views of an object can be obtained in an image. The camera intrinsic parameters can be obtained by capturing three images from different positions. According to the symmetry of plane-mirror imaging, the circles fitted by some of the CPSs are coaxial parallel. Using two CPSs, two infinity points can be obtained, whose connection provides an infinity line. Using the property that the direction of the infinity point is orthogonal to the direction of its polar line with respect to a circle, two pairs of infinity points in the orthogonal direction can be obtained. Moreover, the infinity line and circles fitted by the CPSs intersect at the circular points. The orthogonal vanishing points or the image of the circular points can be obtained on the image plane. In addition, based on the pole-polar relationship and inference of the Laguerre theorem, the image of the common axis of the coaxial parallel circles can be obtained; therefore, the vanishing point can be determined in the axis direction, which is in the direction of the normal vectors of the projection plane including the parallel circles. Both of them have a pole-polar relationship with respect to the IAC. Therefore, if two pairs of ICPSs are extracted on the image plane, the three calibration methods can be applied to solve the camera intrinsic parameters. Analysis of the features of the neighboring point set indicate that the centers of the circles comprising the corresponding catadioptric five-point sets are on a straight line; hence, the image of the circle fitted on the image can be optimized. The determined camera parameters are considered as the initial values and further optimized through binding adjustment. In simulation, the absolute errors of our methods are slightly lower than those obtained using Ying et al.’s [21] and Zhao and Li’s [25] methods, and the optimized absolute errors converge rapidly, resulting in estimated values that are stable around exact values with a small degree of jitter. Comparison of our results with the ground truth value obtained using real data indicated that the methods have low absolute error. In performing 3D reconstruction, the two methods for example produced similar norms of deviation from ground truth.
Our experimental approach led to several limitations in terms of results. First, to ensure that five points within each catadioptric mirror would be on the same plane, our experiments were limited to the use of rectangular mirrors. As a result, we were unable to examine cases in which the five points were not on the same plane, as would occur with the use of nonrectangular mirrors. Second, we did not consider cases in which the shape features of an imaged object were partially occluded. Finally, no experiments were carried out in overexposed or dark settings.
Appendix
Corresponding Definitions, Propositions, and Proofs
In this section, the constraints on the corresponding image are described in detail using the spatial scenario described in Section 3.
Proposition 5. In the double-mirror catadioptric imaging system, if two ICPSs are given, two sets of orthogonal vanishing points can be determined.
Proof. As shown in Figure 10, the projections of CPSs and on the image plane are the image points and , respectively, and the fitted conics and are the images of circles and , respectively. The line passing through points and intersects with the line passing through points and at vanishing point , which is the image of ; therefore, Similarly, vanishing point is the image of : In the projection plane containing circle , the vanishing line through and is , which is the image of : The polar line of points and with respect to circular image intersects with vanishing line at vanishing points and , respectively, which are the images of infinity points and , respectively. According to the principle of projective invariance, and are two pairs of orthogonal vanishing points.
Proposition 6. In the double-mirror catadioptric imaging system, if the vanishing line of the planes containing two coaxial parallel circles is determined, a group of ICPs and a group of pole and polar with respect to the IAC can be determined.
Proof. In Figure 10, the intersection points of vanishing line and conic are ICPs and , respectively, on the projection plane containing circle . The equations of the simultaneous vanishing line and and , respectively, are as follows:

Theoretically, the solutions of (A.6) and (A.7) are the same for and . The presence of noise during actual shooting can cause the solutions of the ICPs to differ; therefore, the average of the two solutions is taken as the ICP. In addition, the harmonic conjugate is invariant in the perspective projective and the inference of the Laguerre theorem interprets the relationship between infinite and circular points. Therefore, a pole and polar pair with respect to the IAC can be obtained according to the inference of the Laguerre theorem [28].
In Figure 10, is the image of the intersection line between the two plane mirrors and , which is the connecting line between the images of the two centers of and . The vanishing line and the images of the two centers and have a pole-polar relationship with respect to and , respectively:
Therefore, the connecting line of and is
The intersection of lines and is the vanishing point:
Line is perpendicular to the line passing through the optical center of camera ; inferring from the Laguerre theorem [28] and the projective invariance principle, there exists where is the vanishing point in the normal direction of the plane determining line and the optical center of camera . and have a pole-polar relationship with respect to the IAC [28], given by
Because and are a pair of orthogonal vanishing points,
and the ICPs lie on the IAC:
Because and are a pair of conjugate points, where Re and Im are the real and imaginary parts, respectively.
Equation Equations (A.14), (A.16), and (A.13) provide two independent constraints for the elements of [26]. To completely calibrate the internal parameters, at least three images are needed to estimate based on Propositions 5 and 6.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (NSFC) under Grant Nos. 61663048 and 11861075.