如何使用Python绑定在ITK中创建网格对象？

Question

我正在尝试在Python中创建一个Mesh对象。我使用的是从以下网页安装的Python绑定。关于C++代码，我们可以这样做：

      MeshType::Pointer mesh = MeshType::New();

我对ITK还很陌生，甚至不知道怎么创建它。在C++文档中，Mesh的构造函数提到一个必要的参数，叫做TPixelType。我也找不到这个参数。

请问有人能帮我一下吗？

谢谢

Answer 1

如果我是你，我会看看ITK 4.0自带的Python绑定。你可以通过在cmake中开启WRAP_ITK_PYTHON这个选项来访问它们。

一旦你编译了ITK并开启了Python绑定，你就可以直接创建两种网格类型：

import itk
meshType2D = itk.Mesh.D2Q.New()
meshType3D = itk.Mesh.D3Q.New()

另外，你也可以像下面这样明确地实例化你的类：

import itk
meshType2D = itk.Mesh[itk.D, 2, itk.QuadEdgeMeshTraits.D2BBFF]
meshType3D = itk.Mesh[itk.D, 3, itk.QuadEdgeMeshTraits.D3BBFF]

这样你就能得到2维和3维的网格，像素值是双精度类型，使用的是默认的网格特性。至于ITK中的像素类型，这些其实就是基本的C++变量类型，比如：double、float、unsigned int等等。这些基本类型在Python中被封装了，可以在ITK的命名空间里找到：itk.D、itk.F、itk.UI等等。