如何在rgl::surface3d中增大网格线间距？适配小步长数据

问题场景 #

当使用rgl包绘制3D曲面时，如果x、y序列的步长极小（比如length=100），直接用surface3d绘制前后网格线会导致线条过密，影响可视化效果。我们需要在不更换包的前提下，实现类似plot3D::persp3D中resfac参数的功能，控制网格线的间距。

解决方案思路 #

核心是分离曲面绘制和网格线绘制：用完整的原始数据绘制曲面保证平滑度，对数据进行抽稀后再绘制网格线，通过抽稀间隔控制网格线的密度。

修改后的代码实现 #

# 带网格线密度控制的绘图函数
my_plot <- function(x, y, z, grid_res = 10){
  # 绘制完整曲面，保证细节
  colours <- heat.colors(length(unique(z)))
  minval <- min(z)
  maxval <- max(z) + 10
  col <- colours[(z - minval)/(maxval - minval)*99 + 1]
  rgl::persp3d(x, y, z, col = col) 
  
  # 对数据抽稀，控制网格线密度
  # 生成抽稀的索引间隔
  x_idx <- seq(1, length(x), by = grid_res)
  y_idx <- seq(1, length(y), by = grid_res)
  
  # 抽取对应的子数据
  x_sub <- x[x_idx]
  y_sub <- y[y_idx]
  z_sub <- z[x_idx, y_idx]
  
  # 用抽稀后的数据绘制网格线
  rgl::surface3d(x_sub, y_sub, z_sub, front = "lines")
  rgl::surface3d(x_sub, y_sub, z_sub, back = "lines")
}

# 测试小步长数据（length=100）
x <- seq(-2, 2, length = 100) 
y <- seq(-3, 3, length = 100)
z <- outer(x, y, function(x, y) x + y^2)

# 绘图，grid_res=10表示每10个点取一个绘制网格线
my_plot(x, y, z, grid_res = 10)

参数说明 #

• grid_res：控制网格线的抽稀间隔，数值越大，网格线越稀疏（间距越宽）。比如设为5表示每5个原始点取一个来绘制网格线，可根据实际数据的步长调整。
• 曲面依然使用完整的原始数据绘制，不会丢失曲面的细节，仅网格线进行抽稀处理。

效果验证 #

当grid_res设置合适时，即使x、y的步长极小，网格线也能保持合适的间距，既清晰展示网格结构，又不会因为过密导致画面杂乱。

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