9

我想创建一个球体的 3d 散点图,它们的颜色是第四维。我在 csv 文件中有数据,其中每行表示粒子的 x、y、z 位置,并且有一列告诉我粒子的值(1,2 或 3)。如果它们的值为 1,我想用一种颜色为球着色,否则用另一种颜色着色。

编辑:

我创建了以下代码:

library(rgl)
m <- read.csv(file="mem0.csv", sep = ",", head=TRUE)
mcol = m$val
i = 1 
mdim = dim(m)

while (i <= mdim[1] ){
   if (mcol[i] == 1){
      mcol[i] = "red"
   }else {
      mcol[i] = "blue"
   }
   i = i +1
}

plot3d(m$x, m$y, m$z, col = mcol, type='s', size=0.1)

编辑编号 2:

我使用 rgl.snapshot() 导出到 svg 文件:

我的 rgl.shanpshot 的快照

数据应该再次显示一层红球、4层蓝球和一层红球。

4

1 回答 1

17

包的plot3d()功能rgl允许很容易地做这样的事情。你甚至可以交互式地探索你的情节:

R> library(rgl)

R> df <- data.frame(x=runif(10,0,1),
+                  y=runif(10,0,1),
+                  z=runif(10,0,1),
+                  color=round(runif(10,1,3)))
R> df
            x         y          z color
1  0.73518229 0.1385970 0.69053482     2
2  0.88789302 0.6872121 0.54734176     2
3  0.79402546 0.5771570 0.89613292     1
4  0.19922140 0.2117405 0.25116078     1
5  0.31825325 0.7449661 0.01174593     2
6  0.64614521 0.4704698 0.68905621     1
7  0.15242295 0.6461338 0.77896858     1
8  0.32698024 0.4548752 0.33969754     3
9  0.00793849 0.6557488 0.75901935     2
10 0.20460232 0.9302882 0.23413984     3

你可以plot3d()这样调用:

R> plot3d(df$x, df$y, df$z, col=df$color, size=2, type='s')

这会给你类似的东西:

plot3d 结果

于 2011-03-07T15:27:02.387 回答