医学 医学图像的三维可视化学习课件.pptVIP

精品文档*三维表面重建重建数据的采集边界轮廓曲线表面绘制设置图像的颜色及阴影效果设置图像光照效果设置图像的显示效果精品文档*重建数据的采集运用傅立叶级数的系数，求出边界上若干个点x，y向坐标值，并为其加上适当的z坐标值xo=[0:pi/180:2*pi];%x的值在[0,2π]中选取yo=yo+a(i)*cos((i-1)*xo)+b(i)*sin((i-1)*xo);%通过傅立叶系数求y值，其中yo初始值为a0consx=[consx;yo.*cos(xo)];%将x，y值从极坐标系转换到直角坐标系consy=[consy;yo.*sin(xo)];consz=[consz;ones(1,length(xo))*iLayer*(-4.0)];%为每一切片层赋予z坐标值，iLayer为层数精品文档*边界轮廓曲线表面绘制surf(consx,consy,consz);%利用surf()函数进行三维表面绘制。精品文档*设置图像的颜色及阴影效果。colormap(gray);%利用colormap()函数为图像定义颜色集shadingflat;%利用shading定义显示图像的颜色阴影精品文档*设置图像光照效果light(Position,[-80,-262,-200],style,infinite);%利用light()函数为图像设置光照效果light(Position,[-500,-0,-4500],style,infinite);light(Position,[5000,100,-300],style,infinite);精品文档*设置图像的显示效果

view(-144,20);%利用view()函数定义观察者视角lightinggouraud;%利用lighting定义显示图像的光线阴影axisequal;%利用axis定义显示图像的轴精品文档*精品文档*运用MATLAB程序在进行CT图像边界轮廓提取的基础上得到三维表面重建图像。重建速度快、效果好；但是面绘制的缺点是信息的丢失比较大，运算量与景物和物体形状有关。精品文档*CT图像三维体重建体绘制通过计算所有体素对光线的作用得到二维投影图像，基于体绘制的三维体重建方法计算量不依赖于景物的复杂程度和物体形状的复杂程度，也不需要对切片的边界轮廓进行提取，其计算过程不依赖于视点，处理三维采样信号方便，便于显示物体的内部结构。但是，三维体重建所需数据量大，运算速度较慢。精品文档*重建数据的采集重建数据预处理计算数据集在显示平面累计投影构造三维体重建碎片设置图像的颜色、阴影及显示效果精品文档*重建数据的采集对现有的n幅头部CT图像数据进行三维数据集D的构造，得到的数据集D为一个x×y×n的矩阵image1=imread(01.bmp);%使用imread()函数读入现有的n幅图像image2=imread(02.bmp);﹕﹕imagen=imread(n.bmp);D＝cat(3,image1,image2,image3,……imagen);%使用cat()函数创建三维矩阵D精品文档*重建数据预处理采用上述方法构造的三维数据集D，数据量大，在体重建中速度慢，并且可能在计算中超出内存。因而，可以根据实际情况，对数据集D进行预处理，减少数据量。[xyzD]=reducevolume(D,[abc]);%使用reducevolume()函数减少数据量，其中a，b，c为x，y和z轴数据抽取的比例，根据数据情况自行定义。D=smooth3(D);%使用smooth()函数对数据进行平滑处理精品文档*计算数据集在显示平面累计投影fv=isosurface(x,y,z,D,isovalue);%使用isosurface()函数计算数据集在显示平面累计投影，isovalue根据实际情况自行定义精品文档*构造三维体重建碎片p=patch(fv,FaceColor,yellow,EdgeColor,none);%使用patch()函数对碎片进行构造，并对图像的颜色，光线进行定义，其中fv是第(3)步中得到的。精品文档*设置图像的颜色、阴影及显示效果colormap(gray);%利用colormap()函数为图像定义颜色集view(3);%利用view()函数定义观察者视角lightinggouraud;%利用lighting定义显示图像的光线阴影axisequal;%利用axis定义显示图像的轴daspect([xyz]);%使用daspect()定义x、y、z轴的显示比例精品