WO2019148739A1 - Comprehensive processing method and system for blurred image - Google Patents

Comprehensive processing method and system for blurred image Download PDF

Info

Publication number
WO2019148739A1
WO2019148739A1 PCT/CN2018/091164 CN2018091164W WO2019148739A1 WO 2019148739 A1 WO2019148739 A1 WO 2019148739A1 CN 2018091164 W CN2018091164 W CN 2018091164W WO 2019148739 A1 WO2019148739 A1 WO 2019148739A1
Authority
WO
WIPO (PCT)
Prior art keywords
edge
image
blur
blurred
optimal
Prior art date
Application number
PCT/CN2018/091164
Other languages
French (fr)
Chinese (zh)
Inventor
白海玲
Original Assignee
上海康斐信息技术有限公司
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 上海康斐信息技术有限公司 filed Critical 上海康斐信息技术有限公司
Publication of WO2019148739A1 publication Critical patent/WO2019148739A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4053Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20061Hough transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • G06T2207/20192Edge enhancement; Edge preservation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • G06T2207/20201Motion blur correction

Definitions

  • the present invention relates to the field of image processing technologies, and in particular, to a fuzzy image integration processing method and system.
  • Image acquisition may result in inaccurate focus, or relative motion of the camera and the camera, camera distortion, air diffraction, etc., resulting in unsatisfactory images, losing useful information in the image, that is, image degradation. . It hinders the subsequent processing analysis, especially in the face of some non-replicable scenes. Therefore, for such cases, restoring the original appearance of the scene is an urgent problem to be solved.
  • the actual design can add optical devices, improve transmission equipment and other hardware means such as electronic image shift compensation, optical image shift compensation, mechanical image shift compensation to improve the quality of the captured image, but because the device process is too complicated, The impact of many factors such as high cost and long adjustment time is not universal. Therefore, using image restoration technology to improve image quality from image compensation is a correct, low-cost, low-cost option. .
  • PSF point spread function
  • the degradation model can be used to reverse the degradation process by algorithm to achieve the restoration of the blurred image.
  • the most suitable restoration algorithm can not be selected to help the image restoration.
  • the point spread function estimation is not accurate, and the selected restoration algorithm is not appropriate, which will result in low quality.
  • Restoration not only can not get better fuzzy image restoration effect, but also make the blurred image more fuzzy, so in addition to the accurate estimation of the point spread function is important, it also needs to obtain accurate prior knowledge, that is, the image Causes of degradation and types of blurring, etc.
  • Image restoration algorithms are mathematically defined as a class of ill-posed inverse problems. Knowing the cause, the result is a positive problem, and the known result is the inverse problem; the image restoration algorithm is an algorithm that knows the reason for the reverse of the result.
  • the ill-posed inverse problem means that the reverse push process is very unstable, that is, it is affected by slight noise, which will cause very large interference to the final speculation, resulting in incorrect results. Therefore, the image is restored to minimize the impact of interference on the reverse process.
  • the image restoration algorithms commonly used in the prior art include RL filtering, constrained least squares filtering, Wiener filtering, and regular filtering.
  • Common fuzzy image types include motion blur and defocus blur.
  • the patent document disclosed in the publication No. CN104331871A discloses "an image deblurring method and apparatus", comprising: performing blur region detection on a to-be-processed image, determining a blur region image, and determining a blur type of the blur region, if the blur region If the blur type of the image is defocus blur, the defocus blur parameter estimation algorithm based on differential image autocorrelation is used to determine the defocus radius. If the blur type of the blur region image is motion blur, the motion blur based on cepstrum analysis The parameter estimation algorithm determines the blur direction and the blur scale, and substitutes the estimated parameters into the classical image restoration algorithm to obtain the restored image.
  • the patent uses different methods to estimate the parameters necessary for restoration according to different types of blur.
  • the disc model has a large limitation.
  • the motion blur parameter estimation algorithm based on cepstrum analysis is used to determine the blur direction and blur.
  • the scale does not provide perfect and realistic prior knowledge for image restoration; the classic image restoration algorithm includes Wiener filtering and LR filtering. Wiener filtering is not ideal when the peak signal-to-noise ratio of the blurred image is small, and the LR filtering is not ideal. It is sensitive to noise, and the restored image has a significant ringing effect, that is, the adaptation range is not wide enough.
  • the technical problem to be solved by the present invention is to provide a fuzzy image comprehensive processing method suitable for more real fuzzy image restoration, simple restoration method, strong pertinence, accurate restoration, high speed and high quality, in view of the above-mentioned deficiencies of the prior art. And system.
  • a method for comprehensively processing blurred images comprising the following steps:
  • S10 Identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
  • the edge kernel method is used to obtain the blur kernel; if the blurred image is motion blur, the adaptive preset algorithm obtains the blur kernel, and the preset algorithm includes at least one algorithm;
  • the method further includes the following steps:
  • S40 processing the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image.
  • acquiring the fuzzy kernel by using the edge method includes the following steps:
  • S202 Acquire a point spread equation according to the optimal edge edge image, where the point spread equation is the blur kernel.
  • step S201 includes the following steps:
  • the gradient values of the edge image are calculated one by one, and the optimal edge image is extracted based on the gradient values.
  • the calculating the gradient value of the edge image and extracting the optimal edge image according to the gradient value comprises the following steps:
  • the absolute value of the difference between the average pixel values of the two sides of the new edge edge point is counted as the gradient value of the edge edge image, and the edge edge image having the largest gradient value is extracted as the optimal edge edge image.
  • step S202 includes the following steps:
  • the point spread equation is calculated using the edge spread function.
  • step S20
  • the adaptive regularization method based on sparse prior is used to estimate the fuzzy kernel.
  • a fuzzy image integrated processing system comprising:
  • a fuzzy type identification module configured to identify a blur type of the blurred image according to an elongation of the blurred image spectrogram, the blur type including defocus blur and motion blur;
  • An estimation module configured to acquire a fuzzy kernel by using a blade edge method if the blurred image is out-of-focus blur; and obtain a fuzzy kernel by using an adaptive preset algorithm if the blurred image is motion blur, the preset algorithm includes at least one algorithm;
  • An image restoration module is configured to obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
  • system further includes:
  • a high resolution processing module that processes the sharp image using a super-resolution reconstruction technique to obtain a high resolution image
  • the estimation module includes:
  • a defocusing unit for extracting an edge-edge image based on a gradient criterion and acquiring a point-diffusion equation according to the optimal edge-edge image if the blurred image is defocused, and the point-diffusion equation is The blur kernel;
  • a motion blur unit is used to acquire a fuzzy kernel based on a sparse priori regularization method if the blurred image is motion blur.
  • the defocus blur unit includes:
  • An edge detection subunit configured to perform Canny edge detection on the blurred image to obtain an edge detection image
  • a line detection subunit configured to perform a Hough transform on the edge detection image to obtain a step edge image
  • An image intercepting subunit for intercepting a blade edge image of each edge centering on a center point of each edge of the step edge image, the size of the edge edge image being according to the defocused image size and The degree of blur selects a predetermined size to be formed;
  • An image extraction subunit configured to calculate a gradient value of the edge edge image one by one, and extract an optimal edge edge image according to the gradient value;
  • a linear fitting sub-unit for linearly fitting edges in the optimal edge-edge image by using a least squares method to obtain an optimal edge line
  • a scatter plot sub-unit configured to obtain a vertical distance of each pixel point in the optimal edge-edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate, forming a dispersion Dot map
  • An edge diffusion function sub-unit for linearly fitting the scatter plot with a Fermi function to obtain an edge spread function
  • a calculation subunit is configured to calculate a point spread equation using the edge spread function.
  • Determining the type of fuzzy image according to the elongation is not only simple, but also suitable for the clear definition of the fuzzy image type in the actual real scene, which can ensure the accuracy of the fuzzy type identification and help to accurately estimate the fuzzy kernel with different degradation mechanisms. .
  • the edge edge method is used to adaptively extract the optimal edge image of the defocused image, and the optimal edge image is obtained based on the gradient criterion.
  • it is beneficial to improve the speed of defocus image restoration, on the other hand, it is beneficial to Avoid blindly extracting the edge image, reduce human interference, and improve the accuracy of image restoration.
  • the fuzzy kernel estimated from the optimal edge-edge image is also closer to the true degradation model, and the recovery speed is faster and more accurate.
  • the adaptive sparse priori-based regularization method estimates the fuzzy kernel of the motion blurred image, and iterates the target multiple times until convergence, which results in more stable and ideal results, thus solving the ill-posed problem in the image restoration problem.
  • the super Laplacian prior deconvolution algorithm is used to obtain clear images.
  • the super Laplace a priori is used as a regular term.
  • the accuracy of the fuzzy kernel is relatively low, and the heavy tail of the natural image gradient can be well satisfied. Distribution, reducing the ringing effect and quickly recovering high quality images.
  • super-resolution reconstruction technology to process the clear image can restore more image details of real objects, help to improve image quality, and help image recognition and image data acquisition and analysis after a series of image restoration operations.
  • Embodiment 1 is a flowchart of a method for comprehensively processing a blurred image according to Embodiment 1 of the present invention
  • FIG. 2 is a flowchart of a method for comprehensively processing a blurred image according to Embodiment 2 of the present invention
  • FIG. 3 is a block diagram of a fuzzy image integrated processing system according to Embodiment 3 of the present invention.
  • FIG. 4 is an overall block diagram of a defocus blur unit according to Embodiment 3 of the present invention.
  • FIG. 5 is a front and rear image display of a defocused image using the present invention, wherein (a) is a blurred image, and (b) is a restored image;
  • Figure 6 is a front and rear image display of a motion blurred image using the present invention, wherein (a) is a blurred image and (b) is a restored image.
  • Image restoration is the process of restoring a clear, high-quality original image from a blurred, noisy, low-quality, poorly-resolution, degraded image.
  • Image restoration firstly analyzes the cause of image distortion according to the image distortion phenomenon.
  • the fuzzy type of fuzzy image is generally divided into two categories: defocus blur and motion blur. Among them, defocus blur is caused by the fact that the image plane is not on the corresponding focal plane; the motion blur is caused by the relative motion between the imaging system and the target. Understand the causes of image distortion and then create different distortion models for different distortion causes, also known as degenerate models, and finally invert them to restore the original clear image.
  • the invention distinguishes the fuzzy type of the blurred image according to the elongation of the blurred image spectrogram, and based on different blur type images, different distortion causes, and performs different algorithms for defocus blur and motion blur to establish respective degradation models, which are more targeted.
  • this embodiment provides a method for comprehensive processing of a blurred image, which includes the following steps:
  • S10 Identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
  • the blurred image does not necessarily have defocus blur or motion blur in the complete sense, or it may be a mixture of the two, accompanied by some other types of noise, so it is more convenient to extract from the spectrogram.
  • Related geometric features generally perform preprocessing such as smoothing, image enhancement, binarization, etc., and these preprocessing are only used to eliminate noise interference and do not affect the differentiation of true fuzzy image degradation types.
  • the fuzzy type of the blurred image is identified according to the elongation of the blurred image spectrogram, and the elongation is specifically defined as:
  • A is the area of the strip or circle in the spectrogram
  • W and L are the width and length of the smallest rectangle surrounding the target respectively.
  • the method of distinguishing fuzzy types by elongation is simple, the calculation is convenient and fast, and the accuracy of identification is high, which can lay a good foundation for the later use of fuzzy kernel estimation methods for fuzzy kernel estimation.
  • the edge kernel method is used to obtain the blur kernel; if the blurred image is motion blur, the adaptive preset algorithm obtains the blur kernel, and the preset algorithm includes at least one algorithm;
  • the estimation of the fuzzy kernel plays a vital role in the restoration process of the blurred image. If the accuracy of the acquired fuzzy kernel is high, the subsequent restoration can adopt a simpler processing method.
  • the edge is a kind of image feature, which is the most uncertain place in the image, and the most concentrated image information.
  • the edge is also an important basis for image segmentation, and also an important reference for texture analysis and image recognition.
  • the blur of the defocused image mainly refers to the blur of the edge.
  • the edge edge method is used to obtain the step edge first, and the step edge is derived to calculate the point spread function, so that the fuzzy kernel closer to the real data can be obtained.
  • the preset algorithm may be a fuzzy kernel based on two-dimensional discrete wavelet transform and cepstrum analysis. Estimation algorithm, maximum likelihood method, Bayesian estimation algorithm, regularization algorithm, maximum entropy method, partial differential equation based algorithm and many other fuzzy kernel estimation algorithms.
  • the super Laplacian prior deconvolution algorithm is used to obtain a clear image, that is, using super Laplac as the image first. Knowledge, modeling, and rapid restoration of images for high-quality, clear images.
  • Embodiment 2 As shown in FIG. 2, the difference from Embodiment 1 is that the present embodiment provides a method for comprehensively processing a blurred image, and the method further includes the following steps:
  • S40 processing the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image.
  • acquiring the fuzzy kernel by using the edge method includes the following steps:
  • S202 Acquire a point spread equation according to the optimal edge edge image, where the point spread equation is the blur kernel, and the point spread equation is a point spread function, referred to as PSF, which is a spatial function representation of the fuzzy kernel.
  • the step S201 includes the following steps:
  • the Canny edge detection is performed on the blurred image to obtain the edge detection image.
  • the reason why the Canny edge detection operator is selected is that the algorithm has a good suppression effect on noise, single line response, high positioning accuracy, and some parameters involved in the algorithm have Versatility, which can be used for post-image restoration, Canny edge detection operator can detect the general edges in the image.
  • the quasi-determined position of the edge plays a role in the effective estimation of the subsequent transfer function (MTF) and point spread function (PSF).
  • MTF transfer function
  • PSF point spread function
  • the edge detection image is subjected to Hough transform to acquire a step edge image.
  • the optimal edge edge image extraction based on the gradient criterion is provided.
  • a blade edge image of each edge is taken centering on a center point of each edge in the step edge image, and the size of the edge edge image is preferably formed by a predetermined size according to the size of the defocused image and the degree of blur.
  • edge-edge image size should be appropriate, and should include all the important information of the blurred image point spread function; therefore, it is necessary to consider the defocused image when determining the edge-edge image size.
  • the size and blur degree, the edge image block is too small may not contain sufficient information, too large may cause the selected edge to be too curved, the calculation deviation is large, and the calculation process is time consuming, therefore, in actual application, Different images should be properly sized to the edge of the image.
  • the gradient value r of the edge image is calculated one by one, and the optimal edge image is extracted based on the gradient value r.
  • the gradient value r of the edge image is calculated one by one, and an optimal edge image is extracted in a plurality of the edge images according to the magnitude of the gradient value r.
  • the range of the line is larger than the range of the number of lines in which it is located. If the range of the number of columns in the original edge image is smaller than the range of the number of lines in which it is located, the angle is rotated so that the range of the number of columns in which the edge point is located is larger than the range of the number of columns.
  • said calculating the gradient value r of said edge image and extracting the optimal edge image based on the gradient value r comprises the steps of:
  • the edge in the edge image is linearly fitted by the least squares method to obtain the edge line.
  • the edge in the selected edge image is likely not a straight line in the strict sense, or the edge edge point distribution model caused by the edge detection error is not a straight line, so the edge is assumed in this embodiment.
  • the point obeys the straight line model, and the edge point of the edge in the edge image is straight-line fitted by the least square method.
  • the expression of the fitted line is:
  • n is the number of edge points of the edge
  • x k is the number of edge points
  • y k represents the relative position of the edge line.
  • the absolute value of the difference between the average pixel values of the two sides of the new edge edge point is counted as the gradient value r of the edge edge image, and the edge edge image with the largest gradient value r is extracted as the optimal edge edge image.
  • the gradient value r is calculated as follows:
  • the gradient value r calculated by using the edge edge information clearly reflects the distribution of the gray value of the bright and dark areas on both sides of the edge.
  • the step S202 includes the following steps:
  • the optimal edge line here is consistent with the calculation formula of the edge line in step S201. If there is a calculation before the edge line is saved, it can be read directly for calculation.
  • the vertical distance d of each pixel in the optimal edge image to the optimal edge line is obtained as the abscissa, and the gray value of each pixel is the ordinate. Scatter plot.
  • the value is used as the ordinate, so that the ESF scatter plot can be formed smoothly.
  • the distance from the point to the line is calculated as:
  • the scatter plot is linearly fitted using the Fermi function to obtain an edge spread function.
  • the edge spread function ESF can be fitted.
  • the present embodiment selects to be efficient and robust to random noise.
  • the improved Fermi function to fit the ESF the expression is as follows:
  • a point spread equation is calculated using the edge spread function, which is also a point spread function.
  • the line spread function LSF in the y direction can also be obtained in the above manner.
  • the MTF After obtaining the LSF in the x and y directions, the MTF can be obtained by Fourier transform, and then the point spread equation is obtained after convolving the MTF.
  • the PSF model is considered to be isotropic, so the two-dimensional PSF is separable, so the PSF can also be quickly calculated by:
  • the adaptive edge-edge method based on Gaussian model is used to estimate the fuzzy kernel, and the research shows that the degradation type of defocus blur is also coincident with the Gaussian model. Therefore, the method is active compared with the traditional edge-edge method.
  • the optimal edge-edge image is selected as the basic parameter of fuzzy kernel estimation, and the estimated fuzzy kernel energy is closer to the real fuzzy kernel.
  • the reason for the motion blur is mainly caused by the relative displacement change of the pixel, and in actual use, when the motion blur is generated, the relative motion between the target scene and the imaging device is not a uniform motion, so preferably, the In step S20:
  • the adaptive sparse prior method is used to estimate the fuzzy kernel, and the sparsity of the image gradient domain is used as the regular constraint. Firstly, the image pyramid of the blurred image is established, and then the fuzzy kernel and the clear image optimal value of each layer of image are calculated layer by layer by the method of alternating iteration until the last layer calculates the best fuzzy kernel as the parameter used for image restoration.
  • the filtered image is reprocessed by the impulse filter, which can suppress the noise well and strengthen the edge information of the image;
  • IRLS unconstrained iterative Re-weighted Least Squares
  • the embodiment provides a fuzzy image synthesis processing system, which is used to provide a physical implementation basis of the method in Embodiment 2, including:
  • the fuzzy type identification module 100 is configured to identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
  • the estimation module 200 is configured to obtain a fuzzy kernel by using a blade edge method if the blurred image is out-of-focus blur, and obtain a blur kernel by using an adaptive preset algorithm if the blurred image is motion blur, and the preset algorithm includes at least one Algorithm
  • the image restoration module 300 is configured to obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
  • system further includes:
  • the high-resolution processing module 400 processes the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image
  • the estimating module 200 includes:
  • the defocusing unit 210 is configured to extract an edge edge image based on a gradient criterion and obtain a point spread equation according to the optimal edge edge image if the blurred image is out of focus, and the point diffusion equation That is, the fuzzy kernel;
  • the motion blur unit 220 is configured to obtain a fuzzy kernel based on a sparse priori regularization method if the blurred image is motion blur.
  • the defocus blur unit 210 includes:
  • the edge detection sub-unit 211 is configured to perform Canny edge detection on the blurred image to obtain an edge detection image.
  • a line detection sub-unit 212 configured to perform a Hough transform on the edge detection image to acquire a step edge image
  • An image intercepting sub-unit 213, configured to intercept a blade edge image of each edge centering on a center point of each edge of the step edge image, the size of the edge edge image being according to the size of the defocused image And the degree of blurring is preferably formed by a predetermined size;
  • An image extraction sub-unit 214 configured to calculate a gradient value r of the edge-edge image one by one, and extract an optimal edge-edge image according to the gradient value r;
  • a linear fitting sub-unit 215, configured to perform linear fitting on the edge in the optimal edge-edge image by using a least square method to obtain an optimal edge line;
  • a scatter plot sub-unit 216 configured to obtain a vertical distance d of each pixel point in the optimal edge-edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate, Forming a scatter plot;
  • An edge diffusion function sub-unit 217 configured to perform a linear fit on the scatter plot by using a Fermi function to obtain an edge spread function
  • a calculation subunit 218 is configured to calculate a point spread equation using the edge spread function.
  • the fuzzy type recognition module 100 firstly identifies the blur type of the blurred image according to the elongation of the blurred image spectrogram; and then estimates the blur type to defocus blur by the estimation module 200 using different estimation methods. Or a blurred kernel of the blurred image of the motion blur; and then the image restoration module 300 obtains a clear image according to the restoration of the fuzzy kernel, and finally the clear image is processed by the high-resolution processing module 400 to obtain richer details and higher resolution, which is more advantageous. Identify and analyze images.
  • the estimation module 200 provided by the embodiment has a more accurate and fast fuzzy kernel calculation process for the blur type image with defocus blur, and can reduce human interference, avoid blindly extracting relevant parameters, and the possibility of reducing the reliability of the fuzzy kernel estimation.
  • the evaluation accuracy and speed of the image evaluation factor modulation transfer function (MTF) can be improved, which is more suitable for the restoration of real blurred images than the prior art.
  • FIG. 5 and FIG. 6 are front and rear comparison images of the blurred image after the above processing, and FIG. 5 shows the integrated processing of the defocused image, wherein (a) is a blurred image, and (b) is a restored image; 6 shows the image processing of motion blur, in which (a) is a blurred image and (b) is a restored image.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

