CN106918302A - A kind of spatial digitizer self-adapting calibration method - Google Patents
A kind of spatial digitizer self-adapting calibration method Download PDFInfo
- Publication number
- CN106918302A CN106918302A CN201710147917.8A CN201710147917A CN106918302A CN 106918302 A CN106918302 A CN 106918302A CN 201710147917 A CN201710147917 A CN 201710147917A CN 106918302 A CN106918302 A CN 106918302A
- Authority
- CN
- China
- Prior art keywords
- point
- points
- harris angle
- proj
- collection
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Image Processing (AREA)
Abstract
The invention provides a kind of spatial digitizer self-adapting calibration method, it is characterised in that comprise the following steps:The point that Harris angle points are concentrated is divided three classes;It is every point of Harris angle points concentration, sets up its Laplce's coordinate;If 3 D scanning system there are n platform cameras, take on the basis of the Harris angle point collection that the 1st camera is gathered, the Harris angle points collection of other each camera carries out Automatic adjusument.Using after the method that the present invention is provided, the angle point overlap problem in three-dimensional scanner system marked point process can be solved, realize the self-adapting calibration to three-dimensional scanner system.
Description
Technical field
The present invention relates to a kind of method that self-adapting calibration is carried out to three-dimensional scanner system.
Background technology
When three-dimensional scanner system is built using multiple depth cameras, it is necessary to the three-dimensional point cloud that single camera is gathered
Piece is stitched together.Its splicing principle is that, using same scaling board, there is 16 × 10 chequered with black and white square dices above,
The angle point for being constituted, referred to as Harris angle points.According to the image-forming principle of camera, if the purpose of camera calibration is by cadre's phase
Several Harris angle point figures that machine is seen overlap exactly.However, because the processing technology of camera is not quite similar,
And when three-dimensional scanner system is built, (such as camera case is under different screw tension force for the objective reality of camera rigging error
Different degreeof tortuositys), angle point excessive effects can be caused.
The content of the invention
The purpose of the present invention is:Overcome angle point excessive effects, realize the self-adapting calibration to three-dimensional scanner system.
In order to achieve the above object, the technical scheme is that providing a kind of spatial digitizer self-adapting calibration side
Method, it is characterised in that comprise the following steps:
Step 1, by Harris angle points concentrate point be divided three classes, I types point be Harris angle points concentrate non-boundary point
And non-four angle points, Type-II point is the boundary point of non-four angle points that Harris angle points are concentrated, and type III point is Harris angles
Four angle points that point is concentrated;
Step 2, it is point that every Harris angle point is concentrated, sets up its Laplce's coordinate, comprises the following steps:
Step 2.1, calculating projection plane P=nxx+nyy+nzZ+d, its normal direction n are
Wherein:
viIt is the nearest-neighbors point of point v, for I types point, possesses four nearest-neighbors points, k=4;For Type-II point,
Possess three nearest-neighbors points, k=3;For type III point, possess two nearest-neighbors points, k=2;
D is average distance,
nx、ny、nzIt is components of the normal direction n on tri- directions of xyz;
Step 2.2, by point v and its nearest-neighbors spot projection to projection plane P, the locus expression of its subpoint is such as
Under:
vproj=v- (d+ (vn)) n;
vi_proj=vi-(d+(viN)) n, in formula, vprojIt is projections of the point v on projection plane P, vi_projFor its is nearest
Projection of neighbours' point on projection plane P;
Step 2.3, calculating point v points are relative to viThe tangential weight w of pointiAnd normal direction weight bi;
If step 3,3 D scanning system have n platform cameras, on the basis of taking the Harris angle point collection that the 1st camera is gathered,
The Automatic adjusument process of the Harris angle point collection of other each camera is:
Step 3.1, for jth platform camera, j ≠ 1, Harris angle point collection for, according to the method in step 2.1, meter
Calculate its projection plane being currently located, and average distance d ' and normal direction n ';
Step 3.2, for jth platform camera, j ≠ 1, Harris angle point collection in i-th point for, by its neighbours' point
v′iProject on new projection plane, gained projected position v 'i_projFor:
v′i_proj=v 'i-[d′+(v′i·n′)]n′
Step 3.3, the tangential weight w according to step 2.3i, calculate jth platform camera, j ≠ 1, Harris angle point collection in
I-th point of projected position v 'proj:
Step 3.4, calculate jth platform camera, j ≠ 1, Harris angle point collection in i-th point of final position
Preferably, tangential weight wiCalculating process be:
In formula,Wherein:
Preferably, normal direction weight biCalculating process be:
In formula,
Asked using the overlap of the angle point in three-dimensional scanner system marked point process after the method that the present invention is provided, can be solved
Topic, realizes the self-adapting calibration to three-dimensional scanner system.
Brief description of the drawings
Fig. 1 is using the effect contrast figure after the present invention.
Specific embodiment
With reference to specific embodiment, the present invention is expanded on further.It should be understood that these embodiments are merely to illustrate the present invention
Rather than limitation the scope of the present invention.In addition, it is to be understood that after the content for having read instruction of the present invention, people in the art
Member can make various changes or modifications to the present invention, and these equivalent form of values equally fall within the application appended claims and limited
Scope.
The invention provides a kind of spatial digitizer self-adapting calibration method, comprise the following steps:
Step 1, by Harris angle points concentrate point be divided three classes, I types point be Harris angle points concentrate non-boundary point
And non-four angle points, Type-II point is the boundary point of non-four angle points that Harris angle points are concentrated, and type III point is Harris angles
Four angle points that point is concentrated;
Step 2, it is point that every Harris angle point is concentrated, sets up its Laplce's coordinate, comprises the following steps:
Step 2.1, calculating projection plane P=nxx+nyy+nzZ+d, its normal direction n are Wherein:
viIt is the nearest-neighbors point of point v, for I types point, possesses four nearest-neighbors points, k=4;For Type-II point,
Possess three nearest-neighbors points, k=3;For type III point, possess two nearest-neighbors points, k=2;
D is average distance,
nx、ny、nzIt is components of the normal direction n on tri- directions of xyz;
Step 2.2, by point v and its nearest-neighbors spot projection to projection plane P, the locus expression of its subpoint is such as
Under:
vproj=v- (d+ (vn)) n;
vi_proj=vi-(d+(viN)) n, in formula, vprojIt is projections of the point v on projection plane P, vi_projFor its is nearest
Projection of neighbours' point on projection plane P;
Step 2.3, calculating point v points are relative to viThe tangential weight w of pointiAnd normal direction weight bi;
Tangential weight wiCalculating process be:
In formula,Wherein:
Normal direction weight biCalculating process be:
In formula,
If step 3,3 D scanning system have n platform cameras, on the basis of taking the Harris angle point collection that the 1st camera is gathered,
The Automatic adjusument process of the Harris angle point collection of other each camera is:
Step 3.1, for jth platform camera, j ≠ 1, Harris angle point collection for, according to the method in step 2.1, meter
Calculate its projection plane being currently located, and average distance d ' and normal direction n ';
Step 3.2, for jth platform camera, j ≠ 1, Harris angle point collection in i-th point for, by its neighbours' point
v′iProject on new projection plane, gained projected position v 'i_projFor:
v′i_proj=v 'i-[d′+(v′i·n′)]n′
Step 3.3, the tangential weight w according to step 2.3i, calculate jth platform camera, j ≠ 1, Harris angle point collection in
I-th point of projected position v 'proj:
Step 3.4, calculate jth platform camera, j ≠ 1, Harris angle point collection in i-th point of final position
Fig. 1 is shown and is compared using the effect of camera collection three-dimensional point cloud before and after the technology of the present invention, it can be seen that by this
After inventive technique is demarcated to camera, high-quality three-dimensional point cloud can be generated.
Claims (3)
1. a kind of spatial digitizer self-adapting calibration method, it is characterised in that comprise the following steps:
Step 1, the point that Harris angle points are concentrated is divided three classes, I types point is the non-boundary point and non-that Harris angle points are concentrated
Four angle points, Type-II point is the boundary point of non-four angle points that Harris angle points are concentrated, and type III point is Harris angle point collection
In four angle points;
Step 2, it is point that every Harris angle point is concentrated, sets up its Laplce's coordinate, comprises the following steps:
Step 2.1, calculating projection plane P=nxx+nyy+nzZ+d, its normal direction n are
Wherein:
viIt is the nearest-neighbors point of point v, for I types point, possesses four nearest-neighbors points, k=4;For Type-II point, possess
Three nearest-neighbors points, k=3;For type III point, possess two nearest-neighbors points, k=2;
D is average distance,
nx、ny、nzIt is components of the normal direction n on tri- directions of xyz;
Step 2.2, by point v and its nearest-neighbors spot projection to projection plane P, the locus of its subpoint is expressed as follows:
vproj=v- (d+ (vn)) n;
vi_proj=vi-(d+(viN)) n, in formula, vprojIt is projections of the point v on projection plane P, vi_projIt is its nearest-neighbors
Projection of the point on projection plane P;
Step 2.3, calculating point v points are relative to viThe tangential weight w of pointiAnd normal direction weight bi;
If step 3,3 D scanning system have n platform cameras, on the basis of taking the Harris angle point collection that the 1st camera is gathered, other
The Automatic adjusument process of each Harris angle point collection of camera is:
Step 3.1, for jth platform camera, j ≠ 1, Harris angle point collection for, according to the method in step 2.1, calculate it
The projection plane being currently located, and average distance d ' and normal direction n ';
Step 3.2, for jth platform camera, j ≠ 1, Harris angle point collection in i-th point for, by its neighbours' point v 'iProjection
Onto new projection plane, gained projected position v 'i_projFor:
v′i_proj=v 'i-[d′+(v′i·n′)]n′
Step 3.3, the tangential weight w according to step 2.3i, calculate jth platform camera, j ≠ 1, Harris angle point collection in i-th
The projected position v ' of pointproj:
Step 3.4, calculate jth platform camera, j ≠ 1, Harris angle point collection in i-th point of final position v ',
2. a kind of spatial digitizer self-adapting calibration method as claimed in claim 1, it is characterised in that in the step 1.3
In, tangential weight wiCalculating process be:
In formula,Wherein:
3. a kind of spatial digitizer self-adapting calibration method as claimed in claim 1, it is characterised in that in the step 1.3
In, normal direction weight biCalculating process be:
In formula,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710147917.8A CN106918302B (en) | 2017-03-13 | 2017-03-13 | A kind of spatial digitizer self-adapting calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710147917.8A CN106918302B (en) | 2017-03-13 | 2017-03-13 | A kind of spatial digitizer self-adapting calibration method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106918302A true CN106918302A (en) | 2017-07-04 |
CN106918302B CN106918302B (en) | 2019-09-27 |
Family
ID=59461115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710147917.8A Expired - Fee Related CN106918302B (en) | 2017-03-13 | 2017-03-13 | A kind of spatial digitizer self-adapting calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106918302B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020101438A1 (en) * | 2001-01-31 | 2002-08-01 | Harris Corporation | System and method for identifying tie point collections used in imagery |
CN101321303A (en) * | 2008-07-17 | 2008-12-10 | 上海交通大学 | Geometric and optical correction method for non-plane multi-projection display |
CN101763643A (en) * | 2010-01-07 | 2010-06-30 | 浙江大学 | Automatic calibration method for structured light three-dimensional scanner system |
JP2011197120A (en) * | 2010-03-17 | 2011-10-06 | Toppan Printing Co Ltd | Pattern evaluation method and pattern evaluation device |
CN103236064A (en) * | 2013-05-06 | 2013-08-07 | 东南大学 | Point cloud automatic registration method based on normal vector |
CN104517316A (en) * | 2014-12-31 | 2015-04-15 | 中科创达软件股份有限公司 | Three-dimensional modeling method and terminal equipment |
CN104990515A (en) * | 2015-06-02 | 2015-10-21 | 江苏科技大学 | Three-dimensional shape measurement system and method for large-size object |
CN105551039A (en) * | 2015-12-14 | 2016-05-04 | 深圳先进技术研究院 | Calibration method and calibration device for structured light 3D scanning system |
-
2017
- 2017-03-13 CN CN201710147917.8A patent/CN106918302B/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020101438A1 (en) * | 2001-01-31 | 2002-08-01 | Harris Corporation | System and method for identifying tie point collections used in imagery |
CN101321303A (en) * | 2008-07-17 | 2008-12-10 | 上海交通大学 | Geometric and optical correction method for non-plane multi-projection display |
CN101763643A (en) * | 2010-01-07 | 2010-06-30 | 浙江大学 | Automatic calibration method for structured light three-dimensional scanner system |
JP2011197120A (en) * | 2010-03-17 | 2011-10-06 | Toppan Printing Co Ltd | Pattern evaluation method and pattern evaluation device |
CN103236064A (en) * | 2013-05-06 | 2013-08-07 | 东南大学 | Point cloud automatic registration method based on normal vector |
CN104517316A (en) * | 2014-12-31 | 2015-04-15 | 中科创达软件股份有限公司 | Three-dimensional modeling method and terminal equipment |
CN104990515A (en) * | 2015-06-02 | 2015-10-21 | 江苏科技大学 | Three-dimensional shape measurement system and method for large-size object |
CN105551039A (en) * | 2015-12-14 | 2016-05-04 | 深圳先进技术研究院 | Calibration method and calibration device for structured light 3D scanning system |
Non-Patent Citations (2)
Title |
---|
何欣荣等: "基于差分形态分解的多尺度Harris角点检测器", 《应用科技》 * |
李端等: "基于骨架重合的真实人体模型动态仿真", 《纺织学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106918302B (en) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10972630B2 (en) | Method for flattening laser-based image of curved book page | |
WO2020135446A1 (en) | Target positioning method and device and unmanned aerial vehicle | |
KR100966592B1 (en) | Method for calibrating a camera with homography of imaged parallelogram | |
EP3435282A3 (en) | Laser speckle system and method for an aircraft | |
JP2021016153A5 (en) | ||
US9430865B2 (en) | Real-time dynamic non-planar projection apparatus and method | |
CN104732539A (en) | Projector calibration method | |
JP2018514783A5 (en) | Method for calculating the distance to an object | |
CN106331527A (en) | Image splicing method and device | |
CN105096283A (en) | Panoramic image acquisition method and device | |
CN102402785B (en) | Camera self-calibration method based on quadratic curves | |
JP2011128117A5 (en) | Information processing apparatus, pattern data generation apparatus, information processing method, pattern data generation method, and program | |
WO2020147574A1 (en) | Deep-learning-based stereo matching method for binocular dynamic vision sensor | |
CN103852060A (en) | Visible light image distance measuring method based on monocular vision | |
CN102103746A (en) | Method for calibrating parameters in camera through solving circular ring points by utilizing regular tetrahedron | |
CN106709865A (en) | Depth image synthetic method and device | |
TW201243764A (en) | Orthoimage color correction method using multiple aerial images | |
CN111784648A (en) | Soft material fitting precision detection method, device, equipment and storage medium | |
CN104123725B (en) | A kind of computational methods of single line array camera homography matrix H | |
CN105139336B (en) | A kind of method of multichannel full-view image conversion ball curtain flake film | |
CN103994779A (en) | Panorama camera calibrating method based on three-dimensional laser-point cloud | |
CN110363801A (en) | The corresponding point matching method of workpiece material object and workpiece three-dimensional CAD model | |
CN106918302A (en) | A kind of spatial digitizer self-adapting calibration method | |
US10281265B2 (en) | Method and system for scene scanning | |
CN207133863U (en) | For eliminating the device of the overlapping deviation of 360 degree of shooting images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190927 |