CN102867319B - Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images - Google Patents

Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images Download PDF

Info

Publication number
CN102867319B
CN102867319B CN201210361585.0A CN201210361585A CN102867319B CN 102867319 B CN102867319 B CN 102867319B CN 201210361585 A CN201210361585 A CN 201210361585A CN 102867319 B CN102867319 B CN 102867319B
Authority
CN
China
Prior art keywords
image
subimage
original image
degraded image
discrete
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.)
Active
Application number
CN201210361585.0A
Other languages
Chinese (zh)
Other versions
CN102867319A (en
Inventor
谭久彬
赵烟桥
刘俭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201210361585.0A priority Critical patent/CN102867319B/en
Publication of CN102867319A publication Critical patent/CN102867319A/en
Application granted granted Critical
Publication of CN102867319B publication Critical patent/CN102867319B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

The invention discloses a method for constructing an integral discrete degraded image in one-way movement by weighted superposition of sub-images, belongs to a section using more than one image in the field of general image data processing or generation, and particularly relates to a method for constructing a blurred image in discrete movement. The method includes constructing n sub-images fig i (i=1, 2, ..., n) according to n-pixel movement distances of an original image with the resolution of M*N along a row or column direction of the original image; and then carrying out weighted linear superposition for the obtained n sub-images according to a formula shown in the description. The n is smaller than the N if the original image moves along the row direction, and the n is smaller than the M if the original image moves along the column direction. In the formula, wi is a weighting coefficient, and a fig is the constructed discrete degraded image. The method for constructing the discrete degraded image has the advantages that operation time is short, a degradation process is visible and is easy to understand, and the image does not need to be readjusted.

Description

The discrete degraded image building method of view picture one-way movement of subimage weighted stacking
Technical field
The view picture one-way movement discrete degraded image building method of subimage weighted stacking to belong in general image real time transfer or generation field by using the part more than piece image, particularly relates to a kind of discrete motion blurred picture building method.
Background technology
If there is relative motion between imageing sensor and target in imaging process, gained image will produce motion blur phenomenon.At daily life, commercial production, aerospace field, this phenomenon is very general.Although motion blur image can present the aesthetic feeling of art in some special dimension, but in most field such as communications and transportation, commercial production, motion blur image but can only make troubles to the identification of target in image and to the acquisition of target detail information to us.As the electronic eyes in traffic and transport field, if the fog-level photographing image brings difficulty to the identification of license plate number, be so just difficult to press chapter punishment to vehicles peccancy, be unfavorable for the conventional maintenance of traffic order, cause potential safety hazard to the life of people.
Avoid the generation of motion blur phenomenon in image, generally adopt stabilization technology, stabilization technology comprises optical anti-vibration and electronic flutter-proof, optical anti-vibration is divided into again camera lens stabilization and imaging stabilization, camera lens stabilization refers to and arrange special stabilization compensating glass group in camera lens, according to jitter direction and the degree of camera, and the corresponding adjustment position of compensating glass group and angle, light path is made to keep stable, as Canon EF IS series camera lens, Nikon VR series camera lens, suitable horse OS series camera lens; Imaging stabilization refers to image device after perception camera shake, and the position or the angle that change image device are held in the stable of picture, and this technology was widely used in the digital camera epoch.Electronic flutter-proof refers to by analyzing become image, then utilizes the stabilization technology that algorithm compensates image, and this technology carrys out compensate for jitter indeed through reduction image quality, attempts between image quality and float, obtain an equilibrium point.Electronic flutter-proof, compared with optical anti-vibration, has cost low, the feature of weak effect, and therefore electronic flutter-proof is only used in low side camera.But the algorithm about electronic flutter-proof is but more subject to the concern of academia than optical anti-vibration technology.
For the algorithm of electronic flutter-proof, it is exactly the restoration algorithm of corresponding academia motion blur image, present stage, restoration algorithm was very many, there are traditional liftering algorithm, Wiener filtering algorithm, also has numerous blind restoration algorithm such as Kalman filtering algorithm and projections onto convex sets, up to now, the new algorithm of improvement is still had to continue to bring out out.In order to verify the adaptability of these new restoration algorithms, needing to restore the different degraded image of only degradation parameter, and comparing with original image.When each collection image, although we artificially can set the degradation parameter of image as required, but the impact of random noise cannot be avoided, the image sequence that actual acquisition is arrived is except degradation parameter difference, the impact of random noise will inevitably be subject to, therefore cannot actual acquisition to the different degraded image of only degradation parameter.
The way overcoming this problem is very simple, manually degenerates exactly by the mode of software simulation to the different degenrate function of original non degenerate imagery exploitation.By people's works such as Gonzalez, the people such as Ruan Qiuqi translation, and " Digital Image Processing " of being published by China Machine Press summarizes the method for two kinds of artificial degraded images that prior art adopts in book:
The first is spatial domain convolution degeneration method, if original image is f (x, y), degenrate function is h (x, y), then degraded image g (x, y) is expressed as:
g(x,y)=f(x,y)*h(x,y)
In formula, " * " represents convolution algorithm; For the discrete picture of M × N, the process that the first spatial domain convolution way of degeneration obtains degraded image can be write as further:
g ( x , y ) = Σ m = 1 M Σ n = 1 N f ( m , n ) h ( x - m , y - n )
In formula, x=1,2 ..., M; Y=1,2 ..., N.Can know according to formula above, calculate discrete degraded image g (x, y), need x, y, m, n completes quadruple loop computation and could realize, quadruple loop computation makes the computation process of discrete degraded image g (x, y) very consuming time, and this is the shortcoming of spatial domain convolution way of degeneration.
The second is that frequency domain Fourier degenerates method, if the frequency spectrum of original image f (x, y) is F (u, v), the frequency spectrum of degenrate function h (x, y) is H (u, v), then the frequency spectrum designation of degraded image g (x, y) is:
G(u,v)=F(u,v)H(u,v)
In formula, u=1,2 ..., M; V=1,2 ..., N.Due to the existence of Fast Fourier Transform (FFT) method, frequency domain Fourier method of degenerating is made to compare spatial domain convolution degeneration method and have on operation time and significantly promote, but, this method also has himself shortcoming: first, whole degenerative process completes in a frequency domain, degenerative process is neither directly perceived, not easily understands again; Secondly, by the frequency spectrum G (u of degraded image g (x, y), v) carry out inverse Fourier transform and obtain degraded image g (x, y) in process, also need image border to move to center, otherwise correspondingly with real image not go up.
Summary of the invention
The present invention is exactly the shortcoming for spatial domain convolution degeneration method length operation time, and frequency domain Fourier transform degeneration method intuitively, is not easily understood, and after inverse Fourier transform, also need shortcoming image being moved to operation, propose a kind of discrete degraded image building method of view picture one-way movement of subimage weighted stacking; The method not only operation time short, and degenerative process is directly perceived, is convenient to understand, and without the need to adjusting image again.
The object of the present invention is achieved like this:
The discrete degraded image building method of view picture one-way movement of subimage weighted stacking comprises the following steps:
A. be that the original image of M × N moves along its row or column direction the distance of n pixel according to resolution, construct n subimage fig i(i=1,2 ..., n), wherein:
If original image moves in the row direction, then n < N;
If original image moves along column direction, then n < M;
B. n the subimage obtained by step a is weighted linear superposition according to following formula:
fig = &Sigma; i = 1 n w i &CenterDot; fig i
In formula, w ifor weighting coefficient, fig is the discrete degraded image constructed.
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, the 1st subimage fig 1be expressed as:
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, original image moves upward, i-th subimage fig i(2≤i≤n) is expressed as:
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, original image moves downward, i-th subimage fig i(2≤i≤n) is expressed as:
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, original image to left movement, i-th subimage fig i(2≤i≤n) is expressed as:
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, original image moves right, i-th subimage fig i(2≤i≤n) is expressed as:
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step b, weighting coefficient w iratio be expressed as:
w 1 : w 2 : &CenterDot; &CenterDot; &CenterDot; : w n - 1 : w n = 1 v 1 : 1 v 2 : &CenterDot; &CenterDot; &CenterDot; : 1 v n - 1 : 1 v n
In formula, v ii-th subimage fig icorresponding image motion speed, and have: v i≠ 0 (i=1,2 ..., n).
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, the discrete degraded image fig obtained by step b carries out gray-scale value adjustment according to following formula:
fig_improve=k·fig
In formula, k is regulation coefficient, and fig_improve is the discrete degraded image after adjustment.
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, described regulation coefficient k is after the described discrete degraded image fig mean value divided by all grey scale pixel values of this image, then is multiplied by the 1st subimage fig 1the mean value of all grey scale pixel values.
The invention has the beneficial effects as follows:
1) due to discrete degraded image building method of the present invention be according to original image motion n pixel distance, construct n subimage, and linear, additive computing is weighted to these subimages, the operation time of the method is made up of three parts, Part I is the computing time of n weight, Part II is the computing time of the numeral that is multiplied with image of n respective weights and matrix multiple, Part III is the computing time that n the matrix representing n subimage is added, avoid four of spatial domain convolution degeneration method due to method of the present invention to recirculate, and each step computing is all the simplest arithmetic, therefore the method has short beneficial effect operation time,
2) due to n locus in n of the present invention subimage respectively corresponding imaging process residing for target, linear, additive computing is weighted to this n subimage and represents that the subimage by this n locus has different motion speed is recorded, whole degenerative process completes in spatial domain, the direct corresponding imaging process of calculating process, therefore to have degenerative process directly perceived for the method, be convenient to understand, and the beneficial effect without the need to adjusting again image.
Accompanying drawing explanation
Fig. 1 is original image.
Fig. 2 is the degraded image based on subimage stacking method.
Fig. 3 is the degraded image based on convolution method.
Fig. 4 is the degraded image based on Fourier transformation method.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the invention is described in further detail.
Fig. 1 does not have the resolution of degenerating to be the original image of 256 × 256, and in the present embodiment, original image is from its initial position along the distance of column direction upwards uniformly accelerated motion 10 pixels.And at its initial position, the movement velocity of original image is:
v 1 = 2 a &CenterDot; 2 L
In formula, v 1for the movement velocity of original image, a is acceleration, and L is pel spacing.
The discrete degraded image building method of view picture one-way movement of subimage weighted stacking comprises the following steps:
A. according to resolution be 256 × 256 original image to move along its column direction the distance of 10 pixels, construct 10 subimage fig 1~ fig 10, and meet 10 < 256;
Wherein, fig 1there is the function distribution that same original image is identical:
Element value in matrix represents the gray-scale value of this image correspondence position pixel.
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, in described step a, original image moves upward, i-th subimage fig i(2≤i≤10) are expressed as:
B. 10 subimages obtained by step a are weighted linear superposition according to following formula:
fig = &Sigma; i = 1 10 w i &CenterDot; fig i
In formula, w ifor weighting coefficient, fig is the discrete degraded image constructed.Can obtain by calculating, comprising original image 10 sub-image motion speed and being respectively:
v i = 2 a ( i + 1 ) L , ( i = 1,2 , &CenterDot; &CenterDot; &CenterDot; , 10 )
So weighting coefficient w iratio be expressed as:
w 1 : w 2 : &CenterDot; &CenterDot; &CenterDot; : w 10 = 1 v 1 : 1 v 2 : &CenterDot; &CenterDot; &CenterDot; : 1 v 10 = 1 2 : 1 3 : &CenterDot; &CenterDot; &CenterDot; : 1 11
In formula, v ii-th subimage fig icorresponding image motion speed, and have: v i≠ 0 (i=1,2 ..., 10).
The discrete degraded image building method of view picture one-way movement of above-mentioned subimage weighted stacking, the discrete degraded image fig obtained by step b carries out gray-scale value adjustment according to following formula:
fig_improve=k·fig
In formula, k is regulation coefficient, and fig_improve is the discrete degraded image after adjustment.Described regulation coefficient k is after the described discrete degraded image fig mean value divided by all grey scale pixel values of this image, then is multiplied by the 1st subimage fig 1the mean value of all grey scale pixel values, the degraded image obtained after adjustment as shown in Figure 2.
Short in order to verify that the view picture one-way movement discrete degraded image building method of subimage weighted stacking of the present invention has not only operation time further, and degenerative process is directly perceived, be convenient to understand, and the beneficial effect without the need to adjusting again image, method of spatial domain convolution degeneration method and frequency domain Fourier in method therefor of the present invention and prior art being degenerated contrasts.
No matter spatial domain convolution degeneration method or frequency domain Fourier degenerate method, and its spatial domain degenrate function h is the matrix of 256 × 256, and 119th ~ 128 row in this matrix, 10 element ratios of the 128th row are: rest of pixels is 0.According to spatial domain convolution degeneration method, utilize degraded image g (x, y) that g (x, y)=f (x, y) * h (x, y) obtains as shown in Figure 3; To degenerate method according to frequency domain Fourier, utilize G (u, v)=F (u, v) H (u, v) first obtains the frequency spectrum of degraded image g (x, y), again through inverse Fourier transform, and image border is moved to center, the degraded image g (x, y) obtained is as shown in Figure 4.
Fig. 2, Fig. 3, Fig. 4 are contrasted, the maximum difference of Fig. 2 and Fig. 3 respective pixel gray-scale value is only 2.5580 × 10 -13, the maximum difference of Fig. 2 and Fig. 4 respective pixel gray-scale value is only 3.4106 × 10 -13, the small difference between each method rounds off due to computer-internal mathematical operation and causes, completely negligible.This conclusion illustrates that the view picture one-way movement discrete degraded image building method of subimage weighted stacking of the present invention has same prior art spatial domain convolution degeneration method and frequency domain Fourier and to degenerate the identical deteriorating effect of method.
In addition, on operation time, the view picture one-way movement discrete degraded image building method used time 0.2500s of subimage weighted stacking of the present invention, spatial domain convolution degeneration method used time 343.6720s, frequency domain Fourier degenerates method used time 0.2030s, can find out, the present invention compares with spatial domain convolution degeneration method, has short beneficial effect operation time, compare with frequency domain Fourier method of degenerating, there is degenerative process directly perceived, be convenient to understand, and the beneficial effect without the need to adjusting again image.

