STL 三维颗粒模型几何形状特征计算 Python 代码

要计算 STL 三维颗粒模型的圆度、球形度、内径比、矩形度和棱角性，需要进行一系列的几何计算。以下是一种可能的实现方式：\n\npython\nimport numpy as np\nfrom scipy.spatial import ConvexHull\n\ndef calculate_roundness(vertices):\n # 计算圆度\n hull = ConvexHull(vertices)\n volume = hull.volume\n surface_area = hull.area\n roundness = (36 * np.pi * volume) / (surface_area ** 2)\n return roundness\n\ndef calculate_sphericity(vertices):\n # 计算球形度\n hull = ConvexHull(vertices)\n volume = hull.volume\n sphere_volume = (4/3) * np.pi * (hull.equations[0, 3] ** 3)\n sphericity = volume / sphere_volume\n return sphericity\n\ndef calculate_inner_diameter_ratio(vertices):\n # 计算内径比\n hull = ConvexHull(vertices)\n max_distance = 0\n for i, vertex_i in enumerate(hull.vertices[:-1]):\n for vertex_j in hull.vertices[i+1:]:\n distance = np.linalg.norm(vertices[vertex_i] - vertices[vertex_j])\n if distance > max_distance:\n max_distance = distance\n inner_diameter = max_distance\n diameter = np.sqrt(hull.equations[0, 3])\n inner_diameter_ratio = inner_diameter / diameter\n return inner_diameter_ratio\n\ndef calculate_rectangularity(vertices):\n # 计算矩形度\n hull = ConvexHull(vertices)\n volume = hull.volume\n rectangular_volume = hull.area * (3/2 * np.mean(hull.equations[:, 3]))\n rectangularity = volume / rectangular_volume\n return rectangularity\n\ndef calculate_edge_sharpness(vertices):\n # 计算棱角性\n hull = ConvexHull(vertices)\n num_edges = len(hull.simplices)\n num_vertices = len(hull.vertices)\n edge_sharpness = num_edges / num_vertices\n return edge_sharpness\n\n# 示例使用\n# 读取STL文件模型并提取顶点信息\nvertices = np.array([[0, 0, 0], [1, 0, 0], [1, 1, 0], [0, 1, 0], [0, 0, 1], [1, 0, 1], [1, 1, 1], [0, 1, 1]])\nroundness = calculate_roundness(vertices)\nsphericity = calculate_sphericity(vertices)\ninner_diameter_ratio = calculate_inner_diameter_ratio(vertices)\nrectangularity = calculate_rectangularity(vertices)\nedge_sharpness = calculate_edge_sharpness(vertices)\n\nprint("圆度:", roundness)\nprint("球形度:", sphericity)\nprint("内径比:", inner_diameter_ratio)\nprint("矩形度:", rectangularity)\nprint("棱角性:", edge_sharpness)\n\n\n需要注意的是，这里假设STL文件模型中的顶点信息以Numpy数组的形式提供。如果实际情况不同，可以根据实际情况进行调整。此外，代码中的计算方法仅提供了一种可能的实现方式，具体的计算公式可以根据需求进行调整。