open3d体积计算 体素计算

# -*- 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, 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) # 计算三角网格的体积 volume = mesh.get_volume() # 输出体积结果 print("volume1:", volume) # 改变两个底面顶点顺序 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) volume = mesh.get_volume() print("volume2:", volume) # 改变两个底面和一个侧面顶点顺序 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) volume = mesh.get_volume() print("volume3:", volume) # 删掉一个底面 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) volume = mesh.get_volume() print("volume4:", volume) # 再删掉一个底面 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) volume = mesh.get_volume() print("volume5:", volume)