Purpose
To obtain three‐dimensional (3D), quantitative and motion‐robust imaging with magnetic resonance fingerprinting (MRF).
Methods
Our acquisition is based on a 3D spiral projection k‐space ...scheme. We compared different orderings of trajectory interleaves in terms of rigid motion‐correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7‐s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole‐brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k‐space data, and the aligned data were matched with the dictionary to obtain motion‐corrected maps.
Results
A significant improvement on the motion‐affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T1 and T2 estimations after motion correction. In addition, the average motion‐induced quantification bias of 70 ms for T1 and 18 ms for T2 values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.
Conclusion
We established a method that allows correcting 3D rigid motion on a 7‐s timescale during the reconstruction of MRF data using self‐navigators, improving the image quality and the quantification robustness.
Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and ...reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories.
Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions.
After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site.
Our analysis demonstrated a high repeatability (CVs 0.7–1.3% for T1, 2.0–7.8% for T2, 1.4–2.5% for normalized PD) and reproducibility (CVs of 2.0–5.8% for T1, 7.4–10.2% for T2, 5.2–9.2% for normalized PD) in gray and white matter.
Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.
Magnetic resonance imaging (MRI) nowadays plays an important role in the identification of brain underpinnings in a wide range of neuropsychiatric disorders, including Autism Spectrum Disorders ...(ASD). Characterizing the hallmarks in these pathologies is not a straightforward task and machine learning (ML) is certainly one of the most promising tools for addressing complex and non-linear problems. ML algorithms and, in particular, deep neural networks (DNNs), need large datasets in order to be properly trained and thus ensure generalization capabilities on new data. Large datasets can be obtained by collecting images from different centers, thus bringing unavoidable biases in the analysis due to differences in hardware and scanning protocols between different centers. In this work, we dealt with the issue of multicenter MRI data harmonization by comparing two different approaches: the analytical ComBat-GAM procedure, whose effectiveness is already documented in the literature, and an originally developed site-adversarial deep neural network (ad-DNN). The latter aims to perform a classification task while simultaneously searching for site-relevant patterns in order to make predictions free from site-related biases. As a case study, we implemented DNN and ad-DNN classifiers to distinguish subjects with ASD with respect to typical developing controls based on functional connectivity measures derived from data of the multicenter ABIDE collection. The classification performance of the proposed ad-DNN, measured in terms of the area under the ROC curve (AUC), achieved the value of AUC = 0.70±0.03, which is comparable to that obtained by a DNN on data harmonized according to the analytical procedure (AUC = 0.71±0.01). The relevant functional connectivity alterations identified by both procedures showed an agreement between each other and with the patterns of neuroanatomical alterations previously detected in the same cohort of subjects.
Numerous publications and commercial systems are available that deal with automatic detection of pulmonary nodules in thoracic computed tomography scans, but a comparative study where many systems ...are applied to the same data set has not yet been performed. This paper introduces ANODE09 (
http://anode09.isi.uu.nl), a database of 55 scans from a lung cancer screening program and a web-based framework for objective evaluation of nodule detection algorithms. Any team can upload results to facilitate benchmarking. The performance of six algorithms for which results are available are compared; five from academic groups and one commercially available system. A method to combine the output of multiple systems is proposed. Results show a substantial performance difference between algorithms, and demonstrate that combining the output of algorithms leads to marked performance improvements.
Machine Learning (ML) is nowadays an essential tool in the analysis of Magnetic Resonance Imaging (MRI) data, in particular in the identification of brain correlates in neurological and ...neurodevelopmental disorders. ML requires datasets of appropriate size for training, which in neuroimaging are typically obtained collecting data from multiple acquisition centers. However, analyzing large multicentric datasets can introduce bias due to differences between acquisition centers. ComBat harmonization is commonly used to address batch effects, but it can lead to data leakage when the entire dataset is used to estimate model parameters. In this study, structural and functional MRI data from the Autism Brain Imaging Data Exchange (ABIDE) collection were used to classify subjects with Autism Spectrum Disorders (ASD) compared to Typical Developing controls (TD). We compared the classical approach (external harmonization) in which harmonization is performed before train/test split, with an harmonization calculated only on the train set (internal harmonization), and with the dataset with no harmonization. The results showed that harmonization using the whole dataset achieved higher discrimination performance, while non-harmonized data and harmonization using only the train set showed similar results, for both structural and connectivity features. We also showed that the higher performances of the external harmonization are not due to larger size of the sample for the estimation of the model and hence these improved performance with the entire dataset may be ascribed to data leakage. In order to prevent this leakage, it is recommended to define the harmonization model solely using the train set.
No univocal and reliable brain-based biomarkers have been detected to date in Autism Spectrum Disorders (ASD). Neuroimaging studies have consistently revealed alterations in brain structure and ...function of individuals with ASD; however, it remains difficult to ascertain the extent and localization of affected brain networks. In this context, the application of Machine Learning (ML) classification methods to neuroimaging data has the potential to contribute to a better distinction between subjects with ASD and typical development controls (TD). This study is focused on the analysis of resting-state fMRI data of individuals with ASD and matched TD, available within the ABIDE collection. To reduce the multiple sources of heterogeneity that impact on understanding the neural underpinnings of autistic condition, we selected a subgroup of 190 subjects (102 with ASD and 88 TD) according to the following criteria: male children (age range: 6.5–13 years); rs-fMRI data acquired with open eyes; data from the University sites that provided the largest number of scans (KKI, NYU, UCLA, UM). Connectivity values were evaluated as the linear correlation between pairs of time series of brain areas; then, a Linear kernel Support Vector Machine (L-SVM) classification, with an inter-site cross-validation scheme, was carried out. A permutation test was conducted to identify over-connectivity and under-connectivity alterations in the ASD group. The mean L-SVM classification performance, in terms of the area under the ROC curve (AUC), was 0.75 ± 0.05. The highest performance was obtained using data from KKI, NYU and UCLA sites in training and data from UM as testing set (AUC = 0.83). Specifically, stronger functional connectivity (FC) in ASD with respect to TD involve (p < 0.001) the angular gyrus with the precuneus in the right (R) hemisphere, and the R frontal operculum cortex with the pars opercularis of the left (L) inferior frontal gyrus. Weaker connections in ASD group with respect to TD are the intra-hemispheric R temporal fusiform cortex with the R hippocampus, and the L supramarginal gyrus with L planum polare. The results indicate that both under- and over-FC occurred in a selected cohort of ASD children relative to TD controls, and that these functional alterations are spread in different brain networks.
We propose and evaluate a procedure for the explainability of a breast density deep learning based classifier. A total of 1662 mammography exams labeled according to the BI-RADS categories of breast ...density was used. We built a residual Convolutional Neural Network, trained it and studied the responses of the model to input changes, such as different distributions of class labels in training and test sets and suitable image pre-processing. The aim was to identify the steps of the analysis with a relevant impact on the classifier performance and on the model explainability. We used the grad-CAM algorithm for CNN to produce saliency maps and computed the Spearman’s rank correlation between input images and saliency maps as a measure of explanation accuracy. We found that pre-processing is critical not only for accuracy, precision and recall of a model but also to have a reasonable explanation of the model itself. Our CNN reaches good performances compared to the state-of-art and it considers the dense pattern to make the classification. Saliency maps strongly correlate with the dense pattern. This work is a starting point towards the implementation of a standard framework to evaluate both CNN performances and the explainability of their predictions in medical image classification problems.
Background:
The integration of the information encoded in multiparametric MRI images can enhance the performance of machine-learning classifiers. In this study, we investigate whether the combination ...of structural and functional MRI might improve the performances of a deep learning (DL) model trained to discriminate subjects with Autism Spectrum Disorders (ASD) with respect to typically developing controls (TD).
Material and methods
We analyzed both structural and functional MRI brain scans publicly available within the ABIDE I and II data collections. We considered 1383 male subjects with age between 5 and 40 years, including 680 subjects with ASD and 703 TD from 35 different acquisition sites. We extracted morphometric and functional brain features from MRI scans with the Freesurfer and the CPAC analysis packages, respectively. Then, due to the multisite nature of the dataset, we implemented a data harmonization protocol. The ASD vs. TD classification was carried out with a multiple-input DL model, consisting in a neural network which generates a fixed-length feature representation of the data of each modality (FR-NN), and a Dense Neural Network for classification (C-NN). Specifically, we implemented a
joint fusion
approach to multiple source data integration. The main advantage of the latter is that the loss is propagated back to the FR-NN during the training, thus creating informative feature representations for each data modality. Then, a C-NN, with a number of layers and neurons per layer to be optimized during the model training, performs the ASD-TD discrimination. The performance was evaluated by computing the Area under the Receiver Operating Characteristic curve within a nested 10-fold cross-validation. The brain features that drive the DL classification were identified by the SHAP explainability framework.
Results
The AUC values of 0.66±0.05 and of 0.76±0.04 were obtained in the ASD vs. TD discrimination when only structural or functional features are considered, respectively. The
joint fusion
approach led to an AUC of 0.78±0.04. The set of structural and functional connectivity features identified as the most important for the two-class discrimination supports the idea that brain changes tend to occur in individuals with ASD in regions belonging to the Default Mode Network and to the Social Brain.
Conclusions
Our results demonstrate that the multimodal
joint fusion
approach outperforms the classification results obtained with data acquired by a single MRI modality as it efficiently exploits the complementarity of structural and functional brain information.
Purpose
This study aims at exploiting artificial intelligence (AI) for the identification, segmentation and quantification of COVID-19 pulmonary lesions. The limited data availability and the ...annotation quality are relevant factors in training AI-methods. We investigated the effects of using multiple datasets, heterogeneously populated and annotated according to different criteria.
Methods
We developed an automated analysis pipeline, the
LungQuant
system, based on a cascade of two U-nets. The first one (U-net
1
) is devoted to the identification of the lung parenchyma; the second one (U-net
2
) acts on a bounding box enclosing the segmented lungs to identify the areas affected by COVID-19 lesions. Different public datasets were used to train the U-nets and to evaluate their segmentation performances, which have been quantified in terms of the Dice Similarity Coefficients. The accuracy in predicting the CT-Severity Score (CT-SS) of the
LungQuant
system has been also evaluated.
Results
Both the volumetric DSC (vDSC) and the accuracy showed a dependency on the annotation quality of the released data samples. On an independent dataset (COVID-19-CT-Seg), both the vDSC and the surface DSC (sDSC) were measured between the masks predicted by
LungQuant
system and the reference ones. The vDSC (sDSC) values of 0.95±0.01 and 0.66±0.13 (0.95±0.02 and 0.76±0.18, with 5 mm tolerance) were obtained for the segmentation of lungs and COVID-19 lesions, respectively. The system achieved an accuracy of 90% in CT-SS identification on this benchmark dataset.
Conclusion
We analysed the impact of using data samples with different annotation criteria in training an AI-based quantification system for pulmonary involvement in COVID-19 pneumonia. In terms of vDSC measures, the U-net segmentation strongly depends on the quality of the lesion annotations. Nevertheless, the CT-SS can be accurately predicted on independent test sets, demonstrating the satisfactory generalization ability of the
LungQuant
.
Nowadays, MRI is focused on using ultra-high static magnetic fields (>7 T) to increase the signal-to-noise ratio. The use of high fields, on the other hand, requires novel technical solutions as well ...as more stringent design criteria for specific absorption rate levels, reducing radiative effect and coil resistance. In this paper, two flexible RF coils for 7 T human magnetic resonance, and 298 MHz ultra-high frequency operations were analyzed and characterized. Imaging of lower human limbs is regarded as a case study. The lumped element theory and subsequent numerical simulations were used to fine-tune the single-coil element and the dual-coil array design, respectively. Here, we demonstrate how the shape, size, configuration, and presence of the sample influence the coil performance. The penetration depth of the <inline-formula> <tex-math notation="LaTeX">B_{1} </tex-math></inline-formula>-field and the specific absorption rate values have been determined numerically using two numerical surface phantoms: saline and a multilayer human tissue. A preliminary study in the presence of a saline solution phantom has been carried out to develop and validate the dual-coil system. The frequency response of the dual-coil array was measured to assess its robustness when coupled to twelve human volunteers. We found that our design is robust to variations in the anatomical properties of the human thighs, and hence to coil bending. The presented approach can be useful for the implementation of flexible devices with high sensitivity levels and low specific absorption rate.