步骤 1:加载和查看水平 MRI
此示例使用 MATLAB® 附带的用在 montage 和 immovie 的帮助示例中的 MRI 数据集。加载 mri.mat 会向工作区添加两个变量:D(128×128×1×27,uint8 类)和灰度颜色图 map(89×3,double 类)。
D 包含人体头盖骨 MRI 扫描数据的 27 个 128×128 水平切片。D 中的值的范围是从 0 到 88,因此需要颜色图来生成实用的可视图形。D 的维度与 montage 兼容。前两个维度是空间维度。第三个维度是颜色维度,大小为 1,因为它是颜色图的索引。对于 RGB 图像序列,size(D,3) 为 3。第四个维度是时序维度(和任何图像序列一样),但在本例中,它也是空间维度。因此,D 中有三个空间维度,我们可以使用 imtransform 或 tformarray 将水平切片转换为矢状切片(显示头部侧面的视图)或冠状(正面)切片(显示头部正面或背面的视图)。
D 的空间维度排序如下:
维度 1:头部从前到后(从吻突/前侧到尾椎/后侧)
维度 2:头部从左到右
维度 4:头部从底到顶(从下到上)。
一个重要因素是三个维度上的采样间隔不同:沿垂直维度 (4) 的采样间隔比沿水平维度的采样间隔宽 2.5 倍。
加载 MRI 数据集,将 27 个水平切片以蒙太奇形式进行查看。
load mri;
montage(D,map)
title('Horizontal Slices');