数字图像处理实验结果.doc

1、 注意：打印后实验三第一个代码要解释代码试验一1.1 灰度图像到索引图像转换RGB=imread(2.jpg); figure(3); imshow(RGB); figure(1); RGB1,map1=gray2ind(RGB,128);imshow(RGB1,map1) ; figure(2); RGB2,map2=gray2ind(RGB,16); imshow(RGB2,map2) ; imwrite(RGB1,map1,3.bmp); imwrite(RGB2,map2,4.bmp); info=imfinfo(3.bmp)中中2.jpg中中中中1 3.bmp中中中中2 4.bmp1

2、.1.2 索引到灰度转换参考代码X,map=imread(3.bmp); figure(1); imshow(X); figure(2); I=ind2gray(X,map);imshow(I) ; imwrite(I,new.jpg); info=imfinfo(new.jpg)中中中中中中1.2 RGB 图像到索引图像的转换RGB=imread(1.jpg); figure(1); imshow(RGB); figure(2); y=rgb2ind(RGB,256); imwrite(y,1s.bmp); imshow(y); info=imfinfo(1s.bmp) RGB中中中中中中1

3、.2.2 索引转换 RGB参考代码及结果X,MAP=imread(1s.bmp); figure(1); imshow(X); figure(2); RGB=ind2rgb(X,MAP) imwrite(RGB,RGB.jpg); imshow(RGB); info=imfinfo(RGB.jpg)中中中中中中中中RGB1.3.1RGB 图像到二值图像的转换RGB=imread(1.jpg); % 将图像“1.jpg ”的数据赋给变量RGB figure(3); % 设定显示窗口 3 imshow(RGB); % 显示原“1.jpg ”文件 bw1=im2bw(RGB,0.5); % 阀值为

4、0.5 进行二值转换 bw2=im2bw(RGB,0.6); % 阀值为 0.6 进行二值转换 figure(1); % 设定显示窗口 1 imshow(bw1); % 显示索引图像 bw1 figure(2); % 设定显示窗口 2 imshow(bw2); % 显示索引图像 bw2 中RGB中中中中0.5中中中中中中0.6中中中中1.3.2 选取灰度图像进行二值转换参考代码如下（阈值要修改）RGB=imread(3.bmp); figure(3); imshow(RGB); bw1=im2bw(RGB,0.25); bw2=im2bw(RGB,0.35); figure(1); imsho

