US20020186897A1  Camera calibration system using planar concentric circles and method thereof  Google Patents
Camera calibration system using planar concentric circles and method thereof Download PDFInfo
 Publication number
 US20020186897A1 US20020186897A1 US09919086 US91908601A US2002186897A1 US 20020186897 A1 US20020186897 A1 US 20020186897A1 US 09919086 US09919086 US 09919086 US 91908601 A US91908601 A US 91908601A US 2002186897 A1 US2002186897 A1 US 2002186897A1
 Authority
 US
 Grant status
 Application
 Patent type
 Prior art keywords
 camera
 circle
 concentric
 calibration
 points
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Abandoned
Links
Images
Classifications

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
 G06T7/00—Image analysis
 G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration

 H—ELECTRICITY
 H04—ELECTRIC COMMUNICATION TECHNIQUE
 H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
 H04N17/00—Diagnosis, testing or measuring for television systems or their details
 H04N17/002—Diagnosis, testing or measuring for television systems or their details for television cameras
Abstract
The invention relates to a camera calibration system and method thereof, which is capable of easily performing camera calibration using a concentric circle pattern. According to the invention, a method of calibrating a camera calibrates the internal parameters of the camera required to measure geometric information of an object using projection invariable characteristic of concentric circles. The method comprises a first step of taking images of the calibration pattern consisting of two or more concentric circles located at the same plane and having different radius at different angles to obtain projection images of ellipses, a second step of calculating the central point of the projected concentric circles using a given algorithm, and a third step of calculating the principal point and focal point of camera using a nonlinear minimization algorithm based on the central point thus obtained.
Description
 [0001]The invention relates generally to a method of calibrating a camera for obtaining the camera parameters required when the camera is used to obtain geometric information. More particularly, the invention relates to a camera calibration system and method thereof, which use a concentric circle pattern as an artificial calibration pattern in the process of calibrating the camera and obtain camera parameters by analyzing two images of the concentric circle pattern taken at different angles.
 [0002]To obtain geometric information of an object using a camera, it is required to calibrate the camera by estimating parameters between the image information obtained from the camera and the actual geometric information of the object. FIG. 1 illustrates an image of a 3D artificial calibration object model in a conventional imethod of calibrating a camera using 3D calibration object and FIG. 2 is an image of a plane model in a conventional method of calibrating a camera using coordinates of points on a twodimensional plane.
 [0003]Conventional methods for calibrating a camera includes {circle over (1)} a calibration method using a 3D calibration object, {circle over (2)} self calibration method, {circle over (3)} a calibration method using coordinates of points on a twodimensional plane.
 [0004]First, a calibration method widely used so far is one to calibrate the camera using a 3D artificial calibration object, such as a rectangular parallelepiped calibration object, as shown in FIG. 1. As can been seen from FIG. 1, a photograph of the rectangular calibration object is taken to obtain a geometric relation of the rectangular parallelepiped object.
 [0005]However, in the method of calibrating a camera using a 3D artificial calibration object, it is difficult to manufacture and maintain the calibration object of a rectangular parallelepiped. This is because the calibration object should have the characteristic of a normal rectangular parallelepiped in order to calculate the parameters of the camera from its image. In other words, it must be a rectangular parallelepiped in which three plane and twelve edges must maintain a right angle with respect to a vertex. Otherwise, exact calculation of the camera parameter becomes difficult and the reliability of calibration is degraded accordingly.
 [0006]Second, the self calibration method is one to calculate the parameters of the camera using only information of each of corresponding points from several sheets of images. Although this self calibration method can be widely applied without the limitation of using an artificial calibration object, it is difficult to exactly define the corresponding points. Due to this, this method has problems that the process of finding a solution is very complicated and it is also difficult to find a correct solution.
 [0007]Third, as shown in FIG. 2, the method of calibrating the coordinates of the points on a plane is relates to the calibration method using artificial calibration object and the self calibration method. As shown in FIG. 2, this method estimates the parameters of the camera by taking several images of the plane pattern to exactly find the coordinates of the points on a plane and the coordinates of the corresponding points on a corresponding image plane.
 [0008]For this method of using the coordinates of the points on a plane, one can manufacture and maintain a 2D plane pattern easier than the method of using a 3D artificial calibration object. However, this method is very complicated since it must use a plurality of calibration points of the images of a plane and compare them to each other.
 [0009]It is therefore an object of the present invention to provide a system for calibrating internal parameters using projection invariable characteristics of a concentric circle and a method thereof, which is able to use a plane pattern and increase the reliability of calibration, and calculate parameters of a camera without identifying the corresponding points.
 [0010]In order to accomplish the above object, a method of calibrating a camera by calculating camera parameters required to obtain geometric information of an object using projective invariant characteristic of a concentric circle is characterized in that it comprises the steps of; a) taking a plurality of images of a calibration pattern consisting of two or more concentric circles located at a same plane and having different radiuses at different angles to obtain images of ellipses as a result of the projection of said concentric circles; b) calculating the central points of the projected concentric circles using a given algorithm; and c) calculating the principal point and focal length of said camera for tracing the location of a circle in a 3D space using a nonlinear minimizing method based on the central points calculated in said step b).
 [0011]Also, a system for calibrating internal parameters of a camera for calibrating parameters between an actual object and images of the object is characterized in that it comprises a camera; a calibration pattern located on a same plane and consisting of two or more concentric circles having different radiuses, the image of which is taken by said camera; and a controller for calculating a straight line connecting the central points of the images of ellipses obtained by projecting said calibration pattern of concentric circles to said camera at a given angle, and finding crossing points of said straight line and said projected ellipses to obtain the central point coordinate of the concentric circle using a cross ratio (Cr), and tracing the location of the circle in a 3D space based on the coordinate of the central point of the concentric circle to calculate a principal point and focal length of the camera.
 [0012]The aforementioned aspects and other features of the present invention will be explained in the following description, taken in conjunction with the accompanying drawings, wherein:
 [0013][0013]FIG. 1 illustrates an image of a 3D artificial calibration object model taken in a conventional method of calibrating a camera using a 3D calibration object;
 [0014][0014]FIG. 2 is an image of a plane model taken in a conventional method of calibrating a camera using coordinates of points on a twodimensional plane;
 [0015][0015]FIG. 3 is a conceptual view showing that images of a concentric circle pattern is taken at different angles in a method of calibrating camera parameters using a concentric circle pattern according to one embodiment of the present invention;
 [0016][0016]FIG. 4a illustrates two images showing a concentric circle taken at different angles in FIG. 3, and
 [0017][0017]FIG. 4b illustrates images where boundary lines are extracted from the ellipses in FIG. 4a;
 [0018][0018]FIG. 5 is an image showing the central points of the ellipses in FIG. 4b;
 [0019][0019]FIG. 6 is an image showing a straight line connecting the central points of the ellipses in FIG. 5; and
 [0020][0020]FIG. 7 illustrates a position relation of the points on the straight line in FIG. 6.
 [0021]The present invention will be described in detail by way of a preferred embodiment with reference to accompanying drawings, in which like reference numerals are used to identify the same or similar parts.
 [0022][0022]FIG. 3 is a conceptual view showing that images of a concentric circle pattern is taken at different angles in a method of calibrating camera parameters using a concentric circle pattern according to one embodiment of the present invention; FIG. 4a illustrates two images showing a concentric circle taken at different angles in FIG. 3, and FIG. 4b illustrates images where boundary lines are extracted from the ellipses in FIG. 4a; FIG. 5 is an image showing the central points of the ellipses in FIG. 4b; FIG. 6 is an image showing a straight line connecting the central points of the ellipses in FIG. 5; and FIG. 7 illustrates a position relation of the points on the straight line in FIG. 6.
 [0023]As shown in FIG. 3, a concentric circle located in a 3D space has the following geometric conditions.
 [0024]1) Two circles are located at the same coordinate, 2) two circles are located at the same plane, 3) the radiuses of two circles are different. An image of the concentric circle having this geometric structure is taken using a camera.
 [0025]Meanwhile, a concentric circle consisting of a circle having a smaller radius and a circle having a larger radius including the smaller circle has the same normal vector and a vertical distance according to above condition 2).
 [0026]Two images of this concentric circle pattern are taken by one camera at different angles to obtain two sheets of image information. A concentric circle projected at a given angle has a shape of an ellipse, as shown in FIG. 4a to FIG. 5. Generally, the central points of two ellipses do not match. This is because two ellipses obtained by projecting concentric circles at different angles have different shapes according to the angle of projection. The slant of the straight line including the central points of the two ellipses varies depending on the projection angle of the camera. Meanwhile, in the drawings, the central point of an ellipse is a middle point of the straight line connecting two focuses, that is, the point at which the long axis and short axis of the ellipse meet.
 [0027]Meanwhile, as mentioned above, the straight line connecting the central points of the ellipses projected from a 3D concentric circle can be represented by the following Equation 1.
 [0028]First, a projection transformation (T) into a homogenous coordinate system widely used in image processing and projection geometry can be represented into a projection matrix as Equation 1.
