CN101957175B  Threepoint microplanebased normal detection method  Google Patents
Threepoint microplanebased normal detection method Download PDFInfo
 Publication number
 CN101957175B CN101957175B CN 201010281290 CN201010281290A CN101957175B CN 101957175 B CN101957175 B CN 101957175B CN 201010281290 CN201010281290 CN 201010281290 CN 201010281290 A CN201010281290 A CN 201010281290A CN 101957175 B CN101957175 B CN 101957175B
 Authority
 CN
 China
 Prior art keywords
 coordinate
 measuring machine
 laser
 normal
 displacement sensor
 Prior art date
Links
 238000001514 detection method Methods 0.000 title claims abstract 4
 238000006073 displacement reactions Methods 0.000 claims abstract description 20
 238000000034 methods Methods 0.000 claims abstract description 11
 238000005457 optimization Methods 0.000 claims abstract description 5
 238000006243 chemical reactions Methods 0.000 claims description 14
 230000001131 transforming Effects 0.000 claims description 6
 210000001503 Joints Anatomy 0.000 claims description 3
 239000011159 matrix materials Substances 0.000 claims description 3
 238000004364 calculation methods Methods 0.000 claims description 2
 238000006467 substitution reactions Methods 0.000 claims description 2
 238000005516 engineering processes Methods 0.000 description 10
 235000013350 formula milk Nutrition 0.000 description 8
 238000004519 manufacturing process Methods 0.000 description 4
 238000010586 diagrams Methods 0.000 description 3
 241001074085 Scophthalmus aquosus Species 0.000 description 1
 230000000694 effects Effects 0.000 description 1
