3810

A joint-community effort to standardize quantitative MRI data: Updates from the BIDS extension proposal
Agah Karakuzu1,2, Gilles Hollander3,4, Stefan Appelhof5, Tibor Auer6, Mathieu Boudreau1,2, Franklin Feingold7, Ali R. Khan8, Alberto Lazari9, Christophe Phillips10, Nikola Stikov1,2, and Kirstie Whitaker11,12
1NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada, 2Montreal Heart Institute, Montreal, QC, Canada, 3Laboratory for Social and Neural Systems Research (SNS Lab), Department of Economics, University of Zurich, Zürich, Switzerland, 4Spinoza Centre for Neuroimaging, Amsterdam, Netherlands, 5Center for Adaptive Rationality, Max Planck Institute for Human Development, Berlin, Germany, 6University of Surrey, Guildford, United Kingdom, 7Stanford University, Stanford, CA, United States, 8Department of Medical Biophysics, Robarts Research Institute, University of Western Ontario, London, ON, Canada, 9Wellcome Centre for Integrative Neuroimaging, FMRIB, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom, 10GIGA Institute, University of Liège, Liège, Belgium, 11Department of Psychiatry, University of Cambridge, Cambridge, United Kingdom, 12The Alan Turing Institute, London, United Kingdom

Synopsis

The Brain Imaging Data Structure (BIDS) aims to develop a standard for organizing and describing neuroimaging data (https://bids.neuroimaging.io/). It has rapidly gained traction across neuroimaging disciplines through community-driven development of BIDS extension proposals (BEPs; http://bit.ly/bids_bep). Here were present such an extension proposal, that captures a wide range of structural MRI contrasts and parametric mapping protocols that are of interest to the broader ISMRM community (https:/bit.ly/bep001). This abstract reports the current state of the proposal and illustrates the developments made since the 2018 ISMRM virtual meeting Bringing BIDS closer to quantitative MRI.

INTRODUCTION

As neuroimaging research is becoming increasingly data-intensive and with the surgeance of dataset sharing within the last decade, proper data standardization and annotation is of paramount importance1. The Brain Imaging Data Structure2,3 (BIDS) aims to develop a standard for organizing and describing neuroimaging data.

Different vendors and labs have different conventions for organizing neuroimaging data, which makes it difficult to standardize processing (Figure 1). BIDS standardizes this organisation, which eases the development of open-source turn-key BIDS apps and sharing of public datasets. Its applicability has steadily been increased by the community-driven development of BIDS extension proposals (BEP).

Over the last two years, we have worked to bring MRI physicists and neuroimaging researchers together and collaborate on a BIDS specification to facilitate the standardization and use of quantitative MRI (qMRI) methods by neuroscientists. In this work, we introduce our specification, share example datasets, and illustrate how multiple BIDS apps4 can easily work on the same BIDS-compliant qMRI data.

METHODS

Community-driven open development

BEP001 contributors have been holding monthly video-conferences over the last two years, making consensus decisions for improvement and sharing action items. The standard operational procedure to advance our proposal is summarized in Figure 2.

qMRI within the BIDS framework

BIDS was originally developed by functional neuroimaging researchers who laid out a directory structure covering anatomical, functional, diffusion-weighted and field map images. BEP001 focuses on structural (anatomical) acquisitions that include multiple contrasts and sensitivity maps. Therefore has placed all qMRI data in the anat and fmap folders.

In the BIDS specification the suffix of the file name indicates the type of data captured in the image, for example _T1w and _bold. In this extension, we added 3 categories of suffixes to annotate the use-case of more advanced anatomical imaging data: i) conventional images, ii) quantitative maps and iii) grouped scan collections (Figure 3). Grouped scan collections identify protocols where multiple images are acquired with largely identical parameters, except for a few key parameters such as inversion time, echo time or flip angle.

We introduce the indexable_metadata-<index> entity to capture such key parameters. For example a multi-echo gradient echo-recalled sequence will include images with the key pairs _echo-1 and _echo-2 in their respective filenames. If the parameters that vary between the constituents of a grouped scan are categorical rather than numerical, the _acq-<label> entity is used. BEP001 suggests a list of labels to be used with the _acq key tag for the sake of creating a coherent naming entity for grouped scans. Currently available _acq labels are MTon, MToff, and T1w, creating sensible filenames for e.g., magnetization transfer saturation6 (MTsat) and multi-parameter mapping7 .

BIDS makes a strong distinction between "raw" images that are directly imported from the scanner and "derived" images that have been processed in some way. However, it is important to emphasize that "raw" in the BIDS context does not refer to MR k-space data before reconstruction, which means that BIDS is compatible with the ISMRM raw data format8. To maintain consistency and ensure compatibility between the two standards, we propose that quantitative parameter maps that are estimated with a model are stored in the /derivatives folder, whereas unprocessed are stored in the raw /anat folder. To aid BIDS processing pipelines that only take raw data as an input, researchers can use file links to make parameter maps also available in the raw data folder.

Metadata for grouped scan collections

In addition to a file naming system, BIDS also prescribed how to store relevant metadata in a "sidecar" json-file. This is especially relevant for qMRI applications, where key imaging parameters are essential for further processing. Our proposal standardizes which parameters are required and optional for specific grouped acquisition protocols (e.g., MP2RAGE, VFA; Figure 4a). It also standardizes optional parameters that are relevant for specific qMRI applications (Figure 4b).

