数字图像处理——实验.ppt

数字图像处理——实验;实验一、计算图象统计参数;实验二;实验二——参考答案;实验二——参考答案;实验三;i=imread(‘test1.tif); figure,imshow(i); title(‘test1.tif); f=fft2(i); f1=fftshift(f); rr=real(f1); ii=imag(f1); a1=abs(f1); a2=sqrt(rr.^2+ii^2); aa1=log(1+a1); aa2=log(1+a2); figure,imshow(aa1,[]); title(a1 abs); figure,imshow(aa2,[]); title(a2 sqrt);;实验三;A=zeros(99,99); A(49:51,49:51)=1; B=fft2(A); B1=fftshift(B); figure,imshow(A); title(冲击函数imsow); figure,mesh(A); title(冲击函数mesh); figure,imshow(log(1+abs(B1))); title(冲击函数频域模mesh); figure,mesh(log(1+abs(B1))) title(冲击函数频域模mesh);;实验三; g=imread(test2.tif); G=fft2(g); w=fspecial(gaussian,9,2); fi1=imfilter(g,w); W=freqz2(w,256,256); W1=ifftshift(W); Ff=G.*W1; fi2=real(ifft2(Ff)); figure,imshow(w); title(w); figure,mesh(W); title(W); figure,imshow(g); title(‘test2.tif(g)); figure,imshow(fi1,[]); title(空间滤波); figure,imshow(fi2,[]); title(频域滤波);实验四;实验四;实验五;实验五——参考答案;f=imread(‘strawberries_fullcolor.tif); [x1,map1]=rgb2ind(f,8,nodither); figure,imshow(x1,map1); title(ind nodither); [x2,map1]=rgb2ind(f,8,dither); figure,imshow(x2,map1); title(ind dither); g=rgb2gray(f); g1=dither(g); figure,imshow(g); title(gray nodither); figure,imshow(g1); title(gray dither); ;2、彩色图像平滑与锐化 对图像“iris.tif ” 在 RGB 空 间 实 现 彩色图像的平滑 （ w=ones（25）./(25*25)）与锐化（w= [1 1 1 1 1; 1 1 1 1 1; 1 1 -24 1 1 ; 1 1 1 1 1; 1 1 1 1 1];）（函数 imfilter）; 在 HSI 空间实现彩色图像的平滑（w=ones（25）./(25*25)），（函数 imfilter, rgb2hsi, cat），观察仅平滑亮度分量（intensity）和平滑全部三个分量结果的差 别。 ;RGB;i=imread(‘iris.tif ); figure,imshow(i); w1=ones(25)./(25*25); w2=[1 1 1 1 1;1 1 1 1 1; 1 1 -24 1 1;1 1 1 1 1; 1 1 1 1 1]; i1=imfilter(i,w1,replicate); i2=imfilter(i,w2,replicate); figure,imshow(i1); title(rgb平滑); figure,imshow(i2) title(rgbr锐化); (1);3、彩色图像分割 用 colorseg 函 数 实 现 基 于 肤 色 的 人 脸 分 割 。 对 图 像“woman_baby_original.tif”在 HSI 空间，用函数 colorseg 分别用参数“euclidean”和“mahalanobis”实现人脸分割。