CN104899890A - Detection method based on polarized light imaging technology - Google Patents
Detection method based on polarized light imaging technology Download PDFInfo
- Publication number
- CN104899890A CN104899890A CN201510351426.6A CN201510351426A CN104899890A CN 104899890 A CN104899890 A CN 104899890A CN 201510351426 A CN201510351426 A CN 201510351426A CN 104899890 A CN104899890 A CN 104899890A
- Authority
- CN
- China
- Prior art keywords
- parameter
- image
- frequency coefficient
- degree
- high frequency
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 60
- 238000001514 detection method Methods 0.000 title claims abstract description 56
- 238000005516 engineering process Methods 0.000 title abstract description 3
- 230000010287 polarization Effects 0.000 claims abstract description 104
- 238000000354 decomposition reaction Methods 0.000 claims description 13
- 238000000034 method Methods 0.000 claims description 9
- 230000015556 catabolic process Effects 0.000 claims description 6
- 230000007547 defect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 3
- 230000002708 enhancing effect Effects 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 6
- 238000012935 Averaging Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000001303 quality assessment method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 208000031968 Cadaver Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 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/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/0008—Industrial image inspection checking presence/absence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
- G01N21/88—Investigating the presence of flaws or contamination
- G01N21/8851—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
- G01N2021/8887—Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges based on image processing techniques
-
- 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
-
- 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/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
Abstract
The invention relates to the technical field of image processing, and especially relates to a detection method and a detection device based on a polarized light imaging technology. According to the invention, polarized light imaging and intensity imaging principles are combined to fuse multi-source images, an enhanced fusion image is generated by utilizing complementation information in polarization image information and intensity image information, the details of the fusion image is more outstanding than that of any of polarization images and intensity images, and a new enhanced image deriving from the multi-source images is formed. The contrast of objective discrimination characteristics in the fusion image is far higher than the contrast of the objective discrimination characteristics in the intensity images, so that describing of relative scenes which is unable to obtained by a single image is realized, the limitation of the single image is overcome, the detection efficiency of objects to be detected is improved, and the accuracy of the detection efficiency is enhanced.
Description
Technical field
The present invention relates to a kind of technical field of image processing, particularly relate to a kind of detection method based on polarized light imaging technique.
Background technology
Dynamic Non-Destruction Measurement and non-damaged data, utilize the characteristics such as the sound of thing to be detected, light, magnetic and electricity, do not destroying under the prerequisites such as the original state of thing to be detected, chemical property, for obtaining the inspection method that physics, the chemical corps intelligences such as content, character or the composition relevant with article matter to be detected adopt, the existing lossless detection method being applied to industrial detection mainly contains mechanical type contact detection method and photo-electric non-contact detection method.Wherein mechanical type contact detection method is measured in the slippage of measured surface by the contact pilotage of instrument, this kind of method needs contact pilotage to carry out comprehensive engagement on the surface of thing to be detected, thus detection speed is slow, and easily scratch thing to be detected surface, cause secondary cut, photo-electric non-contact detection method is now widely used method, by analyzing the two dimensional image feature of testee, can realize quick, non-contact detecting.Existing non-contact detection method mainly adopts light intensity image checking and polarized light image checking, light intensity image checking is by light source irradiation thing to be detected, obtain the radiation intensity information on thing surface to be detected, and according to the flaw of thing to be detected in this radiation intensity information combining image Processing Algorithm detected image.Because thing surface blemish feature to be detected (as cut, surface is uneven, and smoothness is low) is usually more hidden, not easily present in intensity imaging system, this just for later image process brings larger pressure, not only can increase algorithm complexity, can reduce the detection efficiency of system equally.Polarized light image checking is that the target reflecting light by obtaining different polarization direction obtains target polarization information, target polarization parameter information is obtained by carrying out parsing to target polarization information, as degree of polarization, polarization angle etc., carry out Inversion Calculation according to target polarization parameter information and obtain object reconstruction figure, object reconstruction image comprises geometric configuration, surfaceness, texture, the polarization information of the physicochemical properties such as conductance, but there is certain difference in the surface of object reconstruction image and background on some polarization parameter and original image, compare intensity image, polarization parameter image readability is not good, need later stage decipher.
Summary of the invention
Thing surface blemish to be detected is detected in order to better utilize polarization characteristic, improve detection efficiency, the invention provides a kind of detection method based on polarized light imaging technique, by intensity imaging combine with technique polarized light imaging technique, intensity image basis utilizes carry out of polarization characteristic to thing surface blemish to be detected strengthen, improve detection efficiency and detection accuracy.
Based on a detection method for polarized light imaging technique, be applied to thing surface blemish to be detected and detect, wherein, comprise the steps,
Step S1, obtain the light intensity parameter on described thing surface to be detected in three predetermined angulars by a light intensity imaging device;
Step S2, obtain described thing surface image to be detected in three described predetermined angulars by a polarized light image-forming detecting system, forms the described surface image output that three light intensities are different;
Step S3, a computing unit described surface image different according to three light intensities obtains Stokes parameter and degree of polarization parameter P in conjunction with described light intensity parameter;
Step S4, described surface image is divided into several little images, calculates the average energy value of each described little image according to described Stokes parameter and degree of polarization parameter P,
Step S5, to choose the highest described little image of the average energy value be image to be fused;
Step S6, wavelet decomposition is carried out to the light intensity parameter I of described image to be fused and degree of polarization parameter P, obtain the low frequency coefficient of the light intensity parameter I of described image to be fused, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P respectively;
Step S7, the low frequency coefficient according to described light intensity parameter I, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P obtain low frequency coefficient to be fused, high frequency coefficient to be fused,
Step S8, carry out Image Reconstruction according to described high frequency coefficient to be fused and described low frequency coefficient to be fused, obtain fused images;
Step S9, Defect Detection is carried out to fused images.
The above-mentioned detection method based on polarized light imaging technique, wherein, three predetermined angulars are defined as α 1, α 2, α 3 respectively; Three light intensity parameters form I respectively
0(α 1), I
0(α 2), I
0(α 3), I
0the light intensity parameter on the thing surface described to be detected that (α 1) obtains for α 1 polarization direction, I
0the light intensity parameter on the thing surface described to be detected that (α 2) obtains for α 2 polarization direction, I
0the light intensity parameter on the thing surface described to be detected that (α 3) obtains for α 3 polarization direction.
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step S3, wherein, the computing formula of described Stokes parameter is:
Wherein, I is light intensity parameter, and Q is the first polarization parameter, and U is the second polarization parameter, and V is circular polarization parameter, and because in natural light, V is less, value is zero;
The computing formula of described degree of polarization parameter P is:
Wherein, P is degree of polarization parameter.
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step S4,
The computing formula of wherein said the average energy value is:
Wherein: E is the average energy value, M is the length of each described little image, N is the width of each described little image, x be each described little image in the picture breakdown coordinate points in horizontal ordinate direction, y is that each described little image is in the picture breakdown coordinate points in ordinate direction.
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step S6, also specifically comprises,
Step S61, carry out wavelet decomposition to the light intensity parameter I of described image to be fused and degree of polarization parameter P, Decomposition order is 1;
Step S62, obtain the low frequency coefficient (I of the light intensity parameter I of described image to be fused
lL), the high frequency coefficient (I of light intensity parameter I
lH, I
hL, I
hH),
Step S63, obtain the low frequency coefficient (P of the described degree of polarization parameter P of described image to be fused
lL), the high frequency coefficient (P of described degree of polarization parameter P
lH, P
hL, P
hH).
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step S7, specifically comprises the steps,
Step S71, low frequency coefficient (I according to described light intensity parameter I
lL) with the low frequency coefficient (P of described degree of polarization parameter P
lL) obtain low frequency coefficient to be fused, described low frequency coefficient F to be fused
lLfusion formula be:
F
LL=log(P
LL)·I
LL n
Wherein, n is I
lLindex, span is 0 ~ 1;
Step S72, high frequency coefficient (I according to described light intensity parameter I
lH, I
hL, I
hH) with the high frequency coefficient (P of described degree of polarization parameter P
lH, P
hL, P
hH) obtain high frequency coefficient F to be fused
h,
Wherein specifically comprise the steps,
High frequency coefficient (the I of step S721, definition light intensity parameter I
lH, I
hL, I
hH) be I
m; High frequency coefficient (the P of definition degree of polarization parameter P
lH, P
hL, P
hH) be P
m;
Step S722, set up light intensity parameter high frequency coefficient I
mwith described degree of polarization parameter P
mcorrelation coefficient M
i,P,
M
i,Pcomputing formula be:
Wherein, E
hPfor the high frequency coefficient energy value of degree of polarization parameter P, P
mfor the high frequency coefficient of degree of polarization parameter P; E
hIfor the high frequency coefficient energy value of intensive parameter I, I
mfor the high frequency coefficient of light intensity parameter I;
Step S722, a default relative degree threshold value δ,
Step S723, calculate high frequency coefficient F to be fused according to described relative degree threshold value δ
h
Work as M
i,Pduring > δ,
Work as M
i,Pduring < δ, F
h=ε
hIi
h+ ε
hPp
h,
Wherein ε
hPfor the weighting coefficient of degree of polarization parameter P, ε
hIfor the weighting coefficient of light intensity parameter I.
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step 8, the reconstructing method of its image to be fused is:
Wherein F
lLfor view data to be fused,
h
r, H
c, G
r, G
cassociate matrix, F
lHbe the 1st layer of vertical direction high frequency imaging, F
hLbe the 1st layer of horizontal direction high frequency imaging, F
hHbe the 1st layer to angular direction high frequency imaging.
The above-mentioned detection method based on polarized light imaging technique, wherein, it is 0 ° that described predetermined angular is defined as α 1 respectively, and α 2 is 60 °, and α 3 is 120 °.
The above-mentioned detection method based on polarized light imaging technique, wherein, in described step 6, described wavelet decomposition is Mallat Wavelet Transformation Algorithm.
Compared with prior art, advantage of the present invention is:
Polarized light imaging and intensity imaging principle combine by the present invention, multi-source image (image that polarized light imaging obtains and the image that intensity imaging obtains) is merged, complementary information in polarization image information and intensity image information is utilized to generate the fused images of a width enhancing, the details of fused images is all more outstanding than a width any in polarization image and intensity image, what form a width enhancing has drawn the new images of multi-source image, in fused images, the contrast of target distinguishing characteristics is far above the contrast of target distinguishing characteristics in intensity image, thus obtain the description of the associated scenario that single image cannot obtain, overcome the limitation existing for single image, improve the detection efficiency of thing to be detected, enhance the accuracy of detection efficiency.
Accompanying drawing explanation
Fig. 1 is a kind of detection method schematic flow sheet based on polarized light imaging technique.
Embodiment
Below in conjunction with the drawings and specific embodiments, the invention will be further described, but not as limiting to the invention.
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, the every other embodiment that those of ordinary skill in the art obtain under the prerequisite of not making creative work, all belongs to the scope of protection of the invention.
It should be noted that, when not conflicting, the embodiment in the present invention and the feature in embodiment can combine mutually.
Below in conjunction with the drawings and specific embodiments, the invention will be further described, but not as limiting to the invention.
As shown in Figure 1, a kind of detection method based on polarized light imaging technique, is applied to thing surface blemish to be detected and detects, wherein, comprise the steps,
Step S1, obtain the light intensity parameter on thing surface to be detected in three predetermined angulars by a light intensity imaging device, predetermined angular is defined as α 1, α 2, α 3 respectively; Light intensity parameter comprises I
0(α 1), I
0(α 2), I
0(α 3), wherein, I
0the light intensity parameter on the thing surface to be detected that (α 1) obtains for α 1 polarization direction, I
0the light intensity parameter on the thing surface to be detected that (α 2) obtains for α 2 polarization direction, I
0the light intensity parameter on the thing surface to be detected that (α 3) obtains for α 3 polarization direction;
Step S2, obtain thing surface image to be detected in three predetermined angulars by a polarized light image-forming detecting system, the surface image forming three light intensities different exports;
Step S3, the computing unit surface image different according to three light intensities obtains Stokes parameter and degree of polarization parameter P in conjunction with described light intensity parameter;
Step S4, surface image is divided into several little images, calculates the average energy value of each little image according to Stokes parameter and degree of polarization parameter P,
Step S5, to choose the highest little image of the average energy value be image to be fused;
Step S6, the light intensity parameter I treating fused images and degree of polarization parameter P carry out wavelet decomposition, obtain the low frequency coefficient of the light intensity parameter I of image to be fused, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P respectively;
Step S7, the low frequency coefficient according to light intensity parameter I, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P obtain low frequency coefficient to be fused, high frequency coefficient to be fused,
Step S8, carry out Image Reconstruction according to high frequency coefficient to be fused and low frequency coefficient to be fused, obtain fused images;
Step S9, Defect Detection is carried out to fused images.
Polarized light imaging and intensity imaging principle combine by the present invention, multi-source image (image that polarized light imaging obtains and the image that intensity imaging obtains) is merged, complementary information in polarization image information and intensity image information is utilized to generate the fused images of a width enhancing, the details of fused images is all more outstanding than a width any in polarization image and intensity image, what form a width enhancing has drawn the new images of multi-source image, in fused images, the contrast of target distinguishing characteristics is far above the contrast of target distinguishing characteristics in intensity image, thus obtain the description of the associated scenario that single image cannot obtain, overcome the limitation existing for single image, improve the detection efficiency of thing to be detected, enhance the accuracy of detection efficiency.
The above-mentioned detection method based on polarized light imaging technique, wherein, in step S3, wherein, the computing formula of Stokes parameter is:
Wherein, I is light intensity parameter, and Q is the first polarization parameter, and U is the second polarization parameter, and V is circular polarization parameter, and because in natural light, V is less, value is zero;
The computing formula of degree of polarization parameter P is:
Wherein, P is degree of polarization parameter.
The above-mentioned detection method based on polarized light imaging technique, wherein, in step S4,
Wherein the computing formula of the average energy value is:
Wherein: E is the average energy value, M is the length of each little image, and N is the width of each little image, x be each little image in the picture breakdown coordinate points in horizontal ordinate direction, y is that each little image is in the picture breakdown coordinate points in ordinate direction.
The above-mentioned detection method based on polarized light imaging technique, wherein, in step S6, also specifically comprises,
Step S61, the light intensity parameter I treating fused images and degree of polarization parameter P carry out wavelet decomposition, and Decomposition order is 1;
Step S62, obtain the low frequency coefficient (I of the light intensity parameter I of image to be fused
lL),
High frequency coefficient (the I of light intensity parameter I
lH, I
hL, I
hH),
Step S63, obtain the low frequency coefficient (P of the degree of polarization parameter P of image to be fused
lL), the high frequency coefficient (P of degree of polarization parameter P
lH, P
hL, P
hH),
The above-mentioned detection method based on polarized light imaging technique, wherein, in step S7, specifically comprises the steps,
Step S71, low frequency coefficient (I according to light intensity parameter I
lL) with the low frequency coefficient (P of degree of polarization parameter P
lL) obtain low frequency coefficient to be fused, low frequency coefficient F to be fused
lLfusion formula be:
F
LL=log(P
LL)·I
LL n
Wherein, n is I
lLindex, span is 0 ~ 1;
Low frequency coefficient F to be fused
lLreflect the approximate of original image and evenness, concentrated the most information of original image, adopt modulation fusion method to obtain to low frequency coefficient.
Step S72, high frequency coefficient (I according to light intensity parameter I
lH, I
hL, I
hH) with the high frequency coefficient (P of degree of polarization parameter P
lH, P
hL, P
hH) obtain high frequency coefficient F to be fused
h,
Wherein specifically comprise
High frequency coefficient (the I of step S721, definition light intensity parameter I
lH, I
hL, I
hH) be I
m; High frequency coefficient (the P of definition degree of polarization parameter P
lH, P
hL, P
hH) be P
m;
Step S722, high frequency coefficient F to be fused reflect the catastrophe characteristics of original image, to HFS respectively along horizontal direction, vertical direction with adopt coefficient weighted mean fusion rule to merge to angular direction.The high frequency coefficient of image after multi-resolution decomposition contains the detailed information of such as edge, region contour etc. in image.When fusion treatment, consider the correlativity between adjacent coefficient, set up light intensity parameter high frequency coefficient I
mwith described degree of polarization parameter P
mcorrelation coefficient M
i,P,
M (
i,
p) computing formula be:
Wherein, E
hPfor the high frequency coefficient energy value of degree of polarization parameter P, P
mfor the high frequency coefficient of degree of polarization parameter P;
E
hIfor the high frequency coefficient energy value of intensive parameter I, I
mfor the high frequency coefficient of light intensity parameter I;
Step S722, default a relative degree threshold value δ, relative degree threshold value δ according to actual selection, can not do concrete restriction herein.
Step S723, calculate high frequency coefficient F to be fused according to relative degree threshold value δ
h
Work as M
i,Pduring > δ,
Work as M
i,Pduring < δ, F
h=ε
hIi
h+ ε
hPp
h,
ε
HP=1-ε
HI
Wherein ε
hPfor the weighting coefficient of degree of polarization parameter P, ε
hIfor the weighting coefficient of light intensity parameter I.
The above-mentioned detection method based on polarized light imaging technique, wherein, in step 8, the reconstructing method of its image to be fused is:
Wherein F
lLfor view data to be fused,
h
r, H
c, G
r, G
cassociate matrix, F
lHbe the 1st layer of vertical direction high frequency imaging, F
hLbe the 1st layer of horizontal direction high frequency imaging, F
hHbe the 1st layer to the relational expression between the high frequency imaging of angular direction.
The above-mentioned detection method based on polarized light imaging technique, wherein, it is 0 ° that predetermined angular is defined as α 1 respectively, and α 2 is 60 °, and α 3 is 120.
The above-mentioned detection method based on polarized light imaging technique, wherein, in step 6, wavelet decomposition is Mallat Wavelet Transformation Algorithm.
Enumerate a specific embodiment to explain:
CD is as thing to be detected, and the cut that optical disc surface has two places less, the degree of depth is about 0.5mm, and length is about 12mm.For obtaining stronger polarization information, select to test under fixing special light source incidence angle conditions.Gathered the original polarization frame of optical disc surface flaw respectively in three angles (0 °, 60 °, 120 °) by polarized light image-forming detecting system, and the image after gathering is handled as follows: obtain strength parameter I and degree of polarization parameter P, treat fused images and use the above-mentioned detection method based on polarized light imaging technique to obtain fused images F.
Effect is strengthened in order to verify to merge, using strength parameter I and degree of polarization parameter P as with reference to image, to the image of high light imaging, degree of polarization image and adopt the application to obtain fused images respectively computed image quality assessment parameter (comprising image information entropy, image averaging gradient, sharpness, contrast and Gauss three rank details average statistics amount) carry out image quality evaluation.
The computing formula of image information entropy is
Wherein, pi is the probability that the gray level of i occurs, i is gray level, and H is image information entropy.
Image information entropy weighs image information to enrich an important indicator of degree, the size of the entropy of fused images represent the average information that image comprises number.By the details expressive ability that relatively can draw image to image information entropy.The size of entropy reflects the number of the quantity of information that image carries, and entropy is larger, illustrates that the quantity of information of carrying is larger.
The computing formula of image averaging gradient is
Wherein,
image X-direction derivative,
it is image Y-direction derivative.
Average gradient is that the index of responsive reflection image to minor detail contrast and texture variations feature representation ability is general, and g is larger, and image level is more, and fused image is more clear, merges to reach to put forward high-resolution object.
Sharpness computation formula is
Wherein, I here
xand I
ybe respectively the difference of image I in x and y direction, i is X-direction coordinate figure in image, and j is Y-direction coordinate figure in image.The sharpness of image represents with average gradient, and sharpness shows that more greatly the details such as image texture are more obvious, is more easy to observe.
The computing formula of the contrast gray-scale value contrast of object and background (in the image) is:
Wherein W
t(x, y), W
b(x, y) is respectively target and background image intensity value, and R is selected zone, and E is mathematical expectation, and C is the gray-scale value contrast of object and background in image.
The concrete formula of Gauss three rank details average statistics amount is:
Wherein, c
3ffor Gauss three rank details average statistics amount, m
1f(i, j) is the local mean value in the picture in (i, j) the little field of putting, m
2f(i, j) is the second moment of image f (i, j), m
3f(i, j) is the third moment of image f (i, j), and K is the departure of coordinate i, and Q is image template size, and i is X-direction coordinate figure in image, and j is Y-direction coordinate figure in image, and l is the departure of coordinate j.
Gauss three rank details average statistics metering method does not need original image for referencial use, insensitive to white Gaussian noise, can reflect the signal to noise ratio (S/N ratio) of image preferably, wherein defines the larger picture quality of numerical value higher.
Table 1 is image, the degree of polarization image of high light imaging and adopts the fused images of the application's acquisition to distinguish the picture appraisal result of computed image quality assessment parameter.
Table 1 picture appraisal result
Can find from table 1 picture appraisal result:
(1) information entropy comparing strength parameter image, degree of polarization parametric image and fused images can be found out, fused images is compared degree of polarization parametric image information entropy and improved a lot.Compared with degree of polarization parametric image, the quantity of information of fused images obtains enhancing;
(2) average gradient of the fused images obtained through blending algorithm is greatly improved, and improves overall sharpness, and the detail section reflecting image obtains enhancing;
(3) from the result of the contrast of optical disc surface flaw and intact optical disc surface, in fused images, target contrast is far above strength parameter image;
(4) from visual effect, fused images remains the scene information in strength parameter image, highlights the detailed information of optical disc surface flaw in degree of polarization parametric image simultaneously;
(5) visible according to Gauss three rank details average statistics amount result of calculation, the picture quality after fusion is all better than the strength parameter image before merging and degree of polarization parametric image;
(6) detection and indentification that enhancing result is easier to optical disc surface flaw is merged.
These are only preferred embodiment of the present invention; not thereby embodiments of the present invention and protection domain is limited; to those skilled in the art; should recognize and all should be included in the scheme that equivalent replacement done by all utilizations instructions of the present invention and diagramatic content and apparent change obtain in protection scope of the present invention.
Claims (9)
1., based on a detection method for polarized light imaging technique, be applied to thing surface blemish to be detected and detect, it is characterized in that, comprise the steps,
Step S1, obtain the light intensity parameter on described thing surface to be detected in three predetermined angulars by a light intensity imaging device;
Step S2, obtain described thing surface image to be detected in three described predetermined angulars by a polarized light image-forming detecting system, forms the described surface image output that three light intensities are different;
Step S3, a computing unit described surface image different according to three light intensities obtains Stokes parameter and degree of polarization parameter P in conjunction with described light intensity parameter;
Step S4, described surface image is divided into several little images, calculates the average energy value of each described little image according to described Stokes parameter and degree of polarization parameter P,
Step S5, to choose the highest described little image of the average energy value be image to be fused;
Step S6, wavelet decomposition is carried out to the light intensity parameter I of described image to be fused and degree of polarization parameter P, obtain the low frequency coefficient of the light intensity parameter I of described image to be fused, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P respectively;
Step S7, the low frequency coefficient according to described light intensity parameter I, the high frequency coefficient of light intensity parameter I, the low frequency coefficient of degree of polarization parameter P, the high frequency coefficient of degree of polarization parameter P obtain low frequency coefficient to be fused, high frequency coefficient to be fused,
Step S8, carry out Image Reconstruction according to described high frequency coefficient to be fused and described low frequency coefficient to be fused, obtain fused images;
Step S9, Defect Detection is carried out to fused images.
2. the detection method based on polarized light imaging technique according to claim 1, is characterized in that, three predetermined angulars are defined as α 1, α 2, α 3 respectively; Three light intensity parameters form I respectively
0(α 1), I
0(α 2), I
0(α 3), wherein I
0the light intensity parameter on the thing surface described to be detected that (α 1) obtains for α 1 polarization direction, I
0the light intensity parameter on the thing surface described to be detected that (α 2) obtains for α 2 polarization direction, I
0the light intensity parameter on the thing surface described to be detected that (α 3) obtains for α 3 polarization direction.
3. the detection method based on polarized light imaging technique according to claim 2, is characterized in that, in described step S3, wherein, the computing formula of described Stokes parameter is:
Wherein, I is light intensity parameter, and Q is the first polarization parameter, and U is the second polarization parameter, and V is circular polarization parameter, and because in natural light, V is less, value is zero;
The computing formula of described degree of polarization parameter P is:
Wherein, P is degree of polarization parameter.
4. the detection method based on polarized light imaging technique according to claim 1, is characterized in that, in described step S4,
The computing formula of wherein said the average energy value is:
Wherein: E is the average energy value, M is the length of each described little image, N is the width of each described little image, x be each described little image in the picture breakdown coordinate points in horizontal ordinate direction, y is that each described little image is in the picture breakdown coordinate points in ordinate direction.
5. the detection method based on polarized light imaging technique according to claim 1, is characterized in that, in described step S6, also specifically comprises,
Step S61, carry out wavelet decomposition to the light intensity parameter I of described image to be fused and degree of polarization parameter P, Decomposition order is 1;
Step S62, obtain the low frequency coefficient (I of the light intensity parameter I of described image to be fused
lL), the high frequency coefficient (I of light intensity parameter I
lH, I
hL, I
hH),
Step S63, obtain the low frequency coefficient (P of the degree of polarization parameter P of described image to be fused
lL), the high frequency coefficient (P of degree of polarization parameter P
lH, P
hL, P
hH).
6. the detection method based on polarized light imaging technique according to claim 5, is characterized in that, in described step S7, specifically comprises the steps,
Step S71, low frequency coefficient (I according to described light intensity parameter I
lL) with the low frequency coefficient (P of described degree of polarization parameter P
lL) obtain low frequency coefficient to be fused, described low frequency coefficient F to be fused
lLfusion formula be:
F
LL=log(P
LL)·I
LL n
Wherein, n is I
lLindex, span is 0 ~ 1;
Step S72, high frequency coefficient (I according to described light intensity parameter I
lH, I
hL, I
hH) with the high frequency coefficient (P of described degree of polarization parameter P
lH, P
hL, P
hH) obtain high frequency coefficient F to be fused
h,
Wherein specifically comprise
High frequency coefficient (the I of step S721, definition light intensity parameter I
lH, I
hL, I
hH) be I
m; High frequency coefficient (the P of definition degree of polarization parameter P
lH, P
hL, P
hH) be P
m;
Step S722, set up light intensity parameter high frequency coefficient I
mwith described degree of polarization parameter P
mcorrelation coefficient M
i,P,
M
i,Pcomputing formula be:
Wherein, E
hPfor the high frequency coefficient energy value of degree of polarization parameter P, P
mfor the high frequency coefficient of degree of polarization parameter P;
E
hIfor the high frequency coefficient energy value of intensive parameter I, I
mfor the high frequency coefficient of light intensity parameter I;
Step S722, a default relative degree threshold value δ,
Step S723, calculate high frequency coefficient F to be fused according to described relative degree threshold value δ
h
Work as M
i,Pduring > δ,
Work as M
i,Pduring < δ, F
h=ε
hIi
h+ ε
hPp
h,
Wherein ε
hPfor the weighting coefficient of degree of polarization parameter P, ε
hIfor the weighting coefficient of light intensity parameter I.
7. the detection method based on polarized light imaging technique according to claim 6, is characterized in that, in described step 8, the reconstructing method of its image to be fused is:
Wherein F
lLfor view data to be fused,
h
r, H
c, G
r, G
cassociate matrix, F
lHbe the 1st layer of vertical direction high frequency imaging, F
hLbe the 1st layer of horizontal direction high frequency imaging, F
hHbe the 1st layer to angular direction high frequency imaging.
8. the detection method based on polarized light imaging technique according to claim 1, is characterized in that, it is 0 ° that described predetermined angular is defined as α 1 respectively, and α 2 is 60 °, and α 3 is 120 °.
9. the detection method based on polarized light imaging technique according to claim 1, is characterized in that, in described step 6, described wavelet decomposition is Mallat Wavelet Transformation Algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510351426.6A CN104899890A (en) | 2015-06-23 | 2015-06-23 | Detection method based on polarized light imaging technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510351426.6A CN104899890A (en) | 2015-06-23 | 2015-06-23 | Detection method based on polarized light imaging technology |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104899890A true CN104899890A (en) | 2015-09-09 |
Family
ID=54032535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510351426.6A Pending CN104899890A (en) | 2015-06-23 | 2015-06-23 | Detection method based on polarized light imaging technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104899890A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110044931A (en) * | 2019-04-23 | 2019-07-23 | 华中科技大学 | A kind of detection device on bend glass surface and internal flaw |
CN110412035A (en) * | 2019-06-13 | 2019-11-05 | 合刃科技(武汉)有限公司 | A kind of high reflecting material surface inspecting method and system |
CN111328370A (en) * | 2019-04-23 | 2020-06-23 | 合刃科技(深圳)有限公司 | Surface defect detection system and method |
CN114324369A (en) * | 2022-03-11 | 2022-04-12 | 北京新研创能科技有限公司 | System and method for detecting scratches on surface of bipolar plate |
CN114460104A (en) * | 2022-04-13 | 2022-05-10 | 清远市通用皮具配件有限公司 | Zipper tooth arrangement injection molding quality detection system and method based on machine vision |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530853A (en) * | 2013-10-17 | 2014-01-22 | 中北大学 | Infrared intensity image and infrared polarization image enhancement and fusion method |
CN104574296A (en) * | 2014-12-24 | 2015-04-29 | 长春理工大学 | Multi-wavelet fusion image processing method for polarization haze removal |
-
2015
- 2015-06-23 CN CN201510351426.6A patent/CN104899890A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103530853A (en) * | 2013-10-17 | 2014-01-22 | 中北大学 | Infrared intensity image and infrared polarization image enhancement and fusion method |
CN104574296A (en) * | 2014-12-24 | 2015-04-29 | 长春理工大学 | Multi-wavelet fusion image processing method for polarization haze removal |
Non-Patent Citations (3)
Title |
---|
张晶晶 等: "基于小波的偏振图像融合算法及性能评价", 《合肥工业大学学报(自然科学版)》 * |
朱瑞磊: "一种关于纺织品图像带状瑕疵的检测方法", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
田光 等: "相位相关法及其在偏振图像配准中的应用研究", 《大气与环境光学学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110044931A (en) * | 2019-04-23 | 2019-07-23 | 华中科技大学 | A kind of detection device on bend glass surface and internal flaw |
CN111328370A (en) * | 2019-04-23 | 2020-06-23 | 合刃科技(深圳)有限公司 | Surface defect detection system and method |
CN110044931B (en) * | 2019-04-23 | 2021-03-26 | 华中科技大学 | Detection apparatus for curved surface glass surface and internal defect |
CN111328370B (en) * | 2019-04-23 | 2023-05-02 | 合刃科技(深圳)有限公司 | Surface defect detection system and method |
CN110412035A (en) * | 2019-06-13 | 2019-11-05 | 合刃科技(武汉)有限公司 | A kind of high reflecting material surface inspecting method and system |
CN114324369A (en) * | 2022-03-11 | 2022-04-12 | 北京新研创能科技有限公司 | System and method for detecting scratches on surface of bipolar plate |
CN114324369B (en) * | 2022-03-11 | 2022-06-07 | 北京新研创能科技有限公司 | System and method for detecting scratches on surface of bipolar plate |
CN114460104A (en) * | 2022-04-13 | 2022-05-10 | 清远市通用皮具配件有限公司 | Zipper tooth arrangement injection molding quality detection system and method based on machine vision |
CN114460104B (en) * | 2022-04-13 | 2022-07-12 | 清远市通用皮具配件有限公司 | Zipper tooth arrangement injection molding quality detection system and method based on machine vision |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104899890A (en) | Detection method based on polarized light imaging technology | |
Li et al. | No-reference image quality assessment with shearlet transform and deep neural networks | |
Zalama et al. | Road crack detection using visual features extracted by Gabor filters | |
CN102288613B (en) | Surface defect detecting method for fusing grey and depth information | |
Huang et al. | A pavement crack detection method combining 2D with 3D information based on Dempster‐Shafer theory | |
Guo et al. | Defect Detection in Tire X‐Ray Images Using Weighted Texture Dissimilarity | |
Chen et al. | Structure-aware-based crack defect detection for multicrystalline solar cells | |
CN103759676A (en) | Non-contact type workpiece surface roughness detecting method | |
CN104574296B (en) | A kind of method for polarizing the m ultiwavelet fusion treatment picture for removing haze | |
CN103729842A (en) | Fabric defect detection method based on local statistical characteristics and overall significance analysis | |
CN204903413U (en) | Detection device based on polarized light imaging technique | |
CN103839243B (en) | Multi-channel satellite cloud picture fusion method based on Shearlet conversion | |
Hoseini et al. | Fabric defect detection using auto-correlation function | |
CN106157232A (en) | A kind of general steganalysis method of digital picture characteristic perception | |
CN102651074A (en) | Texture feature-based printed paper identification method | |
Chen et al. | Multidirectional edge detection based on gradient ghost imaging | |
CN106023150A (en) | Real-time objective quality evaluation method for digital hologram noise interference | |
CN116030056A (en) | Detection method and system for steel surface cracks | |
CN112508862B (en) | Method for enhancing magneto-optical image of crack by improving GAN | |
Teeravech et al. | Discovering repetitive patterns in facade images using a RANSAC-style algorithm | |
Pei et al. | Image enhancement for on-site X-ray nondestructive inspection of reinforced concrete structures | |
CN113092487A (en) | Rapid detection method for full-surface flaws of safety belt roller | |
Cui et al. | Image enhancement-based detection of concrete cracks under turbid water bodies | |
Guo et al. | A novel laser stripe center extraction method for pavement rut detection | |
Xu et al. | Crack identification algorithm based on MASK Dodging principle and deep learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150909 |