Disclosed in the present invention are a comprehensive processing method and system for a blurred image. The method comprises: S10, identifying the blur type of a blurred image according to the elongation degree of a blurred image spectrogram, the blur type comprising defocus blur and motion blur; S20, if the blur type of the blurred image is defocus blur, obtaining a blur kernel by using a knife-edge method, and if the blur type of the blurred image is motion blur, obtaining the blur kernel by an adaptive preset algorithm; and S30, obtaining a clear image by a hyper-Laplace prior deconvolution algorithm on the basis of the blur kernel. According to the present invention, the blur types of blurred images are distinguished, and corresponding blur kernel estimation algorithms are set according to the difference between degradation mechanisms of different blur types, so that estimated blur kernels are more in line with actual situations; and especially for images of defocus blur, accurate blur kernel estimation based on a Gaussian model is performed. The present invention is suitable for restoration of more real blurred images, with the restoration method being simple and more targeted and the restoration being accurate, quick, and high-quality.

Description

一种模糊图像综合处理方法和系统Fuzzy image comprehensive processing method and system
本申请要求2018年01月31日提交的申请号为:201810094023.1、发明名称为“一种模糊图像综合处理方法和系统”的中国专利申请的优先权,其全部内容合并在此。The present application claims the priority of the Chinese Patent Application No. 201810094023.1, the entire disclosure of which is incorporated herein by reference.
技术领域Technical field
本发明涉及图像处理技术领域,尤其涉及一种模糊图像综合处理方法和系统。The present invention relates to the field of image processing technologies, and in particular, to a fuzzy image integration processing method and system.
背景技术Background technique
图像获取时会由于对焦不准,或者目标与相机产生相对运动、相机畸变、空气衍射等一些不确定的因素,导致拍摄的图像不尽如人意,丢失了图像中有用信息,也就是图像退化现象。阻碍了后续的处理分析,尤其是面对一些不可复制的场景,故对于此类情况,恢复出场景的原貌是一件迫切需要解决的问题。尽管现实设计中可以增设光学器件、改善传输设备等硬件手段如电子式像移补偿、光学式像移补偿、机械式像移补偿的方法来提高拍摄图像的质量,但因其器件工艺太复杂、成本造价较高、装调时间较长等诸多因素的影响,不具有普遍性,因此,利用图像复原技术,从图像补偿方面来提高图像的质量当是一个正确、低成本、低耗时的选择。Image acquisition may result in inaccurate focus, or relative motion of the camera and the camera, camera distortion, air diffraction, etc., resulting in unsatisfactory images, losing useful information in the image, that is, image degradation. . It hinders the subsequent processing analysis, especially in the face of some non-replicable scenes. Therefore, for such cases, restoring the original appearance of the scene is an urgent problem to be solved. Although the actual design can add optical devices, improve transmission equipment and other hardware means such as electronic image shift compensation, optical image shift compensation, mechanical image shift compensation to improve the quality of the captured image, but because the device process is too complicated, The impact of many factors such as high cost and long adjustment time is not universal. Therefore, using image restoration technology to improve image quality from image compensation is a correct, low-cost, low-cost option. .
图像复原需要了解图像退化的类型、机理和过程,只有知道了上述这些先验知识,就能够有针对性的确定准确的点扩散函数(PSF),也称模糊核,从而建立起退化模型,有了退化模型就可以通过算法将退化过程逆转过去,实现对模糊图像的复原。Image restoration needs to understand the type, mechanism and process of image degradation. Only by knowing these a priori knowledge, it is possible to determine the exact point spread function (PSF), also called fuzzy kernel, to establish a degradation model. The degradation model can be used to reverse the degradation process by algorithm to achieve the restoration of the blurred image.
图像复原时,如果没有合理的方法帮助了解图像退化的原因和模糊类型,就不能选取出最适合的复原算法帮助图像复原,点扩散函数估计不准确,选取的复原算法不恰当,都会造成低质量复原,更甚者不仅不能得到较好了模糊图像复原效果,还会使模糊图像更加模糊,因此除了准确估计点扩散函数十分重要外,还需要有针对性的获取准确的先验知识,即图像退化原因和模糊类型等。When the image is restored, if there is no reasonable way to understand the cause of the image degradation and the type of the blur, the most suitable restoration algorithm can not be selected to help the image restoration. The point spread function estimation is not accurate, and the selected restoration algorithm is not appropriate, which will result in low quality. Restoration, not only can not get better fuzzy image restoration effect, but also make the blurred image more fuzzy, so in addition to the accurate estimation of the point spread function is important, it also needs to obtain accurate prior knowledge, that is, the image Causes of degradation and types of blurring, etc.
图像复原算法在数学上被定义为一类不适定的反问题。已知原因求结果, 是正问题,已知结果反推原因是反问题;图像复原算法就是已知结果反推原因的算法。不适定的反问题是指,反推过程非常不稳定,即受到轻微的噪声影响,对最终的推测会造成非常大的干扰,导致结果错误。因此图像复原时尽量减少干扰对反推过程的影响。Image restoration algorithms are mathematically defined as a class of ill-posed inverse problems. Knowing the cause, the result is a positive problem, and the known result is the inverse problem; the image restoration algorithm is an algorithm that knows the reason for the reverse of the result. The ill-posed inverse problem means that the reverse push process is very unstable, that is, it is affected by slight noise, which will cause very large interference to the final speculation, resulting in incorrect results. Therefore, the image is restored to minimize the impact of interference on the reverse process.
现有技术中常用的图像复原算法有RL滤波、约束最小二乘方滤波、维纳滤波、正则滤波,常见的模糊图像类型有运动模糊和离焦模糊,复原清晰图像时,通过区分图像模糊类型,选择有针对性的图像复原方法进行复原不失为一个较佳的方法。区分的图像模糊类型的标准是否准确,图像复原方法对实际的图像复原是否恰当都是复原过程中需要反复考量和设计的。The image restoration algorithms commonly used in the prior art include RL filtering, constrained least squares filtering, Wiener filtering, and regular filtering. Common fuzzy image types include motion blur and defocus blur. When recovering clear images, by distinguishing image blur types It is a better method to select a targeted image restoration method for recovery. Whether the criteria for distinguishing image blur types are accurate, and whether the image restoration method is appropriate for the actual image restoration is repeated consideration and design during the restoration process.
如公开号为CN104331871A的专利文献公开了“一种图像去模糊方法及装置”,包括:对待处理图像进行模糊区域检测,确定模糊区域图像,判断所述模糊区域的模糊类型,若所述模糊区域图像的模糊类型为离焦模糊,则利用基于微分图像自相关的离焦模糊参数估计算法来确定离焦半径,若所述模糊区域图像的模糊类型为运动模糊,则基于倒谱分析的运动模糊参数估计算法来确定模糊方向和模糊尺度,将估计的参数代入经典图像复原算法中,得到复原图像。The patent document disclosed in the publication No. CN104331871A discloses "an image deblurring method and apparatus", comprising: performing blur region detection on a to-be-processed image, determining a blur region image, and determining a blur type of the blur region, if the blur region If the blur type of the image is defocus blur, the defocus blur parameter estimation algorithm based on differential image autocorrelation is used to determine the defocus radius. If the blur type of the blur region image is motion blur, the motion blur based on cepstrum analysis The parameter estimation algorithm determines the blur direction and the blur scale, and substitutes the estimated parameters into the classical image restoration algorithm to obtain the restored image.
该专利根据模糊类型不同选用不同的方法来进行复原时所必须的参数的估计,但是基于圆盘模型估计出离焦半径时,圆盘模型局限性较大,对于真实拍摄的离焦模糊图像,很难复原出清晰的图像;另现实生活当中,很多复杂的运动是无法确定方向及尺度的,即不是成线性运动的,故采用基于倒谱分析的运动模糊参数估计算法来确定模糊方向和模糊尺度并不能为图像复原提供完善的符合实际的先验知识;经典图像复原算法包括维纳滤波和LR滤波,维纳滤波在模糊图像峰值信噪比很小时,复原效果并不理想;LR滤波则是对噪声较为敏感,复原的图像具有明显的振铃效应,即适应范围不够广。The patent uses different methods to estimate the parameters necessary for restoration according to different types of blur. However, when the defocus radius is estimated based on the disc model, the disc model has a large limitation. For the defocused image of the real shot, It is difficult to recover a clear image; in real life, many complex motions cannot determine the direction and scale, that is, they are not linear motion, so the motion blur parameter estimation algorithm based on cepstrum analysis is used to determine the blur direction and blur. The scale does not provide perfect and realistic prior knowledge for image restoration; the classic image restoration algorithm includes Wiener filtering and LR filtering. Wiener filtering is not ideal when the peak signal-to-noise ratio of the blurred image is small, and the LR filtering is not ideal. It is sensitive to noise, and the restored image has a significant ringing effect, that is, the adaptation range is not wide enough.
发明内容Summary of the invention
本发明要解决的技术问题是针对上述现有技术的不足,提供一种适应于更多真实模糊图像复原,复原方法简单,针对性强,复原精准、速度快、质量高的模糊图像综合处理方法和系统。The technical problem to be solved by the present invention is to provide a fuzzy image comprehensive processing method suitable for more real fuzzy image restoration, simple restoration method, strong pertinence, accurate restoration, high speed and high quality, in view of the above-mentioned deficiencies of the prior art. And system.
为了实现上述目的,本发明采用的技术方案为:In order to achieve the above object, the technical solution adopted by the present invention is:
一种模糊图像综合处理方法,包括以下步骤:A method for comprehensively processing blurred images, comprising the following steps:
S10:根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;S10: Identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
S20:若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;S20: if the blurred image is out-of-focus blur, the edge kernel method is used to obtain the blur kernel; if the blurred image is motion blur, the adaptive preset algorithm obtains the blur kernel, and the preset algorithm includes at least one algorithm;
S30:基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。S30: Obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
进一步地,还包括以下步骤:Further, the method further includes the following steps:
S40:利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像。S40: processing the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image.
进一步地,所述步骤S20中若模糊图像为离焦模糊,则采用刃边法获取模糊核包括以下步骤:Further, if the blurred image is defocused in the step S20, acquiring the fuzzy kernel by using the edge method includes the following steps:
S201:提取基于梯度准则的最优刃边图像;S201: extract an optimal edge edge image based on a gradient criterion;
S202:根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核。S202: Acquire a point spread equation according to the optimal edge edge image, where the point spread equation is the blur kernel.
进一步地,所述步骤S201包括以下步骤:Further, the step S201 includes the following steps:
对模糊图像进行Canny边缘检测获取边缘检测图像;Performing an edge detection image by performing Canny edge detection on the blurred image;
对所述边缘检测图像进行霍夫变换获取阶跃边缘图像;Performing a Hough transform on the edge detection image to obtain a step edge image;
以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据所述离焦模糊图像大小和模糊程度选择一预设尺寸形成;Cutting a blade edge image of each edge centering on a center point of each edge in the step edge image, the size of the edge edge image selecting a preset size according to the size and blur degree of the defocused image form;
逐一计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像。The gradient values of the edge image are calculated one by one, and the optimal edge image is extracted based on the gradient values.
进一步地,所述计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像包括以下步骤:Further, the calculating the gradient value of the edge image and extracting the optimal edge image according to the gradient value comprises the following steps:
利用最小二乘法对所述刃边图像中的边缘进行线性拟合获取刃边直线;Using a least squares method to linearly fit the edges in the edge image to obtain a edge line;
提取所述刃边图像中的边缘与所述刃边直线重叠形成的交叉点作为新刃边边缘点;Extracting an intersection formed by an edge of the edge image and a straight line of the edge as a new edge edge point;
统计所述新刃边边缘点两侧区域的平均像素值之差的绝对值作为所述刃边图像的梯度值,提取梯度值最大的所述刃边图像作为最优刃边图像。The absolute value of the difference between the average pixel values of the two sides of the new edge edge point is counted as the gradient value of the edge edge image, and the edge edge image having the largest gradient value is extracted as the optimal edge edge image.
进一步地,所述步骤S202包括以下步骤:Further, the step S202 includes the following steps:
利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;Using an least square method to linearly fit the edges in the optimal edge image to obtain an optimal edge line;
获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离作为横坐标,每个像素点的灰度值为纵坐标,形成散点图;Obtaining a vertical distance of each pixel point in the optimal edge edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate to form a scatter plot;
利用费米函数对所述散点图进行线性拟合获取边缘扩散函数;Using the Fermi function to linearly fit the scatter plot to obtain an edge spread function;
利用所述边缘扩散函数计算点扩散方程。The point spread equation is calculated using the edge spread function.
进一步地,所述步骤S20中:Further, in the step S20:
若模糊图像为运动模糊,则自适应基于稀疏先验的正则化方法估计模糊核。If the blurred image is motion blur, the adaptive regularization method based on sparse prior is used to estimate the fuzzy kernel.
一种模糊图像综合处理系统,包括:A fuzzy image integrated processing system, comprising:
模糊类型识别模块,用于根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;a fuzzy type identification module, configured to identify a blur type of the blurred image according to an elongation of the blurred image spectrogram, the blur type including defocus blur and motion blur;
估计模块,用于若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;An estimation module, configured to acquire a fuzzy kernel by using a blade edge method if the blurred image is out-of-focus blur; and obtain a fuzzy kernel by using an adaptive preset algorithm if the blurred image is motion blur, the preset algorithm includes at least one algorithm;
图像复原模块,用于基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。An image restoration module is configured to obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
进一步地,所述系统还包括:Further, the system further includes:
高分辨率处理模块,利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像;a high resolution processing module that processes the sharp image using a super-resolution reconstruction technique to obtain a high resolution image;
所述估计模块包括:The estimation module includes:
离焦模糊单元,用于若模糊图像为离焦模糊时,采用刃边法提取基于梯度准则的最优刃边图像并根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核;a defocusing unit for extracting an edge-edge image based on a gradient criterion and acquiring a point-diffusion equation according to the optimal edge-edge image if the blurred image is defocused, and the point-diffusion equation is The blur kernel;
运动模糊单元,用于若模糊图像为运动模糊时,自适应基于稀疏先验的正则化方法获取模糊核。A motion blur unit is used to acquire a fuzzy kernel based on a sparse priori regularization method if the blurred image is motion blur.
进一步地,所述离焦模糊单元包括:Further, the defocus blur unit includes:
边缘检测子单元,用于对模糊图像进行Canny边缘检测获取边缘检测图像;An edge detection subunit, configured to perform Canny edge detection on the blurred image to obtain an edge detection image;
直线检测子单元,用于对所述边缘检测图像进行霍夫变换获取阶跃边缘图像;a line detection subunit, configured to perform a Hough transform on the edge detection image to obtain a step edge image;
图像截取子单元,用于以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据所述离焦模糊图像大小和模糊程度选择一预设尺寸形成;An image intercepting subunit for intercepting a blade edge image of each edge centering on a center point of each edge of the step edge image, the size of the edge edge image being according to the defocused image size and The degree of blur selects a predetermined size to be formed;
图像提取子单元,用于逐一计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像;An image extraction subunit, configured to calculate a gradient value of the edge edge image one by one, and extract an optimal edge edge image according to the gradient value;
线性拟合子单元,用于利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;a linear fitting sub-unit for linearly fitting edges in the optimal edge-edge image by using a least squares method to obtain an optimal edge line;
散点图子单元,用于获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离作为横坐标,每个像素点的灰度值为纵坐标,形成散点图;a scatter plot sub-unit, configured to obtain a vertical distance of each pixel point in the optimal edge-edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate, forming a dispersion Dot map
边缘扩散函数子单元,用于利用费米函数对所述散点图进行线性拟合获取边缘扩散函数;An edge diffusion function sub-unit for linearly fitting the scatter plot with a Fermi function to obtain an edge spread function;
计算子单元,用于利用所述边缘扩散函数计算点扩散方程。A calculation subunit is configured to calculate a point spread equation using the edge spread function.
采用上述技术方案后,本发明的有益效果是:After adopting the above technical solution, the beneficial effects of the present invention are:
根据伸长度确定模糊图像的类型不仅方法简单,适用于实际真实场景中的模糊图像类型清晰界定,能够保证模糊类型鉴别的精度,有助于有针对性地进行不同退化机理的模糊核的准确估计。Determining the type of fuzzy image according to the elongation is not only simple, but also suitable for the clear definition of the fuzzy image type in the actual real scene, which can ensure the accuracy of the fuzzy type identification and help to accurately estimate the fuzzy kernel with different degradation mechanisms. .
采用刃边法自适应地提取离焦模糊图像的最优刃边图像,且最优刃边图像是基于梯度准则获取的,一方面有利于提高离焦模糊图像复原的速度,另一方面有利于避免盲目提取刃边图像,减少人为干扰,提高图像复原的精度。The edge edge method is used to adaptively extract the optimal edge image of the defocused image, and the optimal edge image is obtained based on the gradient criterion. On the one hand, it is beneficial to improve the speed of defocus image restoration, on the other hand, it is beneficial to Avoid blindly extracting the edge image, reduce human interference, and improve the accuracy of image restoration.
依据最优刃边图像估计出的模糊核也更加接近于真实的退化模型,同时复原速度更快精度更高。The fuzzy kernel estimated from the optimal edge-edge image is also closer to the true degradation model, and the recovery speed is faster and more accurate.
自适应基于稀疏先验的正则化方法估计运动模糊图像的模糊核,对目标进行多次迭代直至收敛,得到更加稳定和理想的结果,从而解决了图像复原问题中的不适定的问题。The adaptive sparse priori-based regularization method estimates the fuzzy kernel of the motion blurred image, and iterates the target multiple times until convergence, which results in more stable and ideal results, thus solving the ill-posed problem in the image restoration problem.
采用超拉普拉斯先验去卷积算法获得清晰图像,以超拉普拉斯先验作为正则项,对模糊核精确性的要求相对较低,能够较好的满足自然图像梯度的重尾分布,减少振铃效应,快速的恢复出高质量的图像。The super Laplacian prior deconvolution algorithm is used to obtain clear images. The super Laplace a priori is used as a regular term. The accuracy of the fuzzy kernel is relatively low, and the heavy tail of the natural image gradient can be well satisfied. Distribution, reducing the ringing effect and quickly recovering high quality images.
采用超分辨率重建技术处理所述清晰图像,能够还原真实物体的更多图像细节,有助于提高图像质量,帮助图像识别和图像数据采集分析等一系列图像复原后的操作。The use of super-resolution reconstruction technology to process the clear image can restore more image details of real objects, help to improve image quality, and help image recognition and image data acquisition and analysis after a series of image restoration operations.
附图说明DRAWINGS
为了更清楚地说明本发明实施例或现有技术的技术方案,附图如下:In order to more clearly illustrate the embodiments of the present invention or the prior art, the drawings are as follows:
图1为本发明实施例1提供的一种模糊图像综合处理方法流程图;1 is a flowchart of a method for comprehensively processing a blurred image according to Embodiment 1 of the present invention;
图2为本发明实施例2提供的一种模糊图像综合处理方法流程图;2 is a flowchart of a method for comprehensively processing a blurred image according to Embodiment 2 of the present invention;
图3为本发明实施例3提供的一种模糊图像综合处理系统框图;3 is a block diagram of a fuzzy image integrated processing system according to Embodiment 3 of the present invention;
图4为本发明实施例3提供的离焦模糊单元整体框图;4 is an overall block diagram of a defocus blur unit according to Embodiment 3 of the present invention;
图5为离焦模糊图像利用本发明处理前后图像展示,其中(a)为模糊图像,(b)为复原图像;5 is a front and rear image display of a defocused image using the present invention, wherein (a) is a blurred image, and (b) is a restored image;
图6为运动模糊图像利用本发明处理前后图像展示,其中(a)为模糊图像,(b)为复原图像。Figure 6 is a front and rear image display of a motion blurred image using the present invention, wherein (a) is a blurred image and (b) is a restored image.
具体实施方式Detailed ways
以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。The technical solutions of the present invention are further described below with reference to the accompanying drawings, but the present invention is not limited to the embodiments.
图像复原是从模糊的、含噪声的、质量低的、分辨率差的退化图像中复原出清晰的、高质量的原始图像的过程。图像复原首先要根据图像失真现象,分析引起图像失真的原因,如模糊图像的模糊类型一般分为离焦模糊和运动模糊两大类。其中,离焦模糊是因为像面没有在对应的焦平面上而引起的一种模糊;运动模糊是由于成像系统和目标物间具有相对运动而引起的模糊。了解引起图像失真的原因后接着针对不同的失真原因建立不同失真模型,也称退化模型,最后对其求逆恢复出原来的清晰图像。Image restoration is the process of restoring a clear, high-quality original image from a blurred, noisy, low-quality, poorly-resolution, degraded image. Image restoration firstly analyzes the cause of image distortion according to the image distortion phenomenon. For example, the fuzzy type of fuzzy image is generally divided into two categories: defocus blur and motion blur. Among them, defocus blur is caused by the fact that the image plane is not on the corresponding focal plane; the motion blur is caused by the relative motion between the imaging system and the target. Understand the causes of image distortion and then create different distortion models for different distortion causes, also known as degenerate models, and finally invert them to restore the original clear image.
本发明根据模糊图像频谱图的伸长度区分模糊图像的模糊类型,基于不同模糊类型图像产生失真原因的不同,对离焦模糊和运动模糊执行不同的算法建立各自的退化模型,更有针对性的、准确的估计模糊核,并采用振铃效应少、计算速度快的超拉普拉斯先验去卷积算法复原图像。The invention distinguishes the fuzzy type of the blurred image according to the elongation of the blurred image spectrogram, and based on different blur type images, different distortion causes, and performs different algorithms for defocus blur and motion blur to establish respective degradation models, which are more targeted. Accurately estimate the fuzzy kernel, and use the super Laplacian prior deconvolution algorithm with less ringing effect and fast calculation speed to recover the image.
实施例1Example 1
如图1所示,本实施例提供一种模糊图像综合处理方法,包括以下步骤:As shown in FIG. 1 , this embodiment provides a method for comprehensive processing of a blurred image, which includes the following steps:
S10:根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;S10: Identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
尽管离焦模糊和运动模糊失真原因不同,但光从肉眼上很难准确分辨出这两种模糊类型,所以区分离焦模糊和运动模糊还需要寻找其他突破点,经过研究发现,离焦模糊和运动模糊两者的模糊图像的频谱图存在较大的差异,从频谱图中提取出相关的几何特征,设置阈值,就能够区分出模糊类型,方法简便容易实现;需要说明的是,此处区分模糊类型不是单纯为了区分模糊退化类型,而是希望依据类型归纳区分,建立更吻合其退化过程的退化模型,估计出更符合实际的模糊核。Although the causes of defocus blur and motion blur distortion are different, it is difficult for the light to accurately distinguish the two types of blur from the naked eye. Therefore, it is necessary to find other breakthrough points for the separation of focal blur and motion blur. After research, it is found that defocus blur and There is a big difference in the spectrogram of the blurred image of the motion blur. The relevant geometric features are extracted from the spectrogram, and the threshold is set to distinguish the fuzzy type. The method is simple and easy to implement. It should be noted that the distinction is made here. The fuzzy type is not simply to distinguish the type of fuzzy degradation, but it is hoped to differentiate according to the type, to establish a degradation model that is more consistent with its degradation process, and to estimate a more realistic fuzzy kernel.
实际场景中,模糊图像不一定是完全意义上只存在离焦模糊或者运动模糊的,也可能是这两者的混合,同时伴有些其他类型的噪声,故为了能从频谱图更方便的提取出相关的几何特征,一般会对频谱图进行平滑滤波、图像增强、二值化等预处理,且这些预处理仅用于排除噪声干扰,不影响真实模糊图像退化类型的区分。本实施例根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,伸长度具体定义为:In the actual scene, the blurred image does not necessarily have defocus blur or motion blur in the complete sense, or it may be a mixture of the two, accompanied by some other types of noise, so it is more convenient to extract from the spectrogram. Related geometric features generally perform preprocessing such as smoothing, image enhancement, binarization, etc., and these preprocessing are only used to eliminate noise interference and do not affect the differentiation of true fuzzy image degradation types. In this embodiment, the fuzzy type of the blurred image is identified according to the elongation of the blurred image spectrogram, and the elongation is specifically defined as:
T=W*L/AT=W*L/A
其中,A是频谱图中条状或圆形的面积,W、L分别是包围目标的最小矩形的宽度和长度,根据上述的计算公式,可见伸长度能够容易的区分出圆形和条状目标,即越接近于圆形,伸长度的值越小。而实验研究表明运动模糊的频谱图呈条状,离焦模糊的频谱图呈圆形,根据伸长度的计算公式可得:离焦模糊的伸长度小于运动模糊的伸长度,设定一个阈值,当伸长度大于阈值则表示该图像的模糊类型为运动模糊,反之则是离焦模糊。Where A is the area of the strip or circle in the spectrogram, and W and L are the width and length of the smallest rectangle surrounding the target respectively. According to the above formula, it can be seen that the elongation can easily distinguish the circular and strip targets. That is, the closer to the circle, the smaller the value of the elongation. The experimental research shows that the spectrogram of motion blur is strip-shaped, and the speckle spectrum is circular. According to the formula of elongation, the elongation of defocus blur is smaller than the elongation of motion blur, and a threshold is set. When the elongation is greater than the threshold, the blur type of the image is motion blur, and vice versa.
通过伸长度来区分模糊类型的方法简单,计算方便快速,且鉴别的精度高,能够为后期有针对性地利用模糊核估计方法进行模糊核估计打下良好的基础。The method of distinguishing fuzzy types by elongation is simple, the calculation is convenient and fast, and the accuracy of identification is high, which can lay a good foundation for the later use of fuzzy kernel estimation methods for fuzzy kernel estimation.
S20:若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;S20: if the blurred image is out-of-focus blur, the edge kernel method is used to obtain the blur kernel; if the blurred image is motion blur, the adaptive preset algorithm obtains the blur kernel, and the preset algorithm includes at least one algorithm;
模糊核的估计在模糊图像复原过程中,起着至关重要的作用,如果获取的模糊核准确度高,后续复原就能够采用更简单的处理方式。The estimation of the fuzzy kernel plays a vital role in the restoration process of the blurred image. If the accuracy of the acquired fuzzy kernel is high, the subsequent restoration can adopt a simpler processing method.
边缘是图像特征的一种,是图像中不确定性最大的地方,也是图像信息最集中的地方,边缘也是图像分割的重要基础,也是纹理分析和图像识别的重要参考。离焦模糊图像的模糊主要指的就是边缘的模糊,采用刃边法先获取阶跃边缘,对阶跃边缘推导计算点扩散函数,就能够获取较为接近真实数据的模糊核。The edge is a kind of image feature, which is the most uncertain place in the image, and the most concentrated image information. The edge is also an important basis for image segmentation, and also an important reference for texture analysis and image recognition. The blur of the defocused image mainly refers to the blur of the edge. The edge edge method is used to obtain the step edge first, and the step edge is derived to calculate the point spread function, so that the fuzzy kernel closer to the real data can be obtained.
另外,本实施例中,若模糊图像为运动模糊时,则自适应选取一种预设定算法估计模糊核,所述预设定算法可以是基于二维离散小波变换和倒频谱分析的模糊核估计算法、最大似然法、贝叶斯估计算法、正则化算法、最大熵法、基于偏微分方程的算法等等众多模糊核估计算法中的任一一种或多种。In addition, in this embodiment, if the blurred image is motion blur, adaptively selecting a preset algorithm to estimate the fuzzy kernel, the preset algorithm may be a fuzzy kernel based on two-dimensional discrete wavelet transform and cepstrum analysis. Estimation algorithm, maximum likelihood method, Bayesian estimation algorithm, regularization algorithm, maximum entropy method, partial differential equation based algorithm and many other fuzzy kernel estimation algorithms.
S30:基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。S30: Obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
不管是离焦模糊的模糊图像还是运动模糊的模糊图像,计算出模糊核后,均利用超拉普拉斯先验去卷积算法获得清晰图像,即采用将超拉普拉斯作为图像的先验知识,进行建模,快速的复原图像获得高质量清晰图像。Regardless of whether it is a defocused blur image or a motion blur blur image, after calculating the blur kernel, the super Laplacian prior deconvolution algorithm is used to obtain a clear image, that is, using super Laplac as the image first. Knowledge, modeling, and rapid restoration of images for high-quality, clear images.
实施例2Example 2
如图2所示,与实施例1的区别在于,本实施例提供一种模糊图像综合处理方法,所述方法还包括以下步骤:As shown in FIG. 2, the difference from Embodiment 1 is that the present embodiment provides a method for comprehensively processing a blurred image, and the method further includes the following steps:
S40:利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像。S40: processing the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image.
具体地,所述步骤S20中若模糊图像为离焦模糊,则采用刃边法获取模糊核包括以下步骤:Specifically, if the blurred image is defocused in the step S20, acquiring the fuzzy kernel by using the edge method includes the following steps:
S201:提取基于梯度准则的最优刃边图像;S201: extract an optimal edge edge image based on a gradient criterion;
S202:根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核,点扩散方程就是点扩散函数,简称PSF,是模糊核的空间函数表示。S202: Acquire a point spread equation according to the optimal edge edge image, where the point spread equation is the blur kernel, and the point spread equation is a point spread function, referred to as PSF, which is a spatial function representation of the fuzzy kernel.
所述步骤S201包括以下步骤:The step S201 includes the following steps:
对模糊图像进行Canny边缘检测获取边缘检测图像,之所以选择Canny边缘检测算子其原因在于该算法对噪声有很好的抑制作用,单线响应,定位精度高,且算法里涉及到的一些参数具有通用性,能够为后期图像复原所利用,Canny边缘检测算子对于图像中一般的边缘都可以检测出来。边缘的准确定 位,对于后面的传递函数(MTF)及点扩散函数(PSF)的有效估计起到锦上添花的作用。The Canny edge detection is performed on the blurred image to obtain the edge detection image. The reason why the Canny edge detection operator is selected is that the algorithm has a good suppression effect on noise, single line response, high positioning accuracy, and some parameters involved in the algorithm have Versatility, which can be used for post-image restoration, Canny edge detection operator can detect the general edges in the image. The quasi-determined position of the edge plays a role in the effective estimation of the subsequent transfer function (MTF) and point spread function (PSF).
对所述边缘检测图像进行霍夫变换获取阶跃边缘图像。The edge detection image is subjected to Hough transform to acquire a step edge image.
对边缘检测后的边缘检测图像执行一次霍夫(Hough)变换,此时会出现一些很明显的Hough变换峰值点,根据这些峰值点,即可寻找并链接到峰值点所对应的直线,且可以标记出这些直线在图像中的具体空间位置,这些直线也就是阶跃边缘,常规阶跃边缘图像中会存在多个阶跃边缘,为了避免人为盲目随意的选取一阶跃边缘作为模糊核估计的依据,本实施例中提供基于梯度准则的最优刃边图像提取。Perform a Hough transform on the edge detection image after edge detection. At this time, some obvious Hough transform peak points appear. According to these peak points, you can find and link to the line corresponding to the peak point, and Mark the specific spatial positions of these lines in the image. These lines are also step edges. There are multiple step edges in the regular step edge image. In order to avoid artificially blindly selecting a step edge as the fuzzy kernel estimation. According to the embodiment, the optimal edge edge image extraction based on the gradient criterion is provided.
以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据离焦模糊图像大小和模糊程度优选一预设尺寸形成。A blade edge image of each edge is taken centering on a center point of each edge in the step edge image, and the size of the edge edge image is preferably formed by a predetermined size according to the size of the defocused image and the degree of blur.
具体地,确定阶跃边缘图像中每一边缘中心点坐标,并以其为中心点,截取一定尺寸大小的含有部分边缘的刃边图像,即阶跃边缘图像中有多少条独立的阶跃边缘,实际就有多少刃边图像,需要说明的是,刃边图像尺寸大小应合适,应尽量包含模糊图像点扩散函数的所有重要信息;故确定刃边图像大小的时候需要考虑离焦模糊图像的大小及模糊程度,刃边图像块太小可能包含不了充分的信息,太大可能造成选取的刃边太弯曲,计算偏差较大,同时计算的过程中较耗时,因此,实际运用中,针对不同的图像要做适当的刃边图像尺寸大小调整。Specifically, determining the coordinates of the center point of each edge in the step edge image, and taking the center point as a center point, intercepting a certain edge image with a partial edge, that is, how many independent step edges are in the step edge image Actually, there are many edge-edge images. It should be noted that the edge-edge image size should be appropriate, and should include all the important information of the blurred image point spread function; therefore, it is necessary to consider the defocused image when determining the edge-edge image size. The size and blur degree, the edge image block is too small may not contain sufficient information, too large may cause the selected edge to be too curved, the calculation deviation is large, and the calculation process is time consuming, therefore, in actual application, Different images should be properly sized to the edge of the image.
逐一计算所述刃边图像的梯度值r,并根据梯度值r提取最优刃边图像。The gradient value r of the edge image is calculated one by one, and the optimal edge image is extracted based on the gradient value r.
具体地,逐一计算所述刃边图像的梯度值r,并根据梯度值r的大小在诸多所述刃边图像中提取最优刃边图像。需要说明的是,计算所述刃边图像的梯度值r时,为了减小梯度准则计算误差,提高最优刃边图像提取的准确率,需要确保刃边图像中边缘的边缘点所在的列数的范围大于其所在的行数范围,若原始刃边图像中边缘点所在的列数的范围小于其所在的行数范围,则旋转一定角度,使得边缘点所在的列数的范围大于其所在的行数范围后,作为新的刃边图像,再来执行梯度值r计算。Specifically, the gradient value r of the edge image is calculated one by one, and an optimal edge image is extracted in a plurality of the edge images according to the magnitude of the gradient value r. It should be noted that, when calculating the gradient value r of the edge edge image, in order to reduce the gradient criterion calculation error and improve the accuracy of the optimal edge edge image extraction, it is necessary to ensure the number of columns of the edge edge of the edge edge image. The range of the line is larger than the range of the number of lines in which it is located. If the range of the number of columns in the original edge image is smaller than the range of the number of lines in which it is located, the angle is rotated so that the range of the number of columns in which the edge point is located is larger than the range of the number of columns. After the number of rows, as a new edge image, the gradient value r calculation is performed.
作为优选地,所述计算所述刃边图像的梯度值r,并根据梯度值r提取最 优刃边图像包括以下步骤:Preferably, said calculating the gradient value r of said edge image and extracting the optimal edge image based on the gradient value r comprises the steps of:
首先利用最小二乘法对所述刃边图像中的边缘进行线性拟合获取刃边直线。First, the edge in the edge image is linearly fitted by the least squares method to obtain the edge line.
通常情况下,所选取的刃边图像中的边缘很有可能不是一条严格意义上的直线,或者是边缘检测误差等原因造成的刃边边缘点分布模型不是一条直线,故本实施例中假设边缘点服从直线模型,借助最小二乘法对所述刃边图像中边缘的边缘点进行直线拟合,拟合直线的表达式为:In general, the edge in the selected edge image is likely not a straight line in the strict sense, or the edge edge point distribution model caused by the edge detection error is not a straight line, so the edge is assumed in this embodiment. The point obeys the straight line model, and the edge point of the edge in the edge image is straight-line fitted by the least square method. The expression of the fitted line is:
y=ax+by=ax+b
Figure PCTCN2018091164-appb-000001
Figure PCTCN2018091164-appb-000001
Figure PCTCN2018091164-appb-000002
Figure PCTCN2018091164-appb-000002
其中,n是刃边边缘点数目,x k是边缘点列数,y k代表边缘点行的相对位置,先求得直线系数a和b,从而确定拟合的刃边直线。 Where n is the number of edge points of the edge, x k is the number of edge points, and y k represents the relative position of the edge line. First, the linear coefficients a and b are obtained to determine the fitted edge line.
其次,提取所述刃边图像中的边缘与所述刃边直线重叠形成的交叉点作为新刃边边缘点。Next, an intersection formed by overlapping the edge in the edge image with the edge of the edge is extracted as a new edge edge point.
最后,统计所述新刃边边缘点两侧区域的平均像素值之差的绝对值作为所述刃边图像的梯度值r,提取梯度值r最大的所述刃边图像作为最优刃边图像。Finally, the absolute value of the difference between the average pixel values of the two sides of the new edge edge point is counted as the gradient value r of the edge edge image, and the edge edge image with the largest gradient value r is extracted as the optimal edge edge image. .
具体地,所述梯度值r计算公式如下:Specifically, the gradient value r is calculated as follows:
Figure PCTCN2018091164-appb-000003
Figure PCTCN2018091164-appb-000003
Figure PCTCN2018091164-appb-000004
Figure PCTCN2018091164-appb-000004
r=|gc_1-gc_2|r=|gc_1-gc_2|
其中,m和n分别是刃边图像的长和宽,r j(j=1,2,...,n)代表每列上刃边边缘点所在的行数,e i,j代表位置(i,j)处的灰度值。 Where m and n are the length and width of the edge image, respectively, r j (j=1, 2, . . . , n) represents the number of lines in the edge of each edge of the edge, and e i,j represents the position ( Gray value at i, j).
根据上式利用刃边边缘点的信息计算出的梯度值r较清晰地反映了刃边两侧亮暗区域的灰度值的分布情况,r值越大,表示刃边两侧的灰度值对比度越大以及同一侧的灰度值相似度越高,越有利于模糊核的估计及后续的图像恢复。故提取梯度值r最大的所述刃边图像作为最优刃边图像。According to the above formula, the gradient value r calculated by using the edge edge information clearly reflects the distribution of the gray value of the bright and dark areas on both sides of the edge. The larger the r value, the gray value on both sides of the edge. The greater the contrast and the higher the similarity of the gray values on the same side, the more favorable the estimation of the fuzzy kernel and the subsequent image restoration. Therefore, the edge edge image having the largest gradient value r is extracted as an optimum edge edge image.
所述步骤S202包括以下步骤:The step S202 includes the following steps:
利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;此处的最优刃边直线与步骤S201中的刃边直线的计算公式相一致,若刃边直线之前有过计算就会进行保存,则此处可以直接读取用于计算即可。Using the least square method to linearly fit the edge in the optimal edge image to obtain the optimal edge line; the optimal edge line here is consistent with the calculation formula of the edge line in step S201. If there is a calculation before the edge line is saved, it can be read directly for calculation.
得到最优刃边直线后,获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离d作为横坐标,每个像素点的灰度值为纵坐标,形成散点图。After obtaining the optimal edge line, the vertical distance d of each pixel in the optimal edge image to the optimal edge line is obtained as the abscissa, and the gray value of each pixel is the ordinate. Scatter plot.
通常,边缘扩散函数(ESF)散点图的准确获取是拟合ESF曲线的必要过程,以最优刃边图像块中每个像素(i,j)(i=1,2,…,m;j=1,2,…,n)到拟合的最优刃边直线ax+by+c=0的垂直距离d作为横坐标(以像素为单位),每个像素点的归一化灰度值作为纵坐标,从而就可以顺利形成ESF散点图。其中,点到直线的距离计算公式为:In general, the accurate acquisition of the edge spread function (ESF) scatter plot is a necessary process for fitting the ESF curve to each pixel (i, j) in the optimal edge-edge image block (i = 1, 2, ..., m; j=1,2,...,n) to the fitted edge edge straight line ax+by+c=0 vertical distance d as the abscissa (in pixels), normalized grayscale of each pixel The value is used as the ordinate, so that the ESF scatter plot can be formed smoothly. Among them, the distance from the point to the line is calculated as:
Figure PCTCN2018091164-appb-000005
Figure PCTCN2018091164-appb-000005
利用费米函数对所述散点图进行线性拟合获取边缘扩散函数。The scatter plot is linearly fitted using the Fermi function to obtain an edge spread function.
形成距离-灰度的ESF散点图后,就可以拟合边缘扩散函数ESF,鉴于ESF散点图所反应的原始模糊图像的边缘信息的局限性,本实施例选择对随机噪声具有高效鲁棒性的改进后的费米(Fermi)函数来拟合ESF,其表达式如下:After forming the distance-gray ESF scatter plot, the edge spread function ESF can be fitted. In view of the limitation of the edge information of the original blurred image reflected by the ESF scatter plot, the present embodiment selects to be efficient and robust to random noise. The improved Fermi function to fit the ESF, the expression is as follows:
Figure PCTCN2018091164-appb-000006
Figure PCTCN2018091164-appb-000006
利用所述边缘扩散函数计算点扩散方程(PSF),所述点扩散方程也就是点扩散函数。A point spread equation (PSF) is calculated using the edge spread function, which is also a point spread function.
对边缘扩散函数进行求导,即可获得x方向上的线扩散函数LSF,Deriving the edge spread function to obtain the line spread function LSF in the x direction,
Figure PCTCN2018091164-appb-000007
Figure PCTCN2018091164-appb-000007
此外,y方向上的线扩散函数LSF也可通过上述方式获得。Further, the line spread function LSF in the y direction can also be obtained in the above manner.
在得到x和y方向上的LSF后,可通过傅里叶变换获得MTF,再对MTF卷积后获得点扩散方程。After obtaining the LSF in the x and y directions, the MTF can be obtained by Fourier transform, and then the point spread equation is obtained after convolving the MTF.
一般的,认为PSF模型为各向同性分布,故二维PSF具有可分离性,故也可以通过下式快速计算PSF:In general, the PSF model is considered to be isotropic, so the two-dimensional PSF is separable, so the PSF can also be quickly calculated by:
PSF(x,y)=LSF(x)×LSF(y)PSF(x,y)=LSF(x)×LSF(y)
综上,本实施例采用了基于高斯模型的自适应刃边法估计模糊核,而且研究表明离焦模糊的退化类型也恰好与高斯模型相吻合,故较传统的刃边法,该方法通过主动选取较优的刃边图像,作为模糊核估计的基础参数,估计出的模糊核能更加接近于真实模糊核。In summary, the adaptive edge-edge method based on Gaussian model is used to estimate the fuzzy kernel, and the research shows that the degradation type of defocus blur is also coincident with the Gaussian model. Therefore, the method is active compared with the traditional edge-edge method. The optimal edge-edge image is selected as the basic parameter of fuzzy kernel estimation, and the estimated fuzzy kernel energy is closer to the real fuzzy kernel.
由于运动模糊产生的原因主要是像素点的相对位移变化造成的,且实际 运用中,运动模糊产生时目标场景和成像设备之间的相对运动不是匀速运动的情况较多,故优选地,所述步骤S20中:The reason for the motion blur is mainly caused by the relative displacement change of the pixel, and in actual use, when the motion blur is generated, the relative motion between the target scene and the imaging device is not a uniform motion, so preferably, the In step S20:
若模糊图像为运动模糊,则自适应基于稀疏先验的正则化方法估计模糊核,利用图像梯度域的稀疏性来作为正则约束项。首先建立模糊图像的图像金字塔,然后利用交替迭代的方法逐层计算每层图像的模糊核以及清晰图像最优值,直至最后一层计算出最佳模糊核作为图像复原所使用的参数。If the blurred image is motion blur, the adaptive sparse prior method is used to estimate the fuzzy kernel, and the sparsity of the image gradient domain is used as the regular constraint. Firstly, the image pyramid of the blurred image is established, and then the fuzzy kernel and the clear image optimal value of each layer of image are calculated layer by layer by the method of alternating iteration until the last layer calculates the best fuzzy kernel as the parameter used for image restoration.
具体步骤如下:Specific steps are as follows:
(1)输入模糊图像g,模糊核尺寸m,其中模糊核尺寸为给定值;(1) input a blurred image g, a fuzzy kernel size m, wherein the fuzzy kernel size is a given value;
(2)根据模糊核尺寸确定分解层数,对每层图像进行双边滤波;(2) determining the number of decomposition layers according to the size of the fuzzy kernel, and performing bilateral filtering on each layer of the image;
(3)采用冲激滤波器对滤波后的图像进行再处理,双重滤波后能够很好的抑制噪声,并且强化图像的边缘信息;(3) The filtered image is reprocessed by the impulse filter, which can suppress the noise well and strengthen the edge information of the image;
(4)一阶求导得到的梯度图像y;(4) Gradient image y obtained by first order derivation;
(5)采用迭代收缩阈值算法(Iteration Shrinkage Thresholding Algorithm,ISTA)算法更新求解f子问题,其中f就是每层图像的清晰图像最优值;(5) Using the Iteration Shrinkage Thresholding Algorithm (ISTA) algorithm to update the solution to the f-sub-problem, where f is the clear image optimal value of each layer of image;
(6)采用无约束迭代重加权最小二乘法(Iterative Re-weighted Least Squares,IRLS)算法求解h子问题,其中h就是每层图像的模糊核;(6) The unconstrained iterative Re-weighted Least Squares (IRLS) algorithm is used to solve the h sub-problem, where h is the fuzzy kernel of each layer of image;
(7)将本层求得的模糊核h加入到下一层的迭代过程中直至最后一层求得较为稳定的模糊核,作为最佳模糊核。(7) Add the fuzzy kernel h obtained by this layer to the iterative process of the next layer until the last layer obtains a relatively stable fuzzy kernel as the best fuzzy kernel.
实施例3Example 3
如图3所示,本实施例提供一种模糊图像综合处理系统,用于提供实施例2中所述方法的物理实现基础,包括:As shown in FIG. 3, the embodiment provides a fuzzy image synthesis processing system, which is used to provide a physical implementation basis of the method in Embodiment 2, including:
模糊类型识别模块100,用于根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;The fuzzy type identification module 100 is configured to identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
估计模块200,用于若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;The estimation module 200 is configured to obtain a fuzzy kernel by using a blade edge method if the blurred image is out-of-focus blur, and obtain a blur kernel by using an adaptive preset algorithm if the blurred image is motion blur, and the preset algorithm includes at least one Algorithm
图像复原模块300,用于基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。The image restoration module 300 is configured to obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
可选地,所述系统还包括:Optionally, the system further includes:
高分辨率处理模块400,利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像;The high-resolution processing module 400 processes the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image;
所述估计模块200包括:The estimating module 200 includes:
离焦模糊单元210,用于若模糊图像为离焦模糊时,采用刃边法提取基于梯度准则的最优刃边图像并根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核;The defocusing unit 210 is configured to extract an edge edge image based on a gradient criterion and obtain a point spread equation according to the optimal edge edge image if the blurred image is out of focus, and the point diffusion equation That is, the fuzzy kernel;
运动模糊单元220,用于若模糊图像为运动模糊时,自适应基于稀疏先验的正则化方法获取模糊核。The motion blur unit 220 is configured to obtain a fuzzy kernel based on a sparse priori regularization method if the blurred image is motion blur.
可选地,如图4所示,所述离焦模糊单元210包括:Optionally, as shown in FIG. 4, the defocus blur unit 210 includes:
边缘检测子单元211,用于对模糊图像进行Canny边缘检测获取边缘检测图像;The edge detection sub-unit 211 is configured to perform Canny edge detection on the blurred image to obtain an edge detection image.
直线检测子单元212,用于对所述边缘检测图像进行霍夫变换获取阶跃边缘图像;a line detection sub-unit 212, configured to perform a Hough transform on the edge detection image to acquire a step edge image;
图像截取子单元213,用于以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据所述离焦模糊图像大小和模糊程度优选一预设尺寸形成;An image intercepting sub-unit 213, configured to intercept a blade edge image of each edge centering on a center point of each edge of the step edge image, the size of the edge edge image being according to the size of the defocused image And the degree of blurring is preferably formed by a predetermined size;
图像提取子单元214,用于逐一计算所述刃边图像的梯度值r,并根据梯度值r提取最优刃边图像;An image extraction sub-unit 214, configured to calculate a gradient value r of the edge-edge image one by one, and extract an optimal edge-edge image according to the gradient value r;
线性拟合子单元215,用于利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;a linear fitting sub-unit 215, configured to perform linear fitting on the edge in the optimal edge-edge image by using a least square method to obtain an optimal edge line;
散点图子单元216,用于获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离d作为横坐标,每个像素点的灰度值为纵坐标,形成散点图;a scatter plot sub-unit 216, configured to obtain a vertical distance d of each pixel point in the optimal edge-edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate, Forming a scatter plot;
边缘扩散函数子单元217,用于利用费米函数对所述散点图进行线性拟合获取边缘扩散函数;An edge diffusion function sub-unit 217, configured to perform a linear fit on the scatter plot by using a Fermi function to obtain an edge spread function;
计算子单元218,用于利用所述边缘扩散函数计算点扩散方程。A calculation subunit 218 is configured to calculate a point spread equation using the edge spread function.
本实施例提供的模糊图像综合处理系统,首先由模糊类型识别模块100根据模糊图像频谱图的伸长度识别模糊图像的模糊类型;然后通过估计模块 200采用不同的估计方法估计模糊类型为离焦模糊或运动模糊的模糊图像的模糊核;再接着由图像复原模块300,根据模糊核复原获得清晰图像,最后清晰图像经过高分辨率处理模块400处理得到细节更丰富、分辨率更高,更有利于识别和分析的图像。In the fuzzy image integration processing system provided by the embodiment, the fuzzy type recognition module 100 firstly identifies the blur type of the blurred image according to the elongation of the blurred image spectrogram; and then estimates the blur type to defocus blur by the estimation module 200 using different estimation methods. Or a blurred kernel of the blurred image of the motion blur; and then the image restoration module 300 obtains a clear image according to the restoration of the fuzzy kernel, and finally the clear image is processed by the high-resolution processing module 400 to obtain richer details and higher resolution, which is more advantageous. Identify and analyze images.
本实施例提供的估计模块200对模糊类型为离焦模糊的图像,具有更加精准快速的模糊核计算过程,且能够减少人为干扰,避免盲目提取相关参数,导致模糊核估算可靠性降低的可能性,同时基于推导计算过程参数的计算和提取,能够提高图像评价因子调制传递函数(MTF)的评价精度和速度,与现有技术相比,更适用于真实模糊图像的复原。具体地,图5、图6中为经过上述处理后模糊图像的前后对比图,图5显示了离焦模糊的图像综合处理情况,其中(a)为模糊图像,(b)为复原图像;图6显示了运动模糊的图像综合处理情况,其中(a)为模糊图像,(b)为复原图像。The estimation module 200 provided by the embodiment has a more accurate and fast fuzzy kernel calculation process for the blur type image with defocus blur, and can reduce human interference, avoid blindly extracting relevant parameters, and the possibility of reducing the reliability of the fuzzy kernel estimation. At the same time, based on the calculation and extraction of the derivation calculation process parameters, the evaluation accuracy and speed of the image evaluation factor modulation transfer function (MTF) can be improved, which is more suitable for the restoration of real blurred images than the prior art. Specifically, FIG. 5 and FIG. 6 are front and rear comparison images of the blurred image after the above processing, and FIG. 5 shows the integrated processing of the defocused image, wherein (a) is a blurred image, and (b) is a restored image; 6 shows the image processing of motion blur, in which (a) is a blurred image and (b) is a restored image.
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。The specific embodiments described herein are merely illustrative of the spirit of the invention. A person skilled in the art can make various modifications or additions to the specific embodiments described or in a similar manner, without departing from the spirit of the invention or as defined by the appended claims. The scope.

