没有什么特殊的算法,直接用包,开箱即用
from stl import mesh import numpy as np # 读取stl文件 filename = './text.stl' mesh_data = mesh.Mesh.from_file(filename) xyz = (mesh_data.max_ - mesh_data.min_) sizel = round(xyz[0] / 10, 2) sizew = round(xyz[1] / 10, 2) sizeh = round(xyz[2] / 10, 2) # 计算体积 volume, _, _ = mesh_data.get_mass_properties() # 计算表面积 surface_area = 0.0 for triangle in mesh_data.vectors: # 获取三角形的三个顶点 a = triangle[0] b = triangle[1] c = triangle[2] # 计算三角形的面积并累加 surface_area += 0.5 * np.linalg.norm(np.cross(b - a, c - a)) print("Surface Area:", surface_area) print("Volume:", volume)
下面是PHP计算体积和表面积 应该是顶点问题 计算有bug 等有空再来修复 仅做参考
<?php // 从STL文件中解析顶点信息 // 从二进制STL文件中解析顶点信息 function parseBinarySTLFile($filePath) { $vertices = array(); // 读取STL文件内容 $stlContent = file_get_contents($filePath); //var_dump($stlContent); // 假设STL文件头部80字节为文件头信息,跳过这部分 $stlContent = substr($stlContent, 80); // 之后的4字节为三角形数量,按little-endian解析 $numTriangles = unpack('V', substr($stlContent, 0, 4))[1]; $stlContent = substr($stlContent, 4); // 依次解析每个三角形的顶点信息 for ($i = 0; $i < $numTriangles; $i++) { // 每个三角形占据50个字节,包括面法向量和三个顶点坐标 $data = unpack('f*', substr($stlContent, 0, 50)); $vertices[] = array( [$data[3], $data[4], $data[5]], [$data[6], $data[7], $data[8]], [$data[9], $data[10], $data[11]] ); $stlContent = substr($stlContent, 50); } return $vertices; } // 计算模型体积 function calculateModelVolume($vertices) { $volume = 0; // 对于每个三角形,计算有向体积并累加 foreach ($vertices as $triangle) { // 以第一个顶点为起点,计算有向体积 $v321 = $triangle[2][0] * $triangle[1][1] * $triangle[0][2]; $v231 = $triangle[1][0] * $triangle[2][1] * $triangle[0][2]; $v312 = $triangle[2][0] * $triangle[0][1] * $triangle[1][2]; $v132 = $triangle[0][0] * $triangle[2][1] * $triangle[1][2]; $v213 = $triangle[1][0] * $triangle[0][1] * $triangle[2][2]; $v123 = $triangle[0][0] * $triangle[1][1] * $triangle[2][2]; $volume += (1.0 / 6.0) * (-$v321 + $v231 + $v312 - $v132 - $v213 + $v123); //return (-v321 + v231 + v312 - v132 - v213 + v123)/6.0; } return abs($volume); } // 计算模型表面积 function calculateModelSurfaceArea($vertices) { $surfaceArea = 0; // 对于每个三角形,计算面积并累加 foreach ($vertices as $triangle) { // 计算三角形的两条边的向量 $v1 = array( $triangle[1][0] - $triangle[0][0], $triangle[1][1] - $triangle[0][1], $triangle[1][2] - $triangle[0][2] ); $v2 = array( $triangle[2][0] - $triangle[0][0], $triangle[2][1] - $triangle[0][1], $triangle[2][2] - $triangle[0][2] ); // 计算三角形的面积(使用向量叉乘) $crossProduct = array( $v1[1] * $v2[2] - $v1[2] * $v2[1], $v1[2] * $v2[0] - $v1[0] * $v2[2], $v1[0] * $v2[1] - $v1[1] * $v2[0] ); $triangleArea = 0.5 * sqrt(pow($crossProduct[0], 2) + pow($crossProduct[1], 2) + pow($crossProduct[2], 2)); // 累加三角形的面积 $surfaceArea += $triangleArea; } return $surfaceArea; } function calculateModelDimensions($vertices) { $minX = $maxX = $vertices[0][0][0]; $minY = $maxY = $vertices[0][0][1]; $minZ = $maxZ = $vertices[0][0][2]; foreach ($vertices as $triangle) { foreach ($triangle as $vertex) { $minX = min($minX, $vertex[0]); $maxX = max($maxX, $vertex[0]); $minY = min($minY, $vertex[1]); $maxY = max($maxY, $vertex[1]); $minZ = min($minZ, $vertex[2]); $maxZ = max($maxZ, $vertex[2]); } } $width = $maxX - $minX; $height = $maxY - $minY; $depth = $maxZ - $minZ; return array('width' => $width, 'height' => $height, 'depth' => $depth); } // STL文件路径 $stlFilePath = '111.stl'; // 解析STL文件并获取顶点信息 $vertices = parseBinarySTLFile($stlFilePath); //var_dump($vertices); // 计算模型体积 $modelVolume = calculateModelVolume($vertices); // 表面积 //$area = calculateModelSurfaceArea($vertices); // 计算模型的长宽高尺寸 //$modelDimensions = calculateModelDimensions($vertices); // 输出模型的长宽高尺寸 //echo "Model Dimensions: Width = " . $modelDimensions['width'] . ", Height = " . $modelDimensions['height'] . ", Depth = " . $modelDimensions['depth']; //echo $area; // 输出模型体积 echo "The volume of the STL model is: " . $modelVolume; ?>
标签:表面积,area,STL,vertices,python,mesh,计算,data From: https://www.cnblogs.com/we-jack/p/17817627.html