数字图像处理从3维MRI数据集中提取切片图小论文.doc

2012-2013年第一学期《数字图像处理》科目考查卷 专业：通信工程 班 级： 任课教师： 姓名： 学 号： 成 绩： Extracting Slices from a 3-Dimensional MRI Data Set 如何从3维MRI数据集中提取切片图 从3维MRI数据集中提取切片图是一种计算机系统和方法，通过收集在一个集2D计算机体层摄影术（CT）的切片的病人的骨骼和一个同轴集2D CT球探图像，数字二维透视图像，可以是形式的三维数据的通过CT扫描仪产生的从每个这些三维数据中提取的二维轮廓设置一个相应的堆栈构造对象的3D几何模型 在matlab中，imtransform和tformarray功能可用于MRI数据集的三维内插和重新切割，提供一种方便的方式来查看的数据量 。 内容 第1步：加载和水平MRI 第2步：提取矢状切片水平切片IMTRANSFORM 第3步：提取矢状切片水平切片TFORMARRAY 第4步：创建和显示矢状切片 第5步：创建和显示冠状切片 第1步：加载和水平MRI 该演示程序使用MRI的数据集，与MATLAB?，用于为蒙太奇和immovie都在帮助例子。装载mri.mat添加到工作区中的两个变量：e（128-128-1-27级UINT8）和一个灰度色彩对照表，地图（89-3，双级）。 e包括27 128-128水平切片从MRI数据扫描的人类头骨。e范围从0到88的值，所以colormap的需要产生一个有用的视觉范围内的数字。D的维数，使得它适合用蒙太奇。前两个维度空间。第三个维度是颜色尺寸，大小为1到彩色地图，因为它的索引。（尺寸（D，3）将3的RGB图像序列。）第四维是时间（与任何图像序列），但在这个特殊的情况下，它也是空间。因此，有在 D中的三个空间维度，我们可以使用imtransform或tformarray矢状切片（从侧面示出了视图头部）或冠状面（正面）切片（示出了视图从正面或背面的水平切片转换头）的空间尺寸的 D 顺序如下： 尺寸1：从前到后的头（头端/前向尾侧/后） 维度2：左到右的头 尺寸：底部到头顶（劣至优）。 一个重要因素是，采样间隔是不一样的，沿着三个维度：沿垂直维度（4）的样品的间距的2.5倍，更广泛地比沿水平尺寸。 MRI数据集加载和查看蒙太奇的27水平切片。 加载MRI ;蒙太奇（D，地图）标题（“水平切片”） 第2步：提取矢状切片水平切片IMTRANSFORM 我们可以构造一个中期矢状切片从MRI数据采取的 D的一个子集，将其转化为不同的采样间隔和空间方向的尺寸的 D到账户 。 下面的语句提取物的矢状切片所需的所有数据。 M1 = D（64，:);尺寸（M1） ANS = 128 1 1 27 然而，我们无法查看M1作为图像，因为它是128-1-1-27。重塑（或挤压）可以转换到一个128-27的图像是可见的与imshow M1。 M2 =重塑（M1，[128 27]），尺寸（M2）图，imshow（M2，地图）;标题（矢状面-原始数据“）; ANS = 128 27 在M2中的尺寸如下排列： 尺寸1：从前到后的头（吻端到尾鳍） 尺寸：底部到头顶（劣至优）。 通过变换M2以改变其方向，并增加了2.5倍的沿垂直方向的采样（劣质高级）维度-使得采样间隔等于在所有三个空间维度，我们可以得到一个令人满意得多的视图。为此，我们可以在转置开始的步骤，但以下的仿射变换使一个单步的变换和更经济地利用存储器。 T0 = maketform（ “ 仿射，[0 -2.5 0 1 0 0]）; 上部的2×2的矩阵块传递给maketform，0 -2.5 1 0] ，结合旋转和缩放。改造后，我们有： 尺寸1：排名最前的底部头（优到劣）。 维度2：从前到后的头（吻端到尾鳍） imtransform（M2，T0，“立方”） 就足够了申请T到M2，并提供良好的分辨率，而内插沿上下方向的顶部。然而，也没有必要为三次插值中的前后方向，由于没有重采样会发生沿（输出）尺寸2。因此，我们在这个维度指定最近邻重采样，以更高的效率和相同的结果。 R2 = makeresampler（{ 立方，最近 }， “ 填补”）;M3 = imtransform（M2，T0，R2）;图，imshow（M3，地图）;标题（的矢状面- IMTRANSFORM“ ） 第3步：提取矢状切片水平切片TFORMARRAY 在此步骤中，我们获得同样的结果作为第2步中，，但使用tformarray去在单个操作中从三个空间维度两个。步骤2从一个数组具有三个空间维度和结尾的一个阵列，其具有两个空间维度，但中间的二维图像（M1和M2）的方式为的调用imtransform创建M3