CN107886485A - 基于Logistic隶属度函数的图像模糊增强方法及系统 - Google Patents
基于Logistic隶属度函数的图像模糊增强方法及系统 Download PDFInfo
- Publication number
- CN107886485A CN107886485A CN201711249581.2A CN201711249581A CN107886485A CN 107886485 A CN107886485 A CN 107886485A CN 201711249581 A CN201711249581 A CN 201711249581A CN 107886485 A CN107886485 A CN 107886485A
- Authority
- CN
- China
- Prior art keywords
- mrow
- image
- fuzzy
- msub
- msup
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000002708 enhancing effect Effects 0.000 claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 230000002146 bilateral effect Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 4
- 238000005728 strengthening Methods 0.000 claims description 4
- 230000007547 defect Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000033228 biological regulation Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000004321 preservation Methods 0.000 abstract 1
- 230000008859 change Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000013507 mapping Methods 0.000 description 4
- 230000001629 suppression Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000001931 thermography Methods 0.000 description 2
- 238000007630 basic procedure Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
- G06T2207/20028—Bilateral filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于Logistic隶属度函数的图像模糊增强方法,其包括步骤:S10:获取待增强图像,并将其表达为空间域上的待增强图像;S20:采用Logistic函数作为模糊隶属度函数,将所述空间域上的待增强图像从空间域映射至模糊域,得到模糊域上的待增强图像;S30:对模糊域上的待增强图像进行增强处理,得到模糊域上的增强图像;S40:对模糊域上的增强图像进行逆变换,得到空间域上的增强图像。此外,本发明还公开了相应的系统。本发明能简单有效地对待增强图像进行增强,在保存图像边缘信息完整性、调节的适度性、算法的复杂度上有相对的优势,可用于电力设备红外热像增强,以实现设备的故障或热缺陷的正确检出。
Description
技术领域
本发明涉及图像增强技术,尤其涉及一种图像模糊增强方法及系统。
背景技术
红外热像检测技术是保障电力设备安全的重要巡检手段之一。然而红外热成像技术往往受外界环境影响较大,存在设备图像边缘模糊、对比度差等问题,同时还存在各种各样的噪声,这往往会导致设备热缺陷的正确识别、分割和提取存在较大的困难。因此,有必要进行电力设备红外热像增强以实现设备的故障或热缺陷的正确检出。一般而言,电力设备红外热像增强方法可以分为两类,其一为基于空间域处理方法,即直接对灰度级进行处理,如Laxmikant Dash,Chatterji,B.N.首先提出利用亮度变化进行图像增强,该法相比较传统的增强方法更优,但受制于噪声,而且存在过增强和过抑制现象;其二为基于频域处理方法,即在变换域内改变图像参数的间接增强方法,如基于小波变换、金字塔变换等的图像增强方法,该类方法虽可有效抑制噪声,但算法的复杂程度较高,图像本身的不确定性因素也较多。
为解决上述问题,模糊理论逐渐应用到电力设备的红外图像增强处理之中:首先对给定图像选取适当的隶属度函数,通过该函数将原图像映射到对应的模糊集上,得到模糊特征平面,经模糊处理增强结果,再通过其反函数去模糊化输出结果。模糊增强理论常采用经典S型函数、Pal函数等作为隶属函数,但是在计算程度、去模糊化等方面效果不甚理想。Pal算法则将处理后的像素负值硬性规定为0,这样将导致被强行削弱为0值部分像素信息的缺失;经典S型函数与经典Pal函数相比,不存在过调或过抑制现象,而且计算复杂度较低,但是因其为分段函数,对其参数进行控制相对较难。
发明内容
本发明的目的之一是提供一种基于Logistic隶属度函数的图像模糊增强方法,其能简单有效地对待增强图像进行增强。
基于上述目的,本发明提供了一种基于Logistic隶属度函数的图像模糊增强方法,其包括步骤:
S10:获取待增强图像,并将其表达为空间域上的待增强图像;
S20:采用Logistic函数作为模糊隶属度函数,将所述空间域上的待增强图像从空间域映射至模糊域,得到模糊域上的待增强图像;
S30:对模糊域上的待增强图像进行增强处理,得到模糊域上的增强图像;
S40:对模糊域上的增强图像进行逆变换,得到空间域上的增强图像。
本发明所述的基于Logistic隶属度函数的图像模糊增强方法,其运用基于Logistic隶属度函数的模糊理论对待增强图像进行模糊化,并对所得作为所述模糊域上的待增强图像的模糊特征平面进行增强,然后再通过逆变换进行去模糊化处理,得到效果明显增强的图像。
Logistic函数为模糊理论隶属函数在定义域上为连续函数,逆变换不存在无解情况,作为模糊算子可以保证增强的模糊特征平面逆变到空间域时,不存在信息丢失。当自变量区域无穷时,该函数变化率趋于0,可以使值域在一个有限区间内,避免过增强或过抑制现象。同时,与经典S函数相比,该方法在参数控制上也易于实现。
此外,与经典pal函数、经典S型函数进行分析比较,可以发现该函数在保存图像边缘信息完整性、调节的适度性、算法的复杂度上有相对的优势,在增强故障部分和抑制缓冲区域具有可佳的效果,并且运行效率较高,符合图像处理的实际要求。
因此,相对于现有的常规隶属度函数等处理算法,本发明方法处理的图像边缘信息更为完整,图像对比度大幅提高,并且执行效率高,可批量处理,因而能简单有效地对待增强图像进行增强。
进一步地,本发明所述的基于Logistic隶属度函数的图像模糊增强方法中,所述步骤S10中采用双边滤波器对空间域上的待增强图像进行双边滤波。
上述方案中,双边滤波的作用是抑制设备红外图像中的噪声干扰,为后面的图像增强提供较高质量的红外图像原图。
更进一步地,上述基于Logistic隶属度函数的图像模糊增强方法中,所述双边滤波包括步骤:
计算距离权重G和强度权重H;
按照式F(i,j)=H×G计算双边滤波器的权重F(i,j);
基于所述双边滤波器的权重F(i,j)计算滤波响应。
更进一步地,上述基于Logistic隶属度函数的图像模糊增强方法中,所述距离权重G按照下式计算:
其中,和为核距离矩阵,w为双边滤波器的窗口半径,δd为空间域距离标准方差,e为自然指数。
更进一步地,上述基于Logistic隶属度函数的图像模糊增强方法中,所述强度权重H按照下式计算:
其中,x(i,j)为空间域上的坐标为(i,j)的待增强图像的像素的灰度值,δr为强度标准方差,e为自然指数。
更进一步地,上述基于Logistic隶属度函数的图像模糊增强方法中,按照下式计算滤波响应B(i,j):
其中,I(i,j)为以x(i,j)为窗口中心,w为窗口半径提取的图像,F(i,j)为双边滤波器的权重,imax、imin分别为坐标i最大值和坐标i最小值,jmax、jmin分别为坐标j最大值和坐标j最小值。
进一步地,本发明所述的基于Logistic隶属度函数的图像模糊增强方法中,所述步骤S20中通过以下Logistic隶属度函数将所述空间域上的待增强图像从空间域映射至模糊域:
其中,u(x)表示模糊域上的待增强图像;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值;x(i,j)表示空间域上的待增强图像的坐标为(i,j)的像素的灰度值;k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
进一步地,本发明所述的基于Logistic隶属度函数的图像模糊增强方法中,所述步骤S30中通过以下函数对模糊域上的待增强图像进行增强处理:
其中u(x)′表示模糊域上的增强图像;uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值,a为取值在区间[0,1]内的阈值。
上述方案中,通常热缺陷区域为感兴趣区域,因此取模糊域上的待增强图像中一像素值a作为阈值将模糊域上的待增强图像分为两个亮度区域,大于阈值部分进行增强,小于阈值部分则进行抑制,从而实现模糊域上的待增强图像的增强。
进一步地,本发明所述的基于Logistic隶属度函数的图像模糊增强方法中,所述步骤S40通过以下函数进行所述逆变换,对模糊域上的增强图像进行逆变换,得到空间域上的增强图像:
其中,x′(i,j)为空间域上的增强图像中坐标为(i,j)的像素的灰度值,uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值,k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
本发明的另一目的是提供一种基于Logistic隶属度函数的图像模糊增强系统,其能简单有效地对待增强图像进行增强。
基于上述目的,本发明提供了一种基于Logistic隶属度函数的图像模糊增强系统,其采用上述任一方法对待增强图像进行增强。
本发明所述的基于Logistic隶属度函数的图像模糊增强系统,由于其采用了本发明所述的方法,同样能简单有效地对待增强图像进行增强。具体原理前已描述,在此不再赘述。
所述系统可以是具有对应本发明方法的软件的计算机。
本发明所述的基于Logistic隶属度函数的图像模糊增强方法,其具有以下优点和有益效果:
1)Logistic函数为模糊理论隶属函数在定义域上为连续函数,逆变换不存在无解情况,作为模糊算子可以保证增强的模糊特征平面逆变到空间域时,不存在信息丢失。当自变量区域无穷时,该函数变化率趋于0,可以使值域在一个有限区间内,避免过增强或过抑制现象。同时,与经典S函数相比,该方法在参数控制上也易于实现。
2)与经典pal函数、经典S型函数进行了分析比较,可以发现该函数在保存图像边缘信息完整性、调节的适度性、算法的复杂度上有相对的优势,在增强故障部分和抑制缓冲区域具有可佳的效果,并且运行效率较高,符合图像处理的实际要求。
3)可用于电力设备红外热像增强,以实现设备的故障或热缺陷的正确检出。
4)相对于现有的常规隶属度函数等处理算法,本发明处理后的红外图像边缘信息更为完整、图像对比度大幅提高。
5)执行效率高,可实现电力设备红外热像的批量处理以及设备的故障诊断。
6)克服现有技术中存在的噪声缺陷和复杂程度缺陷。
本发明所述的基于Logistic隶属度函数的图像模糊增强系统,其同样具有上述优点和有益效果。
附图说明
图1为本发明所述的基于Logistic隶属度函数的图像模糊增强方法的基本流程示意图。
图2为本发明实例中空间域上的待增强图像。
图3为本发明实例中模糊域上的待增强图像。
图4为本发明实例中模糊域上的增强图像。
图5为本发明实例中空间域上的增强图像。
具体实施方式
下面结合说明书附图及实施例进一步说明本发明所述的技术方案。
图1示意了本发明所述的基于Logistic隶属度函数的图像模糊增强方法的基本流程。
如图1所示,本发明的基于Logistic隶属度函数的图像模糊增强方法包括步骤:
S10:获取待增强图像,并将其表达为空间域上的待增强图像。
在某些实施方式下,该步骤还采用双边滤波器对空间域上的待增强图像进行双边滤波。
其中,双边滤波包括步骤:
首先,按照下式计算距离权重G:
其中,和为核距离矩阵,w为双边滤波器的窗口半径,δd为空间域距离标准方差。
其次,按照下式计算强度权重H:
其中,x(i,j)为空间域上的坐标为(i,j)的待增强图像的像素的灰度值,δr为强度标准方差。
再次,按照式F(i,j)=H×G计算双边滤波器的权重F(i,j);
最后,基于双边滤波器的权重F(i,j)计算滤波响应:
按照下式计算滤波响应B(i,j):
其中,I(i,j)为以x(i,j)为窗口中心,w为窗口半径提取的图像,F(i,j)为双边滤波器的权重,imax、imin分别为坐标i最大值和坐标i最小值,jmax、jmin分别为坐标j最大值和坐标j最小值。
S20:采用Logistic函数作为模糊隶属度函数,将空间域上的待增强图像从空间域映射至模糊域,得到模糊域上的待增强图像。
在某些实施方式下,该步骤通过以下函数将空间域上的待增强图像从空间域映射至模糊域:
其中,u(x)表示模糊域上的待增强图像;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值;x(i,j)表示空间域上的待增强图像的坐标为(i,j)的像素的灰度值;k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
S30:对模糊域上的待增强图像进行增强处理,得到模糊域上的增强图像。
在某些实施方式下,该步骤通过以下函数对模糊域上的待增强图像进行增强处理:
其中u(x)′表示模糊域上的增强图像;uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值,a为取值在区间[0,1]内的阈值。
S40:对模糊域上的增强图像进行逆变换,得到空间域上的增强图像。
在某些实施方式下,该步骤通过以下函数进行逆变换,对模糊域上的增强图像进行逆变换,得到空间域上的增强图像:
其中,x′(i,j)为空间域上的增强图像中坐标为(i,j)的像素的灰度值,uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值,k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
下面通过一个具体实例进一步说明本发明,该实例采用上述方法对作为待增强图像的电压互感器的红外图像进行增强。
该实例系统是具有对应上述方法的软件的计算机。
该实例通过上述计算机运行软件实现以下步骤:
步骤10:获取电压互感器的红外图像,将其转换到灰度级下,得到如图2所示的N×M的像素矩阵x(i,j),从而将电压互感器的红外图像表达为空间域上的待增强图像,此外还采用双边滤波器对其进行双边滤波。
该步骤中,导入分辨率为240×320的电压互感器的红外图像并灰度化处理得到240×320的灰度化像素矩阵x(i,j),双边滤波器的窗口半径w取5,空间域距离标准方差δd取3,强度标准方差δr取0.1,得到核距离矩阵 距离权重强度权重 然后令x(i,j)为窗口中心,以w为半径提取该区域内原图像赋给I(i,j),滤波响应为:其中双边滤波器的权重F(i,j)=H(i,j)×G(i,j)。
步骤20:采用Logistic函数作为模糊隶属度函数,将空间域上的待增强图像x从空间域映射至模糊域,得到如图3所示的模糊域上的待增强图像u(x)。
该步骤通过以下函数将空间域上的待增强图像从空间域映射至模糊域:
其中,u(x)表示模糊域上的待增强图像;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值;x(i,j)表示空间域上的待增强图像的坐标为(i,j)的像素的灰度值;k为曲线的倾斜度,本实例中取12;L为曲线的最值,本实例中取1;a为取值在区间[0,1]内的阈值,本实例中取0.5;xmin、xmax分别为最小像素灰度值和最大像素灰度值。
步骤30:对模糊域上的待增强图像u(x)进行增强处理,得到如图4所示的模糊域上的增强图像u(x)′。
该步骤通过以下函数对模糊域上的待增强图像进行增强处理:
其中u(x)′表示模糊域上的增强图像;uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值,a为取值在区间[0,1]内的阈值,本实例中取0.5。
步骤40:对模糊域上的增强图像u(x)′进行逆变换,得到如图5所示的空间域上的增强图像x′(i,j)。
该步骤通过以下函数进行逆变换,对模糊域上的增强图像进行逆变换,得到空间域上的增强图像:
其中,x′(i,j)为空间域上的增强图像中坐标为(i,j)的像素的灰度值,uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值,k为曲线的倾斜度,本实例中取12;L为曲线的最值,本实例中取1;a为取值在区间[0,1]内的阈值,本实例中取0.5;xmin、xmax分别为最小像素灰度值和最大像素灰度值。
需要说明的是,本发明的保护范围中现有技术部分并不局限于本申请文件所给出的实施例,所有不与本发明的方案相矛盾的现有技术,包括但不局限于在先专利文献、在先公开出版物,在先公开使用等等,都可纳入本发明的保护范围。
另外,还需要说明的是,本案中各技术特征的组合方式并不限本案权利要求中所记载的组合方式或是具体实施例所记载的组合方式,本案所记载的所有技术特征可以以任何方式进行自由组合或结合,除非相互之间产生矛盾。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (10)
1.一种基于Logistic隶属度函数的图像模糊增强方法,其特征在于,包括步骤:
S10:获取待增强图像,并将其表达为空间域上的待增强图像;
S20:采用Logistic函数作为模糊隶属度函数,将所述空间域上的待增强图像从空间域映射至模糊域,得到模糊域上的待增强图像;
S30:对模糊域上的待增强图像进行增强处理,得到模糊域上的增强图像;
S40:对模糊域上的增强图像进行逆变换,得到空间域上的增强图像。
2.如权利要求1所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述步骤S10中采用双边滤波器对空间域上的待增强图像进行双边滤波。
3.如权利要求2所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述双边滤波包括步骤:
计算距离权重G和强度权重H;
按照式F(i,j)=H×G计算双边滤波器的权重F(i,j);
基于所述双边滤波器的权重F(i,j)计算滤波响应。
4.如权利要求3所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述距离权重G按照下式计算:
<mrow>
<mi>G</mi>
<mo>=</mo>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>Y</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<msub>
<mi>&delta;</mi>
<mi>d</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
其中,和为核距离矩阵,w为双边滤波器的窗口半径,δd为空间域距离标准方差。
5.如权利要求4所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述强度权重H按照下式计算:
<mrow>
<mi>H</mi>
<mo>=</mo>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<mi>x</mi>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msup>
<msub>
<mi>&delta;</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
其中,x(i,j)为空间域上的坐标为(i,j)的待增强图像的像素的灰度值,δr为强度标准方差。
6.如权利要求5所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,按照下式计算滤波响应B(i,j):
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mi>min</mi>
<mo>,</mo>
<mi>j</mi>
<mi>min</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>,</mo>
<mi>j</mi>
<mi>max</mi>
</mrow>
</msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mi>min</mi>
<mo>,</mo>
<mi>j</mi>
<mi>min</mi>
</mrow>
<mrow>
<mi>i</mi>
<mi>max</mi>
<mo>,</mo>
<mi>j</mi>
<mi>max</mi>
</mrow>
</msubsup>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,I(i,j)为以x(i,j)为窗口中心,w为窗口半径提取的图像,F(i,j)为双边滤波器的权重,imax、imin分别为坐标i最大值和坐标i最小值,jmax、jmin分别为坐标j最大值和坐标j最小值。
7.如权利要求1所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述步骤S20中通过以下函数将所述空间域上的待增强图像从空间域映射至模糊域:
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<msub>
<mi>i</mi>
<mi>j</mi>
</msub>
</msub>
<mo>=</mo>
<mfrac>
<mi>L</mi>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>x</mi>
<mi>max</mi>
</msub>
</mfrac>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
其中,y(x)表示模糊域上的待增强图像;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值;x(i,j)表示空间域上的待增强图像的坐标为(i,j)的像素的灰度值;k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
8.如权利要求1所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述步骤S30中通过以下函数对模糊域上的待增强图像进行增强处理:
<mrow>
<mi>u</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</msup>
<mo>=</mo>
<msup>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
</msup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mo>&times;</mo>
<msup>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo><</mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mo>&times;</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中u(x)′表示模糊域上的增强图像;uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值;uij表示模糊域上的待增强图像的坐标为(i,j)的像素的灰度值,a为取值在区间[0,1]内的阈值。
9.如权利要求1所述的基于Logistic隶属度函数的图像模糊增强方法,其特征在于,所述步骤S40通过以下函数进行所述逆变换,对模糊域上的增强图像进行逆变换,得到空间域上的增强图像:
<mrow>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<mfrac>
<mi>L</mi>
<mrow>
<msup>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中,x′(i,j)为空间域上的增强图像中坐标为(i,j)的像素的灰度值,uij′表示模糊域上的增强图像的坐标为(i,j)的像素的灰度值,k为曲线的倾斜度,L为曲线的最值,a为取值在区间[0,1]内的阈值,xmin、xmax分别为最小像素灰度值和最大像素灰度值。
10.一种基于Logistic隶属度函数的图像模糊增强系统,其采用如权利要求1-9中任意一项权利要求所述方法对待增强图像进行增强。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711249581.2A CN107886485A (zh) | 2017-12-01 | 2017-12-01 | 基于Logistic隶属度函数的图像模糊增强方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711249581.2A CN107886485A (zh) | 2017-12-01 | 2017-12-01 | 基于Logistic隶属度函数的图像模糊增强方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107886485A true CN107886485A (zh) | 2018-04-06 |
Family
ID=61776550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711249581.2A Pending CN107886485A (zh) | 2017-12-01 | 2017-12-01 | 基于Logistic隶属度函数的图像模糊增强方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107886485A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113643219A (zh) * | 2021-08-03 | 2021-11-12 | 武汉三江中电科技有限责任公司 | 一种基于三光融合的图像成像方法及装置 |
CN113763265A (zh) * | 2021-08-09 | 2021-12-07 | 云南北方光电仪器有限公司 | 红外图像对比度增益调节方法、存储介质及红外探测器 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020923A (zh) * | 2011-09-28 | 2013-04-03 | 中国航天科工集团第二研究院二0七所 | 模糊增强的光电搜索跟踪系统图像预处理方法 |
CN103632342A (zh) * | 2013-11-12 | 2014-03-12 | 华南理工大学 | 一种用于集成电路封装中的x射线图像模糊增强方法 |
CN103871029A (zh) * | 2014-01-28 | 2014-06-18 | 西安科技大学 | 一种图像增强及分割方法 |
CN106204470A (zh) * | 2016-07-01 | 2016-12-07 | 湖南源信光电科技有限公司 | 基于模糊理论的低照度成像方法 |
-
2017
- 2017-12-01 CN CN201711249581.2A patent/CN107886485A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020923A (zh) * | 2011-09-28 | 2013-04-03 | 中国航天科工集团第二研究院二0七所 | 模糊增强的光电搜索跟踪系统图像预处理方法 |
CN103632342A (zh) * | 2013-11-12 | 2014-03-12 | 华南理工大学 | 一种用于集成电路封装中的x射线图像模糊增强方法 |
CN103871029A (zh) * | 2014-01-28 | 2014-06-18 | 西安科技大学 | 一种图像增强及分割方法 |
CN106204470A (zh) * | 2016-07-01 | 2016-12-07 | 湖南源信光电科技有限公司 | 基于模糊理论的低照度成像方法 |
Non-Patent Citations (7)
Title |
---|
AKIRA TAGUCHI 等: "Image Enhancement of Noisy Images by Using Fuzzy Data-Dependent Bilateral Filter", 《2009 SIXTH INTERNATIONAL CONFERENCE ON FUZZY SYSTEMS AND KNOWLEDGE DISCOVERY》 * |
M. M. RIAZ 等: "PRINCIPLE COMPONENT ANALYSIS AND FUZZY LOGIC BASED THROUGH WALL IMAGE ENHANCEMENT", 《PROGRESS IN ELECTROMAGNETICS RESEARCH》 * |
余胜威 等编著: "《MATLAB图像滤波去噪分析及其应用》", 30 June 2015, 北京:北京航空航天大学出版社 * |
张玺: "基于稳态特性的柔性作业车间重调度决策方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑(月刊)》 * |
李红光 等: "基于双边滤波的变电站设备图像增强算法", 《电子世界》 * |
赵军 等: "基于Logistic统计的模糊法识别储层流体", 《地球物理学进展》 * |
赵永珠: "船舶焊缝图像缺陷识别系统研究", 《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113643219A (zh) * | 2021-08-03 | 2021-11-12 | 武汉三江中电科技有限责任公司 | 一种基于三光融合的图像成像方法及装置 |
CN113643219B (zh) * | 2021-08-03 | 2023-11-24 | 武汉三江中电科技有限责任公司 | 一种基于三光融合的图像成像方法及装置 |
CN113763265A (zh) * | 2021-08-09 | 2021-12-07 | 云南北方光电仪器有限公司 | 红外图像对比度增益调节方法、存储介质及红外探测器 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108932697B (zh) | 一种失真图像的去失真方法、装置及电子设备 | |
Li et al. | Adaptive fractional differential approach and its application to medical image enhancement | |
CN109151475A (zh) | 一种视频编码方法、解码方法、装置及电子设备 | |
CN103942758A (zh) | 基于多尺度融合的暗通道先验图像去雾方法 | |
WO2019001095A1 (zh) | 一种视频编码方法、解码方法、装置及电子设备 | |
CN108038833B (zh) | 一种梯度相关性检测的图像自适应锐化方法和存储介质 | |
Zhao et al. | Deep reconstruction of least significant bits for bit-depth expansion | |
CN107886486A (zh) | 基于暗通道先验与变分Retinex水下图像增强方法 | |
CN101207697B (zh) | 一种融合多模糊算子的数字图像模糊增强方法 | |
CN105046658A (zh) | 一种低照度图像处理方法和装置 | |
CN111105357B (zh) | 一种失真图像的去失真方法、装置及电子设备 | |
Liu et al. | Learning hadamard-product-propagation for image dehazing and beyond | |
CN101188671A (zh) | 视频图像的锐化方法 | |
Chen et al. | Saliency-directed image interpolation using particle swarm optimization | |
CN107516302A (zh) | 一种基于OpenCV的混合图像增强的方法 | |
CN107886485A (zh) | 基于Logistic隶属度函数的图像模糊增强方法及系统 | |
CN107146202B (zh) | 基于l0正则化和模糊核后处理的图像盲去模糊的方法 | |
CN117252773A (zh) | 基于自适应颜色校正和导向滤波的图像增强方法及系统 | |
Lei et al. | Low-light image enhancement using the cell vibration model | |
CN103700077A (zh) | 一种基于人眼视觉特性的自适应图像增强方法 | |
CN116681627B (zh) | 一种跨尺度融合的自适应水下图像生成对抗增强方法 | |
Yadav et al. | Underwater image enhancement using convolutional neural network | |
CN113763268A (zh) | 人脸图像盲修复方法及系统 | |
Hanumantharaju et al. | Adaptive color image enhancement based geometric mean filter | |
CN114565543A (zh) | 一种基于uv直方图特征的视频色彩增强方法及系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180406 |