如何平滑pyvista.StructuredGrid，消除霍普夫环面渲染着色断层

问题描述 #

以下是使用Python的PyVista生成的霍普夫环面（Hopf torus）：

现有生成代码如下：

import numpy as np
import pyvista as pv

A = 0.44
n = 3
def Gamma(t):
    alpha = np.pi/2 - (np.pi/2-A)*np.cos(n*t)
    beta = t + A*np.sin(2*n*t)
    return np.array([
      np.sin(alpha) * np.cos(beta),
      np.sin(alpha) * np.sin(beta),
      np.cos(alpha)
    ])

def HopfInverse(p, phi):
    return np.array([
      (1+p[2])*np.cos(phi),
      p[0]*np.sin(phi) - p[1]*np.cos(phi), 
      p[0]*np.cos(phi) + p[1]*np.sin(phi),
      (1+p[2])*np.sin(phi)
    ]) / np.sqrt(2*(1+p[2]))

def Stereo(q):
    return 2*q[0:3] / (1-q[3])

def F(t, phi):
    return Stereo(HopfInverse(Gamma(t), phi))

angle = np.linspace(0, 2 * np.pi, 300)
theta, phi = np.meshgrid(angle, angle)
x, y, z = F(theta, phi)

# Display the mesh
grid = pv.StructuredGrid(x, y, z)
grid.plot(smooth_shading=True)

当前渲染效果的色彩过渡不够平滑，在右下角的波瓣处可见一条分隔浅灰色与深灰色的线条，需要消除这条分隔线，实现pyvista.StructuredGrid的平滑渲染。

解决方案 #

出现该硬边的核心原因是环面属于闭合周期性曲面，现有代码生成网格时包含了theta/phi的重复端点，且默认法线计算没有处理周期边界的连续性，导致平滑着色时相邻重合点的法线不一致，出现明暗分界。可按以下步骤修复：

• 移除重复端点：生成角度序列时设置endpoint=False，避免2π和0位置的重复顶点，保证网格天然连续：

  # 替换原角度生成代码
  angle = np.linspace(0, 2 * np.pi, 300, endpoint=False)
  

• 手动计算连续法线：调用compute_normals方法，关闭边拆分、开启法线方向统一，保证周期边界处法线过渡平滑：

  # 创建网格后、渲染前添加法线计算逻辑
  grid = pv.StructuredGrid(x, y, z)
  grid = grid.compute_normals(splitting=False, consistent_normals=True)
  

• 可选优化渲染效果：开启抗锯齿、选择过渡更平滑的配色，进一步弱化边界感：

  grid.plot(smooth_shading=True, cmap='viridis', anti_aliasing='msaa')
  

内容的提问来源于stack exchange，提问作者Stéphane Laurent