第3章数字图像处理.pdf

《第3章数字图像处理.pdf》由会员分享，可在线阅读，更多相关《第3章数字图像处理.pdf（62页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章图像增强 图像增强是采用一系列技术去改善图像的视觉效果, 或将图像转换成一种更适合于人或机器进行分析和处理 的形式。例如采用一系列技术有选择地突出某些感兴趣 的信息，同时抑制一些不需要的信息，提高图像的使用 价值。 图像增强方法从增强的作用域出发，可分为空间域 增强和频率域增强两种。 空间域增强是直接对图像各像素进行处理； 频率域增强是对图像经傅立叶变换后的频谱成分进 行处理，然后逆傅立叶变换获得所需的图像。 讲解内容 图像的代数运算 彩色变换及应用 伪彩色增强 假彩色增强 彩色增强 同态滤波增强 低通滤波 高通滤波 频率域 图像锐化 图像平滑 局部运算 局部统计法 规定化 均衡化 直方

2、图修正法 灰度变换 点运算 空间域 图像增强 目的 1. 熟悉并掌握本章基本概 念、空间域图像增强的原 理、方法及其特点； 2. 了解频率域图像增强的 方法及其实现过程； 3.重点掌握直方图修正方 法、特点及其应用；空间 域平滑、锐化和彩色增强 技术。 4.1图像增强的点运算 4.1.2 灰度变换 灰度变换可调整图像的灰度动态范围或图像对比度，是图 像增强的重要手段之一。 ) 51 . 4 (), (), ( +=ajif ab ab ajig 黑白 1线性变换 令图像f(i，j)的灰度范围为a,b， 线性变换后图像g(i，j)的范围为 a,b，如图，g(i，j)与f(i，j)之间 的关系式为

3、： 在曝光不足或过度的情况下，图像灰度可能会局限 在一个很小的范围内。这时在显示器上看到的将是一个 模糊不清、似乎没有灰度层次的图像。 下图是对曝光不足的图像采用线性变换对图像每一 个像素灰度作线性拉伸。可有效地改善图像视觉效果。 2分段线性变换 为了突出感兴趣目标所在 的灰度区间，相对抑制那些不 感兴趣的灰度区间，可采用分 段线性变换。 设原图像f(x,y)在0，Mf,感 兴趣目标的灰度范围在a,b,欲 使其灰度范围拉伸到c,d,则对 应的分段线性变换表达式为 + + = ffg MyxfbdbyxfbMdM byxfacayxfabcd ayxfyxfac yxg ),(),()/()(

4、),(),()/()( ),(0),()/( ),( 通过细心调整折线拐点的位置及控制分段直线的 斜率，可对任一灰度区间进行拉伸或压缩。 3非线性灰度变换 当用某些非线性函数如对数函数、指数函数等，作为 映射函数时，可实现图像灰度的非线性变换。 对数变换 对数变换的一般表达式为 ) 71 . 4 ( ln 1), (ln ), ( + += cb jif ajig 这里a,b,c是为了调整曲线 的位置和形状而引入的参数。 当希望对图像的低灰度区较 大的拉伸而对高灰度区压缩 时，可采用这种变换，它能 使图像灰度分布与人的视觉 特性相匹配。 f (i,j) g(i,j) 指数变换 指数变换的一般表

5、达式为 这里参数a,b,c用来调整曲线的位置和形状。这种变 换能对图像的高灰度区给予较大的拉伸。 ) 81 . 4 (1), ( ), ( = ajifc bjig g (i,j) f (i,j) 4.1.3 直方图修整法 灰度直方图反映了数字图像中每一灰度级与其出现频率 间的关系,它能描述该图像的概貌。通过修改直方图的方法 增强图像是一种实用而有效的处理技术。 直方图修整法包括直方图均衡化直方图均衡化及直方图规定化直方图规定化两类。 1.直方图均衡化直方图均衡化 直方图均衡化是将原图像通过某种变换，得到一幅灰度 直方图为均匀分布的新图像的方法。 直方图均衡化 下面先讨论连续变化图像的均衡化问

6、题，然后推广到 离散的数字图像上。 设r和s分别表示归一化了的原图像灰度和经直方图修 正后的图像灰度。即 （4.1-9） 在0,1区间内的任一个r值，都可产生一个s值，且 (4.1-10) 1,0sr )(rTs = T(r)作为变换函数，满足下列条件： 在0r1内为单调递增函数，保证灰度级从黑到白 的次序不变； 在0r1内，有0T(r)1，确保映射后的像素灰 度在允许的范围内。 反变换关系为 （4.1-11） T-1(s)对s同样满足上述两个条件。 由概率论理论可知，如果已知随机变量r的概率密度为 pr(r),而随机变量s是r的函数，则s的概率密度ps(s)可以由 pr(r)求出。 假定随机

7、变量s的分布函数用Fs(s)表示，根据分布函数 定义 )( 1 sTr = = = r r s sS drrpdsspsF)121 . 4()()()( 利用密度函数是分布函数的导数的关系，等式两边对s 求导，有： （4.1-13） 可见，输出图像的概率密度函数可以通过变换函数T(r) 控制原图像灰度级的概率密度函数得到,因而改善原图像的 灰度层次,这就是直方图修改技术的基础。 从人眼视觉特性来考虑，一幅图像的直方图如果是均匀 分布的，即Ps(s)=k(归一化时k=1)时，该图像色调给人的感 觉比较协调。因此将原图像直方图通过T(r)调整为均匀分布 的直方图，这样修正后的图像能满足人眼视觉要求

8、。 因为归一化假定 由（4.1-13）则有 )141 . 4(1)(=sP s drrpds r )(= )()()( 1 sT ds d p ds dr pdrrp ds d sP rr r rs = = 两边积分得 上式表明，当变换函数为r的累积直方图函数时，能达 到直方图均衡化的目的。 对于离散的数字图像，用频率来代替概率,则变换函数 T(rk)的离散形式可表示为: 上式表明，均衡后各像素的灰度值sk可直接由原图像的 直方图算出。 )151 . 4()()( 0 = r r drrprTs = = k j j k j jrkk n n rprTs 00 )()( 一幅图像的sk与rk之间

9、的关系称为该图像的累积灰 度直方图。 rk Pr(rk) rk S(rk) 1.0 1.0 1.0 下面举例说明直方图均衡过程。 rknkpr(rk)=nk/nsk计sk并sknskpk(s) r0=07900.190.191/7s0=1/77900.19 r1=1/710230.250.443/7s1=3/710230.25 r2=2/78500.210.655/7s2=5/78500.21 r3=3/76560.160.816/7 r4=4/73290.080.896/7s3=6/79850.24 r5=5/72450.060.951 r6=6/71220.030.981 r7=1810.

10、021.001s4=14480.11 例例假定有一幅总像素为n=6464的图像，灰度级数为8，各灰 度级分布列于表中。对其均衡化计算过程如下： ？ 若在原图像一行上连续8个像素的灰度值分别为：0、1、 2、3、4、5、6、7，则均衡后，他们的灰度值为多少？ 原图像的直方图均衡后图像的直方图 直方图均衡化示例 2.2.直方图规定化直方图规定化 在某些情况下，并不一定需要具有均匀直方图的图像， 有时需要具有特定的直方图的图像，以便能够增强图像中某 些灰度级。直方图规定化方法就是针对上述思想提出来的。 直方图规定化是使原图像灰度直方图变成规定形状的直方图 而对图像作修正的增强方法。 可见，它是对直方

11、图均衡化处理的一种有效的扩展。直 方图均衡化处理是直方图规定化的一个特例。 对于直方图规定化，下面仍从灰度连续变化的概率密度 函数出发进行推导，然后推广出灰度离散的图像直方图规定 化算法。 假设pr(r)和pz(z)分别表示已归一化的原始图像灰度分 布的概率密度函数和希望得到的图像的概率密度函数。 首先对原始图像进行直方图均衡化，即求变换函数： 假定已得到了所希望的图像，对它也进行均衡化处理，即 它的逆变换是 这表明可由均衡化后的灰度得到希望图像的灰度。 若对原始图像和希望图像都作了均衡化处理，则二者 均衡化的ps(s)和pv(v)相同，即都为均匀分布的密度函数。 由s代替v 得z=G-1(s

12、) = r r drrprTs 0 )171.4()()( )181.4()()( 0 = z z drrpzGv )191. 4()( 1 = vGz 这就是所求得的变换表达式。根据上述思想，可总 结出直方图规定化增强处理的步骤如下： 对原始图像作直方图均衡化处理； 按照希望得到的图像的灰度概率密度函数pz(z)，求得 变换函数G(z)； 用步骤得到的灰度级s作逆变换z= G-1(s)。 经过以上处理得到的图像的灰度级将具有规定的概 率密度函数pz(z)。 采用与直方图均衡相同的原始图像数据（6464像 素且具有8级灰度），其灰度级分布列于表中。给定的 直方图的灰度分布列于表中。 对应的直方

13、图如下： 原图像的直方图规定化直方图 rjsknkps(sk)zkpz(zk)vkzk并nkpz(zk) r0s0=1/77900.19z0=00.000.00z000.00 r1s1=3/710230.25z1=1/70.000.00z100.00 r2s2=5/78500.21z2=2/70.000.00z200.00 r3s3=6/7z3=3/70.150.15z3s0=1/77900.19 r4s3=6/79850.24z4=4/70.200.35z4s1=3/71023 0.25 r5s4=1z5=5/70.300.65z5s2=5/78500.21 r6s4=1z6=6/70.20

14、0.85z6s3=6/79850.24 r7s4=14480.1110.151.00z7s4=14480.11 117/67/3 17/67/57/2 17/57/47/ 1 7/67/47/30 7763 7652 7541 6430 = = = = zrzr zrzr zrzr zrzr 原图像的直方图规定的直方图规定化后图像的直方图 ？ 若在原图像一行上连续8个像素的灰度值分别为：0、1、 2、3、4、5、6、7，则规定化后，他们的灰度值为多少？ 利用直方图规定化方法进行图像增强的主要困难在于 要构成有意义的直方图。图像经直方图规定化，其增强效 果要有利于人的视觉判读或便于机器识别。 下

15、面是一个直方图规定化应用实例。 图(C)、(c)是将图像(A)按图(b)的直方图进行规定化得 到的结果及其直方图。通过对比可以看出图(C)的对比度同 图（B)接近一致，对应的直方图形状差异也不大。这样有利 于影像融合处理，保证融合影像光谱特性变化小。 111 111 111 9 1 模 板 4.2 图像的空间域平滑 任何一幅原始图像，在其获取和传输等过程中，会受 到各种噪声的干扰，使图像恶化，质量下降，图像模糊， 特征淹没，对图像分析不利。 为了抑制噪声改善图像质量所进行的处理称图像平滑 或去噪。它可以在空间域和频率域中进行。本节介绍空间 域的几种平滑法。 4.2.14.2.1局部平滑法局部平

