CT图像三维体绘制

MATLAB实现三维重建，对CT图像实现三维体绘制。

 %1.重建数据的采集
 image1=imread（‘1.jpg‘）;%使用imread（）函数读入现有的n幅图像
 image2=imread（‘2.jpg‘）;
 image3=imread（‘3.jpg‘）;
 image4=imread（‘4.jpg‘）;
 image5=imread（‘5.jpg‘）;
 image6=imread（‘6.jpg‘）;
 image7=imread（‘7.jpg‘）;
 image8=imread（‘8.jpg‘）;
 image9=imread（‘9.jpg‘）;
 image10=imread（‘10.jpg‘）;
 D=cat（3image1image2image3image4image5image6image7image8image9image10）;%使用cat（）函数创建三维矩阵D
 %2.重建数据项预处理
 
 [x y z D]=reducevolume（D[5 5 1]）;%使用reducevolume（）函数减少数据量，其中a，b，c为
                           %x，y，z轴数据抽取的比例，根据数据情况自行定义
 figure                         
 D=smooth3（D）;%使用smooth（）函数对数据惊醒平滑处理
 %3.计算数据集在显示平面累计投影
 
 fv=isosurface（xyzD‘noshare‘）;%使用isosurface（）函数计算数据集在显示平面累计投影，（从体数据中提取等值面数据.
                                %isovalue根据实际情况自行定义
 %4.构造三维体重建碎片
 p=patch（fv‘FaceColor‘[1.75.65]‘EdgeColor‘‘none‘）;%使用patch（）函数对碎片进行构造
                                            %并对图像的颜色，光线进行定义，其中fv是
                                           %第三步中得到的
 %5.设置图像的颜色，阴影及显示效果
 colormap（gray）;%利用colormap（）函数为图像定义颜色集
 view（3）;%利用view（）函数定义观察者视角
 lighting gouraud;%利用lighting定义显示图像的光线阴影
 axis equal;%利用axis定义显示图像的轴
 daspect（[11.08]）;%使用daspect（）定义x，y，z轴的显示比例

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-01-08 11:35  实验\
     目录           0  2013-01-08 11:35  实验\1\
     文件       20383  2012-07-13 19:13  实验\1\1.jpg
     文件       10340  2012-07-13 19:14  实验\1\10.jpg
     文件       20938  2012-07-13 19:14  实验\1\2.jpg
     文件       19431  2012-07-13 19:13  实验\1\3.jpg
     文件       19241  2012-07-13 19:14  实验\1\4.jpg
     文件       19091  2012-07-13 19:13  实验\1\5.jpg
     文件       18306  2012-07-13 19:13  实验\1\6.jpg
     文件       16884  2012-07-13 19:14  实验\1\7.jpg
     文件       15362  2012-07-13 19:13  实验\1\8.jpg
     文件       13203  2012-07-13 19:14  实验\1\9.jpg
     文件        1456  2012-07-19 17:10  实验\1\CT.m
     文件        2482  2012-07-19 17:57  实验\1\CT1.asv
     文件        2493  2012-07-19 17:58  实验\1\CT1.m
     文件        1371  2012-07-19 12:45  实验\1\CT程序说明.txt
     文件       24465  2012-07-21 11:44  实验\1\untitled.jpg
     文件       19055  2012-07-21 11:44  实验\1\untitled1.jpg
     目录           0  2013-01-08 11:35  实验\2\
     文件        2572  2012-07-24 18:10  实验\2\CT1.m