$\begin{array}{cc}T=[\text{\hspace{1em}}\ue89e\begin{array}{cccc}{p}_{11}& {p}_{12}& {p}_{13}& {p}_{14}\\ {p}_{21}& {p}_{22}& {p}_{23}& {p}_{24}\\ {p}_{31}& {p}_{32}& {p}_{33}& {p}_{34}\end{array}\ue89e\text{\hspace{1em}}]& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e1\right]\end{array}$  [0029]According to invariable characteristic of a circle which is a trace of points located at a constant distance (radius) from a central point, two concentric circles satisfying the above geometric condition can be represented by following Equation 2.
 (X−x _{0})^{2}+(Y−y _{0})^{2} =r ^{2} , Z=0 [Equation 2]
 [0030]where, x_{0 }and y_{0 }are the central point coordinates of the circle, r is the radius of the circle and Z is a value of the vertical axis with respect to the plane.
 [0031]At this time, using the condition Z=0, a general projection equation is converted into a twodimentional transformation (II) between an image plane and a plane where Z=0 as in Equation 3.
$\begin{array}{cc}\left[\begin{array}{c}s\ue89e\text{\hspace{1em}}\ue89ex\\ s\ue89e\text{\hspace{1em}}\ue89ey\\ s\end{array}\right]=[\text{\hspace{1em}}\ue89e\begin{array}{cccc}{p}_{11}& {p}_{12}& {p}_{13}& {p}_{14}\\ {p}_{21}& {p}_{22}& {p}_{23}& {p}_{24}\\ {p}_{31}& {p}_{32}& {p}_{33}& {p}_{34}\end{array}\ue89e\text{\hspace{1em}}]\ue8a0\left[\begin{array}{c}X\\ Y\\ 0\\ 1\end{array}\right]=[\text{\hspace{1em}}\ue89e\begin{array}{ccc}{p}_{11}& {p}_{12}& {p}_{14}\\ {p}_{21}& {p}_{22}& {p}_{24}\\ {p}_{31}& {p}_{32}& {p}_{34}\end{array}\ue89e\text{\hspace{1em}}]\ue8a0\left[\begin{array}{c}X\\ Y\\ 1\end{array}\right]& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e3\right]\end{array}$  [0032]where, X,Y indicate actual 3D coordinates in projection geometry theory, x, y indicate coordinates in the projected image, s is a scale factor, and p is a value of the projection transformation matrix.
 [0033]Also, rearranging Equation 3 in order to obtain an equation of the projected circle Equation 4 is obtained.
