医学图像处理实验一指导书2013.pdfVIP

实验一：图像的傅里叶变换与空域滤波 实验学时：4 实验类型：设计 实验要求：必修 一、实验目的 通过本实验，掌握数字图像的表示、傅里叶变换的原理、快速傅里叶变换的 实现方法、图像的常见噪声类型及基本的空域滤波方法，在此基础上熟悉 MATLAB 编程语言，熟悉编程方法。 二、实验内容 1. 用 MATLAB 读取、显示图像 square.bmp 和 lena.jpg，并获取图像的全部 信息及矩阵参数； 2. 根据二维离散傅里叶变换（DFT ）公式和快速傅里叶变换（FFT ）算法分 别编写程序，计算图像 square.bmp 的傅里叶变换；用MATLAB 提供的二维 DFT 函数作相同的计算，比较计算结果、所需计算时间，并与理论分析作对比； 3.用 MABLAB 提供的函数分别给图像 lena.bmp 加入高斯噪声、椒盐噪声及 乘性噪声，采用均值滤波、中值滤波对加噪声后的图像进行处理，比较不同滤波 方法的处理结果，分析不同参数设置对滤波效果的影响。 三、实验原理、方法和手段 1. MATLAB 图像读取、显示、信息获取的基本函数 I=imread (filename,format); %读取图像 imshow(I) ； %显示灰度图像 whos I %显示出图像一个数组的附加信息 imfinfo (filename, format) %显示出图像的所有详细信息 [M N]=size(I); %确定一幅图像的行数和列数 No=numel(I); %确定一幅图像的像素数目，即矩阵（数组）的元素数目 Y1=fft(x); %计算一维快速傅里叶变换，x 为一维数组 Y2=fft2(I); %计算二维快速傅里叶变换，I 为图像矩阵 t=dec2bin(a,b); %将十进制数 a 转换成至少 b 位的二进制数 a=bin2dec(t);%将二进制数 t 转换成十进制数 c=bitshift(a,b)%将整数a 以二进制形式左移b位，输出整数c tic; toc; %分别为计时开始和计时结束语句，调用计时器测算其间的程序运行时间 J = imnoise(I,type,parameters); %给图像加入制定类型的噪声； H=fspecial(type,parameter) ；%产生指定类型的二维滤波器算子； A= imfilter(J,H) ；%用指定算子对图像进行滤波； B = medfilt2(A, [m n]) ；%用指定尺寸的模板对图像进行中值滤波； 1 2. 二维 DFT 公式与性质 设f (x , y ) ( x 0,1, , N 1; y 0,1, , N 1) 是一幅N N 图像，其二维离散傅里 叶变换F (u,v) 定义为： N 1 N 1 uxvy j 2 F (u,v) f (x , y )e N x 0 y 0 （1） u 0,1, , N 1; v 0,1, , N 1; 基于二维 DFT 的可分离性，二维 FFT 算法可用两个一维 FFT 算法来实现，如图 1 所 示 N-1 v N-1 y N-1 v 0 0 行变换 0 列变换