Image restoration method of blurred image, storage medium and terminal

A technology for blurring images and images, which is applied in the field of image restoration and can solve problems such as noise amplification

Pending Publication Date: 2021-06-08
UNIV OF ELECTRONIC SCI & TECH OF CHINA
0 Cites 0 Cited by

AI-Extracted Technical Summary

Problems solved by technology

[0004] The purpose of the present invention is to overcome the problem of noise amplification when the noise increases in the e...
View more

Method used

The present invention carries out double-tree complex wavelet decomposition twice to be restored image, and the image after decomposing has time-frequency local analysis characteristic, has translation invariance and multi-direction selectivity approximately simultaneously, can better keep image detail, avoids Gibuss phenomenon.
The present invention utilizes wavelet decomposition to separate eac...
View more

Abstract

The invention discloses an image restoration method of a blurred image, a storage medium and a terminal, and belongs to the technical field of image restoration, and the method comprises the following steps: carrying out the wavelet decomposition of a to-be-restored image, and obtaining a multi-band decomposition image; calculating an autocorrelation function of the decomposed image of each frequency band so as to obtain an initial point spread function of the decomposed image; carrying out R-L iteration processing on the decomposed image according to the initial point spread function; and carrying out wavelet inverse transformation on the decomposed image subjected to R-L iteration processing to realize image restoration processing. According to the method, wavelet decomposition is utilized to separate each frequency band component in the to-be-restored image, and R-L iteration processing is performed on the decomposed image based on the point spread function of each decomposed image, so that the image noise amplification phenomenon can be effectively reduced, the image detail information can be well stored, the precision of the deconvolution algorithm is improved, and a clear restored image is obtained.

Application Domain

Image enhancement

Technology Topic

Frequency bandImage based +8

Image

  • Image restoration method of blurred image, storage medium and terminal
  • Image restoration method of blurred image, storage medium and terminal
  • Image restoration method of blurred image, storage medium and terminal

Examples

  • Experimental program(3)

Example Embodiment

