在绘图表面上使用条件着色

Question

我第一次通过 R 使用 plotly 并尝试从网格创建一个表面并根据计算为其着色。

例如，我想使用来自data(volcano)的表面，如

library(plotly)
plot_ly(z = ~volcano) %>% add_surface()

但不是基于 z 值（高度）的颜色，假设我想根据在 (20,60) 处的小台地上与我家的距离进行颜色。

house_loc <- c(20,60,150) # (x,y,z) of my house
dist_to_house <- Vectorize(function(x,y,z){sqrt(sum( (c(x,y,z)-house_loc)^2 ))})

到目前为止我已经尝试过：

color_me <-function(x){
  colorRampPalette(c('tan','blue')
  )(24L)[findInterval(x,seq(0,1,length.out=25),
                      all.inside=TRUE)]
} 

library(dplyr)
library(reshape2)
volcano %>% 
     melt( varnames=c('y','x'),value.name='z' ) %>%
     mutate( d = dist_to_house(x, y, z) ,
             d_rel = d/max(d),
             d_color = color_me(d_rel) 
     ) -> df

plot_ly(df,
        type='scatter3d',
        mode='none', # no markers, just surface
        x=~x,
        y=~y,
        z=~z,
        surfaceaxis=2,
        surfacecolor=~d_color) # last argument seems not to work

刚刚返回：

理想的结果是将房屋区域的景观染成棕褐色，并在远离房屋的区域逐渐褪色为蓝色。

related question 使用了在其他地方找到的 mesh3d 代码，并没有解释如何计算 (i, j, k)

Answer 1

您的代码几乎拥有您需要的一切，只需使用 surface 绘图并将距离数组用作 color。

library(plotly)
library(dplyr)
library(reshape2)

house_loc <- c(20,60,150)
dist_to_house <- Vectorize(function(x,y,z){sqrt(sum( (c(x,y,z)-house_loc)^2 ))})

volcano %>% 
  melt( varnames=c('y','x'),value.name='z' ) %>%
  mutate( d = dist_to_house(x, y, z) ,
          d_rel = d/max(d)
  ) -> df

color <- df$d_rel
dim(color) <- dim(volcano)

plot_ly(df,
        type='surface',
        z=volcano,
        surfacecolor=color,
        colors=c('tan','blue'))

Answer 2

除了surface 绘图（请参阅已接受的答案）之外，我们还可以执行mesh3d 绘图并避免绘图所需的重塑（回到网格）步骤。

但是，比例尺仍然不正确（显示z 的范围，而不是d_rel）

plot_ly(df,
        type='mesh3d',
        x = ~x,
        y = ~y,
        z = ~z,
        intensity=~d_rel,
        colors = colorRamp(c("tan", "blue"))
    )

与直觉相反，似乎控制条件着色的是 intensity= 而不是 color=。

我最初避免使用mesh3d，因为我认为我必须创建一个三角形网格（Delaunay 之类的）并且不知道如何做到这一点，但在这种情况下它似乎是自动处理的。