实验一 图像的二维离散傅立叶变换资料.docx

1、实验一 图像的二维离散傅立叶变换资料实验一 图像的二维离散傅立叶变换一、 实验目的掌握图像的二维离散傅立叶变换以及性质二、 实验要求 1) 建立输入图像，在6464的黑色图像矩阵的中心建立1616的白色矩形图像点阵，形成图像文件。对输入图像进行二维傅立叶变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上。2) 调整输入图像中白色矩形的位置，再进行变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上，比较变换结果。3) 调整输入图像中白色矩形的尺寸（4040，44），再进行变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上，比较变换结果。三、 实验仪器设备及软件HP D538、M

2、ATLAB四、 实验原理在二维情况下，定义 f(x,y)的傅立叶变换(u,v) ：它表明了空间频率成分与二维图像信号之间的相互关系对于我们要处理的实际二维图像，其傅氏变换一般是在频率域上有界的，亦即有用成分总是落在一定的频率域范围之内上述的频率域性质的依据在于：一是图像中景物的复杂性具有一定的限度，其中大部分内容是变化不大的区域完全像“雪花”点似的图像没有任何实际意义。二是人眼对空间复杂性（频率）的分辨率以及显示器的分辨能力都是具有一定限度。若实变量函数f(x)是绝对可积的，即：且F(u)是可积的，则傅立叶变换对一定存在。如果f(x)考虑为实函数，它的傅立叶变换通常是复数形式，即： 也可表为：