16、滑法 局部平滑法是一种直接在空间域上进行平滑处理的技 术。假设图像是由许多灰度恒定的小块组成，相邻像素间 存在很高的空间相关性，而噪声则是统计独立的。因此， 可用邻域内各像素的灰度平均值代替该像素原来的灰度值， 实现图像的平滑。 设有一幅NN的图像f(x,y)，若平滑图像为 g(x,y),则有 式中x,y=0,1,N-1； s为（x,y）邻域内像素坐标的集合； M表示集合s内像素的总数。 可见邻域平均法就是将当前像素邻域内各像 素的灰度平均值作为其输出值的去噪方法。 ) 12 . 4(),( 1 ),( , = sji jif M yxg (m-1,n-1) (m-1,n)(m-1,n+1)

17、(m,n-1)(m,n)(m,n+1) (m+1,n-1)(m+1,n)(m+1,n+1) 例如，对图像采用33的邻域平均法，对于像素 （m,n)，其邻域像素如下： 则有： ),(),( 9 1 jnimfnmg ZiZj += 其作用相当于用这样的模板同图像卷积。 设图像中的噪声是随机不相关的加性噪声，窗口内 各点噪声是独立同分布的，经过上述平滑后，信号与噪 声的方差比可望提高M倍。 这种算法简单，但它的主要缺点是在降低噪声的同 时使图像产生模糊，特别在边缘和细节处。而且邻域越 大，在去噪能力增强的同时模糊程度越严重。如图 4.2.1(c)和(d)。 = 111 111 111 9 1 H

18、(a)原图像(b) 对(a)加椒盐噪声的图像 (c)33邻域平滑(d) 55邻域平滑 为克服简单局部平均法的弊病，目前已提出许多保 边缘、细节的局部平滑算法。它们的出发点都集中在如 何选择邻域的大小、形状和方向、参加平均的点数以及 邻域各点的权重系数等，下面简要介绍几种算法。 4.2.2 4.2.2 超限像素平滑法超限像素平滑法 对邻域平均法稍加改进，可导出超限像素平滑法。它 是将f(x,y)和邻域平均g(x,y)差的绝对值与选定的阈值进 行比较，根据比较结果决定点（x,y）的最后灰度g(x,y)。 其表达式为 这算法对抑制椒盐噪声比较有效，对保护仅有微小灰 度差的细节及纹理也有效。可见随着邻

19、域增大，去噪能力 增强，但模糊程度也大。 同局部平滑法相比，超限像元平 滑法去椒盐噪声效果更好。 (a)原图像 (b)对(a)加椒盐噪声的图像 (c)33邻域平滑 (d) 55邻域平滑 (e)33超限像素平滑(T=64) (f)55超限像素平滑(T=48) 4.2.3 4.2.3 灰度最相近的灰度最相近的K K个邻点平均法个邻点平均法 该算法的出发点是：在nn的窗口内，属于同一集合 体的像素，它们的灰度值将高度相关。因此，可用窗口内 与中心像素的灰度最接近的K个邻像素的平均灰度来代替 窗口中心像素的灰度值。这就是灰度最相近的K个邻点平 均法。 较小的K值使噪声方差下降较小，但保持细节效果较 好

20、；而较大的K值平滑噪声较好，但会使图像边缘模糊。 实验证明，对于33的窗口，取K=6为宜。 4.2.4 4.2.4 最大均匀性平滑最大均匀性平滑 为避免消除噪声引起边缘模糊，该 算法先找出环绕图像中每像素的最均匀 区域，然后用这区域的灰度均值代替该 像素原来的灰度值。 4.2.5 4.2.5 有选择保边缘平滑法有选择保边缘平滑法 该方法对图像上任一像素 (x,y)的55邻域，采用9个掩模， 其中包括一个33正方形、4个 五边形和4个六边形。计算各个 掩模的均值和方差，对方差进行 排序，最小方差所对应的掩模区 的灰度均值就是像素（x,y) 的 输出值。 该方法以方差作为各个区域灰度均匀性的测度。

