R中的体积空间中的切片

时间:2018-01-25 19:23:53

标签: r 3d plotly r-plotly

我试图制作一个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)"
               )
             ))

1 个答案:

答案 0 :(得分:1)

在Trace2中,Get http://IP:17428/update?id=sigfoxtest&data={data} 是一个整数。而且我相信,两条轨迹中的z都需要是一个矩阵才能工作。

z

因此,使用两个> class(Trace2$y) [1] "matrix" > class(Trace2$z) [1] "integer" 矩阵的简单示例,它可以工作。在第一个矩阵中,所有z值都相同(matrix2 = 6)。但在第二个中,z值是常量,但y值不同。

z

给出这个:

enter image description here