5、w(bw1); figure(2); imshow(bw2); 中中中中中中0.25中中中中中中0.35中中中中2.1 卷积实验A=magic(3); B=ones(3); A(8,8)=0; B(8,8)=0; C=ifft2(fft2(A).*fft2(B); C=C(1:5,1:5); C=real(C) C =8.0000 9.0000 15.0000 7.0000 6.000011.0000 17.0000 30.0000 19.0000 13.000015.0000 30.0000 45.0000 30.0000 15.00007.0000 21.0000 30.0000 23.0

6、000 9.00004.0000 13.0000 15.0000 11.0000 2.0000参考代码C=conv2(A,B); C=C(1:5,1:5) C =8 9 15 7 611 17 30 19 1315 30 45 30 157 21 30 23 94 13 15 11 22.2 对一副灰度图片进行二维傅立叶变换参考代码figure(1); RGB=imread(2.jpg); imshow(RGB); figure(2); B=fftshift(fft2(RGB); imshow(log(abs(B),),colormap(jet(64),colorbar;-5052.3创建大小

7、为 128*128,背景为黑色，中间开出一个 88 白色的窗口图像参考代码clearx=zeros(128,128);x(60:68,60:68)=1;imshow(x)imwrite(x,new.bmp)figure(1);RGB=imread(new.bmp);imshow(RGB);figure(2);B=fftshift(ifft2(RGB);imshow(log(abs(B),),colormap(jet(64),colorbar;3.1 进行时域灰度变换X1=imread(image.tif); figure,imshow(X1) f0=0;g0=0; f1=70;g1=30; f

8、2=180;g2=230; f3=255;g3=255; r1=(g1-g0)/(f1-f0); b1=g0-r1*f0; r2=(g2-g1)/(f2-f1); b2=g1-r2*f1; r3=(g3-g2)/(f3-f2); b3=g2-r3*f2; m,n=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; if(f=0) end end end figure,imshow(mat2gray(g) -10-503.2 实现 Butterworth 低通滤波器的参考代码I1=imread(3.bmp); fig

9、ure,imshow(I1) I2=imnoise(I1,salt); figure,imshow(I2) f=double(I2); g=fft2(f); g=fftshift(g); N1,N2=size(g); n=2; d0=50; n1=fix(N1/2); n2=fix(N2/2); for i=1:N1 for j=1:N2 d=sqrt(i-n1)2+(j-n2)2); h=1/(1+ (d/d0)(2*n); result(i,j)=h*g(i,j); end end result=ifftshift(result); X2=ifft2(result); X3=uint8(r

10、eal(X2); figure,imshow(X3) 4.1均值滤波3*3滤波代码I=imread(2.jpg); I=imnoise(I, gaussian,0,0.02); % 添加高斯噪声 figure(1) imshow(I); h=1 1 1 1 1 1 1 1 1 ; % 产生滤波模板, 不同模板对应的滤波形 式也不同 h=h/25; % 滤波模板归一化 J=conv2(I,h) ; % 滤波并显示图像 figure(2);imshow(J,) 5*5,滤波代码I=imread(2.jpg); I=imnoise(I, gaussian,0,0.02); % 添加高斯噪声 figu

11、re(1) imshow(I); h=1 1 1 1 11 1 1 1 11 1 1 1 11 1 1 1 11 1 1 1 1 ; % 产生滤波模板,不同模板对应的滤波形 式也不同 h=h/25; % 滤波模板归一化 J=conv2(I,h) ; % 滤波并显示图像 figure(2);imshow(J,) 7*7滤波参考代码I=imread(2.jpg); I=imnoise(I, gaussian,0,0.02); % 添加高斯噪声 figure(1) imshow(I); h=1 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1

12、 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 1 ; % 产生滤波模板,不同模板对应的滤波形 式也不同 h=h/49; % 滤波模板归一化 J=conv2(I,h) ; % 滤波并显示图像 figure(2);imshow(J,)中中中中3*3中中5*5 中中7*7中中4.2 中值滤波参考代码I=imread(2.jpg); % 将图像的数据赋给变量 I J=imnoise(I,gaussian,0.02); % 加入噪声 K=medfilt2(J); % 进行中值滤波 subplot(211),imshow(J); % 分别显示原图像和变换后的图像 subplot(2

13、12),imshow(K) 中中中中中中中中中4.3加椒盐噪声后中值滤波与均值滤波对比参考代码I=imread(b.jpg); N=rgb2gray(I);subplot(2,2,1),imshow(N); title(灰度图)J=imnoise(J,salt K=medfilt2(J); subplot(2,2,2),imshow(J); title(加椒盐噪声图) subplot(2,2,3),imshow(K);title(中值滤波图)h= 1 1 1 1 1 1 1 1 1 ; h=h/9; M=conv2(N,h) ; subplot(2,2,4),imshow(M);title(均

14、值滤波图)处理图像结果中中中 中中中中中中中中中中中 中中中中中加高斯噪声后中值滤波与均值滤波对比参考代码I=imread(b.jpg); N=rgb2gray(I);subplot(2,2,1),imshow(N); title(灰度图)J=imnoise(J,gaussian,0,0.02); K=medfilt2(J); subplot(2,2,2),imshow(J); title(加高斯噪声图) subplot(2,2,3),imshow(K);title(中值滤波图)h=1 1 1 1 h=h/4; M=conv2(N,h) ; subplot(2,2,4),imshow(M);t

15、itle(均值滤波图)中 中 中 中 中 中 中 中 中中 中 中 中 中 中 中 中 中 中4.3 边缘检测I=imread(2.jpg); subplot(411),imshow(I);title(原图)BW1=edge(I,Roberts); subplot(412),imshow(BW1);title(Roberts)BW2=edge(I,sobel) ; subplot(413),imshow(BW2);title(sobel)BW3=edge(I,Prewitt)subplot(414),imshow(BW3);title(Prewitt)中中RobertssobelPrewitt