$\begin{array}{cc}s\ue8a0\left[\begin{array}{c}X\\ Y\\ 1\end{array}\right]={[\text{\hspace{1em}}\ue89e\begin{array}{ccc}{p}_{11}& {p}_{12}& {p}_{14}\\ {p}_{21}& {p}_{22}& {p}_{24}\\ {p}_{31}& {p}_{32}& {p}_{34}\end{array}\ue89e\text{\hspace{1em}}]}^{1}\ue8a0\left[\begin{array}{c}x\\ y\\ 1\end{array}\right]=[\text{\hspace{1em}}\ue89e\begin{array}{ccc}{\alpha}_{11}& {\alpha}_{12}& {\alpha}_{13}\\ {\alpha}_{21}& {\alpha}_{22}& {\alpha}_{23}\\ {\alpha}_{31}& {\alpha}_{32}& {\alpha}_{33}\end{array}\ue89e\text{\hspace{1em}}]\ue8a0\left[\begin{array}{c}x\\ y\\ 1\end{array}\right]& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e4\right]\end{array}$  [0034]Considering the homogeneous characteristic of the coordinate system in Equation 4, multiplying terms of Equation 4 are deleted and Equation 2 is substituted to obtain the following Equation 5.
$\begin{array}{cc}{\left(\frac{{\alpha}_{11}\ue89ex+{\alpha}_{12}\ue89ey+{\alpha}_{13}}{{\alpha}_{31}\ue89ex+{\alpha}_{32}\ue89ey+{\alpha}_{33}}{x}_{0}\right)}^{2}+{\left(\frac{{\alpha}_{21}\ue89ex+{\alpha}_{22}\ue89ey+{\alpha}_{23}}{{\alpha}_{31}\ue89ex+{\alpha}_{32}\ue89ey+{\alpha}_{33}}{y}_{0}\right)}^{2}={r}^{2}& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e5\right]\end{array}$  [0035]Meanwhile, Equation 5 may be rearranged to obtain an equation of an ellipse. Then, the coordinate of the central point of the projected ellipse can be obtained by following Equation 6.
$\begin{array}{cc}\left({X}^{\prime},{Y}^{\prime}\right)=\left(\frac{{C}^{\prime}+{r}^{2}\ue89e{D}^{\prime}}{{A}^{\prime}+{r}^{2}\ue89e{B}^{\prime}},\frac{{E}^{\prime}+{r}^{2}\ue89e{F}^{\prime}}{{A}^{\prime}+{r}^{2}\ue89e{B}^{\prime}}\right)& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e6\right]\end{array}$  [0036]Where, A′, B′, C′, D′, E′ and F′ consist of a combinations of factors α_{ij }in Equation 5 and r is a radius constant of the circle.
 [0037]As can be seen from Equation 6, A′, B′, C′, D′, E′ and F′ are independent from r.
 [0038]A straight line connecting the central points of the ellipses is defined as Equation 7 using the coordinates of the central points of the ellipses calculated by Equation 6.
 (B′E′−A′F′)X′+(A′D′−B′C′)Y′+(D′E′+F′C′)=0 [Equation 7]
 [0039]Equation 7 is an equation of the straight line including only constant coefficients with respect to r. This means that the central points of the ellipses obtained as a result of the projection is located on the straight line defined independently from the radius of the circle, where X′, Y′ are coordinates of the centers of the ellipses. As the center of the projected concentric circle is same to the center of the ellipses when r=0 in Equation 2, the center of the circle is also located on the straight line defined by Equation 7.
 [0040]From the relation that the central points of the projected ellipses are located on one straight line, an algorithm to find the coordinate of the central point of the circle will be now described in detail.
 [0041]It is defined that the crossing points of the projected internal ellipse and the straight line are A, B, the crossing points of the external ellipse the straight line are A′, B′, the central point of the circle is O and a given point separated infinitely on the straight line is M_{∞}. Location of these points can be found in FIG. 6 and FIG. 7 and this positional relation is same in 3D and 2D spaces.
 [0042]Meanwhile, a given point infinitely separated from the center of the circle O on the straight line is a point virtually set and is frequently used in projection theory. This point is typically used because it is easy to serve as the boundary value of a condition and it has a characteristic limiting the cross ratio.
 [0043]Therefore, a cross ratio Cr (A, O, B and M.) and Cr (A′, O, B′ and M_{∞}) are same, which is calculated in the same manner in a 3D space by Equation 8.
