ISRN Machine Vision

Volume 2013 (2013), Article ID 874084, 7 pages

http://dx.doi.org/10.1155/2013/874084

Research Article

Vision Measurement Scheme Using Single Camera Rotation

Shidu Dong

College of Computer Science and Engineering, Chongqing University of Technology, Chongqing 400050, China

Received 10 April 2013; Accepted 8 May 2013

Academic Editors: A. Nikolaidis and J. P. Siebert

Copyright © 2013 Shidu Dong. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We propose vision measurement scheme for estimating the distance or size of the object in static scene, which requires single camera with 3-axis accelerometer sensor rotating around a fixed axis. First, we formulate the rotation matrix and translation vector from one coordinate system of the camera to another in terms of the rotation angle, which can be figured out from the readouts of the sensor. Second, with the camera calibration data and through coordinate system transformation, we propose a method for calculating the orientation and position of the rotation axis relative to camera coordinate system. Finally, given the rotation angle and the images of the object in static scene at two different positions, one before and the other after camera rotation, the 3D coordinate of the point on the object can be determined. Experimental results show the validity of our method.

1. Introduction

Nowadays, digital camera or mobile phone with camera is very popular. It is appealing and convenient, if they are utilized to estimate the distance or size of an object. For this purpose, stereo images with disparity should be taken [1, 2]. One obvious method for the stereo image acquisition is using two cameras with different view angles. With two images of an object from two cameras and the relative orientation and position of the two different viewpoints, through the correspondence between image points in the two views, the 3D coordinates of the points on the object can be determined [3, 4]. But, in general, mobile phone or professional camera has only one camera and cannot acquire two images from different views simultaneously. Fortunately, there have been many methods for stereo vision system with single camera. The methods may be broadly divided into three categories. First, to obtain virtual images from different viewpoints, additional optical devices are introduced, such as two planar mirrors [5], a biprism [6], convex mirrors [1, 7], or the double lobed hyperbolic mirrors [8]. But these optical devices are expensive and space consuming. Second, 3D information of an object is inferred directly from a still image under the knowledge of some geometrical scene constraints such as planarity of points and parallelism of lines and planes [9–11] or prior knowledge about the scene obtained from the supervised learning [12]. Nevertheless, these methods require constrained scenes or extra computation for training the depth models. Third, 3D information is extracted from sequential images with respect to camera movement, which is often adopted in robot area. Due to the uncertainties in the sequential camera position, however, it is difficult to get the accurate 3D information in that method [1].

In this paper, we propose a novel vision measurement method for estimating the distance or size of the object in static scene, which requires single camera with 3-axis accelerometer rotating around a fixed axis. Through the 3-axis accelerometer sensor, the slope angle of the camera relative to gravity direction can be obtained [13]. The angle can uniquely determine the position of the camera if the camera is rotated around a fixed axis which is not parallel to gravity. Moreover, the relative position and orientation of the camera between two positions, one before and the other after camera rotation, can be determined by the rotation angle. Therefore, at the given two positions, if the camera is calibrated, the classical binocular view methods [4] can be adopted to extract the 3D coordinates of the points on objects. Unfortunately, it is very difficult for the user to rotate the camera to the same slope angle as the calibrated one.

To deal with this problem, we firstly formulate the rotation matrix and translation vector from one coordinate system of the camera to another in terms of the rotation angle. Then, with camera calibration data at various positions, we calculate the orientation and position of the rotation axis relative to the camera coordinate system. Accordingly, at given two positions, the rotation matrix and translation vector from one coordinate system of the camera to another can be calculated by rotation angle, and with the collected two images at the two positions, the 3D coordinate of the points on the object can be determined.

The paper is organized as follows: Section 2 provides a stereo vision system from single camera through rotation, and Section 3 indicates calculation method of the rotation matrix and translation vector from one coordinate system of the camera to another by the rotation angle. The calculation method for the position and orientation of rotation axis relative to the camera coordination system is proposed in Section 4, and the complete calibration and the 3D coordinate calculation of the point on an object are presented in Section 5. Experimental results are given in Section 6, and some conclusions and discussions are drawn in Section 7.

2. Stereo Vision through Camera Rotation

In what follows, we adopt the following hypotheses.

 (H1) The camera is provided with a 3-axis accelerometer sensor, whose readouts, , , and , are the divisions of gravity on the , , and axes of the sensor coordinate system, respectively. (H2) In calibration and measurement processes, the camera rotates around the fixed rotation axis which is parallel to the axis of the sensor coordinate system and is not parallel to the direction of gravity.

