ggseg是款包含ggplot2几何对象(geom)的R包,可利用简单要素(simple features)绘制脑图谱。该包的核心组成部分是2个内置脑图谱的数据。通常,在脑区或脑网络上可视化分析结果时,需要在R等统计工具与脑成像软件之间切换以确保呈现准确。而ggseg包旨在实现直接通过R绘制分析结果,简化可视化流程。图谱- dk:Desikan-Killiany图谱(aparc版本,皮层图谱)
- 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/ggsegExtrainstall.packages("ggseg")
# 启用 r-universe 仓库options(repos = c( ggseg = 'https://ggseg.r-universe.dev', CRAN = 'https://cloud.r-project.org'))# 安装包install.packages('ggseg')
install.packages("remotes")remotes::install_github("ggseg/ggseg")
安装完成后,可在需要时加载包。所有函数均遵循标准R文档规范。library(ggseg)library(ggplot2)plot(dk) # 绘制 Desikan-Killiany 皮层图谱
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/