解决Vedo中网格法向量计算偏移问题

在使用Vedo库进行3D可视化时，当网格对象不在坐标系原点时，可能会遇到法向量计算不正确的问题。本文将详细分析这一现象的原因，并提供正确的解决方案。

问题现象

当用户尝试为一个偏移了原点的球体网格计算法向量并绘制箭头时，发现箭头方向并不垂直于球体表面。具体表现为：球体中心位于(100,0,0)位置时，计算得到的法向量箭头方向异常。

原因分析

这种现象的根本原因在于计算法向量时没有考虑网格对象的实际位置偏移。在3D图形学中，法向量是相对于对象局部坐标系定义的，而直接使用顶点坐标计算时，如果不考虑对象的变换矩阵，就会得到错误的结果。

解决方案

Vedo库提供了内置的compute_normals()方法来正确处理这种情况。正确用法如下：

1. 首先创建网格对象并设置其位置：

s1 = Sphere(pos=[100,0,0], r=10, res=8).wireframe().c('white')

2. 调用compute_normals()方法计算顶点法向量：

s1.compute_normals(points=True, cells=False)

3. 获取顶点坐标和法向量：

start = s1.vertices
end = start + s1.vertex_normals * (5 + np.random.randn(len(start), 3))

4. 创建箭头对象并可视化：

a1 = Arrows(start, end, c="viridis", s=2, res=24).lighting("plastic")
show(s1, a1)

技术要点

1. 顶点法向量：vertex_normals属性存储了每个顶点的法向量，这些向量已经考虑了对象的变换矩阵。

2. 箭头可视化：使用Arrows类可以直观地显示法向量方向，通过设置颜色映射可以更好地观察法向量的大小变化。

3. 光照效果：添加lighting("plastic")可以使箭头具有更好的3D视觉效果。

最佳实践

• 对于任何变换过的网格对象，都应使用内置的compute_normals()方法而非手动计算
• 可视化法向量时，可以添加随机分量来避免箭头重叠
• 使用颜色映射可以同时展示法向量的大小信息

通过这种方法，可以确保无论网格对象位于坐标系中的任何位置，都能正确计算和显示其法向量。