数字图像处理 许录平 实验 MATLABPPT.pptVIP

邹媛;实验一 MATLAB编程基础练习;x=0:0.05:5; y=sin(x.^2) plot(x,y);;x=0:0.1:4; y=sin(x.^2).*exp(-x); stem(x,y);;x=-2.9:0.2:2.9; bar(x,exp(-x.*x)) ;surf(40*membrane(1,25)) shading interp daspect([1 1 0.9]) camlight axis off tight title(MATLAB.the language of technical computing.) ;z=peaks(25)mesh(z) ;z=peaks(25) surf(z) colormap(jet) ;z=peaks(25) surf(z) shading interp colormap(pink) ;clearclose allI=imread(rice.tif)imshow(I) ;Background=imopen(I,strel(disk,15)) I2=imsubtract(I,Background) figure,imshow(I2) I3=imadjust(I2,stretchlim(I2),[0 1]) figure,imshow(I3);t=0:0.1:2*pi y1=sin(t) h_line1=plot(t,y1,:) hold on y2=cos(t) h_line2=plot(t,y2,*) h_label=xlabel(Time(0~2\pi),FontWeigh,bold) h_text1=text(pi,0,\leftarrowsin wave) h_text2=text(pi/2,0,\leftarrowcos wave) ;实验二 中值滤波的MATLAB实现;I=imread(eight.tif) J=imnoise(I,salt pepper,0.02) K=medfilt2(J) imshow(J),figure,imshow(K) ;I=imread(moon.tif) I1=imnoise(I,salt pepper,0.06) I2=double(I1)/255 h1=[1/9 1/9 1/9;1/9 1/9 1/9;1/9 1/9 1/9]; J1=conv2(I2,h1,same) J2=medfilt2(I2,[3 3]) figure,imshow(I) figure,imshow(J1) figure,imshow(J2);;实验三 基于MATLAB的二值图像操作;load trees figure(1) imshow(X,map) ff0=im2bw(X,map,0.3) figure(2) imshow(ff0) ff1=im2bw(X,map,0.5) figure(3) imshow(ff1) ff2=im2bw(X,map,0.7) figure(4) imshow(ff2) ;;BW1=imread(circlesm.tif) SE=eye(5) BW2=erode(BW1,SE) BW3=dilate(BW1,SE) subplot(2,2,1),imshow(BW1),title(原始图像) subplot(2,2,2),imshow(BW2),title(腐蚀后图像) subplot(2,2,3),imshow(BW3),title(膨胀后图像) ;I=imread(circbw.tif)subplot(2,2,1),imshow(I),title(原始图像)I8=bwperim(I)subplot(2,2,2),imshow(I8),title(识别边界（8-连接边沿）)I4=bwperim(I,4)subplot(2,2,3),imshow(I4),title(识别边界（4-连接边沿）) ;A=imread(testpat2.tif) imshow(A) S=bwarea(A) E8=bweuler(A) E4=bweuler(A,4) ;I=imread(cell.tif) figure,imshow(I),title(original image) BWs=edge(I,sobel,(graythresh(I) * .1)) figure,imshow(BWs),title(binary gradient mask) se90=strel(line,3,90) se0=strel(line,3,0) BWsdil=imdilate(BWs,[se90 se0]) figure,imshow(BWsdil),title(dilated gradient mask) BWdfill=imf