第六次课_图像增强_直方图矫正_最后

1、图像增强图像增强（Image EnhancementImage Enhancement） 概述概述一、处理原因：一、处理原因： 图像形成、传输、转换、显示产生降质图像形成、传输、转换、显示产生降质二、改善方法：二、改善方法：图像增强：目的是为了改善图像的视觉效果，或图像增强：目的是为了改善图像的视觉效果，或者是为了更便于人或机器的分析和处理，提高者是为了更便于人或机器的分析和处理，提高图像的可懂度。在不考虑降质原因的情况下，图像的可懂度。在不考虑降质原因的情况下，用试探的方式对图像进行加工，力求改善图像用试探的方式对图像进行加工，力求改善图像的质量，如突出了一部分信息，同时可能压制的质量，如突

2、出了一部分信息，同时可能压制另一部分的信息。另一部分的信息。 图像复原：图像复原：考虑降质原因，分析降质模型，考虑降质原因，分析降质模型，试图利用退化现象的某种先验知识，把已试图利用退化现象的某种先验知识，把已经退化了的图像加以重建或恢复。经退化了的图像加以重建或恢复。 三、三、增强的方法：增强的方法：1.1.空域法空域法 2.2.频域法频域法 四、四、图像增强技术：图像增强技术：主要可分为主要可分为1 1点运算增强算子：如图像灰度倒置、对比点运算增强算子：如图像灰度倒置、对比度伸缩、灰度动态范围的伸缩、灰度级分片、度伸缩、灰度动态范围的伸缩、灰度级分片、图像减影、直方图修正等；图像减影、直方

3、图修正等；2 2区域（模板）运算增强算子：如平滑、中区域（模板）运算增强算子：如平滑、中值滤波、值滤波、 锐化等；锐化等；3 3变换增强算子：如低通滤波、高通滤波、变换增强算子：如低通滤波、高通滤波、带通滤波、同态增晰等；带通滤波、同态增晰等；4 4色彩算子：如伪彩色处理。色彩算子：如伪彩色处理。灰度修正（空域法）灰度修正（空域法）如如K=1K=1，即为点增强处理，即为点增强处理特点：特点：1 1） 输出图像在像素点输出图像在像素点(m, n)(m, n)的灰度的灰度值值 g(m, n)g(m, n)仅取决于输入图像在像素点仅取决于输入图像在像素点(m, (m, n)n)的灰度值的灰度值f(m

4、, n) f(m, n) ，与像素点，与像素点(m, n)(m, n)的的邻近点无关。邻近点无关。2 2） 我们通常写成我们通常写成 s = T(r) ,s = T(r) ,其中其中s s是输出像素点值，是输出像素点值，r r是输出像素点值。是输出像素点值。3 3） T T可以是任一从可以是任一从0,10,1到到0,10,1映射映射的递增函数。的递增函数。灰度变换灰度变换( (一一) )线性灰度变换线性灰度变换 当图象成象时曝光不足或过度当图象成象时曝光不足或过度, , 或由于或由于成象设备的非线性和图象记录设备动态范围成象设备的非线性和图象记录设备动态范围太窄等因素太窄等因素, ,都会产生对

