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

问题描述

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

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

要转换为体素坐标,数组的代码和结果形状如下:

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

理想情况下,我希望生成的数组具有以下形状:(160,4) 而不是 (160,2086604)

你们中的任何人都可以帮我解决这个问题吗?

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)