Claims (10)

  1. 一种模糊图像综合处理方法,其特征在于,包括以下步骤:A method for comprehensively processing a blurred image, comprising the steps of:
    S10:根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;S10: Identify a blur type of the blurred image according to the elongation of the blurred image spectrogram, where the blur type includes defocus blur and motion blur;
    S20:若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;S20: if the blurred image is out-of-focus blur, the edge kernel method is used to obtain the blur kernel; if the blurred image is motion blur, the adaptive preset algorithm obtains the blur kernel, and the preset algorithm includes at least one algorithm;
    S30:基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。S30: Obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
  2. 根据权利要求1所述的一种模糊图像综合处理方法,其特征在于,还包括以下步骤:A method for comprehensively processing blurred images according to claim 1, further comprising the steps of:
    S40:利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像。S40: processing the clear image by using a super-resolution reconstruction technique to obtain a high-resolution image.
  3. 根据权利要求1所述的一种模糊图像综合处理方法,其特征在于,所述步骤S20中若模糊图像为离焦模糊,则采用刃边法获取模糊核包括以下步骤:The method of processing a blurred image according to claim 1, wherein if the blurred image is defocused in the step S20, acquiring the fuzzy kernel by using the edge method comprises the following steps:
    S201:提取基于梯度准则的最优刃边图像;S201: extract an optimal edge edge image based on a gradient criterion;
    S202:根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核。S202: Acquire a point spread equation according to the optimal edge edge image, where the point spread equation is the blur kernel.
  4. 根据权利要求3所述的一种模糊图像综合处理方法,其特征在于,所述步骤S201包括以下步骤:The method of processing a blurred image according to claim 3, wherein the step S201 comprises the following steps:
    对模糊图像进行Canny边缘检测获取边缘检测图像;Performing an edge detection image by performing Canny edge detection on the blurred image;
    对所述边缘检测图像进行霍夫变换获取阶跃边缘图像;Performing a Hough transform on the edge detection image to obtain a step edge image;
    以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据所述离焦模糊图像大小和模糊程度选择一预设尺寸形成;Cutting a blade edge image of each edge centering on a center point of each edge in the step edge image, the size of the edge edge image selecting a preset size according to the size and blur degree of the defocused image form;
    逐一计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像。The gradient values of the edge image are calculated one by one, and the optimal edge image is extracted based on the gradient values.
  5. 根据权利要求4所述的一种模糊图像综合处理方法,其特征在于,所述计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像包括以下步骤:A fuzzy image synthesis processing method according to claim 4, wherein said calculating a gradient value of said edge image and extracting an optimal edge image based on the gradient value comprises the following steps:
    利用最小二乘法对所述刃边图像中的边缘进行线性拟合获取刃边直线;Using a least squares method to linearly fit the edges in the edge image to obtain a edge line;
    提取所述刃边图像中的边缘与所述刃边直线重叠形成的交叉点作为新刃边边缘点;Extracting an intersection formed by an edge of the edge image and a straight line of the edge as a new edge edge point;
    统计所述新刃边边缘点两侧区域的平均像素值之差的绝对值作为所述刃边图像的梯度值,提取梯度值最大的所述刃边图像作为最优刃边图像。The absolute value of the difference between the average pixel values of the two sides of the new edge edge point is counted as the gradient value of the edge edge image, and the edge edge image having the largest gradient value is extracted as the optimal edge edge image.
  6. 根据权利要求3所述的一种模糊图像综合处理方法,其特征在于,所述步骤S202包括以下步骤:The method of processing a blurred image according to claim 3, wherein the step S202 comprises the following steps:
    利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;Using an least square method to linearly fit the edges in the optimal edge image to obtain an optimal edge line;
    获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离作为横坐标,每个像素点的灰度值为纵坐标,形成散点图;Obtaining a vertical distance of each pixel point in the optimal edge edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate to form a scatter plot;
    利用费米函数对所述散点图进行线性拟合获取边缘扩散函数;Using the Fermi function to linearly fit the scatter plot to obtain an edge spread function;
    利用所述边缘扩散函数计算点扩散方程。The point spread equation is calculated using the edge spread function.
  7. 根据权利要求1所述的一种模糊图像综合处理方法,其特征在于,所述步骤S20中:A method for comprehensively processing blurred images according to claim 1, wherein in step S20:
    若模糊图像为运动模糊,则自适应基于稀疏先验的正则化方法获取模糊核。If the blurred image is motion blur, the adaptive sparse priori-based regularization method obtains the fuzzy kernel.
  8. 一种模糊图像综合处理系统,其特征在于,包括:A fuzzy image integrated processing system, comprising:
    模糊类型识别模块,用于根据模糊图像频谱图的伸长度识别模糊图像的模糊类型,所述模糊类型包括离焦模糊和运动模糊;a fuzzy type identification module, configured to identify a blur type of the blurred image according to an elongation of the blurred image spectrogram, the blur type including defocus blur and motion blur;
    估计模块,用于若模糊图像为离焦模糊,则采用刃边法获取模糊核;若模糊图像为运动模糊,则自适应预设定算法获取模糊核,所述预设定算法包括至少一种算法;An estimation module, configured to acquire a fuzzy kernel by using a blade edge method if the blurred image is out-of-focus blur; and obtain a fuzzy kernel by using an adaptive preset algorithm if the blurred image is motion blur, the preset algorithm includes at least one algorithm;
    图像复原模块,用于基于所述模糊核,利用超拉普拉斯先验去卷积算法获得清晰图像。An image restoration module is configured to obtain a clear image by using a super Laplacian prior deconvolution algorithm based on the fuzzy kernel.
  9. 根据权利要求8所述的一种模糊图像综合处理系统,其特征在于,所述系统还包括:A fuzzy image synthesis processing system according to claim 8, wherein the system further comprises:
    高分辨率处理模块,利用超分辨率重建技术处理所述清晰图像,获得高分辨率图像;a high resolution processing module that processes the sharp image using a super-resolution reconstruction technique to obtain a high resolution image;
    所述估计模块包括:The estimation module includes:
    离焦模糊单元,用于若模糊图像为离焦模糊时,采用刃边法提取基于梯度准则的最优刃边图像并根据所述最优刃边图像获取点扩散方程,所述点扩散方程即为所述模糊核;a defocusing unit for extracting an edge-edge image based on a gradient criterion and acquiring a point-diffusion equation according to the optimal edge-edge image if the blurred image is defocused, and the point-diffusion equation is The blur kernel;
    运动模糊单元,用于若模糊图像为运动模糊时,自适应基于稀疏先验的正则化方法获取模糊核。A motion blur unit is used to acquire a fuzzy kernel based on a sparse priori regularization method if the blurred image is motion blur.
  10. 根据权利要求8所述的一种模糊图像综合处理系统,其特征在于,所述离焦模糊单元包括:A blurred image synthesis processing system according to claim 8, wherein the defocus blur unit comprises:
    边缘检测子单元,用于对模糊图像进行Canny边缘检测获取边缘检测图像;An edge detection subunit, configured to perform Canny edge detection on the blurred image to obtain an edge detection image;
    直线检测子单元,用于对所述边缘检测图像进行霍夫变换获取阶跃边缘图像;a line detection subunit, configured to perform a Hough transform on the edge detection image to obtain a step edge image;
    图像截取子单元,用于以所述阶跃边缘图像中每一边缘的中心点为中心,截取每一边缘的刃边图像,所述刃边图像的尺寸大小根据所述离焦模糊图像大小和模糊程度选择一预设尺寸形成;An image intercepting subunit for intercepting a blade edge image of each edge centering on a center point of each edge of the step edge image, the size of the edge edge image being according to the defocused image size and The degree of blur selects a predetermined size to be formed;
    图像提取子单元,用于逐一计算所述刃边图像的梯度值,并根据梯度值提取最优刃边图像;An image extraction subunit, configured to calculate a gradient value of the edge edge image one by one, and extract an optimal edge edge image according to the gradient value;
    线性拟合子单元,用于利用最小二乘法对所述最优刃边图像中的边缘进行线性拟合获取最优刃边直线;a linear fitting sub-unit for linearly fitting edges in the optimal edge-edge image by using a least squares method to obtain an optimal edge line;
    散点图子单元,用于获取所述最优刃边图像中每个像素点到所述最优刃边直线的垂直距离作为横坐标,每个像素点的灰度值为纵坐标,形成散点图;a scatter plot sub-unit, configured to obtain a vertical distance of each pixel point in the optimal edge-edge image to the optimal edge line as an abscissa, and a gray value of each pixel point is an ordinate, forming a dispersion Dot map
    边缘扩散函数子单元,用于利用费米函数对所述散点图进行线性拟合获取边缘扩散函数;An edge diffusion function sub-unit for linearly fitting the scatter plot with a Fermi function to obtain an edge spread function;
    计算子单元,用于利用所述边缘扩散函数计算点扩散方程。A calculation subunit is configured to calculate a point spread equation using the edge spread function.
