基于水平集的SAR图像分割 源代码

function [phi,mean_1,mean_2]=EVOLUTION(I,phi0,nu,step,numIter) [m,n]=size(phi0); p=zeros(size(phi0)); for k=1:numIter for i=1:m for j=1:n if phi0(i,j)>0 p(i,j)=1; else p(i,j)=0; end end end area_1=sum(sum(p)); I0=double(I).*p; Sarea_1=sum(sum(I0)); mean_1=Sarea_1/(area_1+eps); ER_1=log(mean_1+eps)+I/(mean_1+eps); area_2=sum(sum(~p)); I0=double(I).*(~p); Sarea_2=sum(sum(I0)); mean_2=Sarea_2/(area_2+eps); ER_2=log(mean_2+eps)+I/(mean_2+eps); di=gradient(double(phi0)); sum_phi0=sum(abs(di)); gm=max(sum_phi0); phi0=BoundMirrorExpand(phi0); curv=curvature_1(phi0); phi0=BoundMirrorShrink(phi0); curv=BoundMirrorShrink(curv); phi0=-1*step*(double(ER_1)-double(ER_2)+1*nu*curv)*gm+phi0; end for i=1:size(phi0,1) for j=1:size(phi0,2) if phi0(i,j)>0 phi0(i,j)=1; else phi0(i,j)=0; end end end phi0=bwareaopen(phi0,10); phi0=~phi0; phi0=bwareaopen(phi0,10); phi=double(phi0);