[go: up one dir, main page]

CN107481197A - Image Fuzzy Complementary characterizing method - Google Patents

Image Fuzzy Complementary characterizing method Download PDF

Info

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
Application number
CN201710558860.0A
Other languages
Chinese (zh)
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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710558860.0A priority Critical patent/CN107481197A/en
Publication of CN107481197A publication Critical patent/CN107481197A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering 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

The invention discloses a kind of image Fuzzy Complementary characterizing method, for solving the ropy technical problem of conventional images complementarity characterizing method image restoration.Technical scheme is to extract the complementary feature of blurred picture first, by carrying out High frequency filter to multiframe blurred picture, obtains detail of the high frequency corresponding to image;Gradient convolution operation is carried out for image high-frequency information, extracts the gradient information of image;Binarization operation is carried out to the high frequency gradient map of blurred picture, obtains the high frequency gradient binary map of blurred picture.The complementary feature of blurred picture extracted is then based on, complementary quantitative calculating is carried out using complementary calculation formula, obtains the complementary numerical values recited between blurred picture.This method is weighed to the complementarity of different blurred pictures, the image provided definitely, which is restored, for blurred picture multiframe selects frame foundation, for the sequence multiple image that picture quality is uneven, quality that the image that is characterized based on complementarity selects frame mode to improve image restoration.

Description

图像模糊互补性表征方法Image Fuzzy Complementarity Representation Method

技术领域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和h2Step 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和竖直方向梯度图像hdiyPerform 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和h2For 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和竖直方向梯度图像hdiyPerform 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)

