图像处理计算机作业

《图像处理计算机作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第一章上机作业使用函数读取图像lena512_gray.bmp，存取到矩阵中clear;A=imread('lena512_gray.bmp');B=imread('lena512_color.bmp');C=rgb2gray(B);figure(1);subplot(1,2,1);imshow(A);title('读取灰度图像');subplot(1,2,2);imshow(B);title('读取彩色图像');figure(2);subplot(1,2,1);imshow(B);title('原彩色图像');subplot(1,2,2);i

2、mshow(C);title('把彩色图像转换为灰度图像');fori=1:256forj=1:256D(i,j)=A(i,j);endendfigure(3);subplot(1,2,1);imshow(A);title('显示原图像');subplot(1,2,2);imshow(D);title('显示图像的局域部分');………………………………………………………………………………………………………………………………………………………………………………………………………………第二章计算机作业(1)图像的采样%下采样A=imread('ba

3、rb1.bmp');[m,n]=size(A);B=zeros(m/2,n);fori=1:m/2B(i,1:end)=A(2*i,1:end);endC=zeros(m/2,n/2);forj=1:n/2C(1:end,j)=uint8(B(1:end,2*j));endD=uint8(C);imshow(D);%下采样A=imread('barb1.bmp');[m,n]=size(A);B=zeros(m/2,n);fori=1:m/2B(i,1:end)=A(2*i,1:end);endC=zeros(m/2,n/2);forj=1:n

4、/2C(1:end,j)=uint8(B(1:end,2*j));endD=uint8(C);subplot(1,2,1);imshow(D);%上采样[m,n]=size(D);F=zeros(2*m,n);fori=1:2*m-1F(i,1:end)=D(floor(i/2)+1,1:end);endG=zeros(2*m,2*n);forj=1:2*n-1G(1:end,j)=uint8(F(1:end,floor(j/2)+1));endE=uint8(G);subplot(1,2,2);imshow(E);>>B=imresize(A

5、,[row/2,column/2]);>>imshow(B);>>A=imread('barb1.bmp');>>[row,column]=size(A);>>B=imresize(A,[row/2,column/2]);>>imshow(B);（2）>>G=im2bw(A,0.4);>>imshow(G);>>G=im2bw(A,0.48);>>imshow(G);>>G=im2bw(A,0.5);>>imshow(G);>>G=im2bw(A,0.55);>>imshow(G);对比可知阈值设置在0.48左右比较合适。(3)>>A=imre

6、ad('barb2.bmp');>>noise=uint8(randn(size(A)));>>B=A+noise;>>imshow(B);第三章上机作业(1)图像的直方图I=imread('malight.bmp');m1=zeros(1,256);[m,n]=size(I);fori=1:256forj=1:256m1(I(i,j))=m1(I(i,j))+1;endendx=1:256;y=m1;subplot(1,2,1);imshow(I);title('原图像');subplot(1,2,2);bar(x,y,'g');title(

7、'直方图');调用imhist的效果更好一些。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。(1)图像的直方图均衡I=imread('malight.bmp');m1=zeros(1,256);[m,n]=size(I);fori=1:256forj=1:256m1(I(i,j))=m1(I(i,j))+1;endendm2=zeros(1,256);fori=1:256m2(i)=m1(i)./(m*n);endm3=zeros(1,2

8、56);fori=1:256forj=1:im3(i)=floor(m2(j)*255+0.5)+m3(i);endendm4=zeros(m,n);