vtk-py求3d模型表面积
模型格式：.obj
环境：python3.6+vtk7.1
vtk版：
1 import vtk 2 3 filename = "XXXX.obj" 4 reader = vtk.vtkOBJReader() 5 reader.SetFileName(filename) 6 reader.Update() 7 triangleFilter = vtk.vtkTriangleFilter() 8 triangleFilter.SetInputData(reader.GetOutput()) 9 triangleFilter.Update() 10 11 polygonProperties = vtk.vtkMassProperties() 12 polygonProperties.SetInputData(triangleFilter.GetOutput()) 13 polygonProperties.Update() 14 15 area = polygonProperties.GetSurfaceArea() 16 print(area)
自己写的版本：
先将文件数据读入数组face，和vertices
1 #计算3d 面片面积 2 def getarea(a, b, c): 3 ab = np.sqrt((a[0] - b[0]) ** 2 + (a[1] - b[1]) ** 2 + (a[2] - b[2]) ** 2) 4 ac = np.sqrt((a[0] - c[0]) ** 2 + (a[1] - c[1]) ** 2 + (a[2] - c[2]) ** 2) 5 bc = np.sqrt((b[0] - c[0]) ** 2 + (b[1] - c[1]) ** 2 + (b[2] - c[2]) ** 2) 6 s = (ab + ac + bc) / 2.0 7 area = np.sqrt(s * (s - ab) * (s - ac) * (s - bc)) 8 return area 9 def getarea(self): #总片面面积 10 for index in self.faces: 11 self.area = self.area + getarea(self.vertices[index[0]-1], self.vertices[index[1]-1], self.vertices[index[2]-1])
经测试，两个算出来的误差非常小，可以使用。
但是如果用vtk算体积的话，会有误差，可能适合模型封闭性有关。
