CN106918302A - A kind of spatial digitizer self-adapting calibration method - Google Patents

A kind of spatial digitizer self-adapting calibration method Download PDF

Info

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
Application number
CN201710147917.8A
Other languages
Chinese (zh)
Other versions
CN106918302B (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.)
Donghua University
Original Assignee
Donghua 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 Donghua University filed Critical Donghua University
Priority to CN201710147917.8A priority Critical patent/CN106918302B/en
Publication of CN106918302A publication Critical patent/CN106918302A/en
Application granted granted Critical
Publication of CN106918302B publication Critical patent/CN106918302B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring 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

A kind of spatial digitizer self-adapting calibration method
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
v p r o j ′ = Σ i = 1 k w i v i _ p r o j ′
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,
CN201710147917.8A 2017-03-13 2017-03-13 A kind of spatial digitizer self-adapting calibration method Expired - Fee Related CN106918302B (en)

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)

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

Patent Citations (8)

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

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