基于matlab的用各向异性扩散的尺度空间和边缘检测算法仿真

function [x,y] = snakeinit(II,delta) BW = edge(II,'canny',0.6,0.3); hold on; x = []; y = []; z =0; [n,m] = size(BW); for q = 5:1:n-5 for p = 5:1:m-5 if(BW(q,p) == 1) z = z+1; x(z,1)=p; y(z,1)=q; end end end plot(x,y,'b.'); xx = []; yy = []; number1 = 0; for i = 1:z if(x(i,1) == min(x(:))) number1 = number1 + 1; xx(number1,1) = min(x(:))-3; yy(number1,1) = y(i,1); end end for i = 1:z if(y(i,1) == max(y(:))) number1 = number1 + 1; yy(number1,1) = max(y(:))+3; xx(number1,1) = x(i,1); end end number2 = 0; for i = 1:z if(x(i,1) == max(x(:))) number2 = number2 + 1; number1 = number1 + 1; end end %number1 = number1+1; for i = 1:z if(x(i,1) == max(x(:))) number1 = number1 - 1; yy(number1,1) = y(i,1); xx(number1,1) = max(x(:))+3; end end number1 = number1+ number2; number2 = 0; for i = 1:z if(y(i,1) == min(y(:))) number2 = number2 + 1; number1 = number1 + 1; end end %number1 = number1 + 1; for i = 1:z if(y(i,1) == min(y(:))) number1 = number1 - 1; xx(number1,1) = x(i,1); yy(number1,1) = min(y(:))-3; end end x = [xx;xx(1,1)]; y = [yy;yy(1,1)]; plot(x, y, '-'); hold off % sampling and record number to N z=size(x); x = [x;x(1,1)]; y = [y;y(1,1)]; t = 1:z+1; ts = [1:delta:z+1]'; xi = interp1(t,x,ts); yi = interp1(t,y,ts); n = length(xi); x = xi(1:z-1); y = yi(1:z-1);