评分
大理大学实验报告
课程名称 医学图像处理 实验名称 图像分割
专业班级 13级生物医学工程1班 学 号 2013166143 2013166141 实验日期 2015/12/17 实验地点 三教204
2015—2016学年度第 2 学期
姓 名 范小欢 杨飞
《医学图像处理》
一、 实验目的 掌握常用的边缘提取算法，哦那个图像中提取感兴趣的区域，实现图像分割。在图像中，寻找灰度相同或相似的区域，区分图像中的背景区域和目标区域，利用matlab实现图像的边缘检测，进行图像分割。 二、实验环境 1、硬件配置：处理器：Intel（R） celeron(R) CPU N2940 @1.83GHz 1.83GHz 安装内存（RAM）：4.00GB 系统类型：64位操作系统 2、软件环境：MATLAB R2012b 三、实验内容 (包括本实验要完成的实验问题及需要的相关知识简单概述) 一、实验原理 1、 函数imadjust是对灰度图像进行灰度变换的基本工具。其语法为 g=imadjust（f，[low_in high_in],[low_out high_out],gamma） 2、直方图均衡化 直方图生成函数为imhist，其基本语法为 g=imhist（f，b） 直方图均衡化用工具箱中的函数histeq histeq（f，nlev） 3、直方图规定化 函数histeq实现直方图匹配的形式为 g=histeq（f，hspec） 4、图像正交变换 傅里叶频谱可以使用函数abs来获得 F=abs（Y） 傅里叶逆变换用函数ifft2实现，基本语法为 F=ifft2（F） 第2页 共16页
《医学图像处理》
四、实验结果与分析 （包括实验原理、数据的准备、运行过程分析、源程序（代码）、图形图象界面等） 【例题】 例1手动阈值分割 ① 源代码 [I,map]=imread('eight.tif'); %读取原图像eight imshow(I); figure;J=imhist(I);imhist(I);[M,N]=size(I); %显示图像并生成直方图，图像大小为M*N像素 for i=1:1:M %用for语句依次读取图像像素点，若超过了指定的阈值200，取像素为0，否则为1 for j=1:1:N if I(i,j)>200 %若像素点大于200，则取0，反之 g(i,j)=0; else g(i,j)=1; %否则取0 end %if循环的end end end figure;imshow(g); %显示输出图像g ② 图像输出 图1原图像 图3 迭代阈值二值图像 图2 原图像的灰度直方图 ③ 结果分析 从图1可看出原图像灰度级单一，可选择一个合适的阈值来进行阈值分割，由图2可知，其灰度直方图中出现双峰，即利用双峰法选取阈值为200（双峰之间波谷的灰度值即作为阈值），将目标和背景分离，产生相应的二值图像（如图3）。 第3页 共16页