资源描述：
《数字图像处理试验程序清单》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1、灰度阈值分割i=imread('frog512.bmp');subplot(2,2,1),imshow(i),title('原始图像')ic=imcomplement(i);bw=im2bw(ic,graythresh(ic));subplot(2,2,2),imshow(bw),title('阈值截取分割后图像')se=strel('disk',6);bwc=imclose(bw,se);bwco=imopen(bwc,se);subplot(2,2,3),imshow(bwco),title('对小
2、图像进行删除后图像')mask=bw&bwco;subplot(2,2,4),imshow(mask),title('检测结果的图像')2、直方图分割clear;I=imread('Boat.bmp');imshow(I);figure;imhist(I);Inew=im2bw(I,140/255);figure;imshow(Inew);3、全局阈值分割%MATLAB对图像基于全局阈值的分割程序%functiontest2(a)clearall;a=150;I=imread('cameraman.bmp'
3、);figure(1),imshow(I);fori=1:256forj=1:256ifdouble(I(i,j))>aI(i,j)=255;endifdouble(I(i,j))<=aI(i,j)=0;endendendfigure(2),imshow(I)4.调用qtdecomp函数实现图像的四叉数分解i=imread('girl.bmp');s=qtdecomp(i,27);sblocks=repmat(uint8(0),size(s));fordim=[5122561286432168421];nu
4、mblocks=length(find(s==dim));if(numblocks>0)values=repmat(uint8(1),[dimdimnumblocks]);values(2:dim,2:dim,:)=0;blocks=qtsetblk(blocks,s,dim,values);endendblocks(end,1:end)=1;blocks(1:end,end)=1;imshow(i);figure,imshow(blocks,[])5、最大熵阈值clear;closeall;I=imrea
5、d('256gray.BMP');figure,imshow(I);iflength(size(I))==3I=rgb2gray(I);end[X,Y]=size(I);V_max=max(max(I));V_min=min(min(I));T0=(V_max+V_min)/2;h=imhist(I);figure,plot(h);grayp=imhist(I)/numel(I);I=double(I);H0=-sum(grayp(find(grayp(1:end)>0)).*log(grayp(find(
6、grayp(1:end)>0))));cout=100;while(cout>0)Tmax=0;grayPd=0;grayPb=0;Hd=0;Hb=0;T1=T0;A1=0;A2=0;B1=0;B2=0;fori=1:Xforj=1:Yif(I(i,j)<=T1)A1=A1+1;B1=B1+I(i,j);elseA2=A2+1;B2=B2+I(i,j);endendendM1=B1/A1;M2=B2/A2;T2=(M1+M2)/2;TT=round(T2);grayPd=sum(grayp(1:TT));i
7、fgrayPd==0grayPd=eps;endgrayPb=1-grayPd;ifgrayPb==0grayPb=eps;endHd=-sum((grayp(find(grayp(1:TT)>0))/grayPd).*log((grayp(find(grayp(1:TT)>0))/grayPd)));Hb=-sum(grayp(TT+(find(grayp(TT+1:end)>0)))/grayPb.*log(grayp(TT+(find(grayp(TT+1:end)>0)))/grayPb));H1=
8、Hd+Hb;cout=cout-1;if(abs(H0-H1)<0.0001)
9、(cout==0)Tmax=T2;break;elseT0=T2;H0=H1;endendTmaxcoutfori=1:Xforj=1:Yif(I(i,j)<=Tmax)I(i,j)=0;elseI(i,j)=1;endendendfigure,imshow(I);6.编程实现逆滤波、维纳滤波和约束最小二乘滤波复原。I=imread(
