CN106023190A - Automatic focusing image focusing degree evaluation method - Google Patents
Automatic focusing image focusing degree evaluation method Download PDFInfo
- Publication number
- CN106023190A CN106023190A CN201610325972.7A CN201610325972A CN106023190A CN 106023190 A CN106023190 A CN 106023190A CN 201610325972 A CN201610325972 A CN 201610325972A CN 106023190 A CN106023190 A CN 106023190A
- Authority
- CN
- China
- Prior art keywords
- focusing
- image
- sigma
- phase
- gradient
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012634 optical imaging Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research 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/0002—Inspection of images, e.g. flaw detection
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses an automatic focusing image focusing degree evaluation method. The method adopts a phase consistency image and gradient image combined method, solves the problem of peak collapse caused when a phase consistency method is adopted as a focusing evaluation function, has good expression performance no matter in a normal scene, or a scene, where signal-to-noise ratio is low or brightness is changed, is high in practical value, and has the advantages of good unbiasedness, unimodality, sensitivity and stability.
Description
Technical field
The present invention relates to image quality evaluating method, particularly relate to the evaluation side of a kind of auto-focusing image quasi-Jiao degree
Method.
Background technology
Along with digital picture and the fast development of multimedia technology, various types of optical imaging systems also get more and more,
Autofocus Technology plays vital effect in optical imaging systems.Autofocus Technology mainly has two kinds, is respectively
Active focusing and passive type focusing.Actively focusing needs to add extra range-measurement system in systems, adds the one-tenth of system
Basis and the complexity of system, simultaneously for having glass or having the scene of shelter, range finding there will be problem, therefore, active
Focusing technology is replaced by self-action focusing technology at leisure.
Passive focusing is to judge out of focus degree by the fog-level of image obtained in analysis focus process, thus refers to
Lead focus process.This Atomatic focusing method is to utilize certain Digital Image Processing algorithm, obtains image definition eigenvalue,
And control motor according to this eigenvalue, to drive eyeglass or imageing sensor to change position, do corresponding adjustment, until this
Till one eigenvalue meets a certain condition made an appointment.Therefore, the auto focus evaluation function of superperformance is selected to judge
The fog-level of image thus to guide autofocus system to focus exactly be the key of problem.Weight due to focusing evaluation function
The property wanted, the evaluation of image definition has become a popular research field.Preferably focusing evaluation function should have unbiased
Property, unimodality, susceptiveness, stability, some the most also should to scene light and shade change and noise have preferably
Robustness.Existing auto focus evaluation function is very sensitive to the brightness and contrast of scene, because its principle is utilization figure
The intensity of picture judges out of focus degree.And along with imaging system popularizing in people live, people like at various environment
Lower use focusing system, uses the situation of focusing system also to get more and more under or scene that signal to noise ratio is relatively low inconsistent in brightness.
In consideration of it, still have the evaluation function of good behaviour performance to become the urgent needs of people this kind of scene.
Summary of the invention
Present invention aims to use phase equalization method to have what peak value subsided as focusing evaluation function
Problem, proposes the evaluation methodology of a kind of auto-focusing image quasi-Jiao degree.
It is an object of the invention to be achieved through the following technical solutions: the evaluation side of a kind of auto-focusing image quasi-Jiao degree
Method, the method comprises the following steps:
(1) image I is obtained;
(2) from image I, focusing window region I is chosenc;
(3) utilize log-Gabor small echo to calculate focusing window region IcPhase equalization collection of illustrative plates PC, each pixel
Phase equalization is defined as:
Wherein, m and n represents the direction yardstick number of Log-Gabor wavelet filter, A respectivelyji(x, y) and φji(x,y)
The amplitude being illustrated respectively in the jth logarithmic scale of i-th direction and local phase deviation.TiIt is estimating on i-th direction
Meter noise, Wi(x y) is the weighting function on i-th direction, fp(u)=(u > 0) u.
Δ φ (x, y)=cos [φr(x,y)-φ(x,y)]-|sin[φr(x,y)-φ(x,y)]| (2)
(x is y) at (x, y) the local phase weighted mean of position, φ to φr(x, y) represent local phase (x, y)
The r level Fourier expansion series component of position;
(4) utilize gradient calculation method obtain phase equalization collection of illustrative plates PC gradient magnitude spectrum G:
(5) suing for peace all gradient magnitudes in G, as final focusing function evaluation of estimate, gradient magnitude summation is the highest,
Current image location is closer to focus position.
The beneficial effects of the present invention is: 1. solve scene by phase equalization collection of illustrative plates PC and there is brightness flop and make an uproar
Problem in terms of sound, further combined with gradient calculation method overcome phase equalization method as during evaluation function in accurate burnt position
Putting neighbouring curve and there is the problem subsided, the evaluation function curve finally given has good unbiasedness, unimodality, susceptiveness
And stability, practical value is higher, can be widely used in and include brightness flop, and the situations such as signal to noise ratio is low are interior various
Actual photographed scene.
Accompanying drawing explanation
Fig. 1 is the FB(flow block) of inventive method.
Fig. 2 is the schematic diagram of optional 5 width in one group of emulation out of focus sequence chart.
Fig. 3 is as the criterion burnt image.
Fig. 4 is the focusing window image chosen.
Fig. 5 is the phase equalization collection of illustrative plates of focusing window area image.
Fig. 6 is the gradient magnitude spectrum G of phase equalization collection of illustrative plates PC.
Fig. 7 is the evaluation function curve utilizing phase equalization method to obtain.
Fig. 8 is the evaluation function curve utilizing phase equalization and gradient associated methods to obtain.
Fig. 9 is the evaluation function curve utilizing existing method to obtain.
Detailed description of the invention
The evaluation methodology of the present invention a kind of auto-focusing image quasi-Jiao degree, the method uses phase equalization collection of illustrative plates and ladder
The degree method that combines of collection of illustrative plates, solves asking of having when using phase equalization method as focusing evaluation function that peak value subsides
Topic, the scene that or brightness relatively low for normal scene or signal to noise ratio changes all has good performance performance, practical valency
Being worth higher, range of application is relatively wide, has good unbiasedness, unimodality, susceptiveness, stability.
It is described in detail with example below in conjunction with the accompanying drawings:
Fig. 1 is the simple process block diagram of the inventive method.Below in conjunction with embodiment, the invention will be further described.
The present embodiment is as a example by one group of emulation out of focus sequence chart, and Fig. 2 is for arbitrarily to choose 5 width figures in one group of out of focus sequence chart
As schematic diagram, wherein (a)-(e) is respectively the 1st width, the 9th width, the 17th width, the 25th width, the 33rd width out of focus sequence analogous diagram, wherein
17th width is as the criterion burnt figure.Wherein, taper into from the 1st width~the 17th width defocusing amount, from the 17th width~the 33rd width defocusing amount gradually
Become big.
The present embodiment carries out defocusing amount evaluation initially with existing method to 33 width figures, owing to brightness is inconsistent, existing
Method can not truly reflect the defocusing amount of each figure, as shown in Figure 9.
The present embodiment have employed the method for simple employing phase equalization method and the present invention further to above-mentioned 33 width
Image is evaluated as follows (as a example by the 17th width image):
(1) image I is obtained as shown in Figure 3;
(2) from image I, choose the focusing window region I of 200 × 200 pixels as shown in Figure 4c;
(3) utilize log-Gabor small echo to calculate focusing window area image IcPhase equalization collection of illustrative plates PC, 2D scheme
In Xiang, the phase equalization of each pixel is defined as:
Wherein, m and n represents the direction yardstick number of Log-Gabor wavelet filter, A respectivelyji(x, y) and φji(x,y)
The amplitude being illustrated respectively in the jth logarithmic scale of i-th direction and local phase deviation.TiIt is estimating on i-th direction
Meter noise, Wi(x y) is the weighting function on i-th direction, fp(u)=(u > 0) u.
Δ φ (x, y)=cos [φr(x,y)-φ(x,y)]-|sin[φr(x,y)-φ(x,y)]| (2)
(x is y) at (x, y) the local phase weighted mean of position, φ to φr(x, y) represent local phase (x, y)
The r level Fourier expansion series component of position, phase equalization collection of illustrative plates PC is as shown in Figure 5;
(4) utilize gradient calculation method obtain phase equalization collection of illustrative plates PC gradient magnitude spectrum G, as shown in Figure 6,
Wherein, GxWith GyRepresent phase equalization collection of illustrative plates PC partial derivative in the horizontal and vertical directions respectively;
(5) all gradient magnitudes in G are sued for peace, as final focusing function evaluation of estimate.
Repeat the above steps carries out definition evaluation of estimate calculating to the 33 width out of focus sequence analogous diagram obtained, and obtains this group
The evaluation function curve of out of focus sequence focusing window, as shown in Figure 8.Preferably solve as shown in Figure 7 use merely phase place
The problem that peak value subsides is there is in coherence method as focusing evaluation function.
The present invention proposes the auto focus evaluation function computational methods that a kind of phase equalization is combined, the method with gradient
Solve employing phase equalization method and there is, as focusing evaluation function, the problem that peak value subsides, for normal scene the most also
Being that relatively low or brightness change the scene of signal to noise ratio all has good performance performance, practical value is higher, and range of application is relatively wide, tool
There are good unbiasedness, unimodality, susceptiveness, stability.
Claims (1)
1. the evaluation methodology of auto-focusing image quasi-Jiao degree, it is characterised in that the method comprises the following steps:
(1) image I is obtained;
(2) from image I, focusing window region I is chosenc;
(3) utilize log-Gabor small echo to calculate focusing window region IcPhase equalization collection of illustrative plates PC, the phase place of each pixel
Conformance definition is:
Wherein, m and n represents the direction yardstick number of Log-Gabor wavelet filter, A respectivelyji(x, y) and φji(x, y) respectively
Represent the amplitude in the jth logarithmic scale of i-th direction and local phase deviation.TiIt is that the estimation on i-th direction is made an uproar
Sound, Wi(x y) is the weighting function on i-th direction, fp(u)=(u > 0) u.
Δ φ (x, y)=cos [φr(x,y)-φ(x,y)]-|sin[φr(x,y)-φ(x,y)]| (2)
(x is y) at (x, y) the local phase weighted mean of position, φ to φr(x y) represents local phase in (x, y) position
R level Fourier expansion series component;
(4) utilize gradient calculation method obtain phase equalization collection of illustrative plates PC gradient magnitude spectrum G:
(5) suing for peace all gradient magnitudes in G, as final focusing function evaluation of estimate, gradient magnitude summation is the highest, currently
Picture position is closer to focus position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610325972.7A CN106023190B (en) | 2016-05-16 | 2016-05-16 | A kind of evaluation method of the quasi- burnt degree of auto-focusing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610325972.7A CN106023190B (en) | 2016-05-16 | 2016-05-16 | A kind of evaluation method of the quasi- burnt degree of auto-focusing image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106023190A true CN106023190A (en) | 2016-10-12 |
CN106023190B CN106023190B (en) | 2019-05-07 |
Family
ID=57097213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610325972.7A Expired - Fee Related CN106023190B (en) | 2016-05-16 | 2016-05-16 | A kind of evaluation method of the quasi- burnt degree of auto-focusing image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106023190B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108513115A (en) * | 2018-03-14 | 2018-09-07 | 深圳市点睛创视技术有限公司 | A kind of dynamic focusing window selection method and apparatus in optical projection system |
CN109801289A (en) * | 2019-01-25 | 2019-05-24 | 淮阴师范学院 | A kind of image Focus field emission array implementation method based on Gabor transformation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140247870A1 (en) * | 2011-04-28 | 2014-09-04 | Koninklijke Philips N.V. | Apparatuses and methods for hdr image encoding and decodng |
CN104680541A (en) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | Remote sensing image quality evaluation method based on phase congruency |
CN105243660A (en) * | 2015-09-16 | 2016-01-13 | 浙江大学 | Quality evaluation method of light source scene-containing automatic focusing image |
CN105354835A (en) * | 2015-10-16 | 2016-02-24 | 浙江工业大学 | Method for evaluating medical image quality in combination with phase consistency, gradient magnitude and structural prominence |
-
2016
- 2016-05-16 CN CN201610325972.7A patent/CN106023190B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140247870A1 (en) * | 2011-04-28 | 2014-09-04 | Koninklijke Philips N.V. | Apparatuses and methods for hdr image encoding and decodng |
CN104680541A (en) * | 2015-03-15 | 2015-06-03 | 西安电子科技大学 | Remote sensing image quality evaluation method based on phase congruency |
CN105243660A (en) * | 2015-09-16 | 2016-01-13 | 浙江大学 | Quality evaluation method of light source scene-containing automatic focusing image |
CN105354835A (en) * | 2015-10-16 | 2016-02-24 | 浙江工业大学 | Method for evaluating medical image quality in combination with phase consistency, gradient magnitude and structural prominence |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108513115A (en) * | 2018-03-14 | 2018-09-07 | 深圳市点睛创视技术有限公司 | A kind of dynamic focusing window selection method and apparatus in optical projection system |
CN108513115B (en) * | 2018-03-14 | 2020-08-18 | 深圳市点睛创视技术有限公司 | Dynamic focusing window selection method and device in projection system |
CN109801289A (en) * | 2019-01-25 | 2019-05-24 | 淮阴师范学院 | A kind of image Focus field emission array implementation method based on Gabor transformation |
Also Published As
Publication number | Publication date |
---|---|
CN106023190B (en) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9438816B2 (en) | Adaptive image acquisition for multiframe reconstruction | |
CN102314683B (en) | Computational imaging method and imaging system based on nonplanar image sensor | |
CN113140011B (en) | Infrared thermal imaging monocular vision distance measurement method and related components | |
CN102053339B (en) | Automatic focusing method | |
CN106249325A (en) | A kind of bionical quick focus adjustment method of vision based on liquid lens | |
CN104618664A (en) | Adaptive integral time based embedded infrared image superframe processing method | |
US20070286517A1 (en) | Method and apparatus for multifocus digital image restoration using image integration technology | |
CN102968800B (en) | A kind of evaluation method of image definition | |
CN101943839A (en) | Integrated automatic focusing camera device and definition evaluation method | |
CN104914446A (en) | Three-dimensional distance image time domain real-time denoising method based on photon counting | |
CN112019751B (en) | Calibration information based automatic focusing method | |
JP2013531268A (en) | Measuring distance using coded aperture | |
CN106526839B (en) | Synchronous wavefront-free self-adaptive optical system based on mode | |
CN103487927A (en) | Automatic focusing method of microscope | |
CN103916599A (en) | Rapid focusing method for remote sensing camera exterior imaging | |
CN109064508A (en) | A kind of detecting a laser beam method and system | |
CN106023190A (en) | Automatic focusing image focusing degree evaluation method | |
CN108765505A (en) | The coding of infrared imaging perceives matrix estimation method | |
CN102883103B (en) | Auto-focusing method and photographic device | |
CN102170521A (en) | Non-uniform-sampling-window-based automatic focusing method for digital still camera | |
CN104813217B (en) | For the method for the passive single channel imager for being designed to estimate the depth of field | |
CN106888344A (en) | Camera module and its inclined acquisition methods of image planes and method of adjustment | |
CN112113661B (en) | Deep learning type snapshot spectrum imaging device and detection method thereof | |
CN103873773B (en) | Primary-auxiliary synergy double light path design-based omnidirectional imaging method | |
CN101561267B (en) | Distance-measuring device and 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190507 |