CT头颅三维重建

在MATLAB R2016a版本下运行，里面包含头颅的CT图像，运行结果可以显示出头颅的立体图像以及三视图。

function   [xcycRa] = circfit（xy）
%
%   [xc yx R] = circfit（xy）
%
%   fits a circle  in xy plane in a more accurate
%   （less prone to ill condition ）
%  procedure than circfit2 but using more memory
%  xy are column vector where （x（i）y（i）） is a measured point
%
%  result is center point （ycxc） and radius R
%  an optional output is the vector of coeficient a
% describing the circle‘s equation
%
%   x^2+y^2+a（1）*x+a（2）*y+a（3）=0
%
%  By:  Izhak bucher 25/oct /1991 
    x=x（:）; y=y（:）;
   a=[x y ones（size（x））]\[-（x.^2+y.^2）];
   xc = -.5*a（1）;
   yc = -.5*a（2）;
   R  =  sqrt（（a（1）^2+a（2）^2）/4-a（3））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-05-28 17:59  SkullReconstruct-master\
     文件         101  2015-05-28 17:59  SkullReconstruct-master\Contents.m
     文件         326  2015-05-28 17:59  SkullReconstruct-master\README.md
     文件         633  2015-05-28 17:59  SkullReconstruct-master\circfit.m
     文件        2879  2015-05-28 17:59  SkullReconstruct-master\ellipsefit.m
     文件         450  2015-05-28 17:59  SkullReconstruct-master\extractContour.m
     文件         304  2015-05-28 17:59  SkullReconstruct-master\plotcircle.m
     文件        2499  2015-05-28 17:59  SkullReconstruct-master\plotellipse.m
     文件        2995  2015-05-28 17:59  SkullReconstruct-master\solution.m
     文件        6371  2015-05-28 17:59  SkullReconstruct-master\solution2.m
     文件        8237  2015-05-28 17:59  SkullReconstruct-master\solution3.m
     目录           0  2015-05-28 17:59  SkullReconstruct-master\tiff\
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00001.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00002.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00003.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00004.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00005.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00006.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00007.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00008.tif
     文件      535453  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00009.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00010.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00011.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00012.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00013.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00014.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00015.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00016.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00017.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00018.tif
     文件      535461  2015-05-28 17:59  SkullReconstruct-master\tiff\Image00019.tif
............此处省略234个文件信息