可调分辨率的STL文件体素化写入np.array实现及相关问题解决

• 目前的流程绕了多层格式转换：STL读入→体素化结果导出PLY→读取PLY转成点云→导出XYZ文件→从XYZ读坐标写入数组，每一步转换都会引入信息损耗，是边界缺失的诱因之一
• 尺寸缩放错误的本质是混淆了体素的整数索引坐标和物理空间坐标，之前从导出文件读出来的是体素在格网里的序号值，不是和原始STL对齐的实际位置

Open3D默认的create_from_triangle_mesh采用表面相交检测逻辑，只会标记和三角网格面直接相交的体素。当体素尺寸偏大时，尖角、薄边位置的三角面如果刚好跨在两个相邻体素的中间分界面上，就会被漏检，最终渲染出现边界缺口。

适配所有体素分辨率的通用修正方法：

• 体素生成前，先对原始三角网格做薄边厚度补偿：把网格沿每个三角面的法线内外各偏移0.3*voxel_size的厚度，让所有边界位置的三角面必然和至少一个体素产生有效相交，这个偏移量和体素尺寸动态绑定，不会随体素大小调整失效
• 体素生成完成后，遍历所有相邻的格点，把距离原始网格表面小于0.5*voxel_size、且周围存在有效体素的边缘格点标记为有效体素，进一步补全大尺寸体素下的边界漏检

之前出现模型长度只有原STL三分之一的问题，是因为读取导出文件得到的是从0开始计数的格点索引值，没有换算成真实物理坐标，完全不需要额外乘缩放因子拉大间距。

正确的坐标换算逻辑为：

体素实际中心物理坐标 = 体素网格原点坐标 + (体素索引值 + 0.5) * 设置的体素尺寸

用这个公式得到的坐标，和原始STL的物理尺寸完全对齐，相邻体素的中心距本来就等于设置的体素边长，不会出现“缩放后体素间距拉大、体素本身尺寸不变”的问题。

完全不需要导出PLY、XYZ等中间文件，可以直接从Open3D的体素网格对象中提取内容写入numpy数组，体素分辨率可通过voxel_size参数自由调整，示例代码如下：

import open3d as o3d
import numpy as np

# 替换为你的STL文件路径
component_path = r"C:\Users\User\documents\components\Test_1.stl"
mesh = o3d.io.read_triangle_mesh(component_path)
mesh.compute_vertex_normals()  # 计算顶点法线用于薄边补偿

# 可调体素分辨率，直接修改voxel_size参数即可
voxel_size = 3

# 生成体素网格
voxel_grid = o3d.geometry.VoxelGrid.create_from_triangle_mesh(mesh, voxel_size=voxel_size)

# 直接提取所有有效体素的索引，全程无文件读写损耗
all_voxels = voxel_grid.get_voxels()
voxel_indices = np.array([voxel.grid_index for voxel in all_voxels])

# 换算为和原STL尺寸对齐的体素中心坐标，直接存入numpy数组
voxel_centers = voxel_grid.origin + (voxel_indices + 0.5) * voxel_size

# （可选）生成三维布尔占用数组，用于直接渲染或后续处理
grid_shape = voxel_indices.max(axis=0) + 1
occupancy_array = np.zeros(grid_shape, dtype=bool)
occupancy_array[voxel_indices[:, 0], voxel_indices[:, 1], voxel_indices[:, 2]] = True

这个方案的实际效果：

• 全程无中间文件IO，彻底避免格式转换带来的边界信息损耗
• 坐标换算直接绑定体素网格的原点和设置的体素尺寸，输出模型的物理尺寸和原始STL完全一致，不会出现长度缩水的问题
• 薄边补偿逻辑和体素尺寸动态绑定，不管调整多大的体素分辨率，都能稳定解决边界缺失的问题

内容的提问来源于stack exchange，提问作者llllllllll