$\begin{array}{cc}\mathrm{Cr}\ue8a0\left(A,O,B,{M}_{\infty}\right)=\frac{\stackrel{\_}{{\mathrm{BAM}}_{\infty}\ue89eO}}{\stackrel{\_}{{\mathrm{BOM}}_{\infty}\ue89eA}}=2\ue89e\text{}\ue89e\mathrm{Cr}\ue8a0\left({A}^{\prime},O,{B}^{\prime},{M}_{\infty}\right)=\frac{\stackrel{\_}{{B}^{\prime}\ue89e{A}^{\prime}\ue89e{M}_{\infty}\ue89eO}}{\stackrel{\_}{{B}^{\prime}\ue89e{\mathrm{OM}}_{\infty}\ue89e{A}^{\prime}}}=2& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e8\right]\end{array}$  [0044]Where, the points that can be found in the image are A, B, A′ and B′, and Equation 8 has two variables of points O and Mo. These two points are commonly used in the two Equations of Cr (A, O, B and M.) and Cr (A′, O, B and M_{∞}). Therefore, two variables and two equations can be obtained. As a result, the coordinate of point O can be obtained by solving the two Equations.
 [0045]All the items of below Equation 9 can be found from thus obtained coordinate of the central point of the circle. Also, calibration is performed by finding internal parameters u_{0}, v_{0 }and f which minimize the value of Equation 9. Calibration of the internal parameters uses below Equation 9.
$\begin{array}{cc}F\ue8a0\left({u}_{0},{v}_{0},f\right)={\alpha \ue8a0\left(\stackrel{\_}{{n}_{1}}\stackrel{\_}{{n}_{2}}\right)}^{2}+{\beta \ue8a0\left(\stackrel{\_}{{n}_{3}}\stackrel{\_}{{n}_{4}}\right)}^{2}+{\gamma \ue8a0\left({d}_{1}{d}_{2}\right)}^{2}+{\lambda \ue8a0\left({d}_{3}{d}_{4}\right)}^{2}+{\rho}_{1}\ue89e\sum \uf603{R}_{{C}_{1}}\ue8a0\left({n}_{1},{d}_{1}\right){R}_{1}\uf604+{\rho}_{2}\ue89e\sum \uf603{R}_{{C}_{2}}\ue8a0\left({n}_{2},{d}_{2}\right){R}_{2}\uf604+{\rho}_{3}\ue89e\sum \uf603{R}_{{C}_{3}}\ue8a0\left({n}_{3},{d}_{3}\right){R}_{3}\uf604+{\rho}_{4}\ue89e\sum \uf603{R}_{{C}_{4}}\ue8a0\left({n}_{4},{d}_{4}\right){R}_{4}\uf604& \left[\mathrm{Equation}\ue89e\text{\hspace{1em}}\ue89e9\right]\end{array}$  [0046]In Equation 9, u_{0}, v_{0 }and f, respectively indicate a principal point and a focal length of the camera. The first two items (α({overscore (n)}_{1}{overscore (n)}_{2})+β({overscore (n)}_{3}−{overscore (n)}_{4})) indicate that two concentric circles have the same normal vector and next two items (γ(d_{1}−d_{2})+λ(d_{3}−d_{4})) indicate that two concentric circles are located on the plane having the same distance. In addition, next fcur items (ρ_{1}ΣR_{C} _{ l }(n_{l},d_{l})−R_{l}+ρ_{2}ΣR_{C} _{ l }(n_{2},d_{2})−R_{2}+ρ_{3}ΣR_{C} _{ ll }(n_{3},d_{3})−R_{3}+ρ_{4}ΣR_{C} _{ll}(n_{4},d_{4})−R_{4}) indicate that a radius obtained from given internal parameters is same as the actual radius.
 [0047]A concentric circle is drawn on the plane so that the central points of two circles having a radius of 40 mm and 80 mm, respectively are located at one point. Also, as shown in FIG. 4a and FIG. 4b, images of the circle is taken at two different angles to obtain two images. At this time, a XC003 camera manufactured by SONY and with focal radius of 8 mm was used.
 [0048]After obtaining two projection images (containing the boundary lines), a calibration value as shown in Table 1 is obtained by means of a camera calibration method using the concentric circle pattern according to the present invention.
TABLE 1 Class u_{O} v_{O} α_{u} α_{v} Comparative example 306.70 240.44 1104.2 1112.5 (Z. Zhang) Embodiment of the 309.37 246.01 1104.7 1110.7 present invention  [0049]The comparative example (Z.Zhang) in Table 1 shows the calibrated values according to the camera calibration method taught in an article Z. Zhang, “Flexible Camera Calibration by Viewing a Plane From Unknown Orientations”, In Proc. 7th International Conference Computer Vision, Corfu, Greece, pp.666673, September 1999. Table 1 shows that the calibration value calculated using a concentric circle pattern is similar to the calibration value in the comparative example of Z. Zhang.
 [0050]As mentioned above, the method of calibrating internal parameters of a camera using a concentric circle pattern according to the present invention provides the calibration parameters of the camera from images of a 2D concentric circle taken from different angles. Thus, the present invention has an advantage that a concentric circle pattern which is easy to manufacture and maintain.
 [0051]Also, the method of calibrating internal parameters of a camera using a concentric circle pattern according to the present invention calculates camera calibration parameters using a plurality of Equations. Therefore, the present invention has an advantage that it can easily calibrate a camera by calculating the calibration parameters.
 [0052]The present invention has been described with reference to a particular embodiment in connection with a particular application. Those having ordinary skill in the art and access to the teachings of the present invention will recognize additional modifications and applications within the scope thereof. It is therefore intended by the appended claims to cover any and all such applications, modifications, and embodiments within the scope of the present invention.