1.一种图像模糊互补性表征方法,其特征在于包括以下步骤:1. An image fuzzy complementarity characterization method is characterized in that comprising the following steps: 步骤一、针对大小为256×256的序列图像中的任意两幅模糊图像g1和g2,采用高斯滤波器G进行图像滤波分解,计算得到模糊图像的高频分量h1和h2Step 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: <mrow> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>G</mi> <mo>&amp;CircleTimes;</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>h</mi><mi>i</mi></msub><mo>=</mo><msub><mi>g</mi><mi>i</mi></msub><mo>-</mo><mi>G</mi><mo>&amp;CircleTimes;</mo><msub><mi>g</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow> 其中,为卷积算子,二维高斯滤波卷积算子G如公式(2)所示:in, is the convolution operator, and the two-dimensional Gaussian filter convolution operator G is shown in formula (2): <mrow> <mi>G</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>16</mn> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>2</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>2</mn> </mtd> <mtd> <mn>4</mn> </mtd> <mtd> <mn>2</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>2</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow><mi>G</mi><mo>=</mo><mfrac><mn>1</mn><mn>16</mn></mfrac><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>2</mn></mtd><mtd><mn>4</mn></mtd><mtd><mn>2</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow> 分别在水平方向和竖直方向上对各高频分量图像hi做一阶梯度微分,即水平和竖直梯度算子卷积操作,得到水平方向梯度图像hdix和竖直方向梯度图像hdiyPerform 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 : <mrow> <msub> <mi>hd</mi> <mrow> <mi>i</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>&amp;CircleTimes;</mo> <msub> <mi>d</mi> <mi>x</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>hd</mi><mrow><mi>i</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>h</mi><mi>i</mi></msub><mo>&amp;CircleTimes;</mo><msub><mi>d</mi><mi>x</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow> <mrow> <msub> <mi>hd</mi> <mrow> <mi>i</mi> <mi>y</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>h</mi> <mi>i</mi> </msub> <mo>&amp;CircleTimes;</mo> <msub> <mi>d</mi> <mi>y</mi> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>hd</mi><mrow><mi>i</mi><mi>y</mi></mrow></msub><mo>=</mo><msub><mi>h</mi><mi>i</mi></msub><mo>&amp;CircleTimes;</mo><msub><mi>d</mi><mi>y</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow> 其中,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 image in the horizontal direction and the vertical direction, further calculate and extract the image The global gradient image hd i : <mrow> <msub> <mi>hd</mi> <mi>i</mi> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <msub> <mi>hd</mi> <mrow> <mi>i</mi> <mi>x</mi> </mrow> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>hd</mi> <mrow> <mi>i</mi> <mi>y</mi> </mrow> </msub> <mn>2</mn> </msup> </mrow> </msqrt> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>hd</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><msup><msub><mi>hd</mi><mrow><mi>i</mi><mi>x</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>hd</mi><mrow><mi>i</mi><mi>y</mi></mrow></msub><mn>2</mn></msup></mrow></msqrt><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow> 对计算得到的梯度图像hdi进行图像二值化操作,采用公式(6)提取模糊图像的互补性特征Ti(x,y);Carry out image binarization operation on the calculated gradient image hd i , and use formula (6) to extract the complementary feature T i (x, y) of the blurred image; <mrow> <msub> <mi>T</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>1</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>hd</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <msub> <mi>k</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>0</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <msub> <mi>hd</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>&lt;</mo> <msub> <mi>k</mi> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>T</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>hd</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&amp;GreaterEqual;</mo><msub><mi>k</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>hd</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>k</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow> 其中,二值化阈值ki为:Among them, the binarization threshold k i is: <mrow> <msub> <mi>k</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>(</mo> <mrow> <msub> <mi>hd</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>-</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mo>(</mo> <mrow> <msub> <mi>hd</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> <mrow><msub><mi>k</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mi>m</mi><mi>a</mi><mi>x</mi><mo>(</mo><mrow><msub><mi>hd</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>-</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>(</mo><mrow><msub><mi>hd</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow> 式中,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; <mrow> <mi>c</mi> <mi>o</mi> <mi>m</mi> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>g</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>&amp;cup;</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>&amp;cap;</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>|</mo> <msub> <mo>|</mo> <mn>1</mn> </msub> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>&amp;cup;</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>|</mo> <msub> <mo>|</mo> <mn>1</mn> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow><mi>c</mi><mi>o</mi><mi>m</mi><mrow><mo>(</mo><msub><mi>g</mi><mn>1</mn></msub><mo>,</mo><msub><mi>g</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msub><mi>T</mi><mn>1</mn></msub><mo>&amp;cup;</mo><msub><mi>T</mi><mn>2</mn></msub><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub><mo>&amp;cap;</mo><msub><mi>T</mi><mn>2</mn></msub><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub></mrow><mrow><mo>|</mo>mo><mo>|</mo><msub><mi>T</mi><mn>1</mn></msub><mo>&amp;cup;</mo><msub><mi>T</mi><mn>2</mn></msub><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow> 公式(8)表示模糊图像g1和g2之间的互补性大小;其中,T1∪T2表示计算二值图像T1与T2的并集,具体方法是图像对应位置像素进行点对点或运算,计算结果是一个256×256的二值图像,T1∩T2表示计算二值图像T1与T2的交集,具体方法是对应位置像素进行点对点与运算,计算结果也是一个256×256的二值图像;||T||1表示对二值图像T计算其l1范数:Formula (8) represents the complementarity between blurred images g 1 and g 2 ; 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, the calculation result is a 256×256 binary image, T 1 ∩ T 2 means to calculate the intersection of the binary image 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: <mrow> <mo>|</mo> <mo>|</mo> <mi>T</mi> <mo>|</mo> <msub> <mo>|</mo> <mn>1</mn> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>x</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>y</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <msub> <mi>t</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow> 1 <mrow><mo>|</mo><mo>|</mo><mi>T</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>=</mo><munderover><mo>&amp;Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mo>&amp;Sigma;</mo><mrow><mi>y</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>|</mo><msub><mi>t</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow><mo>.</mo></mrow> 1
CN201710558860.0A 2017-07-11 2017-07-11 Image Fuzzy Complementary characterizing method Pending CN107481197A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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