图像处理实验参考程序

《图像处理实验参考程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验一matlab基础任务一：用记事本在系统work@录下建立一datl.txt数据文件，然后使用loaddatl.txt命令将其载入工作空间。任务二：_1111、输入矩阵Xl=[l122],X2=11,X3=,11JL」「123]X4=1232、计算各数组的维数及最长维数大小3、计算X3*X2（矩阵相乘）和X3.*X4（数组相乘），并比较二者的差异4、计算X3&X4任务三：编程计算1!+2!+3!+…+20!,并将程序保存为M文件参考程序：sum=0;for1=1:20prd=1;fork=l:iprd=prd*k;endsum二sum+pnl;end

2、sum实验二、对灰度图像作傅立叶变换并显示其频谱实验要求:K生成一个30X30的黑白矩形图像G,_II5

3、直方图实验要求:统计并显示图像pou.tif的灰度直方图，并与系统函数imhist的计算结果相比较I=imread(，pout,tif');[m,n]二size(I);x二class(I);switchxcase1uintS1max1=256;case'uintWmax1=65536;case'double"max1=1;endr=zeros(maxi,1);total=l;fori=l:mforj=l:nc二double(I(i,j))+l;r(c)=r(c)+l/total;endendx=1:max1;subplot(2,2,1);imshow(I

4、);subplot(2,2,3);bar(x,r),axis([0maxi0max(r(:))])subplot(2,2,4);imhist(I),axis([0maxi0max(r(:))]);实验四.对灰度图象进行直方图均衡化并显示对比效果实验要求：对灰度图象pou.tif进行直方图均衡化，并与系统函数histeq的处理结果进行对比。I=imrcad(，pout・tif');[m,n]二size(I);x=class(I);switchxcase1uint8"max1=256;case'uint16’max1=65536;case'double'ma

5、x1=1;endr=zeros(max1,1);total=m*n;%total=l;fori=l:mforj=l:nc二doubled(i,j))+l;r(c)=r(c)+l/total;endencl%计算累积直方图summ(l)=r(1);forx=2:maxlsumin(x)=suitiiti(x-1)+r(x~l);encl%计算变换后的灰度值，得到映射关系z=round((max1-1).*sunun);%图像直方图均衡G二z(I+l);subplot(2,2,1);imshow(I);subplot(2,2,3);imshow(G,[025

6、5]),subplot(2,2,4);histeq(I);实验五、设计频域理想低通滤波器，并加以验证实验要求：K对图象cameraman.tif添加高斯噪声2、设计一个理想低通滤波器对1得到的含噪图象进行处理3、依次显示原始图象、含噪图象、滤波器特性和去噪图象%读入图像，并变换到频域I=imread('cameraman.tif);Il=imnoise(I/gaussian,0.02);fftl=fftshift(fft2(ll,256,256));[m,n]=size(I);%设计理想低通滤波器[fl,f2]=freqspace([m,n],'me

7、shgrid');Hd=zeros(m,n);d=sqrt(fl.A2+f2.A2)<0.5;Hd(d)=1;%频域滤波jl=ifft2(Hd.*fftI);j2=wcodemat(abs(jl),256/matl);subplot(2,2,1)Jmshow(I)subplot(2,2,2),imshow(Il)subplot(2,2,3),mesh(fl,f2,Hd)subplot(2,2,4),imshow(J2,[0255])f0.5-1-1