Claims (8)
 1. A method of calibrating a camera by calculating camera parameters required to obtain geometric information of an object using projection invariable characteristic of a concentric circle comprising the steps of;a) taking a plurality of images of a calibration pattern consisting of two or more concentric circles located at a same plane and having different radiuses at different angles to obtain images of ellipses as a result of the projection of said concentric circles;b) calculating the central points of the projected concentric circles using a given algorithm; andc) calculating the principal point and focal length of said camera for tracing the location of a circle in a 3D space using a nonlinear minimizing method based on the central points calculated in said step b).
 2. The method of calibrating a camera according to
claim 1 , wherein said concentric circle includes two circles and said step b) includes the steps of:d) calculating the central points of the two ellipses to obtain a straight line passing through said central points;e) calculating crossing points of said straight line and said two ellipses; andf) calculating the coordinate of the central point of said concentric circle using a cross ratio depending on the central point (0) of the concentric circle, a given point (M_{∞}) infinitely separated from the central point of the concentric circle on the straight line passing through the central points of the ellipses, and the crossing points of the straight line and the ellipses.  3. The method of calibrating a camera according to
claim 2 , wherein said step d) calculates the straight line passing through the centers of the ellipses by Equation 1:(B′E′−A′F′)X′+(A′D′−B′C′)Y′+(D′E′+F′C′)=0 (1)(where A′, B′, C′, D′, E′ and F′ are constants with respect to the radius of the circle and, X′,Y′ are the coordinates of the ellipse center)  4. The method of calibrating a camera according to
claim 3 , wherein said step f) calculates the coordinate of the center of the circle by Equation 2:$\begin{array}{cc}\mathrm{Cr}\ue8a0\left(A,O,B,{M}_{\infty}\right)=\frac{\stackrel{\_}{{\mathrm{BAM}}_{\infty}\ue89eO}}{\stackrel{\_}{{\mathrm{BOM}}_{\infty}\ue89eA}}=2\ue89e\text{}\ue89e\mathrm{Cr}\ue8a0\left({A}^{\prime},O,{B}^{\prime},{M}_{\infty}\right)=\frac{\stackrel{\_}{{B}^{\prime}\ue89e{A}^{\prime}\ue89e{M}_{\infty}\ue89eO}}{\stackrel{\_}{{B}^{\prime}\ue89e{\mathrm{OM}}_{\infty}\ue89e{A}^{\prime}}}=2& \left(2\right)\end{array}$ where A,B are the coordinates of crossing points of an ellipse generated by projection of smaller circle and the straight line, and A′,B′ are the coordinates of crossing points of an ellipse generated by projection of an external circle surrounding the smaller circle and the straight line)  5. The method of calibrating a camera according to
claim 3 , wherein said step c) includes the steps of:g) assuming that the central point and a approximate focal length obtained in said step b) as a principal point focal length;h) estimating the location of the circle in a 3D space using said principle point and focal length; andi) repeating said steps g) and h) using a given nonlinear minimization algorithm using a normal vector of said concentric circle, a vertical distance from an original point and the radius of the concentric circle based on the location of the circle obtained in said step h), so that below Equation 3 is minimized.F(u _{0},ν_{0} ,f)=α({overscore (n)} _{1} −{overscore (n)} _{2})^{2}+β({overscore (n)} _{3} −{overscore (n)} _{4})^{2}+γ(d _{1} −d _{2})^{2}+λ(d _{3} −d _{4})^{2}+ρ_{1} ΣR _{C} _{ 1 }(n _{1} ,d _{1})−R _{1}+ρ_{2} ΣR _{C} _{ 2 }(n _{2} ,d _{2})−R_{2}+ρ_{3} ΣR _{C} _{ 1 }(n _{3} ,d _{3})−R _{3}+ρ_{4}ρ_{4} ΣR _{C} _{ 4 }(n _{4} ,d _{4})−R _{4} (3)where u_{0}, v_{0 }and f indicate a principle point and a focal length of the camera, the first two items in Equation 3 indicate that two concentric circles have the same normal vector, next two items in Equation 3 indicate that two concentric circles are located on a plane having the same length, and next four items indicate that a radius obtained from given internal parameters is same as the actual radius.  6. A system for calibrating internal parameters of a camera for calibrating parameters between an actual object and images of the object, comprising;a camera;a calibration pattern located on a same plane and consisting of two or more concentric circles having different radiuses, the image of which is taken by said camera; anda controller for calculating a straight line connecting the central points of the images of ellipses obtained by projecting said calibration pattern of concentric circles to said camera at a given angle, and finding crossing points of said straight line and said projected ellipses to obtain the central point coordinate of the concentric circle using a cross ratio (Cr), and tracing the location of the circle in a 3D space based on the coordinate of the central point of the concentric circle to calculate a principal point and focal length of the camera.
 7. The system for calibrating internal parameters of a camera according to
claim 6 , wherein said principal point and focal length are obtained using a given nonlinear minimization algorithm so that below Equation 3 is minimized.F(u _{0},ν_{0} ,f)=α({overscore (n)} _{1} −{overscore (n)} _{2})^{2}+β({overscore (n)} _{3} −{overscore (n)} _{4})^{2}+γ(d _{1} −d _{2})^{2}+λ(d _{3} −d _{4})^{2}+ρ_{1} ΣR _{C} _{ 1 }(n _{1} ,d _{1})−R _{1}+ρ_{2} ΣR _{C} _{ 2 }(n _{2} ,d _{2})−R_{2}+ρ_{3} ΣR _{C} _{ 1 }(n _{3} ,d _{3})−R _{3}+ρ_{4}ρ_{4} ΣR _{C} _{ 4 }(n _{4} ,d _{4})−R _{4} (3)where u_{0}, v_{0 }and f indicate a principle point and a focal length of the camera, the first two items in Equation 3 indicate that two concentric circles have the same normal vector, next two items in Equation 3 indicate that two concentric circles are located on a plane having the same length, and next four items indicate that a radius obtained from given internal parameters is same as the actual radius.  8. The system for calibrating internal parameters of a camera according to
claim 6 , wherein said camera is a CCD camera.
Priority Applications (2)
Application Number  Priority Date  Filing Date  Title 

