Matlab:修改每个三维单元矩阵的立方体

Question

我使用mat2cell创建了一个12x13x5单元矩阵，其中每个单元格包含一个29x29x29立方体。如何生成一个新的三维单元格矩阵，其中只包含每个29x29x29立方体的中心21x21x21立方体？也就是说，新的矩阵仍然是12x13x5，但是每个矩阵都包含一个2x21x21的原始立方体中心体素？

Answer 1

要裁剪一个独立块，可以使用一个简单的索引表达式：

c=4 % cut of c rows on each side.
for ix=1:numel(X)
   X{ix}=X{ix}(1+c:end-c,1+c:end-c,1+c:end-c);
end

如果您想要一个性能更好的解决方案，我建议切换到多维矩阵：

c=4
Y=reshape(X,29,12,29,13,29,5);
%now the first block is squeeze(Y(:,1,:,1,:,1))
%now cut to smaller blocks:
Y=Y(1+c:end-c,:,1+c:end-c,:,1+c:end-c,:);