Code：

本代码完成的功能是：

对3维数据进行thresholding，之后以此作为水平集的初始化轮廓，最后在3D中显示Lsf的0水平集

clear all;

close all;

load lsf_half.mat;

im = lsf;

BW = (lsf>100);

lsf = 4.*(0.5-BW);

% using mr li

[m,n,p] = size(lsf);

[xx,yy,zz] = meshgrid(1:n,1:m,1:p);

v = lsf;

a = 0;

aspect_ratio = [1,1,1];