KR20010017418A KR100386090B1 (en)  20010402  20010402  Camera calibration system and method using planar concentric circles 
KR200117418  20010402 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US11074562 US7155030B2 (en)  20010402  20050308  Camera calibration system using planar concentric circles and method thereof 
Publications (1)
Publication Number  Publication Date 

US20020186897A1 true true US20020186897A1 (en)  20021212 
Family
ID=19707748
Family Applications (2)
Application Number  Title  Priority Date  Filing Date 

US09919086 Abandoned US20020186897A1 (en)  20010402  20010731  Camera calibration system using planar concentric circles and method thereof 
US11074562 Expired  Fee Related US7155030B2 (en)  20010402  20050308  Camera calibration system using planar concentric circles and method thereof 
Family Applications After (1)
Application Number  Title  Priority Date  Filing Date 

US11074562 Expired  Fee Related US7155030B2 (en)  20010402  20050308  Camera calibration system using planar concentric circles and method thereof 
Country Status (2)
Country  Link 

US (2)  US20020186897A1 (en) 
KR (1)  KR100386090B1 (en) 
Cited By (21)
Publication number  Priority date  Publication date  Assignee  Title 

US20050228555A1 (en) *  20030820  20051013  Samsung Electronics Co., Ltd.  Method of constructing artificial mark for autonomous driving, apparatus and method of determining position of intelligent system using artificial mark and intelligent system employing the same 
US20090066643A1 (en) *  20070907  20090312  Samsung Electronics Co., Ltd  Touch screen panel to input multidimension values and method for controlling touch screen panel 
WO2009056050A1 (en) *  20071024  20090507  Shenzhen Huawei Communication Technologies Co. , Ltd.  Video camera calibration method and device thereof 
CN101303768B (en)  20080617  20100623  东南大学  Method for correcting circle center error of circular index point when translating camera perspective projection 
CN101813465A (en) *  20100330  20100825  中北大学  Monocular vision measuring method of noncontact precision measuring corner 
CN102622744A (en) *  20120117  20120801  浙江大学  Telephoto camera calibration method based on polynomial projection model 
CN102853903A (en) *  20110628  20130102  中国科学院西安光学精密机械研究所  Method for demarcating linearity indexes of scientific grade Chargecoupled Device (CCD) 
US20130033596A1 (en) *  20110803  20130207  The Boeing Company  Projection aided feature measurement using uncalibrated camera 
CN102103746B (en)  20091221  20130213  云南大学  Method for calibrating parameters in camera through solving circular ring points by utilizing regular tetrahedron 
CN102982550A (en) *  20121214  20130320  云南大学  Method of solving intrinsic parameters of camera with regular pentagonal prismatic table 
CN102982549A (en) *  20121214  20130320  云南大学  Method of solving intrinsic parameters of camera with two centersharing and principal axis orthonormal intersected identical ellipses 
CN103592285A (en) *  20131126  20140219  长春新产业光电技术有限公司  Nonlinear calibration method of Raman spectrum wavelength and CCD picture element 
CN104155003A (en) *  20140623  20141119  中国科学院光电研究院  Highstability rotating mirror interferometer 
CN104165692A (en) *  20140725  20141126  华东师范大学  Portable spectrum rapid detector based on quantum effect photodetector 
CN104165693A (en) *  20140729  20141126  中国科学院西安光学精密机械研究所  Largeaperture staticstate interference spectrum imaging instrument interferogram optical splicing method 
CN104331900A (en) *  20141125  20150204  湖南科技大学  Corner subpixel positioning method in CCD (charge coupled device) camera calibration 
CN104380064A (en) *  20120615  20150225  株式会社日立高新技术  Optical signal detecting circuit, light volume detecting device, and charged particle beam device 
CN104848944A (en) *  20150528  20150819  上海理工大学  Z scanning optical polarization measurement device and measurement method thereof 
CN105241427A (en) *  20151013  20160113  中北大学  Singlecamera vision measurement method using concentric circle target for measuring position posture and rotating angle of object 
CN106017669A (en) *  20160513  20161012  中国科学院紫金山天文台  Multifunctional reading circuit system for KID detector array 
CN106017667A (en) *  20160719  20161012  苏州市职业大学  Multifunctional ultraviolet ray detector powered by solar energy 
Families Citing this family (21)
Publication number  Priority date  Publication date  Assignee  Title 