Abstract
Description
Technical field
The present invention relates to a kind of normal direction detection technique in the field of measurement, specifically is a kind of based on design and the correlation parameter scaling method of Surface Method 3 little planar processes, the employing noncontact measurement to measuring method.
Background technology
In commercial production, measuring technique is being brought into play very important effect, along with development of technology; The particularly development of laser technology, sensor technology; Measurement means is variation and evolution earlier more, and the measurand that requires in the production run simultaneously is also more and more, many times need measure the normal vector on a certain curved surface; Like the brill riveting in aircraft manufacturing process technology, the method vector coincidence at point to be processed place on the processing axis that just requires boring and riveting machine and the aircraft skin; Domestic technology and document about this respect is fewer at present, therefore proposes Surface Method to measurement, and designs a kind of measuring method and correlation parameter scaling method, and the field of measurement in commercial production has the significant values meaning.
Summary of the invention
The object of the present invention is to provide the scaling method of a kind of Surface Method to measurement and correlation parameter.The present invention has utilized laser sensor technology and data acquisition technology and based on the design philosophy of 3 little planar processes, can carry out the measurement of curved surface normal direction accurately, quickly and easily.
The present invention realizes that through following technical scheme implementation step is following:
1) set up model: 3 laser displacement sensor annulars are distributed in the head of three coordinate measuring machine, and the spatial attitude of sensor can make that laser is tapered to be penetrated on curved surface, controls 3 formed little plane sizes of laser point according to accuracy requirement adjustment tapering;
2) set up space coordinates: at first, set up a world coordinate system 1 with three coordinate measuring machine; Select for use the terminal joint coordinate system of three coordinate measuring machine as intermediate conversion coordinate system 2 then; On each laser displacement sensor, setting up one in addition penetrates direction with laser rays and is z axle negative sense, is the local coordinate system 3 of true origin with laser displacement sensor work RP;
3) image data: utilize data collecting card and LabVIEW software collection laser displacement sensor voltage signal and convert shift value L into _{i}
4) calculate the coordinate of laser point: the displacement that utilizes laser displacement sensor to record promptly know the coordinate of laser point under local coordinate system 3 (0,0, w _{i}), w _{i}=L _{i}, according to the Eulerian coordinates conversion formula, try to achieve the coordinate of 3 laser point under coordinate system 1 respectively, the coordinate conversion formula is following:
P in the formula _{i}The coordinate of expression laser point; R _{0}Be the terminal pose of three coordinate measuring machine, from three coordinate measuring machine software, obtain; R _{i}Be the homogeneous transformation matrix of local coordinate system 3 with respect to intermediate conversion coordinate system 2, wherein r _{Z1}, r _{X2}Be respectively coordinate system 3 around the z of coordinate system 2 axle, the x axle anglec of rotation, (a _{i}, b _{i}, c _{i}) be the translation vector of 3 pairs of coordinate systems 2 of coordinate system, to r _{Z1}, r _{X2}, (a _{i}, b _{i}, c _{i}) demarcate;
5) computing method vector: 3 laser point coordinates of trying to achieve with step 4) can solve 3 vectors, and wherein any 2 vectorial multiplication crosses get final product to such an extent that 3 little planar processes are vectorial again, and formula is following:
The normal vector that calculates is applied to normal error correction, the adjustment of normal direction attitude, normal direction detection.
Described laser displacement sensor attitude parameter scaling method: implementation step is following:
1) according to step 1) in (1) laser displacement sensor is installed on the three dimensional coordinate measuring machine head and sets up following measurement space: according to step 2 in (1)) use three coordinate measuring machine to set up a world coordinate system 1 and set up an intermediate conversion coordinate system 2, in world coordinate system 1, draft a spherical equation (xa) then with standard ball ^{2}+ (yb) ^{2}+ (zc) ^{2}=r ^{2}, wherein, (a, b, c) expression sphere centre coordinate, r representes the radius of a ball, and sphere centre coordinate is demarcated with three coordinate measuring machine, and the radius of a ball is known;
2) on sphere, evenly gather 30 measuring points with one of them sensor operated three coordinate measuring machine, obtain 30 shift values;
3) the shift value substitution spherical equation that records is used fminsearch minimum optimization method and find the solution each parameter value;
4) repeating step 2), step 3) finds the solution the attitude parameter of other 2 sensors.
Design concept of the present invention replaces little curved surface based on differential geometric thoughtlittle plane approximation, has adopted laser sensor technology, not only can carry out noncontact measurement, and measuring accuracy height, convenient measurement, fast; Calibration algorithm among the present invention is simple, easy in addition realizes, and when application fminsearch minimum optimization is found the solution, initial value is chosen less demanding.
Description of drawings
Fig. 1 is a principle key diagram of the present invention;
Fig. 2 is a model synoptic diagram of the present invention;
(a) figure is a world coordinate system 1;
(b) figure is that laser displacement sensor is arranged synoptic diagram and intermediate conversion coordinate system 2;
(c) be laser displacement sensor local coordinate system 3;
Below in conjunction with accompanying drawing content of the present invention is done further explain.
Embodiment
With reference to shown in Figure 1, be based on infinitesimal geometry thought, choosing 1 the 3 small plane of point type near certain point on the curved surface, with the approximate normal vector that replaces this some place on the curved surface of the normal vector on this small plane, resultant error receives the influence of little plane sizes.Therefore, little plane sizes is chosen according to accuracy requirement, 1 expression sensor.
With reference to shown in Figure 2,1 expression sensor, 2 expression spheres.Implementation step of the present invention is: set up measurement modelattitude parameter demarcation, and Surface Method is to measurement, and embodiment is following:
(1) 3 laser sensors is installed on the three coordinate measuring machine, and makes 3 laser of ejaculation constitute taper, regulate the size that tapering is controlled little plane according to accuracy requirement.
(2) use three coordinate measuring machine to set up a world coordinate system 1; It is the coordinate system 1 in the accompanying drawing 2; Selecting the coordinate system of the terminal joint of three coordinate measuring machine for use is intermediate conversion coordinate system 2; Being coordinate system 2 in the accompanying drawing 2, is the Z axle along laser rays respectively in addition, is that initial point is set up 3 local coordinate systems 3 with the working sensor RP.The relation of each coordinate system is following: coordinate system 1 is the terminal pose of three coordinate measuring machine with the transformation relation of coordinate system 2, from three coordinate measuring machine software, obtains; Coordinate system 3 is that the laser displacement sensor attitude parameter is to be calibrated with the transformation relation of coordinate system 2.
(3) in world coordinate system 1, place a standard ball, use three coordinate measuring machine to demarcate sphere centre coordinate, and set up spherical equation: (xa) according to the standard ball radius ^{2}+ (yb) ^{2}+ (zc) ^{2}=r ^{2}, be used for parameter calibration.
(4) with the singlesensor be research object, the control three coordinate measuring machine is gathered 30 measuring points on the standard sphere, collect the shift value L of each measuring point through data acquisition equipment _{i}, can obtain the coordinate of measuring point under local coordinate system 3 (0,0, w _{i}), w _{i}=L _{i}
(5) carry out coordinate transform: with the measuring point coordinate (0,0, w _{i}) transform in the world coordinate system 1, transformation for mula does
P in the formula _{i}The coordinate of expression laser point; R _{0}Be the terminal pose of three coordinate measuring machine, from three coordinate measuring machine software, obtain; R _{i}Be the homogeneous transformation matrix of local coordinate system 3 with respect to intermediate conversion coordinate system 2, wherein r _{Z1}, r _{X2}Be respectively coordinate system 3 around the z of coordinate system 2 axle, the x axle anglec of rotation, (a _{i}, b _{i}, c _{i}) be the translation vector of 3 pairs of coordinate systems 2 of coordinate system, each parameter is to be calibrated.Coordinate satisfies spherical equation (x after the conversion _{P}a) ^{2}+ (y _{P}b) ^{2}+ (z _{P}c) ^{2}=r ^{2}Obtain 30 spherical equations.
(6) in matlab, use fminsearch minimum optimization method and find the solution 30 spherical equations that obtain, can solve each parameter in the coordinate transform formula.
(7) after each parametric solution comes out, then can carry out the curved surface normal direction and measure.Implementation algorithm is: through 3 sensor acquisition data; On curved surface, obtain the volume coordinate of 3 points; Can constitute 3 vectors, any 2 vectorial multiplication crosses can be tried to achieve little planar process vector, i.e. Surface Method vector; Algorithm realizes that with programming software the vector that draws is used for field of measurement such as normal error correction, the adjustment of normal direction attitude.
Bore in the riveting technology in the aircraft robotization, require cutter usually, to guarantee the crudy in hole along the normal direction feeding of covering processing stand place.Need carry out automatic posture adjustment in the automation process, must obtain the normal direction at processing stand place, then normal direction fed back to control system and carry out the attitude adjustment.Characteristics such as the present invention has just solved how to realize obtaining of normal direction, and has utilized the laser without contact measurement, and measuring accuracy is high, and is convenient can meet practical application request.
Claims (2)
 One kind based on 3 little planar processes to detection method, it is characterized in that: carry out as follows:1) set up model: 3 laser displacement sensor annulars are distributed in the head of three coordinate measuring machine, and the spatial attitude of sensor can make that laser is tapered to be penetrated on curved surface, controls 3 formed little plane sizes of laser point according to accuracy requirement adjustment tapering;2) set up space coordinates: at first, set up a world coordinate system 1 with three coordinate measuring machine; Select for use the terminal joint coordinate system of three coordinate measuring machine as intermediate conversion coordinate system 2 then; On each laser displacement sensor, setting up one in addition penetrates direction with laser rays and is z axle negative sense, is the local coordinate system 3 of true origin with laser displacement sensor work RP;3) image data: utilize data collecting card and LabVIEW software collection laser displacement sensor voltage signal and convert shift value L into _{i}4) calculate the coordinate of laser point: the displacement that utilizes laser displacement sensor to record promptly know the coordinate of laser point under local coordinate system 3 (0,0, w _{i}), w _{i}=L _{i}, according to the Eulerian coordinates conversion formula, try to achieve the coordinate of 3 laser point under coordinate system 1 respectively, the coordinate conversion formula is following:
