最新医学图像处理实验报告.docx

最新医学图像处理实验报告.docx

《最新医学图像处理实验报告.docx》由会员分享，可在线阅读，更多相关《最新医学图像处理实验报告.docx（11页珍藏版）》请在冰豆网上搜索。

最新医学图像处理实验报告

医学图像处理实验报告

班级专业姓名学号

实验名称：

图像增强

一、实验目的

1：

理解并掌握常用的图像的增强技术。

2：

熟悉并掌握MATLAB图像处理工具箱的使用。

3：

实践几种常用数字图像增强的方法，增强自主动手能力。

二、实验任务

对于每张图像（共三张图片），实现3种图像增强方法。

根据图像的特点，分别选用不用的图像增强算法。

三、实验内容（设计思路）

1、artery_vessel

（1）直方图均衡化

直方图是图像的最基本的统计特征，它反映的是图像的灰度值的分布情况。

直方图均衡化的目的是使图像在整个灰度值动态变化范围内的分布均匀化，改善图像的亮度分布状态，增强图像的视觉效果。

灰度直方图是图像预处理中涉及最广泛的基本概念之一。

图像的直方图事实上就是图像的亮度分布的概率密度函数，是一幅图像的所有象素集合的最基本的统计规律。

直方图反映了图像的明暗分布规律，可以通过图像变换进行直方图调整，获得较好的视觉效果。

直方图均衡化是指：

采用累积分布函数（CDF）变化生成一幅图像，该图像的灰度级较为均衡化，且覆盖了整个范围[0，1]，均衡化处理的结果是一幅扩展了动态范围的图像。

直方图均衡化就是通过灰度变换将一幅图像转换为另一幅具有均衡直方图，即在每个灰度级上都具有相同的象素点数的过程。

主要用途是：

将一幅灰度分布集中在较窄区间，细节不够清晰的图像，修正后使图像的灰度间距增大或灰度分布均匀，令图像的细节清晰，达到图像增强的目的。

（2）中值滤波加直方图均衡化

中值滤波法是一种非线性平滑技术，它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值。

中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围的像素值接近的真实值，从而消除孤立的噪声点。

方法是用某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。

二维中值滤波输出为g（x,y）=med{f（x-k,y-l）,（k,l∈W）}，其中，f（x,y），g（x,y）分别为原始图像和处理后图像。

W为二维模板，通常为3*3，5*5区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。

（3）直方图均衡化加维纳滤波

维纳滤波也称为最小均方误差滤波器，或最小二乘方误差滤波器。

它是建立在人为图像和噪声是随机过程的基础上，而目标是找一个未污染图像f的估计值，使它们之间的均方误差最小。

误差度量由下式给出：

其中，是参数的期望值。

这里假设噪声和图像不相关，其中一个或另一个有零均值，且估计中的灰度级是退化图像中灰度级的线性函数。

基于这些条件，误差函数的最小值在频率域中由如下表达式给出：

式中，为退化函数，为其复共轭，

为噪声的功率谱，为未退化图像的功率谱。

2、Brain_CT

（1）均衡化加巴特沃斯低通滤波

截止频率位于距原点处的n阶巴特沃斯低通滤波器的传递函数定义为：

其中，

巴特沃斯低通滤波器传递函数并没有在通过频率和滤除频率之间给出明显截止的尖锐的不连续性。

对于具有平滑传递函数的滤波器，可在这样一点上定义截止频率，即使得H（u，v）下降为其最大值的某个百分比的点，在上式中，截止频率点是当的点（即H（u，v）从其最大值1下降为50%）。

（2）空域线性平滑滤波

平滑线性空间滤波器的输出（响应）是包含在滤波器模版邻域内的像素的简单平均值，这些滤波器也称为均值滤波器。

它使用滤波器模版确定的邻域内像素的平均灰度值代替图像中每个像素的值，这种处理的结果降低了图像灰度的“尖锐”变化。

即对待处理的当前像素点（x，y），选择一个模板，该模板由其近邻的若干像素组成，求模板中所有像素的均值，再把该均值赋予当前像素点（x，y），作为处理后图像在该点上的灰度值g（x，y），即，m为该模板中包含当前像素在内的像素总个数。

