MRI图像增强和灰度插值.doc

课程设计报告 题 目： 专 业： 班 级： 学 号： 姓 名： 指导教师： 天津医科大学 生物医学工程年 月 日 在学习了数字信号处理和医学图像处理等课程的基础上，以MRI学习医学图像的读取保存等基本操作，图像的灰度映射增强以及 三．课题研究的主要内容 1、原始图像是MRI脑图像，分别为t1加权像，t2加权像和pd质子加权像。根据图像的不同类型，把原始MRI图像转为灰度图像，将图像拆分为大小相等的t1, t2, pd三幅图像 3、绘出t1, t2, pd三幅图像的灰度直方图 4、对t1, t2, pd三幅图像t1, t2, pd三幅图像 3、利用imhist函数绘出三幅图像的灰度直方图，再分别用histeq函数，自定义的 twomodegauss函数，imadjust函数对三幅图像进行直方图均衡化，直方图规定化， Gamma 校正增强处理； 4、对t1, t2, pd三幅图像 五．步骤与结果 1、利用imfinfo函数获取所给MRI图像信息，程序见附录（1），结果如下: Width: 520 Height: 284 ColorType: indexed 由此可知原始MRI图像的宽为520列，高为284行，图像类型是索引图像； 2、利用ind2gray函数将原始图像转换为灰度图像，程序见附录（1），结果如图1: 图1 原始MRI的灰度图像 3、根据得到的图像的矩阵信息将图像拆分为大小相等的t1, t2, pd三幅图像 图2.1 原始t1加权像的灰度图像 图2.2 原始t2加权像的灰度图像 图2.3 原始pd质子加权像的灰度图像 绘出t1, t2, pd三幅图像的灰度直方图 图3.1 t1加权像的灰度直方图 图3.2 t2加权像的灰度直方图 图3.3 pd质子加权像的灰度直方图 结果分析：由三幅加权像的灰度直方图可以看出，脑组织的灰度值主要分布在40-150 之间。 图像增强处理 直方图均衡化 利用histeq函数对三幅加权像的灰度图像进行均衡化处理，程序见附录（3），结果如图4.1~4.6： 图4.1 t1加权像均衡化后图像 图4.2 t1加权像均衡化后图像的灰度直方图 图4.3 t2加权像均衡化后图像 图4.4 t2加权像均衡化后图像的灰度直方图 图4.5 pd加权像均衡化后图像 图4.6 pd加权像均衡化后图像的灰度直方图 结果分析： 直方图均衡化的基本思想是把原始图像的直方图变换成均匀分布的形式，这样就增加了图像灰度值的动态范围，从而达到了增强图像整体对比度的效果。但这种方法的缺点是它对处理的数据不加选择，它可能会增加背景杂讯的对比度，并且降低有用信号的对比度，变换后图像的灰度级减少，某些细节消失。从而导致三幅加权像均衡化效果不佳。 直方图规定化 利用自定义的一个双峰高斯函数作为规定化函数对三幅加权像的灰度图像进行规定化处理，双峰高斯函数如图5.1，程序见附录（4），规定化程序见附录（5），结果如图5.2~5.7：