图像处理实验报告解读.doc

实 验 报 告 实验课程名称： 数字图像处理 班级： 学号： 姓名： 实验名称 报告规范程度（10分） 实验过程与结果（10分） 实验成绩 实验一图的二维傅立叶变换实验图的增强2016年 4 月18日 实验一 图的二维傅立叶变换二维傅立叶变换建立输入图，在6464的黑色图矩阵的中心建立1616的白色矩形图点阵，形成图文件。对输入图进行二维傅立叶变换，将原始图及变换图都显示于屏幕上。调整输入图中白色矩形的，再进行变换，将原始图及变换图都显示于屏幕上比较变换结果调整输入图中白色矩形的尺寸40，44），再进行变换，将原始图及变换图都显示于屏幕上比较变换结果。 …，M-1；y=0，1，…N-1 …，M-1；y=0，1，…N-1 在图像处理中，有事为了讨论上的方便，取M=N，这样二维离散傅里叶变换对就定义为 …，N-1 ，…，N-1 其中，是正变换核，是反变换核。将二维离散傅里叶变换的频谱的平方定义为的功率谱，记为 功率谱反映了二维离散信号的能量在空间频率域上的分布情况。 实验步骤、程序及结果： 1、实验步骤： （1）、编写程序建立输入图像； （2）、对上述图像进行二维傅立叶变换，观察其频谱 （3）、改变输入图像中白框的位置，在进行二维傅里叶变换，观察频谱； （4）、改变输入图像中白框的大小，进行二维傅里叶变换，观察频谱。 2、实验程序： （1）、原始图像及变换程序： %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));%频谱中心化 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) （2）、调整输入图像中白色矩形的位置，再进行变换后的程序 %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(47:63,47:63)=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); mesh(abs(real(F)));title(三维频谱图);%显示三维频谱图 subplot(235); mesh(x,y,F2(x,y)); title(FFT) （3）、整输入图像中白色矩形的尺寸（40*40），再进行变换的程序： clear %原始图象 f=zeros(64,64);%输入64*64的黑色图像矩阵 f(13:52,13:52)=1;%建立40*40的白色矩行图像点阵 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); mesh(abs(real(F)));title(三维频谱图);%显示三维频谱图 subplot(235); mesh(x,y,F2(x,y)); title(FFT) （4）、整输入图像中白色矩形的尺寸（4*4），再进行变换的程序： clear %原