《小论文--三维MRI切片数据提取》毕业学术论文.docVIP

2012-2013年第一学期《数字图像处理》科目考查卷 专业：通信工程 班 级：B090507 任课教师： 姓名： 学 号： 成 绩： 三维MRI切片数据提取 一、基本概况描述 这里imtransform和tformarray函数可以用来插入和reslice三维MRI数据集,提供一种方便的方式来查看的数据量。 该演示程序使用MRI的数据集，与MATLAB中使用的M-文件的帮助都蒙太奇和immovie的例子。装载MRI添加到工作区中的两个变量：D（128-128-1-27级UINT8）和一个灰度色彩对照表，地图（89-3，双级）。 D包括27 128-128水平切片从MRI数据扫描的人类头骨。 D系列中的值从0到88，所以颜色映射需要生成一个数字一个有用的可视范围。 D的维数，兼容immovie。前两个维度空间。第三个维度是颜色尺寸，大小为1到彩色地图，因为它的索引。 （尺寸（D，3）将3的RGB电影。）第四维是时间（与任何电影），但在这个特殊的情况下，它也是空间。因此，有三个空间维度在D和转换的矢状切片（从侧面示出了视图头部）或冠状面（正面）切片（示出了视图从正面或背面的水平切片，我们可以使用imtransform或tformarray头）。空间尺寸的D命令如下：尺寸：从前到后的头（头端/前尾侧/后）尺寸：左到右的头尺寸：底部到头顶（劣至优）。 二、窗口化图像实现与操作说明 1.该演示包括以下步骤： 第1步：加载和显示水平MRI切片影像 D空间维度的命令如下： 维度1：从头部的前到后（头端/前尾侧/后）； 维度2：从头部的左到右； 维度4：从头底到头顶（劣至优）。 沿着三个维度采样间隔不一样是一个重要因素：沿垂直维度（4）的样品的间距的2.5倍，更广泛地比沿水平尺寸。加载MRI数据集，可以看到水平方向的27张影像，然后作为一个蒙太奇混合处理。做蒙太奇混合处理时，为了避免TRUESIZE警告，该示例使用iptgetpref关闭警告。 truesizewarning = iptgetpref(TruesizeWarning); iptsetpref(TruesizeWarning,off); % Turn warning off load mri; figure; immovie(D,map); montage(D,map); title(Horizontal Slices); 第2步：使用Imtransform从水平切片提取矢状面片 我们可以通过D的一个子集构造一个MRI数据并将其转化为下面占不同的采样间隔时间和空间定位尺寸的语句.提取所需的所有数据的正中矢状切片中矢状面片。 下面的语句提取物的矢状切片所需的所有数据。 M1 = D(:,64,:,:); size(M1) 然而，我们不能查看M1作为图像，因为它是128-1-1-27。重塑（或紧缩），可以将M1到一个128-27图像，可视与imshow。 M2 = reshape(M1,[128 27]); size(M2) figure, imshow(M2,map); title(Sagittal - Raw Data); 在M2的尺寸进行排序如下：维度：从前到后的头（吻端到尾鳍）；维度：头底部到头顶（劣至优）。通过变换M2以改变其方向，并增加了2.5倍的沿垂直方向的采样（劣质高级）维度 - 使得采样间隔等于在所有三个空间维度，我们可以得到一个令人满意得多的视图。为此，我们可以在转置开始的步骤，但以下的仿射变换使一个单步的变换和更经济地利用存储器。 上面的2×2的矩阵块传递给maketform[0 -2.5 1 0]结合旋转和缩放。改转换后，我们得到了：维度1：头顶部至底部（优到劣）。维度2：从头部的前到后（吻端到尾鳍）的调用imtransform（M2，T0，“cube”），就能够了申请T到M2，并提供良好的分辨率，而内插沿上下方向的顶部。然而，也没有必要为三次插值中的前后方向，由于没有重采样会发生沿（输出）尺寸2。因此，我们在这个维度指定最近邻重采样，以取得更高的效率和相同的结果。 R2 = makeresampler({cubic,nearest},fill); M3 = imtransform(M2,T0,R2); figure, imshow(M3,map); title(Sagittal - IMTRANSFORM) 第3步：使用Tformarray从水平切片提取矢状面片 在此步骤中，我们获得像第2步中同样的结果，但使用tformarray单个操作中使得空间三个维度变为两个。步骤2从具有三个空间维度一个数组开始并以具有两个空间维度一个阵列结尾。但其中间的二维图像（M1和M