The brainglobe atlas API (brainglobe-atlasapi) provides a common interface for programmers to download and process brain atlas data from multiple sources.
A number of atlases are in development, but those available currently are:
- Allen Mouse Brain Atlas at 10, 25, 50 and 100 micron resolutions
- Allen Human Brain Atlas at 100 micron resolution
- Max Planck Zebrafish Brain Atlas at 1 micron resolution
- Enhanced and Unified Mouse Brain Atlas at 10, 25, 50 and 100 micron resolutions
- Smoothed version of the Kim et al. mouse reference atlas at 10, 25, 50 and 100 micron resolutions
- Gubra's LSFM mouse brain atlas at 20 micron resolution
- 3D version of the Allen mouse spinal cord atlas at 20 x 10 x 10 micron resolution
- AZBA: A 3D Adult Zebrafish Brain Atlas at 4 micron resolution
- Waxholm Space atlas of the Sprague Dawley rat brain at 39 micron resolution
- 3D Edge-Aware Refined Atlases Derived from the Allen Developing Mouse Brain Atlases (E13, E15, E18, P4, P14, P28 & P56)
- Princeton Mouse Brain Atlas at 20 micron resolution
- Kim Lab Developmental CCF (P56) at 10 micron resolution with 8 reference images - STP, LSFM (iDISCO) and MRI (a0, adc, dwo, fa, MTR, T2)
brainglobe-atlasapi works with Python >3.6, and can be installed from PyPI with:
pip install brainglobe-atlasapi
Full information can be found in the documentation
To see a list of atlases use brainglobe_atlasapi.show_atlases
from brainglobe_atlasapi import show_atlases
show_atlases()
# Brainglobe Atlases
# â•â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”¬â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”¬â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”¬â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â”€â•®
# │ Name │ Downloaded │ Local version │ Latest │
# │ │ │ │ version │
# ├──────────────────────────────────┼────────────┼───────────────┼──────────────┤
# │ allen_human_500um │ ✔ │ 0.1 │ 0.1 │
# │ mpin_zfish_1um │ ✔ │ 0.3 │ 0.3 │
# │ allen_mouse_50um │ ✔ │ 0.3 │ 0.3 │
# │ kim_unified_25um │ ✔ │ 0.1 │ 0.1 │
# │ allen_mouse_25um │ ✔ │ 0.3 │ 0.3 │
# │ allen_mouse_10um │ ✔ │ 0.3 │ 0.3 │
# │ example_mouse_100um │ --- │ --- │ 0.3 │
# ╰──────────────────────────────────┴────────────┴───────────────┴──────────────╯
All the features of each atlas can be accessed via the BrainGlobeAtlas
class.
e.g. for the 25um Allen Mouse Brain Atlas:
from brainglobe_atlasapi.bg_atlas import BrainGlobeAtlas
atlas = BrainGlobeAtlas("allen_mouse_25um")
The various files associated with the atlas can then be accessed as attributes of the class:
# reference image
reference_image = atlas.reference
print(reference_image.shape)
# (528, 320, 456)
# annotation image
annotation_image = atlas.annotation
print(annotation_image.shape)
# (528, 320, 456)
# a hemispheres image (value 1 in left hemisphere, 2 in right) can be generated
hemispheres_image = atlas.hemispheres
print(hemispheres_image.shape)
# (528, 320, 456)
There are multiple ways to work with individual brain regions. To see a dataframe of each brain region, with it's unique ID, acronym and full name, use atlas.lookup_df
:
atlas.lookup_df.head(8)
# acronym id name
# 0 root 997 root
# 1 grey 8 Basic cell groups and regions
# 2 CH 567 Cerebrum
# 3 CTX 688 Cerebral cortex
# 4 CTXpl 695 Cortical plate
# 5 Isocortex 315 Isocortex
# 6 FRP 184 Frontal pole, cerebral cortex
# 7 FRP1 68 Frontal pole, layer 1
Each brain region can also be access by the acronym, e.g. for primary visual cortex (VISp):
from pprint import pprint
VISp = atlas.structures["VISp"]
pprint(VISp)
# {'acronym': 'VISp',
# 'id': 385,
# 'mesh': None,
# 'mesh_filename': PosixPath('/home/user/.brainglobe/allen_mouse_25um_v0.3/meshes/385.obj'),
# 'name': 'Primary visual area',
# 'rgb_triplet': [8, 133, 140],
# 'structure_id_path': [997, 8, 567, 688, 695, 315, 669, 385]}
Working with both image coordinates and cartesian coordinates in the same space can be confusing!
In brainglobe-atlasapi
, the origin is always assumed to be in the upper left corner of the image (sectioning along the first dimension), the "ij" convention.
This means that when plotting meshes and points using cartesian systems, you might encounter confusing behaviors coming from the fact that in cartesian plots one axis is inverted with respect to ij coordinates (vertical axis increases going up, image row indexes increase going down).
To make things as consistent as possible, in brainglobe-atlasapi
the 0 of the meshes coordinates is assumed to coincide with the 0 index of the images stack, and meshes coordinates increase following the direction stack indexes increase.
To deal with transformations between your data space and brainglobe-atlasapi
, you might find the brainglobe-space package helpful.
Contributors to bg-atlaspi are absolutely encouraged, whether you want to fix bugs, add/request new features or simply ask questions.
If you would like to contribute to brainglobe-atlasapi
(or any of the downstream tools like brainrender etc.) please get in touch by opening a new issue or pull request on GitHub.
Please also see the developers guide.
Someone might have already asked a question you might have, so if you're not sure where to start, check out the issues (and the issues of the other repositories).
If you find the BrainGlobe Atlas API useful, please cite the paper in your work:
Claudi, F., Petrucco, L., Tyson, A. L., Branco, T., Margrie, T. W. and Portugues, R. (2020). BrainGlobe Atlas API: a common interface for neuroanatomical atlases. Journal of Open Source Software, 5(54), 2668, https://doi.org/10.21105/joss.02668
Don't forget to cite the developers of the atlas that you used!
For full instructions to add a new BrainGlobe atlas, please see here.
The brainglobe_atlasapi.atlas_generation
submodule contains code for the generation of cleaned-up data, for the main brainglobe_atlasapi
module.
This code was previously the bg-atlasgen
module.
- Fork this repo
- Clone your repo
- Run
git clone https://github.com/brainglobe/brainglobe-atlasapi
- Install an editable version of the package; by running
pip install -e .
within the cloned directory - Create a script to package your atlas, and place into
brainglobe_atlasapi/atlas_generation/atlas_scripts
. Please see other scripts for examples.
Your script should contain everything required to run. The raw data should be hosted on a publicly accessible repository so that anyone can run the script to recreate the atlas.
If you need to add any dependencies, please add them as an extra in the pyproject.toml
file, e.g.:
[project.optional-dependencies]
allenmouse = ["allensdk"]
newatlas = ["dependency_1", "dependency_2"]