21、若区域 含有尖锐的边缘，它的灰度方差必定很大，而不含边缘或灰 度均匀的区域，它的方差就小，那么最小方差所对应的区域 就是灰度最均匀区域。因此有选择保边缘平滑法既能够消除 噪声，又不破坏区域边界的细节。另外，五边形和六边形在 （x,y）处都有锐角，这样，即使像素（x,y）位于一个复杂 形状区域的锐角处，也能找到均匀的区域。从而在平滑时既 不会使尖锐边缘模糊，也不会破坏边缘形状。 例如，某像素55邻域的灰度分布如图4.2.4，经计算9 个掩模区的均值和方差为 最小方差为0，对应的灰度均值3， 采用有选择保边缘平滑，该像素的输出值为3。 4.2.6 4.2.6 空间低通滤波法空间低通滤波法 邻域平均

22、法可看作一个掩模作用于图像f(x,y)的低通 空间滤波，掩模就是一个滤波器，它的响应为H(r,s)， 于是滤波输出的数字图像g(x,y)用离散卷积表示为 均值443234233 对应的 方差 54717 17 28 31 23 260 36421 47324 84143 42153 43216 )62 . 4(),(),(),(= = l ls k kr srHsyrxfyxg 常用的掩模有 掩模不同，中心点或邻域的重要程度也不相同，因此， 应根据问题的需要选取合适的掩模。但不管什么样的掩模， 必须保证全部权系数之和为单位值，这样可保证输出图像灰 度值在许可范围内，不会产生“溢出”现象。 =

23、111 111 111 9 1 1 H = 111 121 111 10 1 2 H = 121 242 121 16 1 3 H = 111 101 111 8 1 4 H = 00 1 00 2 1 4 1 4 1 4 1 4 1 5 H 4.2.7 4.2.7 中值滤波中值滤波 中值滤波是对一个滑动窗口内的诸像素灰度值排序，用 中值代替窗口中心像素的原来灰度值，因此它是一种非线性 的图像平滑法。 例：采用13窗口进行中值滤波 原图像为：2 2 6 2 1 2 4 4 4 2 4 处理后为：2 2 2 2 2 2 4 4 4 4 4 它对脉冲干扰及椒盐噪声的抑制效果好，在抑制随机噪 声的同

24、时能有效保护边缘少受模糊。但它对点、线等细节较 多的图像却不太合适。 对中值滤波法来说，正确选择窗口尺寸的大小是很重要 的环节。一般很难事先确定最佳的窗口尺寸，需通过从小窗 口到大窗口的中值滤波试验，再从中选取最佳的。 原图像中值滤波 一维中值滤波的几个例子（N=5） 离散阶跃信号、斜升信号没有受到影响。离散三角信号 的顶部则变平了。对于离散的脉冲信号，当其连续出现的次 数小于窗口尺寸的一半时，将被抑制掉，否则将不受影响。 一维中值滤波的概念很容易推广到二维。一般来说，二 维中值滤波器比一维滤波器更能抑制噪声。 二维中值滤波器的窗口形状可以有多种，如线状、方 形、十字形、圆形、菱形等（见图）。

25、 不同形状的窗口产生不同的滤波效果，使用中必须根据 图像的内容和不同的要求加以选择。从以往的经验看，方形 或圆形窗口适宜于外轮廓线较长的物体图像，而十字形窗口 对有尖顶角状的图像效果好。 图(a)为原图像；图(b)为加椒盐噪声的图像；图(c)和图 (d)分别为33、55模板进行中值滤波的结果。 可见中值滤波法能有效削弱椒盐噪声，且比邻域、超限 像素平均法更有效。 111 191 111 4.3 图像空间域锐化 在图像的识别中常需要突出边缘和轮廓信息。图像锐化 就是增强图像的边缘或轮廓。 图像平滑通过积分过程使得图像边缘模糊，图像锐化则 通过微分而使图像边缘突出、清晰。 4.3.1 梯度锐化法梯

26、度锐化法 图像锐化法最常用的是梯度法。 对于图像f(x，y)，在(x， y)处的梯度定义为 梯度是一个矢量，其大小和方向为 )13.4(),( ),( ),( = = y yxf x yxf y x f f yxgrad )23 . 4()/()/( )()(y)grad(x, ),(),(11 2),(2),( 2 2 = +=+= x yxf y yxf xy y yxf x yxf yx tgfftg ff 对于离散图像处理而言，常用到梯度的大小，因此把 梯度的大小习惯称为“梯度”。并且一阶偏导数采用一阶 差分近似表示，即 fx=f(x +1，y)-f(x，y) fy=f(x，y +1)

27、-f(x，y) 为简化梯度的计算，经常使用 grad(x，y)=Max(|fx|，|fy|)(4.3-4) 或grad（x，y）=|fx|+|f y| (4.3-5) 除梯度算子以外，还可采用Roberts、Prewitt和 Sobel 算子计算梯度，来增强边缘。 Roberts对应的模板如图4.3.2所示。差分计算式如下 fx=|f(x+1,y+1)-f(x,y)| fy=|f(x+1,y)-f(x,y+1)| 为在锐化边缘的同时减少噪声的影响，Prewitt从加大 边缘增强算子的模板大小出发，由2x2扩大到3x3来计算差 分，如图(a)所示。 (a)Prewitt 算子(b)Sobel算子

28、 Sobel在Prewitt算子的基础上，对4-邻域采用带权的 方法计算差分，对应的模板如图(b)。 根据梯度计算式就可以计算Roberts、Prewitt和Sobel 梯度。一旦梯度算出后，就可根据不同的需要生成不同的 梯度增强图像。 -101-1-1-1-101-1-2-1 -101000-202000 -101111-101121 第一种输出形式 g(x,y)=grad(x,y)(4.3-7) 此法的缺点是增强的图像仅显示灰度变化比较徒的边缘轮 廓，而灰度变化比较平缓或均匀的区域则呈黑色。 第二种输出形式 式中T是一个非负的阈值。适当选取T，可使明显的边缘 轮廓得到突出，又不会破坏原来灰

29、度变化比较平缓的背景 第三种输出形式 它将明显边缘用一固定的灰度级LG来表现。 = 其它),( ),(),( ),( yxf Tyxgradyxgrad yxg = 其他， ， ),( ),( ),( yxf TyxgradL yxg G 第四种输出形式 此方法将背景用一个固定的灰度级 LB来表现，便于研究 边缘灰度的变化。 第五种输出形式 这种方法将明显边缘和背景分别用灰度级LG和LB表示， 生成二值图像，便于研究边缘所在位置。 = 其他, ),(,),( ),( B L Tyxgradyxgrad yxg = 其他, ),(, ),( B G L TyxgradL yxg 4.3.2 La

30、placian增强算子增强算子 Laplacian 算子是线性二阶微分算子。即 2f(x，y)= 2 2 2 2 ),(),( y yxf x yxf + 对离散的数字图像而言，二阶偏导数可 用二阶差分近似，可推导出Laplacian算子 表达式为 2f(x，y)= f(x+1，y)+f(x-1，y)+ f(x，y+1)+f(x，y-1)-4f(x，y) Laplacian增强算子为： g(x，y)=f(x，y)- 2f(x，y) =5f(x，y)- f(x+1，y)+ f(x- 1，y)+f(x，y+1)+ f(x，y-1) Laplacian增强算子 Laplacian算子 其特点是： 1

31、、在灰度均匀的区域或斜坡中间2f(x，y)为0，增强 图像上像元灰度不变； 2、在斜坡底或低灰度侧形成“下冲”；而在斜坡顶或高 灰度侧形成“上冲”。 0 -10-1 1 1 H1= -15 1H2= -19 1 0 -10-1 1 1 4.3.3 4.3.3 高通滤波法高通滤波法 高通滤波法就是用高通滤波算子和图像卷积来增强边 缘。常用的算子有： 4.44.4图像的频率域增强图像的频率域增强 图像增强的目的主要包括：消除噪声，改善图像的视 觉效果；突出边缘，有利于识别和处理。前面是关于图像 空间域增强的知识，下面介绍频率域增强的方法。 假定原图像为f(x，y)，经傅立叶变换为F(u，v)。频率

32、域 增强就是选择合适的滤波器H(u,v)对F(u,v)的频谱成分进行 处理，然后经逆傅立叶变换得到增强的图像g(x,y)。 频率域增强的一般过程如下： DFT H(u，v) IDFT f(x，y) F(u，v) F(u，v)H(u，v) g(x，y) 滤波 图像的平滑除了在空间域中进行外，也可以在频率域中 进行。由于噪声主要集中在高频部分，为去除噪声改善图像 质量，滤波器采用低通滤波器H(u,v)来抑制高频成分，通过 低频成分，然后再进行逆傅立叶变换获得滤波图像，就可达 到平滑图像的目的。常用的频率域低滤波器H(u,v)有四种： 1理想低通滤波器 设傅立叶平面上理想低通滤波器离开原点的截止频率

33、为 D0，则理想低通滤波器的传递函数为 由于高频成分包含 有大量的边缘信息，因 此采用该滤波器在去噪 声的同时将会导致边缘信息损失而使图像边模糊。 4.4.1频率域平滑频率域平滑 ) 14 . 4 ( ),(0 ),(1 ),( 0 0 = DvuD DvuD vuH 2Butterworth低通滤波器 n阶Butterworth滤波器的传递函数为： 它的特性是连续性衰减，而不象理想滤波器那样陡峭 变化，即明显的不连续性。因此采用该滤波器滤波在抑制 噪声的同时，图像边缘的模糊程度大大减小，没有振铃效 应产生。 )24 . 4(),( 2 0 ),( 1 1 = + n D vuD vuH 3指

34、数低通滤波器 指数低通滤波器是图像处理中常用的另一种平滑滤波器。 它的传递函数为： 采用该滤波器滤波在抑制噪声的同时，图像边缘的模糊程 度较用Butterworth滤波产生的大些，无明显的振铃效应。 )34 . 4(e v)H(u, - 0 D v)D(u, = n 4. 梯形低通滤波器 梯形低通滤波器是理想低通滤波器和完全平滑滤波器的 折中。它的传递函数为： 它的性能介于理想低通滤波器和指数滤波器之间，滤波 的图像有一定的模糊和振铃效应。 )44 . 4( Dv)D(u,0 D),(D Dv)D(u,1 v)H(u, 1 10DD D-v)D(u, 0 10 1 = vuD 4.4.2 4.

35、4.2 频率域锐化频率域锐化 图像的边缘、细节主要位于高频部分，而图像的模糊 是由于高频成分比较弱产生的。频率域锐化频率域锐化就是为了消除 模糊，突出边缘。因此采用高通滤波器让高频成分通过， 使低频成分削弱，再经逆傅立叶变换得到边缘锐化的图像 。常用的高通滤波器有： 1）理想高通滤波器 二维理想高通滤波器的传递函数为 )54 . 4( ),(1 ),(0 ),( 0 0 = DvuD DvuD vuH 2）巴特沃斯高通滤波器 n阶巴特沃斯高通滤波器的传递函数定义如下 H(u，v)=1/1+( D0/D(u，v)2n 3）指数滤波器 指数高通滤波器的传递函数为 )74 . 4(),( ),( 0

36、 = n vuD D evuH 4）梯形滤波器 梯形高通滤波器的定义为 )84 . 4( Dv)D(u,1 D),(D Dv)D(u,0 v)H(u, 0 01DD D-v)D(u, 1 10 1 = vuD 四种滤波函数的选用类似于低通。理 想高通有明显振铃现象，即图像的边缘有 抖动现象；Butterworth高通滤波效果较 好，但计算复杂，其优点是有少量低频通 过，H(u，v)是渐变的，振铃现象不明显 指数高通效果比Butterworth差些，振铃现象不明显； 梯形高通会产生微振铃效果，但计算简单，较常用。 一般来说，不管在图像空间域还是频率域，采用高频 滤波不但会使有用的信息增强，同时也

37、使噪声增强。因此 不能随意地使用。 4.5 彩色增强技术 人眼的视觉特性 ： 分辨的灰度级介于十几到二十几级之间 ； 彩色分辨能力可达到灰度分辨能力的百倍以上。 彩色增强技术是利用人眼的视觉特性，将灰度图像变 成彩色图像或改变彩色图像已有彩色的分布，改善图像的 可分辨性。彩色增强方法可分为伪彩色增强和假彩色增强 两类。 4 4. .5 5. .1 1 伪彩色增强伪彩色增强 伪彩色增强是把灰度图像的各个不同灰度级按照线性 或非线性的映射函数变换成不同的彩色，得到一幅彩色图 像的技术。使原图像细节更易辨认，目标更容易识别。 伪彩色增强的方法主要有密度分割法、灰度级一彩色 变换和频率域伪彩色增强三种

38、。 1.1.密度分割法密度分割法 密度分割法是把灰度图像的灰度级从0（黑）到M0（白） 分成N个区间Ii(i=1，2，N),给每个区间Ii指定一种彩 色Ci，这样，便可以把一幅灰度图像变成一幅伪彩色图像。 该方法比较简单、直观。缺点是变换出的彩色数目有限。 2 2. .空间域灰度级一彩色变换空间域灰度级一彩色变换 根据色度学原理，将原图像f(x,y)的灰度范围分段， 经过红、绿、蓝三种不同变换TR()、TG()和TB()，变成 三基色分量IR(x,y)、IG(x,y)、IB(x,y)，然后用它们分别 去控制彩色显示器的红、绿、蓝电子枪，便可以在彩色显 示器的屏幕上合成一幅彩色图像。 3.3.频

39、率域伪彩色增强频率域伪彩色增强 频率域伪彩色增强方法的步骤为： 1）把灰度图像经傅立叶变换到频率域，在频率域内用三个不 同传递特性的滤波器分离成三个独立分量； 2）然后对独立分量进行逆傅立叶变换，便得到三幅代表不同 频率分量的单色图像，接着对这三幅图像作进一步的处理 （如直方图均衡化）； 3）最后将它们作为三基色分量分别加到彩色显示器的红、绿、 蓝显示通道，得到一幅彩色图像。 4.5.2 4.5.2 假彩色增强假彩色增强 假彩色增强是对一幅自然彩色图像或同一景物的多 光谱图像，通过映射函数变换成新的三基色分量，彩色 合成使感兴趣目标呈现出与原图像中不同的、奇异的彩 色。 假彩色增强目的：一是使

40、感兴趣的目标呈现奇异的 彩色或置于奇特的彩色环境中，从而更引人注目；一是 使景物呈现出与人眼色觉相匹配的颜色，以提高对目标 的分辨力。 多光谱图像的假彩色增强可表示为 将可见光与非可见光波段结合起来，通过假彩色处 理，就能获得更丰富的信息，便于对地物识别。 ) 15 . 4 ( ,.,., ,.,., ,.,., 21 21 21 = = = iBF iGF iRF gggfB gggfG gggfR 对于自然景色图像，通用的线性假彩色映射可表示为 例如采用以下的映射关系 则原图像中绿色物体会呈红色，蓝色物体会呈绿色，红 色物体则呈兰色。 伪彩色增强与假彩色增强有何区别？ ) 25 . 4 ( 333 222 111 = f f f F F F B G R cba cba cba B G R )35 . 4( 001 100 010 = f f f F F F B G R B G R ？