数字图像处理实验报告(六)
以下为《数字图像处理实验报告(六)》的无排版文字预览，完整格式请下载
下载前请仔细阅读文字预览以及下方图片预览。图片预览是什么样的，下载的文档就是什么样的。
数字图像处理实验报告(六)
学院名称
信息电子技术
专业名称
通信工程
年级/班级
一
姓 名
张某某
学 号
***
指导教师
杜某某
实验日期
2021-04-29
实验地点
1C06-327
实验成绩
实验项目名称
图像频域处理
实
验
目
的
了解图像频域分析的目的、意义和方法；
掌握傅立叶变换和离散余弦变换在图像处理中的应用；
熟悉图像中各种频率分量的分布；
掌握频域处理的编程技巧。
实
验
内
容
2-1）、编程实现对“rec1.bmp”图像的频域处理，在一个窗口下显示原图像、频谱图像、中心平移后的频谱图像和中心平移后的增强型频谱图像；
参考函数 fft2、fftshift、abs、log
2-2）、编程实现对“rec1.bmp”和“rec2.bmp”图像的频域处理，在一个窗口下显示两幅原图像及其对应的中心平移后的增强型频谱图像，分析研究傅立叶变换的旋转不变性；
参考函数 fft2、fftshift、abs、log
2-3）、编程实现对“lina.bmp”图像的频域处理，然后在其中心平移后的频谱图像中截取出中心120×120区域，进行反变换，在一个窗口下显示原图、中心平移频谱图、截取区域图和其反变换效果图；
参考函数 fft2、fftshift、abs、log、ifft2
2-4）、编程实现对“lina.bmp”图像的频域处理，然后将其中心平移后的频谱图像中心120×120区域设为0，进行反变换，在一个窗口下显示原图、中心平移频谱图、去低频频谱图和其反变换效果图；
参考函数 fft2、fftshift、abs、log、ifft2
2-5）、以图像“text.png”为例，编程实现字符“a”在图片中的定位，分别在不同窗口下显示原图、字符a、及定位后的效果图效果图；
参考函数 fft2、ifft2、rot90、
（提示：字符a在图片中的(32:45,88:98)区域，可将其提取出作为模板，然后通过两图的卷积实现定位；显示定位效果时要设定相关阈值，即保留卷积后的峰值部分）
2-6）、编程实现对“autumn.tif”图像的离散余弦变换，并将其中小于10的系数设为0，然后进行反变换，在一个窗口下显示原图、其灰度图、灰度图的余弦变换图和反变换效果图。
参考函数 rgb2gray、dct2、idct2
实
验
要
求
1）、上述实验要以给定图像为目标进行展开；
2）、所有程序采用函数调用方式实现，且程序用独立的M文件来保存，文件分别为exp61.m～exp66.m；
3）、同一个窗口显示时，图片数量≤3时，以1行的方式进行显示，图片数量＞3时，以2行的方式进行显示；
4）、在关键语句和重要步骤上添加必要的注释，以辅助说明该语句和段落的功能。
实
验
程
序
五号宋体，单倍行距（逐个实验阐述相关实验程序，具体形式见下）
实验2-1程序如下：
A=imread('rec1.bmp');
N=100;
F=fft2(A,N,N);
F2=fftshift(abs(F));
F3=log(abs(F));
subplot(2,2,1);imshow(A);title('yuan')
subplot(2,2,2);imshow(F);title('FFT')
subplot(2,2,3);imshow(F2);title('pinyi')
subplot(2,2,4);imshow(F3);title('LOG')
实验2-2程序如下：
A=imread('rec1.bmp');
B=imread('rec2.bmp');
N=100;
F=fft2(A,N,N);
F2=fftshift(abs(F));
E=fft2(B,N,N);
E2=fftshift(abs(E));
subplot(2,2,1);imshow(A);title('Ayuan');
subplot(2,2,2);imshow(F2);title('Abian');
subplot(2,2,3);imshow(B);title('Byuan');
subplot(2,2,4);imshow(E2);title('Abian');
实验2-3程序如下：
A=imread('lina.bmp');
N=256;
a=double(A);
Z=fft2(a);
Z1=fftshift(Z);
B=Z1(68:188,68:188);
F=ifft2(B,N,N);
F1=uint8(abs(F));
subplot(2,2,1);imshow(A);title('原图');
subplot(2,2,2);imshow(Z1);title('中心平移频谱图');
subplot(2,2,3);imshow(B);title('截取区域图');
subplot(2,2,4);imshow(F1);title('反变换效果图');
实验2-4程序如下：
A=imread('lina.bmp');
N=256;
Z=fft2(A);
Z1=fftshift(Z);
Z3=fftshift(Z);
Z3(68:168,68:168)=0;
F=ifft2(Z1,N,N);
F1=uint8(abs(F));
subplot(2,2,1);imshow(A);title('原图');
subplot(2,2,2);imshow(Z1);title('中心平移频谱图');
subplot(2,2,3);imshow(Z3);title('去低频频谱图');
subplot(2,2,4);imshow(F1);title('反变换效果图');
实验2-5程序如下：
bw=imread('text.png');
a=bw(32:45,88:98);
B=real(ifft2(fft2(bw).*fft2(rot90(a,2),256,256)));
subplot(2,2,1);imshow(bw);
subplot(2,2,2);imshow(a);
subplot(2,2,3);imshow(B,[]);
thresh=60;
subplot(2,2,4);imshow(B>thresh);
实验2-6程序如下：
A=imread('autumn.tif');
H=rgb2gray(A);
C=dct2(H);
C(abs(C)请点击下方选择您需要的文档下载。
以上为《数字图像处理实验报告(六)》的无排版文字预览，完整格式请下载
下载前请仔细阅读上面文字预览以及下方图片预览。图片预览是什么样的，下载的文档就是什么样的。