Claims (1)

1. the discrete degraded image building method of the view picture one-way movement of subimage weighted stacking, is characterized in that said method comprising the steps of:
A. be that the original image of M × N moves along its row or column direction the distance of n pixel according to resolution, construct n subimage fig i(i=1,2 ..., n), wherein:
If original image moves in the row direction, then n≤N;
If original image moves along column direction, then n≤M;
B. n the subimage obtained by step a is weighted linear superposition according to following formula:
fig = &Sigma; i = 1 n w i fig i
In formula, w ifor weighting coefficient, fig is the discrete degraded image constructed;
C. the discrete degraded image fig obtained by step b carries out gray-scale value adjustment according to following formula:
fig_improve=k·fig
In formula, k is regulation coefficient, and fig_improve is the discrete degraded image after adjustment, and described regulation coefficient k is after the described discrete degraded image fig mean value divided by all grey scale pixel values of this image, then is multiplied by the 1st subimage fig 1the mean value of all grey scale pixel values;
In described step a, at the 1st subimage fig 1be expressed as:
When:
Original image moves upward, i-th subimage fig i(2≤i≤n) is expressed as:
Original image moves downward, i-th subimage fig i(2≤i≤n) is expressed as:
Original image to left movement, i-th subimage fig i(2≤i≤n) is expressed as:
Original image moves right, i-th subimage fig i(2≤i≤n) is expressed as:
In described step b, weighting coefficient w iratio be expressed as:
w 1 : w 2 : . . . : w n - 1 : w n = 1 v 1 : 1 v 2 : . . . : 1 v n - 1 : 1 v n
In formula, v ii-th subimage fig icorresponding image motion speed, and have: v i≠ 0 (i=1,2 ..., n).
CN201210361585.0A 2012-09-25 2012-09-25 Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images Active CN102867319B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210361585.0A CN102867319B (en) 2012-09-25 2012-09-25 Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210361585.0A CN102867319B (en) 2012-09-25 2012-09-25 Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images

