CN107481197A - Image Fuzzy Complementary characterizing method - Google Patents
Image Fuzzy Complementary characterizing method Download PDFInfo
- Publication number
- CN107481197A CN107481197A CN201710558860.0A CN201710558860A CN107481197A CN 107481197 A CN107481197 A CN 107481197A CN 201710558860 A CN201710558860 A CN 201710558860A CN 107481197 A CN107481197 A CN 107481197A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- image
- mtd
- gradient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种图像互补性表征方法,特别是涉及一种图像模糊互补性表征方法。The invention relates to an image complementarity characterization method, in particular to an image fuzzy complementarity characterization method.
背景技术Background technique
由于运动模糊模型的参数不同,会形成不同模糊程度的降质图像,不同模糊图像之间反映了目标的不同信息,这些信息就称为模糊图像的互补性。图像间互补性能够为模糊图像的清晰化重建提供更多信息。Due to the different parameters of the motion blur model, degraded images with different blur degrees will be formed, and different blurred images reflect different information of the target, which is called the complementarity of blurred images. The complementarity between images can provide more information for the sharp reconstruction of blurred images.
文献“Li W,Zhang J,Dai Q.Exploring aligned complementary image pairfor blind motion deblurring[C],IEEE Conference on Computer Vision and PatternRecognition,2011:273-280”提出了一种基于两幅图像的模糊图像复原算法,其利用了图像之间的互补信息进行多帧复原,但并未明确说明图像之间互补性的大小,并且,在多帧图像复原的选帧过程中,该方法并未提及选帧的方式,对于图像质量参差不齐的序列多帧图像,选帧方式的不同严重影响图像复原的质量。The paper "Li W, Zhang J, Dai Q. Exploring aligned complementary image pair for blind motion deblurring[C], IEEE Conference on Computer Vision and Pattern Recognition, 2011: 273-280" proposed a blurred image restoration algorithm based on two images , which utilizes the complementary information between images to perform multi-frame restoration, but does not clearly specify the size of the complementarity between images, and, in the frame selection process of multi-frame image restoration, this method does not mention the frame selection For a sequence of multi-frame images with uneven image quality, different frame selection methods seriously affect the quality of image restoration.
发明内容Contents of the invention
为了克服现有图像互补性表征方法图像复原质量差的不足,本发明提供一种图像模糊互补性表征方法。该方法首先采用基于图像滤波的方法提取模糊图像的互补性特征,通过对多帧模糊图像进行高频滤波,得到图像对应的高频细节信息;针对图像高频信息进行梯度卷积操作,分别提取水平方向以及竖直方向的图像梯度信息;对模糊图像的高频梯度图进行二值化操作,得到模糊图像的高频梯度二值图,该二值图即为模糊图像的互补性特征。然后基于提取出的模糊图像的互补性特征,利用互补性计算公式进行互补性定量计算,最终得到模糊图像间的互补性数值大小。该方法对不同模糊图像的互补性进行衡量,为模糊图像多帧复原提供了更加明确的图像选帧依据,对于图像质量参差不齐的序列多帧图像,基于互补性表征的图像选帧方式提高了图像复原的质量。In order to overcome the shortcomings of poor image restoration quality in existing image complementarity characterization methods, the present invention provides an image fuzzy complementarity characterization method. This method first uses the method based on image filtering to extract the complementary features of the fuzzy image, and obtains the high-frequency detail information corresponding to the image by performing high-frequency filtering on multiple frames of fuzzy images; performs gradient convolution operations on the high-frequency information of the image to extract Image gradient information in the horizontal direction and vertical direction; binarize the high-frequency gradient map of the blurred image to obtain a high-frequency gradient binary map of the blurred image, which is the complementary feature of the blurred image. Then, based on the extracted complementarity features of the fuzzy images, the complementarity calculation formula is used to carry out quantitative calculation of complementarity, and finally the value of complementarity between fuzzy images is obtained. This method measures the complementarity of different blurred images, and provides a clearer basis for image frame selection for multi-frame restoration of blurred images. quality of image restoration.
本发明解决其技术问题所采用的技术方案是:一种图像模糊互补性表征方法,其特点是包括以下步骤:The technical scheme adopted by the present invention to solve the technical problem is: a kind of image fuzzy complementarity characterization method, which is characterized in that it includes the following steps:
步骤一、针对大小为256×256的序列图像中的任意两幅模糊图像g1和g2,采用高斯滤波器G进行图像滤波分解,计算得到模糊图像的高频分量h1和h2:Step 1. For any two blurred images g 1 and g 2 in the sequence image with a size of 256×256, use the Gaussian filter G to perform image filtering and decomposition, and calculate the high-frequency components h 1 and h 2 of the blurred image:
其中,为卷积算子,二维高斯滤波卷积算子G如公式(2)所示:in, is the convolution operator, and the two-dimensional Gaussian filter convolution operator G is shown in formula (2):
分别在水平方向和竖直方向上对各高频分量图像hi做一阶梯度微分,即水平和竖直梯度算子卷积操作,得到水平方向梯度图像hdix和竖直方向梯度图像hdiy:Perform first-order gradient differentiation on each high-frequency component image h i in the horizontal direction and vertical direction, that is, the horizontal and vertical gradient operator convolution operation, to obtain the horizontal gradient image hd ix and the vertical gradient image hd iy :
其中,dx为水平梯度算子[1,-1],dy为竖直梯度算子[1,-1]T。在计算水平方向和竖直方向梯度图像的基础上,进一步计算提取图像的全局梯度图像hdi:Among them, d x is the horizontal gradient operator [1,-1], d y is the vertical gradient operator [1,-1] T . On the basis of calculating the gradient images in the horizontal direction and the vertical direction, further calculate the global gradient image hd i of the extracted image:
对计算得到的梯度图像hdi进行图像二值化操作,采用公式(6)提取模糊图像的互补性特征Ti(x,y)。Perform image binarization on the calculated gradient image hd i , and use formula (6) to extract the complementary feature T i (x, y) of the blurred image.
其中,二值化阈值ki为:Among them, the binarization threshold k i is:
式中,max(hdi)表示提取梯度图像hdi的最大像素灰度值,min(hdi)表示提取梯度图像hdi的最小像素灰度值。In the formula, max(hd i ) represents the maximum pixel gray value of the extracted gradient image hd i , and min(hd i ) represents the minimum pixel gray value of the extracted gradient image hd i .
步骤二、对于计算得到的二值图像T1和T2,以其作为模糊图像的互补性特征,基于数学集合论的方法,采用公式(8),对模糊图像g1和g2的互补性进行表征。Step 2. For the calculated binary images T 1 and T 2 , use them as the complementary features of the fuzzy image, based on the method of mathematical set theory, use the formula (8), the complementarity of the fuzzy images g 1 and g 2 To characterize.
公式(8)表示模糊图像g1和g2之间的互补性大小。其中,T1∪T2表示计算二值图像T1与T2的并集,具体方法是图像对应位置像素进行点对点或运算,计算结果是一个256×256的二值图像,T1∩T2表示计算二值图像T1与T2的交集,具体方法是对应位置像素进行点对点与运算,计算结果也是一个256×256的二值图像。||T||1表示对二值图像T计算其l1范数:Formula ( 8 ) expresses the magnitude of complementarity between blurred images g1 and g2. Among them, T 1 ∪ T 2 means to calculate the union of binary images T 1 and T 2. The specific method is to perform point-to-point OR operation on the corresponding position pixels of the image. The calculation result is a 256×256 binary image, T 1 ∩ T 2 Indicates to calculate the intersection of binary images T 1 and T 2 , the specific method is to perform point-to-point AND operation on the corresponding position pixels, and the calculation result is also a 256×256 binary image. ||T|| 1 means to calculate its l 1 norm for the binary image T:
本发明的有益效果是:该方法首先采用基于图像滤波的方法提取模糊图像的互补性特征,通过对多帧模糊图像进行高频滤波,得到图像对应的高频细节信息;针对图像高频信息进行梯度卷积操作,分别提取水平方向以及竖直方向的图像梯度信息;对模糊图像的高频梯度图进行二值化操作,得到模糊图像的高频梯度二值图,该二值图即为模糊图像的互补性特征。然后基于提取出的模糊图像的互补性特征,利用互补性计算公式进行互补性定量计算,最终得到模糊图像间的互补性数值大小。该方法对不同模糊图像的互补性进行衡量,为模糊图像多帧复原提供了更加明确的图像选帧依据,对于图像质量参差不齐的序列多帧图像,基于互补性表征的图像选帧方式提高了图像复原的质量。The beneficial effect of the present invention is that: the method first adopts the method based on image filtering to extract the complementary features of the fuzzy image, and obtains the high-frequency detail information corresponding to the image by performing high-frequency filtering on multi-frame fuzzy images; The gradient convolution operation extracts the image gradient information in the horizontal direction and the vertical direction respectively; performs a binarization operation on the high-frequency gradient map of the blurred image to obtain a high-frequency gradient binary map of the blurred image, which is the blurred image. Complementary features of images. Then, based on the extracted complementarity features of the fuzzy images, the complementarity calculation formula is used to carry out quantitative calculation of complementarity, and finally the value of complementarity between fuzzy images is obtained. This method measures the complementarity of different blurred images, and provides a clearer basis for image frame selection for multi-frame restoration of blurred images. quality of image restoration.
下面结合具体实施方式对本发明作详细说明。The present invention will be described in detail below in combination with specific embodiments.
具体实施方式detailed description
本发明图像模糊互补性表征方法具体步骤如下:The specific steps of the image fuzzy complementarity characterization method of the present invention are as follows:
步骤一、提取模糊图像的互补性特征。Step 1, extracting the complementary features of the blurred image.
针对大小为256×256的序列图像中的任意两幅模糊图像g1和g2,采用高斯滤波器G进行图像滤波分解,并计算得到模糊图像的高频分量h1和h2:For any two blurred images g 1 and g 2 in the sequence image with a size of 256×256, the Gaussian filter G is used for image filtering and decomposition, and the high frequency components h 1 and h 2 of the blurred image are calculated:
其中,为卷积算子,二维高斯滤波卷积算子G如下所示:in, is a convolution operator, and the two-dimensional Gaussian filter convolution operator G is as follows:
分别在水平方向和竖直方向上对各高频分量图像hi做一阶梯度微分,即水平和竖直梯度算子卷积操作,得到水平方向梯度图像hdix和竖直方向梯度图像hdiy:Perform first-order gradient differentiation on each high-frequency component image h i in the horizontal direction and vertical direction, that is, the horizontal and vertical gradient operator convolution operation, to obtain the horizontal gradient image hd ix and the vertical gradient image hd iy :
其中,dx为水平梯度算子[1,-1],dy为竖直梯度算子[1,-1]T。在计算水平方向和竖直方向梯度图像的基础上,进一步计算提取图像的全局梯度图像hdi:Among them, d x is the horizontal gradient operator [1,-1], d y is the vertical gradient operator [1,-1] T . On the basis of calculating the gradient images in the horizontal direction and the vertical direction, further calculate the global gradient image hd i of the extracted image:
对计算得到的梯度图像hdi进行图像二值化操作,采用公式(6)提取模糊图像的互补性特征Ti(x,y)。Perform image binarization on the calculated gradient image hd i , and use formula (6) to extract the complementary feature T i (x, y) of the blurred image.
其中,二值化阈值ki为:Among them, the binarization threshold k i is:
式中,max(hdi)表示提取梯度图像hdi的最大像素灰度值,min(hdi)表示提取梯度图像hdi的最小像素灰度值。In the formula, max(hd i ) represents the maximum pixel gray value of the extracted gradient image hd i , and min(hd i ) represents the minimum pixel gray value of the extracted gradient image hd i .
步骤二、模糊图像互补性表征。Step 2: Fuzzy image complementarity characterization.
对于计算得到的二值图像T1和T2,以其作为模糊图像的互补性特征,基于数学集合论的思想,采用公式(8),对模糊图像g1和g2的互补性进行表征。For the calculated binary images T 1 and T 2 , as the complementarity features of blurred images, based on the idea of mathematical set theory, formula (8) is used to characterize the complementarity of blurred images g 1 and g 2 .
上述公式表示模糊图像g1和g2之间的互补性大小。其中,T1∪T2表示计算二值图像T1与T2的并集,具体方法是图像对应位置像素进行点对点“或”运算,计算结果是一个256×256的二值图像,T1∩T2表示计算二值图像T1与T2的交集,具体方法是对应位置像素进行点对点“与”运算,计算结果也是一个256×256的二值图像。||T||1表示对二值图像T计算其l1范数:The above formula expresses the magnitude of complementarity between blurred images g1 and g2. Among them, T 1 ∪ T 2 means to calculate the union of binary images T 1 and T 2. The specific method is to perform a point-to-point "OR" operation on the corresponding pixels of the image. The calculation result is a 256×256 binary image, T 1 ∩ T 2 means to calculate the intersection of binary images T 1 and T 2. The specific method is to perform point-to-point "AND" operation on the corresponding position pixels, and the calculation result is also a 256×256 binary image. ||T|| 1 means to calculate its l 1 norm for the binary image T:
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710558860.0A CN107481197A (en) | 2017-07-11 | 2017-07-11 | Image Fuzzy Complementary characterizing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710558860.0A CN107481197A (en) | 2017-07-11 | 2017-07-11 | Image Fuzzy Complementary characterizing method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107481197A true CN107481197A (en) | 2017-12-15 |
Family
ID=60594951
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710558860.0A Pending CN107481197A (en) | 2017-07-11 | 2017-07-11 | Image Fuzzy Complementary characterizing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107481197A (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867289A (en) * | 2012-08-20 | 2013-01-09 | 清华大学 | Complementary blurred image acquisition system and blurred image recovery method using complementary blurred image acquisition system |
US20130230244A1 (en) * | 2012-03-02 | 2013-09-05 | Chintan Intwala | Continuously Adjustable Bleed for Selected Region Blurring |
CN104282003A (en) * | 2014-10-29 | 2015-01-14 | 哈尔滨工业大学 | Digital blurred image blind restoration method based on gradient screening |
-
2017
- 2017-07-11 CN CN201710558860.0A patent/CN107481197A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130230244A1 (en) * | 2012-03-02 | 2013-09-05 | Chintan Intwala | Continuously Adjustable Bleed for Selected Region Blurring |
CN102867289A (en) * | 2012-08-20 | 2013-01-09 | 清华大学 | Complementary blurred image acquisition system and blurred image recovery method using complementary blurred image acquisition system |
CN104282003A (en) * | 2014-10-29 | 2015-01-14 | 哈尔滨工业大学 | Digital blurred image blind restoration method based on gradient screening |
Non-Patent Citations (1)
Title |
---|
王珮等: "基于互补性的运动模糊图像复原方法", 《中国科技论文在线》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103606132B (en) | Based on the multiframe Digital Image Noise method of spatial domain and time domain combined filtering | |
CN101639938B (en) | Image restoration method based on double-edge wave filter and margin deconvolution | |
CN103942758B (en) | Dark channel prior image dehazing method based on multiscale fusion | |
CN103426182B (en) | The electronic image stabilization method of view-based access control model attention mechanism | |
CN102446352B (en) | Method of video image processing and device | |
WO2018214671A1 (en) | Image distortion correction method and device and electronic device | |
CN104282003B (en) | Digital blurred image blind restoration method based on gradient screening | |
CN103996178A (en) | Sand and dust weather color image enhancing method | |
CN101996406A (en) | No-reference structure sharpness image quality assessment method | |
Shen et al. | Convolutional neural pyramid for image processing | |
CN102819837B (en) | Method and device for depth map processing based on feedback control | |
CN112163994B (en) | Multi-scale medical image fusion method based on convolutional neural network | |
CN107464217A (en) | A kind of image processing method and device | |
CN107146209A (en) | A Single Image Dehazing Method Based on Gradient Domain | |
CN108492245A (en) | Low light images based on wavelet decomposition and bilateral filtering are to fusion method | |
CN104504676A (en) | Full-reference image quality evaluation method based on multi-vision sensitive feature similarity | |
CN103886553A (en) | Method and system for non-local average value denoising of image | |
CN105184743A (en) | Image enhancement method based on non-linear guiding filtering | |
CN101504769B (en) | Self-adaptive noise intensity estimation method based on encoder frame work | |
CN108596843A (en) | A kind of underwater picture color recovery algorithm based on bright channel | |
Xu et al. | Image deblurring with blur kernel estimation in RGB channels | |
CN103313068B (en) | White balance corrected image processing method and device based on gray edge constraint gray world | |
CN104917933B (en) | The low-light video denoising method estimated based on phase equalization and scene motion | |
CN104754183B (en) | A kind of real-time monitor video adaptive filter method and its system | |
CN119784648A (en) | A multi-stage turbulence dynamic video restoration method based on physical model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171215 |