BRIEF RESEARCH REPORT article
Beware (Surprisingly Common) Left-Right Flips in Your MRI Data: An Efficient and Robust Method to Check MRI Dataset Consistency Using AFNI
- 1Scientific and Statistical Computing Core, NIMH/NIH/DHHS, Bethesda, MD, United States
- 2Rotman Research Institute at Baycrest, Toronto, ON, Canada
- 3Department of Psychology, University of Toronto, Toronto, ON, Canada
Knowing the difference between left and right is generally assumed throughout the brain MRI research community. However, we note widespread occurrences of left-right orientation errors in MRI open database repositories where volumes have contained systematic left-right flips between subject EPIs and anatomicals, due to having incorrect or missing file header information. Here we present a simple method in AFNI for determining the consistency of left and right within a pair of acquired volumes for a particular subject; the presence of EPI-anatomical inconsistency, for example, is a sign that dataset header information likely requires correction. The method contains both a quantitative evaluation as well as a visualizable verification. We test the functionality using publicly available datasets. Left-right flipping is not immediately obvious in most cases, so we also present visualization methods for looking at this problem (and other potential problems), using examples from both FMRI and DTI datasets.
As part of the NIFTI dataset standard (Cox et al., 2004), orientation and location information were included in the file headers, to be able to reduce uncertainty in interpretation across software and systems. However, it is possible for mistakes to occur while recording information at the scanner, while interpreting DICOM fields, when converting to NIFTI or another format, or during a subsequent processing step. This can lead to seriously erroneous results (or to catastrophic results in clinical surgical cases).
While some header mistakes can be easily spotted visually (e.g., having incorrect voxel dimensions recorded), there are more subtle changes, particularly involving the positioning of the data within the acquisition field of view (FOV). The header contains the “orientation” parameters for how the data matrix is stored on disk, so that the rows and slices appear in the correct locations and view planes; an incorrect value can lead to axes that are switched (e.g., axial slices are interpreted as coronal ones) or flipped (e.g., the anterior part of the brain is labeled as posterior). Some of these problems can be recognized instantaneously, such as a flip in the anterior-posterior or inferior-superior axes in a whole brain acquisition; but others are much more subtle, such as left-right flips, due to the large-scale structural symmetry of the brain. Consider Figure 1, which shows an axial slice of a T1w anatomical volume with both the original version of that subject’s EPI and a left-right flipped version (created by altering the orientation field of the file header). The correct orientation is often not visually obvious even when looking for this problem.
Figure 1. Axial slice of a T1w anatomical volume (left) with two versions of the same subject’s EPI volume (right): one is the original, and one is left-right flipped (through changing the orientation value stored in the file header). It is not immediately apparent which EPI image matches correctly with the anatomical, which is a problem for trusting results of analyzes.
Here, we propose a simple method to detect relative left-right flips within pairs of MRI datasets (though it may also be more broadly applicable), available within the AFNI software package (Cox, 1996). The method applies to datasets with multiple acquisition types, such as EPI and T1w datasets acquired together for FMRI protocols, or DWI and T2w datasets for DTI protocols. How do we know that this problem actually occurs in practice? By using our approach, we have found systematic left-right flips in datasets submitted to each large public repository in which we have looked: the Functional Connectome Project (FCP) (Biswal et al., 2010), OpenFMRI (Poldrack et al., 2013), and ABIDE (Di Martino et al., 2014). These problems were verified by each notified consortium and subsequently fixed1. We note this not to criticize these projects (indeed, the openness of these databases has made it possible to address these issues), but instead to point out how important it is for researchers to be able to check the basic properties of their data. This problem also exists in clinical imaging, with several studies reporting on rates of “laterality errors” in radiology (e.g., Bernstein, 2003; Sangwaiya et al., 2009; Landau et al., 2015; Digumarthy et al., 2018). If this fundamental problem can occur within large, public dataset collections that have been looked at by many people, as well as in clinical settings that affect patient outcome directly, then surely this is an issue that all researchers and neuroimagers should be aware of.
We present a simple diagnosis for this data issue. This efficient and robust check is available within the AFNI software package, along with several ways to visually verify the quantitative results. The visualization methods are shown using AFNI, but some of these methods are likely available in other software packages. In the Discussion, we note some of the potential causes of the flip issue.
The primary method introduced here is to check for left-right flipping through volumetric alignment: comparing the alignment cost function value between a pair of original datasets with the cost function of alignment when one dataset has been intentionally left-right flipped. This relies on having an appropriate cost function for the alignment. For instance, in FMRI studies, one generally wants to align a subject’s T1w anatomical with an EPI volume, which has a very different contrast pattern; T1w volumes in adult human datasets show decreasing brightness from white matter (WM) to gray matter (GM) and then to cerebrospinal fluid (CSF); in an adult human T2w EPI, the intensities are reordered relative to the T1w volumes, with CSF the brightest, then GM, and then WM. The reversal of the CSF relative intensities (i.e., the CSF is bright in EPI data yet dark in T1w datasets) is a key property that can be used for the advantage of alignment, in conjunction with aligning sulcal and gyral features. For these cases, AFNI uses the “local Pearson correlation (LPC)” cost function (Saad et al., 2009) or a variant “LPC+ZZ,” which have been shown to be robust in such cases. The LPC cost function is the negative of the sum of correlations computed over local regions (“patches”). Alignment then proceeds by optimizing for the minimal cost (negative correlation). The related “LPC+ZZ” cost function operates similarly but uses a combination of cost functions for its initial estimate and then finalizes the alignment parameters with the standard LPC cost2; this refinement has shown improvement in stability when the volumes being aligned have greater initial differences (e.g., large relative rotation). For aligning brains with similar tissue contrast (e.g., T2w anatomical with the b = 0 s/mm2 volume in DWI datasets), AFNI typically uses the similarly robust “local Pearson absolute (LPA)” or related “LPA+ZZ” cost function, both using the absolute value of the local Pearson correlation.
The inputs for the check are simply the EPI and either anatomical volume; because the check is for relative flips, neither dataset is tagged as “correct” a priori. AFNI’s align_epi_anat.py program performs linear affine (12 degrees of freedom, by default) alignment between the volumes3. There is no requirement for the datasets to be in particular coordinate or orientation systems; typically, they are in native/acquired coordinates. One simply provides an option flag “-check_flip,” and the program will generate cost function results and aligned volumetric datasets for the two cases: for the original volumes, and with one volume (e.g., the anatomical) left-right flipped. If the alignment is improved (lower cost function metric) for the intentionally left-right flipped data, then the flipped dataset is likely better, and an error is probably present in one of the original volume headers (NB: this method detects inconsistencies between volumes in a dataset, with further follow-up needed to determine which file is incorrect; see the Discussion for more on this). The cost function results and an ensuing recommendation about whether a dataset needs to be flipped (“NO_FLIP” or “DO_FLIP”), are saved in a text file. The aligned volume from each case is also output, for visual comparison and verification.
The left-right flip check can also be performed conveniently as part of the FMRI processing pipeline constructed by AFNI’s afni_proc.py. The option “-align_opts_aea -check_flip” can be added so that the same check is performed (via align_epi_anat.py) during the “align” block of processing. The quantitative results of the flip check are automatically parsed and presented to the researcher—along with image snapshots for visual verification—as part of the automatically-generated quality control (QC) HTML output, which afni_proc.py also creates.
To test the efficacy of this left-right flip test on a range of data, we downloaded publicly available FMRI datasets, both human and non-human (macaque). For variety, datasets were downloaded from a range of locations and projects, with a range of data quality and acquisition parameters; data was only used from subjects who had both an anatomical and FMRI volume in the same session. Human datasets were: OpenFMRI (ds0000003, ds000114, ds000172), FCP (Beijing-Zang, part 1; Cambridge-Buckner, part 1; New York, a/ADHD; Taipei, a) and ABIDE-II (UPSM_Long). Macaque datasets were downloaded from the PRIMate Data Exchange (PRIME-DE; Milham et al., 2018), from the following sites: Institute of Neuroscience (ION), Shanghai; Netherlands Institute of Neuroscience (NIN); Newcastle University; Stem Cell and Brain Research Institute (SBRI); and University of Minnesota (UMinn).
For each study, a subject was randomly assigned to be either “flipped” or “unflipped,” with the anatomical volume left-right flipped for the former group. Each subject’s data were processed using a brief afni_proc.py command with alignment blocks, including the left-right flip check option (see Supplementary Figure S1 for the full command). The flip results were tallied, and visual verification using the QC HTML output was performed.
An example4 of output from running the left-right flip check as part of FMRI processing with afni_proc.py is shown in Figure 2. The figure shows the relevant part of afni_proc.py’s QC HTML doc, contained within the “warns” block that reports on potential warnings during the processing. The quantitative results are reported (they are also stored in a text or JSON file, for any later use), and in the QC document they are also parsed for immediate identification; also, axial slice montages of each alignment’s result are shown for verification, in the form of the edges of the EPI volume displayed over the original or flipped anatomical. Here, indeed, matching of sulcal and gyral patterns confirm the quantitative flip-check results.
Figure 2. An example of left-right flip checking during FMRI processing with afni_proc.py, in the automatically generated QC HTML doc. This information is displayed in the “warns” block of the HTML page. The name of the cost function, its values for both cases, and a resulting “guess” evaluation by the program are reported at the top: the “none” in green denotes that NO_FLIP was the result here; a warning in red would be reported there and at the top of the QC page, instead, if DO_FLIP were the result (i.e., if evidence of flipping were found). For visual verification, montages of the results of each alignment are shown as edges of the EPI volume overlaid on each case: original anatomical at top, and flipped version at bottom. In particular, the structures in the superior slices of the cortex provide clear verification of the “flip guess” alignment results.
The results from running the left-right flip check on several publicly available human datasets are summarized in Table 1. Out of 178 subjects analyzed, 100% received the correct left-right flip recommendation: 90/90 subjects in the “flipped” group received DO_FLIP, and 88/88 in the “unflipped” group received NO_FLIP. Visual checks of each subject’s automatically-produced QC HTML (see example in Figure 2) verified the results.
In the case of macaque datasets, 26 subjects contained both an anatomical and a functional MRI from the same session. Compared to the human datasets, several EPIs had lower tissue contrast, brightness inhomogeneity across the volume and greater distortion (e.g., EPI distortion); additionally, the FOV of EPI and anatomical datasets were shifted in many cases, making alignment much more difficult. In total, either anatomical skullstripping or EPI-anatomical alignment failed in seven subjects. From the remaining 19 subjects, 7/9 “flipped” and 10/10 “unflipped” were correctly categorized (89% correct total); see Table 2. However, it should be noted that the cost function values in seven of these cases were extremely close and difficult to differentiate. This might be expected due to the aforementioned contrast, brightness and distortion effects, as well as to the inherently high left-right symmetry of anatomy in macaques (much higher than in humans); this symmetry also renders visual verification of the flip check more difficult.
In addition to the QC images, one can view the outputs of align_epi_anat.py (run either directly or via afni_proc.py) interactively in the AFNI graphical user interface (GUI), and there are several features to aid in the visual comparison of overlaps. While many viewer software packages allow for control of transparency to show one dataset over another, we recommend additional methods for interactive comparisons. Several of these visualization methods are shown in Supplementary Figure S1. Supplementary Figure S1 displays several of these for cases of EPI-T1w overlap (differing tissue contrasts) and T2w-DWI comparison (similar tissue contrasts).
Left-right flipping is an unfortunately too common occurrence in MRI analysis. As noted in the Introduction, several of the largest public repositories of publicly available data have had systematically flipped datasets, even after an initial curation process by some of the field’s most experienced researchers. This problem can occur in any sized dataset collection, however. In all cases it is incumbent upon the researchers using the data to guard against such problems. We have presented here a simple method for doing so in AFNI, using alignment. This method can be integrated seamlessly into a processing stream (via afni_proc.py) or checked separately (via align_epi_anat.py). This approach is efficient, as it takes of order 1 min to perform the check. The combination of both quantitative evaluation (using LPC/LPA cost function results) and qualitative verification (using overlaid images of structural features) has performed robustly.
We note that the percent correct (100%) in the tests of human public datasets here was higher than expected. Due to the varying quality of tissue contrast, noise and artifacts that can occur in MRI data, we would generally expect some false positives and negatives in the left-right flip checks (which would be verified with visual checks). We also note that no “hidden” flips were found in the data – i.e., datasets which turned out to have mismatched header information. Based on previous use of the check flip tool introduced here, we had already alerted several repositories, and they have fixed many of the problems. We hope that such curation can continue prospectively in future uploads, and we would recommend the tools presented here for doing so. Such tests are also similarly useful for any center or institute acquiring data, even before it might be analyzed or made public.
It should be noted that this alignment-based method can identify the presence of a relative left-right flip between volumes in a dataset. However, it cannot identify an absolute left-right flip – that is, it cannot determine which of the volumes has incorrect header information. That kind of information is more difficult to assess, and likely it must be investigated using the original data: checking scan parameters, DICOM conversion and other processing steps. As a corollary, this method cannot detect if both volumes in the pair are left-right flipped. However, it may be possible to extend this approach to using an asymmetric MRI template for reference (e.g., the ICBM 2009a/b Non-linear Asymmetric MNI templates; Fonov et al., 2009). The degree of variability between subjects and the template would likely reduce the certainty of “absolute flip” detection greatly, particularly across a wide age range or in the case of pathology. Potentially, one could determine left-rightedness based on noted population differences like the petalia, often found in the human brain where the right hemisphere protrudes anteriorly and the left hemisphere protrudes posteriorly (Toga and Thompson, 2003). While it is unlikely that the difference would be suitably reliable on an individual basis, such a biologically-based method could be used as a basis for determining left-right flips over a group of subjects (typically, a group of subjects with the same acquisition and conversion steps will all have the same flip properties). We note here that we have started to look at this approach for determination of the likelihood of an entire group of subjects being flipped and hope to develop a useful implementation of this absolute flip check in the future.
There can be many causes for left-right flipping (and other header problems) in datasets. As noted above, in cases where flipping is detected, a researcher will likely have to backtrack through the provenance and processing of the datasets to find the root cause of the misinformation. For example, DICOM data from scanners can be wrong or ambiguous; the mosaic format of Siemens can have slices stacked in a reversed order, which is in a non-standard part of the file header and hence may not be read properly during conversion. The position of the subject in the scanner must be correctly recorded during acquisition as supine or prone (or “sphinx,” for animal studies). DICOM Conversion PACS systems and custom scripts can also misinterpret information that may have changed in the format, particularly in vendor-specific tags. Older Analyze-7.5 format datasets lack accurate orientation information, and such files may have been used or passed along in studies. NIFTI format datasets have stricter definitions of orientation in their headers, but conflicting sform_matrix and qform_matrix information can still occur, leading to incorrect conversion. Some processing scripts rely on read/write functions in tools that apply either a default orientation or no orientation information, missing consistency checks (e.g., in Matlab, ImageJ, or other software); because the NIFTI formats are not part of these codebases’ native format, maintaining correct orientation information in the header is non-trivial and can easily result in mistakes. Some software packages also make assumptions about coordinate storage order and orientation and will either assume the input is in a specific orientation or that multiple input datasets all match each other. Analytical results will have passed through a multitude of steps of conversion, processing and/or regridding in most FMRI and DTI pipelines; maintaining consistency through all steps is a challenge within a software package, and combining processing across packages, each with their own assumptions (e.g., some packages ignore such header information), can potentially result in errors.
While the validation of the flip check results of both the human and animal data were consistently high, we note some important differences approaching each of these datasets. Firstly, it should be expected that the flip check results would be much more sensitive in humans, due to the greater left-right asymmetry of structure; indeed, the symmetry of macaque brains made it more difficult to determine consistency with a high degree of confidence (i.e., clearly differentiated cost function and clear visual verification). For any group of subjects with high left-right symmetry, this check becomes less reliable. The quality of data (e.g., amount of artifact, relative tissue contrast, consistent coordinates) also affects the reliability of the check; in the present case, the macaque datasets presented a greater challenge due to such issues. While there are any number of species (human, macaque, rodent, etc.) and types of data (EPI-anatomical, longitudinal anatomical, etc.) to compare, the present results suggest that: for datasets with reasonable quality and left-right asymmetry, this left-right flip check should provide a useful consistency check.
The method presented here provides a simple, fast and verifiable method for an automatic determination of potential left-right flipping problems. In our opinion, it should be included in basic processing streams as a standard feature to detect potential problems in the data. For researchers making data public, it would be particularly beneficial to the community to have run this test. In particular, it may be possible to integrate this with BIDS (Brain Imaging Data Structure) (Gorgolewski et al., 2016) organizational framework and include some basic header validation for coordinate system and orientation.
It should be noted that this left-right flip problem was first observed by an AFNI user (author BRB) while reviewing data, and this highlights an important point: visualizing data remains extremely important in neuroimaging. While we can devise new methods to automatically find some problems (like the left-right flipping issue), there will always be another unforeseen problem that requires carefully looking at the data, and in different ways. We suggest methods that include overlay opacity control, layer toggling, vertical and horizontal curtains, layer blending, checkerboard, and edge displays. These kinds of visualization methods help researchers to identify unexpected problems.
The presented method is not foolproof. As noted above, if datasets have poor structural contrast, which may occur fairly often in EPI datasets, for example, then the alignment costs are less reliable. EPI with large flip angles can be a source of the lack of structural information (Gonzalez-Castillo et al., 2013). Also for data where there is partial coverage, similar lack of structure can occur and potentially decrease the stability of this method.
We have found even very basic properties of MRI data like left and right can be confused. Here we presented a simple method to determine consistency among datasets and visualization methods for other unforeseen issues. Even in the era of big data, details still matter – some even more than before, because curating large datasets across multisite studies can pose many new challenges. In the end using a definitive, physical marker while scanning (such as a vitamin E capsule, with a recorded side of placement) is the most robust method to recognize the presence of left-right flipping in a dataset; such a method would also have the benefit of determining flips absolutely. However, to date this practice has not been widely adopted across the neuroimaging community. The presented left-right flip check method in AFNI is simple, efficient and robust. Under the guiding principle of “caveat emptor” when using public data, and good practice when using self-acquired data, we strongly recommend the inclusion of this check in all MRI processing pipelines.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://openneuro.org/datasets/ds000114/versions/1.0.1/download, https://openneuro.org/datasets/ds000172/versions/1.0.1/download, https://openneuro.org/datasets/ds000003/versions/00001/download, https://www.nitrc.org/frs/?group_id=404, https://www.nitrc.org/frs/?group_id=296.
DG: concept and scripting. PT, DG, RR: scripting, testing, analysis, and writing. BB: concept and initial observations. RC: underlying software. RR: scripting and analysis.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The research and writing of the manuscript were supported by the NIMH and NINDS Intramural Research Programs (ZICMH002888) of the NIH (HHS, United States), and by the NIH’s Brain Initiative (1R24MH117467-01). This work utilized the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov).
We thank O. Esteban, G. A. Devenyi and an anonymous reviewer for suggestions to clarify this work.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2020.00018/full#supplementary-material
- ^ Discussion of orientation issues in FCP datasets: https://www.nitrc.org/forum/forum.php?thread_id=1310&forum_id=1243; Note about dataset orientation fixes in OpenFMRI (after email communication): https://openfmri.org/dataset-orientation-issues/ Note about dataset orientation fixes in ABIDE (after email communication): fcon_1000.projects.nitrc.org/indi/abide/updates/ABIDEII-Usernotes_Updates_ Fixes_9_25_16.pdf fcon_1000.projects.nitrc.org/indi/abide/updates/ABIDEII-Usernotes_Updates_ Fixes_3_27_17.pdf
- ^ The initial estimate of LPC+ZZ uses the following formula of cost functions: LPC + 0.4∗(HEL + CRA + OV) + 0.2∗(NMI + MI), where Hellinger metric, normalized mutual information (NMI), additively symmetrized correlation ratio (CRA), overlap (OV), mutual information (MI) and normalized mutual information (NMI); using multiple cost functions here stabilizes the initial matching. The final stages of alignment are still evaluated using only the LPC cost function, which provides good detail matching due to its designed localness.
- ^ Its name reflects its original purpose – aligning a subject’s EPI to their T1w anatomical – but it is actually more broadly applicable to any linear affine alignment (typically, to any volumes belonging to the same subject); other programs are recommended for nonlinear alignment.
- ^ This dataset is publicly available as part of the AFNI Bootcamp demo package (https://afni.nimh.nih.gov/pub/dist/edu/data/CD.tgz), located in the “AFNI_data6/FT_analysis/” directory and processed with afni_proc.py using the accompanying s05.ap.uber script.
Cox, R. W., Ashburner, J., Breman, H., Fissell, K., Haselgrove, C., Holmes, C. J., et al. (2004). “A (sort of) new image data format standard: NIfTI-1,” in 10th Annual Meeting of the Organization for Human Brain Mapping, Budapset, 2004.
Di Martino, A., Yan, C. G., Li, Q., Denio, E., Castellanos, F. X., Alaerts, K., et al. (2014). The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. Mol. Psychiatry 19, 659–667. doi: 10.1038/mp.2013.78
Digumarthy, S. R., Vining, R., Tabari, A., Nandimandalam, S., Otrakji, A., Shepard, J. O., et al. (2018). Process improvement for reducing side discrepancies in radiology reports. Acta Radiol. Open 7:2058460118794727. doi: 10.1177/2058460118794727
Gonzalez-Castillo, J., Duthie, K. N., Saad, Z. S., Chu, C., Bandettini, P. A., and Luh, W. M. (2013). Effects of image contrast on functional MRI image registration. NeuroImage 67, 163–174. doi: 10.1016/j.neuroimage.2012.10.076
Gorgolewski, K. J., Auer, T., Calhoun, V. D., Craddock, R. C., Das, S., Duff, E. P., et al. (2016). The brain imaging data structure, a format for organizing and describing outputs of neuroimaging experiments. Sci. Data 3:160044. doi: 10.1038/sdata.2016.44
Poldrack, R. A., Barch, D. M., Mitchell, J. P., Wager, T. D., Wagner, A. D., Devlin, J. T., et al. (2013). Toward open sharing of task-based fMRI data: the OpenfMRI project. Front. in Neuroinformatics 7:12. doi: 10.3389/fninf.2013.00012
Saad, Z. S., Glen, D. R., Chen, G., Beauchamp, M. S., Desai, R., and Cox, R. W. (2009). A new method for improving functional-to-structural MRI alignment using local Pearson correlation. Neuroimage 44, 839–848. doi: 10.1016/j.neuroimage.2008.09.037
Sangwaiya, M. J., Saini, S., Blake, M. A., Dreyer, K. J., and Kalra, M. K. (2009). Errare humanum est: frequency of laterality errors in radiology reports. Am. J. Roentgenol. 192, W239–W244. doi: 10.2214/AJR.08.1778
Keywords: symmetry, asymmetry, DICOM, Nifti, fMRI, DTI, T1-weighted brain image, EPI (echo planar imaging)
Citation: Glen DR, Taylor PA, Buchsbaum BR, Cox RW and Reynolds RC (2020) Beware (Surprisingly Common) Left-Right Flips in Your MRI Data: An Efficient and Robust Method to Check MRI Dataset Consistency Using AFNI. Front. Neuroinform. 14:18. doi: 10.3389/fninf.2020.00018
Received: 18 October 2019; Accepted: 14 April 2020;
Published: 25 May 2020.
Edited by:Pierre Bellec, Université de Montréal, Canada
Reviewed by:Oscar Esteban, Stanford University, United States
Gabriel Allan Devenyi, Douglas Hospital Research Centre, McGill University, Canada
Copyright © 2020 Glen, Taylor, Buchsbaum, Cox and Reynolds. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel R. Glen, firstname.lastname@example.org