image - Matlab DICOM 切片

Question

我有一个 DICOM 图像作为矩阵加载到 matlab 中。我的问题是，如何在每个正交方向上显示该图像的特定切片？

比如在位置 100 处查看切片 x，y=0，z=0

Answer 1

如果您的矩阵是 M，并且具有 d 维(3，或者您有什么)并且您想要绘制其中一个维的一维“切片”，则:

  plot(squeeze(M(n1,n2, ...,:,...));

其中 n1,n2,... 是要切片的维度 x,y,... 的位置，运算符 (:) 是您要绘制的维度。

例如，给定一个 5 维矩阵 M=rand(10,10,10,10,10)，让我们围绕某些点(x=n1、y=n2 等...)对第 4 维进行切片

   M=rand(10,10,10,10,10);
   n1=4; n2=7; n3=3; n5=5; 
   plot(squeeze(M(n1,n2, n3, :, n5)));

如果切片是二维的，那么您可以使用 imshow 或 imagesc 来显示二维切片，例如显示第 2 维和第 4 维:

 imagesc(squeeze(M(n1,:,n3,:,n5)))