医学影像图像处理2.ppt

1、第2章 医学影像图像数字化及其运算图像数字化 图像数字化是将一幅画面转化成计算机能处理的形式数字图像的过程。模拟图像 数字图像 正方形点阵 具体来说，就是把一幅图画分割成如图所示的一个个小区域（像元或像素），并将各小区域 灰度用整数来表示，形成一幅数字图像。它包括采样和量化两个过程。小区域的位置和灰度就是像素的属性。采样 将空间上连续的图像变换成离散点的操作称为采样。采样间隔和采样孔径的大小是两个很重要的参数。当对图像进行实际的抽样时，怎样选择各抽样点的间隔是个非常重要的问题。关于这一点，图像包含何种程度的细微的浓淡变化，取决于希望忠实反映图像的程度。经采样图像被分割成空间上离散的像素，但其灰

2、度是连续的，还不能用计算机进行处理。将像素灰度转换成离散的整数值的过程叫量化。一幅数字图像中不同灰度值的个数称为灰度级数，用G表示。黑白图像黑白图像 是指图像的每个像素只能是黑或者白，没有中间的过渡，故又称为值图像。2值图像的像素值为0、1。2灰度图像灰度图像 灰度图像是指每个像素的信息由一个量化的灰度级来描述的图像，没有彩色信息。灰度图像描述示例彩色图像彩色图像 彩色图像是指每个像素的信息由RGB三原色构成的图像，其中RBG是由不同的灰度级来描述的。量化参数与数字化图像间的关系数字化方式可分为均匀采样、量化和非均匀采样、量化。所谓“均匀”，指的是采样、量化为等间隔。图像数字化一般采用均匀采样

3、和均匀量化方式。非均匀采样是根据图象细节的丰富程度改变采样间距。细节丰富的地方，采样间距小，否则间距大。非均匀量化是对像素出现频度少的间隔大，而频度大的间隔小。采用非均匀采样与量化，会使问题复杂化，因此很少采用。一般来说，采样间隔越大，所得图像像素数越少，空间分辨率低，质量差，严重时出现像素呈块状的国际棋盘效应；采样间隔越小，所得图像像素数越多，空间分辨率高，图像质量好，但数据量大。对于数字图像，对于数字图像，且且 1234566432211664663456661466231364661234565456214灰度直方图一幅图像对应唯一的灰度直方图，反之不成立。不同的图像可对应相同的直方图。

4、下图给出了一个不同的图像具有相同直方图的例子。一幅图像分成多个区域，多个区域的直方图之和即为原图像的直方图。图像的直方图H(i)=区域的直方图H1(i)+区域的直方图H2(i)三、直方图的应用1.数字化参数（判断量化是否恰当）直方图给出了一个简单可见的指示，用来判断一幅图象是否合理的利用了全部被允许的灰度级范围。一般一幅图应该利用全部或几乎全部可能的灰度级，否则等于增加了量化间隔。丢失的信息将不能恢复。(a)恰当量化 (b)未能有效利用动态范围 (c)超过了动态范围.边界阈值选取（确定图像二值化的阈值）假设某图象的灰度直方图具有 二峰性，则表明这个图象的较亮的区域和较暗的区域可以较好地分离，以

5、这一点为阈值点，可以得到好的值处理的效果。3.当物体部分的灰度值比其它部分灰度值大时，可利用直方图统计图像中物体的面积。A=式中n为图像像素总数，vi是图像灰度级为i的像素出现的频率。4.计算图像信息量H（熵）Pi是图像灰度级为i的像素出现的频率，图像的灰度范围在0，L-1。彩色图的灰度直方图灰度图的灰度直方图灰度图具有二峰性具有二峰性的灰度图的2值化例：设图象有例：设图象有64*64*6464=4096=4096个象素，有个象素，有8 8个灰度级，灰度个灰度级，灰度分布如表所示分布如表所示。进行进行直直方图均衡化方图均衡化。rkr0=0r1=1/7r2=2/7r3=3/7r4=4/7r5=5

6、/7r6=6/7r7=1 nk 790102385065632924512281p(rk)0.190.250.210.160.080.060.030.02rkr0=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.00sk舍入舍入 1/73/75/76/76/7111nsk 7901023850985448p(sk)0.190.250.210.240.11均衡化

7、前后均衡化前后直直方图比较方图比较直直方图均衡化方图均衡化均衡化直方图均衡化 直方图均衡化实质上是减少图象的灰度级以换取对比度的加大。在均衡过程中，原来的直方图上频数较小的灰度级被归入很少几个或一个灰度级内，故得不到增强。若这些灰度级所构成的图象细节比较重要，则需采用局部区域直方图均衡。医学图像灰度直方图均衡化平移变换公式yx放大、缩小变换公式使图像在x轴方向放大(缩小)c倍，在y轴方向上放大(缩小)d倍 由于放大(缩小)算子运算不是一一映射，只是简单的重复放大，因此将产生所谓的“方块”效应，为改善这种可视效果，需要进行插值运算 yx旋转变换公式变换公式对图像绕原点顺时针角旋转变换 x x yyXY令f(x,y)为两个变量的函数，其在单位正方形顶点的值已知，假设我们希望通过插值得到正方形内任意点的f(x,y)的灰度值，则可由双曲线方程从a到d这四个系数需由已知的四个顶点的f(x,y)灰度值拟合 xy123xy123步骤1对上端的两个顶点进行线性插值，可得：步骤2对底端的两个顶点进行线性插值，可得 步骤3进行垂直方向的线性插值 合并