Thus, in the course of the rotation, the readout keeps steady. As a result, the division of gravity on the plane of the sensor coordinate system also keeps steady (see Figure 1). Therefore, the slope angle with respect to gravity division (position, for short) can be determined by

Figure 1: Position and rotation angle of the camera. denotes gravity division on the - plane of the sensor coordinate system. and denote the divisions of gravity on the and axes of the sensor coordinate system, respectively, and and denote the camera positions. indicates the rotation angle.

From position to , the rotation angle of device is governed by

Two images, and , of an object are collected at the positions and , respectively. and denote the camera coordinate systems at and , respectively. () and () denote the coordinates of a point on the object relative to - and -, respectively, and () and () denote image coordinates of the point of the object on and (image coordinate for short), respectively.

The projection of object coordinates relative to - and - into image coordinates is summarized by the following forms:

Let

where and denote rotation matrix and translation vector between - and -.

Substituting (3) into (4), we get

Thus, the 3D coordinate of the point on the object relative to - (object coordinate) can be determined by

where

From (6), we can see that the object coordinate can be found provided that the intrinsic parameters , , , and , the rotation matrix , and translation vector are available.

With the camera calibration method proposed by Zhang [14], the camera intrinsic parameters and extrinsic parameters describing the camera motion around a static scene can be figured out. Let and denote the extrinsic parameters at and and the extrinsic parameters at . For simplicity, = (), = (), and = (, and ) stand for the coordinates of the point on the object relative to the world coordinate system, - and -, respectively. Thus,

From (8), we get

Substituting (9) into (8), we get

Equation (4) can be rewritten as

Thus,

Based on the previous discussions, the object coordinate can be determined by the following process.

(1)The camera is calibrated, and the intrinsic parameters, , , , and , are obtained.(2)At the positions and , the camera is calibrated, and the corresponding extrinsic parameters, , , , and , are obtained. Then, and can be acquired from (12).(3)In the course of measurement, the camera is rotated carefully to the positions and . At these two positions, the two images of the object, and , are collected. Once the image coordinates of the object, () and (), are known, the object coordinate can be figured out by (6).

It should be noticed that it is rather difficult for the user to rotate the camera accurately to the positions and , where camera is calibrated.

3. Calculation Method of Rotation Matrix and Translation Vector

Intuitively, the orientation and position of - relative to -, and , depend on the camera rotation angle, . Thus, and may be figured out by .

Let - be the coordinate system associated with the rotation axis, which is unit vector on the axis of the coordinate system. Let and denote the orientation and position of - relative to - (orientation and position of the rotation axis). Its homogeneous matrix can be written as

The rotation matrix for the camera rotation around the axis of - from to (device rotation matrix) can be modeled as

Its homogeneous matrix can be written as

The coordinate system transformation can be represented as a graph (Figure 2). A directed edge represents a relationship between two coordinate systems and is associated with a homogeneous transformation. From Figure 2, we can get

where

Figure 2: Coordinate system transformation represented by a graph.

Substituting (13), (15), and (17) into (16), we get

From (20) and (21), we can see that the and can be calculated by the rotation angle provided that and are available.

4. Calculation Method of Position and Orientation of Rotation Axis

4.1. Orientation of Rotation Axis

Suppose that the camera be calibrated at different positions. At position , we get camera calibration data. Let , and let denote the device rotation matrix with respect to rotation angle . and denote rotation matrix and translation vector between camera coordinate systems at the positions and . Equation (20) can be written as

When the rotation axis is fixed, would be constant. Thus, given the values of , , and , we can solve (22) for .

Using a normalized quaternion to define the rotation between two coordinate systems provides a simple and elegant way to formulate successive rotations [15–17]. Given rotation matrix

it can be transformed to the quaternion with the following equation [18]:

Similarly, the quaternion can be transformed to the rotation matrix with the following equation [18]:

Let , , and denote the quaternion of , , and , respectively. With quaternion, the sequence of rotation can be formulated as an equation without involving rotation matrices [15]. As a result, the problem of solving can be transformed into an equivalent problem involving the corresponding quaternion as follows [15]:

Since the quaternion multiplication can be written in matrix form and with notations introduced in [19], we have the following [16]:

where, letting ,

Moreover, these two matrices are orthogonal [16], that is,

Thus,

where

Thus, the total error function allowing us to compute becomes

where is the number of positions of the camera:

is the unit quaternion. Therefore, can be obtained by solving the following problem:

4.2. Position of Rotation Axis

At position , (21) can be written as

When the rotation axis is fixed, would be constant. Thus, given the values of and , we can solve (35) for .

