《基于MATLAB的MRI及CT图像增强处理》由会员分享，可在线阅读，更多相关《基于MATLAB的MRI及CT图像增强处理（15页珍藏版）》请在金锄头文库上搜索。
1、共17 页第1 页课 程 设 计 说 明 书题目： 基于 MATLAB 的 MRI及 CT图像增强处理学院（系）：年级专业：学号：学生姓名：指导教师：教师职称：共17 页第2 页目录一、课程设计目的4二、课程设计要求三 、课程设计内容,四 、 题目分析五、总体设计5 六、程序设计 ,1、灰度转换5 2、灰度变换6 3、空余滤波 , ,7 4、频域滤波8 5、伪彩色增强11七、心得体会12 八、附录,13 附： 燕山大学课程设计评审意见表共17 页第3 页一、 课程设计目的学习使用 matlab 软件，并使用 matlab 软件编写相应的程序完成对图像的增强处理。二、 课程设计的要求1 熟悉和掌
2、握MATLAB 程序设计方法2 掌握 matlab 软件在图像增强技术的处理方法3 编写相应的matlab 程序进行图像处理三、 课程设计的内容学习MATLAB程序设计，利用MATLAB图像处理工具箱，设计和实现自己的matlab 图像增强设计要求。要求：按照软件工程方法，根据需求进行程序编程，给出设计详细说明。然后按照自己拟定的功能要求进行程序设计和调试。以下几点是程序必须实现的功能。1）图像的读取和保存。2）能够采用各种方法对图像进行处理，显示和对比变换前后的图像。3）编写程序并进行相应的程序说明。比较各种方法的好坏。4）图像直方图统计和直方图均衡，要求显示直方图统计，比较直方图均衡后的效
3、果。5）能对图像加入各种噪声，并通过几种滤波算法实现去噪并显示结果。比较去噪 效果。四、 题目分析信息化社会中， 计算机在各种信息处理中发挥着重要的作用。我们可以借助计算机，对数字图像进行处理，以达到不同的效果。根据题目的要求，除了实现要求的功能外，还有很多的功能需要用到。（1） 、由于医学图像采集到的都是灰色图像，而有许多图像格式并不符合要求，因此要通过编程将一个RGB 图像转换为灰度图像。（2）图像增强可以从以下几方面入手：灰色增强、空运增强、频域增强、彩色增强。（3）必须对图像读取和显示这种基本处理以及相应的对比。（4） 、选择 MA TLAB 计算编程可视化,功能强大 ,学、易用 ;
4、工具箱中包括图像处理函数 ,涵盖了几乎所有的图像处理技术方法; 具有开放性 ,用户可以根据需要书写自己的函数 ,以满足特定的需要;能够准确读取DICOM 格式文件 ,并有丰富的处理DICOM 格式图片的各种函数。医学影像设备产生的文件为DICOM 格式 ,普通图像处理软件无法读取共17 页第4 页此格式 ,若将其转换为 JPG或BMP 又会造成某些信息的丢失,MATLAB 可以克服这一缺点。（5） 、分析一个图像的频谱特征，利用傅里叶变换，将图像从空间域变换到频域，然后进行各种处理，经过高通滤波器或是低通滤波器。（6） 、由于医学图像采集到的都是灰色图像，可以通过matlab 进行伪彩色增强得
5、到相应的彩色图像。 、五、 总体设计本课程设计主要从以下四个方面对图像进行处理：灰度增强、空域增强、频域增强、彩色增强。1、 灰度增强：灰度变换的方法包括直方图灰度变换、直方图均衡化、直方图规定化。本课程设计主要采用直方图均衡化完成图像增强处理。2、 空域增强： 本课程设计采用平滑滤波器来完成空域滤波，其中采用均平滤波器filter2 、中值滤波器 medfilt2 、维纳滤波器wiener 。3、 频域滤波：在医学图像里采用低通滤波器进行滤波，采用高通滤波器进行边缘检出。4、 彩色增强：由于MRI 图像得到的都是灰色图像，而人眼对彩色敏感度更高，因此可以采用伪彩色增强的方法对图像进行增强处理
6、以得到彩色图像六、 程序设计1、灰度转换由于RGB 图像是三维图像，所以图像数据是一个三维数组，为了显示灰度图像，把三维图像降为二维，可以只取其中的二维数据。我在设计中采用rgb2gray 函数来完成数据的转换，但是其必须指出其开始是否灰色图像，否则容易出错。具体程序如下：I=imread(naogan.jpg); if isrgb(I) I2=rgb2gray(I); else msgbox(该图像是灰色图像，无需转换); end figure(1),imshow(I),title(原图 ) figure(2), imshow(I2),title(处理后 ) 处理后的图像如下图一所示：共17
7、 页第5 页2、灰度变换利用图像灰度级的分布可以看出图像灰度分布的特性。如果大部分像素集中在低灰度区则图像呈现暗特性,反之则呈现亮特性。 灰度变换的目的是通过改善直方图的灰度分布特性 ,进而改善图像的质量。灰度变换的方法包括直方图灰度变换、直方图均衡化、直方图规定化。本研究以直方图均衡化为例说明该模块的设计功能。直方图均衡化的基本思想是把原始图像的直方图变换成均匀分布的形式，这样就增加了图像灰度值的动态范围，从而达到了增强图像整体对比度的效果。均匀量化的自然图像的灰度直方图通常在低值灰度区间上频率较大，使得图像中较暗区域中的细节常常看不清楚。为了使图像清晰，可将图像的灰度范围拉开，让灰度频率较
8、小的灰度级变大，即让灰度直方图在较大的动态范国内趋于一致。理想情况下：使一输入图象转换为在每一灰度级上都有相同的像素点数（即输出的直方图是平的）。具体方法是：（1）列出原始图像的灰度级Sk，k 0,1，, ， L1；其中 L 是灰度级的个数。（2）统计原始图像各灰度级的像素数目nk。（3）计算原始图像直方图各灰度级的频率数。（4）计算原始图像的累计直方图。（5）取整计算， tk=int （N1）tkk N（6）确定映像关系，Sknk （7）统计新直方图各个灰度级的像素数目nk。（8）计算新的直方图，Pt（tk） nk / N。如要处理的 mri图像 ,通过直方图分析,原图像大面积为暗色,并且层
9、次不清,经过直方图均衡化后,直方图的灰度间隔被拉大,显得较为“平坦”,灰度层次等级增加,然后用此均衡得到理想的图像。其程序代码如下：I = imread (mri.tif) ; J = histeq( I) ; subplot (121) ;imshow( I) ; title ( 原图 ) subplot (122) ;imshow(J ) ; title ( 直方图均衡化后) ; figure ,subplot (121) ; imhist ( I ,64) ; title ( 原始图像直方图) ; subplot (122) ;imhist (J ,64) title ( 均衡变换后直方
10、图) 其图像处理结果如下：共17 页第6 页均衡化处理前后脑部mri图片均衡化前后的脑膜炎图片直方图比较3、空域增强空域滤波在处理效果上来分,可以分为平滑滤波器和锐化滤波器。下面以平滑滤波器为例介绍课件功能及实现。平滑滤波器的目的在于消除混在图像中的干扰,常使用的滤波器有均平滤波器filter2 ( ) 、中值滤波器 medfilt2 ( ) 、维纳滤波器 wiener () 。由于医学成像设备CT机自身的原因 ,在CT图像中容易产生高斯白噪声,严重影像了图像的质量 ,为了去除此类干扰可采用维纳滤波器;如果影像本身带有盐椒噪声,可采用中值滤波器。 通过给图像增加盐椒噪声( SAL TPEPPE
11、T) 或高斯噪 ( GUASSIAN) 来模拟医院的 ct及存在噪声的 mri 图像 ,然后采用不同滤波器进行处理,对比其效果。其对应的的 matlab程序如下：共17 页第7 页I=imread(mri.tif); J=imnoise ( I ,salt %添加盐椒噪声J=imnoise(I,gaussian,0,0.005); %添加高斯噪声K=filter2 (fspecial (average,3),I)/255; %均平过滤器L=medfilt2(I,3 3 ); %中值过滤器M=wiener2(I,5 5 ); %维纳滤波器subplot(1,5,1),imshow(I),titl
12、e( 原图 ) % 显示图像subplot(1,5,2),imshow(J),title( 添加噪声后的图片) subplot(1,5,3),imshow(K),title(均平过滤后的图片) subplot(1,5,4),imshow(L),title(中值过滤的图片) subplot(1,5,5),imshow(M),title(维纳滤波后的图片) 其相应的对比图如下：经过对比得知经中值滤波得到的图像效果最佳。4、频域滤波由于图像的边缘、跳跃部分及噪声颗粒是图像信号的高频分量,大面积的背景区域代表了图像信号的低频分量。因此 ,在频域内采用低通滤波器可以进行滤波,使用高通滤波器可以进行边缘检
13、出。其中,常用低通滤波器有理想低通滤波器、Butterworth 低通滤波器、指数低通滤波器和梯形低通滤波器等;常用高通滤波器有理想高通滤波器、Butterworth 高通滤波器、指数高通滤波器和梯形低通滤波器等。本次课程设计滤波器我采用二阶巴特沃斯滤波器进行滤波。为了得到图像的频谱图，先要对数据进行傅里叶变换，用fft2函数对二维数据进行快速傅里叶变换，同时为了更好的观察频谱图，需要把fft2 变换后的数据进行平移，利用fftshift 函数，把快速傅里叶变换的DC组件移到光谱中心。这样图像能量的低频成分将集中到频谱中心，图像上的边缘、线条细节信息等高频成分将分散在图像频谱的边缘。（1）高通
14、滤波器通过高通滤波器进行边缘检出，其程序代码如下：共17 页第8 页I= imread (mri.tif); %读入图像 : f=double(I); %数据类型转换k=fft2(f); %傅里叶转换g=fftshift(k); %转换数据矩阵M,N=size(g); %获取图像大小nn=2; %二健巴特沃斯高通滤波器d0=4; %截止频率为 4 m=fix(M/2);n=fix(N/2); for i=1:M for j=1:N d=sqrt(i-m)2+(j-n)2); %计算低通滤波器的传递函数if d=d0 h=0; else h=1; end result(i,j)=h*g(i,j)
15、; end end result=ifftshift(result); %转换数据矩阵J=ifft2(result); %傅里叶反转换K=uint8(real(J); %转换数据类型subplot(1,2,1),imshow(I),title( 原图 ) %显示图像subplot(1,2,2),imshow(K);title(高通滤波后的图像) 其相应的对比图像如下：共17 页第9 页（2）通过低通滤波器采用低通滤波器可以进行滤波，其相应的图像如下：I= imread (mri.tif); %读入图像 : f=double(I); %数据类型转换k=fft2(f); %傅里叶转换g=fftshift(k); %转换数据矩阵M,N=size(g); %获取图像大小nn=2; %二健巴特沃斯低通滤波器d0=10; %截止频率是 10 m=fix(M/2);n=fix(N/2); for i=1:
