图像处理实验报告.doc

《图像处理实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验报告实验课程名称：数字图像处理班级：学号：姓名：实验名称报告规范程度（10分）实验过程与结果（10分）实验成绩实验一：图像的二维离散傅立叶变换实验二：图像的增强实验三圆形物体的图像分割与分析平均成绩注：1、每个实验中各项成绩按照10分制评定，每个实验成绩为两项总和20分。2、平均成绩取三个实验平均成绩。2016年4月18日实验一图像的二维离散傅立叶变换一、实验目的掌握图像的二维离散傅立叶变换以及性质二、实验要求1)建立输入图像，在6464的黑色图像矩阵的中心建立1616的白色矩形图像点阵，形成图像文件。对输入图像进行二维傅立

2、叶变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上。2)调整输入图像中白色矩形的位置，再进行变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上，比较变换结果。3)调整输入图像中白色矩形的尺寸（4040，44），再进行变换，将原始图像及变换图像（三维、中心化）都显示于屏幕上，比较变换结果。三、实验仪器设备及软件HPD538、MATLAB四、实验原理傅里叶变换作为分析数字图像的有利工具，因其可分离性、平移性、周期性和共轭对称性可以定量地方分析数字化系统，并且变换后的图像使得时间域和频域间的联系能够方便直观地解决许多问题

3、。实验通过MATLAB实验该项技能。设是在空间域上等间隔采样得到的M×N的二维离散信号，x和y是离散实变量，和为离散频率变量，则二维离散傅里叶变换对一般地定义为…，M-1；y=0，1，…N-1…，M-1；y=0，1，…N-1在图像处理中，有事为了讨论上的方便，取M=N，这样二维离散傅里叶变换对就定义为…，N-1，…，N-1其中，是正变换核，是反变换核。将二维离散傅里叶变换的频谱的平方定义为的功率谱，记为功率谱反映了二维离散信号的能量在空间频率域上的分布情况。五、实验步骤、程序及结果：1、实验步骤：（1）、编写程序建立输入图像；（

4、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);im

5、show(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）、调整输入图像中白色矩形的位置，再进行变换后的程序%

6、原始图象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('中心化傅里叶频谱图');%显

7、示中心化傅里叶频谱图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);titl

8、e('原始图像')%显示原图像F=fft2(f);%傅立叶变换subplot(232);imshow(abs(F));title('傅里叶变换图像');%显示傅里叶变换图像F2=fftshift(abs(F));%频谱中心化subplot(233);imshow