5、比度不足的弊病，都会产生对比度不足的弊病，使图象中的细节分辨不清。这时可将灰度范使图象中的细节分辨不清。这时可将灰度范围线性扩展。围线性扩展。 ayxfcbyxfacayxfabcdbyxfdyxg),(),(),(),(),( 设设f(x,y)f(x,y)灰度范围为灰度范围为a,ba,b，g(x,y)g(x,y)灰度范灰度范围为围为c,dc,d。0f(x,y)g(x,y)abcd( (二二) )分段线性灰度变换分段线性灰度变换 将感兴趣的灰度范围线性扩展，相对抑将感兴趣的灰度范围线性扩展，相对抑制不感兴趣的灰度区域。制不感兴趣的灰度区域。 设设f(x,y)f(x,y)灰度范围为灰度范围为0,

6、Mf0,Mf，g(x,y)g(x,y)灰度范围为灰度范围为0,Mg0,MgayxfyxfacbyxfacayxfabcdMyxfbdbyxfbMdMyxgffg),(0),(),(),(),(),(),(0f(x,y)g(x,y)abcdMfMg( (三三) )非线性灰度变换非线性灰度变换 （1 1）对数变换）对数变换 低灰度区扩展，高灰度区压缩。低灰度区扩展，高灰度区压缩。 （2 2）指数变换）指数变换 高灰度区扩展，低灰度区压缩。高灰度区扩展，低灰度区压缩。cbyxfayxgln 1),(ln),(对数变换对数变换a,b,ca,b,c是按需要可以调整的参数。是按需要可以调整的参数。对数变换

7、对数变换低灰度区扩展，高灰度区压缩低灰度区扩展，高灰度区压缩1),(),(ayxfcbyxg指数变换指数变换a,b,ca,b,c是按需要可以调整的参数。是按需要可以调整的参数。指数变换指数变换高灰度区扩展，低灰度区压缩。高灰度区扩展，低灰度区压缩。（四）一些基本的灰度变换（四）一些基本的灰度变换1 1） 灰度倒置灰度倒置负像（负像（Image NegativesImage Negatives） s = L - 1 s = L - 1 r r2 2）对数变换对数变换 （Log TransformationsLog Transformations） ，c c是常数，且是常数，且 )1log(rcs

8、0r3 3）幂规律变换（）幂规律变换（Power-Law TransformationsPower-Law Transformations） 或或 ，c c和和是正数，是正数，有时也称有时也称GammaGamma校正。校正。crs )( rcsa)a)原图原图 b) =3.0 c)=4.0 d)=5.0b) =3.0 c)=4.0 d)=5.04 4）分段线性变换（）分段线性变换（Piecewise-Linear Piecewise-Linear Transformation FunctionsTransformation Functions） 对比度伸缩（对比度伸缩（Contrast str

9、etchingContrast stretching） 灰度级分片（灰度级分片（Gray-level slicingGray-level slicing） 使指定的灰度值范围高亮度使指定的灰度值范围高亮度左图转换函数相当于二值化（没有背景）左图转换函数相当于二值化（没有背景） 位平面分片（位平面分片（Bit-plane slicingBit-plane slicing）5 5）图像减影）图像减影 在这种情况，通过计算两帧相似图像在这种情况，通过计算两帧相似图像之间的不同，来增强（高亮度化）两帧间之间的不同，来增强（高亮度化）两帧间的区别的区别 这种方法已经应用于图像分割和增强，这种方法已经应用

10、于图像分割和增强，如如X X光线成像，医学上的运用光线成像，医学上的运用n 逻辑操作(二进制掩膜，binary masking)基于点运算，对两副图像的单个象素进行操作（此时每个象素的值都被看成逻辑值），基本包括与、或、非三者，其他任何逻辑操作都可通过三者之间的组合来完成。逻辑操作通常用于选择ROI (region of interest)，也常与形态学处理相结合。直方图修正直方图修正n一、定义一、定义q灰度直方图（灰度直方图（histogram）是）是灰度级的函数灰度级的函数，描述，描述的是图像中每种灰度级像素的个数，反映图像中每的是图像中每种灰度级像素的个数，反映图像中每种灰度出现的频率。

11、横坐标是灰度级，纵坐标是灰种灰度出现的频率。横坐标是灰度级，纵坐标是灰度级出现的频率。度级出现的频率。q图像及其灰度直方图的例图像及其灰度直方图的例（512像素像素*512像素）像素）n一、定义一、定义q对于连续图像，平滑地从中心的高灰度级变化到边对于连续图像，平滑地从中心的高灰度级变化到边缘的低灰度级。其直方图可定义为：缘的低灰度级。其直方图可定义为：q其中其中A(D)为为阈值面积函数阈值面积函数：为一幅连续图像中被具有灰为一幅连续图像中被具有灰度级度级D的所有轮廓线所包围的面积。的所有轮廓线所包围的面积。q对于离散函数，固定对于离散函数，固定D为为1，则，则H（D）=A（D）-A（D+1）

