Python 3d分析(vedo)：找到适合网格内部的最大球体？

from vedo import *
import numpy as np

mesh = Mesh("your_mesh_file.obj")

bbox = mesh.bounds()

def is_sphere_inside(mesh, center, radius):
    # 将球体的每个顶点与网格进行碰撞检测
    for vertex in mesh.points():
        if np.linalg.norm(vertex - center) > radius:
            return False
    return True

def find_max_sphere_inside(mesh):
    # 初始化最大半径和球体中心
    max_radius = 0
    max_center = None
    
    # 遍历网格的边界框
    for x in np.linspace(bbox[0], bbox[1], 100):
        for y in np.linspace(bbox[2], bbox[3], 100):
            for z in np.linspace(bbox[4], bbox[5], 100):
                center = np.array([x, y, z])
                radius = 0
                
                # 逐渐增加球体的半径，直到不再适合网格内部
                while is_sphere_inside(mesh, center, radius):
                    radius += 0.1
                
                # 更新最大半径和球体中心
                if radius > max_radius:
                    max_radius = radius
                    max_center = center
    
    return max_center, max_radius

center, radius = find_max_sphere_inside(mesh)

show(mesh, at=0)
show(Sphere(center, radius), at=1)