请给出python读取stl文件并计算文件stl三维颗粒模型的perimeter和 surface_area的代码及注释

以下是一种使用Python读取STL文件并计算其周长和表面积的代码：

import math

def calculate_perimeter(vertices):
    perimeter = 0
    for i in range(len(vertices)):
        # 计算当前顶点与下一个顶点之间的距离
        dx = vertices[i][0] - vertices[(i+1) % len(vertices)][0]
        dy = vertices[i][1] - vertices[(i+1) % len(vertices)][1]
        dz = vertices[i][2] - vertices[(i+1) % len(vertices)][2]
        distance = math.sqrt(dx**2 + dy**2 + dz**2)
        perimeter += distance
    return perimeter

def calculate_surface_area(vertices):
    surface_area = 0
    for i in range(len(vertices)):
        # 计算当前顶点与相邻两个顶点之间组成的两个向量
        vector1 = [
            vertices[i][0] - vertices[(i+1) % len(vertices)][0],
            vertices[i][1] - vertices[(i+1) % len(vertices)][1],
            vertices[i][2] - vertices[(i+1) % len(vertices)][2]
        ]
        vector2 = [
            vertices[i][0] - vertices[(i+2) % len(vertices)][0],
            vertices[i][1] - vertices[(i+2) % len(vertices)][1],
            vertices[i][2] - vertices[(i+2) % len(vertices)][2]
        ]
        # 计算两个向量的叉积
        cross_product = [
            vector1[1]*vector2[2] - vector1[2]*vector2[1],
            vector1[2]*vector2[0] - vector1[0]*vector2[2],
            vector1[0]*vector2[1] - vector1[1]*vector2[0]
        ]
        # 计算叉积的模长并累加到表面积
        surface_area += math.sqrt(cross_product[0]**2 + cross_product[1]**2 + cross_product[2]**2) / 2
    return surface_area

def read_stl_file(file_path):
    vertices = []
    with open(file_path, 'r') as file:
        for line in file:
            if line.startswith('vertex'):
                # 提取顶点坐标并添加到列表中
                vertex = [float(num) for num in line.strip().split()[1:]]
                vertices.append(vertex)
    return vertices

# 读取STL文件
stl_file_path = 'example.stl'
vertices = read_stl_file(stl_file_path)

# 计算周长和表面积
perimeter = calculate_perimeter(vertices)
surface_area = calculate_surface_area(vertices)

# 输出结果
print('Perimeter:', perimeter)
print('Surface Area:', surface_area)

注释：

calculate_perimeter函数用于计算给定顶点列表的周长。使用欧几里得距离公式计算顶点之间的距离，并将所有距离相加得到周长。
calculate_surface_area函数用于计算给定顶点列表的表面积。对于每个顶点，计算其与相邻两个顶点之间组成的两个向量，然后计算这两个向量的叉积。叉积的模长除以2即为该三角形的面积，将所有三角形的面积相加得到表面积。
read_stl_file函数用于读取STL文件，并提取顶点坐标。该函数打开STL文件，逐行读取文件内容，当行以"vertex"开头时，提取后面的数字作为顶点坐标，并将其添加到顶点列表中。
使用read_stl_file函数读取STL文件，获取顶点列表。
使用calculate_perimeter函数计算顶点列表的周长。
使用calculate_surface_area函数计算顶点列表的表面积。
输出周长和表面积的结果