Matlab：修改3D单元矩阵的每个立方体

Question

我使用mat2cell创建了一个12 x 13 x 5的单元格矩阵，其中每个单元格包含一个29 x 29 x 29的立方体。 如何生成仅包含每个29 x 29 x 29立方体的中心21 x 21 x 21立方体的新3D单元矩阵？ 也就是说，因此新矩阵仍然是12 x 13 x 5，但每个矩阵都包含21 x 21 x 21原始立方体的中央体素？

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,:);