我是 3D 图像处理的新手。我想知道如何在python中查看dicom系列。我尝试使用 matplotlib 和 VTK。在 matplot 中,我无法像使用 volViewer 在 matlab 中查看那样查看体积。关于 VTK,我无法导入 VTKRAyCASt 来查看 3D。我使用的版本是 8.2.0。
我正在使用 scipy.ndimages 进行处理
请为我的音量 dicom 文件推荐一些资源
最佳答案
您可以尝试使用 ipyvolume https://github.com/maartenbreddels/ipyvolume 进行交互式绘图,我发现它非常有用。
此外,您可以使用 matplotlib 绘制它们,方法是使用行进立方体来获取表面网格,但速度很慢:
from mpl_toolkits.mplot3d.art3d import Poly3DCollection
import numpy as np
from skimage import measure
def plot_3d(image, threshold=-300):
p = image.transpose(2,1,0)
verts, faces, normals, values = measure.marching_cubes_lewiner(p, threshold)
fig = plt.figure(figsize=(10, 10))
ax = fig.add_subplot(111, projection='3d')
mesh = Poly3DCollection(verts[faces], alpha=0.1)
face_color = [0.5, 0.5, 1]
mesh.set_facecolor(face_color)
ax.add_collection3d(mesh)
ax.set_xlim(0, p.shape[0])
ax.set_ylim(0, p.shape[1])
ax.set_zlim(0, p.shape[2])
plt.show()
-300 HU 的阈值适用于可视化胸部 CT 扫描,但如果您打算使用 MRI(检查您的强度值分布)或二进制体积(阈值 = 0),请更改它。
有一些可视化示例:
关于python-3.x - Python 中的 3D Dicom 可视化,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/56035562/