Note
Go to the end to download the full example code
8.4. Computer Parcel Aggregation.#
This example uses the ParcelAggregation
marker to compute the mean of each
parcel using the Schaefer parcellations (100 rois, 7 Yeo networks) for both 3D
and 4D NIfTI.
Authors: Federico Raimondo, Synchon Mandal
License: BSD 3 clause
from junifer.testing.datagrabbers import (
OasisVBMTestingDataGrabber,
SPMAuditoryTestingDataGrabber,
)
from junifer.datareader import DefaultDataReader
from junifer.markers import ParcelAggregation
from junifer.utils import configure_logging
Set the logging level to info to see extra information
configure_logging(level="INFO")
2024-04-30 13:12:22,446 - JUNIFER - INFO - ===== Lib Versions =====
2024-04-30 13:12:22,446 - JUNIFER - INFO - numpy: 1.26.4
2024-04-30 13:12:22,446 - JUNIFER - INFO - scipy: 1.11.4
2024-04-30 13:12:22,446 - JUNIFER - INFO - pandas: 2.1.4
2024-04-30 13:12:22,446 - JUNIFER - INFO - nilearn: 0.10.2
2024-04-30 13:12:22,446 - JUNIFER - INFO - nibabel: 5.2.1
2024-04-30 13:12:22,446 - JUNIFER - INFO - junifer: 0.0.5.dev48
2024-04-30 13:12:22,446 - JUNIFER - INFO - ========================
Perform parcel aggregation on VBM GM data (3D) from OASIS dataset
with OasisVBMTestingDataGrabber() as dg:
# Get the first element
element = dg.get_elements()[0]
# Read the element
element_data = DefaultDataReader().fit_transform(dg[element])
# Initialize marker
marker = ParcelAggregation(parcellation="Schaefer100x7", method="mean")
# Compute feature
feature = marker.fit_transform(element_data)
# Print the output
print(feature.keys())
print(feature["VBM_GM"]["data"].shape) # Shape is (1 x parcels)
/opt/hostedtoolcache/Python/3.11.9/x64/lib/python3.11/site-packages/nilearn/datasets/struct.py:852: UserWarning: `legacy_format` will default to `False` in release 0.11. Dataset fetchers will then return pandas dataframes by default instead of recarrays.
warnings.warn(_LEGACY_FORMAT_MSG)
2024-04-30 13:12:22,451 - JUNIFER - INFO - Getting element sub-01
2024-04-30 13:12:22,452 - JUNIFER - INFO - Reading VBM_GM from /home/runner/nilearn_data/oasis1/OAS1_0001_MR1/mwrc1OAS1_0001_MR1_mpr_anon_fslswapdim_bet.nii.gz
2024-04-30 13:12:22,452 - JUNIFER - INFO - VBM_GM is type NIFTI
2024-04-30 13:12:22,452 - JUNIFER - INFO - Computing VBM_GM
2024-04-30 13:12:22,453 - JUNIFER - INFO - Fetching one of Schaefer parcellations.
2024-04-30 13:12:22,453 - JUNIFER - INFO - Parcellation parameters:
2024-04-30 13:12:22,453 - JUNIFER - INFO - resolution: 2.0
2024-04-30 13:12:22,453 - JUNIFER - INFO - n_rois: 100
2024-04-30 13:12:22,453 - JUNIFER - INFO - yeo_networks: 7
2024-04-30 13:12:22,453 - JUNIFER - INFO - At least one of the parcellation files are missing. Fetching using nilearn.
Downloading data from https://raw.githubusercontent.com/ThomasYeoLab/CBIG/v0.14.3-Update_Yeo2011_Schaefer2018_labelname/stable_projects/brain_parcellation/Schaefer2018_LocalGlobal/Parcellations/MNI/Schaefer2018_100Parcels_7Networks_order.txt ...
...done. (0 seconds, 0 min)
Downloading data from https://raw.githubusercontent.com/ThomasYeoLab/CBIG/v0.14.3-Update_Yeo2011_Schaefer2018_labelname/stable_projects/brain_parcellation/Schaefer2018_LocalGlobal/Parcellations/MNI/Schaefer2018_100Parcels_7Networks_order_FSLMNI152_2mm.nii.gz ...
...done. (0 seconds, 0 min)
2024-04-30 13:12:22,749 - JUNIFER - INFO - Loading parcellation /home/runner/junifer/data/parcellations/schaefer_2018/Schaefer2018_100Parcels_7Networks_order_FSLMNI152_2mm.nii.gz
2024-04-30 13:12:22,779 - JUNIFER - INFO - Found existing xfm file for MNI152NLin6Asym to MNI152Lin at /home/runner/junifer/data/xfms/MNI152NLin6Asym_to_MNI152Lin/MNI152NLin6Asym_to_MNI152Lin_Composite.h5
2024-04-30 13:12:23,831 - JUNIFER - INFO - Downloading template MNI152Lin in resolution 2
2024-04-30 13:12:23,928 - JUNIFER - INFO - antsApplyTransforms command to be executed:
antsApplyTransforms -d 3 -e 3 -n 'GenericLabel[NearestNeighbor]' -i /tmp/tmp1q1papv3/parcellationsfgqbi3lh/Schaefer100x7.nii.gz -r /tmp/tmp1q1papv3/parcellationsfgqbi3lh/MNI152Lin_T1w_2.0.nii.gz -t /home/runner/junifer/data/xfms/MNI152NLin6Asym_to_MNI152Lin/MNI152NLin6Asym_to_MNI152Lin_Composite.h5 -o /tmp/tmpt4t2mbgo/parcellations93d7w814/Schaefer100x7_warped_from_MNI152NLin6Asym_to_MNI152Lin.nii.gz
2024-04-30 13:12:25,322 - JUNIFER - INFO - antsApplyTransforms command succeeded with the following output:
b''
2024-04-30 13:12:26,225 - JUNIFER - INFO - No storage specified, returning dictionary
dict_keys(['VBM_GM'])
(1, 100)
Perform parcel aggregation on BOLD data (4D) from SPM Auditory dataset
with SPMAuditoryTestingDataGrabber() as dg:
# Get the first element
element = dg.get_elements()[0]
# Read the element
element_data = DefaultDataReader().fit_transform(dg[element])
# Initialize marker
marker = ParcelAggregation(
parcellation="Schaefer100x7", method="mean", on="BOLD"
)
# Compute feature
feature = marker.fit_transform(element_data)
# Print the output
print(feature.keys())
print(feature["BOLD"]["data"].shape) # Shape is (timepoints x parcels)
2024-04-30 13:12:26,226 - JUNIFER - INFO - Getting element sub001
Dataset created in /home/runner/nilearn_data/spm_auditory
Data absent, downloading...
Downloading data from https://www.fil.ion.ucl.ac.uk/spm/download/data/MoAEpilot/MoAEpilot.zip ...
Downloaded 3768320 of 34212021 bytes (11.0%, 8.3s remaining)
Downloaded 8814592 of 34212021 bytes (25.8%, 5.9s remaining)
Downloaded 14106624 of 34212021 bytes (41.2%, 4.4s remaining)
Downloaded 19357696 of 34212021 bytes (56.6%, 3.1s remaining)
Downloaded 24666112 of 34212021 bytes (72.1%, 2.0s remaining)
Downloaded 30072832 of 34212021 bytes (87.9%, 0.8s remaining) ...done. (8 seconds, 0 min)
Extracting data from /home/runner/nilearn_data/spm_auditory/sub001/MoAEpilot.zip..... done.
2024-04-30 13:12:37,897 - JUNIFER - INFO - Reading BOLD from /tmp/tmpajv_922o/sub001_bold.nii.gz
2024-04-30 13:12:37,897 - JUNIFER - INFO - BOLD is type NIFTI
2024-04-30 13:12:37,898 - JUNIFER - INFO - Reading T1w from /tmp/tmpajv_922o/sub001_T1w.nii.gz
2024-04-30 13:12:37,898 - JUNIFER - INFO - T1w is type NIFTI
2024-04-30 13:12:37,899 - JUNIFER - INFO - Computing BOLD
2024-04-30 13:12:37,899 - JUNIFER - INFO - Fetching one of Schaefer parcellations.
2024-04-30 13:12:37,899 - JUNIFER - INFO - Parcellation parameters:
2024-04-30 13:12:37,899 - JUNIFER - INFO - resolution: 3.0
2024-04-30 13:12:37,900 - JUNIFER - INFO - n_rois: 100
2024-04-30 13:12:37,900 - JUNIFER - INFO - yeo_networks: 7
2024-04-30 13:12:37,901 - JUNIFER - INFO - Loading parcellation /home/runner/junifer/data/parcellations/schaefer_2018/Schaefer2018_100Parcels_7Networks_order_FSLMNI152_2mm.nii.gz
2024-04-30 13:12:37,930 - JUNIFER - INFO - Found existing xfm file for MNI152NLin6Asym to MNI152Lin at /home/runner/junifer/data/xfms/MNI152NLin6Asym_to_MNI152Lin/MNI152NLin6Asym_to_MNI152Lin_Composite.h5
2024-04-30 13:12:38,972 - JUNIFER - INFO - Downloading template MNI152Lin in resolution 2
2024-04-30 13:12:39,069 - JUNIFER - INFO - antsApplyTransforms command to be executed:
antsApplyTransforms -d 3 -e 3 -n 'GenericLabel[NearestNeighbor]' -i /tmp/tmp1q1papv3/parcellations0z9zo087/Schaefer100x7.nii.gz -r /tmp/tmp1q1papv3/parcellations0z9zo087/MNI152Lin_T1w_3.0.nii.gz -t /home/runner/junifer/data/xfms/MNI152NLin6Asym_to_MNI152Lin/MNI152NLin6Asym_to_MNI152Lin_Composite.h5 -o /tmp/tmpt4t2mbgo/parcellationsy1lc0hy7/Schaefer100x7_warped_from_MNI152NLin6Asym_to_MNI152Lin.nii.gz
2024-04-30 13:12:40,466 - JUNIFER - INFO - antsApplyTransforms command succeeded with the following output:
b''
2024-04-30 13:12:43,131 - JUNIFER - INFO - No storage specified, returning dictionary
dict_keys(['BOLD'])
(96, 100)
Total running time of the script: (0 minutes 20.690 seconds)