Let , ], , and thus

Let denote the number of positions of the camera. The total error function is

Since the camera rotation axis is approximately vertical to , the value of approaches zero. Thus, can be obtained by solving the following problem:

5. Calculation of Position and Orientation of Rotation Axis and 3D Coordinate

5.1. Calculation of Position and Orientation of Rotation Axis

Based on the previous discussions, the complete process of calculation of position and orientation of rotation axis is outlined below:

(1)The chessboard (see Figure 3) is printed and plastered on plane. By rotating and moving the camera properly, a set of chessboard images is collected. Then, the values of the intrinsic parameters of the camera, , , , and , are obtained by calling the method proposed in [14].(2)The camera is fixed on a camera tripod, and the rotation axis of the camera, which is parallel to the axis of the sensor coordinate system, lies in an approximately horizontal plane. (3)By rotating the camera around the fixed axis to the different positions, another set of chessboard images are collected.(4)The extrinsic parameters of the camera at the positions , , and are obtained by calling the function cvFindExtrinsicCameraParams2 in Opencv [20].(5)The rotation matrix and translation vector between the camera coordinate systems at the positions and , , and are figured out by (12). The rotation angle and its corresponding rotation matrix are also calculated by (2) and (14).(6)The and are converted into quaternions and , respectively, by using the method proposed by Bar-Itzhack [18].(7) is found by solving problem (34). As a result, the can be obtained.(8) is obtained by solving problem (38).

Figure 3: Chessboard.

5.2. 3D Coordinate Calculation

Based on the previous discussions, we present the complete process of 3D coordinate calculation as follows:

(1)The camera is fixed on a camera tripod whose rotation axis lies in an approximately horizontal plane.(2)At certain position , the image of the object is collected. (3)By rotating the camera on the fixed axis to another position , we get another image of the object .(4)The rotation angle and its corresponding rotation matrix are figured out by (2) and (14).(5)With , and are calculated by (20) and (21). (6)The image coordinate of the point of the object, (), on the image is appointed manually.(7)The corresponding image coordinate in the image , (), can be determined by stereo correspondence method, for example, the function FindStereoCorrespondenceBM in Opencv [20] or by manual. (8)The 3D coordinate of the point on the object relative to , (), can be figured out by   (6).

6. Experiments

Since digital camera with 3-axis accelerometer sensor is not available for us, IPhone 4, which has the sensor, is adopted. To simulate digital camera, the phone is placed in a box which is fixed on a camera tripod. And it is ensured that the axis of the sensor (experiment results show that the axis is parallel to the optical axis of the camera) is parallel to the rotation axis of the camera tripod, so that the value of keeps steady in the course of the rotation. In calibration course, the distance between the camera and the chessboard is about 1000 mm.

Figure 4 illustrates the curve of the quaternion of the rotation matrix between the camera coordinate systems at the th position and 1st position, with respect to the rotation angle. The quaternion was calculated by the proposed method with rotation angle, while ] was converted directly from rotation matrix , which was from calibration data. From the graph, one can see that the proposed method can calculate the rotation matrix by rotation angle.

Figure 4: Curve of the quaternion of rotation matrix with respect to rotation angle; is converted from rotation matrix, and is calculated by the proposed method.

Figure 5 plots the curve of translation vector between the camera coordinate systems at the th position and 1st position, with respect to rotation angle. The vector was calculated by the proposed method with rotation angle, while was directly from calibration data . From the graphs, one can see that the proposed method can calculate effectively translation vector by rotation angle.

Figure 5: Curve of translation vector with respect to rotation angle; was directly from calibration data, and was calculated by (23) with rotation angle.

In order to estimate the accuracy of 3D coordinate calculated by the proposed method, the chessboard which has a bigger block than the one for calibration is printed. The width of the block is 46.6 mm, and the distance between the chessboard and the camera is about 1200 mm. The distance between two neighbor corners rather than the distance between the camera and the chessboard is calculated, because the measurement of the former by manual is easier. For simplicity, the corners of the blocks in images are automatically detected by calling the Opencv function “cvFindChessboardCorners” [19]. Figure 6 depicts the measurement error ratio of the distance between two corners with respect to rotation angle.

Figure 6: Measurement error ratio (%) of the distance between two corners with respect to rotation angle.

7. Conclusions

This paper proposed a stereo vision system with single camera, which requires digital camera with 3-axis accelerometer sensor rotating around the fixed axis which is parallel to the axis of the sensor. Under these conditions, the slope angle relative to gravity, which can be figured out from the readouts of the sensor, can determine the camera position, and the rotation angle between two positions can determine the rotation matrix and translation vector from one coordinate system of the camera to another. Accordingly, given the rotation angle and the images of the object at two different positions, one before and the other after camera rotation, the 3D coordinates of the points on the object can be determined. Theoretical analysis and experimental results show the validity of our method.

