图像处理实验参考程序.docVIP

实验一 matlab基础 任务一：用 记事本在系统work目录下建立一dat1.txt数据文件，然后使用load dat1.txt命令将其载入工作空间。 任务二： 1、 输入矩阵，，， 2、 计算各数组的维数及最长维数大小 3、计算（矩阵相乘）和X3.*X4（数组相乘），并比较二者的差异 4、计算 任务三：编程计算1！+2！+3！+…+20！，并将程序保存为M文件 参考程序： sum=0; for I=1:20 prd=1; for k=1:i prd=prd*k; end sum=sum+prd; end sum 实验二、对灰度图像作傅立叶变换并显示其频谱 实验要求： 1、生成一个30×30的黑白矩形图像G， 2、计算并显示该图像的频谱 %傅立叶变换的图形表示 f=zeros(30,30); f(5:24,13:17)=1; F=fft2(f,256,256); F2=fftshift(F); F3=log(1+abs(F2)); subplot(2,2,1); imshow(f); title(原图象); subplot(2,2,2); imshow(F3); title(傅立叶幅度谱); 实验三、统计并显示图像的灰度直方图 实验要求： 统计并显示图像pou.tif的灰度直方图，并与系统函数 imhist的计算结果相比较 I=imread(pout.tif); [m,n]=size(I); x=class(I); switch x case uint8 max1=256; case uint16 max1=65536; case double max1=1; end r=zeros(max1,1); total=1; for i=1:m for j=1:n c=double(I(i,j))+1; r(c)=r(c)+1/total; end end x=1:max1; subplot(2,2,1); imshow(I); subplot(2,2,3); bar(x,r),axis([0 max1 0 max(r(:))]) subplot(2,2,4); imhist(I),axis([0 max1 0 max(r(:))]); 实验四、对灰度图象进行直方图均衡化并显示对比效果 实验要求：对灰度图象pou.tif进行直方图均衡化，并与系统函数histeq的处理结果进行对比。 I=imread(pout.tif); [m,n]=size(I); x=class(I); switch x case uint8 max1=256; case uint16 max1=65536; case double max1=1; end r=zeros(max1,1); total=m*n; %total=1; for i=1:m for j=1:n c=double(I(i,j))+1; r(c)=r(c)+1/total; end end %计算累积直方图 summ(1)=r(1); for x=2:max1 summ(x)=summ(x-1)+r(x-1); end %计算变换后的灰度值，得到映射关系 z=round((max1-1).*summ); %图像直方图均衡 G=z(I+1); subplot(2,2,1); imshow(I); subplot(2,2,3); imshow(G,[0 255]), subplot(2,2,4); histeq(I); 实验五、设计频域理想低通滤波器，并加以验证 实验要求： 对图象cameraman.tif添加高斯噪声 设计一个理想低通滤波器对1得到的含噪图象进行处理 依次显示原始图象、含噪图象、滤波器特性和去噪图象 %读入图像,并变换到频域 I=imread(cameraman.tif); I1=imnoise(I,gaussian,0,0.02); fftI=fftshift(fft2(I1,256,256)); [m,n]=size(I); %设计理想低通滤波器 [f1,f2]=freqspace([m,n