数字图像处理实验程序MATLAB.

实验一

内容（一）

（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）彩色平均值

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变换、反变换

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);

imshow(A);

title('原图');

subplot(1,3,2);

imshow(B);

title('二维离散余弦变换频谱图');

subplot(1,3,3);

imshow(C);

title('压缩后图像');

实验三图像增强

（一）灰度图像增强

（1）线性变换法

clc;

clear all;

I=imread('19.jpg');

A=rgb2gray(I);

colormap;

imshow(A);%设置图像倒数参数

j=imadjust(A,[0 1],[1 0],1.5);

figure;subimage(j)