12、q一幅连续图像的直方图定义的示意图。一幅连续图像的直方图定义的示意图。)D(AdDdD)DD(A)D(Alim)DD(D)DD(A)D(Alim)D(H0D0D A2A1D=D2D=D1D2=D1+DXYD(X,Y)二、直方图的计算和性质二、直方图的计算和性质n1 1）计算）计算q依依据定义，若图像具有据定义，若图像具有L（通常（通常L=256,即即8位灰度级）位灰度级）级灰度，则大小为级灰度，则大小为MN的灰度图像的灰度图像f(x,y)的灰度直方图的灰度直方图hist0L-1可用如下计算获得：可用如下计算获得：q初始化初始化histk=0; k=0,L-1q统计统计histf(x,y)+;

13、x=0,M-1, y =0,N-1q归一化归一化histf(x,y)/M*N n2）直方图的性质）直方图的性质q不表示图像的空间信息；不表示图像的空间信息；q任一特定图像都有唯一直方图，但反之并不成立；任一特定图像都有唯一直方图，但反之并不成立；二、直方图的计算和性质二、直方图的计算和性质具有同样直方图的图像对应的直方图图像图像 直方图直方图 : : 确定确定 直方图直方图 图像图像: : 不确定不确定q归一化灰度直方图和面积函数可得到图像的归一化灰度直方图和面积函数可得到图像的概率密度函概率密度函数数PDF和和累积分布函数累积分布函数CDF。DDDDDdppHDCDFADHDPDFNSNLD

14、HAdppHDDAdppHApAdppHDDDAdDdDH)()()0(/)()()()0()(,0)()(0)()()()(25500对于离散图象，图象的面积则若令所以因为进行积分到，并等式两端从替换）（因为二、直方图的计算和性质二、直方图的计算和性质q若一幅图像包含一个灰度均匀一致，且背景与物体对比度很若一幅图像包含一个灰度均匀一致，且背景与物体对比度很强，假设物体的边界由灰度级强，假设物体的边界由灰度级D1定义的轮廓线，则定义的轮廓线，则nn500像素像素*546像素像素=n273000，最左侧直方，最左侧直方图峰值的灰度级为图峰值的灰度级为33 1DdD)D(H物物体体的的面面积积二、

15、直方图的计算和性质二、直方图的计算和性质n从灰度从灰度54到到255级级n约占图像总面积的约占图像总面积的60%2 25 55 55 54 4H H( (D D) )d dD D = =1 16 63 30 00 01 1二、直方图的计算和性质二、直方图的计算和性质q直方图的可相加性直方图的可相加性例如一副图像由若干个不相交的区域构成，则整幅图像例如一副图像由若干个不相交的区域构成，则整幅图像的直方图是这若干个区域直方图之和。的直方图是这若干个区域直方图之和。二、直方图的计算和性质二、直方图的计算和性质直方图修正是灰度级变换的常用方法直方图修正是灰度级变换的常用方法n直方图归一化直方图归一化原

16、始图像灰度级原始图像灰度级 归一化归一化 原始图像灰度原始图像灰度 归一化归一化 变换后图像灰度变换后图像灰度 归一化归一化 1,kzz1 , 0izr1 , 0rizs1 , 0sn直方图不反映灰度值的像素在图像中位置方面的任直方图不反映灰度值的像素在图像中位置方面的任何信息。何信息。 直方图反映的总体性质：明暗程度、细节是否清晰、动态范围大小等直方图修正的技术基础直方图修正的技术基础设设 应满足下列条件：应满足下列条件： 1 1）在区间内）在区间内T T( (r r) )为单值单调为单值单调增加；增加；2 2）对于，有。）对于，有。条件条件1 1）使灰度级保持从黑到白的次序；）使灰度级保持

17、从黑到白的次序；条件条件2 2）保证映射变换后的像素灰度值在）保证映射变换后的像素灰度值在允许范围内允许范围内 从从s s到到r r的反变换用下式表示的反变换用下式表示同样假设对于变量同样假设对于变量s s也要满足条件也要满足条件1 1）和）和2 2）。）。在一幅图像中，在在一幅图像中，在00，11区间内的区间内的灰度级是随机变量，假定对每一个灰度级是随机变量，假定对每一个瞬间它们是连续变量，那么可以用瞬间它们是连续变量，那么可以用概率密度函数和分别表示概率密度函数和分别表示原始图像和变换图像的灰度级。原始图像和变换图像的灰度级。 )(rpr)(sps由概率论知道，如果，T(r)是已知的，满足

18、条件1），那么变换图像灰度级的概率密度函数可以由下式得到：这说明：通过T(r)控制图像灰度级的概率函数，从而改善图像的外貌。关键是 如何选择？)(rpr 三、直方图均衡直方图均衡：把原始图的直方图变换为均匀分布的形式，增加像素灰度值的动态范围，提高图像对比度。T(r)Pr(r)Ps(s)S1S1rr在0,1区间内的任一个r值，都可以产生一个s值，且 ， 为变换函数。为使这种灰度变换具有实际意义， 应满足下列条件：在 区间， 为单调递增函数；在 区间，有 。用r和s分别表示归一化了的原始图像灰度和变换后的图像灰度。即 10 r10 s（0代表黑，1代表白）)(rTs )(rT)(rT10 r)(

19、rT10 r1)(0rT这里，条件保证灰度级从黑到白的次序，条件保证变换后的像素灰度仍在原来的动态范围内。 由s到r的反变换为:)(1sTr10 s)(1sT对s也满足条件和。对于连续图像，当直方图均衡化（并归一化）后有 ，即：)()(rdTdrrpdsr两边取积分得：rrdrrPrTs0)()(就是所求的变换函数，它表明变换函数是原图像的累计分布函数，是一个非负的递增函数。 1)(sPs) 1, 1 , 0, 10(/)(LkrNnrPkkkr对于离散图像，假定数字图像中的总像素为N，灰度级总数为个L，第k个灰度级的值为 ，图像中具有灰度级 的像素数目为 ，则该图像中灰度级 的像素出现的概率

20、（或称频数）为：krkrknkr对其进行均匀化处理的变换函数为：kjjkjjrkkNnrPrTS00)()(逆变换：10)(1kkkSsTr例：设图象有例：设图象有6464* *64=409664=4096个象素，有个象素，有8 8个灰度级，个灰度级，灰度分布如表所示。进行灰度分布如表所示。进行直直方图均衡化方图均衡化。rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1 nk 790102385065632924512281p(rk) 0.190.250.210.160.080.060.030.02步骤：步骤：rkr0=0r1=1/7r2=2/7r3=

21、3/7r4=4/7r5=5/7r6=6/7r7=1 nk 790102385065632924512281p(rk) 0.190.250.210.160.080.060.030.021. 1. 由（由（2-22-2）式计算）式计算s sk k。rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1 nk 790102385065632924512281p(rk) 0.190.250.210.160.080.060.030.02sk计算计算 0.190.440.650.810.890.950.981.00rkr0=0r1=1/7r2=2/7r3=3/7r4=

22、4/7r5=5/7r6=6/7r7=1 nk 790102385065632924512281p(rk) 0.190.250.210.160.080.060.030.02sk计算计算 0.190.440.650.810.890.950.981.00sk舍入舍入 1/73/75/76/76/71112. 2. 把计算的把计算的s sk k就近安排就近安排到到8 8个灰度级中。个灰度级中。rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5/7r6=6/7r7=1 nk 790102385065632924512281p(rk) 0.190.250.210.160.080.060.

23、030.02sk计算计算 0.190.440.650.810.890.950.981.00sk舍入舍入 1/73/75/76/76/7111sk nsk 7901023850985448p(sk) 0.190.250.210.240.113. 3. 重新命名重新命名sksk，归并相，归并相同灰度级的象素数。同灰度级的象素数。均衡化前后均衡化前后直直方图比较方图比较DAnkPr(rk)f取成整数倍均衡后直方图07900.190.191/7(0.14)0.191/710230.250.443/7(0.428)0.252/78500.210.655/7(0.714)0.213/76560.160.8

24、16/7(0.857)0.16+0.08=0.244/73290.080.896/7(0.857)5/72450.060.957/7(1.00)0.06+0.03+0.02=0.116/71220.030.987/7(1.00)7/7810.021.007/7(1.00)8个灰级，总计64*64=4096点例题：原图1/73/75/77/70.20.40.60.81.0DA变换图均衡后1/73/75/77/70.190.250.210.240.110.250.050.250.150.051/73/75/77/7仅存5个灰级，宏观拉平，微观不可能平，层次减少，对比度提高。直方图均衡化处理应注意的问题直方图均衡化处理应注意的问题直方图均衡化实质上是减少图象的灰度级直方图均衡化实质上是减少图象的灰度级以换取对比度的加大。在均衡过程中，原来的以换取对比度的加大。在均衡过程中，原来的直方图上频数较小的灰度级被归入很少几个或直方图上频数较小的灰度级被归入很少几个或一个灰度级内，故得不到增强。若这些灰度级一个灰度级内，故得不到增强。若这些灰度级所构成的图象细节比较重要，则需采用局部区所构成的图象细节比较重要，则需采用局部区域直方图均衡。域直方图均衡。结论：直方图均衡是一种非线性变换，以牺牲图像等级为代