CN100357701C - Grid type sub-pixel extraction of calibrated points - Google Patents
Grid type sub-pixel extraction of calibrated points Download PDFInfo
- Publication number
- CN100357701C CN100357701C CNB2005100827640A CN200510082764A CN100357701C CN 100357701 C CN100357701 C CN 100357701C CN B2005100827640 A CNB2005100827640 A CN B2005100827640A CN 200510082764 A CN200510082764 A CN 200510082764A CN 100357701 C CN100357701 C CN 100357701C
- Authority
- CN
- China
- Prior art keywords
- point
- grid
- pixel
- sub
- location
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
The invention relates to the modification to the extracting method for the grid-point sub image element of the grid-type index point. The process is: first to determine the image element position of the grid point by the Hessian matrix; then to determine the sub image element position of the grid point by the decked search. So it can get the sub image element position quickly and the process is simple. Even the picture has the aberrance and big yawp; it also has the high extracting precision.
Description
Technical field
The invention belongs to Computerized three-dimensional vision measurement technology, relate to improvement grid target point grid point sub-pixel extraction.
Background technology
In the Computerized three-dimensional vision measurement, usually need utilize the three dimensional space coordinate of target calibration point and the model parameter that the two dimensional image coordinate comes calibration system thereof, wherein, the target calibration point pattern that is adopted and the extraction algorithm of respective image coordinate will directly influence the precision and the stability of calibration result.The grid target point (see figure 1) is that the target pattern is used in a kind of vision system demarcation relatively more commonly used, in order to guarantee the stated accuracy of system, need reach sub-pixel to the extraction result of calibration point image coordinate usually.Traditional grid point sub-pixel extraction mainly is to extract near the lines central point of grid point earlier, carries out fitting a straight line then, determines the sub-pixel location of grid point again by the intersection point that calculates two straight lines.The shortcoming of this algorithm mainly contains: need to extract earlier the lines central point on the one hand, leaching process is loaded down with trivial details; When there was distortion in camera lens, near the lines central point the grid point was not to distribute in strict accordance with straight line on the other hand, therefore, utilized the method for fitting a straight line find intersection bring error can for actual extraction, and it is not high to extract precision.
Summary of the invention
The objective of the invention is:, propose a kind of grid target point sub-pixel extraction that extract, that leaching process is simple, precision is high that is used at the deficiency of existing method.
Technical scheme of the present invention is: a kind of grid type sub-pixel extraction of calibrated points, it is characterized in that,
1, utilize the Hessian matrix to judge the location of pixels of grid point;
Decision condition obtains according to the eigenwert of Hessian matrix, the expression formula of Hessian matrix as shown in Equation 1:
Wherein, f
Xx, f
Xy, f
YyBe respectively the second-order partial differential coefficient of gradation of image with respect to x, y, Gauss operator convolution by gradation of image and corresponding differential form obtains, for grid target point, two eigenwerts of Hessian matrix just are, and judge that then the shape operator S of grid type point location of pixels is expressed as:
Wherein, λ
1, λ
2Be the eigenwert of Hessian matrix, ideally the positive extreme point of the S value of grid point near zone distribution is the location of pixels of grid point;
2, determine the sub-pixel location of grid point by hierarchical search;
2.1, near the grid point in the ideal image ± sub-pixel location point in 1 pixel coverage carries out gray-level interpolation, make the curved surface that obtained behind the gray-level interpolation be zero at the grid point place along the first order derivative of all directions, the constraint condition that obtains searching for the grid point sub-pixel location is as follows:
min(r
x 2(x
0+s,y
0+t)+r
y 2(x
0+s,y
0+t))|
(s,t)∈[-1,1]×[-1,1](3)
Wherein, (x
0, y
0) be the location of pixels of grid point, r
x(x
0+ s, y
0+ t), r
y(x
0+ s, y
0+ t) be respectively point (x
0+ s, y
0+ t) locate the single order partial derivative of gradation of image with respect to x, y;
2.2, be the center with the location of pixels that obtains in 1.1, to its ± sub-pixel location point in 1 pixel coverage carries out gray-level interpolation, interpolation is elected 0.1 pixel at interval as, further judges the grid central point according to constraint condition (3) then; Be the center with this point again, to its ± sub-pix point in 0.1 pixel coverage carries out gray-level interpolation, interpolation is spaced apart 0.01 pixel, and then judges grid point according to constraint condition (3); By that analogy, carry out hierarchical search, the search number of plies is 2~6 layers, till satisfying the predetermined precision requirement.
Advantage of the present invention is:
(1) with respect to the extracting method of traditional fitting a straight line find intersection, the present invention is according to the gray feature of topography, can obtain the sub-pixel location of grid point very soon by hierarchical search, and leaching process is simple.
(2) the present invention just utilizes the gradation of image in the grid point close region, therefore, even exist at image under the situation of distortion and big noise, still has higher extraction precision.
Description of drawings
Fig. 1 is the target synoptic diagram with grid target point.
Fig. 2 is that the S value of ideally grid point near zone distributes, and its positive extreme point is a grid point.
Fig. 3 be near the grid point in the ideal image ± sub-pixel location point in 1 pixel coverage carries out the result after the gray-level interpolation.
Fig. 4 is the hierarchical search synoptic diagram in the inventive method.
Fig. 5 is the virtual plane target that the embodiment of the invention is used.
Fig. 6 is the target image of Fig. 5 target after the virtual video camera imaging.
Embodiment
The present invention mainly is the variation of image grayscale according to the grid point near zone, utilizes the location of pixels of the eigenwert judgement grid point of Hassian matrix, on this basis, determines the sub-pixel location of grid point by hierarchical search.Concrete scheme is as follows:
(1) utilize the Hessian matrix to judge the location of pixels of grid point
Before the sub-pixel location of extracting grid point, need to judge earlier the location of pixels of this point, decision condition can obtain according to the eigenwert of Hessian matrix.The expression formula of Hessian matrix is as shown in Equation 1:
Wherein, f
Xx, f
Xy, f
YyBe respectively the second-order partial differential coefficient of gradation of image, can obtain by the Gauss operator convolution of gradation of image and corresponding differential form with respect to x, y.For grid target point, two eigenwerts of Hessian matrix just are, and judge that then the shape operator of grid type point location of pixels is expressed as:
Wherein, λ
1, λ
2Eigenwert for the Hessian matrix.Fig. 2 is that the S value of ideally grid point near zone distributes, and its positive extreme point is a grid point.
(2) determine the sub-pixel location of grid point by hierarchical search
After location of pixels, can determine that the sub-pixel location of grid point promptly is positioned at this pixel adjacent domain according to shape operator preliminary judgement grid point.Fig. 3 be near the grid point in the ideal image ± sub-pixel location point in 1 pixel coverage carries out the result after the gray-level interpolation, this curved surface is zero at the grid point place along the first order derivative of all directions, can search for the constraint condition of grid point sub-pixel location:
min(r
x 2(x
0+s,y
0+t)+r
y 2(x
0+s,y
0+t))|
(s,t)∈[-1,1]×[-1,1](3)
Wherein, (x
0, y
0) be the location of pixels of grid point, r
x(x
0+ s, y
0+ t), r
y(x
0+ s, y
0+ t) be respectively point (x
0+ s, y
0+ t) locate the single order partial derivative of gradation of image with respect to x, y.
In order to improve search speed, can adopt the method for hierarchical search.As shown in Figure 4, utilize shape operator (2) to judge the location of pixels of grid point earlier, be the center with this point then, to its ± sub-pixel location point in 1 pixel coverage carries out gray-level interpolation, interpolation is elected 0.1 pixel at interval as, further judge the grid central point according to constraint condition (3) then, be the center with this point again, to its ± sub-pix point in 0.1 pixel coverage carries out gray-level interpolation, interpolation is spaced apart 0.01 pixel, and then according to constraint condition (3) judgement grid point, by that analogy.When carrying out actual grid point extraction, the search number of plies is made as 3, and last interpolation is made as 0.001 pixel at interval, can satisfy actual calibration request.
Simulation example
It below is one group of simulation example.Fig. 5 is a virtual plane target that is generated by computing machine, has 144 grid points on the target, every distance between two points 16mm.The image resolution ratio of virtual video camera is made as 512 * 512; Inner parameter is set at: α=1000, β=1000, γ=0, u
0=256, v
0=256; External parameter is set at: r
1=[0.951-0.174 0.255]
T, r
2=[0.168 0.9850.045]
T, t=[00500]
TThe lens distortion parameter setting is: k
1=-0.22, k
2=0.20.Fig. 6 is the target image after the virtual video camera imaging.Table 1 is the contrast that classic method and new method are extracted the result, and extraction precision of the present invention is better than classic method.
The extraction precision of two kinds of methods of table 1
Noise level | Extract precision | |
Classic method | |
|
0 | 0.105 | 0.048 |
0.02 | 0.106 | 0.054 |
0.04 | 0.108 | 0.069 |
0.06 | 0.108 | 0.089 |
0.08 | 0.112 | 0.111 |
Claims (2)
1, a kind of grid type sub-pixel extraction of calibrated points is characterized in that,
1.1, utilize the Hessian matrix to judge the location of pixels of grid point;
Decision condition obtains according to the eigenwert of Hessian matrix, the expression formula of Hessian matrix as the formula (1):
Wherein, f
Xx, f
Xy, f
YyBe respectively the second-order partial differential coefficient of gradation of image with respect to x, y, Gauss operator convolution by gradation of image and corresponding differential form obtains, for grid target point, two eigenwerts of Hessian matrix just are, and judge that then the shape operator S of grid type point location of pixels is expressed as:
Wherein, λ
1, λ
2Be the eigenwert of Hessian matrix, ideally the positive extreme point of the S value of grid point near zone distribution is the location of pixels of grid point;
1.2, determine the sub-pixel location of grid point by hierarchical search;
1.2.1, near the grid point in the ideal image ± sub-pixel location point in 1 pixel coverage carries out gray-level interpolation, make the curved surface that obtained behind the gray-level interpolation be zero at the grid point place along the first order derivative of all directions, the constraint condition that obtains searching for the grid point sub-pixel location is as follows:
Wherein, (x
0, y
0) be the location of pixels of grid point, r
x(x
0+ s, y
0+ t), r
y(x
0+ s, y
0+ t) be respectively point (x
0+ s, y
0+ t) locate the single order partial derivative of gradation of image with respect to x, y;
1.2.2, be the center with the location of pixels that obtains in 1.1, to its ± sub-pixel location point in 1 pixel coverage carries out gray-level interpolation, interpolation is elected 0.1 pixel at interval as, further judges the grid central point according to constraint condition (3) then; Be the center with this point again, to its ± sub-pix point in 0.1 pixel coverage carries out gray-level interpolation, interpolation is spaced apart 0.01 pixel, and then judges grid point according to constraint condition (3); By that analogy, carry out hierarchical search, the search number of plies is 2~6 layers, till satisfying the predetermined precision requirement.
2, grid type sub-pixel extraction of calibrated points according to claim 1 is characterized in that, when carrying out hierarchical search, the search number of plies is 3 layers, the 3rd layer of when search, near the central point ± sub-pix point in 0.01 pixel coverage carries out gray-level interpolation, interpolation is spaced apart 0.001 pixel.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100827640A CN100357701C (en) | 2005-07-12 | 2005-07-12 | Grid type sub-pixel extraction of calibrated points |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100827640A CN100357701C (en) | 2005-07-12 | 2005-07-12 | Grid type sub-pixel extraction of calibrated points |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1896678A CN1896678A (en) | 2007-01-17 |
CN100357701C true CN100357701C (en) | 2007-12-26 |
Family
ID=37609256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100827640A Expired - Fee Related CN100357701C (en) | 2005-07-12 | 2005-07-12 | Grid type sub-pixel extraction of calibrated points |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100357701C (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100491903C (en) * | 2007-09-05 | 2009-05-27 | 北京航空航天大学 | Method for calibrating structural parameter of structure optical vision sensor |
CN101867815B (en) * | 2010-04-30 | 2011-09-14 | 西北工业大学 | Rapid fractional pixel hierarchical searching method |
JP5833958B2 (en) * | 2012-03-14 | 2015-12-16 | 富士フイルム株式会社 | Image processing apparatus and method, and program |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6192156B1 (en) * | 1998-04-03 | 2001-02-20 | Synapix, Inc. | Feature tracking using a dense feature array |
US6208769B1 (en) * | 1998-05-28 | 2001-03-27 | Acuity Imaging, Llc | Method of accurately locating the fractional position of a template match point |
-
2005
- 2005-07-12 CN CNB2005100827640A patent/CN100357701C/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6192156B1 (en) * | 1998-04-03 | 2001-02-20 | Synapix, Inc. | Feature tracking using a dense feature array |
US6208769B1 (en) * | 1998-05-28 | 2001-03-27 | Acuity Imaging, Llc | Method of accurately locating the fractional position of a template match point |
Non-Patent Citations (2)
Title |
---|
一种应用于视觉测量标定的亚像素级特征点提取算法 袁萍,袁证鹏.现代机械,第4期 2003 * |
圆结构光光条中心亚像素级提取方法 杨宪铭等.光电工程,第31卷第4期 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN1896678A (en) | 2007-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104359403B (en) | Planar part dimension measurement method based on sub-pixel edge algorithm | |
CN106340044B (en) | Join automatic calibration method and caliberating device outside video camera | |
CN109000557B (en) | A kind of nuclear fuel rod pose automatic identifying method | |
CN107851314A (en) | For the method for the optimization for analyzing surface of tyre uniformity | |
CN105160702A (en) | Stereoscopic image dense matching method and system based on LiDAR point cloud assistance | |
CN104574415B (en) | Target space positioning method based on single camera | |
CN105783779B (en) | The real-time form identification of rail profile and distortion calibration method based on three layers of matching | |
CN101178812A (en) | Mixed image processing process of structure light striation central line extraction | |
CN100468457C (en) | Method for matching depth image | |
CN103955927A (en) | Fillet weld automatic tracking method based on laser vision | |
CN106530334B (en) | A kind of airborne Interference synthetic aperture radar complex image corregistration method and complex image corregistration system | |
CN103168219A (en) | Method for identifying and defining basic patterns forming the tread design of a tyre | |
CN101900536A (en) | Method for measuring object surface appearance based on digital picture method | |
CN107702663A (en) | A kind of point cloud registration method based on the rotation platform with index point | |
CN105938619A (en) | Visual odometer realization method based on fusion of RGB and depth information | |
CN106600647A (en) | Binocular visual multi-line projection structured light calibration method | |
CN104751458B (en) | A kind of demarcation angular-point detection method based on 180 ° of rotation operators | |
CN107240073A (en) | A kind of 3 d video images restorative procedure merged based on gradient with clustering | |
CN103886642B (en) | A kind of surface of steel plate three-dimensional reconstruction Fast implementation | |
CN104089599B (en) | Quasi morphological filtering method for extracting two-dimensional contour in contact measuring head measurement | |
CN101178811A (en) | Image characteristic extracting method of structure light straight-line striation of three steps method | |
CN100357701C (en) | Grid type sub-pixel extraction of calibrated points | |
CN111008602B (en) | Scribing feature extraction method combining two-dimensional vision and three-dimensional vision for small-curvature thin-wall part | |
CN101719287B (en) | Method for rebuilding shape of hemisphere three-dimensional surface with control point information | |
CN109325298A (en) | Tyre tread designing system based on three dimensional design platform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071226 Termination date: 20110712 |