vtk-py求3d模型表面积

模型格式：.obj

环境：python3.6+vtk7.1

vtk版：

 import vtk

 filename = "XXXX.obj"
 reader = vtk.vtkOBJReader()
 reader.SetFileName(filename)
 reader.Update()
 triangleFilter = vtk.vtkTriangleFilter()
 triangleFilter.SetInputData(reader.GetOutput())
 triangleFilter.Update()

 polygonProperties = vtk.vtkMassProperties()
 polygonProperties.SetInputData(triangleFilter.GetOutput())
 polygonProperties.Update()

 area = polygonProperties.GetSurfaceArea()
 print(area)

自己写的版本：

先将文件数据读入数组face，和vertices

 #计算3d 面片面积
 def getarea(a, b, c):
     ab = np.sqrt((a[0] - b[0]) ** 2 + (a[1] - b[1]) ** 2 + (a[2] - b[2]) ** 2)
     ac = np.sqrt((a[0] - c[0]) ** 2 + (a[1] - c[1]) ** 2 + (a[2] - c[2]) ** 2)
     bc = np.sqrt((b[0] - c[0]) ** 2 + (b[1] - c[1]) ** 2 + (b[2] - c[2]) ** 2)
     s = (ab + ac + bc) / 2.0
     area = np.sqrt(s * (s - ab) * (s - ac) * (s - bc))
     return area
 def getarea(self):   #总片面面积
     for index in self.faces:
         self.area = self.area + getarea(self.vertices[index[0]-1], self.vertices[index[1]-1], self.vertices[index[2]-1])

经测试，两个算出来的误差非常小，可以使用。

但是如果用vtk算体积的话，会有误差，可能适合模型封闭性有关。