CN102509324A - Rotational stereovision rotary axis determining method based on quadratic curve fitting - Google Patents
Rotational stereovision rotary axis determining method based on quadratic curve fitting Download PDFInfo
- Publication number
- CN102509324A CN102509324A CN201110338048XA CN201110338048A CN102509324A CN 102509324 A CN102509324 A CN 102509324A CN 201110338048X A CN201110338048X A CN 201110338048XA CN 201110338048 A CN201110338048 A CN 201110338048A CN 102509324 A CN102509324 A CN 102509324A
- Authority
- CN
- China
- Prior art keywords
- point
- turning axle
- image
- stereovision
- projection
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
A rotational stereovision rotary axis determining method based on quadratic curve fitting includes steps: firstly, acquiring trace of points on an image; secondly, acquiring some discrete points on a projection image of a spatial circle formed by a certain rotating point and obtaining a corresponding quadratic curve by fitting the discrete points; thirdly, solving an image of a line at infinity; fourthly, defining the projection of a spatial circle formed by a certain spatial rotating point as a quadratic curve Ci, accordingly the corresponding center oi=Cil infinity, being the projection of a certain point on the rotary axis; and fifthly, resetting up the center point to obtain a spatial point on the rotary axis after the center point of the quadratic curve is obtained; and sixthly, obtaining the rotational stereovision rotary axis by fitting the points on the rotary axis. The rotational stereovision rotary axis determining method based on quadratic curve fitting is capable of simplifying calculation and improving precision and has high practicality.
Description
Technical field
The present invention relates to fields such as computer vision, Flame Image Process, graphics, how much, mathematics, especially a kind of method of rotating stereovision.
Background technology
Turning axle commonly used at present confirms that method is a discrete point coordinate position on the left track when perhaps some put the rotation on universal stage through the demarcation thing; Obtain space tracking by discrete point then, after calculate the character calculating turning axle of confirming that the center of circle is perhaps intersected through Surface Method vector and turning axle.The technological deficiency that existing definite method exists: calculation of complex, precision is lower, practicality is relatively poor.
Summary of the invention
For the deficiency that calculation of complex, precision are lower, practicality is relatively poor of definite method of overcoming existing rotating stereovision turning axle, the present invention provides a kind of and simplifies calculating, promote precision, have the definite method based on the rotating stereovision turning axle of conic fitting of good practicality.
The technical solution adopted for the present invention to solve the technical problems is:
A kind of definite method of the rotating stereovision turning axle based on conic fitting, definite method of said rotating stereovision turning axle may further comprise the steps:
(1) obtains the tracking of putting on the image;
(2) after some discrete points on the projected image that has obtained the formed space circle of certain some rotation, obtain corresponding quafric curve according to said discrete point match, the equation of quafric curve is:
ax
2+by
2+cx+dy+e=0
Wherein a, b, c, d and e are the coefficients of quafric curve, are written as matrix form to be:
As n discrete point, then matrix of coefficients is n * 5 matrixes, goes out quafric curve with least square fitting;
(3) ask the image of line at infinity: focoid are on line at infinity, and the space focoid are [1,0, ± i, 0]
T, after the projection matrix projection, can form two points is c
iAnd c
j, obtain the image l of line at infinity by these two points
∞=c
i* c
j
(4) make space point rotate formed space circle and be projected as quafric curve C
i, its corresponding center o then
i=C
il
∞, promptly be the projection of certain point on the turning axle;
(5) after obtaining the central point of quafric curve, rebuild said central point and promptly obtain the spatial point on the turning axle;
Get wherein two width of cloth images and promptly can rebuild certain point on the turning axle, suppose that certain unknown point is X on the turning axle
w, its projection calculates through said method and is o
i, the projection matrix that wherein two width of cloth images are corresponding is P
1And P
2, then obtain by linear camera model:
o
i=P
1X
w
o
i=P
2X
w
Get a spatial point on the turning axle by top solution of equation;
(6) according to the point on the turning axle, match obtains the rotating stereovision turning axle.
Further, in the said step (3), be the image of line at infinity with respect to the hachure that disappears that the limit of some imaging planes simulates this Plane of rotation through every other image.
Beneficial effect of the present invention mainly shows: simplify calculating, promote precision, have good practicality.
Description of drawings
Fig. 1 is the synoptic diagram of each limit of relative first width of cloth image.
Fig. 2 is the parameter synoptic diagram of discrete point rotational trajectory.
Fig. 3 is the synoptic diagram of line at infinity.
Fig. 4 is the turning axle synoptic diagram of rebuilding.
Embodiment
Below in conjunction with accompanying drawing the present invention is further described.
With reference to Fig. 1~Fig. 4, a kind of definite method of the rotating stereovision turning axle based on conic fitting, definite method of said rotating stereovision turning axle may further comprise the steps:
(1) obtains the tracking of putting on the image;
(2) after some discrete points on the projected image that has obtained the formed space circle of certain some rotation, obtain corresponding quafric curve according to said discrete point match, the equation of quafric curve is:
ax
2+by
2+cx+dy+e=0
Wherein a, b, c, d and e are the coefficients of quafric curve, are written as matrix form to be:
As n discrete point, then matrix of coefficients is n * 5 matrixes, goes out quafric curve with least square fitting;
(3) ask the image of line at infinity: focoid are on line at infinity, and the space focoid are [1,0, ± i, 0]
T, after the projection matrix projection, can form two points is c
iAnd c
j, obtain the image l of line at infinity by these two points
∞=c
i* c
j
(4) make space point rotate formed space circle and be projected as quafric curve C
i, its corresponding center o then
i=C
il
∞, promptly be the projection of certain point on the turning axle;
(5) after obtaining the central point of quafric curve, rebuild said central point and promptly obtain the spatial point on the turning axle;
Get wherein two width of cloth images and promptly can rebuild certain point on the turning axle, suppose that certain unknown point is X on the turning axle
w, its projection calculates through said method and is o
i, the projection matrix that wherein two width of cloth images are corresponding is P
1And P
2, then obtain by linear camera model:
o
i=P
1X
w
o
i=P
2X
w
Get a spatial point on the turning axle by top solution of equation;
(6) according to the point on the turning axle, match obtains the rotating stereovision turning axle.
In the present embodiment; Under the constant situation of the intrinsic parameter of video camera; When rotating with rotation platform, the rotary motion trace of any point on the rotation platform on the object (removing the point on the turning axle) is a space circle, and such space circle will become a quafric curve through projection.Confirm that based on the turning axle of conic fitting the method concrete steps are following:
(1) obtains the tracking of certain point.The tracking of point can realize automatically, but for for the purpose of accurately, we by hand method obtain the tracking of certain point.Also possibly have some error owing to manually confirm certain point on the image, can be auxiliary with program, promptly searching for angle point in the scope around the determined point of mouse, the some tracking meeting that obtains like this is more accurate.
(2) match quafric curve.After some discrete points on the projected image that has obtained the formed space circle of certain some rotation, obtain corresponding quafric curve according to these discrete point matches, the equation of quafric curve is:
ax
2+by
2+cx+dy+e=0
Being written as matrix form is:
As n discrete point, then matrix of coefficients is n * 5 matrixes, can go out quafric curve with least square fitting.
(3) ask the image of line at infinity.The most directly and simply asking the method for line at infinity image for the image of having demarcated is that image by focoid calculates.Focoid are on line at infinity, and the space focoid are [1,0, ± i, 0]
T, after the projection matrix projection, can form two points is c
iAnd c
j, obtain the image l of line at infinity by these two points
∞=c
i* c
j
The another kind of method of asking the line at infinity image is to find the solution through each limit.Under the single shaft rotation; Every other image with respect to the limit of some imaging planes all on same straight line; And this straight line is exactly the hachure that disappears of this Plane of rotation, so another kind of method then is to be the image of line at infinity through every other image with respect to the hachure that disappears that the limit of some imaging planes simulates this Plane of rotation.As shown in Figure 1, the center of all video cameras on same imaging plane, o
i, o
jAnd o
kProjection on first width of cloth image is limit e
I1, e
J1And e
K1, all on same straight line L, L promptly is the image l of line at infinity to these limits
∞
(4) center of calculating quafric curve, the i.e. image of certain point on the turning axle.Make space point rotate formed space circle and be projected as quafric curve C
i, its corresponding center o then
i=C
il
∞, promptly be the projection of certain point on the turning axle.
(5) point on the turning axle of reconstruction space.After obtaining the central point of some quafric curves, rebuild these central points and promptly obtain the spatial point on the turning axle.The image of turning axle is a not moving-wire in all images, so the coordinate of the image of certain point in any image all is identical on the turning axle that we ask, gets wherein two width of cloth images and promptly can rebuild certain point on the turning axle.Suppose that certain unknown point is X on the turning axle
w, its projection calculates through said method and is o
i, the projection matrix that wherein two width of cloth images are corresponding is P
1And P
2, then have linear camera model to obtain:
o
i=P
1X
w
o
i=P
2X
w
Can solve a spatial point on the turning axle by top equation.
(6) match space turning axle.As long as the point on enough turning axles is arranged, the space turning axle promptly can be tried to achieve in match.
Claims (2)
1. based on definite method of the rotating stereovision turning axle of conic fitting, it is characterized in that:
Definite method of said rotating stereovision turning axle may further comprise the steps:
(1) obtains the tracking of putting on the image;
(2) after some discrete points on the projected image that has obtained the formed space circle of certain some rotation, obtain corresponding quafric curve according to said discrete point match, the equation of quafric curve is:
ax
2+by
2+cx+dy+e=0
Wherein a, b, c, d and e are the coefficients of quafric curve, are written as matrix form to be:
As n discrete point, then matrix of coefficients is n * 5 matrixes, goes out quafric curve with least square fitting;
(3) ask the image of line at infinity: focoid are on line at infinity, and the space focoid are [1,0, ± i, 0]
T, after the projection matrix projection, can form two points is c
iAnd c
j, obtain the image l of line at infinity by these two points
∞=c
i* c
j
(4) make space point rotate formed space circle and be projected as quafric curve C
i, its corresponding center o then
i=C
il
∞, promptly be the projection of certain point on the turning axle;
(5) after obtaining the central point of quafric curve, rebuild said central point and promptly obtain the spatial point on the turning axle;
Get wherein two width of cloth images and promptly can rebuild certain point on the turning axle, suppose that certain unknown point is X on the turning axle
w, its projection calculates through said method and is o
i, the projection matrix that wherein two width of cloth images are corresponding is P
1And P
2, then obtain by linear camera model:
o
i=P
1X
w
o
i=P
2X
w
Get a spatial point on the turning axle by top solution of equation;
(6) according to the point on the turning axle, match obtains the rotating stereovision turning axle.
2. definite method of the rotating stereovision turning axle based on conic fitting as claimed in claim 1; It is characterized in that: in the said step (3), be the image of line at infinity with respect to the hachure that disappears that the limit of some imaging planes simulates this Plane of rotation through every other image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110338048XA CN102509324A (en) | 2011-10-31 | 2011-10-31 | Rotational stereovision rotary axis determining method based on quadratic curve fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110338048XA CN102509324A (en) | 2011-10-31 | 2011-10-31 | Rotational stereovision rotary axis determining method based on quadratic curve fitting |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102509324A true CN102509324A (en) | 2012-06-20 |
Family
ID=46221402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110338048XA Pending CN102509324A (en) | 2011-10-31 | 2011-10-31 | Rotational stereovision rotary axis determining method based on quadratic curve fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102509324A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663742A (en) * | 2012-03-22 | 2012-09-12 | 浙江工业大学 | Determination method of rotary stereo visual rotation axis based on quadratic curve fitting |
CN103072302A (en) * | 2013-01-10 | 2013-05-01 | 浙江工业大学 | Braking curve self-learning method for numerical control press |
CN107244165A (en) * | 2017-06-09 | 2017-10-13 | 浙江新盛蓝科技有限公司 | A kind of compass teaching aid and its drawing implementation method |
WO2018076211A1 (en) * | 2016-10-26 | 2018-05-03 | 中国科学院自动化研究所 | Method for quadratic curve fitting in image based on geometric error optimization |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6606404B1 (en) * | 1999-06-19 | 2003-08-12 | Microsoft Corporation | System and method for computing rectifying homographies for stereo vision processing of three dimensional objects |
CN101320474A (en) * | 2008-06-25 | 2008-12-10 | 浙江工业大学 | Exterior parameter self-calibration method for camera with rotating stereovision |
CN101329764A (en) * | 2008-07-31 | 2008-12-24 | 上海交通大学 | Method for positioning video camera using two arbitrary coplane circles |
-
2011
- 2011-10-31 CN CN201110338048XA patent/CN102509324A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6606404B1 (en) * | 1999-06-19 | 2003-08-12 | Microsoft Corporation | System and method for computing rectifying homographies for stereo vision processing of three dimensional objects |
CN101320474A (en) * | 2008-06-25 | 2008-12-10 | 浙江工业大学 | Exterior parameter self-calibration method for camera with rotating stereovision |
CN101329764A (en) * | 2008-07-31 | 2008-12-24 | 上海交通大学 | Method for positioning video camera using two arbitrary coplane circles |
Non-Patent Citations (2)
Title |
---|
BEILI GUO ET AL: "a self-calibration method for rotational stereo vision", 《2009 INTERNATIONAL CONFERENCE ON ARTIFICIAL INTELLIGENCE AND COMPUTATIONAL INTELLIGENCE》, 8 November 2009 (2009-11-08), pages 128 - 132, XP031598326 * |
GUAND JIANG ET AL: "Geometry of single axis motions using conic fitting", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 》, vol. 25, no. 10, 31 October 2003 (2003-10-31), pages 1343 - 1348 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663742A (en) * | 2012-03-22 | 2012-09-12 | 浙江工业大学 | Determination method of rotary stereo visual rotation axis based on quadratic curve fitting |
CN103072302A (en) * | 2013-01-10 | 2013-05-01 | 浙江工业大学 | Braking curve self-learning method for numerical control press |
CN103072302B (en) * | 2013-01-10 | 2015-06-03 | 浙江工业大学 | Braking curve self-learning method for numerical control press |
WO2018076211A1 (en) * | 2016-10-26 | 2018-05-03 | 中国科学院自动化研究所 | Method for quadratic curve fitting in image based on geometric error optimization |
CN107244165A (en) * | 2017-06-09 | 2017-10-13 | 浙江新盛蓝科技有限公司 | A kind of compass teaching aid and its drawing implementation method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108053450B (en) | High-precision binocular camera calibration method based on multiple constraints | |
Pereira et al. | Self calibration of multiple LIDARs and cameras on autonomous vehicles | |
CN100562707C (en) | Binocular vision rotating axis calibration method | |
CN107146254A (en) | The Camera extrinsic number scaling method of multicamera system | |
CN103729841B (en) | A kind of based on side's target model and the camera distortion bearing calibration of perspective projection | |
CN103278138B (en) | Method for measuring three-dimensional position and posture of thin component with complex structure | |
CN104517291B (en) | Pose measuring method based on target coaxial circles feature | |
CN106441242B (en) | A kind of interactive plotting method based on laser point cloud and full-view image | |
CN103175470B (en) | Reference sphere positioning and measuring method based on line-structured light vision sensor | |
CN107588721A (en) | The measuring method and system of a kind of more sizes of part based on binocular vision | |
CN112884902B (en) | Point cloud registration-oriented target ball position optimization method | |
CN103440639A (en) | Camera internal parameter calibration method based on regular triangle template | |
CN110030926B (en) | Calibration method for laser beam space pose | |
CN105427361B (en) | The display methods of moving-target track in a kind of three-dimensional scenic | |
CN104880205B (en) | Nonopiate shafting laser transit measuring system scaling method | |
CN102509324A (en) | Rotational stereovision rotary axis determining method based on quadratic curve fitting | |
CN104655106B (en) | Autonomous positioning based on GPS RTK and full-view image orients plotting method | |
CN103247048A (en) | Camera mixing calibration method based on quadratic curve and straight lines | |
CN105513063A (en) | Calibration of parabolic refraction and reflection camera through Veronese mapping and checkerboard | |
CN105374067A (en) | Three-dimensional reconstruction method based on PAL cameras and reconstruction system thereof | |
CN106546172A (en) | Three-dimensional coordinate measurement method based on nonopiate shafting laser total station | |
CN106023207A (en) | City component collecting method based on double panoramas of mobile measuring system | |
CN104807405A (en) | Three-dimensional coordinate measurement method based on light ray angle calibration | |
CN113870366A (en) | Calibration method and calibration system of three-dimensional scanning system based on pose sensor | |
Zexiao et al. | Study on a full field of view laser scanning system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120620 |