cnsegm.rar_MATLAB 复杂网络_图像复杂度_复杂网络_复杂网络 matlab_复杂网络 图像

function [Threshold]=cnsegm(F) [m,n,k]=size(F); if k==3 F=rgb2gray(F); end Fd=double(F); [R]=size_adj(F); [m,n,k]=size(R); if(k==3) U=rgb2gray(R); else U=R; end U=double(U); E=double(U); E=E'; sig=mean(mean(Fd))/127.5/((max(max(Fd))-min(min(Fd)))/255) % /255 Tn=2*variance(F)/100000*sig [m,n]=size(E); M=m*n; W=zeros(M); for i=1:M for j=i:M W(i,j)=abs(E(j)-E(i)); W(j,i)=W(i,j); end end W=W/255; W=double(W<=Tn&W~=0);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%51*2 Wp=W; for i=1:M W(i,i)=0; end nodemax=[]; for N=1:3 Maxi=find(sum(W)==max(sum(W))); t=Maxi(1); a=find(W(t,:)==1); W(a,:)=0; W(:,a)=0; xt=floor((t-1)/m)+1; yt=mod(t-1,m)+1; nodemax=[nodemax,U(xt,yt)]; x=floor((a-1)/m)+1; y=mod(a-1,m)+1; end nodemax=ceil(nodemax) nd=sum(Wp); pixd=zeros(1,255); for pix=1:255 Eind=(find(E==pix)); se=size(Eind); if se(1)~=0 pixd(pix)=nd(Eind(1)); end end fzero=find(pixd==0); flag=1; for fi=2:length(fzero)-1 if (fzero(1)==1|fzero(length(fzero))==255)&flag==1 flag=0; continue; end pixd(fzero(fi))=(pixd(fzero(fi)-1)+pixd(fzero(fi)+1))/2; end D=sort(nodemax); t1=D(1); t2=D(2); t3=D(3); figure,stem(pixd); wv=min(pixd(t1:t2));%avoid zero index; subpixd=find(pixd==wv); Threshold1=mean(subpixd(find(subpixd>=t1&subpixd<=t2))); wv=min(pixd(t2:t3)); subpixd=find(pixd==wv); Threshold2=mean(subpixd(find(subpixd>=t2&subpixd<=t3))); u=mean(mean(Fd)); u1=(t1+t3)/2; di1=Threshold1-u; di2=Threshold2-u; if abs(di1)<abs(di2) Threshold=Threshold1 else Threshold=Threshold2 end figure,im2bw(F,Threshold/255); figure,im2bw(F,graythresh(F)); end