數(shù)字圖像處理試驗程序清單

1、灰度閾值分割i=imread('frog512.bmp');subplot(2,2,l),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(對小圖像進行刪除后圖像')mask=bw&bwco;subplot(2,2,4),imshow(mask),title(檢測結(jié)果的圖像')2、直方圖分割clear;I=imread('Boat.bmp');imshow(I);figure;imhist(I);Inew=im2bw(I,l40/255);figure;imshow(Inew);3、全局閾值分割%MATLAB對圖像基于全局閾值的分割程序%functiontest2(a)clearall;a=l50;I=imread('cameraman.bmp');figure(l),imshow(I);fori=l:256forj=l:256ifdouble(I(i,j))>aI(i,j)=255;endifdouble(I(i,j))<=aI(i,j)=0;endendendfigure(2),imshow(I)4、調(diào)用qtdecomp函數(shù)實現(xiàn)圖像的四叉數(shù)分解i=imread('girl.bmp');s=qtdecomp(i,27);sblocks=repmat(uint8(0),size(s));fordim=[5l2256l286432l6842l];numblocks=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=imread('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(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));ifgrayPd==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=Hd+Hb;cout=cout-1;if(abs(H0-H1)<0.0001)|(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.編程實現(xiàn)逆濾波、維納濾波和約束最小二乘濾波復原。I=imread('Girl.bmp');figure(1);imshow(I,[]);PSF=fspecial('motion',40,75);MF=imfilter(I,PSF,'circular');noise=imnoise(zeros(size(I)),'gaussian',0,0.01);MFN=imadd(MF,im2uint8(noise));figure(2);imshow(MFN,[]);NSR=sum(noise(:).人2)/sum(MFN(:).人2);figure(3);imshow(deconvwnr(MFN,PSF),[]);figure(4);imshow(deconvwnr(MFN,PSF,NSR),[]);NP=0.02*prod(size(I));[reg1LAGRA]=deconvreg(MFN,PSF,NP/3.0);figure(5);imshow(reg1);edged=edgetaper(MFN,PSF);reg2=deconvreg(edged,PSF,[],LAGRA/300);figure(6);imshow(reg2);1．實驗四4基于邊緣檢測的圖像分割方法。clearx=-2:0.1:2;y=-2:0.1:2;sigma=0.5;y=y;fori=1:(4/0.1+1)xx(i,:)=x;yy(:,i)=y;endr=1/(pi*sigmaA4)*((xx.A2+yy.A2)/(2*sigmaA2)-1)?*exp(-(xx.A2+yy.A2)/(2*sigmaA2));colormap(jet(16));mesh(xx,yy,r)I=imread('barbara?bmp');figure;imshow(I)bw=edge(I,'log');figure;imshow(bw)實驗五1建立由于運動所造成的退化圖像，并去除模糊I=imread('Girl.bmp');figure(1);imshow(I);LEN=30;THETA=75;PSF=fspecial('motion',LEN,THETA);MF=imfilter(I,PSF,'circular','conv');figure(2);imshow(MF);wnr1=deconvwnr(MF,PSF);figure(3);imshow(wnr1);3.利用deconvwnr函數(shù)對有噪聲模糊圖像進行復原重建。cleari=imread('girl,bmp');i=i(5+[1:256],40+[1:256],:);subplot(3,3,1),imshow(i);len=31;theta=1l;psf=fspecial('motion',len,theta);blurred=imfilter(i,psf,'circular','conv');subplot=(3,3,2),imshow(blurred);wnr1=deconvwnr(blurred,psf);subplot(3,3,3),imshow(wnr1);noise=0.1*rand(size(blurred));j=im2uint8(noise);blurrednoise=imadd(blurred,j);subplot(3,3,4),imshow(blurrednoise)wnr2=deconvwnr(blurred,psf);subplot(3,3,5),imshow(wnr2);nsr二sum(noise(:)「2/sum(im2double(i(:))."2);%(blurred)(:))wnr3=deconvwnr(blurrednoise,psf,nsr);subplot(3,3,6),imshow(wnr3)wnr31=deconvwnr(blurrednoise,psf,nsr/2);subplot(3,3,7),imshow(wnr3)np=abs(fftn(noise))."2;ncorr=fftshift(real(ifftn(np)));ip=abs(fftn(im2double(blurred)))「2;icorr=fftshift(real(ifftn(ip)));wnr4=deconvwnr(blurrednoise,psf,ncorr,icorr);subplot(3,3,8),imshow(wnr4);實驗五、4、求出圖像的七階矩clcI=imread('Boat.bmp');I1=I;imshow(I1);I2=imrotate(I,5,'bilinear');figure,imshow(I2);I3=fliplr(I);figure,imshow(I3)I4=imresize(I,0.5,'bilinear');figure,imshow(I4)display(源圖像')qijieju(I1);display('旋轉(zhuǎn)變化')qijieju(I2);display('鏡像變化')qijieju(I3);display('尺度變化')qijieju(I4);functionqijieju(I0)A=double(I0);[nc,nr]=size(A);[x,y]=meshgrid(1:nr,1:nc);x=x(:);y=y(:);A=A(:);m00=sum(A);ifm00==0m00=eps;endm10=sum(x.*A);m01=sum(y.*A);xmean=m10/m00;ymean=m01/m00;cm00=m00;cm02=sum((y-ymean).A2.*A))/(m00A2);cm03=sum((y-ymean).A3.*A))/(m00A2.5);cm11=sum((x-xmean).*(y-ymean).*A))/(mO0A2);cm12=sum((x-xmean).*(y-ymean).A2.*A))/(m00A2.5);cm20=sum((x-xmean).A2.*A))/(m00A2);cm21=sum((x-xmean).A2.*(y-ymean).*A))/(m00A2.5);cm30=sum((x-xmean).A3.*A))/(m00A2.5);ju(1)=cm20+cm02;ju(2)=(cm20-cm02)A2+4*cm11A2;ju(3)=(cm30-3*cm12)A2+(3*cm21-cm03)A2;ju(4)=(cm30+cm12)A2+(cm21+cm03)A2;ju(5)=(cm30-3*cm12)*(cm30+cm12)*(((cm30+cm12)A2-3*(cm21+cm03)A2)+(3*cm21-cm03)*(cm21+cm03)*