读取CSV



f=csvread('x0y0z0.csv',1,1)
t=0.5*(double(min(f(:)))+double(max(f(:))))
done=false;
while ~done
    g=f>t
    tnext=0.5*(mean(f(g))+mean(f(~g)))
    done=abs(t-tnext)<0.5
    t=tnext
end
line1=(line>=t)*255
%plot(1:5000,line,1:5000,line1,'*')
line2=line1(2:5000)
det=abs(line2(1:4999)-line1(1:4999))
plot(1:5000,line,1:4999,det,'r')