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];