信息论在图像处理中的应用
信息论在图像处理中的应用
摘要:把信息论的基本原理应用到图像处理中具有十分重要的价值。本文主要从评估图像捕捉部分性能的评估、图像分割算法这两个个方面阐述信息论在图像处理中的应用。 通过理论分析来说明使用信息论的基本理论对图像处理的价值。
关键字:信息论;图像捕捉;图像分割
第1章 引言
随着科学技术的不断发展,人们对图形图像认识越来越广泛,图形图像处理的应用领域也将随之不断扩大。为了寻找快速有效的图像处理方法,信息理论越来越多地渗透到图像处理技术中。文章介绍了信息论基本理论在图像处理中的应用,并通过理论分析说明其价值。把通信系统的基本理论信息论应用于采样成像系统,对系统作端到端的系统性能评价,从而优化采样成像系统的设计,是当前采样成像系统研究的分支之一。有些图像很繁杂,而我们只需要其中有意义的一部分,图像分割就是将图像分为一些有意义的区域,然后对这些区域进行描述,就相当于提取出某些目标区域图像的特征,随后判断这些图像中是否有感兴趣的目标。
第2章 图像捕捉部分性能评估
2.1 图像捕捉的数学模型
图像捕捉过程如图1所示。G 为系统的稳态增益,p (x , y ) 是图像捕捉设备的空间响应
(x , y ) 代表采样网格函数,o (x , y ), s (x , y ) 分别为函数,n p (x , y ) 是光电探索的噪声。comb
输入、输出信号。
在这种模型下的输出信号
s (x , y ) =[Go (x , y ) *p (x , y )]comb (x , y ) +n p (x , y )
其中,comb (x , y ) =
设备的采样函数。
输出信号的傅立叶变换为:
S (u , v ) =GO (u , v ) P (u , v ) +N (u , v )
其中:O (u , v ) 是输入信号的傅立叶变换,N (u , v ) 是欠采样噪声和光电探测器噪声和,∑δ(x -m , y -n ) ,代表在直角坐标系下,具有单位采样间隔的采样m , n
P (u , v ) 是图像捕捉设备的空间频率响应。
图1 图像捕捉过程
2.2 性能评估
设信源X 通过系统后输出Y 。根据信息论知,X 与Y 之间的互信息量I (X , Y ) 定义为:
I (X , Y ) =H (Y ) -H (Y /X )
I (X , Y ) =H (X ) -H (X /Y )
H (X ), H (Y ) 分别为X,Y 的熵,H (X /Y ), H (Y /X ) 为条件熵。互信息量I (X , Y ) 的物理意义是输出Y 中得到的关于X 的平均信息量的大小。显然I (X , Y ) 越大,Y 得到的关于X 的越多。采样成像系统的目标是一致的。把输入图像O (u , v ) 看着信源,我们希望输出S (u , v ) 中包含多的关于O (u , v ) 的信息,即希望二者的互信息量尽量大。
根据互信息量的定义,o (x , y ), s (x , y ) 之间的互信息量为:
I =H [s (x , y )]-H [s (x , y ) /o (x , y )]=H [S (u , v )]-H [S (u , v ) /O (u , v )]
H [S (u , v ) /O (u , v )]可以称为噪声熵。因此可以用噪声熵H [N (u , v )]代替,则上式可以改写为:
I =H [S (u , v )]-H [N (u , v )]
又I 可改写为:
φo (u , v ) P (u , v ) 1]dudv (1) I =⎰⎰log[1+2G -2φp (u , v )
2,则有 φ0(u , v ) 是输入信号O (u , v ) 的功率谱。可令φp (u , v ) =σp
22 φn (u , v ) =G φo (u , v ) P (u , v ) *COMB (u , v ) +σp (2) 22
式(1),(2)是图像捕捉系统的基础。
第2章 图像分割算法
2.1 图像分割定义及方法概述
图像分割可借助集合的概念做如下的定义:令集合R 代表整个图像区域,对R 的分割可以看成将R 分成若干个满足以下条件的非空子集R 1, R 2, , Rn ;
(1) R
i =1n i =R (2) 对所有的i 和j ,i ≠j , 有R i R j =Θ;
(3) 对i =1, 2, 3, , n , 有P (R i )=TURE ; (4) 对 i ≠j , 有P R i R j =FALSE ;
(5) i =1, 2, 3, , n , R i 是连通区域
其中P (R i ) 是对所有在集合R i 中元素的逻辑谓词,Θ是空集。
对图像的分割可基于相邻像素在像素值方面的两个性质:不连续性和相似性。区域内部的像素一般具有某种相似性,而在区域之间的边界上一般具有某种不连续性。所以分割算法可据此分为利用区域间特性不连续性的基于边界的算法和利用区域内特性相似性的基于区域的算法。基于区域的算法又主要分为区域生长法和特征向量聚类法。另外根据分割过程中处理策略的不同,分割算法又可分为并行算法和串行算法。在并行算法中,所有判断和决定都可独立地和同时地做出,而在串行算法中,早期处理的结果可被其后的处理过程所利用。一般串行分割算法所需的计算时间常比并行分割算法要长,过程控制较复杂,但抗噪声能力也常较强。 ()
2.2 基于最大熵原理的图像分割
图像的信息熵反映了图像的总体概貌。若图像中包含目标,则在目标与背景可分割的交界处信息量(即熵) 最大。将Shannon 熵概念应用于图像分割时,依据是使图像中目标与背景分布的信息量最大,通过分析图像灰度直方图的熵,找到最佳阂值。
应用二值化进行图像分割,最终目的是通过对图像的处理把目标从背景中分割出来,即将图像分成两个区域:目标区和背景区,这样就得到了一幅二值图像。在一幅含有目标的多灰度图像中,必然存在一个灰度t ,以它作为阈值,可使图像得到最佳二值化分割。设t 将多灰度图像分成两个区域:一个区域的灰度值为0~t ,概率分布为F(t)=Pi/Pl ,其中i =0, 1, 2, , t , p l =∑p i 。该区域的熵为“一F(t)lnF(t)”;另一个区域的灰度值为
i =0t
t +1~L -1。概率分布为1一F(t),该区域的熵为“一(1一F(t))In(1一F(t))”。则总的熵为H(F(t))=一F(t)1nF(t)一(1一F(t))In(1一F(t))。根据信息论,当把目标从背景中最佳分割出来时,熵应最大。所以,使H(F(t))最大的t 即为最佳阈值。
设由最大类间方差法得到的分割阈值为t 1,则H(F(t 1))=一F(t 1)ln F(t 1) 一(1一F(t 1))ln(1一F(t 1) ;由一致性准则法得到的分割阈值为t 2,则H(F(t 2))=一F (t 2)lnF(t 2) 一(1一F(t 2))In(1一F(t 2) 。为了使分割后的二值图像同时具有最大类间方差和最大一致性,选择的阈值t 应满足min(t 1,t 2) ≤t ≤maX(,t 2) ,由于F 是t 的增函数,因此可以得到min(F(t 1) ,F(t 2))≤F(t)≤maX(F(t 1) ,F(t 2)) 。根据最大熵法,最佳阈值t 应为t =arg maxH(F(t)),即满足H(F(t))最大的t 的取值就是所要得到的最佳阈值t 。这就是基于最大熵原理的图像分割算法。 ***
第3章 结论
随着信息技术的迅猛发展,信息理论在通信领域中发挥了越来越重要的作用,由于信息理论解决问题的思路和方法独特、新颖和有效,信息论已经渗透到其他科学领域。随着计算机技术的发展,信息熵已经不仅仅在通信领域中使用。将信息熵应用到图像处理中,对图像处理技术的发展有着重大的作用,使得图像处理的算法更加优良。