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 PDF

Info

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
Application number
CN201110338048XA
Other languages
Chinese (zh)
Inventor
陈胜勇
刘大伟
汪晓妍
刘盛
王万良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201110338048XA priority Critical patent/CN102509324A/en
Publication of CN102509324A publication Critical patent/CN102509324A/en
Pending legal-status Critical Current

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

Definite method based on the rotating stereovision turning axle of conic fitting
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:
x 2 y 2 x y 1 a b c d e = 0
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:
x 2 y 2 x y 1 a b c d e = 0
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:
x 2 y 2 x y 1 a b c d e = 0
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:
x 2 y 2 x y 1 a b c d e = 0
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.
CN201110338048XA 2011-10-31 2011-10-31 Rotational stereovision rotary axis determining method based on quadratic curve fitting Pending CN102509324A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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