JP4147059B2 (en) *  20020703  20080910  株式会社トプコン  Calibration data measuring device, measuring method and a program, and a computerreadable recording medium, the image data processing device 
JP4598372B2 (en) *  20030526  20101215  株式会社トプコン  Threedimensional field for calibration, method of photographing a threedimensional field for calibration 
KR100690857B1 (en) *  20051202  20070227  엘지전자 주식회사  Apparatus for controlling camera alignment of mobile communications terminal 
KR100793838B1 (en)  20060927  20080111  한국전자통신연구원  Appratus for findinng the motion of camera, system and method for supporting augmented reality in ocean scene using the appratus 
US8238648B2 (en) *  20061204  20120807  Electronics And Telecommunications Research Institute  Apparatus and method for estimating focal length of camera 
KR100966592B1 (en)  20071217  20100629  한국전자통신연구원  Method for calibrating a camera with homography of imaged parallelogram 
KR101694294B1 (en)  20101222  20170109  한국전자통신연구원  Apparatus and method for calibrating camera with imaged parallelograms and multiview constraints 
US8743214B2 (en) *  20110511  20140603  Intel Corporation  Display screen for camera calibration 
US8872897B2 (en)  20110511  20141028  Intel Corporation  Camera calibration using an easily produced 3D calibration pattern 
KR101638173B1 (en) *  20110906  20160712  한국전자통신연구원  Method and apparatus for providing automated detection of calibration 
US9841367B2 (en)  20110916  20171212  The University Of North Carolina At Charlotte  Methods and devices for optical sorting of microspheres based on their resonant optical properties 
CN102930551B (en) *  20121122  20150506  云南大学  Camera intrinsic parameters determined by utilizing projected coordinate and epipolar line of centres of circles 
US9519968B2 (en)  20121213  20161213  HewlettPackard Development Company, L.P.  Calibrating visual sensors using homography operators 
CN103116892A (en) *  20130315  20130522  云南大学  Linear solving of video camera internal parameters by utilizing two intersected same circles and common tangent 
CN103247048A (en) *  20130510  20130814  东南大学  Camera mixing calibration method based on quadratic curve and straight lines 
CA2819956A1 (en)  20130702  20150102  Guy Martin  High accuracy camera modelling and calibration method 
JP6217227B2 (en) *  20130812  20171025  株式会社リコー  Calibration apparatus, method, and program 
CN104637042A (en) *  20131106  20150520  北京计算机技术及应用研究所  Camera calibration method based on circular reference object and system thereof 
CN104318555B (en) *  20141013  20170208  江南大学  In the center of the target image projected point precise positioning method 
KR20160098543A (en)  20150209  20160819  한국전자통신연구원  Device and method for multiview image calibration 
CN104794704B (en) *  20150327  20171117  华为技术有限公司  Calibration Method template, detection method, and a terminal device 
Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US5559334A (en) *  19950522  19960924  General Electric Company  Epipolar reconstruction of 3D structures 
US5918196A (en) *  19961129  19990629  Cognex Corporation  Vision system for analyzing solidofrevolution radius profile 
US6088098A (en) *  19980117  20000711  Robert Bosch Gmbh  Calibration method for a laserbased splitbeam method 
US6195455B1 (en) *  19980701  20010227  Intel Corporation  Imaging device orientation information through analysis of test images 
US6437823B1 (en) *  19990430  20020820  Microsoft Corporation  Method and system for calibrating digital cameras 
US6671400B1 (en) *  20000928  20031230  Tateyama R & D Co., Ltd.  Panoramic image navigation system using neural network for correction of image distortion 
US6750860B1 (en) *  19981228  20040615  Microsoft Corporation  Rendering with concentric mosaics 
Family Cites Families (4)
Publication number  Priority date  Publication date  Assignee  Title 

US5054907A (en) *  19891222  19911008  Phoenix Laser Systems, Inc.  Ophthalmic diagnostic apparatus and method 
US5886767A (en) *  19961009  19990323  Snook; Richard K.  Keratometry system and method for measuring physical parameters of the cornea 
US6995762B1 (en) *  20010913  20060207  Symbol Technologies, Inc.  Measurement of dimensions of solid objects from twodimensional image(s) 
JP3624288B2 (en) *  20010917  20050302  株式会社日立製作所  Store management system 
Patent Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US5559334A (en) *  19950522  19960924  General Electric Company  Epipolar reconstruction of 3D structures 
US5918196A (en) *  19961129  19990629  Cognex Corporation  Vision system for analyzing solidofrevolution radius profile 
US6088098A (en) *  19980117  20000711  Robert Bosch Gmbh  Calibration method for a laserbased splitbeam method 
US6195455B1 (en) *  19980701  20010227  Intel Corporation  Imaging device orientation information through analysis of test images 
US6750860B1 (en) *  19981228  20040615  Microsoft Corporation  Rendering with concentric mosaics 
US6437823B1 (en) *  19990430  20020820  Microsoft Corporation  Method and system for calibrating digital cameras 
US6671400B1 (en) *  20000928  20031230  Tateyama R & D Co., Ltd.  Panoramic image navigation system using neural network for correction of image distortion 
Cited By (27)
Publication number  Priority date  Publication date  Assignee  Title 