Publications (2)

Publication Number Publication Date
CN102867319A CN102867319A (en) 2013-01-09
CN102867319B true CN102867319B (en) 2015-04-01

Family

ID=47446175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210361585.0A Active CN102867319B (en) 2012-09-25 2012-09-25 Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images

Country Status (1)

Country Link
CN (1) CN102867319B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102651134A (en) * 2012-03-17 2012-08-29 哈尔滨工业大学 Constant-speed blurred image construction method and device based on splicing of two frames of static images

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7885483B2 (en) * 2006-11-24 2011-02-08 National Synchrotron Radiation Research Center Image alignment method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102651134A (en) * 2012-03-17 2012-08-29 哈尔滨工业大学 Constant-speed blurred image construction method and device based on splicing of two frames of static images

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于相位相关的匀速直线运动模糊图像位移参数估计;孙辉等;《中国光学》;20120430;第5卷(第2期);174-180 *
运动模糊图像恢复的参数估算与算法研究;黎玮;《中国优秀硕士学位论文全文数据库 信息科技辑》;20120415(第4期);正文第14-15页、第20页 *

Also Published As

Publication number Publication date
CN102867319A (en) 2013-01-09

Similar Documents

Publication Publication Date Title
Wang et al. Fast image dehazing method based on linear transformation
CN113837938B (en) Super-resolution method for reconstructing potential image based on dynamic vision sensor
CN102369556B (en) Imaging device
CN109272014B (en) Image classification method based on distortion adaptive convolutional neural network
CN110996080A (en) Video presentation device, video presentation method, data structure, and program
CN103501401A (en) Real-time video de-noising method for super-loud noises based on pre-filtering
CN102024252B (en) Matrix rank minimization-based reconstruction method for recovering underwater distorted image
CN101872112B (en) Three-dimensional camera shooting automatic collecting system
CN102930565B (en) Discrete degraded image construction method for retracing moving target in static background
WO2019228450A1 (en) Image processing method, device, and equipment, and readable medium
CN113077505A (en) Optimization method of monocular depth estimation network based on contrast learning
CN112270691A (en) Monocular video structure and motion prediction method based on dynamic filter network
CN104574381A (en) Full reference image quality evaluation method based on LBP (local binary pattern)
CN112819697A (en) Remote sensing image space-time fusion method and system
CN112991167A (en) Aerial image super-resolution reconstruction method based on layered feature fusion network
CN114677479A (en) Natural landscape multi-view three-dimensional reconstruction method based on deep learning
CN102867318B (en) The discrete degraded image building method of view picture fold return motion of subimage weighted stacking
Li et al. Space–time super-resolution with patch group cuts prior
CN102867319B (en) Method for constructing integral discrete degraded image in one-way movement by weighted superposition of sub-images
CN102930566B (en) Discrete degraded image construction method for unidirectional movement object in static background
CN102930568B (en) Discrete degraded image construction method for uniform motion object in static background
CN110365966B (en) Video quality evaluation method and device based on window
CN102867320B (en) Method for constructing integral discrete degraded image in uniform motion on basis of superposition of sub-images
CN115311145B (en) Image processing method and device, electronic equipment and storage medium
US11900564B2 (en) Storage medium storing program, image processing apparatus, and training method of machine learning 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