Matlab DICOM切片

时间:2012-10-13 21:01:38

标签: image matlab slice dicom

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

与位置100处的视图切片x类似,y = 0,z = 0

1 个答案:

答案 0 :(得分:4)

如果您的矩阵是M,并且有d维(3,或者你有什么),并且你想绘制其中一个维度的1-D“切片”,那么:

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

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

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

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

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

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