CN109523565A - A kind of diffraction light-free Moire fringe center positioning method and system - Google Patents
A kind of diffraction light-free Moire fringe center positioning method and system Download PDFInfo
- Publication number
- CN109523565A CN109523565A CN201811360167.3A CN201811360167A CN109523565A CN 109523565 A CN109523565 A CN 109523565A CN 201811360167 A CN201811360167 A CN 201811360167A CN 109523565 A CN109523565 A CN 109523565A
- Authority
- CN
- China
- Prior art keywords
- image
- center
- diffraction light
- circle
- free
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 230000000877 morphologic effect Effects 0.000 claims abstract description 14
- 230000002159 abnormal effect Effects 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 230000000694 effects Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 8
- 239000003446 ligand Substances 0.000 claims description 7
- 150000003839 salts Chemical class 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000003708 edge detection Methods 0.000 claims description 5
- 238000005530 etching Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000007689 inspection Methods 0.000 claims 1
- 238000007781 pre-processing Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- 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/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G06T5/70—
-
- 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/12—Edge-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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- 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/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
Abstract
The present invention relates to a kind of diffraction light-free Moire fringe center positioning method and systems.Method includes the following steps: Step 1: carrying out image preprocessing, the local concentric annular regions of extraction to noise-containing diffraction light-free Moire fringe；Step 2: carrying out a series of Morphological scale-spaces for the image of step 1 processing；Step 3: carrying out local concentric loop detection for the image of step 2 processing；Step 4: carrying out clustering to initial center of circle collection for the image of step 3 processing, the abnormal point after clustering in two groups of center of circle collection is deleted, the positioning of diffraction light-free Moire fringe can be realized in the central point that iteration seeks two groups of center of circle collection.The present invention can simultaneously position two spot center of diffraction light-free Moire fringe, and the degree of automation is strong, and positioning accuracy is high, applied widely, there is good application prospect.
Description
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.
Claims (10)
1. a kind of diffraction light-free Moire fringe center positioning method, which comprises the following steps:
Step 1, pretreatment operation is carried out to noise-containing diffraction light-free Moire fringe image, extracts local concentric circles ring region
Domain；
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, uses error sum of squares
Objective function of the SSE as cluster, when SSE minimum, Clustering Effect is best, to obtain two beam diffraction light-free spot centers circle
Heart 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 used as in initial
The heart；Then using initial center as the center of circle, certain length R is radius, and radius R is formed by the concentration that circle will include respective center of circle collection
Point, that is, abnormal point other than circle is deleted in part；Central point is sought to new center of circle collection again, obtains O iteration w times3(a,
And O b)4(a, b), as more accurate two beams diffraction light-free center, so that it is determined that diffraction light-free Moire fringe position.
2. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: pre-processed in step 1
Operation includes following sub-step,
Step 1.1, it is denoised using median filtering；The neighborhood that size is S is chosen, in the big minispread neighborhood of intensity value
Pixel value selects new value of the median in arrangement pixel set as abnormal point (m, n), retains while realizing denoising effect
The marginal information of image；
Step 1.2, to treated, image carries out histogram equalization；Enhance image pair by stretching image pixel intensity range
Than degree, the tonal range under different experiments environment 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 I is found in whole picture figure againmax,
According to salt free ligands light intensity distribution character under experimental situation, suitable Δ V value is chosen, chooses ImaxΔ V carries out image as threshold value
Binary conversion treatment obtains diffraction light-free close to its spot center part concentric annular regions.
3. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: the specific reality of step 2
Existing mode is as follows,
Firstly, under the premise of not influencing the local concentric loop overall profile of two spot centers of diffraction light-free Moire fringe,
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 member
Plain 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, and ⊕ is expansive working；The area that finally deleting again generates image because of opening operation is less than S2Picture
The new connected region of element.
4. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: the specific reality of step 3
Existing mode is as follows,
Edge detection is carried out using Sobel operator first, enabling by the image that step 2 is handled is I3, selecting size respectively is H, side
To for horizontal kernel GxIt is vertical kernel G with directionyConvolution algorithm is carried out with image, then in the gray scale G of the every bit of image
Size can be approximated to be | I3*Gx|+|I3*Gy|；Then Hough circle transformation is used, image space is transformed into parameter space, is selected
Suitable radius step-length step_r, angle step step_angle, minimum radius of circle r_min, maximum radius of circle r_max are taken, it is right
Point on edge carries out statistics calculating, determines circle that may be present, obtains initial center of circle collection C0(x,y)。
5. diffraction light-free Moire fringe center positioning method as claimed in claim 2, it is characterised in that:
The concentric loop number that diffraction light-free spot center part concentric annular regions are included in step 1.3 is 3~6.
6. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: equal using K in step 4
It is worth clustering algorithm to initial center of circle collection C0(x, y) carries out clustering, and K takes 2.
7. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: radius R institute in step 4
The circle of formation will include the concentrated part of respective center of circle collection using 85% pixel integrated containing the respective center of circle as standard.
8. diffraction light-free Moire fringe center positioning method as described in claim 1, it is characterised in that: iteration time in step 4
Number w should be no more than 3 times.
9. a kind of diffraction light-free Moire fringe center location system, which is characterized in that comprise the following modules:
Local concentric annular regions extraction module, for carrying out pretreatment behaviour to noise-containing diffraction light-free Moire fringe image
Make, extracts local concentric annular regions；
Morphological scale-space module carries out a series of forms for the image for local concentric annular regions extraction module processing
Processing；
Local concentric loop detection module carries out local concentric loop inspection for being directed to the image of Morphological scale-space resume module
It surveys, 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) clustering is carried out, uses error sum of squares SSE as the objective function of cluster, when SSE minimum, Clustering Effect is best,
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 used as in initial
The heart；Then using initial center as the center of circle, certain length R is radius, and radius R is formed by the concentration that circle will include respective center of circle collection
Point, that is, abnormal point other than circle is deleted in part；Central point is sought to new center of circle collection again, obtains O iteration w times3(a,
And O b)4(a, b), as more accurate two beams diffraction light-free center, so that it is determined that diffraction light-free Moire fringe position.
10. diffraction light-free Moire fringe center location system as claimed in claim 9, it is characterised in that: local concentric loop
The specific process flow of region extraction module is as follows,
Step 1.1, it is denoised using median filtering；The neighborhood that size is S is chosen, in the big minispread neighborhood of intensity value
Pixel value selects new value of the median in arrangement pixel set as abnormal point (m, n), retains while realizing denoising effect
The marginal information of image；
Step 1.2, to treated, image carries out histogram equalization；Enhance image pair by stretching image pixel intensity range
Than degree, the tonal range under different experiments environment 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 I is found in whole picture figure againmax,
According to salt free ligands light intensity distribution character under experimental situation, suitable Δ V value is chosen, chooses ImaxΔ V carries out image as threshold value
Binary conversion treatment 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, under the premise of not influencing the local concentric loop overall profile of two spot centers of diffraction light-free Moire fringe,
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 member
Plain 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, and ⊕ is expansive working；The area that finally deleting again generates image because of opening operation is less than S2Picture
The new connected region of element；
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, selecting size respectively is H, side
To for horizontal kernel GxIt is vertical kernel G with directionyConvolution algorithm is carried out with image, then in the gray scale G of the every bit of image
Size can be approximated to be | I3*Gx|+|I3*Gy|；Then Hough circle transformation is used, image space is transformed into parameter space, is selected
Suitable radius step-length step_r, angle step step_angle, minimum radius of circle r_min, maximum radius of circle r_max are taken, it is right
Point on edge carries out statistics calculating, determines circle that may be present, obtains initial center of circle collection C0(x,y)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811360167.3A CN109523565B (en) | 2018-11-15 | 2018-11-15 | Method and system for positioning center of diffraction-free moire fringe |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811360167.3A CN109523565B (en) | 2018-11-15 | 2018-11-15 | Method and system for positioning center of diffraction-free moire fringe |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109523565A true CN109523565A (en) | 2019-03-26 |
| CN109523565B CN109523565B (en) | 2020-07-10 |
Family
ID=65777815
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811360167.3A Active CN109523565B (en) | 2018-11-15 | 2018-11-15 | Method and system for positioning center of diffraction-free moire fringe |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109523565B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111811784A (en) * | 2020-07-02 | 2020-10-23 | 中航华东光电（上海）有限公司 | Laser spot center coordinate determination method, device and equipment |
| CN112116581A (en) * | 2020-09-23 | 2020-12-22 | 中国科学院物理研究所 | Method and apparatus for acquiring atomic position in atomic imaging |
| CN112231939A (en) * | 2020-01-03 | 2021-01-15 | 郑州轻工业大学 | Ant colony sequencing positioning method for circular layout in cable processing |
| CN112344905A (en) * | 2020-10-28 | 2021-02-09 | 孝感华中精密仪器有限公司 | Roadbed settlement detection device and method for prolonging battery replacement period of roadbed settlement detection device |
| CN112710632A (en) * | 2021-03-23 | 2021-04-27 | 四川京炜交通工程技术有限公司 | Method and system for detecting high and low refractive indexes of glass beads |
| CN112782203A (en) * | 2021-02-22 | 2021-05-11 | 长江存储科技有限责任公司 | Judging method of crystalline phase structure and crystalline phase calibration template |
| CN113100942A (en) * | 2021-04-12 | 2021-07-13 | 中国科学院苏州生物医学工程技术研究所 | Laser point identification method and SS-OCT operation navigation system using same |
| CN113532320A (en) * | 2021-07-20 | 2021-10-22 | 武汉华工激光工程有限责任公司 | Image-based light spot diffraction ring analysis method, storage medium and chip |
| CN113793282A (en) * | 2021-09-17 | 2021-12-14 | 中国科学院长春光学精密机械与物理研究所 | Method for removing moire fringe blur of space camera transfer function test image |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101216895A (en) * | 2007-12-26 | 2008-07-09 | 北京航空航天大学 | An automatic extracting method for ellipse image features in complex background images |
| CN101359400A (en) * | 2008-09-17 | 2009-02-04 | 湖南大学 | Process for positioning spatial position of pipe mouth based on vision |
| JP2017093496A (en) * | 2015-11-18 | 2017-06-01 | キヤノン株式会社 | Imaging device |
| CN106919792A (en) * | 2017-02-24 | 2017-07-04 | 天津大学 | Vortex center automatic identifying method based on high accuracy numerical value Wind Data |
| CN108230388A (en) * | 2018-02-06 | 2018-06-29 | 广西艾盛创制科技有限公司 | A kind of recognition positioning method of white body weld point image |
| CN108600414A (en) * | 2018-05-09 | 2018-09-28 | 中国平安人寿保险股份有限公司 | Construction method, device, storage medium and the terminal of device-fingerprint |
-
2018
- 2018-11-15 CN CN201811360167.3A patent/CN109523565B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101216895A (en) * | 2007-12-26 | 2008-07-09 | 北京航空航天大学 | An automatic extracting method for ellipse image features in complex background images |
| CN101359400A (en) * | 2008-09-17 | 2009-02-04 | 湖南大学 | Process for positioning spatial position of pipe mouth based on vision |
| JP2017093496A (en) * | 2015-11-18 | 2017-06-01 | キヤノン株式会社 | Imaging device |
| CN106919792A (en) * | 2017-02-24 | 2017-07-04 | 天津大学 | Vortex center automatic identifying method based on high accuracy numerical value Wind Data |
| CN108230388A (en) * | 2018-02-06 | 2018-06-29 | 广西艾盛创制科技有限公司 | A kind of recognition positioning method of white body weld point image |
| CN108600414A (en) * | 2018-05-09 | 2018-09-28 | 中国平安人寿保险股份有限公司 | Construction method, device, storage medium and the terminal of device-fingerprint |
Non-Patent Citations (3)
| Title |
|---|
| MA HONG .ETC: "Research on clustering algoriths of data streams", 《2010 2ND IEEE INTERNATIONAL CONFERENCE ON INFORMATION MANAGEMENT AND ENGINEERING》 * |
| 李宁: "莫尔圆的特征参数提取方法", 《光电技术应用》 * |
| 黄孝川: "定弧长多重循环无衍射光莫尔条纹计数方法", 《湖北工业大学学报》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112231939A (en) * | 2020-01-03 | 2021-01-15 | 郑州轻工业大学 | Ant colony sequencing positioning method for circular layout in cable processing |
| CN112231939B (en) * | 2020-01-03 | 2024-02-02 | 郑州轻工业大学 | Ant colony sequencing and positioning method for circular layout in cable processing |
| CN111811784A (en) * | 2020-07-02 | 2020-10-23 | 中航华东光电（上海）有限公司 | Laser spot center coordinate determination method, device and equipment |
| CN112116581B (en) * | 2020-09-23 | 2023-09-08 | 中国科学院物理研究所 | Method and device for acquiring atomic position in atomic imaging |
| CN112116581A (en) * | 2020-09-23 | 2020-12-22 | 中国科学院物理研究所 | Method and apparatus for acquiring atomic position in atomic imaging |
| CN112344905A (en) * | 2020-10-28 | 2021-02-09 | 孝感华中精密仪器有限公司 | Roadbed settlement detection device and method for prolonging battery replacement period of roadbed settlement detection device |
| CN112782203A (en) * | 2021-02-22 | 2021-05-11 | 长江存储科技有限责任公司 | Judging method of crystalline phase structure and crystalline phase calibration template |
| CN112782203B (en) * | 2021-02-22 | 2024-02-20 | 长江存储科技有限责任公司 | Crystal phase structure judging method and crystal phase calibration template |
| CN112710632A (en) * | 2021-03-23 | 2021-04-27 | 四川京炜交通工程技术有限公司 | Method and system for detecting high and low refractive indexes of glass beads |
| CN113100942A (en) * | 2021-04-12 | 2021-07-13 | 中国科学院苏州生物医学工程技术研究所 | Laser point identification method and SS-OCT operation navigation system using same |
| CN113532320A (en) * | 2021-07-20 | 2021-10-22 | 武汉华工激光工程有限责任公司 | Image-based light spot diffraction ring analysis method, storage medium and chip |
| CN113793282A (en) * | 2021-09-17 | 2021-12-14 | 中国科学院长春光学精密机械与物理研究所 | Method for removing moire fringe blur of space camera transfer function test image |
| CN113793282B (en) * | 2021-09-17 | 2023-10-24 | 中国科学院长春光学精密机械与物理研究所 | Moire fringe blurring removal method for space camera transfer function test image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109523565B (en) | 2020-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109523565A (en) | A kind of diffraction light-free Moire fringe center positioning method and system | |
| CN109060836B (en) | Machine vision-based high-pressure oil pipe joint external thread detection method | |
| CN109325935B (en) | Power transmission line detection method based on unmanned aerial vehicle image | |
| CN109859226B (en) | Detection method of checkerboard corner sub-pixels for graph segmentation | |
| Hussain et al. | A comparative analysis of edge detection techniques used in flame image processing | |
| CN109785245A (en) | A kind of light spot image dressing method | |
| CN109685045A (en) | A kind of Moving Targets Based on Video Streams tracking and system | |
| Sahoo et al. | Design and simulation of various edge detection techniques using Matlab Simulink | |
| CN109434251A (en) | A kind of weld image tracking based on particle filter | |
| CN103295227B (en) | A kind of circle detection method based on gradient direction segmentation | |
| CN106709500A (en) | Image feature matching method | |
| CN106951860A (en) | A kind of three-dimensional data intelligent identification Method based on a cloud | |
| CN104282026A (en) | Distribution uniformity assessment method based on watershed algorithm and minimum spanning tree | |
| Masuda | Leaf area estimation by semantic segmentation of point cloud of tomato plants | |
| CN112669280B (en) | Unmanned aerial vehicle inclination aerial photography right-angle image control point target detection method based on LSD algorithm | |
| KR20110094957A (en) | Apparatus and method for object segmentation from range image | |
| CN105678737A (en) | Digital image corner point detection method based on Radon transform | |
| CN109636785A (en) | A kind of visual processing method identifying particles of silicon carbide | |
| CN107330881A (en) | A kind of method for recognizing wheel hub valve port | |
| Fatima et al. | An extensive survey on edge detection techniques | |
| CN113554688B (en) | O-shaped sealing ring size measurement method based on monocular vision | |
| CN109544513A (en) | A kind of steel pipe end surface defect extraction knowledge method for distinguishing | |
| Wang et al. | Global contrast of superpixels based salient region detection | |
| Janoriya et al. | Critical review on edge detection techniques in spatial domain on low illumination images | |
| Hou | Analysis of visual navigation extraction algorithm of farm robot based on dark primary colour. |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230815 Address after: 430060 Block B, Building 10, Huijin Center, Financial Port Fourth Road, Donghu New Technology Development Zone, Wuhan City, Hubei Province Patentee after: Wuhan Puchuang Data Technology Co.,Ltd. Address before: 430068 1, Lijia 1 village, Nanhu, Wuchang District, Wuhan, Hubei Patentee before: HUBEI University OF TECHNOLOGY |