CN102385750B - Line matching method and line matching system on basis of geometrical relationship - Google Patents

Line matching method and line matching system on basis of geometrical relationship Download PDF

Info

Publication number
CN102385750B
CN102385750B CN 201110169453 CN201110169453A CN102385750B CN 102385750 B CN102385750 B CN 102385750B CN 201110169453 CN201110169453 CN 201110169453 CN 201110169453 A CN201110169453 A CN 201110169453A CN 102385750 B CN102385750 B CN 102385750B
Authority
CN
China
Prior art keywords
straight line
image
coordinate
point
same image
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.)
Active
Application number
CN 201110169453
Other languages
Chinese (zh)
Other versions
CN102385750A (en
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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN 201110169453 priority Critical patent/CN102385750B/en
Publication of CN102385750A publication Critical patent/CN102385750A/en
Application granted granted Critical
Publication of CN102385750B publication Critical patent/CN102385750B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

The invention discloses a line matching method and a line matching system on the basis of a geometrical relationship. The line matching method comprises the following steps of: a step 1 of carrying out line detection on a plurality of images; a step 2 of carrying out projection transformation on lines detected from the same image; a step 3 of obtaining an intersection point of circular arcs formed by the lines detected from the same image on a projection plane; a step 4 of generating a matching feature of the lines detected from the same image; and a step 5 of carrying out line matching by the matching features corresponding to the lines detected from each image. In the method, a fundamental matrix does not need to be obtained in advanced, the influence of the point matching precision on the precision of the line matching is avoided, and the efficiency of a matching algorithm is improved.

Description

Straight line matching process and system based on geometric relationship
Technical field
The present invention relates to image processing field, relate in particular to straight line matching process and system based on geometric relationship.
Background technology
Straight line is to constitute important geometric element in the object in the three-dimensional world.The linear feature of object still keeps in the two dimensional image that obtains from three-dimensional scenic.In the image of under the different visual angles Same Scene being taken a part of straight line remain unchanged or straight line between have certain correlativity, find the straight line of the correspondence in the different images, utilize the corresponding relation between the straight line to determine that coupling and transformation relation between two width of cloth images have important meaning.The recovery of three-dimensional scenic mainly utilizes a little and line comes out the reproductions such as object target of three-dimensional, and the straight line in the coupling different images is the important foundation of three-dimensional reconstruction.
The method of straight line coupling mainly adopts utmost point geometrical constraint to carry out coupling just, and the recycling correlation method carries out the essence coupling.Utmost point geometrical constraint method is to utilize the geometric relationship of the image that obtains under two visual angles, set up the corresponding geometric transformation relation between picture point, the geometrical constraint of point in the image at another visual angle in piece image can be represented by fundamental matrix, the fundamental matrix that obtains between image utilizes two end points of straight line can obtain the zone of meeting geometric constraint in another visual angle image, belong to the regional straight line of geometrical constraint as candidate's straight line, finish thick coupling.Smart coupling adopts relevant method, and the grey scale change of choosing pixel in the straight line neighborhood realizes coupling as similarity.
Adopt the method for utmost point geometrical constraint need ask for fundamental matrix in advance, asking for of fundamental matrix is to adopt match point that the method for asking linear equation is realized, before the straight line coupling, need to carry out a coupling, the precision receptor site matching precision influence of straight line coupling, matching algorithm efficiency is not high.
Summary of the invention
At the above-mentioned problems in the prior art, the invention provides straight line matching process and system based on geometric relationship.
The invention provides the straight line matching process based on geometric relationship, comprising:
Step 1 is carried out straight-line detection to several images;
Step 2, to from same image detection to straight line carry out projective transformation;
Step 3, ask for same image detection to the intersection point of the circular arc that forms on the projecting plane of straight line;
Step 4 generates the matching characteristic of the straight line that same image detection arrives;
Step 5, utilize from each image detection to the matching characteristic of straight line correspondence carry out the straight line coupling.
In one example, in the step 1, thereby edge of image detected the straight line that detects in the image.
In one example, in the step 2, with image
Figure GDA00003011493100021
Geometric center as the central point of coordinate transform, projection centre is the central point of unit sphere, wherein L 1, L 2, L 3..., L nBe detected straight line in the same image.
In one example, in the step 2, establishing the coordinate of coordinate transform central point in image is (q X0q Y0), in the same image on the detected straight line coordinate of any point be (q xq y), then The Transformation Relation of Projection is:
q x - q x 0 q y - q y 0 - 1 | | q x - q x 0 q y - q y 0 - 1 | | .
In one example, in the step 3, establish on first circular arc of detected two straight line correspondences in the same image and 2 coordinate vector on second circular arc is respectively
Figure GDA00003011493100023
Figure GDA00003011493100024
With
Figure GDA00003011493100026
Then the intersection point of circular arc is
± ( CP 11 1 × CP 12 1 ) × ( CP 21 1 × CP 22 1 ) | | ( CP 11 1 × CP 12 1 ) × ( CP 21 1 × CP 22 1 ) | | .
In one example, step 4 comprises:
Step 41 is chosen the corresponding distance between two points on straight line of the intersection point coordinate points p farthest of camber line Max1, p Max2, and with its mid point As new coordinate transform central point, and detected straight line in the same image carried out new projective transformation;
Figure GDA00003011493100032
Be coordinate points p Max1, p Max2At the mid point of X-axis,
Figure GDA00003011493100033
Be coordinate points p Max1, p Max2Mid point in Y-axis;
Step 42 is asked for coordinate points p Max1, p Max2The unit normal vector of place circular arc under new projection, and this unit normal vector is set to reference vector;
Step 43, the intersection point that calculates detected straight line in the same image under new projection coordinate vector and the angle of reference vector;
Step 44 sorts the back as the matching characteristic of straight line with angle.
In one example, in the step 5, utilize nearest criterion that the matching characteristic of straight line is mated.
The invention provides a kind of system that realizes based on the straight line matching process of geometric relationship.
The present invention need not to ask in advance fundamental matrix, has avoided the influence of some matching precision to the precision of straight line coupling, has improved the efficient of matching algorithm.
Description of drawings
Come the present invention is described in further detail below in conjunction with accompanying drawing, wherein:
Fig. 1 is the straight line matching process process flow diagram based on geometric relationship of the present invention;
Fig. 2 is the projective transformation illustraton of model;
Fig. 3 is that the circular arc intersection point is asked for schematic diagram.
Embodiment
The present invention proposes a kind of new straight line matching process, utilize the geometric relationship between the image cathetus to carry out matching judgment.Have parallel and overlapping relation between straight line in the image and the straight line, crossing or parallel relation between the image straight line that obtains under the different visual angles is constant under certain condition, determine the geometric relationship of straight line and other straight lines from straight line, utilize whether a straightforward differentiation of geometric relationship is same straight line.After linear projection is to the sphere, there are two intersection points in per two straight lines at sphere, this two intersection point is two limits of sphere great circle, choose limit as the expression amount of geometric relationship, with the geometric relationship employing vector representation of each bar straight line with other straight lines, utilize vector correlation to mate as matching criterior.
Method provided by the invention as shown in Figure 1, it has used 2 images (being not limited to 2 images), i.e. image 1 and image 2, but identical to each treatment of picture flow process.Method provided by the invention comprises:
Step 101, image 1 and image 2 are carried out straight-line detection.When carrying out straight-line detection, by to Image Edge-Detection, and then from the image border, detect straight line.Former two width of cloth images are designated as respectively: f 1(x y) and f 2(x y) detects image f 1Straight line n bar in (x y) is designated as respectively
Figure GDA00003011493100041
Image f 2Straight line m bar in (x y) is designated as respectively
Figure GDA00003011493100042
Belong on the straight line L i k k = 1 , i = 1,2 , . . . n k = 2 , i = 1,2 , . . . m On j the point be expressed as
Figure GDA00003011493100044
Only there is the image of rectilinear point to be designated as ∏ 1, ∏ 2, be expressed as
Figure GDA00003011493100046
Step 102, projective transformation.Be on the same orthodrome on the projecting plane through the point on the same straight line on the projective transformation image.Certainly exist two intersection points between different circular arcs, and the line of two intersection points passes through the centre of sphere.The position of the intersection point that the geometry of image cathetus and relation are different corresponding is just different, for example in the image two straight line parallels then the circular arc intersection point on the sphere is near the equator of projection model, the circular arc intersection point that two straight lines intersect in image on the sphere then departs from the equator.Position of intersecting point information response between spherical arc the geometric relationship between straight line.The projective transformation model as shown in Figure 2, with image f 1Straight line in (x y)
Figure GDA00003011493100047
The geometric relationship proper vector is asked for the process into example explanation projective transformation. Be image ∏ 1In straight line, choose ∏ 1Geometric center
Figure GDA00003011493100049
Be the central point of coordinate transform, projection centre C elects the center of sphere as, chooses unit sphere, connects point on the straight line
Figure GDA000030114931000410
C gets projection line with projection centre, the intersection point of projection line and sphere
Figure GDA000030114931000411
For Subpoint on sphere is with ∏ 1On point adopt The Transformation Relation of Projection to project on the sphere.Note: image ∏ 1In Value is at 1,
Figure GDA000030114931000414
The coordinate of coordinate transform central point in image is
Figure GDA000030114931000415
Subpoint cp then 1Coordinate is:
q x 1 - q x 0 1 q y 1 - q y 0 1 - 1 | | q x 1 - q x 0 1 q y 1 - q y 0 1 - 1 | |
Step 103, ask for how much intersection points.Same orthodrome and projection centre can be determined unique plane, and the line that intersects the intersection point of circular arc in twos was the intersection of two Plane intersects of circular arc and projection centre.Can obtain the direction vector of intersection by two normal vectors that intersect the plane.Intersection is crossed the centre of sphere as previously shown, and intersection point is positioned on the sphere, and the coordinate of intersection point can be obtained by the direction vector of intersection.With straight line
Figure GDA00003011493100052
And straight line
Figure GDA00003011493100053
Be example, image ∏ 1Cathetus
Figure GDA00003011493100054
Last 2 subpoints on sphere
Figure GDA00003011493100055
With
Figure GDA00003011493100056
Coordinate vector be
Figure GDA00003011493100057
Straight line
Figure GDA00003011493100059
Last 2 subpoints on sphere
Figure GDA000030114931000510
With
Figure GDA000030114931000511
Coordinate vector Then
Figure GDA000030114931000514
The projection circular arc with The intersection point of projection circular arc is:
Figure GDA000030114931000516
Two intersection points are about centre of sphere symmetry, and amount of orientation first component is not less than zero value and is designated as Obtain other circular arcs together
Figure GDA000030114931000518
The intersection point of projection circular arc,
Figure GDA000030114931000519
At ∏ 1On corresponding point
Figure GDA000030114931000520
In like manner ask for the intersection point between each circular arc successively, i.e. image ∏ 1Cathetus With
Figure GDA000030114931000522
The intersection point of projection circular arc is
Figure GDA000030114931000523
Corresponding point in original image are
Figure GDA000030114931000524
Step 104, generation matching characteristic.In this step, utilize the inverse transformation of projective transformation to obtain the picture point of circular arc intersection point correspondence, these points are the straight-line intersection in the image.Choose the straight-line intersection middle distance farthest 2 mid point as new transform center all intersection points are projected on the sphere again with new transform center, distance farthest the unit normal vector of the circular arc at 2 places as with reference to vector.The angle (or inner product absolute value) of the coordinate vector of intersection point and reference vector is as the expression amount of intersection point under new projection, and this expression amount is as an element of linear feature.To belong to the expression amount of intersection point on the same straight line by arranging the generating feature vector from small to large.
As example, select straight-line intersection middle distance 2 points farthest,
Figure GDA000030114931000525
Figure GDA000030114931000526
2 mid points:
p mid 1 : p mid x 1 q mid y 1 = p max 1 x 1 + p max 2 x 1 2 p max 1 y 1 + p max 2 y 1 2 . With
Figure GDA000030114931000528
For transform center with intersection point
Figure GDA00003011493100061
Project to the projection coordinate that looks for novelty on the sphere.
Straight-line intersection
Figure GDA00003011493100062
Projection coordinate's vector is under new projection:
Figure GDA00003011493100063
Account for
Figure GDA00003011493100064
Figure GDA00003011493100065
Coordinate vector at the subpoint of newly casting:
Figure GDA00003011493100066
Figure GDA00003011493100067
The unit normal vector of crossing these 2 circular arcs is: C ‾ V max 1 = C ‾ p max 1 1 × C _ p max 1 1 | | C ‾ p max 1 1 × C ‾ p max 1 1 | | , Note
Figure GDA00003011493100069
Be reference vector
Intersection point projection coordinate vector And reference vector
Figure GDA000030114931000611
Inner product:
Figure GDA000030114931000612
Straight line
Figure GDA000030114931000613
Proper vector generates:
Figure GDA000030114931000614
Go up each intersection point projection coordinate vector
Figure GDA000030114931000615
Each intersection point projection coordinate vector is with the angle of reference vector:
ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; . . . . . . ; C ‾ P 1 n - 1 1 ] T ) .
Right Sort by size the result
sort ( ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; . . . . . . ; C ‾ P 1 n - 1 1 ] T ) )
Then matching characteristic is expressed as:
MG 1 1 = sort ( ar cos ( C ‾ V max 1 · [ C ‾ P 12 1 ; C ‾ P 13 1 ; · · · · · · ; C ‾ P 1 n - 1 1 ] T ) ) .
Generate the geometric relationship matching characteristic of all straight lines,
Figure GDA000030114931000620
Adopt same step to generate image ∏ 2Feature
Figure GDA000030114931000621
Step 105, straight line coupling.Utilize the proper vector of nearest criteria match straight line.The above only is preferred implementation of the present invention, but protection domain of the present invention is not limited thereto.Any those skilled in the art all can carry out suitable change or variation to it in technical scope disclosed by the invention, and this change or variation all should be encompassed within protection scope of the present invention.