PCT/CN2018/091164 2018-01-31 2018-06-13 Comprehensive processing method and system for blurred image WO2019148739A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810094023.1A CN108305230A (en) 2018-01-31 2018-01-31 A kind of blurred picture integrated conduct method and system
CN201810094023.1 2018-01-31

Publications (1)

Publication Number Publication Date
WO2019148739A1 true WO2019148739A1 (en) 2019-08-08

Family

ID=62867331

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/091164 WO2019148739A1 (en) 2018-01-31 2018-06-13 Comprehensive processing method and system for blurred image

Country Status (2)

Country Link
CN (1) CN108305230A (en)
WO (1) WO2019148739A1 (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112082915A (en) * 2020-08-28 2020-12-15 西安科技大学 Plug-and-play type atmospheric particulate concentration detection device and detection method
CN112132758A (en) * 2020-08-05 2020-12-25 浙江工业大学 Image restoration method based on asymmetric optical system point spread function model
CN113177889A (en) * 2021-04-27 2021-07-27 深圳市慧鲤科技有限公司 Image processing method and device, electronic equipment and storage medium
CN114998282A (en) * 2022-06-16 2022-09-02 平安科技(深圳)有限公司 Image detection method, image detection device, electronic equipment and storage medium
US11704777B2 (en) 2021-08-27 2023-07-18 Raytheon Company Arbitrary motion smear modeling and removal
CN114998282B (en) * 2022-06-16 2024-05-31 平安科技(深圳)有限公司 Image detection method, device, electronic equipment and storage medium

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110148147B (en) * 2018-11-07 2024-02-09 腾讯大地通途(北京)科技有限公司 Image detection method, image detection device, storage medium and electronic device
CN109410152A (en) * 2018-11-26 2019-03-01 Oppo广东移动通信有限公司 Imaging method and device, electronic equipment, computer readable storage medium
CN109544577B (en) * 2018-11-27 2022-10-14 辽宁工程技术大学 Improved straight line extraction method based on edge point grouping
CN109521556A (en) * 2018-12-07 2019-03-26 歌尔科技有限公司 A kind of electron microscopic wearable device
CN109365877B (en) * 2018-12-20 2020-02-07 长安大学 Automatic positioning method for machining of parts of numerical control machine tool
CN109816620B (en) * 2019-01-31 2021-01-05 深圳市商汤科技有限公司 Image processing method and device, electronic equipment and storage medium
CN110211105B (en) * 2019-05-24 2021-09-07 厦门美柚股份有限公司 Blurred image detection method, device and system
CN110782413B (en) * 2019-10-30 2022-12-06 北京金山云网络技术有限公司 Image processing method, device, equipment and storage medium
CN111160486B (en) * 2019-12-31 2023-05-02 三峡大学 Fuzzy image classification method based on support vector machine and wavelet decomposition
CN111311562B (en) * 2020-02-10 2023-10-10 浙江华创视讯科技有限公司 Ambiguity detection method and device for virtual focus image
CN111652811B (en) * 2020-04-14 2022-09-16 中国人民解放军国防科技大学 Motion blurred image restoration method based on edge function and optimal window wiener filtering
CN113643193A (en) * 2020-05-11 2021-11-12 上海耕岩智能科技有限公司 Image deblurring method and device, imaging equipment and storage medium
US20220015622A1 (en) * 2020-07-15 2022-01-20 Tectus Corporation On-eye image processing
CN111899196A (en) * 2020-08-05 2020-11-06 沈阳工业大学 Blade defect motion blurred image restoration method based on classical restoration algorithm
US11615510B2 (en) * 2020-10-21 2023-03-28 Samsung Electronics Co., Ltd. Kernel-aware super resolution
CN113327206B (en) * 2021-06-03 2022-03-22 江苏电百达智能科技有限公司 Image fuzzy processing method of intelligent power transmission line inspection system based on artificial intelligence
CN115861048A (en) * 2021-09-24 2023-03-28 腾讯科技(深圳)有限公司 Image super-resolution method, device, equipment and storage medium
CN114842366B (en) * 2022-07-05 2022-09-16 山东中宇航空科技发展有限公司 Stability identification method for agricultural plant protection unmanned aerial vehicle
CN114972333B (en) * 2022-07-19 2022-10-25 淄博市淄川区市政环卫服务中心 Road crack detection method and system based on artificial intelligence
CN117058118B (en) * 2023-08-28 2024-04-05 北京航空航天大学 Method for evaluating image deterioration effect caused by plane shielding glass silk screen
CN117170083B (en) * 2023-10-24 2024-02-09 合肥埃科光电科技股份有限公司 Automatic focusing method, system and storage medium for superposition dynamic image quality compensation
CN117557626B (en) * 2024-01-12 2024-04-05 泰安大陆医疗器械有限公司 Auxiliary positioning method for spray head installation of aerosol sprayer

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104021529A (en) * 2014-06-17 2014-09-03 中国电子科技集团公司第二十八研究所 Blurred image non-blind restoration method
CN104331871A (en) * 2014-12-02 2015-02-04 苏州大学 Image de-blurring method and image de-blurring device
CN106339996A (en) * 2016-09-09 2017-01-18 江南大学 Image blind defuzzification method based on hyper-Laplacian prior
US20170104933A1 (en) * 2015-06-16 2017-04-13 Chung-Ang University-Academic Cooperation Foundation Method and apparatus for restoring motion blurred image

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102708551A (en) * 2012-05-17 2012-10-03 浙江大学 Image deconvolution method based on super laplace apriori constraint
US9542736B2 (en) * 2013-06-04 2017-01-10 Paypal, Inc. Evaluating image sharpness
CN104091314B (en) * 2014-07-22 2017-02-01 西北工业大学 Turbulence-degraded image blind restoration method based on edge prediction and sparse ratio regular constraints
CN104166961B (en) * 2014-07-28 2018-03-16 西安电子科技大学 The fuzzy core method of estimation that low-rank for blindly restoring image approaches
CN105913392A (en) * 2016-04-08 2016-08-31 西安电子科技大学昆山创新研究院 Degraded image overall quality improving method in complex environment
US20170316552A1 (en) * 2016-04-27 2017-11-02 Ramot At Tel-Aviv University Ltd. Blind image deblurring via progressive removal of blur residual
CN107084794B (en) * 2017-04-10 2021-06-22 东南大学 Flame three-dimensional temperature field measuring system and method based on light field layered imaging technology
CN107133948B (en) * 2017-05-09 2020-05-08 电子科技大学 Image blurring and noise evaluation method based on multitask convolution neural network
CN107316289B (en) * 2017-06-08 2020-05-08 华中农业大学 Method for dividing rice ears in field based on deep learning and superpixel division
CN107369134A (en) * 2017-06-12 2017-11-21 上海斐讯数据通信技术有限公司 A kind of image recovery method of blurred picture

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104021529A (en) * 2014-06-17 2014-09-03 中国电子科技集团公司第二十八研究所 Blurred image non-blind restoration method
CN104331871A (en) * 2014-12-02 2015-02-04 苏州大学 Image de-blurring method and image de-blurring device
US20170104933A1 (en) * 2015-06-16 2017-04-13 Chung-Ang University-Academic Cooperation Foundation Method and apparatus for restoring motion blurred image
CN106339996A (en) * 2016-09-09 2017-01-18 江南大学 Image blind defuzzification method based on hyper-Laplacian prior

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112132758A (en) * 2020-08-05 2020-12-25 浙江工业大学 Image restoration method based on asymmetric optical system point spread function model
CN112082915A (en) * 2020-08-28 2020-12-15 西安科技大学 Plug-and-play type atmospheric particulate concentration detection device and detection method
CN112082915B (en) * 2020-08-28 2024-05-03 西安科技大学 Plug-and-play type atmospheric particulate concentration detection device and detection method
CN113177889A (en) * 2021-04-27 2021-07-27 深圳市慧鲤科技有限公司 Image processing method and device, electronic equipment and storage medium
CN113177889B (en) * 2021-04-27 2022-03-01 深圳市慧鲤科技有限公司 Image processing method and device, electronic equipment and storage medium
US11704777B2 (en) 2021-08-27 2023-07-18 Raytheon Company Arbitrary motion smear modeling and removal
CN114998282A (en) * 2022-06-16 2022-09-02 平安科技(深圳)有限公司 Image detection method, image detection device, electronic equipment and storage medium
CN114998282B (en) * 2022-06-16 2024-05-31 平安科技(深圳)有限公司 Image detection method, device, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN108305230A (en) 2018-07-20

Similar Documents

Publication Publication Date Title
WO2019148739A1 (en) Comprehensive processing method and system for blurred image
Zhu et al. Single image dehazing using color attenuation prior.
US10755390B2 (en) Image deblurring method based on light streak information in an image
Hu et al. Deblurring low-light images with light streaks
Aslantas et al. A pixel based multi-focus image fusion method
Tiwari et al. Review of motion blur estimation techniques
EP3132418B1 (en) Non local image denoising
US9589328B2 (en) Globally dominant point spread function estimation
Yu et al. Efficient patch-wise non-uniform deblurring for a single image
CN108154491B (en) Image reflection eliminating method
Kumar Deblurring of motion blurred images using histogram of oriented gradients and geometric moments
WO2017100971A1 (en) Deblurring method and device for out-of-focus blurred image
Pertuz et al. Reliability measure for shape-from-focus
CN111640077A (en) Simple and efficient fuzzy text picture sharpening processing method
CN107945119B (en) Method for estimating correlated noise in image based on Bayer pattern
CN105574823A (en) Deblurring method and device for out-of-focus blurred image
CN117132503A (en) Method, system, equipment and storage medium for repairing local highlight region of image
Ahn et al. Segmenting a noisy low-depth-of-field image using adaptive second-order statistics
CN115147613A (en) Infrared small target detection method based on multidirectional fusion
Zhu Measuring spatially varying blur and its application in digital image restoration
Bansal et al. Digital image enhancement by improving contrast, removal of noise and motion blurring
CN117788306A (en) Multithreading-based multi-focal-length tab image fusion method
Gong et al. Image Correction
Lu et al. Parameters optimization in blind motion deblurring of UAV images
Feng et al. A restoration method for the turbulent degraded images based on the salient edge selection and the L0 norm constraint

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18903378

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18903378

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 25.03.2021)

122 Ep: pct application non-entry in european phase

Ref document number: 18903378

Country of ref document: EP

Kind code of ref document: A1