资源描述：
《基于数字图像处理上机的实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、实验一一、实验目的1．熟悉有关图像增强的直方图方法；2．掌握傅立叶变换和离散余弦变换等常用的正交变换函数；二、实验内容1．实现对图像IM2.BMP进行直方图均衡化，从而达到增强图像的目的，观察图像直方图的前后变化；程序如下：I=imread('D:IM2.BMP');[I1,T]=histeq(I);imshow(I1);figure,imhist(I1);figure,plot((0:255)/255,T);结果如下：2．实现对图像IM3.BMP加入某种噪声（如椒盐噪声等），运用低通滤波和中值滤波分别进行图像平滑；程序如下：clear;closeall;I=imread('D:IM3.
2、BMP');imshow(I);I1=imnoise(I,'salt&pepper',0.02);figure,imshow(I1);f=double(I1);g=fft2(f);g=fftshift(g);[M,N]=size(g);d0=50;nn=2;m=round(M/2);n=round(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);h=1/(1+0.414*(d/d0)^(2*nn));result(i,j)=h*g(i,j);endendresult=ifftshift(result);I2=ifft2(result);I3=uin
3、t8(real(I2));figure,imshow(I3);结果如下：clear;closeall;I=imread('D:IM3.BMP');imshow(I);I1=imnoise(I,'salt&pepper',0.02);figure,imshow(I1);I2=double(I1);[M,N]=size(I2);I3=ones(size(I2));fori=2:M-1forj=2:N-1I3(i,j)=median([I2(i-1,j-1)I2(i-1,j)I2(i-1,j+1)...I2(i,j-1)I2(i,j)I2(i,j+1)...I2(i+1,j+1)I2(i+1,
4、j)I2(i+1,j+1)]);endendfori=2:M-1I3(i,1)=I3(i,2);I3(i,N)=I3(i,N-1);endI3(1,:)=I3(2,:);I3(M,:)=I3(M-1,:);figure,imshow(uint8(I3));3．实现对图像IM4.BMP的同态增强：对照明分量和反射分量进行同态滤波。程序如下：I=imread('D:IM5.BMP');S=qtdecomp(I,0.27);J=full(S);figure,imshow(J);blocks=repmat(uint8(0),size(J));fordim=[5122561286432168421
5、]numblocks=length(find(S==dim));if(numblocks>0)values=repmat(uint8(255),[dimdimnumblocks]);values(2:dim,2:dim,:);blocks=qtsetblk(blocks,S,dim,values);endendfigure,imshow(blocks);实验二图像处理一、实验目的1．掌握图像分割中的四叉树区域分割与合并法；2．掌握图像分割中的阈值分割法；3．了解图像纹理分析和特征提取方法；4．了解图像目标识别方法及其应用；二、实验内容1．使用四叉树(quadtree)的方法，找出图像IM5
6、.BMP中物体的轮廓；2．选择一个合适的初始阈值，将图像IM5.BMP分为两个区域；程序如下：I0=imread('D:IM5.BMP');imshow(I0);I=double(I0);T=(min(I(:))+max(I(:)))/2;done=false;while~doneI1=I>=T;T1=(mean(I(I1))+mean(I(~I1)))/2;done=abs(T-T1)<0.5;T=T1;endfigure,imshow(I1);figure,imhist(I0);axis([025501000]);结果如下：阈值分割法二clc;clearall;f=imread('D
7、:IM5.bmp');f=im2double(f);T=0.5*(min(f(:))+max(f(:)));done=false;while~doneg=f>=T;Tn=0.5*(mean(f(g))+mean(f(~g)));done=abs(T-Tn)<0.1;T=Tn;endTr=im2bw(f,T);figure,imshow(f),title('原图');figure,imshow(r),title('迭代法
