[计算机]数字图像处理实验程序MATLAB.doc

实验一 内容（一） （1）彩色图像变灰度图像 A=imread(1.jpg); B=rgb2gray(A); figure subplot(1,2,1), imshow(A) title(原图) subplot(1,2,2), imshow(B) title(原图灰度图像) （2）彩色图像变索引图像 A=imread(1.jpg); figure subplot(1,2,1), imshow(A) title(原图) [X,map]=rgb2ind(A,128); subplot(1,2,2), imshow(X,map) title(原图索引图像) （3）彩色图像变二值图像 A=imread(1.jpg); figure subplot(1,2,1), imshow(A) title(原图) C=im2bw(A,0.2); subplot(1,2,2), imshow(C) title(原图二值图像) （4）灰度图像变索引图像（一） A=imread(1.jpg); figure B=rgb2gray(A); subplot(1,2,1), imshow(B) title(灰度图像) C=grayslice(B,39); subplot(1,2,2), imshow(C) title(灰度变索引图像) （5）灰度图像变索引图像（二） A=imread(1.jpg); figure B=rgb2gray(A); subplot(1,2,1), imshow(B) title(灰度图像) [X,map]=gray2ind(B,63); subplot(1,2,2), imshow(X,map) title(灰度变索引图像) （6）灰度图像变彩色图像 A=imread(1.jpg); figure B=rgb2gray(A); subplot(1,2,1), imshow(B) title(灰度图像) C=gray2rgb(B,map); subplot(1,2,2), imshow(C) title(灰度变彩色图像) 内容（二） （1）灰度平均值 A=imread(1.jpg); figure B=rgb2gray(A); subplot(1,2,1), imshow(B) title(灰度图像) B=double(B); [m,n]=size(B); sumg=0.0; for i=1:m; for j=1:n; sumg=sumg+B(i,j); end end avg=sumg/(m*n) % 均值 maxg=max(max(B)) % 区域最大灰度 ming=min(min(B)) % 区域最小灰度 （2）彩色平均值 A=imread(1.jpg); figure imshow(A) title(彩色图像) A=double(A); [m,n]=size(A); sumg=0.0; for i=1:m; for j=1:n; sumg=sumg+A(i,j); end end avg=sumg/(m*n) squre=m*n maxg=max(max(A)) ming=min(min(A)) 内容（三） 采样量化 实验二 图像变换 傅里叶变换、反变换、 I=imread(19.jpg); A=rgb2gray(I); x1=fft2(A); x2=fftshift(x1); x3=ifft(x1)/10; figure, subplot(1,4,1); imshow(A) title(原图); subplot(1,4,2); imshow(x1) title(频谱图); subplot(1,4,3); imshow(log(abs(x2)+1),[0 10]); title(直流分量移至频谱图中心); subplot(1,4,4); imshow(x3,[0 10]) title(傅里叶反变换); DCT变换、反变换 X=imread(19.jpg); I=rgb2gray(X); subplot(1,3,1); imshow(I); title(原图); subplot(1,3,2); J=dct2(I); imshow(log(abs(J)),[0 20]); title(二维离散余弦变换); subplot(1,3,3); K=idct2(J)/20; imshow(K,[0 20]); title(二维离散反余弦变换); 利用DCT变换压缩图像 I=imread(19.jpg); A=rgb2gray(I); B=DCT2(A); B(abs(B)0.1)=0; C=idct2(B)/255; figure, subplot(1,3,1)