${P}_{i}={R}_{0}*{R}_{i}*\left(\begin{array}{c}0\\ 0\\ {w}_{i}\\ 1\end{array}\right),$ i＝1，2，3${R}_{i}={R}_{z1}*{R}_{x2}+{T}_{i}$ $=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \mathrm{cos}\left({r}_{z1}\right)& \mathrm{sin}\left({r}_{z1}\right)& 0\\ 0& \mathrm{sin}\left({r}_{z1}\right)& \mathrm{cos}\left({r}_{z1}\right)& 0\\ 0& 0& 0& 1\end{array}\right]*\left[\begin{array}{cccc}\mathrm{cos}\left({r}_{x2}\right)& \mathrm{sin}\left({r}_{x2}\right)& 0& 0\\ \mathrm{sin}\left({r}_{x2}\right)& \mathrm{cos}\left({r}_{x2}\right)& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right]+\left[\begin{array}{cccc}0& 0& 0& {a}_{i}\\ 0& 0& 0& {b}_{i}\\ 0& 0& 0& {c}_{i}\\ 0& 0& 0& 0\end{array}\right]$ P in the formula _{i}The coordinate of expression laser point; R _{0}Be the terminal pose of three coordinate measuring machine, from three coordinate measuring machine software, obtain; R _{i}Be the homogeneous transformation matrix of local coordinate system 3 with respect to intermediate conversion coordinate system 2, wherein r _{Z1}, r _{X2}Be respectively coordinate system 3 around the z of coordinate system 2 axle, the x axle anglec of rotation, (a _{i}, b _{i}, c _{i}) be the translation vector of 3 pairs of coordinate systems 2 of coordinate system, to r _{Z1}, r _{X2}, (a _{i}, b _{i}, c _{i}) demarcate;5) computing method vector: 3 laser point coordinates of trying to achieve with step 4) can solve 3 vectors, and wherein any 2 vectorial multiplication crosses get final product to such an extent that 3 little planar processes are vectorial again, and formula is following:$\stackrel{\→}{{P}_{2}{P}_{1}}=({P}_{x1}{P}_{x2},{P}_{y1}{P}_{y2},{P}_{z1}{P}_{z2})$ $\stackrel{\→}{{P}_{2}{P}_{3}}=({P}_{x3}{P}_{x2},{P}_{y3}{P}_{y2},{P}_{z3}{P}_{z2})$ $\stackrel{\→}{N}=\stackrel{\→}{{P}_{2}{P}_{1}}\×\stackrel{\→}{{P}_{2}{P}_{3}}$ The normal vector that calculates is applied to normal error correction, the adjustment of normal direction attitude, normal direction detection.  2. method according to claim 1 is characterized in that, described laser displacement sensor attitude parameter scaling method: implementation step is following:(1) according to step 1) in the claim 1 laser displacement sensor is installed on the three dimensional coordinate measuring machine head and sets up following measurement space: according to step 2 in the claim 1) use three coordinate measuring machine to set up a world coordinate system 1 and set up an intermediate conversion coordinate system 2, in world coordinate system 1, draft a spherical equation (xa) then with standard ball ^{2}+ (yb) ^{2}+ (zc) ^{2}=r ^{2}, wherein, (a, b, c) expression sphere centre coordinate, r representes the radius of a ball, and sphere centre coordinate is demarcated with three coordinate measuring machine, and the radius of a ball is known;(2) on sphere, evenly gather 30 measuring points with one of them sensor operated three coordinate measuring machine, obtain 30 shift values;(3) the shift value substitution spherical equation that records is used fminsearch minimum optimization method and find the solution each parameter value;(4) repeating step (2), step (3) are found the solution the attitude parameter of other 2 sensors.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010281290 CN101957175B (en)  20100914  20100914  Threepoint microplanebased normal detection method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010281290 CN101957175B (en)  20100914  20100914  Threepoint microplanebased normal detection method 
Publications (2)
Publication Number  Publication Date 

