CN102842120A - Image blurring degree detection method based on supercomplex wavelet phase measurement - Google Patents
Image blurring degree detection method based on supercomplex wavelet phase measurement Download PDFInfo
- Publication number
- CN102842120A CN102842120A CN2012102995528A CN201210299552A CN102842120A CN 102842120 A CN102842120 A CN 102842120A CN 2012102995528 A CN2012102995528 A CN 2012102995528A CN 201210299552 A CN201210299552 A CN 201210299552A CN 102842120 A CN102842120 A CN 102842120A
- Authority
- CN
- China
- Prior art keywords
- supercomplex
- coefficient
- wavelet
- image
- phase
- 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
Images
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
The invention discloses an image blurring degree detection method based on supercomplex wavelet phase measurement, and relates to an image resolution distinguishing method in image processing field. In order to overcome the fault of narrow detection range of an existing method, the image blurring degree detection method comprises following steps: first, transforming a noise image in a supercomplex wavelet mode to obtain a supercomplex wavelet phase coefficient; second, converting the supercomplex wavelet phase coefficient in a real part - imaginary part representation to a representation of amplitude - phase; third, counting the distribution character of a phase coefficient, and calculating distribution variance. The distribution variance is used to describe the distribution character of the supercomplex wavelet phase coefficient, so that texture variation between a clear image and a blurred image is effectively described, and the detection of blurring degree between different images can be achieved.
Description
Technical field
The present invention relates to a kind of image definition method of discrimination of image processing field, be specifically related to a kind of image blurring degree detecting method of measuring based on the supercomplex Wavelet Phase.
Background technology
Image definition/fuzzy the detection because its various vision is used is paid close attention to and receive widely, merges or the like such as the self-adaptation of the automatic focusing of digital micro-analysis imaging, remote Sensing Image Analysis, many focus charts picture.From a large amount of images that contains different fog-levels, select image automatically, can help the researcher to reduce very big workload for further graphical analysis with high definition.Because the different depth of field is difficult in the piece image each object is realized focusing on,, obtain the image of better visual effect so we can judge focal zone for many focus charts picture, and merge their clear area.These all need stable image blur operator.
Some article research image blur operators or similar notion have been arranged, measured like fuzzy detection, focusing etc.Most of operators comprise variance, image gradient, RMSE, PSNR etc., suppose that all the radio-frequency component of image has been decayed by point spread function.But these criterion functions have been ignored The noise, and noise can cause the variation of image radio-frequency component, but do not influence the integral image sharpness.The width that the wavelet coefficient of Tian J. and Chen J. discovery blurred picture distributes is narrower than the picture rich in detail, and these characteristics impel them to use the Laplacian mixture model with local auto-adaptive property that the marginal distribution of wavelet coefficient is carried out modeling.All above-mentioned operators of mentioning are all relevant with picture material itself, and for example, for the different images with same fog-level, the result that these operators calculate is different.
The supercomplex wavelet transformation (Hypercomplex Wavelet Transform, HWT), as a kind of new image analysis tool; Than wavelet transform; Have superior character, move constantly such as coefficient is approximate, phase coefficient provides more abundant image texture information etc.1999, B ü low made up for the first time the two-dimensional analysis signal in his PhD dissertation, and the concept extension of supercomplex Fourier transform to supercomplex Gabor conversion, is the predecessor of HWT, and uses it for difference estimation and Texture Segmentation.People such as Chan W. used in 2008 two-dimentional Hilbert transform (Hilbert Transform, HT) with the notion of analytic signal with the concept extension of dual-tree complex wavelet to onlapping number field, and be applied to difference estimation.HWT provides the subband of Quaternion Representation, and its coefficient can convert the form of an amplitude and three phase places through algebraic operation into---and its amplitude has to be similar to moves invariant feature, and three two-dimensional phases have comprised the geological information of describing local two-dimensional structure.
Represent ability because the HWT phase coefficient has texture, a kind of new approach is provided so statistics supercomplex Wavelet Phase coefficient is distributed as the detection of fog-level.
Summary of the invention
The objective of the invention is to overcome the existing narrower deficiency of method sensing range, provide a kind of based on supercomplex wavelet transformation (Hypercomplex Wavelet Transform, HWT) the image blur detection method of phase measurement.It can bring into play the ability that supercomplex Wavelet Phase texture is represented, sets up the corresponding relation between the distribution of different fog-levels and phase coefficient, thereby reaches the sensing range of widening image blurring degree.
The objective of the invention is to realize through following technical scheme:
Blurred picture is carried out HWT, convert the wavelet coefficient that obtains the representation of amplitude-phase place into, calculate the distribution statistics characteristic of phase coefficient, be used for the detected image blur level.
Process flow diagram of the present invention is as shown in Figure 1, and concrete steps are following:
Step 1: noise image is carried out HWT, obtain the supercomplex wavelet coefficient.This step is decomposed the multiresolution that image carries out higher-dimension, can excavate the inherent geometry of image, can effectively suppress noise, keeps the image border.
The supercomplex analytic signal is by signal itself
and its section H T(
;
) and complete HT(
) constitute
Wherein,
;
;
is hypercomplex three imaginary parts, and
is the structure result of supercomplex analytic signal.
Two dimension HT is equivalent to the HT that the row and column of matrix is done one dimension respectively.Can divide small echo for two dimension; The one dimension HT that
considers wavelet function is to
with yardstick equation
, and two-dimentional HWT can remember that work can divide the form of small echo product:
;
Through image and formula (2)-(5) are carried out convolution, just can obtain the HWT coefficient of noise image.Wherein,
is the 2-d wavelet wave filter;
is supercomplex small echo low frequency coefficient;
is supercomplex small echo horizontal component coefficient;
is supercomplex small echo vertical component coefficient,
be supercomplex small echo diagonal components coefficient.
Step 2: the supercomplex wavelet coefficient that real part-imaginary part is represented converts amplitude-phase place representation into; Phase place can be expressed the texture information of image; When handling amplitude coefficient, the image texture information that receives the less influence of noise can be kept fully like this.
The supercomplex wavelet coefficient that is obtained by supercomplex small echo and image convolution can be expressed as:
It is the expression mode of real part-imaginary part; Wherein,
is the supercomplex algebraic symbol;
is supercomplex real part coefficient;
is supercomplex imaginary part i component coefficient;
is supercomplex imaginary part j component coefficient;
is supercomplex imaginary part k component coefficient;
,
is three imaginary parts.
According to the computation rule of supercomplex algebraically, formula (6) can convert following form into:
Wherein:
;
is hypercomplex amplitude;
is hypercomplex phase place, and
is exponent sign.The calculation procedure at three phasing degree
is following:
Step 3: the distribution of statistics phase coefficient, calculate distribution variance.
The distribution of statistics phase coefficient; In its minimum value and peaked interval, is minizone with 0.01 index variation with supercomplex Wavelet Phase coefficient, adds up the number of each minizone supercomplex wavelet coefficient; The total number of image pixel is known, and then the probability that obtains distributing.
Fig. 2-shown in Figure 9 is phase coefficient
distribution histogram; Image behind
process 7*7 window Gaussian Blur is narrower than the picture rich in detail; Because the pixel intensity of blurred picture presents consistance, and the high frequency of texture information has been weakened.Therefore, need an operator to describe the distribution histogram shape directly related with image definition.Variance is being controlled the width of Gaussian distribution, and the present invention utilizes the notion of variance to describe the form parameter of image blurring degree.The image blur operator definitions is:
Wherein,
representes high frequency;
is the variance of QWT coefficient, and subscript
is correspondingly represented the level and the vertical component of phase place.The size of this operator is used for describing the fog-level of image, is worth more for a short time, and representative image is fuzzy more.The computing formula of variance
is:
Wherein,
is supercomplex Wavelet Phase coefficient;
is supercomplex wavelet coefficient phase place average;
is the number of supercomplex Wavelet Phase coefficient, following table
mark coefficient index.
The present invention utilizes variance to describe the distribution characteristics of supercomplex Wavelet Phase coefficient, has described the texture variations between clear and the blurred picture effectively, can realize the detection of fog-level between the different images.Compared with prior art have following advantage:
1) image blur detection method proposed by the invention is utilized the supercomplex wavelet transformation, and compares based on method of wavelet, has measurement stability, wideer sensing range.
2) the present invention introduces the supercomplex wavelet transformation and carries out image denoising; In the supercomplex wavelet field; The present invention proposes phase place high-frequency information measuring method, with data by MoM and MEI, has lower complexity and measures consistance; Blur level result of calculation and image itself is irrelevant, has approximate result of calculation for the image of the identical fog-level of different images.
Description of drawings
Fig. 1 is the image blur detection method process flow diagram of measuring based on the supercomplex Wavelet Phase;
Fig. 2 is the supercomplex small echo low frequency coefficient of the picture rich in detail of phase place
expression;
Fig. 3 is the supercomplex small echo horizontal component coefficient of the picture rich in detail of phase place
expression;
Fig. 4 is the supercomplex small echo vertical component coefficient of the picture rich in detail of phase place
expression;
Fig. 5 is the supercomplex small echo diagonal components coefficient of the picture rich in detail of phase place
expression;
Fig. 6 is the supercomplex small echo low frequency coefficient of the image after the bluring of phase place
expression;
Fig. 7 is the supercomplex small echo horizontal component coefficient of the image after the bluring of phase place
expression;
Fig. 8 is the supercomplex small echo vertical component coefficient of the image after the bluring of phase place
expression;
Fig. 9 is the supercomplex small echo diagonal components coefficient of the image after the bluring of phase place
expression;
Figure 10 is a picture rich in detail;
Figure 11 is the image of Figure 10 after fuzzy;
Figure 12 is that Figure 10 is through the result of calculation after bluring in various degree;
Figure 13 is a test pattern one;
Figure 14 is a test pattern two;
Figure 15 is a test pattern three;
Figure 16 is a test pattern four;
Figure 17 is a test pattern five;
Figure 18 is a test pattern six;
Figure 19 is a test pattern seven;
Figure 20 is a test pattern eight;
Figure 21 is a test pattern nine;
Figure 22 is a test pattern ten;
Figure 23 is the result of calculation of the different fog-levels of ten width of cloth test patterns.
Embodiment
Set forth embodiment of the present invention through standard testing image simulation example below:
Execution in step one: the image to picture rich in detail and its after fuzzy carries out HWT, obtains the supercomplex wavelet coefficient that real part-imaginary part is represented.
With Figure 10-11 is example, and image resolution ratio is 256*256, and blurred picture is the result who picture rich in detail is done Gauss's mask convolution.
Execution in step two: the supercomplex wavelet coefficient image transitions that real part-imaginary part that step 1 is obtained is represented is amplitude-phase place representation.
The phase component that is applied to the blur level detection that Figure 10-11 is corresponding is seen Fig. 2-9.
Execution in step three: the distribution character of statistics phase coefficient, calculate distribution variance.
Because the phase component that supercomplex wavelet amplitude-phase place is represented has the ability that texture is represented, the texture of blurred picture can be affected, so the result of calculation behind two width of cloth image normalizations shown in Figure 10 of the present invention-11 is respectively 1 and 0.5335.Figure 12 is the result of calculation of the image after bluring through the template of different sizes; Horizontal ordinate is a template size, and ordinate is a variance result of calculation, can find out the increase along with fog-level; Result of calculation diminishes, and method of the present invention can active zone be told the image of different fog-levels.
Figure 13-22 is 10 width of cloth test patterns, measures, can obtain the result of calculation of Figure 23 through utilizing method of the present invention to carry out blur level behind same degree fuzzy.Be not more than in 9 in template size, have similar result of calculation between the different images for same fog-level; Between the different fog-level different images, also can pass through the different fog-level of result of calculation differentiate between images.
Claims (4)
1. the image blurring degree detecting method of measuring based on the supercomplex Wavelet Phase is characterized in that said method comprises the steps:
Wherein:
is the supercomplex algebraic symbol;
is supercomplex real part coefficient;
is supercomplex imaginary part i component coefficient;
is supercomplex imaginary part j component coefficient;
is supercomplex imaginary part k component coefficient;
,
is three imaginary parts;
Step 2: the supercomplex wavelet coefficient that real part-imaginary part is represented converts amplitude-phase place representation into:
;
Step 3: the distribution of statistics phase coefficient, calculate distribution variance, utilize variance to describe the form parameter of image blurring degree.
2. the image blurring degree detecting method of measuring based on the supercomplex Wavelet Phase according to claim 1; It is characterized in that in the said step 2 that the calculation procedure of phasing degree
is following:
3. the image blurring degree detecting method of measuring based on the supercomplex Wavelet Phase according to claim 1; It is characterized in that in the said step 3; The image blur operator definitions is:
; Be worth more for a short time, representative image is fuzzy more;
4. the image blurring degree detecting method of measuring based on the supercomplex Wavelet Phase according to claim 1; It is characterized in that in the said step 3 that the computing formula of variance
is:
;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210299552.8A CN102842120B (en) | 2012-08-22 | 2012-08-22 | Image blurring degree detection method based on supercomplex wavelet phase measurement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210299552.8A CN102842120B (en) | 2012-08-22 | 2012-08-22 | Image blurring degree detection method based on supercomplex wavelet phase measurement |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102842120A true CN102842120A (en) | 2012-12-26 |
CN102842120B CN102842120B (en) | 2015-04-01 |
Family
ID=47369446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210299552.8A Active CN102842120B (en) | 2012-08-22 | 2012-08-22 | Image blurring degree detection method based on supercomplex wavelet phase measurement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102842120B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106647A (en) * | 2013-03-06 | 2013-05-15 | 哈尔滨工业大学 | Multi-focal-point image fusion method based on quaternion wavelet and region segmentation |
CN103353985A (en) * | 2013-03-27 | 2013-10-16 | 西华大学 | Measurement method for image Gaussian Blur |
CN104637046A (en) * | 2013-11-13 | 2015-05-20 | 索尼公司 | Image detection method and device |
CN105118057A (en) * | 2015-08-18 | 2015-12-02 | 江南大学 | Image sharpness evaluation method based on quaternion wavelet transform amplitudes and phase positions |
CN106485703A (en) * | 2016-09-30 | 2017-03-08 | 杭州电子科技大学 | Fuzzy detection method based on image gradient dct transform |
CN106504221A (en) * | 2016-10-13 | 2017-03-15 | 江南大学 | Based on the Medical image fusion new method that quaternion wavelet converts context mechanism |
CN106934804A (en) * | 2017-03-13 | 2017-07-07 | 重庆贝奥新视野医疗设备有限公司 | Approach for detecting image sharpness and device |
CN111860057A (en) * | 2019-04-29 | 2020-10-30 | 北京眼神智能科技有限公司 | Face image blurring and living body detection method and device, storage medium and equipment |
CN116847209A (en) * | 2023-08-29 | 2023-10-03 | 中国测绘科学研究院 | Log-Gabor and wavelet-based light field full-focusing image generation method and system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980286A (en) * | 2010-11-12 | 2011-02-23 | 西安电子科技大学 | Method for reducing speckles of synthetic aperture radar (SAR) image by combining dual-tree complex wavelet transform with bivariate model |
WO2012066568A1 (en) * | 2010-11-15 | 2012-05-24 | Indian Institute Of Technology | An improved ultrasound imaging method/technique for speckle reduction/suppression in an improved ultra sound imaging system |
-
2012
- 2012-08-22 CN CN201210299552.8A patent/CN102842120B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980286A (en) * | 2010-11-12 | 2011-02-23 | 西安电子科技大学 | Method for reducing speckles of synthetic aperture radar (SAR) image by combining dual-tree complex wavelet transform with bivariate model |
WO2012066568A1 (en) * | 2010-11-15 | 2012-05-24 | Indian Institute Of Technology | An improved ultrasound imaging method/technique for speckle reduction/suppression in an improved ultra sound imaging system |
Non-Patent Citations (2)
Title |
---|
EDUARDO BAYRO-CORROCHANO: "Multi-resolution image analysis using the quaternion wavelet transform", 《NUMERICAL ALGORITHMS》 * |
PING YANG ET AL.: "19-element sensor-less adaptive optical system based on modified hill-climbing and genetic Algorithms", 《PROC. OF SPIE》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106647B (en) * | 2013-03-06 | 2015-08-19 | 哈尔滨工业大学 | Based on the Multi-focal-point image fusion method of quaternion wavelet and region segmentation |
CN103106647A (en) * | 2013-03-06 | 2013-05-15 | 哈尔滨工业大学 | Multi-focal-point image fusion method based on quaternion wavelet and region segmentation |
CN103353985A (en) * | 2013-03-27 | 2013-10-16 | 西华大学 | Measurement method for image Gaussian Blur |
CN103353985B (en) * | 2013-03-27 | 2016-05-25 | 西华大学 | A kind of Measurement Method of image Gaussian Blur |
CN104637046B (en) * | 2013-11-13 | 2020-03-13 | 索尼公司 | Image detection method and device |
CN104637046A (en) * | 2013-11-13 | 2015-05-20 | 索尼公司 | Image detection method and device |
CN105118057A (en) * | 2015-08-18 | 2015-12-02 | 江南大学 | Image sharpness evaluation method based on quaternion wavelet transform amplitudes and phase positions |
CN106485703B (en) * | 2016-09-30 | 2019-05-21 | 杭州电子科技大学 | Fuzzy detection method based on image gradient dct transform |
CN106485703A (en) * | 2016-09-30 | 2017-03-08 | 杭州电子科技大学 | Fuzzy detection method based on image gradient dct transform |
CN106504221A (en) * | 2016-10-13 | 2017-03-15 | 江南大学 | Based on the Medical image fusion new method that quaternion wavelet converts context mechanism |
CN106504221B (en) * | 2016-10-13 | 2019-01-25 | 江南大学 | Method of Medical Image Fusion based on quaternion wavelet transformation context mechanism |
CN106934804A (en) * | 2017-03-13 | 2017-07-07 | 重庆贝奥新视野医疗设备有限公司 | Approach for detecting image sharpness and device |
CN106934804B (en) * | 2017-03-13 | 2019-12-13 | 重庆贝奥新视野医疗设备有限公司 | Image definition detection method and device |
CN111860057A (en) * | 2019-04-29 | 2020-10-30 | 北京眼神智能科技有限公司 | Face image blurring and living body detection method and device, storage medium and equipment |
CN116847209A (en) * | 2023-08-29 | 2023-10-03 | 中国测绘科学研究院 | Log-Gabor and wavelet-based light field full-focusing image generation method and system |
CN116847209B (en) * | 2023-08-29 | 2023-11-03 | 中国测绘科学研究院 | Log-Gabor and wavelet-based light field full-focusing image generation method and system |
Also Published As
Publication number | Publication date |
---|---|
CN102842120B (en) | 2015-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102842120A (en) | Image blurring degree detection method based on supercomplex wavelet phase measurement | |
CN111062973B (en) | Vehicle tracking method based on target feature sensitivity and deep learning | |
CN104835175B (en) | Object detection method in a kind of nuclear environment of view-based access control model attention mechanism | |
CN102609701B (en) | Remote sensing detection method based on optimal scale for high-resolution SAR (synthetic aperture radar) | |
CN104966085B (en) | A kind of remote sensing images region of interest area detecting method based on the fusion of more notable features | |
CN102999886B (en) | Image Edge Detector and scale grating grid precision detection system | |
CN104408482B (en) | A kind of High Resolution SAR Images object detection method | |
Espinal et al. | Wavelet-based fractal signature analysis for automatic target recognition | |
CN108537787B (en) | Quality judgment method for face image | |
CN104899888B (en) | A kind of image sub-pixel edge detection method based on Legendre squares | |
CN109214439A (en) | A kind of infrared image icing River detection method based on multi-feature fusion | |
CN108257125B (en) | Depth image quality non-reference evaluation method based on natural scene statistics | |
CN110189375A (en) | A kind of images steganalysis method based on monocular vision measurement | |
CN106682678A (en) | Image angle point detection and classification method based on support domain | |
CN104021567B (en) | Based on the fuzzy altering detecting method of image Gauss of first numeral law | |
CN108550145A (en) | A kind of SAR image method for evaluating quality and device | |
CN105809173A (en) | Bionic vision transformation-based image RSTN (rotation, scaling, translation and noise) invariant attributive feature extraction and recognition method | |
CN105550694A (en) | Method for measurement of fuzzy degree of face image | |
CN104217430B (en) | Image significance detection method based on L1 regularization | |
CN103914829B (en) | Method for detecting edge of noisy image | |
CN104637060B (en) | A kind of image partition method based on neighborhood principal component analysis-Laplce | |
Junwu et al. | An infrared and visible image fusion algorithm based on LSWT-NSST | |
Jia et al. | Fabric defect inspection based on lattice segmentation and lattice templates | |
Hu et al. | On the use of joint sparse representation for image fusion quality evaluation and analysis | |
Lijun et al. | Geo‐information mapping improves Canny edge detection method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200326 Address after: 150001 No. 118 West straight street, Nangang District, Heilongjiang, Harbin Patentee after: Harbin University of technology high tech Development Corporation Address before: 150000 Harbin, Nangang, West District, large straight street, No. 92 Patentee before: HARBIN INSTITUTE OF TECHNOLOGY |