将患者坐标转换为DICOM图像的体素坐标

Question

我有一个关于将患者坐标转换为从DICOM获得的3D MRI point cloud模型的voxel coordinates的问题。我有一个从2D MRI图像获得的3D点云文件，这些点位于患者坐标系中。我试图使用我最初创建的仿射矩阵的逆矩阵将云中的这些点转换为voxel coordinates。但是我得到的仿射矩阵是奇异的，因此我不能直接对它应用反函数。有没有一种方法可以构造逆仿射矩阵，而不在仿射矩阵上应用逆函数？或者是否有其他方法来实现此转换？我也试着从仿射矩阵中得到伪逆矩阵，但我不确定这是否有任何帮助

我在{}工作。仿射矩阵的形状是(160, 4, 4)，伪逆矩阵的形状也是(160, 4, 4)，点云阵列的形状是(2086604, 4)

要转换为体素坐标，阵列的代码和结果形状如下所示：

voxels = np.dot(p_inverse_whole, stack_patient.T)  
(160, 4, 2086604)

理想情况下，我希望得到如下形状的数组：(160, 4, 4)而不是(160, 4, 2086604)

你们谁能帮我解决这个问题