Claims (3)

1. based on the straight line matching process of geometric relationship, it is characterized in that, comprising:
Step 1 is carried out straight-line detection to several images;
Step 2, to from same image detection to straight line carry out projective transformation, wherein, with image
Figure FDA00003011493000011
Geometric center as the central point of coordinate transform, projection centre is the central point of unit sphere, wherein L 1, L 2, L 3..., L nBe detected straight line in the same image, establishing the coordinate of coordinate transform central point in image is (q X0q Y0), in the same image on the detected straight line coordinate of any point be (q xq y), then The Transformation Relation of Projection is:
Figure FDA00003011493000012
Step 3, ask for same image detection to the intersection point of the circular arc that forms on the projecting plane of straight line, establish on first circular arc of detected two straight line correspondences in the same image and 2 coordinate vector on second circular arc is respectively
Figure FDA00003011493000014
With
Figure FDA00003011493000016
Then the intersection point of circular arc is
Figure FDA00003011493000017
Step 4 generates the matching characteristic of the straight line that same image detection arrives, and is in described step 4, further comprising the steps of:
Step 41 is chosen the corresponding distance between two points on straight line of the intersection point coordinate points p farthest of camber line Max1, p Max2, and with described coordinate points p farthest Max1, p Max2Mid point
Figure FDA00003011493000018
As new coordinate transform central point, and detected straight line in the same image carried out new projective transformation; Be coordinate points p Max1, p Max2At the mid point of X-axis,
Figure FDA000030114930000110
Be coordinate points p Max1, p Max2Mid point in Y-axis;
Step 42 is asked for coordinate points p Max1, p Max2The unit normal vector of place circular arc under new projection, and this unit normal vector is set to reference vector;
Step 43, calculate detected straight-line intersection in the same image under new projection coordinate vector and the angle of reference vector;
Step 44 sorts the back as the matching characteristic of straight line with angle;
Step 5, utilize from each image detection to the matching characteristic of straight line correspondence carry out the straight line coupling.
2. the straight line matching process based on geometric relationship as claimed in claim 1 is characterized in that, in the step 1, thereby edge of image is detected the straight line that detects in the image.
3. the straight line matching process based on geometric relationship as claimed in claim 1 or 2 is characterized in that, in the step 5, utilizes nearest criterion that the matching characteristic of straight line is mated.
CN 201110169453 2011-06-22 2011-06-22 Line matching method and line matching system on basis of geometrical relationship Active CN102385750B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110169453 CN102385750B (en) 2011-06-22 2011-06-22 Line matching method and line matching system on basis of geometrical relationship

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110169453 CN102385750B (en) 2011-06-22 2011-06-22 Line matching method and line matching system on basis of geometrical relationship

