是的,您可以像处理任何数组一样操作图像/序列。
这是一个简单的例子,其中包含来自 NIH 的数据here.
该数据集为 4D 形式,名为“filtered_func_data.nii”。
让我们加载数据集并访问img
结果结构的字段:
S = load_nii('filtered_func_data.nii')
这里,S是一个具有以下字段的结构体:
S =
hdr: [1x1 struct]
filetype: 2
fileprefix: 'filtered_func_data'
machine: 'ieee-be'
img: [4-D int16]
original: [1x1 struct]
我们可以使用该字段访问图像数据img
(你已经知道了):
A = S.img
如果我们检查尺寸,我们会得到:
size(A)
ans =
64 64 21 180
因此,数据集由 64x64 图像组成,深度/切片数量为 21,帧数量为 180。
此时我们可以操纵A
因为我们喜欢重塑形状、调整大小或任何东西。
下面是一个简单的代码(gif 动画),用于循环遍历 4D 数组中第一个时间点的每个切片:
NumSlices = size(A,3)
figure(1)
filename = 'MRI_GIF.gif';
for k = 1:NumSlices
imshow(A(:,:,k,1),[])
drawnow
frame = getframe(1);
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
if k == 1;
imwrite(imind,cm,filename,'gif', 'Loopcount',inf);
else
imwrite(imind,cm,filename,'gif','WriteMode','append');
end
pause(.1)
end
Output:
这对我来说看起来很不错。
所以在你的情况下,你得到的大小是[39 305 305]
您可以应用与我相同的操作来处理您的 3D 数据集。
EDIT对于你的数据来说是一样的:
S = load_nii('Subject01.nii');
A = S.img;
NumSlices = size(A,3);
然后,如果您想要 2D 图像,则需要在 3D 数组中选择一个切片。
例如,第一个切片的访问方式如下:
A(:,:,1)
其余的依此类推。
要将图像保存为 png,请使用imwrite.
希望有帮助!