【正文】 end kk=1:1:G。 gbest1。 figure(1)。 imshow(J)。 figure(2)。 plot(kk,GG) tt=etime(clock,t1)。endI=imread(39。39。)。tmin=min(I(:))。tmax=max(I(:))。th=(tmin+tmax)/2。ok=true。while ok g1=I=th。 g2=Ith。 u1=mean(I(g1))。 u2=mean(I(g2))。 tnew=(u1+u2)/2。 if abs(thtnew)1 ok=0。 endendth=tnew。th=floor(th)。Inew=im2bw(I,th/255)。subplot(1,2,1)imshow(I)。title(39。原始圖像39。)。subplot(1,2,2)imshow(Inew)。t=[39。迭代法分割后的圖像，閾值=39。 num2str(th)]。title(t)。A0=imread(39。39。)。seed=[100,220]。thresh=15。%相似性選擇閾值A(chǔ)=rgb2gray(A0)。A=imadjust(A,[min(min(double(A)))/255,max(max(double(A)))/255],[])。A=double(A)。 B=A。[r,c]=size(B)。n=r*c。pixel_seed=A(seed(1),seed(2))。q=[seed(1) seed(2)]。top=1。M=zeros(r,c)。M(seed(1),seed(2))=1。count=1。 while top~=0。 r1=q(1,1)。 c1=q(1,2)。 p=A(r1,c1)。 dge=0。 for i=1:1 for j=1:1 if r1+i=r amp。 r1+i0 amp。 c1+j=c amp。 c1+j0 if abs(A(r1+i,c1+j)p)=thresh amp。 M(r1+i,c1+j)~=1 top=top+1。 q(top,:)=[r1+i c1+j]。 M(r1+i,c1+j)=1。 count=count+1。 B(r1+i,c1+j)=1。 end if M(r1+i,c1+j)==0。 dge=1。 end else dge=1。 end end end if dge~=1 B(r1,c1)=A(seed(1),seed(2))。 end if count=n top=1。 end q=q(2:top,:)。 top=top1。endsubplot(1,2,1),imshow(A,[])。subplot(1,2,2),imshow(B,[])。A0=imread(39。39。)。seed=[100,220]。thresh=15。%相似性選擇閾值A(chǔ)=rgb2gray(A0)。A=imadjust(A,[min(min(double(A)))/255,max(max(double(A)))/255],[])。A=double(A)。 B=A。[r,c]=size(B)。n=r*c。pixel_seed=A(seed(1),seed(2))。q=[seed(1) seed(2)]。top=1。M=zeros(r,c)。M(seed(1),seed(2))=1。count=1。 while top~=0。 r1=q(1,1)。 c1=q(1,2)。 p=A(r1,c1)。 dge=0。 for i=1:1 for j=1:1 if r1+i=r amp。 r1+i0 amp。 c1+j=c amp。 c1+j0 if abs(A(r1+i,c1+j)p)=thresh amp。 M(r1+i,c1+j)~=1 top=top+1。 q(top,:)=[r1+i c1+j]。 M(r1+i,c1+j)=1。 count=count+1。 B(r1+i,c1+j)=1。 end if M(r1+i,c1+j)==0。 dge=1。 end else dge=1。 end end end if dge~=1 B(r1,c1)=A(seed(1),seed(2))。 end if count=n top=1。 end q=q(2:top,:)。 top=top1。endsubplot(1,2,1),imshow(A,[])。subplot(1,2,2),imshow(B,[])。A0=imread(39。39。)。seed=[100,220]。thresh=15。%相似性選擇閾值A(chǔ)=rgb2gray(A0)。A=imadjust(A,[min(min(double(A)))/255,max(max(double(A)))/255],[])。A=double(A)。 B=A。[r,c]=size(B)。n=r*c。pixel_seed=A(seed(1),seed(2))。q=[seed(1) seed(2)]。top=1。M=zeros(r,c)。M(seed(1),seed(2))=1。count=1。 while top~=0。 r1=q(1,1)。 c1=q(1,2)。 p=A(r1,c1)。 dge=0。 for i=1:1 for j=1:1 if r1+i=r amp。 r1+i0 amp。 c1+j=c amp。 c1+j0 if abs(A(r1+i,c1+j)p)=thresh amp。 M(r1+i,c1+j)~=1 top=top+1。 q(top,:)=[r1+i c1+j]。 M(r1+i,c1+j)=1。 count=count+1。 B(r1+i,c1+j)=1。 end if M(r1+i,c1+j)==0。 dge=1。 end else dge=1。 end end end if dge~=1 B(r1,c1)=A(seed(1),seed(2))。 end if count=n top=1。 end q=q(2:top,:)。 top=top1。endsubplot(1,2,1),imshow(A,[])。subplot(1,2,2),imshow(B,[])。A0=imread(39。39。)。seed=[100,220]。thresh=15。%相似性選擇閾值A(chǔ)=rgb2gray(A0)。A=imadjust(A,[min(min(double(A)))/255,max(max(double(A)))/255],[])。A=double(A)。 B=A。[r,c]=size(B)。n=r*c。pixel_seed=A(seed(1),seed(2))。q=[seed(1) seed(2)]。top=1。M=zeros(r,c)。M(seed(1),seed(2))=1。count=1。 while top~=0。 r1=q(1,1)。 c1=q(1,2)。 p=A(r1,c1)。 dge=0。 for i=1:1 for j=1:1 if r1+i=r amp。 r1+i0 amp。 c1+j=c amp。 c1+j0 if abs(A(r1+i,c1+j)p)=thresh amp。 M(r1+i,c1+j)~=1 top=top+1。 q(top,:)=[r1+i c1+j]。 M(r1+i,c1+j)=1。 count=count+1。 B(r1+i,c1+j)=1。 end if M(r1+i,c1+j)==0。 dge=1。 end else dge=1。 end end end if dge~=1 B(r1,c1)=A(seed(1),seed(2))。 end if count=n top=1。 end q=q(2:top,:)。 top=top1。endsubplot(1,2,1),imshow(A,[])。subplot(1,2,2),imshow(B,[])。附錄二 labview程序框圖參考文獻36