（3）高斯滤波

高斯滤波器的二维形式为：

，D（u，v）是距频率矩形中心的距离。

是关于中心的扩展度的度量。

通过令，该滤波器可表示为：

其中，是截止频率。

当时，H下降到其最大值的0.607处。

3、BreastCancer_a

（1）同态滤波

同态滤波是把频率过滤和灰度变换结合起来的一种图像处理方法，它依靠图像的照度/反射率模型作为频域处理的基础，利用压缩亮度范围和增强对比度来改善图像的质量。

使用这种方法可以使图像处理符合人眼对于亮度响应的非线性特性，避免了直接对图像进行傅立叶变换处理的失真。

同态滤波的基本原理是：

将像元灰度值看作是照度和反射率两个组份的产物。

由于照度相对变化很小，可以看作是图像的低频成份，而反射率则是高频成份。

通过分别处理照度和反射率对像元灰度值的影响，达到揭示阴影区细节特征的目的。

同态滤波处理的基本流程如下：

f（x,y）---->Log---->DFT---->频域滤波---->IDFT---->Exp---->g（x,y）

其中f（x,y）表示原始图像；g（x,y）表示处理后的图像；Log代表对数运算；DFT代表傅立叶变换（实际操作中运用快速傅立叶变换FFT）；IDFT代表傅立叶逆变换（实际操作中运用快速傅立叶逆变换IFFT）；Exp代表指数运算。

（2）四八领域均值滤波

均值滤波是典型的线性滤波算法，它是指在图像上对目标像素给一个模板，该模板包括了其周围的临近像素（以目标象素为中心的周围8个像素，构成一个滤波模板，即去掉目标像素本身），再用模板中的全体像素的平均值来代替原来像素值。

（3）巴特沃斯高通滤波

截止频率为的n阶巴特沃斯高通滤波器定义为：

式中，。

与巴特沃斯低通滤波器一样，高通滤波器更平滑。

四、实验步骤

1、artery_vessel

（1）直方图均衡化

a、打开计算机，启动MATLAB程序；程序组中“work”文件夹中应有待处理的图像文件； 

b、调入“实验一”中获取的数字图像，并进行计算机均衡化处理； 

c、显示原图像和经过均衡化处理过的图像； 

d、记录和整理实验报告。

（2）中值滤波加直方图均衡化

a、将模板在图中漫游，并将模板中心与图中某个像素位置重合； 

b、读取模板下各对应像素的灰度值； 

c、将这些灰度值从小到大排成1列； 

d、找出这些值中排在中间的1个； 

e、将这个中间值赋给对应模板中心位置的像素；

f、中值滤波之后的像素值进行直方图均衡化处理；

g、显示原图像和处理后的图像，

h、记录和整理实验报告。

（3）直方图均衡化加维纳滤波

a、读取图像；

b、对数字图像进行直方图均衡化处理；

c、选择滤波器窗口为5*5，对均衡化后的图像进行二维自适应维纳滤波；

d、对滤波之后的图像进行灰度变换；

e、显示原图像和处理后的图像。

2、Brain_CT

（1）均衡化加巴特沃斯低通滤波

a、读取图像；

b、对数字图像进行直方图均衡化处理；

c、将图像矩阵转换成double型，对数字图像进行二维傅里叶变换并进行fftshift调整；

d、设计二阶巴特沃斯低通滤波器；

e、频域滤波后，进行傅里叶反变换转换到空域；

f、将数据转换成uint8型，并显示原图像和处理后的图像。

（2）空域线性平滑滤波

a、读入图像；

b、对数字图像进行直方图均衡化处理；

c、利用预定义函数fspecial 命令产生平均（average）滤波器 

d、分别对其进行10*10、5*5、2*2的均值滤波；

e、显示原图像和选用不同大小模版处理后的图像。

（3）高斯滤波

a、读入图像；

b、对数字图像进行直方图均衡化处理；

c、选择高斯滤波参数（标准差）sigma为1.6；

d、选择滤波器尺寸为5*5；

e、创建高斯滤波器进行滤波；

f、显示原图像和处理后的图像。

3、BreastCancer_a