CN101957175A CN101957175A (en)  20110126 
CN101957175B true CN101957175B (en)  20120328 
Family
ID=43484648
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 201010281290 CN101957175B (en)  20100914  20100914  Threepoint microplanebased normal detection method 
Country Status (1)
Country  Link 

CN (1)  CN101957175B (en) 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN108332909A (en) *  20180207  20180727  哈尔滨工业大学  A kind of multiple spot for centroid measurement is weighed the modification method of support construction model 
Families Citing this family (18)
Publication number  Priority date  Publication date  Assignee  Title 

CN102248450B (en) *  20110420  20130320  上海交通大学  Method for quickly detecting normal vector of largecurvatureradius curved surface 
CN102284956B (en) *  20110504  20131120  南京航空航天大学  Normal centering method of automatic drilling and riveting robot 
CN102322794B (en) *  20110512  20130605  浙江大学  Realtime radius compensation method and measuring head device based on LED (lightemitting diode) annular light 
CN102636110B (en) *  20120326  20150325  南京航空航天大学  Reference detecting device of automatic drilling and riveting system of airplane components and detecting method thereof 
CN102768006B (en) *  20120802  20150513  西北工业大学  Method for automatically measuring and adjusting largecurvature curve part in normal direction 
CN103419212A (en) *  20130520  20131204  北京航天金宇电子技术有限公司  Normal alignment method and device for arbitrary curved surface 
CN103307968B (en) *  20130603  20161102  西北工业大学  The detection method of robot carrying platform attitude 
CN103447877B (en) *  20130916  20150819  南京航空航天大学  Autonomous institutionalization is vowed and is detected and posture adjustment movement technique 
CN103528519B (en) *  20131021  20160817  中国航空工业集团公司北京航空精密机械研究所  A kind of measurement scaling method of some projection optics gauge head spatial position vector 
CN104615230B (en) *  20131101  20190426  腾讯科技（深圳）有限公司  A kind of operation recognition methods and terminal 
CN103968767B (en) *  20140514  20160629  哈尔滨工程大学  A kind of laser measurement method of spheroid identification 
CN105345600B (en) *  20151020  20200221  上海拓璞数控科技股份有限公司  Normal measurement and adjustment method for freeform surface drilling of aviation thinwalled part 
CN105783722B (en) *  20160412  20180406  南京航空航天大学  A kind of range sensor method arrow measuring method based on tetrahedral configuration 
CN106123817B (en) *  20160614  20181127  昆明理工大学  A kind of calibration system and method based on any space length 
CN107063086B (en) *  20170313  20190312  浙江优迈德智能装备有限公司  A kind of local surface normal vector based on double infrared sensors determines method 
CN107063089A (en) *  20170410  20170818  北京航空航天大学  A kind of surface normal measurement accuracy verification method 
CN108375352A (en) *  20180131  20180807  歌尔股份有限公司  The method of adjustment of two planpositions based on the depth of parallelism, apparatus and system 
CN110470223B (en) *  20190916  20200825  大连理工大学  Threedimensional space coordinate measuring method based on 1D displacement sensor 
Family Cites Families (5)
Publication number  Priority date  Publication date  Assignee  Title 

