医学图像锐化系统设计汇编.docVIP

医学图像处理课程设计说明书 医学图像锐化系统设计 班 级 生医113班 姓 名 付道远 学 号 111417010308 指导教师 李振伟、何琳 目 录 1 设计目的 1 2 设计内容与原理 1 3 设计思路与方案 5 4 详细设计步骤 6 5 结果与分析 10 6 总结 14 参考文献 15 附录 17 基于灰度变换的医学图像增强系统设计 1 设计目的 1、巩固《医学图像处理》。 2、掌握MATLAB GUI设计基本原理方法，进行医学图像应用处理系统的设计与开发。]。这将导致原始图像在平滑处理之后，图像边缘和图像轮廓模糊的情况出现。 为了减少这类不利效果的影响，就需要利用图像锐化技术，使图像的边缘变得清晰。图像锐化处理的目的是为了使图像的边缘、轮廓线以及图像的细节变得清晰，经过平滑的图像变得模糊的根本原因是因为图像受到了平均或积分运算，因此可以对其进行逆运算(如微分运算)就可以使图像变得清晰。从频率域来考虑，图像模糊的实质是因为其高频分量被衰减，因此可以用高通滤波器来使图像清晰。但要注意能够进行锐化处理的图像必须有较高的性噪比，否则锐化后图像性噪比反而更低，从而使得噪声增加的比信号还要多，因此一般是先去除或减轻噪声后再进行锐化处理[]。 2）图像锐化的目的 加强图像的轮廓边缘或轮廓，使图像看起来比较清晰。与图像平滑不同，图像平滑通过积分过程使得图像边缘模糊，图像锐化则微分而使图像边缘突出、清晰。 由于我们常常无法事先确定轮廓的取向，因而挑选用于轮廓增强的微分算子时，必须选择那些不具备空间方向性的和具有旋转不变的线性微分算子。最基本的一类边缘检测算子是微分类算子。包括：梯度算子、Robel梯度算子、Sobel算子、Prewitt边缘检测算子、Kridch边缘检测算子、Laplacian算子等。 2.1 Sobel微分算子函数 Sobel算子具有一定的噪声抑制能力，在检测阶跃边缘时得到的边缘宽度至少为二像素，它不依赖于边缘方向的二阶微分算子，是一个标准量而不是向量，具有旋转不变即各向同性的性质。在图像处理中常被用于提取图像的边缘。 Sobel微分算子定义如下： （式1） （式2） 如果用图像处理的模板来表示，则有： 　　　　　 （式3） 下面以一个简单的例子来体会该算法。设原图像为 对中模板罩不住的地方其处理结果令为0，其它的像素一次按照式（式1）进行计算。例如，对的模板下的图像子块为： 计算结果为 最终处理结果为 2.2 Priwitt微分算子函数 Priwitt微分算子的思路与Sobel微分算子的思路类似，是在一个奇数大小的模板中定义其微分运算。 Priwitt微分算子定义如下： 2.3拉普拉斯微分算子函数 最简单的二阶各向同性微分算子是拉普拉斯微分算子，二维图像f（x,y）的拉普拉斯微分算子定义为： 写成模板系数形式形式即为Laplacian算子： 二阶微分算子所提取出的细节较一阶微分算子提出的细节多，表明了二阶微分算子在对图像细节更加敏感。 2.4理想高通滤波器的转移函数： 0， D(u,v) ≤ D0 H(u,v)= 1, D(u,v) D0 2.5 Butterworth滤波器的转移函数： H(u,v)=1/{1+[} 2.6 指数高通滤波器的转移函数： H(u,v)= 注：式中，D0为截止频率；D(u,v)=是点(u,v)到频率平面原点的距离。 3 设计思路与方案 利用matlab数字图像处理，编写相关函数，具体过程如下：双击打开MATLAB 7.0→File→New→GUI→单击，调整axes1大小→单击OK，调整按钮大小和颜色，修改名称→再建axes2→单击OK，调整按钮大小和颜色，修改名称→保存→View→M-file Edit→针对每个功能按钮，输入对应的函数。 图像锐化设计流程图： 4 详细设计步骤 4.1 Sobel锐化算子代码 w=fspecial(sobel); k1=imfilter(im,w); axes(handles.axes2); imshow(k1,[]); 4.2 Prewitt锐化算子代码 w=fspecial(sobel); k1=imfilter(im,w); axes(handles.axes2); imshow(k1,[]); 4.3 Log锐化算子代码 w=fspeci