Plotly 3D表面图比例与轴值异常问题排查求助

你的代码存在两个核心问题，导致3D表面图的轴值和比例渲染异常：

1. X/Y轴取值错误 #

你在提取轴数值时，错误地将y和x都赋值为rownames(grid.df_matrix)。实际上acast(grid.df, X~Y, value.var="rkgls")生成的矩阵中，行对应X的取值，列对应Y的取值，正确的取值逻辑应该是：

x <- as.numeric(rownames(grid.df_matrix))
y <- as.numeric(colnames(grid.df_matrix))

之前的错误导致X、Y轴使用了同一组数值，完全不符合数据集的实际结构。

2. 3D轴比例未手动调整 #

你的X/Y轴实际范围是500米，但数值本身是264250这类大整数，而Z轴（高程）的范围仅约90米。Plotly默认的自动缩放会将Z轴压缩得几乎无法分辨地形起伏，必须手动设置轴比例：

自定义比例（推荐地形可视化） #

通过aspectmode = "manual"和aspectratio参数，让Z轴的视觉比例与X/Y轴的实际范围匹配：

layout(
  scene=list(
    xaxis=list(title='x (m)'),
    yaxis=list(title='y (m)'),
    zaxis=list(title='Elevation (m)'),
    # 按高程差与XY范围的比值设置Z轴比例，这里100/500=0.2
    aspectratio = list(x=1, y=1, z=0.2),
    aspectmode = "manual"
  )
)

可以根据实际高程差调整z的数值，确保地形起伏清晰可见。

等比例模式（适合精确对比） #

如果需要三个轴的单位长度视觉完全一致，可使用aspectmode = "cube"，但这会拉伸X/Y轴的显示范围。

完整修正代码 #

# 先对数据集按X、Y排序，避免矩阵行/列顺序混乱
grid.df <- grid.df[order(grid.df$X, grid.df$Y),]
grid.df_matrix <- acast(grid.df, X~Y, value.var="rkgls")

# 正确提取X/Y轴数值
x <- as.numeric(rownames(grid.df_matrix))
y <- as.numeric(colnames(grid.df_matrix))

# 绘制3D表面图
fig <- plot_ly(showscale = FALSE)
fig <- fig %>% add_surface(z = ~grid.df_matrix, x= x, y=y, colors = rainbow(100)) %>%
  layout(
    scene=list(
      xaxis=list(title='x (m)'),
      yaxis=list(title='y (m)'),
      zaxis=list(title='Elevation (m)'),
      aspectratio = list(x=1, y=1, z=0.2),
      aspectmode = "manual"
    )
  )

fig

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