（1）同态滤波

a、读入图像；

b、对数字图像进行直方图均衡化处理；

c、转换图像矩阵为双精度型；

d、取对数；

e、对其做傅里叶变换；

f、选择参数，截止频率为10，锐化系数为2，=1.5，=2.0；

g、进行高斯同态滤波；

h、滤波之后进行傅里叶逆变换；

i、显示原图像和处理后的图像。

（2）四八领域均值滤波

a、读入图像；

b、转换图像矩阵为双精度型；

c、创建4邻域平均滤波模版[010;101;010]；

创建8邻域平均滤波模版[111;101;111];

d、进行滤波；

e、显示原图像和处理后图像。

（3）巴特沃斯高通滤波

a、读取图像；

b、对数字图像进行直方图均衡化处理；

c、将图像矩阵转换成double型，对数字图像进行二维傅里叶变换并进行fftshift调整；

d、设计十阶巴特沃斯低通滤波器；

e、频域滤波后，进行傅里叶反变换转换到空域；

f、将数据转换成uint8型，并显示原图像和处理后的图像。

五、实验结果和分析

1、artery_vessel

（1）直方图均衡化

实验结果：

由图得：

直方图均衡化后，图像对比度增强，血管更加清晰。

（2）中值滤波加直方图均衡化

实验结果：

由图得：

3*3中值滤波效果不明显，5*5和7*7滤波之后对比度增加，但是上方亮区过亮，是的上方血管分辨不清，同时中值滤波之后，细节（即图中较细的血管）变得难以分辨。

（3）直方图均衡化加维纳滤波

实验结果：

由图得：

经过均衡化维纳滤波之后，图像的对比度增强，但上方亮区极度增强，几乎已经难以分辨该区域血管的存在，但细节观察由于中值滤波。

综上：

单纯的直方图均衡化能够对该血管图更好的增加对比度且对于细节的展示得以增强，同时对于上方亮区的血管的分辨影响不大。

表明该图中，不存在太多噪声干扰，仅均衡化增加其对比度即有明显效果。

2、Brain_CT

（1）均衡化加巴特沃斯低通滤波

实验结果：

由图得：

巴特沃斯低通滤波处理后与清晰的原图像相比，变得更平滑、更柔和。

（2）空域线性平滑滤波

手工艺品，它运用不同的材料，通过不同的方式，经过自己亲手动手制作。

看着自己亲自完成的作品时，感觉很不同哦。

不论是01年的丝带编织风铃，02年的管织幸运星，03年的十字绣，04年的星座手链，还是今年风靡一时的针织围巾等这些手工艺品都是陪伴女生长大的象征。

为此，这些多样化的作品制作对我们这一创业项目的今后的操作具有很大的启发作用。

实验结果：

（二）创业优势分析由图得：

空域线性平滑滤波会模糊一副图像，能使得较小物体的灰度与背景混合起来，较大物体变得像“斑点”易于检测。

同时，模版尺寸m越大，模糊效果越明显。

1．www。

cer。

net/artide/2004021313098897。

shtml。

在大学生对DIY手工艺品价位调查中，发现有46%的女生认为在十元以下的价位是可以接受；48%的认为在10-15元；6%的则认为50-100元能接受。

如图1-2所示

我们从小学、中学到大学，学的知识总是限制在一定范围内，缺乏在商业统计、会计，理财税收等方面的知识；也无法把自己的创意准确而清晰地表达出来，缺少个性化的信息传递。

对目标市场和竞争对手情况缺乏了解，分析时采用的数据经不起推敲，没有说服力等。

这些都反映出我们大学生创业知识的缺乏；

（3）高斯滤波

实验结果：

由图得：

高斯低通滤波后，图像同样变得模糊，较大物体得以更好表现。

综上：

低通滤波后，能使得图像模糊，看上去更平滑、更柔和，较小物体得以与背景融合，较大物体更易于检查。

3、BreastCancer_a

（1）同态滤波

调研要解决的问题：

实验结果：

由图得：

直方图均衡化后，对比度明显增加；同态滤波后，图像中暗区更暗，使得肿瘤更为清晰，可以见到更多的细节。

（2）四八领域均值滤