Publications (2)

Publication Number Publication Date
CN102385750A CN102385750A (en) 2012-03-21
CN102385750B true CN102385750B (en) 2013-07-10

Family

ID=45825139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110169453 Active CN102385750B (en) 2011-06-22 2011-06-22 Line matching method and line matching system on basis of geometrical relationship

Country Status (1)

Country Link
CN (1) CN102385750B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345642B (en) * 2013-06-28 2016-05-25 华中科技大学 A kind of image matching method based on dotted line antithesis
CN106296645A (en) * 2015-06-25 2017-01-04 株式会社理光 Image processing method and image processing apparatus
CN105160311B (en) * 2015-08-26 2018-08-28 清华大学 Algorism of Matching Line Segments method and device
CN106934788B (en) * 2015-12-30 2020-11-24 中国科学院沈阳自动化研究所 Rapid extraction method of straight line inclination angle
CN105719309B (en) * 2016-01-27 2018-08-14 大连理工大学 A kind of matching line segments method based on projective invariant
CN107480710B (en) * 2017-08-01 2020-05-22 歌尔股份有限公司 Feature point matching result processing method and device
CN109064440B (en) * 2018-06-19 2022-02-22 广东工业大学 Loudspeaker voice coil bonding wire identification method based on machine vision

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1277079B1 (en) * 2000-04-25 2006-03-29 Rodenstock GmbH Method for calculating a progressive spectacle lens and method for producing a spectacle lens of this type
CN1941850A (en) * 2005-09-29 2007-04-04 中国科学院自动化研究所 Pedestrian tracting method based on principal axis marriage under multiple vedio cameras
CN101561269A (en) * 2009-05-26 2009-10-21 张征宇 Method for automatically matching characteristic lines of close-range photogrammetry
CN101621711A (en) * 2009-07-23 2010-01-06 东南大学 Method for calibrating camera by adopting two same circles
CN101862177A (en) * 2010-04-20 2010-10-20 中山大学中山眼科中心 Method and device for three-dimensionally positioning retinal hole

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1277079B1 (en) * 2000-04-25 2006-03-29 Rodenstock GmbH Method for calculating a progressive spectacle lens and method for producing a spectacle lens of this type
CN1941850A (en) * 2005-09-29 2007-04-04 中国科学院自动化研究所 Pedestrian tracting method based on principal axis marriage under multiple vedio cameras
CN101561269A (en) * 2009-05-26 2009-10-21 张征宇 Method for automatically matching characteristic lines of close-range photogrammetry
CN101621711A (en) * 2009-07-23 2010-01-06 东南大学 Method for calibrating camera by adopting two same circles
CN101862177A (en) * 2010-04-20 2010-10-20 中山大学中山眼科中心 Method and device for three-dimensionally positioning retinal hole

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种基于特征编组的直线立体匹配全局算法;文贡坚;《软件学报》;20061231;第17卷(第12期);2471-2484 *
文贡坚.一种基于特征编组的直线立体匹配全局算法.《软件学报》.2006,第17卷(第12期),2471-2484.

