当前位置:首页>java>「编程虚拟器官」R包ggseg绘制脑图谱

「编程虚拟器官」R包ggseg绘制脑图谱

  • 2026-10-23 12:03:19
「编程虚拟器官」R包ggseg绘制脑图谱
ggseg是款包含ggplot2几何对象(geom)的R包,可利用简单要素(simple features)绘制脑图谱。该包的核心组成部分是2个内置脑图谱的数据。通常,在脑区或脑网络上可视化分析结果时,需要在R等统计工具与脑成像软件之间切换以确保呈现准确。
而ggseg包旨在实现直接通过R绘制分析结果,简化可视化流程。
图谱
目前包中提供4款脑图谱:
  1. dk:Desikan-Killiany图谱(aparc版本,皮层图谱)
  2. aseg:自动皮层下分割图谱(subcortical segmentation)
注意:自1.5.3版本起,ggseg已拆分为2个独立包:
在ggplot中绘制2D多边形图的ggseg(本包)
通过plotly绘制3D网格图的ggseg3d(需前往ggseg3d代码仓库获取)
https://github.com/ggseg/ggseg3d
更多图谱及自定义图谱构建函数可在配套包ggsegExtra中查找
https://github.com/ggseg/ggsegExtra
安装
🍀 从CRAN安装
install.packages("ggseg")
🍀 通过ggseg r-universe安装
# 启用 r-universe 仓库options(repos = c(    ggseg = 'https://ggseg.r-universe.dev',    CRAN = 'https://cloud.r-project.org'))# 安装包install.packages('ggseg')
🍀 安装开发版(需remotes包)
install.packages("remotes")remotes::install_github("ggseg/ggseg")
安装完成后,可在需要时加载包。所有函数均遵循标准R文档规范。
使用
绘制内置图谱
library(ggseg)library(ggplot2)plot(dk)  # 绘制 Desikan-Killiany 皮层图谱
plot(aseg)  # 绘制自动皮层下分割图谱
投射自定义数据至脑图谱
library(dplyr)# 构建示例数据some_data <- tibble(  region = rep(c("transverse temporal", "insula", "precentral", "superior parietal"), 2),  p = sample(seq(0, .5, .001), 8),  groups = c(rep("g1", 4), rep("g2", 4)))# 绘图:按分组 facet,填充色映射数据值 psome_data |>  group_by(groups) |>  ggplot() +  geom_brain(atlas = dk,             position = position_brain(hemi ~ side),             aes(fill = p)) +  facet_wrap(~groups)#> merging atlas and data by 'region'(按“脑区”列合并图谱与数据)

包内置多个教程文档(vignettes)

https://ggseg.github.io/ggseg/

查看包开发者博客的使用入门教程

https://drmowinckels.io/blog/2021/new-ggseg-with-geom/

参考

https://ggseg.github.io/ggseg/

最新文章

随机文章