3、若二变量函数f(x,y) 是绝对可积的，即：且F(u,v)是可积的，则傅立叶变换对一定存在。二维函数的傅立叶谱，振幅谱相位谱和能量谱分别为：五、 实验步骤及程序%clear%原始图象f=zeros(64,64);%输入64*64的黑色图像矩阵f(25:40,25:40)=1;%建立16*16的白色矩行图像点阵figure(1);subplot(231),imshow(f);title(原始图像)%显示原图像F=fft2(f);%傅立叶变换subplot(232)imshow(abs(F);title(傅里叶变换图像);%显示傅里叶变换图像F2=fftshift(abs(F);%频谱中心化sub

4、plot(233);imshow(abs(F2);title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图x=1:64;y=1:64;subplot(234);mesh(abs(real(F);title(三维频谱图);%显示三维频谱图subplot(235)mesh(x,y,F2(x,y);title(FFT)2、调整输入图像中白色矩形的位置，再进行变换后的程序clear%原始图象f=zeros(64,64);%输入64*64的黑色图像矩阵f(47:63,47:63)=1;%建立16*16的白色矩行图像点阵figure(1);subplot(231),imshow(f);title(原始图

5、像)%显示原图像F=fft2(f);%傅立叶变换subplot(232)imshow(abs(F);title(傅里叶变换图像);%显示傅里叶变换图像F2=fftshift(abs(F);%频谱中心化subplot(233);imshow(abs(F2);title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图x=1:64;y=1:64;subplot(234);mesh(abs(real(F);title(三维频谱图);%显示三维频谱图subplot(235)mesh(x,y,F2(x,y);title(FFT)3、整输入图像中白色矩形的尺寸（4040，44），再进行变换的程序4040cl

6、ear%原始图象f=zeros(64,64);%输入64*64的黑色图像矩阵f(13:52,13:52)=1;%建立16*16的白色矩行图像点阵figure(1);subplot(231),imshow(f);title(原始图像)%显示原图像F=fft2(f);%傅立叶变换subplot(232)imshow(abs(F);title(傅里叶变换图像);%显示傅里叶变换图像F2=fftshift(abs(F);%频谱中心化subplot(233);imshow(abs(F2);title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图x=1:64;y=1:64;subplot(234);me

7、sh(abs(real(F);title(三维频谱图);%显示三维频谱图subplot(235)mesh(x,y,F2(x,y);title(FFT)44clear%原始图象f=zeros(64,64);%输入64*64的黑色图像矩阵f(13:52,13:52)=1;%建立16*16的白色矩行图像点阵figure(1);subplot(231),imshow(f);title(原始图像)%显示原图像F=fft2(f);%傅立叶变换subplot(232)imshow(abs(F);title(傅里叶变换图像);%显示傅里叶变换图像F2=fftshift(abs(F);%频谱中心化subplot

8、(233);imshow(abs(F2);title(中心化傅里叶频谱图);%显示中心化傅里叶频谱图x=1:64;y=1:64;subplot(234);mesh(abs(real(F);title(三维频谱图);%显示三维频谱图subplot(235)mesh(x,y,F2(x,y);title(FFT)六、 实验结果与分析图1.1将原始图像及变换图像都显示的实验图像图1.2调整输入图像中白色矩形的位置，再进行变换后的实验图像图1.3调整输入图像中白色矩形的尺寸（4040），再进行变换的实验图像幅度谱决定了一幅图像中含有的各种频率分量的多少相位谱决定了每一种频率分量在图像中的位置。 只要每一

9、种频率分量保持在图像中的正确位置，那么图像的完整性就能得到很好的保持，这也就是为什么在信号或图像处理中通常只对幅度谱进行处理的原因。实验二 图像的增强一、 实验目的1) 掌握在计算机上进行直方图统计，以及直方图均衡化、线性变换的图像增强的方法 2) 掌握在计算机上进行图象平滑、图象锐化特别是中值滤波平滑及拉普拉斯算子锐化的方法 二、 实验要求 1) 显示图像（cameraman.tif）及灰度直方图。2) 对指定图像（cameraman.tif）进行直方图均衡化和线性变换，将原始图像及增强后的图像都显示于屏幕上，比较增强的效果。3) 对指定图像（lena.bmp）加入椒盐噪声，然后进行邻域平滑

10、、中值滤波，将原始图像及平滑后的图像都显示于屏幕上，比较效果。4) 对指定图像（lena.bmp）进行锐化（简单梯度算法、ROBERT算子，Prewitt边缘算子和拉普拉斯算子），将原始图像及锐化后的图像都显示于屏幕上，比较锐化的效果。三、 实验仪器设备及软件HP D538、MATLAB四、 实验原理是一种逐像素点对图像进行变换增强，也称为图像的点运算（一对一变换）。g(x,y) = T f(x,y) 度变换可以选择不同的灰度变换函数，如正比函数和指数函数等 。常用的灰度变换函数主要有： 线性灰度变换 分段线性灰度变换 非线性灰度变换将输入图像（原始图像）灰度值的动态范围按线性关系公式拉伸扩展

11、至指定范围或整个动态范围。线性拉伸采用的变换公式一般为:g(x,y)=f(x,y).C+R C、R的值由输出图像的灰度值动态范围决定。假定原始输入图像的灰度取值范围为fmin, fmax，输出 图像的灰度取值范围gmin, gmax，其变换公式为线性拉伸示意图线性拉伸是将原始输入图像中的灰度值不加区别地 扩展。 而在实际应用中,为了突出图像中感兴趣的研究对象，常常要求局部扩展拉伸某一范围的灰度值，或对不同范围的灰度值进行不同的拉伸处理，即分段线性拉伸。 分段线性拉伸是仅将某一范围的灰度值进行拉伸，而其余范围的灰度值实际上被压缩了常用的几种分段线性拉伸的示意图 :其对应的变换公式如下:给定一个二

12、维连续函数f(x,y) ，它在点(x,y)的梯度是一个矢量，定义为：梯度的重要特点： A. 矢量gradf(x,y)指向f(x,y)的最大增长率的方向。 B. 矢量gradf(x,y)的幅度和方向为：对于离散图像处理f(i,j)，常用到梯度的大小（幅度），因此把梯度的大小习惯称为“梯度”。并且一阶偏导数采用一阶差分近似表示，即 fx =f(i +1，j)-f(i，j) fy=f(i，j +1)-f(i，j)根据梯度计算式就可以计算Roberts、Prewitt和Sobel梯度。一旦梯度算出后，就可根据不同的需要生成不同的梯度增强图像。锐化算子的模板实质上可以看做是一种高通滤波器。因此，图像锐化

13、在增强图像边界和细节的同时，也加强了图像中的噪声。进行图像锐化处理的图像应有较高的信噪比，否则，锐化后信噪比会更低五、 实验步骤及程序（1）I=imread(cameraman.tif); %读入图像subplot(2,2,1);imshow(I);title(cameraman); %显示图像，并命名为cameramansubplot(2,2,2);imhist(I);title(直方图); %显示灰度直方图subplot(2,2,3);J=histeq(I);imshow(J);title(直方图均衡化); %灰度直方图均衡化subplot(2,2,4);K=imadjust(I,0.4,

14、0.6,); %进行线性变化imshow(K);title(线性变换) %显示线性变化的结果（2）I=imread(lena.bmp); %读入图像J=imnoise(I,salt & pepper,0.04); subplot(2,2,1);imshow(I); %显示源图像title(lena);subplot(2,2,2);imshow(J); %显示加入椒盐躁声的图像title(加入椒盐躁声);K=medfilt2(J); %对图像进行中值滤波subplot(2,2,3);imshow(K); %显示中值滤波后的图像title(中值滤波);subplot(2,2,4);H1=ones(3)/9;G=conv2(J,H1); %对图像进行领域平滑imshow(G,); %显示领域平滑后的图形title(邻域平滑)（3）I = imread(lena.bmp); %读入lena图像subplot(321)imshow(I) %显示源图像title(lena) %标题为lenasubplot(322)I=double(I);Gx,Gy=gradient(I);G=sqrt(Gx.*Gx+Gy.*Gy); %利用简单梯度法进行图像锐化O=G