Also Published As

Publication number Publication date
CN102385750A (en) 2012-03-21

Similar Documents

Publication Publication Date Title
CN102385750B (en) Line matching method and line matching system on basis of geometrical relationship
Niu et al. Unsupervised saliency detection of rail surface defects using stereoscopic images
CN103049896B (en) The geometric data of three-dimensional model and data texturing autoregistration algorithm
Tang et al. 3D mapping and 6D pose computation for real time augmented reality on cylindrical objects
CN109671110B (en) Local geometric structure constrained urban wide baseline image feature point matching method
CN113362457B (en) Stereoscopic vision measurement method and system based on speckle structured light
Fehr et al. Covariance based point cloud descriptors for object detection and recognition
CN107886546A (en) Utilize the method for ball picture and public self-polar triangle demarcation parabolic catadioptric video camera
CN112651881B (en) Image synthesizing method, apparatus, device, storage medium, and program product
CN105243137A (en) Draft-based three-dimensional model retrieval viewpoint selection method
CN103945206B (en) A kind of stereo-picture synthesis system compared based on similar frame
CN101383051A (en) View synthesizing method based on image re-projection
CN110111255B (en) Stereo image splicing method
CN1917658B (en) Method for generating sequence of stereo images from monocular image sequence
Lin et al. 3D reconstruction by combining shape from silhouette with stereo
Yang et al. Research on 3D target pose tracking and modeling
Li et al. 3D reconstruction based on stereovision and texture mapping
JP2010237941A (en) Mask image generation device, three-dimensional model information generation device, and program
Xie On 3D reconstruction strategy: a case of conics
Chen et al. Projection ray intersecting location–based multicolour pseudo–random coded projected active vision method
Sun et al. View morphing and interpolation through triangulation
Barrena et al. Wapt-surface normal estimation for improved template matching in visual tracking
Chen et al. 3-D model reconstruction from C-arm images
Wu et al. Towards End-to-End Face Recognition Method Through Spatial Transformer
Liu et al. 3D Street Object Detection from Monocular Images Using Deep Learning and Depth Information

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