如何使用Scipy从3D体积中提取任意2D切片？

Question

我正在使用Scipy从3D数据渲染平面（矢量200x200x200）。 我可以通过2个向量或向量和角度指定所需平面。 我想从这个3D卷中提取这样一个任意切片。 我在Matlab中找到了如何做到这一点： http://www.mathworks.com/help/techdoc/ref/slice.html 我如何在Scipy中做到这一点？

Answer 1

您可以使用scipy.ndimage.interpolation.rotate将您的3D数组旋转到您想要的任何角度（它使用样条插值），然后您可以从中获取切片。

Answer 2

def extract_slice（data，triplet）：
    “”“
    算法：
    1.找到平面与数据框边缘的交点     2.对于这些pts，找到面向轴的b-box
    3.找到从R2到R3的“后”反转（A，T），如X'= AX + T
       使用容易计算的（0,0），（0，h），（w，0）
    4.对2D（w，h）图像中的每个值使用trans（使用三线性插值）     “”“

我会在几个月内正确发布代码，我相信这是该项目的一部分。