JP3466661B2 (en) *  19930629  20031117  キヤノン株式会社  Image processing apparatus and method 
KR100457080B1 (en) *  20020712  20041112  삼성중공업 주식회사  Method for surveying the characteristics of joint on rock slope using image 
DE10313038B4 (en) *  20030324  20050217  Klingelnberg Gmbh  Device for detecting the position of a probe element in a multicoordinate measuring device 
CN100453969C (en) *  20041231  20090121  上海萨克斯动力总成部件系统有限公司  Threecoordinate measuring method and apparatus for measuring complicated curved surface 
WO2008014461A2 (en) *  20060728  20080131  Optimet, Optical Metrology Ltd.  Doublesided measurement of dental objects using an optical scanner 

2010
 20100914 CN CN 201010281290 patent/CN101957175B/en not_active IP Right Cessation
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN108332909A (en) *  20180207  20180727  哈尔滨工业大学  A kind of multiple spot for centroid measurement is weighed the modification method of support construction model 
Also Published As
Publication number  Publication date 

CN101957175A (en)  20110126 
Similar Documents
Publication  Publication Date  Title 

CN103292701B (en)  The online dimension measurement method of accurate device based on machine vision  
CN103895023B (en)  A kind of tracking measurement method of the mechanical arm tail end tracing measurement system based on coding azimuth device  
Stoyanov et al.  Comparative evaluation of range sensor accuracy in indoor environments  
CN101672637B (en)  Digitizing detection method of complicated curved face  
CN105678785B (en)  A kind of laser and the scaling method of camera relative pose relation  
CN100453966C (en)  Spatial threedimensional position attitude measurement method for video camera  
CN102865862B (en)  Measurement method of nonhorizontalization free established station  
CN102460065B (en)  Information processing apparatus and information processing method  
CN104165750B (en)  Measurement method for pose of wind tunnel model combining stereoscopic vision with gyroscope  
CN101943563B (en)  Rapid calibration method of linestructured light vision sensor based on space plane restriction  
CN103198477B (en)  Apple fruitlet bagging robot visual positioning method  
CN102496232B (en)  Transmission facility monitoring method and system  
CN102297658B (en)  Threedimensional information detection method based on dual laser  
CN104219512B (en)  A kind of display device Gamut boundary description method  
CN104458895A (en)  Threedimensional pipeline leakage flux imaging detection method and system  
CN102679925B (en)  Robot localization error measurement method  
CN103278090B (en)  A kind of vision measuring method of irregularly shaped object volume  
CN102809364B (en)  Method for determining complex curved surface profile error  
CN102589429B (en)  Online engine cylinder hexahedral hole set position error detecting method  
CN103236064A (en)  Point cloud automatic registration method based on normal vector  
CN102779345B (en)  A kind of some cloud Precision Registration based on center of gravity Euclidean distance  
CN101871775A (en)  Programming system of three coordinate measuring machine and method thereof  
CN104656097B (en)  Caliberating device based on rotary twodimensional laser threedimensional reconfiguration system and method  
CN104677340B (en)  Point character based monocular vision pose measurement method  
CN103727930B (en)  A kind of laser range finder based on edge matching and camera relative pose scaling method 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
GR01  Patent grant  
C14  Grant of patent or utility model  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20120328 Termination date: 20150914 

EXPY  Termination of patent right or utility model 