A kind of diffraction light-free Moire fringe center positioning method and system
Technical field
The present invention relates to diffraction light-free fields of measurement, more particularly to form the accurate of Moire fringe in two beam diffraction light-frees and survey
In amount, two location of the core methods.
Background technique
Diffraction light-free is concerned since its center spot size, shape do not change with transmission range.With people
Deepen continuously to diffraction light-free research, diffraction light-free application field it is also more and more extensive, as accurate measurement, laser add
Work, field depth extending etc..The accurate position for finding diffraction light-free center spot is an important prerequisite for realizing high measurement accuracy.Often
The diffraction light-free Facula Center Location method seen has light intensity gravity model appoach, Least Square Circle fitting process and Hough circle transformation method, annulus
Filter method and chaos ant colony method etc..
Since diffraction light-free has good collimation and the medium characteristic of energy collection, by diffraction light-free in conjunction with other technologies,
The application field of diffraction light-free can be greatly enriched.It is one in accurate measurement that diffraction light-free is combined with Moire fringe technology
The advantages of new combined method of kind, this method not only remains Moire fringe, but also introduce the characteristic of diffraction light-free.It will be without spreading out
It penetrates light to combine with Moire fringe, have been widely used in field of precision measurement, as space linearity measurement, multifreedom motion miss
Difference measurements etc..Such as in multifreedom motion error measure, two beam diffraction light-frees are subjected to interference and generate diffraction light-free More item
Line, four kinds of different freedom degrees variations that you can get it by calculating the location variation of two spot center of diffraction light-free Moire fringe
Amount.Therefore, the position for how accurately finding two center spot of diffraction light-free Moire fringe is to realize diffraction light-free Moire fringe essence
The premise of close measurement.
The positioning at two centers of diffraction light-free Moire fringe, is different from single beam diffraction light-free, and diffraction light-free Moire fringe exists
Radially fringe distribution on shape feature, there are the regions that two light intensity are most concentrated in light distribution.Because of diffraction light-free
The complexity of Moire fringe, so diffraction light-free Moire fringe positioning difficulty for diffraction light-free positioning is bigger.At present
Mainly have to the method for diffraction light-free Moire fringe positioning: (1) blocking method.Block the first beam salt free ligands first during the experiment
Light calculates the spot center of the second beam diffraction light-free;Then the second beam diffraction light-free is blocked, the light of the first beam diffraction light-free is calculated
Spot center;The positioning of diffraction light-free Moire fringe is realized in the position for finally combining two groups of spot centers.The method will be to diffraction light-free
Moire fringe positioning, which is converted into, positions diffraction light-free, positioning accuracy with higher, but needs respectively to a branch of salt free ligands
Light carries out mechanization and blocks, and cannot calculate two spot center change in location of diffraction light-free Moire fringe, the degree of automation simultaneously
It is not high.(2) the Schlieren method.Concentric circles Moire fringe have it is preferable Targets, theory analysis it is found that two raster centers away from every increase and decrease
One pitch, Moire fringe number increase and decrease four.Equally, diffraction light-free Moire fringe also has this rule.Not to diffraction light-free
That striped is scanned counting in rectangular coordinate system or polar coordinate system, pushes away two light of hot spot by the way that the variable quantity of Moire fringe number is counter
The variable quantity at spot center realizes the positioning of diffraction light-free Moire fringe.The method needs to carry out accurate metering to Moire fringe, still
Since the quality of Moire fringe is very big on count results influence, so Facula Center Location precision is not high.
In short, currently used diffraction light-free Moire fringe localization method is in positioning accuracy, the degree of automation and is applicable in
The problem of range etc. is further studied there are still many needs.
Summary of the invention
The object of the present invention is to provide a kind of diffraction light-free Moire fringe center positioning methods, can be realized diffraction light-free not
The accurate positioning of your two central points of striped.This method can simultaneously position two spot center of diffraction light-free Moire fringe,
The influence that the method for blocking causes the degree of automation low is avoided, the problem that the Schlieren method causes positioning accuracy not high is also reduced.It completes
The technical solution of foregoing invention task is: a kind of diffraction light-free Moire fringe center positioning method, comprising the following steps:
Step 1, pretreatment operation is carried out to noise-containing diffraction light-free Moire fringe image, extracts local concentric loop
Region;
Step 2, for the image of step 1 processing, a series of Morphological scale-spaces are carried out;
Step 3, for the image of step 2 processing, local concentric loop detection is carried out, initial center of circle collection C is obtained0(x,y);
Step 4, for the image of step 3 processing, to initial center of circle collection C0(x, y) carries out clustering, flat using error
The objective function of side and SSE as cluster, when SSE minimum, Clustering Effect is best, to obtain in two beam diffraction light-free hot spots
Heart center of circle collection C1(x, y) and C2(x, y), wherein the calculation formula of SSE be,
Wherein CiIt is i-th of cluster, p is CiIn sample point, diIt is CiIn all samples mean value;
To the two groups of center of circle collection C detected1(x, y) and C2(x, y) seeks its central point O1(a, b) and O2(a, b) is as just
Beginning center;Then using initial center as the center of circle, certain length R is radius, and radius R is formed by circle will be comprising respective center of circle collection
Concentrated part deletes point, that is, abnormal point other than circle;Central point is sought to new center of circle collection again, is obtained for iteration w times
O3(a, b) and O4(a, b), as more accurate two beams diffraction light-free center, so that it is determined that diffraction light-free Moire fringe position.
Further, pretreatment operation includes following sub-step in step 1,
Step 1.1, it is denoised using median filtering;The neighborhood that size is S is chosen, with the big minispread neighborhood of intensity value
Interior pixel value selects new value of the median in arrangement pixel set as abnormal point (m, n), while realizing denoising effect
Retain the marginal information of image;
Step 1.2, to treated, image carries out histogram equalization;By stretching image pixel intensity range enhancing figure
Tonal range under different experiments environment is balanced to 0~255 by image contrast, weakens influence of the external environment to experimental result;
Step 1.3, to treated, image carries out Threshold segmentation and binaryzation;Light intensity maximum value is found in whole picture figure again
Imax, according to salt free ligands light intensity distribution character under experimental situation, suitable Δ V value is chosen, chooses ImaxΔ V is as threshold value to image
Binary conversion treatment is carried out, obtains diffraction light-free close to its spot center part concentric annular regions.
Further, firstly, the local concentric loop in two spot centers for not influencing diffraction light-free Moire fringe is whole
Under the premise of profile, deletes area and be less than S1The connected region of pixel, the image obtained at this time are I1;Then use size for P's
Disc-shaped structure element I2To image I1The opening operation expanded afterwards, the calculation formula of opening operation are first corroded in progress are as follows:
Wherein $ is etching operation,For expansive working;It finally deletes again because the area that opening operation generates image is less than
S2The new connected region of pixel.
Further, the specific implementation of step 3 is as follows,
Edge detection is carried out using Sobel operator first, enabling by the image that step 2 is handled is I3, size is selected respectively
For H, direction is horizontal kernel GxIt is vertical kernel G with directionyConvolution algorithm is carried out with image, then in the every bit of image
Gray scale G size can be approximated to be | I3*Gx|+|I3*Gy|;Then Hough circle transformation is used, image space is transformed into parameter
Suitable radius step-length step_r, angle step step_angle, minimum radius of circle r_min, maximum radius of circle r_ are chosen in space
Max carries out statistics calculating to the point on edge, determines circle that may be present, obtains initial center of circle collection C0(x,y)。
Further, the concentric loop that diffraction light-free spot center part concentric annular regions are included in step 1.3
Number is 3~6.
Further, use K mean cluster algorithm to initial center of circle collection C in step 40(x, y) carries out clustering, and K takes
2。
Further, it will include the concentrated part of respective center of circle collection containing the respective center of circle that radius R, which is formed by circle, in step 4
85% pixel integrated is standard.
Further, the number of iterations w should be no more than 3 times in step 4.
The present invention also provides a kind of diffraction light-free Moire fringe center location systems, comprise the following modules:
Local concentric annular regions extraction module, for being located in advance to noise-containing diffraction light-free Moire fringe image
Local concentric annular regions are extracted in reason operation;
Morphological scale-space module carries out a series of for the image for local concentric annular regions extraction module processing
Morphological scale-space;
Local concentric loop detection module carries out local concentric circles for being directed to the image of Morphological scale-space resume module
Ring detection, obtains initial center of circle collection C0(x,y);
Moire fringe locating module, for the image for local concentric loop detection module processing, to initial center of circle collection
C0(x, y) carries out clustering, uses error sum of squares SSE as the objective function of cluster, when SSE minimum, Clustering Effect
Preferably, to obtain two beam diffraction light-free spot center center of circle collection C1(x, y) and C2(x, y), wherein the calculation formula of SSE be,
Wherein CiIt is i-th of cluster, p is CiIn sample point, diIt is CiIn all samples mean value;
To the two groups of center of circle collection C detected1(x, y) and C2(x, y) seeks its central point O1(a, b) and O2(a, b) is as just
Beginning center;Then using initial center as the center of circle, certain length R is radius, and radius R is formed by circle will be comprising respective center of circle collection
Concentrated part deletes point, that is, abnormal point other than circle;Central point is sought to new center of circle collection again, is obtained for iteration w times
O3(a, b) and O4(a, b), as more accurate two beams diffraction light-free center, so that it is determined that diffraction light-free Moire fringe position.
Further, the specific process flow of local concentric annular regions extraction module is as follows,
Step 1.1, it is denoised using median filtering;The neighborhood that size is S is chosen, with the big minispread neighborhood of intensity value
Interior pixel value selects new value of the median in arrangement pixel set as abnormal point (m, n), while realizing denoising effect
Retain the marginal information of image;
Step 1.2, to treated, image carries out histogram equalization;By stretching image pixel intensity range enhancing figure
Tonal range under different experiments environment is balanced to 0~255 by image contrast, weakens influence of the external environment to experimental result;
Step 1.3, to treated, image carries out Threshold segmentation and binaryzation;Light intensity maximum value is found in whole picture figure again
Imax, according to salt free ligands light intensity distribution character under experimental situation, suitable Δ V value is chosen, chooses ImaxΔ V is as threshold value to image
Binary conversion treatment is carried out, obtains diffraction light-free close to its spot center part concentric annular regions.
The specific process flow of Morphological scale-space module is as follows,
Firstly, the premise of the local concentric loop overall profile in two spot centers for not influencing diffraction light-free Moire fringe
Under, it deletes area and is less than S1The connected region of pixel, the image obtained at this time are I1;Then use size for the disc-shaped structure of P
Element I2To image I1The opening operation expanded afterwards, the calculation formula of opening operation are first corroded in progress are as follows:
Wherein $ is etching operation,For expansive working;It finally deletes again because the area that opening operation generates image is less than
S2The new connected region of pixel.
Local concentric loop detection module detailed process processing is as follows,
Edge detection is carried out using Sobel operator first, enabling by the image that step 2 is handled is I3, size is selected respectively
For H, direction is horizontal kernel GxIt is vertical kernel G with directionyConvolution algorithm is carried out with image, then in the every bit of image
Gray scale G size can be approximated to be | I3*Gx|+|I3*Gy|;Then Hough circle transformation is used, image space is transformed into parameter
Suitable radius step-length step_r, angle step step_angle, minimum radius of circle r_min, maximum radius of circle r_ are chosen in space
Max carries out statistics calculating to the point on edge, determines circle that may be present, obtains initial center of circle collection C0(x,y)。
Compared with prior art, the method for the present invention has the advantage that
(1) image under different experiments environment, strong applicability can preferably be handled.
(2) retain the local concentric loop of two spot centers close to diffraction light-free Moire fringe, the speed of service is fast.
(3) two spot center of diffraction light-free Moire fringe can be positioned simultaneously.The method of blocking is avoided to cause to automate
The low influence of degree also reduces the problem that the Schlieren method causes positioning accuracy not high.
Detailed description of the invention
Fig. 1: process blocks schematic diagram of the invention;
Fig. 2: the noise-containing diffraction light-free Moire fringe image of emulation of the invention;
Fig. 3: the diffraction light-free Moire fringe image after image preprocessing of the invention;
Fig. 4: the Moire fringe image after K mean cluster of the present invention;
Fig. 5: the diffraction light-free Moire fringe image after center of circle exception point deletion of the invention;
Fig. 6: the diffraction light-free Moire fringe image of two center spots positioning of the invention;
Specific embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, with reference to the accompanying drawings and embodiments to this hair
It is bright to be described in further detail, it should be understood that implementation example described herein is merely to illustrate and explain the present invention, not
For limiting the present invention.
Below with reference to Fig. 1, the present invention will be further described:
Step 1, pretreatment operation, extraction office are carried out to noise-containing diffraction light-free Moire fringe image shown in Fig. 2
Portion concentric annular regions, as shown in Figure 3.Specifically:
Step 1.1, it is denoised using median filtering.Selection size S is 3 × 3 neighborhoods, adjacent with the big minispread of intensity value
Pixel value in domain selects new value of the median in arrangement pixel set as point A (m, n).The pixel value of postulated point A is
87, and be outlier pixel values, 3 × 3 neighborhoods are [158 155 151;157 87 156;152 153 150], in this neighborhood
Interior, normal pixel value is 155 or so.After median filtering, the new pixel value of point A is 153, thus obtained new pixel value with
Normal pixel value size is close, and the marginal information of image is retained while realizing denoising effect.
Step 1.2, to treated, image carries out histogram equalization.Image under this experimental situation is united by histogram
Meter gray value of image range is mainly distributed on 0~125 section.Enhance picture contrast using image pixel intensity range is stretched,
Tonal range under this experimental situation is balanced to 0~255, weakens influence of the external environment to experimental result.
Step 1.3, to treated, image carries out Threshold segmentation and binaryzation.Light intensity maximum value is found in whole picture figure
Imax=255, according to salt free ligands light intensity distribution character under this experimental situation, Δ V=195 is chosen, chooses ImaxΔ V=60 conduct
Threshold value carries out binary conversion treatment to image, so that the concentric loop that diffraction light-free spot center part concentric annular regions are included
With 3~6 numbers for standard.
Step 2, for the image of step 1 processing, a series of Morphological scale-spaces are carried out: firstly, not influencing diffraction light-free
Under the premise of the local concentric loop overall profile of two spot centers of Moire fringe, deletes area and be less than S1=100 pixels
Connected region, improves positioning accuracy, and the image obtained at this time is I1.Then the disc-shaped structure element I of size P=1 is used2It is right
Image I1The opening operation expanded afterwards, the calculation formula of opening operation are first corroded in progress are as follows:
Wherein $ is etching operation,For expansive working.While the effect of opening operation guarantees using opening operation smoothed profile
Keep profile close to circle.The area that finally deleting again generates image because of opening operation is less than S2The new connected region of=60 pixels
Domain.
Step 3, for the image of step 2 processing, local concentric loop detection is carried out.It is carried out first using Sobel operator
Edge detection.Enabling by the image that step 2 is handled is I3, selecting size respectively is H=3, and direction is horizontal kernel Gx=[- 1
0 1;-20 2;- 10 1] and direction is vertical kernel Gy=[- 1-2-1;0 0 0;12 1] convolution fortune is carried out with image
Calculate, then can be approximated to be in the gray scale G size of the every bit of image | I3*Gx|+|I3*Gy|.Then Hough circle transformation is used, it will
Image space is transformed into parameter space, chooses radius step-length step_r=1, angle step step_angle=0.1, smallest circle half
Diameter r_min=2, maximum radius of circle r_max=80, carry out statistics calculating to the point on edge, determine circle that may be present, obtain
Initial center of circle collection C0(x,y)。
Step 4, for the image of step 3 processing, to initial center of circle collection C0(x, y) carries out clustering.K mean cluster category
In a kind of unsupervised machine learning, similar object is grouped into same cluster.It is by two beams according to diffraction light-free Moire fringe
Diffraction light-free is interfered, and takes K=2 that can preferably find two beam diffraction light-free spot center center of circle collection, as shown in Figure 4.Make
Use error sum of squares SSE as the objective function of cluster, when SSE minimum, Clustering Effect is best.The calculation formula of SSE is
Wherein CiIt is i-th of cluster, p is CiIn sample point, diIt is CiIn all samples mean value.To two groups detected
Center of circle collection C1(x, y) and C2(x, y) seeks its central point O1(131,128) and O2(190,188) it is used as initial center.Then with
Initial center is the center of circle, and certain length R=10 is radius, and it will include the 85% of respective center of circle collection that radius R=10, which is formed by circle,
Pixel is standard, deletes point, that is, abnormal point other than circle, as shown in Figure 5.Center is finally sought again to new center of circle collection
Point, obtains O iteration 2 times3(128,128) and O4(191,190), as more accurate two beams diffraction light-free center,
So that it is determined that diffraction light-free Moire fringe position, as shown in Figure 6.
The present invention also provides a kind of diffraction light-free Moire fringe center location systems, comprise the following modules:
Local concentric annular regions extraction module, for being located in advance to noise-containing diffraction light-free Moire fringe image
Local concentric annular regions are extracted in reason operation;
Morphological scale-space module carries out a series of for the image for local concentric annular regions extraction module processing
Morphological scale-space;
Local concentric loop detection module carries out local concentric circles for being directed to the image of Morphological scale-space resume module
Ring detection, obtains initial center of circle collection C0(x,y);
Moire fringe locating module, for the image for local concentric loop detection module processing, to initial center of circle collection
C0(x, y) carries out clustering, uses error sum of squares SSE as the objective function of cluster, when SSE minimum, Clustering Effect
Preferably, to obtain two beam diffraction light-free spot center center of circle collection C1(x, y) and C2(x, y), wherein the calculation formula of SSE be,
Wherein CiIt is i-th of cluster, p is CiIn sample point, diIt is CiIn all samples mean value;
To the two groups of center of circle collection C detected1(x, y) and C2(x, y) seeks its central point O1(a, b) and O2(a, b) is as just
Beginning center;Then using initial center as the center of circle, certain length R is radius, and radius R is formed by circle will be comprising respective center of circle collection
Concentrated part deletes point, that is, abnormal point other than circle;Central point is sought to new center of circle collection again, is obtained for iteration w times
O3(a, b) and O4(a, b), as more accurate two beams diffraction light-free center, so that it is determined that diffraction light-free Moire fringe position.
The specific implementation of each module and each step is corresponding, and the embodiment of the present invention is not write.
Above in conjunction with attached drawing, an exemplary description of the invention, it is clear that realization of the invention is not by aforesaid way
Limitation.
It should be understood that the part that this specification does not elaborate belongs to the prior art.
It should be understood that the above-mentioned description for preferred embodiment is more detailed, can not therefore be considered to this
The limitation of invention patent protection range, those skilled in the art under the inspiration of the present invention, are not departing from power of the present invention
Benefit requires to make replacement or deformation under protected ambit, fall within the scope of protection of the present invention, this hair
It is bright range is claimed to be determined by the appended claims.