CN106023190A - Automatic focusing image focusing degree evaluation method - Google Patents

Automatic focusing image focusing degree evaluation method Download PDF

Info

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
Application number
CN201610325972.7A
Other languages
Chinese (zh)
Other versions
CN106023190B (en
Inventor
冯华君
王烨茹
徐之海
李奇
陈跃庭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201610325972.7A priority Critical patent/CN106023190B/en
Publication of CN106023190A publication Critical patent/CN106023190A/en
Application granted granted Critical
Publication of CN106023190B publication Critical patent/CN106023190B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image 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

A kind of evaluation methodology of auto-focusing image quasi-Jiao degree
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:
P C ( x , y ) = Σ i = 1 m Σ j = 1 n W i ( x , y ) f p { A j i ( x , y ) [ Δφ j i ( x , y ) ] - T i } Σ i = 1 m Σ i = 1 n A j i ( x , y ) - - - ( 1 )
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:
P C ( x , y ) = Σ i = 1 m Σ j = 1 n W i ( x , y ) f p { A j i ( x , y ) [ Δφ j i ( x , y ) ] - T i } Σ i = 1 m Σ i = 1 n A j i ( x , y ) - - - ( 1 )
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,
G = G x 2 + G y 2 - - - ( 3 )
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:
P C ( x , y ) = Σ i = 1 m Σ j = 1 n W i ( x , y ) f p { A j i ( x , y ) [ Δφ j i ( x , y ) ] - T i } Σ i = 1 m Σ i = 1 n A j i ( x , y ) - - - ( 1 )
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.
CN201610325972.7A 2016-05-16 2016-05-16 A kind of evaluation method of the quasi- burnt degree of auto-focusing image Expired - Fee Related CN106023190B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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