showbianyuan2.m

来自「水平集分割算法的C-V模型」· M 代码 · 共 13 行

function showbianyuan2(phi1,phi2,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)
                    u1(i,j)=0;
          else
               u1(i,j)=255;
        end
     end
end
subplot(2,2,4);
imshow(u1,[]);
title('边缘');