[理学]数字图像处理指导-测控10版本.docVIP

实验一：图像预处理实验（验证性实验） 实验目的： 要求学生掌握图像读取、显示和保存的方法，掌握空间域增强方法和频率域增强方法。 实验要求： 上机运行，调试通过。 实验内容： 对附录一中的练习1~11、练习13，进行上机运行和调试。 对上述练习，改变一些函数的参数，观察运行结果的变化。 提交试验报告的内容： 实现对一幅图（比如：pout.tif）进行直方图均衡化的源代码； 显示直方图均衡化前后的图像及其直方图； 针对处理前后的图像和对应的直方图进行分析，得出直方图均衡化的作用。 实验二、图像的傅立叶变换（验证性实验） 实验目的 1了解图像变换的意义和手段； 2熟悉傅里叶变换的基本性质； 3熟练掌握FFT的方法与应用； 4理解频率域高低频与时域图像的对应关系； 5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换。 实验原理 1、应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间用在学习和掌握博里叶变换上是很有必要的。 2、傅立叶（Fourier）变换的定义 对于二维信号，二维Fourier变换定义为：二维离散傅立叶变换为： 图像的傅立叶变换与一维信号的傅立叶变换一样，有快速算法，具体参见参考书目，有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片，可以实时实现傅立叶变换。 3、用matlab编程观察图像频域中高低频率与时域图像的对应关系： %%---------------------------------------------------------- % 本代码主要演示图像fft2变换的频谱图中高频和低频分别对应 % 图像时域图中的哪些部分，加深对频域的理解 %----------------------------------------------------------- clear;close all; A=imread(lena.jpg); %读入图像 A=rgb2gray(A); figure(1);imshow(A); %显示图像 A=double(A); %%---------------------------------------------------------- % 对图像进行fft变换，看看频谱图和ifft变换的结果 % %----------------------------------------------------------- F0=fft2(A); %计算二维傅里叶变换 F1=fftshift(F0); %将直流分量移到频谱图的中心 figure(2), imshow(log(abs(F1)+1),[0 10]);%显示变换后的频谱图 F0=fftshift(F1); R0=ifft2(F0); %进行反变换 figure(3);imshow(uint8(R0)); %将变换的结果转换为uint8数据类型，并显示，比较反变换的结果与原图。 %%---------------------------------------------------------- % 把频谱中直流和低频系数变为0，只保留高频系数，然后反变化， % 看看高频系数对应图像时域中的哪一部分？ %----------------------------------------------------------- [m,n]=size(R0); %得到图像大小 i0=fix(m/2); %图像中心点的坐标 j0=fix(n/2); F2=F1; r=50; %设置系数改变区域的半径。 %把中心点附近2W*2W的低频全部设置为0； for i=1:m for j=1:n d=sqrt((i-i0)^2+(j-j0)^2); if(d=r) F2(i,j)=0; end; end; end; F3=fftshift(F2); %将直流分量移到频谱图的四个角 R1=ifft2(F3); %进行反变换 figure(4), imshow(log(abs(F2)+1),[0 10]);%显示去除低频的频谱图 figure(5);imshow(ui