It should be noticed that few digital cameras are provided with 3-axis accelerometer sensor. However, to obtain stereo vision, we believe that the inexpensive sensor embedded in digital camera is worthy. Moreover, due to higher image quality and larger focus range, higher accuracy and larger range of measurement may be obtained. Furthermore, the smart phone which has the sensor is popular. If a mini fixed rotation axis is built in a corner of the phone and it does not move in the course of rotation, with the proposed method, the phone may estimate the size of the object being focused on and distance between the phone and the object.

Acknowledgment

This work is supported by Science and Technology Research Project of Chongqing’s Education Committee (KJ110806).

References

Y. Sooyeong and N. Ahuja, “An omnidirectional stereo vision system using a single camera,” in Proceedings of the 18th International Conference on Pattern Recognition (ICPR '06), pp. 861–865, August 2006.

D. Marr and T. Poggio, “A computational theory of human stereo vision,” Proceedings of the Royal Society of London, vol. 204, no. 1156, pp. 301–328, 1979.

F. Tombari, S. Mattoccia, L. D. Stefano, and E. Addimanda, “Classification and evaluation of cost aggregation methods for stereo correspondence,” in Proceedings of the 26th IEEE Conference on Computer Vision and Pattern Recognition (CVPR '08), pp. 1–8, June 2008.

H. C. Longuet-higgins, “A computer algorithm for reconstructing a scene from two projections,” Nature, vol. 293, no. 5828, pp. 133–135, 1981.

J. Gluckman and S. K. Nayar, “Catadioptric stereo using planar mirrors,” International Journal of Computer Vision, vol. 44, no. 1, pp. 65–79, 2001.

D. Lee and I. Kweon, “A novel stereo camera system by a biprism,” IEEE Transactions on Robotics and Automation, vol. 16, no. 5, pp. 528–541, 2000.

T. Svoboda and T. Pajdla, “Epipolar geometry for central catadioptric cameras,” International Journal of Computer Vision, vol. 49, no. 1, pp. 23–37, 2002.

D. Southwell, A. Basu, M. Fiala, and J. Reyda, “Panoramic stereo,” International Conference on Pattern Recognition, vol. 1, pp. 378–382, 1996.

A. Criminisi, I. Reid, and A. Zisserman, “Single view metrology,” International Journal of Computer Vision, vol. 40, no. 2, pp. 123–148, 2000.

A. Criminisi, “Single-view metrology: algorithms and applications,” in Proceedings of the 24th DAGM Symposium on Pattern Recognition, vol. 2449 of Lecture Notes in Computer Science, pp. 224–239, 2002.

G. Wang, Z. Hu, F. Wu, and H. Tsui, “Single view metrology from scene constraints,” Image and Vision Computing, vol. 23, no. 9, pp. 831–840, 2005.

A. Saxena, S. H. Chung, and A. Y. Ng, “3-D depth reconstruction from a single still image,” International Journal of Computer Vision, vol. 76, no. 1, pp. 53–69, 2008.

J. Rekimoto, “Tilting operations for small screen interfaces,” in Proceedings of the 9th ACM Annual Symposium on User Interface Software and Technology, pp. 167–168, November 1996.

Z. Zhang, “A flexible new technique for camera calibration,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 22, no. 11, pp. 1330–1334, 2000.

J. C. K. Chou and M. Kamel, “Finding the position and orientation of a sensor on a robot manipulator using quaternions,” International Journal of Robotics Research, vol. 10, no. 3, pp. 240–254, 1991.

F. Dornaika and R. Horaud, “Simultaneous robot-world and hand-eye calibration,” IEEE Transactions on Robotics and Automation, vol. 14, no. 4, pp. 617–622, 1998.

K. H. Strobl and G. Hirzinger, “Optimal hand-eye calibration,” in Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS '06), pp. 4647–4653, October 2006.

I. Y. Bar-Itzhack, “New method for extracting the quaternion from a rotation matrix,” Journal of Guidance, Control, and Dynamics, vol. 23, no. 6, pp. 1085–1087, 2000.

M. W. Walker, L. Shao, and R. A. Volz, “Estimating 3-D location parameters using dual number quaternions,” Image Understanding, vol. 54, no. 3, pp. 358–367, 1991.

B. Gary and K. Adrian, Learning OpenCV, O'Reilly Media, 2008.