[0041] Example 1
[0042] like Figure 1-2 As shown, in embodiment 1, a kind of image restoration method of blurred image specifically comprises the following steps:
[0043] S11: performing wavelet decomposition on the image to be restored to obtain a multi-band decomposition image;
[0044] S12: Calculate the autocorrelation function of the decomposed image in each frequency band, and then obtain the initial point spread function of the decomposed image;
[0045] S13: Perform R-L iterative processing on the decomposed image according to the initial point spread function;
[0046] S14: Perform wavelet inverse transform on the decomposed image after R-L iterative processing to realize image restoration processing.
[0047] The present invention uses wavelet decomposition to separate the components of each frequency band in the image to be restored (fuzzy image), and performs R-L iterative processing on the decomposed image based on the point spread function of each decomposed image, which can effectively reduce the phenomenon of image noise amplification and at the same time, better Save image detail information, improve the accuracy of deconvolution algorithm, and then get a clear restored image.
[0048] Further, performing wavelet decomposition on the image to be restored is specifically performing dual-tree complex wavelet decomposition on the image to be restored twice. Specifically, the highest decomposition level of dual-tree complex wavelet decomposition for the image to be restored is M, M≤10, and M=2 in this embodiment. After two decompositions of the image to be restored, multiple high-frequency decomposed images and multiple A low frequency band decomposes the image.
[0049]More specifically, in this embodiment, the blurred image is recorded as image I, and the size of image I is m*n. After image I undergoes the first dual-tree complex wavelet transform (decomposition), the image I is filtered And after downsampling, the size obtained is The real decomposition image I xx , imaginary decomposition image I' xx , based on the first dual-tree complex wavelet decomposition of the image I, the second dual-tree complex wavelet decomposition is performed, that is, after filtering and down-sampling processing, the size is The real decomposition image I xxx , imaginary decomposition image I' xxx ,like Figure 3-4 as shown, image 3 middle h 0 、h 1 is the conjugate orthogonal filter pair of tree A, g 0 , g 1 is the conjugate orthogonal filter pair of tree B, and ↓ represents downsampling. filter pair h 0 、h 1 The corresponding real scaling function φ h (t) and wavelet function for:
[0050]
[0051]
[0052] filter pair g 0 , g 1 The corresponding real scaling function φ g (t) and wavelet function for:
[0053]
[0054]
[0055] This embodiment Figure 4 In, I 01 , I 10 , I 11 The high-frequency detail images in the directions of ±15°, ±45°, and ±75° of the tree A part after the original image I completes a dual-tree complex wavelet transform, I' 01 , I' 10 , I' 11 Respectively, the original image I completes a dual-tree complex wavelet transform and the high-frequency detail images in the directions of ±15°, ±45°, and ±75° of the tree B part, I 01 , I 10 , I 11 and I' 01 , I' 10 , I' 11 The image dimensions of the decomposed image are I 001 , I 010 , I 011 After the original image I completes the first dual-tree complex wavelet transform, the low-frequency image I 00 The high-frequency detail images in the directions of ±15°, ±45°, and ±75° obtained after double-tree complex wavelet decomposition again, I' 001 , I' 010 , I' 011 After the original image I completes the first dual-tree complex wavelet transform, the low-frequency image I' 00 The high-frequency detailed images in the directions of ±15°, ±45°, and ±75° obtained after double-tree complex wavelet decomposition again, I 001 , I 010 , I 011 and I' 001 , I' 010 , I' 011 The image dimensions of the decomposed image are while I' 000 , I 000 is the low-frequency information obtained after two times of dual-tree complex wavelet decomposition, and the image size is
[0056] The invention performs double-tree complex wavelet decomposition on the image to be restored twice, and the decomposed image has time-frequency local analysis characteristics, and has approximate translation invariance and multi-directional selectivity, can better retain image details, and avoid Gibuss phenomenon.
[0057] Further, the calculation formula of the autocorrelation function R of each frequency band decomposition image in step S12 is:
[0058] R=g(x,y)*g(x,y)
[0059] Among them, g(x, y) is the decomposed image of each frequency band, * represents the complex conjugate, and x, y represent the pixel coordinates of the image respectively.
[0060] Further, decompose the initial point spread function (PSF) h of the image 0 The calculation formula is:
[0061] h 0 =R-min(R)+ε[max(R)-min(R)]
[0062] Among them, ε represents the dynamic change percentage of the autocorrelation function R, which is generally a small non-negative integer, and the value here is 0.01.
[0063] Further, performing R-L iterative processing on the decomposed image according to the initial point spread function also includes:
[0064] The corresponding number of R-L iterations is determined according to the size and frequency band of each decomposed image. The present invention determines the corresponding number of R-L iterations according to the size and frequency band of the decomposed image, and can maximize the removal of noise in each decomposed image according to the characteristics of each decomposed image.
[0065] Further, the number of R-L iterations of each decomposed image follows:
[0066] The number of R-L iterations of the large-scale image < the number of R-L iterations of the high-frequency decomposition image < the number of R-L iterations of the low-frequency decomposition image, that is, the number of R-L algorithm iterations is used for the high-frequency detail part, and the low-frequency image part uses the R-L algorithm for multiple iterations. Larger images require fewer iterations of the R-L algorithm.
[0067] In this example, I 000 size is The low-frequency image of I 001 , I 010 , I 011 and I' 001 , I' 010 , I' 011 is the size of The high-frequency detail image of I 01 , I 10 , I 11 and I' 01 , I' 01 , I' 11 is the size of High-frequency detailed images, such as image 3 , Figure 5 As shown, will decompose the image I 01 , I 10 , I 11 , I' 01 , I' 01 , I' 11 The number of iterations is set to N, then the decomposition I 001 , I 010 , I 011 and I' 001 , I' 010 , I' 011 The number of iterations is 2N, decomposing the image I 000 The number of iterations is 3N. Furthermore, the number of R-L iterations of the large-scale image is greater than or equal to 10 (the number of iterations of the large-scale high-frequency decomposition image (I 01 ) 001 ) 000 )), then in this embodiment, the decomposition image I 01 , I 10 , I 11 , I' 01 , I' 01 , I' 11 The number of iterations is set to 10, then the decomposition I 001 , I 010 , I 011 and I' 001 , I' 010 , I' 011 The number of iterations is 20, decomposing the image I 000 The number of iterations is 30.
[0068] Further, the calculation formula for R-L iterative processing of the decomposed image according to the initial point spread function is:
[0069]
[0070]
[0071] Among them, n is the number of iterations, h(x,y) is the point spread function, h 0 is the initial point spread function, and f(x,y) is the restored image.
[0072] Finally, the dual-tree complex wavelet inverse transform is performed on the decomposed image f(x, y) that has completed R-L iterative processing, namely image 3 The inverse operation of , a clear restored image M can be obtained.

Example Embodiment

[0073] Example 2
[0074] This embodiment provides a storage medium having the same inventive concept as Embodiment 1, on which computer instructions are stored, and the steps of an image restoration method for a blurred image in Embodiment 1 are executed when the computer instructions are run.
[0075] Based on this understanding, the technical solution of this embodiment is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium. Several instructions are included to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic disk or optical disk and other various media that can store program codes. .

Example Embodiment

[0076] Example 3
[0077] This embodiment also provides a terminal, which has the same inventive concept as Embodiment 1, and includes a memory and a processor. The memory stores computer instructions that can be run on the processor, and the processor executes the computer instructions in Embodiment 1 when running the computer instructions. The invention relates to steps of an image restoration method for a blurred image. The processor may be a single-core or multi-core central processing unit or a specific integrated circuit, or one or more integrated circuits configured to implement the present invention.
[0078] Each functional unit in the embodiments provided by the present invention may be integrated into one processing unit, or each unit may physically exist separately, or two or more units may be integrated into one unit.

PUM

no PUM

Description & Claims & Application Information

We can also present the details of the Description, Claims and Application information to help users get a comprehensive understanding of the technical details of the patent, such as background art, summary of invention, brief description of drawings, description of embodiments, and other original content. On the other hand, users can also determine the specific scope of protection of the technology through the list of claims; as well as understand the changes in the life cycle of the technology with the presentation of the patent timeline. Login to view more.
Who we serve
  • R&D Engineer
  • R&D Manager
  • IP Professional
Why Eureka
  • Industry Leading Data Capabilities
  • Powerful AI technology
  • Patent DNA Extraction
Social media
Try Eureka
PatSnap group products