CN110458822A - A kind of complex curved surface parts non-contact 3-D matching detection method - Google Patents

A kind of complex curved surface parts non-contact 3-D matching detection method Download PDF

Info

Publication number
CN110458822A
CN110458822A CN201910724278.6A CN201910724278A CN110458822A CN 110458822 A CN110458822 A CN 110458822A CN 201910724278 A CN201910724278 A CN 201910724278A CN 110458822 A CN110458822 A CN 110458822A
Authority
CN
China
Prior art keywords
point
curved surface
striation
measurement
reference line
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
CN201910724278.6A
Other languages
Chinese (zh)
Other versions
CN110458822B (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.)
Hefei College
Original Assignee
Hefei College
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 Hefei College filed Critical Hefei College
Priority to CN201910724278.6A priority Critical patent/CN110458822B/en
Publication of CN110458822A publication Critical patent/CN110458822A/en
Application granted granted Critical
Publication of CN110458822B publication Critical patent/CN110458822B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0004—Industrial image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/136—Segmentation; Edge detection involving thresholding
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/757—Matching configurations of points or features
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10004—Still image; Photographic image
    • G06T2207/10012—Stereo images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30108—Industrial image inspection
    • G06T2207/30164—Workpiece; Machine component

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Geometry (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Multimedia (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The present invention is suitable for curved surface part detection technique field, provides a kind of complex curved surface parts non-contact 3-D matching detection method, comprising the following steps: A, carries out three-dimensional one-point measurement to detection part；B, reference line is drawn；C, real-time region confirms；D, the integral face type profile of free form surface sample must be tested；E, 3 D stereo modeling is carried out to workpiece to be processed using mesh segmentation form to obtained curved surface part integral face type profile.Whereby, the present invention can be realized effectively and be accurately measured, and reduce measurement error, improve measurement accuracy.

Description

A kind of complex curved surface parts non-contact 3-D matching detection method
Technical field
The present invention relates to curved surface part detection technique field more particularly to a kind of complex curved surface parts non-contact 3-Ds With detection method.
Background technique
Currently, for complex curved surface parts detection there are two main classes the method for degree of precision: contact type measurement method and non-connecing Touch mensuration.Contact type measurement samples part using traditional three coordinate measuring machine point by point, and measurement accuracy is higher, but This method measurement efficiency is low, and has particular/special requirement to the material of tested part.With computer vision and mode identification technology Continuous development, non-contact measuring technology obtains with its higher measuring speed and precision in complex curved surface parts detection field It is more and more widely used.
Non-contact measurement obtains the scanning point set of different angle using optical scanner, then spells to these point sets It closes, obtains complete part scanning point set, believed by the error that scanning point set and the analysis of CAD model Point set matching obtain part Breath.In the detection of standard, general main consideration mismachining tolerance but has ignored influence of the measurement error to testing result；Thing In reality, but if the error for matching itself does not reach the ideal order of magnitude, the measurement result of mistake can be obtained, especially More it is preferably minimized matching error during the matching detection of high-precision complex curved surface parts.
In summary, the existing technology has inconveniences and defects in actual use, so it is necessary to be improved.
Summary of the invention
For above-mentioned defect, the purpose of the present invention is to provide a kind of matchings of complex curved surface parts non-contact 3-D to examine Survey method can reduce measurement error, improve measurement precision.
To achieve the goals above, the present invention provides a kind of complex curved surface parts non-contact 3-D matching detection method, The following steps are included:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, to curved surface part section to be measured in measurement process Selection fixation is carried out, is then gradually moved on the same axis, reference line measurement is carried out, it is mobile on the vertical axis after being measured to advise Gradually traverse measurement is carried out to another section after set a distance, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point Away from the uniform velocity carrying out surface scan to part to be measured by scanner, formed in scanning process complete after point is away from being provided with Striation real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image Striation location status be Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+dω/ F, i=0,1,2 ..., n, wherein Carrying out Threshold segmentation for image, treated that striation passes through The minimum component of region in the horizontal direction,For largest component, d is the average survey of scanner to piece surface to be measured Span is from ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, curved surface part progress X is tested to be obtained to straight-line motion accuracy when Y-direction Scanning Detction by laser interferometer measurement To displacement data compensate, by free form surface sample three-dimensional appearance data { D11 (x, y, z), D12 (x, y, z) ..., D12 (x, y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation is used to obtained curved surface part integral face type profile Form carries out 3 D stereo modeling to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, carries out the about fasciculation triangulation network Lattice subdivision extracts the skeleton of 2-d contour, chooses skeletal point and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce two Face angle principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
Complex curved surface parts non-contact 3-D matching detection method according to the present invention when the reference line measures, is surveyed Head is moved to B point from A point, C point is retracted into after B point has been surveyed, then by regulation step pitch to D point, repeat the measurement of next point E, surveyed It during amount, is pinpointed by angle of origin A, according to institute's measured data, draws curved surface part reference line.
Complex curved surface parts non-contact 3-D matching detection method according to the present invention, in the step C, using threshold value The image striation connected region being partitioned into searches for striation width value d by row (column) from top to bottom, obtains striation width value ordered series of numbers such as [d shown in loweri-λ..., di..., di+λ] (λ >=1, i=k, k+1 ..., k+n).
Complex curved surface parts non-contact 3-D matching detection method according to the present invention, the dkIt is wide for row k striation Angle value, k are the first trip of striation connected region, and n is total line number of striation connected region, with the width value d of every row iiAnd its up and down Each λ row element construction calculates ordered series of numbers [di-λ..., di..., di+λ] (λ >=1, i=k, k+1 ..., k+n), calculate its striation width Change rate ψi, with the width of the formal definition striation of discrete type variance, change rate is as follows
Wherein: pjFor jth row (j ∈ [i- λ, i+ λ]) striation width value djThe probability of appearance.
Complex curved surface parts non-contact 3-D matching detection method according to the present invention, the pjThe probability of appearance is pj =1/ (2 λ+1), taking λ=1, μ is the average value of each width element of calculating ordered series of numbers, and therefore, the striation width for obtaining the i-th row becomes Rate is
Wherein, the width value d of part is overflowed at i=1 and i=nk-1And dk+n+1It is handled by 0.
Complex curved surface parts non-contact 3-D matching detection method according to the present invention, the about fasciculation triangle gridding cut open The spatial topotaxy that GIS is utilized is divided to pre-process algorithm input data, the uniform data structure based on triangle is real Existing mesh refinement, reference line, which is drawn, is based on 2-d contour and first ball Modeling Technology using 2-d contour, utilizes central axes Geometrical property carries out data processing to it.
Complex curved surface parts non-contact 3-D matching detection method according to the present invention, the length variation of the reference line Formula: F=X2+Y2-R2, the deviation that starts to walk is preset or operation in quadrant transformation when according to quadrant where starting point and trend be to leave X-axis or the drift correction for carrying out F=F-X+Y or F=F-Y+X towards X-axis are preset, and to the recursion formula in curve motion: working as X When ± 1: F=F ± 2*X+1, when ± 1 Y: F=F ± 2*Y+1 is changed to: when ± 1 X: F=F ± 2*X+2, when ± 1 Y: F= F ± 2*Y+2, to make the correction amount to benchmark R realized in a quadrant by 0.8 through 0.6 to 1.5 variation, reach removal The maximum deviation of controlled point relative datum R in reference axis, keeps reference axis two sides track symmetrical, entire track is centered on R Operation.
The present invention provides a kind of complex curved surface parts non-contact 3-D matching detection methods, comprising the following steps:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, to curved surface part section to be measured in measurement process Selection fixation is carried out, is then gradually moved on the same axis, reference line measurement is carried out, it is mobile on the vertical axis after being measured to advise Gradually traverse measurement is carried out to another section after set a distance, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point Away from the uniform velocity carrying out surface scan to part to be measured by scanner, formed in scanning process complete after point is away from being provided with Striation real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image Striation location status be Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+dω/ F, i=0,1,2 ..., n, wherein Carrying out Threshold segmentation for image, treated that striation passes through The minimum component of region in the horizontal direction,For largest component, d is the average survey of scanner to piece surface to be measured Span is from ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, curved surface part progress X is tested to be obtained to straight-line motion accuracy when Y-direction Scanning Detction by laser interferometer measurement To displacement data compensate, by free form surface sample three-dimensional appearance data { D11 (x, y, z), D12 (x, y, z) ..., D12 (x, y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation is used to obtained curved surface part integral face type profile Form carries out 3 D stereo modeling to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, carries out the about fasciculation triangulation network Lattice subdivision extracts the skeleton of 2-d contour, chooses skeletal point and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce two Face angle principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
Beneficial effects of the present invention: the three-dimensional matching detection technique by designing complex curved surface parts can make full use of Its fast convergence rate, robustness are good and are not easy the characteristics of falling into local optimum, and test shows that high-precision and high efficiency can be obtained Three-dimensional matching result, using non-cpntact measurement, lesion element machined surface, not can determine that whether processing qualified by image, It can be used as online quick, batch detection effective means.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to embodiments, to the present invention It is further elaborated, it should be understood that the specific embodiments described herein are merely illustrative of the present invention, is not used to Limit the present invention.
The present invention provides a kind of complex curved surface parts non-contact 3-D matching detection methods, comprising the following steps:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, to curved surface part section to be measured in measurement process Selection fixation is carried out, is then gradually moved on the same axis, reference line measurement is carried out, it is mobile on the vertical axis after being measured to advise Gradually traverse measurement is carried out to another section after set a distance, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point Away from the uniform velocity carrying out surface scan to part to be measured by scanner, formed in scanning process complete after point is away from being provided with Striation real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image Striation location status be Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+dω/ F, i=0,1,2 ..., n, wherein Carrying out Threshold segmentation for image, treated that striation passes through The minimum component of region in the horizontal direction,For largest component, d is the average survey of scanner to piece surface to be measured Span is from ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, curved surface part progress X is tested to be obtained to straight-line motion accuracy when Y-direction Scanning Detction by laser interferometer measurement To displacement data compensate, by free form surface sample three-dimensional appearance data { D11 (x, y, z), D12 (x, y, z) ..., D12 (x, y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation is used to obtained curved surface part integral face type profile Form carries out 3 D stereo modeling to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, carries out the about fasciculation triangulation network Lattice subdivision extracts the skeleton of 2-d contour, chooses skeletal point and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce two Face angle principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
Preferably, when reference line of the invention measurement, gauge head is moved to B point from A point, is retracted into C after B point has been surveyed Point, then arrive D point by regulation step pitch, repeats the measurement of next point E, in measurement process, using origin A as angle fixed point, according to being surveyed Data draw curved surface part reference line, by repeatedly being moved based on origin, achieve the available point of reference line away from survey Amount, it is further to improve measurement precision.
In addition, the spatial topotaxy of GIS is utilized to algorithm input data in about fasciculation triangular mesh generation of the invention It is pre-processed, the uniform data structure based on triangle realizes mesh refinement, and reference line drafting is based on using 2-d contour 2-d contour and first ball Modeling Technology, carry out data processing to it using the geometrical property of central axes, utilize about fasciculation triangle Grid and 2-d contour and first ball Modeling Technology, to geometrical line away from effectively being measured, to guarantee in subsequent measurement process Precision.
Further, in step C of the invention, using Threshold segmentation go out image striation connected region, from top to bottom by Row (column) searches for striation width value d, obtains striation width value ordered series of numbers [d as followsi-λ..., di... di+λ] (λ >=1, i=k, K+1 ..., k+n), by thresholding method, ranging is split to the surface of curved surface part, and then guarantee its subsequent measurement Precision.
Preferably, the length variation formula of reference line of the invention: F=X2+Y2-R2, starting deviation it is preset or operation X-axis is left according to quadrant where starting point and trend when middle quadrant converts or carries out F=F-X+Y or F=F-Y+X towards X-axis Drift correction is preset, and to the recursion formula in curve motion: when ± 1 X: F=F ± 2*X+1, when ± 1 Y: F=F ± 2*Y + 1, it is changed to: when ± 1 X: F=F ± 2*X+2, when ± 1 Y: F=F ± 2*Y+2, to make the correction amount to benchmark R at one Realize in quadrant by 0.8 through 0.6 to 1.5 variation, reach the maximum deviation for removing the controlled point relative datum R in reference axis, Keep reference axis two sides track symmetrical, entire track is run centered on R, using length variation formula to the position of reference line and For point away from effective interpolation, and then the precision after guarantee measurement is carried out, raising is final to draw precision.
The present invention provides a kind of complex curved surface parts non-contact 3-D matching detection methods, comprising the following steps:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, to curved surface part section to be measured in measurement process Selection fixation is carried out, is then gradually moved on the same axis, reference line measurement is carried out, it is mobile on the vertical axis after being measured to advise Gradually traverse measurement is carried out to another section after set a distance, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point Away from the uniform velocity carrying out surface scan to part to be measured by scanner, formed in scanning process complete after point is away from being provided with Striation real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image Striation location status be Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+dω/ F, i=0,1,2 ..., n, wherein Carrying out Threshold segmentation for image, treated that striation passes through The minimum component of region in the horizontal direction,For largest component, d is the average survey of scanner to piece surface to be measured Span is from ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, curved surface part progress X is tested to be obtained to straight-line motion accuracy when Y-direction Scanning Detction by laser interferometer measurement To displacement data compensate, by free form surface sample three-dimensional appearance data { D11 (x, y, z), D12 (x, y, z) ..., D12 (x, y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation is used to obtained curved surface part integral face type profile Form carries out 3 D stereo modeling to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, carries out the about fasciculation triangulation network Lattice subdivision extracts the skeleton of 2-d contour, chooses skeletal point and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce two Face angle principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
When the reference line measures, gauge head is moved to B point from A point, C point is retracted into after B point has been surveyed, then by regulation step pitch To D point, the measurement of next point E is repeated, be that angle pinpoints using origin A, according to institute's measured data, drafting curved surface zero in measurement process Part reference line.In the step C, the image striation connected region gone out using Threshold segmentation searches for light by row (column) from top to bottom Width value d obtains striation width value ordered series of numbers [d as followsi-λ..., di... di+λ] (λ >=1, i=k, k+1 ..., k+n). The dkFor row k striation width value, k is the first trip of striation connected region, and n is total line number of striation connected region, with every row i Width value diAnd its each λ row element construction calculates ordered series of numbers [d up and downi-λ..., di..., di+λ] (λ >=1, i=k, k+1 ..., k+ N), its striation change width rate ψ is calculatedi, with the width of the formal definition striation of discrete type variance, change rate is as follows
Wherein: pjFor jth row (j ∈ [i- λ, i+ λ]) striation width value djThe probability of appearance.
The pjThe probability of appearance is pj=1/ (2 λ+1), taking λ=1, μ is the average value of each width element of calculating ordered series of numbers, Therefore, the striation change width rate for obtaining the i-th row is
Wherein, the width value d of part is overflowed at i=1 and i=nk-1And dk+n+1It is handled by 0.
The spatial topotaxy that about GIS is utilized in fasciculation triangular mesh generation pre-processes algorithm input data, base Mesh refinement is realized in the uniform data structure of triangle, and reference line, which is drawn, is based on 2-d contour and member using 2-d contour Ball Modeling Technology carries out data processing to it using the geometrical property of central axes.The length variation formula of reference line: F=X2+ Y2-R2, the deviation that starts to walk is preset or operation in quadrant transformation when according to quadrant where starting point and trend be to leave X-axis or towards X The drift correction that axis carries out F=F-X+Y or F=F-Y+X is preset, and to the recursion formula in curve motion: when ± 1 X: F=F ± 2*X+1, when ± 1 Y: F=F ± 2*Y+1 is changed to: when ± 1 X: F=F ± 2*X+2, and when ± 1 Y: F=F ± 2*Y+2, To make the correction amount to benchmark R realized in a quadrant by 0.8 through 0.6 to 1.5 variation, reach removal in reference axis The maximum deviation of controlled point relative datum R keeps reference axis two sides track symmetrical, and entire track is run centered on R.
In conclusion the present invention provides a kind of complex curved surface parts non-contact 3-D matching detection methods, including with Lower step:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, to curved surface part section to be measured in measurement process Selection fixation is carried out, is then gradually moved on the same axis, reference line measurement is carried out, it is mobile on the vertical axis after being measured to advise Gradually traverse measurement is carried out to another section after set a distance, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point Away from the uniform velocity carrying out surface scan to part to be measured by scanner, formed in scanning process complete after point is away from being provided with Striation real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image Striation location status be Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+dω/ F, i=0,1,2 ..., n, wherein Carrying out Threshold segmentation for image, treated that striation passes through The minimum component of region in the horizontal direction,For largest component, d is the average survey of scanner to piece surface to be measured Span is from ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, curved surface part progress X is tested to be obtained to straight-line motion accuracy when Y-direction Scanning Detction by laser interferometer measurement To displacement data compensate, by free form surface sample three-dimensional appearance data { D11 (x, y, z), D12 (x, y, z) ..., D12 (x, y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation is used to obtained curved surface part integral face type profile Form carries out 3 D stereo modeling to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, carries out the about fasciculation triangulation network Lattice subdivision extracts the skeleton of 2-d contour, chooses skeletal point and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce two Face angle principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
Beneficial effects of the present invention: the three-dimensional matching detection technique by designing complex curved surface parts can make full use of Its fast convergence rate, robustness are good and are not easy the characteristics of falling into local optimum, and test shows that high-precision and high efficiency can be obtained Three-dimensional matching result, using non-cpntact measurement, lesion element machined surface, not can determine that whether processing qualified by image, It can be used as online quick, batch detection effective means.
Certainly, the present invention can also have other various embodiments, without deviating from the spirit and substance of the present invention, ripe It knows those skilled in the art and makes various corresponding changes and modifications, but these corresponding changes and change in accordance with the present invention Shape all should fall within the scope of protection of the appended claims of the present invention.

Claims (7)

1. a kind of complex curved surface parts non-contact 3-D matching detection method, which comprises the following steps:
A, three-dimensional one-point measurement is carried out to detection part:
Benchmark line position is measured by Point Measurement method first, curved surface part section to be measured is carried out in measurement process Choose and fix, then gradually move on the same axis, carry out reference line measurement, after being measured on the vertical axis mobile regulation away from Gradually traverse measurement is carried out to another section from rear, forms reference line elementary contour；
B, reference line is drawn:
Reference line elementary contour formed in step A is drawn by mapping software, and sets corresponding point away from point After being provided with, surface scan is at the uniform velocity carried out to part to be measured by scanner, forms complete striation in scanning process Real-time region；
C, real-time region confirms:
Defined nucleotide sequence optical strip image { f ' (x, y, t0), f ' (x, y, t1) ..., f ' (x, y, tn-1) in i-th optical strip image light Location status is Xi, the striation positional relationship of adjacent optical strip image, which passes through following formula reasoning, indicates Xi+1=Xi+ d ω/f, i= 0,1,2 ..., n, wherein Threshold segmentation treated striation is carried out for image to exist by region Minimum component in horizontal direction,For largest component, d is average measurement distance of the scanner to piece surface to be measured, ω is the angular scanning speed of scanner, and f is the acquisition frame frequency of scanner；
D, it is tested curved surface part and carries out X to being obtained by laser interferometer measurement with straight-line motion accuracy when Y-direction Scanning Detction Displacement data compensates, by free form surface sample three-dimensional appearance data D11 (x, y, z), D12 (x, y, z) ..., D12 (x, Y, z), Dij (x, y, z) ..., DMN (x, y, z) } fitting, obtain the integral face type profile of tested free form surface sample；
E, on the reference line drawing basics of step B, mesh segmentation form is used to obtained curved surface part integral face type profile 3 D stereo modeling is carried out to workpiece to be processed, is distributed according to the characteristic point for drawing contour line, is carried out about fasciculation triangle gridding and cut open Point, the skeleton of 2-d contour is extracted, skeletal point is chosen and sampled point projects to three-dimensional space Ellipsoidal Surface, and introduce dihedral angle Principle optimizes the Triangulation Algorithm of spatial spreading data point, finally sutures skeletal point and obtains the three-dimensional grid representation of a surface.
2. complex curved surface parts non-contact 3-D matching detection method according to claim 1, which is characterized in that described When reference line measures, gauge head is moved to B point from A point, C point is retracted into after B point has been surveyed, then by regulation step pitch to D point, under repeating The measurement of one point E in measurement process, is pinpointed by angle of origin A, according to institute's measured data, draws curved surface part reference line.
3. complex curved surface parts non-contact 3-D matching detection method according to claim 1, which is characterized in that described In step C, the image striation connected region gone out using Threshold segmentation is searched for striation width value d by row (column) from top to bottom, obtained Striation width value ordered series of numbers [d as followsi-λ..., di... di+λ] (λ >=1, i=k, k+1 ..., k+n).
4. complex curved surface parts non-contact 3-D matching detection method according to claim 3, which is characterized in that described dkFor row k striation width value, k is the first trip of striation connected region, and n is total line number of striation connected region, with the width of every row i Angle value diAnd its each λ row element construction calculates ordered series of numbers [d up and downi-λ..., di..., di+λ] (λ >=1, i=k, k+1 ..., k+n), meter Calculate its striation change width rate ψi, with the width of the formal definition striation of discrete type variance, change rate is as follows
Wherein: pjFor jth row (j ∈ [i- λ, i+ λ]) striation width value djThe probability of appearance.
5. complex curved surface parts non-contact 3-D matching detection method according to claim 4, which is characterized in that described pjThe probability of appearance is pj=1/ (2 λ+1), taking λ=1, μ is therefore the average value of each width element of calculating ordered series of numbers obtains i-th Capable striation change width rate is
Wherein, the width value d of part is overflowed at i=1 and i=nk-1And dk+n+1It is handled by 0.
6. complex curved surface parts non-contact 3-D matching detection method according to claim 1, which is characterized in that described The spatial topotaxy that about GIS is utilized in fasciculation triangular mesh generation pre-processes algorithm input data, is based on triangle Uniform data structure realize mesh refinement, reference line, which is drawn, is based on 2-d contour and first ball moulding skill using 2-d contour Art carries out data processing to it using the geometrical property of central axes.
7. complex curved surface parts non-contact 3-D matching detection method according to claim 1, which is characterized in that described The length variation formula of reference line: F=X2+Y2-R2, start to walk deviation it is preset or operation in quadrant transformation when according to where starting point Quadrant and trend are that the drift correction leaving X-axis or carry out F=F-X+Y or F=F-Y+X towards X-axis is preset, and to curve motion In recursion formula: when ± 1 X: F=F ± 2*X+1, when ± 1 Y: F=F ± 2*Y+1 is changed to: when ± 1 X: F=F ± 2* X+2, when ± 1 Y: F=F ± 2*Y+2, realize the correction amount to benchmark R by 0.8 in a quadrant through 0.6 to 1.5 Variation, reach remove in reference axis controlled point relative datum R maximum deviation, keep reference axis two sides track symmetrical, entirely Track is run centered on R.
CN201910724278.6A 2019-08-07 2019-08-07 Non-contact three-dimensional matching detection method for complex curved surface part Active CN110458822B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910724278.6A CN110458822B (en) 2019-08-07 2019-08-07 Non-contact three-dimensional matching detection method for complex curved surface part

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910724278.6A CN110458822B (en) 2019-08-07 2019-08-07 Non-contact three-dimensional matching detection method for complex curved surface part

Publications (2)

Publication Number Publication Date
CN110458822A true CN110458822A (en) 2019-11-15
CN110458822B CN110458822B (en) 2022-10-11

Family

ID=68485106

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910724278.6A Active CN110458822B (en) 2019-08-07 2019-08-07 Non-contact three-dimensional matching detection method for complex curved surface part

Country Status (1)

Country Link
CN (1) CN110458822B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118518027A (en) * 2024-07-24 2024-08-20 中国汽车技术研究中心有限公司 Device and method for drawing complex parts of automobile collision dummy

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004079295A2 (en) * 2003-03-06 2004-09-16 Zygo Corporation Profiling complex surface structures using scanning interferometry
CN102305601A (en) * 2011-05-18 2012-01-04 天津大学 High-precision non-contact measurement method and device for three-dimensional profile of optical freeform curved surface
CN104484508A (en) * 2014-11-26 2015-04-01 华中科技大学 Optimizing method for noncontact three-dimensional matching detection of complex curved-surface part
CN105627948A (en) * 2016-01-31 2016-06-01 山东科技大学 Large-scale complex curved surface measurement system and application thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004079295A2 (en) * 2003-03-06 2004-09-16 Zygo Corporation Profiling complex surface structures using scanning interferometry
CN102305601A (en) * 2011-05-18 2012-01-04 天津大学 High-precision non-contact measurement method and device for three-dimensional profile of optical freeform curved surface
CN104484508A (en) * 2014-11-26 2015-04-01 华中科技大学 Optimizing method for noncontact three-dimensional matching detection of complex curved-surface part
CN105627948A (en) * 2016-01-31 2016-06-01 山东科技大学 Large-scale complex curved surface measurement system and application thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杜娟: "应用三坐标测量机检测刀片模具复杂三维曲面形状精度的方法", 《工具技术》 *
闫玉蕾: "基于三维扫描的复杂曲面在机测量系统的研究", 《机械工程师》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118518027A (en) * 2024-07-24 2024-08-20 中国汽车技术研究中心有限公司 Device and method for drawing complex parts of automobile collision dummy

Also Published As

Publication number Publication date
CN110458822B (en) 2022-10-11

Similar Documents

Publication Publication Date Title
CN106643551B (en) A kind of quick scanning means of blade shape and method
CN110014426B (en) Method for grabbing symmetrically-shaped workpieces at high precision by using low-precision depth camera
Lee et al. Data reduction methods for reverse engineering
CN102494657B (en) Measuring head radius compensation method for curve surface profile measuring and detecting
CN112348864B (en) Three-dimensional point cloud automatic registration method for laser contour features of fusion line
CN105654483B (en) The full-automatic method for registering of three-dimensional point cloud
CN107301648B (en) Redundant point cloud removing method based on overlapping area boundary angle
CN110455187B (en) Three-dimensional vision-based box workpiece weld joint detection method
CN104392476A (en) Method of extracting three-dimensional axis of tunnel based on minimum bounding box algorithm
Wu et al. A novel high precise laser 3D profile scanning method with flexible calibration
CN116402866A (en) Point cloud-based part digital twin geometric modeling and error assessment method and system
CN112233249A (en) B spline surface fitting method and device based on dense point cloud
CN108682043A (en) A kind of complex-curved measure planning method based on parameter mapping
CN112033338B (en) Blade curved surface contact type scanning measurement probe radius surface compensation method
CN109448042B (en) Method for calculating thickness of blood vessel wall by using integral curve
Ravishankar et al. Automated inspection of aircraft parts using a modified ICP algorithm
CN109919955A (en) The tunnel axis of ground formula laser radar point cloud extracts and dividing method
CN110458822A (en) A kind of complex curved surface parts non-contact 3-D matching detection method
CN113267122B (en) Industrial part size measurement method based on 3D vision sensor
CN112581521B (en) Method for extracting central line of magnetic suspension track
CN115790440A (en) Profile tolerance measuring method and measuring system based on 3D scanning
CN106407603B (en) A kind of method of virtual gear measuring center last time Rotary Workpiece Fast Collision Detection
CN112365534B (en) Large coal pile volume measurement method based on monocular camera three-dimensional reconstruction
CN112762790B (en) Method for measuring crown width in three-dimensional virtual scene
CN115222893A (en) Three-dimensional reconstruction splicing method for large-size components based on structured light measurement

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