showbianyuan4.m

水平集分割算法的C-V模型

function showbianyuan2(phi1,phi2,phi3,phi4,nrow,ncol)
for i=2:nrow
     for j=2:ncol
        if ((phi1(i,j))*(phi1(i-1,j))<0)|((phi1(i,j))*(phi1(i,j-1))<0| (phi2(i,j))*(phi2(i-1,j))<0)|((phi2(i,j))*(phi2(i,j-1))<0)| ((phi3(i,j))*(phi3(i-1,j))<0)|((phi3(i,j))*(phi3(i,j-1))<0| (phi4(i,j))*(phi4(i-1,j))<0)|((phi4(i,j))*(phi4(i,j-1))<0)
                    u1(i,j)=0;
          else
               u1(i,j)=255;
        end
     end
end
subplot(2,2,4);
imshow(u1,[]);
title('边缘');