在R中实现三维空间体素密度绘图的技术问询

你可以通过数据预处理+rgl包绘图实现需求，没有完全开箱即用的工具，但步骤清晰，具体如下：

1. 数据预处理 #

先把原始数据整理为每行对应一个(x,y,z)记录的长格式，再统计每个体素的观测点数：

# 模拟示例数据（替换为你的真实数据）
set.seed(123)
raw_data <- data.frame(
  x = sample(1:5, 200, replace = TRUE),
  y = sample(1:5, 200, replace = TRUE),
  z = sample(1:3, 200, replace = TRUE),
  obs_data = rnorm(200)
)

# 统计每个(x,y,z)组合的观测次数
voxel_counts <- table(raw_data$x, raw_data$y, raw_data$z)
# 转为数组格式适配rgl
voxel_array <- array(unlist(voxel_counts), dim = dim(voxel_counts))

2. 定义颜色分箱规则 #

根据你的需求创建颜色映射：

library(dplyr)

# 颜色映射函数
get_voxel_color <- function(count) {
  case_when(
    count == 0 ~ "#00000000",  # 透明
    count %in% 1:10 ~ "#00FF00",  # 绿色
    count %in% 11:50 ~ "#FFFF00",  # 黄色
    count > 50 ~ "#FF0000"  # 红色
  )
}

# 生成每个体素的颜色矩阵
voxel_colors <- apply(voxel_array, c(1,2,3), get_voxel_color)

3. 绘制三维体素图 #

使用rgl::voxels()函数实现体素的颜色填充，支持交互式查看：

library(rgl)

# 绘制体素
voxels(
  x = 1:dim(voxel_array)[1],
  y = 1:dim(voxel_array)[2],
  z = 1:dim(voxel_array)[3],
  fill = voxel_colors,
  xlab = "X", ylab = "Y", zlab = "Z"
)

# 导出交互式窗口（在RStudio中可直接查看）
rglwidget()

补充说明 #

• 若你的x/y/z是连续坐标，需要先用cut()函数将空间划分为固定大小的体素，再统计每个体素内的点数
• 透明体素会自动隐藏，无需额外处理
• rgl支持旋转、缩放等交互式操作，方便查看三维空间分布

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