如何实现snake算法的matlab编程,来提取图像中划痕缺陷的边缘，谁有程序不妨发上来

alpha = alpha* ones(1,N)

beta = beta*ones(1,N)

% produce the five diagnal vectors

alpham1 = [alpha(2:N) alpha(1)]

alphap1 = [alpha(N) alpha(1:N-1)]

betam1 = [beta(2:N) beta(1)]

betap1 = [beta(N) beta(1:N-1)]

a = betam1

b = -alpha - 2*beta - 2*betam1

c = alpha + alphap1 +betam1 + 4*beta + betap1

d = -alphap1 - 2*beta - 2*betap1

e = betap1

% generate the parameters matrix

A = diag(a(1:N-2),-2) + diag(a(N-1:N),N-2)

A = A + diag(b(1:N-1),-1) + diag(b(N), N-1)

A = A + diag(c)

A = A + diag(d(1:N-1),1) + diag(d(N),-(N-1))

A = A + diag(e(1:N-2),2) + diag(e(N-1:N),-(N-2))

invAI = inv(A + gamma * diag(ones(1,N)))

for count = 1:ITER,

vfx = interp2(fx,x,y,'*linear')

vfy = interp2(fy,x,y,'*linear')

% deform snake

x = invAI * (gamma* x + kappa*vfx)

y = invAI * (gamma* y + kappa*vfy)

end

利用MATLAB自带的regionprops连通区函数。

记上图为I。

[C num]=bwlabel(I)

s = regionprops(C,'MinorAxisLength','MajorAxisLength')

'MajorAxisLength'：是标量，与区域具有相同标准二阶中心矩的椭圆的长轴长度（像素意义下）。

'MinorAxisLength'：是标量，与区域具有相同标准二隐仿阶中心矩的椭圆的短轴长度（像素意义下）。灶戚纤

num为裂隙的数量，MajorAxisLength为长度仔轿，MinorAxisLength为宽度

附MATLAB连通区函数详解

网页链接