MATALB初步与数字图像处理基础

年级学院专业班级学生姓名学号《数字图像处理》实验报告(1)实验题目数字图像处理初步与数字图像处理基础实验时间2014年11月18日实验地点实验成绩实验性质√应用性√设计性□综合性教师评阅：□实验目的明确；□操作步骤正确；□设计文稿（表格、程序、数据库、网页）符合要求；□保存路径正确；□实验结果正确；□实验分析总结全面；□实验报告规范；□其他：评阅教师签名：一、实验目的1.熟悉MATLAB的界面；2.熟悉MATLAB图像处理工具箱；3．掌握图像矩阵的显示及提取。二、实验主要内容与要求1.对lena图像分别进行4和16倍减采样，查看其减采样效果；2.显示一幅256灰度级的图像，其后转换为64级灰度图像、32级灰度图像、8级灰度图像和2级灰度图像。3.对图像进行模糊并使用维纳滤波器对图像进行复原。4.记录和整理实验报告。三、程序清单（计算机打印出的程序清单）1.a=imread('C:\MATLAB6p5\work\Lena.bmp');imshow(a)[wid,hei]=size(a);%4倍采样quartimg=zeros(wid/2+1,hei/2+1);i1=1;j1=1；fori=1:2:widforj=1:2:heiquartimg(i1,j1)=a(i,j);j1=j1+1;endi1=i1+1;j1=1;endfigureimshow(uint8(quartimg))%16倍采样quartimg=zeros(wid/4+1,hei/4+1);i1=1;j1=1;fori=1:4:widforj=1:4:heiquartimg(i1,j1)=a(i,j);j1=j1+1;endi1=i1+1;j1=1;endimshow(uint8(quartimg))2.a=imread('C:\MATLAB6p5\work\BaboonRGB.bmp');b=rgb2gray(a);figureimshow(b)[wid,hei]=size(b);img64=zeros(wid,hei);img32=zeros(wid,hei);img8=zeros(wid,hei);img2=zeros(wid,hei);%64级灰度fori=1:widforj=1:heiimg64(i,j)=floor(b(i,j)/4);endendfigureimshow(uint8(img64),[0,63])%32级灰度fori=1:widforj=1:heiimg32(i,j)=floor(b(i,j)/8);endendfigureimshow(uint8(img32),[0,31])%8级灰度fori=1:widforj=1:heiimg8(i,j)=floor(b(i,j)/32);endendfigureimshow(uint8(img8),[0,7])%2级灰度fori=1:widforj=1:heiimg2(i,j)=floor(b(i,j)/128);endendfigureimshow(uint8(img2),[0,2])3.a=imread('C:\DocumentsandSettings\maths1\桌面\dip-lb\用于数字图像处理的标准测试影像\数字图像处理标准测试图\PeppersRGB.bmp');imshow(a);%对图像进行运动模糊len=31;theta=11;psf1=fspecial('motion',len,theta);blurred1=imfilter(a,psf1,'circular','conv');subplot(1,2,1);imshow(blurred1);title('blurred1--motion');%对模糊图像加噪声v=0.002;blurredNoisy1=imnoise(blurred1,'gaussian',0,v);figure;subplot(1,2,1);imshow(blurredNoisy1);title('blurredNoisy1');%进行维纳滤波wnr1=deconvwnr(blurred1,psf1);%维纳滤波figure;subplot(1,2,1);imshow(wnr1);title('restored1,Truepsf');四、运行结果（必须是上面的程序清单所对应打印输出的结果）1.2.原图64级灰度图32级灰度图8级灰度图2级灰度图3.原图：模糊：噪声：维纳滤波：五、实验总结通过实验，我熟悉了MATLAB的操作和基于MATLAB进行图像处理的基础。六、思考题如何将一幅彩色图像转换成灰度图像?答：使用rgb2gray（）函数。