CN103824262A - Improved rapid image restoration processing method using R-L iterative algorithm - Google Patents
Improved rapid image restoration processing method using R-L iterative algorithm Download PDFInfo
- Publication number
- CN103824262A CN103824262A CN201410058965.6A CN201410058965A CN103824262A CN 103824262 A CN103824262 A CN 103824262A CN 201410058965 A CN201410058965 A CN 201410058965A CN 103824262 A CN103824262 A CN 103824262A
- Authority
- CN
- China
- Prior art keywords
- image
- iteration
- processing method
- iterative
- iterative algorithm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000012804 iterative process Methods 0.000 claims abstract description 7
- 238000011084 recovery Methods 0.000 claims description 13
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
Disclosed is an improved rapid image restoration processing method using an R-L iterative algorithm. The method includes the steps that an image noise model is established, raw information of an observed image is acquired, and the information is given out in a matrix mode; the Gaussian kernel function is used as a point spread function, and an initial value is initialized to a restoration result of the observed image; numerical value iteration is carried out according to an iterative equation, and a final restored image is obtained after a plurality of times of iteration. According to the method, double regularization parameters are used, the Gaussian kernel function is used as the initial value, and the Alembert convergence criterion is used for conducting numerical value iteration in the iterative process; due to initial value selection, the operation efficiency is high, the iterative computation every time is more effective, and the number of iterations is greatly smaller than the number of iterations in an existing method; the method can be used for image restoration under the condition of high power noise, image fine details can be reserved, the luminance and the edge effect can be improved, the two regularization parameters can be flexibly selected, and the size of the point spread function is far smaller than the size of an original image.
Description
Technical field
The invention belongs to technical field of image processing, particularly a kind of rapid image Recovery processing method of modified R-L iterative algorithm.
Background technology
Image in fields such as medical treatment, astronomy, remote sensing has very high requirement for identification and sharpness, and larger but actual reception image is subject to noise effect, overriding noise type is poisson noise.At present, that relatively more conventional is R-L(Richardson-Lucy) iterative algorithm, but due to the ill characteristic with matrix inversion wherein, so have some technical schemes to propose the method for regularization, what some regularization method adopted is one-parameter regularization, cannot utilize well so each iterative process to realize calculating rapidly and efficiently, what some regularization adopted is the iterative process of partial differential equation, its arithmetic speed is slow, and complexity is high, high to the hardware requirement of calculating process.
Summary of the invention
In order to overcome the shortcoming of above-mentioned prior art, the object of the present invention is to provide a kind of rapid image Recovery processing method of modified R-L iterative algorithm, utilize biregular parameter, adopt gaussian kernel function as initial value, in iterative process, use and reach bright Bel's convergence criterion and carry out numerical value iteration.
To achieve these goals, the technical solution used in the present invention is:
A rapid image Recovery processing method for modified R-L iterative algorithm, sets up image noise model i=h*o+n, and wherein i is the image observing, and h is point spread function, and o is real image, and n is noise, and process comprises the steps:
Step 1. obtains the raw information of the image observing, and provides in the mode of matrix;
Step 2. is used gaussian kernel function for point spread function h and the image restoration result i initialize observing;
Step 3. is carried out numerical value iteration according to following two iterative equations:
In formula, the region at s presentation video place, m is iterations, o
(m)represent the result of the o of the m time iteration, h
(m)the result that represents the h of the m time iteration, * represents convolution algorithm, represents multiplying, o
*and h
*represent respectively the conjugation of o and h, λ
1and λ
2represent regularization parameter, be constant, Σ
so (s) is illustrated in the numerical value summation to o in whole image-region, and div () represents the variable in bracket to ask divergence computing,
represent gradient operator, || || represent norm;
In iterative process, adopt alternately iteration, first use o
(m)and h
(m)iteration obtains o
(m+1), re-use o
(m+1)and h
(m)obtain h
(m+1), or first use o
(m)and h
(m)iteration obtains h
(m+1), re-use h
(m+1)and o
(m)obtain o
(m+1);
Through the o obtaining after several times iteration
(m)the i.e. final image recovering.
In described step 1, if original image is gray scale image, be single matrix; If original image is coloured image, is divided into tri-components of RGB and processes respectively.
In described step 1, described matrix is the picture element matrix of the image that observes.
In described step 2, use gaussian kernel function, provide its size and standard deviation, can complete for the initialized assignment of point spread function, re-use this point spread function and carry out convolution one time with the image receiving, can complete for the initialization assignment of wanting the true picture recovering.
In described step 3, the data scale of h is far smaller than the data scale of i and o, and data scale refers to the size of matrix.
In described step 3, the region at image place is two-dimensional coordinate region.
The numerical value of iterations n is generally no more than 10.
Compared with prior art, beneficial effect of the present invention is that advantage is:
Initial value choose high, the each iterative computation of the operation efficiency causing more effectively, iterations is significantly less than existing method, can use this method to carry out image recovery under compared with very noisy condition, can retain image fine detail, promote brightness and edge effect, two regularization parameters are chosen very flexibly, the size of point spread function is far smaller than the size of original image.
Accompanying drawing explanation
Fig. 1 is the original image of the moonscape of spacecraft camera.
Fig. 2 is the enlarged drawing of oval part in Fig. 1.
Fig. 3 uses the present invention to carry out the result figure after iterative processing Fig. 2.
Embodiment
Describe embodiments of the present invention in detail below in conjunction with drawings and Examples.
As shown in Figure 1, be the original image of the moonscape of spacecraft camera.Process from wherein intercepting a part, as shown in Figure 2, it is very fuzzy as seen for cut-away view picture, luminance shortage, and it is the image i observing.
Utilize image Recovery processing method of the present invention, set up image noise model i=h*o+n, wherein i is the image observing, i.e. Fig. 2, and h is point spread function, is of a size of 15 × 15; O is real image, is the object that finally will solve, and n is noise, and process comprises the steps:
Step 1. obtains the raw information of the image i observing, provides in the mode of matrix, if original image is gray scale image, be single matrix; If original image is coloured image, is divided into tri-components of RGB and processes respectively.The image that the present embodiment observes is for colored, and the line number of this matrix and columns equal respectively the height and width of image pixel, and each element of matrix equals the numerical value of image correspondence position pixel.
Step 2. is used gaussian kernel function, provides its size and standard deviation, is point spread function h and the image restoration result i initialize observing; Re-use this point spread function and carry out convolution one time with the image receiving, can complete for the initialization assignment of wanting the true picture recovering.In the present embodiment, choose gaussian kernel function and be of a size of 15, standard deviation is 2.
The associating deconvolution formula of step 3. based on maximum a posteriori probability is as follows:
The likelihood maximum of taking the logarithm:
It is carried out to the regularization based on total variation:
Above formula is got variation to o and h respectively, can obtain:
To reaching bright Bel's convergence iteration above, there are following two formulas (for iteration time):
In formula, the region at s presentation video place (for example two-dimensional coordinate region), m is iterations, o
(m)represent the result of the o of the m time iteration, h
(m)the result that represents the h of the m time iteration, * represents convolution algorithm, represents multiplying, o* and h* represent respectively the conjugation of o and h, λ
1and λ
2represent regularization parameter, be constant, the selection range of total variation parameter is larger, can choose according to the intensity of actual noise, and the present embodiment is chosen λ
1=0.1, λ
2=0.03, Σ
so (s) is illustrated in the numerical value summation to o in whole image-region, and div () represents the variable in bracket to ask divergence computing,
represent gradient operator, || represent norm;
In iterative process, adopt alternately iteration, first use o
(m)and h
(m)iteration obtains o
(m+1), re-use o
(m+1)and h
(m)obtain h
(m+1), or first use o
(m)and h
(m)iteration obtains h
(m+1), re-use h
(m+1)and o
(m)obtain o
(m+1); Wherein the data scale of h can be far smaller than the data scale of i and o, and data scale refers to the size of matrix.
The o of the present embodiment through obtaining after 7 iteration
(7)the i.e. final image recovering, as shown in Figure 3, its brightness and sharpness have all had very large lifting to the result after iteration as seen.
Claims (7)
1. a rapid image Recovery processing method for modified R-L iterative algorithm, is characterized in that, sets up image noise model i=h*o+n, and wherein i is the image observing, and h is point spread function, and o is real image, and n is noise, and process comprises the steps:
Step 1. obtains the raw information of the image observing, and provides in the mode of matrix;
Step 2. is used gaussian kernel function for point spread function h and the image restoration result i initialize observing;
Step 3. is carried out numerical value iteration according to following two iterative equations:
In formula, the region at s presentation video place, m is iterations, o
(m)represent the result of the o of the m time iteration, h
(m)the result that represents the h of the m time iteration, * represents convolution algorithm, represents multiplying, o
*and h
*represent respectively the conjugation of o and h, λ
1and λ
2represent regularization parameter, be constant, Σ
so (s) is illustrated in the numerical value summation to o in whole image-region, and div () represents the variable in bracket to ask divergence computing,
represent gradient operator, || || represent norm;
In iterative process, adopt alternately iteration, first use o
(m)and h
(m)iteration obtains o
(m+1), re-use o
(m+1)and h
(m)obtain h
(m+1), or first use o
(m)and h
(m)iteration obtains h
(m+1), re-use h
(m+1)and o
(m)obtain o
(m+1);
Through the o obtaining after several times iteration
(m)the i.e. final image recovering.
2. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, is characterized in that, in described step 1, if original image is gray scale image, is single matrix; If original image is coloured image, is divided into tri-components of RGB and processes respectively.
3. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, is characterized in that, in described step 1, described matrix is the picture element matrix of the image that observes.
4. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, it is characterized in that, in described step 2, use gaussian kernel function, provide its size and standard deviation, can complete for the initialized assignment of point spread function, re-use this point spread function and carry out convolution one time with the image receiving, can complete for the initialization assignment of wanting the true picture recovering.
5. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, is characterized in that, in described step 3, the data scale of h is far smaller than the data scale of i and o, and data scale refers to the size of matrix.
6. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, is characterized in that, in described step 3, the region at image place is two-dimensional coordinate region.
7. the rapid image Recovery processing method of modified R-L iterative algorithm according to claim 1, is characterized in that, in described step 3, the numerical value of iterations n is not more than 10.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410058965.6A CN103824262B (en) | 2014-02-21 | 2014-02-21 | The rapid image of a kind of modified model R-L iterative algorithm recovers processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410058965.6A CN103824262B (en) | 2014-02-21 | 2014-02-21 | The rapid image of a kind of modified model R-L iterative algorithm recovers processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103824262A true CN103824262A (en) | 2014-05-28 |
CN103824262B CN103824262B (en) | 2016-07-13 |
Family
ID=50759308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410058965.6A Expired - Fee Related CN103824262B (en) | 2014-02-21 | 2014-02-21 | The rapid image of a kind of modified model R-L iterative algorithm recovers processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103824262B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574302A (en) * | 2014-12-25 | 2015-04-29 | 深圳市一体太赫兹科技有限公司 | Terahertz image restoration method and system |
CN107507135A (en) * | 2017-07-11 | 2017-12-22 | 天津大学 | Image reconstructing method based on coding aperture and target |
CN108230251A (en) * | 2016-12-22 | 2018-06-29 | 深圳清华大学研究院 | Combined type image recovery method and device |
CN109243384A (en) * | 2018-11-09 | 2019-01-18 | 京东方科技集团股份有限公司 | Show equipment and its driving method, driving device and computer-readable medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101742050A (en) * | 2009-12-03 | 2010-06-16 | 浙江大学 | Method for restoring TDICCD image aiming at motion fuzzy core space shift variant |
CN102682437A (en) * | 2012-05-17 | 2012-09-19 | 浙江大学 | Image deconvolution method based on total variation regularization |
CN102708551A (en) * | 2012-05-17 | 2012-10-03 | 浙江大学 | Image deconvolution method based on super laplace apriori constraint |
-
2014
- 2014-02-21 CN CN201410058965.6A patent/CN103824262B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101742050A (en) * | 2009-12-03 | 2010-06-16 | 浙江大学 | Method for restoring TDICCD image aiming at motion fuzzy core space shift variant |
CN102682437A (en) * | 2012-05-17 | 2012-09-19 | 浙江大学 | Image deconvolution method based on total variation regularization |
CN102708551A (en) * | 2012-05-17 | 2012-10-03 | 浙江大学 | Image deconvolution method based on super laplace apriori constraint |
Non-Patent Citations (2)
Title |
---|
YU-WING TAI ET AL: "《Richardson-Lucy Deblurring for Scenes under a Projective Motion Path》", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》, vol. 33, no. 8, 31 August 2011 (2011-08-31), pages 1603 - 1618 * |
赵博等: "《基于Richardson-Lucy 的图像去模糊新算法》", 《计算机工程与应用》, vol. 47, no. 34, 31 December 2011 (2011-12-31) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574302A (en) * | 2014-12-25 | 2015-04-29 | 深圳市一体太赫兹科技有限公司 | Terahertz image restoration method and system |
CN108230251A (en) * | 2016-12-22 | 2018-06-29 | 深圳清华大学研究院 | Combined type image recovery method and device |
CN107507135A (en) * | 2017-07-11 | 2017-12-22 | 天津大学 | Image reconstructing method based on coding aperture and target |
CN107507135B (en) * | 2017-07-11 | 2020-04-24 | 天津大学 | Image reconstruction method based on coding aperture and target |
CN109243384A (en) * | 2018-11-09 | 2019-01-18 | 京东方科技集团股份有限公司 | Show equipment and its driving method, driving device and computer-readable medium |
US11195479B2 (en) | 2018-11-09 | 2021-12-07 | Boe Technology Group Co., Ltd. | Display device and method for driving the same, driving apparatus and computer-readable medium |
Also Published As
Publication number | Publication date |
---|---|
CN103824262B (en) | 2016-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Plug-and-play image restoration with deep denoiser prior | |
Wu et al. | Fast end-to-end trainable guided filter | |
US11715258B2 (en) | Method for reconstructing a 3D object based on dynamic graph network | |
Lai et al. | Blur kernel estimation using normalized color-line prior | |
CN104392243B (en) | Hyperspectral image nonlinear solution mixing method based on the sparse Non-negative Matrix Factorization of core | |
US9171226B2 (en) | Image matching using subspace-based discrete transform encoded local binary patterns | |
CN106408530A (en) | Sparse and low-rank matrix approximation-based hyperspectral image restoration method | |
CN104680491A (en) | Non-uniform image motion blur removing method based on deep neural network | |
CN103150713A (en) | Image super-resolution method of utilizing image block classification sparse representation and self-adaptive aggregation | |
CN110335221B (en) | Multi-exposure image fusion method based on unsupervised learning | |
CN103824262A (en) | Improved rapid image restoration processing method using R-L iterative algorithm | |
Yang et al. | Designing display pixel layouts for under-panel cameras | |
CN104036468A (en) | Super-resolution reconstruction method for single-frame images on basis of pre-amplification non-negative neighbor embedding | |
CN102592129B (en) | Scenario-driven image characteristic point selection method for smart phone | |
CN106529549B (en) | Vision significance detection method based on self-adaptive features and discrete cosine transform | |
US10521918B2 (en) | Method and device for filtering texture, using patch shift | |
Yang et al. | Low‐light image enhancement based on Retinex decomposition and adaptive gamma correction | |
CN109146792A (en) | Chip image super resolution ratio reconstruction method based on deep learning | |
Sun et al. | Tensor Gaussian process with contraction for multi-channel imaging analysis | |
Kaur et al. | QRFODD: Quaternion Riesz fractional order directional derivative for color image edge detection | |
CN106296749A (en) | RGB D image eigen decomposition method based on L1 norm constraint | |
CN112801899A (en) | Internal and external circulation driving image blind deblurring method and device based on complementary structure perception | |
Bertaccini et al. | An adaptive norm algorithm for image restoration | |
CN103413306B (en) | A kind of Harris angular-point detection method of adaptive threshold | |
Shi et al. | Coded diffraction imaging via double sparse regularization model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160713 |