RESULTS

We created an open data repository9 and uploaded example qMRI datasets adhering to the BEP001 proposal. To demonstrate a use case, we used two open-source software packages, pymp2rage10 and qMRLab11,12, for T1 mapping using MP2RAGE data from our open repository (Figure 5). The analyses are version controlled on GitHub and fully reproducible, with no installation, at https://bit.ly/mp2rage.

So far, BEP001 has merged 22 pull requests and hosted 48 issues with the contributions of more than 10 developers around the globe. Currently it lists 5 descriptions for conventional contrasts, 14 for qMRI maps and 9 for grouped scan collections along with their metadata priority level definitions. A sizable multi-center public database for spinal cord MRI13, which attracted tens of thousands of viewers in a few months, has already adopted conventions suggested by BEP001. Members of the ISMRM community are invited to join the development process to make their methods easy to implement for researchers in the cognitive and clinical neurosciences.

CONCLUSION

In this work, we introduced the core components of BEP001, shared real-world example datasets and demonstrated how applications written in different languages can easily work on the same BIDS compliant qMRI data. We hope that these efforts will not only better organize our data, but will also facilitate multi-site collaborative work and encourage neuroscientists to adopt advanced MR techniques.

Acknowledgements

The authors A.K. and G.H. contributed equally to this work.

References

[1] Poldrack, R., et al. “Scanning the horizon: towards transparent and reproducible neuroimaging research.” Nature Reviews Neuroscience 18 (2017): 115-126.

[2] Gorgolewski, K., et al. "The brain imaging data structure, a format for organizing and describing outputs of neuroimaging experiments." Scientific Data 3 (2016): 160044.

[3] https://bids.neuroimaging.io/, Accessed November 5, 2019.

[4] Gorgolewski, K., et al. "BIDS apps: Improving ease of use, accessibility, and reproducibility of neuroimaging data analysis methods." PLoS Comp. Biol. 13.3 (2017): e1005209.

[5] https://github.com/bids-standard/bep001, Accessed November 5, 2019.

[6] Helms, G., et al. "High‐resolution maps of magnetization transfer with inherent correction for RF inhomogeneity and T1 relaxation obtained from 3D FLASH MRI." Magn. Reson. Med. 60 (2008): 1396-1407.

[7] Weiskopf, N., et al. "Quantitative multi-parameter mapping of R1, PD*, MT, and R2* at 3T: a multi-center validation." Frontiers in Neuroscience 7 (2013): 95.

[8] Inati, S., et al. “ISMRM Raw data format: A proposed standard for MRI raw datasets”, Magn. Reson. Med. 7 (2017): 411-421.

[9] Karakuzu, A., et al. (2019). BEP[001]. https://doi.org/10.17605/OSF.IO/K4BS5, Accessed October 28, 2019.

[10] https://github.com/Gilles86/pymp2rage, Accessed November 5, 2019.

[11] Duval, T., et al. “Quantitative MRI made easy with qMRLab”, International Society for Magnetic Resonance in Medicine Meeting (2018): Abstract 2288.

[12] https://qmrlab.org, Accessed November 5, 2019.

[13] Cohen-Adad, J., et al. “Spinal Cord MRI Public Database (multi-subjects)”, https://openneuro.org/datasets/ds001919/versions/1.0.4, DOI: 10.18112/openneuro.ds001919.v1.0.4, Accessed October 26, 2019.


Figures

Figure 1. Despite differences across vendors, the DICOM standard provides a common ground for reconstructed images. However, organizing, annotating and describing these images can easily create idiosyncrasies, challenging interoperability and decreasing efficiency of processing neuroimaging datasets (orange). BIDS resolves this problem through community-driven development of consensus data specifications, establishing a cohesive landscape of neuroimaging data organization (blue).

Figure 2. Transparency and accessibility are central to the standard operational procedure for improving BEP001. Outcomes from the monthly meetings are transferred to a central GitHub repository5, opened for more elaborate public discussions via issues and merged into the proposal through peer-reviewed pull requests. BEP001 is inclusive to all communities who would like to contribute to the proposal or keep themselves up-to-date with the latest developments.

Figure 3. To minimise deviations from the current BIDS specification, qMRI data is stored under the anat folder. Therefore, anatomical imaging data is divided into three subgroups to distinguish different use cases (left panel). New naming entities are created (indexable_metadata) or existing ones are modified (acq, suffix) based on this classification, so that quantitative maps can be differentiated from conventional images, and grouped scan collections are properly banded together to reveal their role (right panel).

Figure 4. a) FlipAngle is defined as a varying REQUIRED field, as it changes between images grouped by the VFA suffix. RepetitionTimeExcitation and SequenceType parameters are REQUIRED for calculations, but the same across the board. The PhaseIncrement is listed under the OPTIONAL columns, indicating that it can be used to derive different qMRI applications. b) The relationship between DESPOT1, DESPOT2 and DESPOT2FM applications and the VFA grouping suffix with respect to the metadata fields.

Figure 5. Interoperable data is one of the key advantages brought by standardization. To demonstrate this, we ran two open-source software (left panel) on a BIDS-formatted MP2RAGE dataset (middle panel) to perform T1 mapping. Despite being written in different languages, outputs were in agreement and followed the exact same data structure (right panel). You can repeat this analysis on your browser with one click: https://bit.ly/mp2rage.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
3810