python三维模型表面积计算方法

# -*- coding: utf-8 -*- ''' CSDN博客：https://blog.csdn.net/suiyingy 乐乐感知学堂公众号 长期分享三维、二维感知、大模型等AI算法和demo以及金融量化分析算法等。 本程序作用：open3d三维模型表面积计算 专栏地址：https://blog.csdn.net/suiyingy/category_12462636.html ''' import numpy as np import open3d as o3d if __name__ == '__main__': vertices = np.array([[0.0, 0.0, 1.0], [3.0, 0.0, 1.0], [0.0, 4.0, 1.0]]) faces = np.array([[1, 2, 3]]) faces = faces - 1 # 减1是因为索引从0开始 mesh = o3d.geometry.TriangleMesh() mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 1: ', mesh.get_surface_area()) vertices = np.array([[0.0, 0.0, 0.0], [3.0, 0.0, 0.0], [0.0, 4.0, 0.0], [0.0, 0.0, 5.0], [3.0, 0.0, 5.0], [0.0, 4.0, 5.0]]) faces = np.array([[1, 3, 2], [2, 3, 5], [5, 3, 6], [1, 4, 3], [3, 4, 6], [1, 2, 4], [2, 5, 4], [4, 5, 6]]) faces = faces - 1 # 减1是因为索引从0开始 # 创建TriangleMesh对象 mesh = o3d.geometry.TriangleMesh() # 设置顶点和面片 mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 2: ', mesh.get_surface_area()) # 改变两个底面顶点顺序 faces = np.array([[1, 2, 3], [2, 3, 5], [5, 3, 6], [1, 4, 3], [3, 4, 6], [1, 2, 4], [2, 5, 4], [4, 6, 5]]) faces = faces - 1 # 减1是因为索引从0开始 mesh = o3d.geometry.TriangleMesh() mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 3: ', mesh.get_surface_area()) # 改变两个底面和一个侧面顶点顺序 faces = np.array([[1, 2, 3], [3, 2, 5], [5, 3, 6], [1, 4, 3], [3, 4, 6], [1, 2, 4], [2, 5, 4], [4, 6, 5]]) faces = faces - 1 # 减1是因为索引从0开始 mesh = o3d.geometry.TriangleMesh() mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 4: ', mesh.get_surface_area()) # 删掉一个底面 faces = np.array([[2, 3, 5], [5, 3, 6], [1, 4, 3], [3, 4, 6], [1, 2, 4], [2, 5, 4], [4, 5, 6]]) faces = faces - 1 # 减1是因为索引从0开始 mesh = o3d.geometry.TriangleMesh() mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 5: ', mesh.get_surface_area()) # 再删掉一个底面 faces = np.array([[2, 3, 5], [5, 3, 6], [1, 4, 3], [3, 4, 6], [1, 2, 4], [2, 5, 4]]) faces = faces - 1 # 减1是因为索引从0开始 mesh = o3d.geometry.TriangleMesh() mesh.vertices = o3d.utility.Vector3dVector(vertices) mesh.triangles = o3d.utility.Vector3iVector(faces) print('area result 6: ', mesh.get_surface_area())