我试图制作一个plot_ly曲面图,其中有两个横截面,即在Z平面上映射一个曲面,在X平面上映射另一个曲面。
z平面上的表面工作正常,但在另一个平面上添加第二个表面让我沮丧。我可能误解了plot_ly如何处理数据。
我正在尝试做类似这样的事情,但我无法从我自己的R(RStudio)中获取此剧情的代码 https://plot.ly/~schippkus/0/slices-in-volumetric-data/#plot
这是一个可重复的例子,单面工作,双面是奇数。此示例中创建的数据与我正在处理的数据具有相同的维度。它最初是一个三维数组,我强迫这种格式。
library(plotly); library(tidyverse);
x = 1:91
y = 1:109
z = 1:91
val = as.numeric(scale(sample(length(x)*length(y)*length(z))))
DATA = expand.grid(X=x,Y=y,Z=z) %>% as.data.frame() %>%
mutate(Val = val)
getTrace = function(data,plane,slice){
require(tidyverse)
tmp = data %>% filter(get(plane) %in% slice)
Trace = list(
x = tmp %>% select(X) %>% unique() %>% unlist() %>% unname(),
y = tmp %>% select(Y) %>% unique() %>% unlist() %>% unname(),
z = tmp %>% select(Z) %>% unique() %>% unlist() %>% unname()
)
switch(plane,
"Z"={
oneplane = expand.grid(x = Trace$x, y = Trace$y)
oneplane$z <- tmp %>% select(Val) %>% unlist() %>% unname()
Trace$SurfaceColor = as.matrix(spread(oneplane, key = x, value = z)[, -1]) %>% unname()
Trace$z <- matrix(slice, ncol=ncol(Trace$SurfaceColor), nrow=nrow(Trace$SurfaceColor))
Trace
},
"Y"={
oneplane = expand.grid(x = Trace$x, z = Trace$z)
oneplane$y <- tmp %>% select(Val) %>% unlist() %>% unname()
Trace$SurfaceColor = as.matrix(spread(oneplane, key = z, value = y)[, -1]) %>% unname()
Trace$y <- matrix(slice, ncol=ncol(Trace$SurfaceColor), nrow=nrow(Trace$SurfaceColor))
Trace
},
"X"={
oneplane = expand.grid(y = Trace$y, z = Trace$z)
oneplane$x <- tmp %>% select(Val) %>% unlist() %>% unname()
Trace$SurfaceColor = as.matrix(spread(oneplane, key = y, value = x)[, -1]) %>% unname()
Trace$x <- matrix(slice, ncol=ncol(Trace$SurfaceColor), nrow=nrow(Trace$SurfaceColor))
Trace
}
)
}
Trace1 = getTrace(DATA,"Z",45)
Trace2 = getTrace(DATA,"Y",45)
# This works as intended (except I can't seem to make the surface greyscale)
plot_ly() %>%
add_trace(x=Trace1$x,
y=Trace1$y,
z=Trace1$z,
type="surface",
cauto=T,
surfacecolor=Trace1$SurfaceColor,
colorscale=c("white","black")
) %>%
layout(title="Ok one-surface plot",
scene = list(
aspectratio = list(
x = 91/109,
y = 1,
z = 91/109
),
xaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)",
range=c(1,91)
),
yaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)",
range=c(1,109)
),
zaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
range = c(1, 91),
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)"
)
))
# This definately does not do as intended
plot_ly() %>%
add_trace(x=Trace1$x,
y=Trace1$y,
z=Trace1$z,
type="surface",
cauto=T,
surfacecolor=Trace1$SurfaceColor,
colorscale=c("white","black")
) %>%
add_trace(x=Trace2$x,
y=Trace2$y,
z=Trace2$z,
type="surface",
cauto=T,
surfacecolor=Trace2$SurfaceColor,
colorscale=c("white","black")
) %>%
layout(title="Bogus plot",
scene = list(
aspectratio = list(
x = 91/109,
y = 1,
z = 91/109
),
xaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)",
range=c(1,91)
),
yaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)",
range=c(1,109)
),
zaxis = list(
backgroundcolor = "rgb(230, 230,230)",
gridcolor = "rgb(255, 255, 255)",
range = c(1, 91),
showbackground = TRUE,
zerolinecolor = "rgb(255, 255, 255)"
)
))