get_dt_vector_kappa.m

来自「图像分割中常用的水平集方法的matlab源代码」· M 代码 · 共 22 行

function [dt] = get_dt_vector_kappa(alpha,dx,dy,u,v,b, dx2, dy2)
%
% Calculate the Euler time step.
%
% Author: Baris Sumengen  [e-mail]
% http://vision.ece.ucsb.edu/~sumengen/
%



if alpha <= 0 | alpha >= 1 
    error('alpha needs to be between 0 and 1!');
end

maxs = max(max(abs(u)/dx + abs(v)/dy + (2*b)/dx2 + (2*b)/dy2));
dt = alpha/(maxs+(maxs==0));