US7437226B2 (en) *  20030820  20081014  Samsung Electronics Co., Ltd.  Method of constructing artificial mark for autonomous driving, apparatus and method of determining position of intelligent system using artificial mark and intelligent system employing the same 
US20050228555A1 (en) *  20030820  20051013  Samsung Electronics Co., Ltd.  Method of constructing artificial mark for autonomous driving, apparatus and method of determining position of intelligent system using artificial mark and intelligent system employing the same 
US8229619B2 (en)  20030820  20120724  Samsung Electronics Co., Ltd.  Method of constructing artificial mark for autonomous driving, apparatus and method of determining position of intelligent system using artificial mark and intelligent system employing the same 
US20090066643A1 (en) *  20070907  20090312  Samsung Electronics Co., Ltd  Touch screen panel to input multidimension values and method for controlling touch screen panel 
US8812992B2 (en) *  20070907  20140819  Samsung Electronics Co., Ltd.  Touch screen panel to input multidimension values and method for controlling touch screen panel 
US8436904B2 (en)  20071024  20130507  Huawei Device Co., Ltd.  Method and apparatus for calibrating video camera 
WO2009056050A1 (en) *  20071024  20090507  Shenzhen Huawei Communication Technologies Co. , Ltd.  Video camera calibration method and device thereof 
US20100259624A1 (en) *  20071024  20101014  Kai Li  Method and apparatus for calibrating video camera 
CN101303768B (en)  20080617  20100623  东南大学  Method for correcting circle center error of circular index point when translating camera perspective projection 
CN102103746B (en)  20091221  20130213  云南大学  Method for calibrating parameters in camera through solving circular ring points by utilizing regular tetrahedron 
CN101813465A (en) *  20100330  20100825  中北大学  Monocular vision measuring method of noncontact precision measuring corner 
CN102853903A (en) *  20110628  20130102  中国科学院西安光学精密机械研究所  Method for demarcating linearity indexes of scientific grade Chargecoupled Device (CCD) 
US8922647B2 (en) *  20110803  20141230  The Boeing Company  Projection aided feature measurement using uncalibrated camera 
US20130033596A1 (en) *  20110803  20130207  The Boeing Company  Projection aided feature measurement using uncalibrated camera 
CN102622744A (en) *  20120117  20120801  浙江大学  Telephoto camera calibration method based on polynomial projection model 
CN104380064A (en) *  20120615  20150225  株式会社日立高新技术  Optical signal detecting circuit, light volume detecting device, and charged particle beam device 
CN102982550A (en) *  20121214  20130320  云南大学  Method of solving intrinsic parameters of camera with regular pentagonal prismatic table 
CN102982549A (en) *  20121214  20130320  云南大学  Method of solving intrinsic parameters of camera with two centersharing and principal axis orthonormal intersected identical ellipses 
CN103592285A (en) *  20131126  20140219  长春新产业光电技术有限公司  Nonlinear calibration method of Raman spectrum wavelength and CCD picture element 
CN104155003A (en) *  20140623  20141119  中国科学院光电研究院  Highstability rotating mirror interferometer 
CN104165692A (en) *  20140725  20141126  华东师范大学  Portable spectrum rapid detector based on quantum effect photodetector 
CN104165693A (en) *  20140729  20141126  中国科学院西安光学精密机械研究所  Largeaperture staticstate interference spectrum imaging instrument interferogram optical splicing method 
CN104331900A (en) *  20141125  20150204  湖南科技大学  Corner subpixel positioning method in CCD (charge coupled device) camera calibration 
CN104848944A (en) *  20150528  20150819  上海理工大学  Z scanning optical polarization measurement device and measurement method thereof 
CN105241427A (en) *  20151013  20160113  中北大学  Singlecamera vision measurement method using concentric circle target for measuring position posture and rotating angle of object 
CN106017669A (en) *  20160513  20161012  中国科学院紫金山天文台  Multifunctional reading circuit system for KID detector array 
CN106017667A (en) *  20160719  20161012  苏州市职业大学  Multifunctional ultraviolet ray detector powered by solar energy 
Also Published As
Publication number  Publication date  Type 

US7155030B2 (en)  20061226  grant 
KR20020077685A (en)  20021014  application 
KR100386090B1 (en)  20030602  grant 
US20050207640A1 (en)  20050922  application 
Similar Documents
Publication  Publication Date  Title 

US6714679B1 (en)  Boundary analyzer  
US6621921B1 (en)  Image processing apparatus  
US5396331A (en)  Method for executing threedimensional measurement utilizing correctively computing the absolute positions of CCD cameras when image data vary  
Stein et al.  Modelbased brightness constraints: On direct estimation of structure and motion  
Fusiello et al.  Quasieuclidean uncalibrated epipolar rectification  
US20090118864A1 (en)  Method and system for finding a tool center point for a robot using an external camera  
US5995639A (en)  Apparatus for identifying person  
US6915008B2 (en)  Method and apparatus for multinodal, threedimensional imaging  
US4893183A (en)  Robotic vision system  
US6985075B2 (en)  Obstacle detection apparatus and method  
US20070127816A1 (en)  Computervision system for classification and spatial localization of bounded 3dobjects  
Smith et al.  Correction of distortion in endoscope images  
US20030144813A1 (en)  Position and orientation determination method and apparatus and storage medium  
US5495535A (en)  Method of inspecting articles  
Wheeler et al.  Consensus surfaces for modeling 3D objects from multiple range images  
US20030144765A1 (en)  Method and apparatus for single camera 3D vision guided robotics  
US6278798B1 (en)  Image object recognition system and method  
US6363168B1 (en)  Measurement position determination on a semiconductor wafer  
US6516099B1 (en)  Image processing apparatus  
Gruen  Adaptive least squares correlation: a powerful image matching technique  
US5926581A (en)  System for topographic mapping from remotely sensed images  
US6751338B1 (en)  System and method of using range image data with machine vision tools  
US20060257051A1 (en)  Method of registering and aligning multiple images  
McIvor et al.  A comparison of local surface geometry estimation methods  
US7313252B2 (en)  Method and system for improving video metadata through the use of frametoframe correspondences 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: KOREA ADVANCED INSTITUTE OF SCIENCE AND TECHNOLOGY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KIM, JUNSIK;KWEON, INSO;REEL/FRAME:012041/0429 Effective date: 20010716 