数字图像处理实验指导书(2014.10.13).doc

数字图像处理 实验指导书（自编） 实验1 空域增强--直接灰度映射 实验2 空域增强—平滑滤波器（邻域平均 vs 中值滤波） 实验3 图像变换--傅里叶变换及离散余弦变换原理及应用 实验4 频域增强—理想低通滤器vs Butterworth低通滤波 实验5 图像恢复—逆滤波 vs 维纳滤波 实验6 图像恢复--交互式恢复去除周期性噪声 实验7 图像编码-算术编码方法 实验8 图像编码--维游程编码方法 实验1 空域增强--直接灰度映射 授课教师 肖红 授课班级 信息11 1,2班 学时 理论48学时+实验16学时 实验环境 Matlab7.1 实验内容： 1．熟悉matlab工作界面 2．直接灰度映射的典型方法 ①图像求反 ②增强对比度 ③动态范围压缩 ④灰度切分 3.注意观察几种方法对图像进行处理后的对比效果 1．直接灰度映射的典型方法 ①图像求反 clear; close all; X1=imread(cameraman.tif); figure,imshow(X1); f0=0;g0=255; f1=255;g1=0 figure,plot([f0,f1],[g0,g1]) axis tight,xlabel(f),ylabel(g) title(intensitytransformation) %绘制变换曲线 r1=(g1-g0)/(f1-f0);%等于-1 b1=g0-r1*f0;%等于255 [m,n]=size(X1); X2=double(X1); for i=1:m for j=1:n f=X2(i,j); g(i,j)=0; g(i,j)=r1*f+b1; %这其实就是用255去减f，255-f end end figure, imshow(mat2gray(g)) figure,imshow(imcomplement(X1)),title(利用imcomplement对图像求反的效果)（1）原始图像 （2）映射曲线 （3）求反后效果 图1-1 图像求反示意图 ②增强对比度 clear; close all; X1=imread(cameraman.tif); figure, imshow(X1) f0=0;g0=0; f1=60;g1=50; f2=150;g2=240; f3=255;g3=255; figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3]) axis tight,xlabel(f),ylabel(g) title(intensitytransformation)%绘制变换曲线 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=f0)(ff1) g(i,j)=r1*f+b1; elseif(f=f1)(f=f2) g(i,j)=r2*f+b2; elseif(f=f2)(f=f3) g(i,j)=r3*f+b3; end end end figure,imshow(mat2gray(g)) （1）原始图像 （2）映射曲线 （3）增强后效果 图1-2 增强对比度示意图 图1-3 imadjdemo调整对比度 ③动态范围压缩 clear close all N=100 f=zeros(50,50); f(15:35,23:28)=1; %这样是一个长方形，很适合观察傅里叶变换的特征 F=fft2(f,N,N); F2=fftshift(abs(F)); %再以灰度图的形式显示一下 figure; subplot(121); imshow(F2,[]);%直接显示的话，图像很黑，看不清楚 subplot(122); imshow(0.25*log10(1+F2),[]);%这条语句用于压缩傅里叶频谱的动态范围。 %因为高频区域衰减的很快，导致频谱中低亮度值的可视细节丢失。通过计算对数，频谱的动态范围会大大减小，便于观察和处理。 %注：此处求的是log10，即以10为底的(1+F2)的对数。 图1-4 傅立叶频谱动态范围压缩 ④灰度切分 clear; close all; X1=imr