STL 三维颗粒模型形状参数计算代码：圆度、球形度、内径比、矩形度和棱角性

以下是一段计算 STL 三维颗粒模型的圆度、球形度、内径比、矩形度和棱角性的示例代码：\n\npython\nimport numpy as np\n\ndef calculate_circularity(triangles):\n total_area = 0\n total_perimeter = 0\n\n for triangle in triangles:\n a = np.linalg.norm(triangle[0] - triangle[1])\n b = np.linalg.norm(triangle[1] - triangle[2])\n c = np.linalg.norm(triangle[2] - triangle[0])\n\n perimeter = a + b + c\n area = np.sqrt(perimeter/2 * (perimeter/2 - a) * (perimeter/2 - b) * (perimeter/2 - c))\n\n total_area += area\n total_perimeter += perimeter\n\n circularity = np.pi * 4 * total_area / (total_perimeter ** 2)\n return circularity\n\ndef calculate_sphericity(triangles):\n total_area = 0\n\n for triangle in triangles:\n a = np.linalg.norm(triangle[0] - triangle[1])\n b = np.linalg.norm(triangle[1] - triangle[2])\n c = np.linalg.norm(triangle[2] - triangle[0])\n\n area = np.sqrt((a + b + c) * (a + b - c) * (b + c - a) * (c + a - b)) / 4\n\n total_area += area\n\n sphericity = (np.pi ** (1/3) * 6 * (total_area ** (2/3))) / (np.sqrt(3) * (total_area ** (1/3)))\n return sphericity\n\ndef calculate_inner_diameter_ratio(triangles):\n total_diameter = 0\n total_inner_diameter = 0\n\n for triangle in triangles:\n a = np.linalg.norm(triangle[0] - triangle[1])\n b = np.linalg.norm(triangle[1] - triangle[2])\n c = np.linalg.norm(triangle[2] - triangle[0])\n\n diameter = (a + b + c) / 3\n inner_diameter = np.sqrt(2 * (a ** 2 + b ** 2 + c ** 2) - (a ** 4 + b ** 4 + c ** 4) / (a ** 2 + b ** 2 + c ** 2))\n\n total_diameter += diameter\n total_inner_diameter += inner_diameter\n\n inner_diameter_ratio = total_inner_diameter / total_diameter\n return inner_diameter_ratio\n\ndef calculate_rectangularity(triangles):\n total_area = 0\n total_perimeter = 0\n\n for triangle in triangles:\n a = np.linalg.norm(triangle[0] - triangle[1])\n b = np.linalg.norm(triangle[1] - triangle[2])\n c = np.linalg.norm(triangle[2] - triangle[0])\n\n perimeter = a + b + c\n area = np.sqrt(perimeter/2 * (perimeter/2 - a) * (perimeter/2 - b) * (perimeter/2 - c))\n\n total_area += area\n total_perimeter += perimeter\n\n rectangularity = total_area / (total_perimeter ** 2)\n return rectangularity\n\ndef calculate_edge_angle(triangles):\n total_angle = 0\n\n for triangle in triangles:\n a = np.linalg.norm(triangle[0] - triangle[1])\n b = np.linalg.norm(triangle[1] - triangle[2])\n c = np.linalg.norm(triangle[2] - triangle[0])\n\n angle_a = np.arccos((b ** 2 + c ** 2 - a ** 2) / (2 * b * c))\n angle_b = np.arccos((c ** 2 + a ** 2 - b ** 2) / (2 * c * a))\n angle_c = np.arccos((a ** 2 + b ** 2 - c ** 2) / (2 * a * b))\n\n total_angle += angle_a + angle_b + angle_c\n\n edge_angle = total_angle / (len(triangles) * 3)\n return edge_angle\n\n# 示例使用：\ntriangles = [[[0, 0, 0], [1, 0, 0], [0, 1, 0]], [[1, 0, 0], [1, 1, 0], [0, 1, 0]]] # 替换为STL模型的三角面片坐标\ncirc = calculate_circularity(triangles)\nspher = calculate_sphericity(triangles)\nidr = calculate_inner_diameter_ratio(triangles)\nrect = calculate_rectangularity(triangles)\nedge = calculate_edge_angle(triangles)\n\nprint("Circularity:", circ)\nprint("Sphericity:", spher)\nprint("Inner Diameter Ratio:", idr)\nprint("Rectangularity:", rect)\nprint("Edge Angle:", edge)\n\n\n请注意，此代码假设三角面片坐标的顺序是逆时针方向的，并且使用了NumPy库进行向量和矩阵计算。在实际使用时，需要根据STL模型的具体格式和数据结构进行适当的修改。