计算机数字图像处理上机作业0524答案.docVIP

北京工商大学 计算机数字图像处理 上机任务书 1、读入图片文件map.png，用变量I1表示； 给图片添加高斯噪声，均值0，方差0.05，用变量I2表示； 用3行3列的均值滤波器对添加噪声后的图片滤波，用变量I3表示； 显示I1,I2,I3 I1=imread(map.png); %读入图片文件map.png，用变量I1表示 imshow(I1) %显示图像I1 I2=imnoise(I1,gaussian,0,0.05); %给图片添加高斯噪声，均值0，方差0.05，用变量I2表示 figure,imshow(I2) %重新打开一个窗口，显示图像I2 h=fspecial(average,3); %用3行3列的均值滤波器对图像进行滤波，生成新模板为h I3=imfilter(I2,h); %用模板h对图像I2进行滤波 figure,imshow(I3) %重新打开一个窗口，显示图像I3 2、读入图片文件lena.jpg，用变量I4表示； 显示I4的直方图； 用imadjust函数将图片调整范围内的灰度值映射到整个灰度空间（用stretchlim函数确定直方图调整范围），调整后的图片用I5表示。 显示I4和I5，以及直方图。 I4=imread(lena.jpg); %读入图片文件lena.jpg，用变量I4表示 imshow(I4) %显示图像I4 imfinfo(lena.jpg) %查找图像信息，lena.jpg为真彩色图像。需将其转化为灰度图像，才能显示直方图 ans = Filename: lena.jpg FileModDate: 20-May-2011 13:35:10 FileSize: 10904 Format: jpg FormatVersion: Width: 256 Height: 256 BitDepth: 24 ColorType: truecolor FormatSignature: Comment: {} I4h=rgb2gray(I4); %将真彩色I4图像转化为灰度图像，用I4h表示 figure,imhist(I4h,64) %重新打开一个窗口，显示I4h的直方图 I5=imadjust(I4h,stretchlim(I4h),[0 1]); %用imadjust函数将图片调整范围内的灰度值映射到整个灰度空间（用stretchlim函数确定直方图调整范围），调整后的图片用I5表示。 figure,imshow(I5) %重新打开一个窗口，显示图像I5 figure,imhist(I5,64) %重新打开一个窗口，显示I5的直方图 3、将I4放大1.5倍，采用双线性插值，放大后的图片用I6表示； 将I6的像素灰度值提高50，用I7表示； 将I7转化为二值图，阈值为0.4. I6=imresize(I4h,1.5,bilinear); %将I4放大1.5倍，采用双线性插值，放大后的图片用I6表示； figure,imshow(I6) %重新打开一个窗口，显示图像I6 I7=imadd(I6,50); %将I6的像素灰度值提高50，用I7表示 figure,imshow(I7) %重新打开一个窗口，显示图像I7 I8=im2bw(I7,0.4); %将I7转化为二值图，阈值为0.4，用I8表示 figure,imshow(I8) %重新打开一个窗口，显示图像I8 4、将真彩色图lena.jpg转化为灰度图，用J1表示； 读入图片rice.tif,用J2表示； 将J1和J2叠加，用J3表示; 读入图片moon.tif，剪切图片至大小[256 256]，用J4表示; 将J3和J4叠加，用J5表示； 将J5顺时针方向旋转60° 显示J1~J5。 J1=rgb2gray(I4); %将真彩色图lena.jpg转化为灰度图，用J1表示 imshow(J1) %显示图像J1 J2=imread(rice.tif); %读入图片rice.tif,用J2表示 figure,imshow(J2) %重新打开一个窗口，显示图像J2 J3=imadd(J1,J2); %将J1和J2叠加，用J3表示 figure,imshow(J3) %重新打开一个窗口，显示图像J3 J40=imread(moon.tif); %读入图片moon.tif，用J40表示 fi