Large medical imaging data sets are becoming increasingly available. A common challenge in these data sets is to ensure that each sample meets minimum quality requirements devoid of significant artefacts. Despite a wide range of existing automatic methods having been developed to identify imperfections and artefacts in medical imaging, they mostly rely on data-hungry methods. In particular, the scarcity of artefact-containing scans available for training has been a major obstacle in the development and implementation of machine learning in clinical research. To tackle this problem, we propose a novel framework having four main components: (1) a set of artefact generators inspired by magnetic resonance physics to corrupt brain MRI scans and augment a training dataset, (2) a set of abstract and engineered features to represent images compactly, (3) a feature selection process that depends on the class of artefact to improve classification performance, and (4) a set of Support Vector Machine (SVM) classifiers trained to identify artefacts. Our novel contributions are threefold: first, we use the novel physics-based artefact generators to generate synthetic brain MRI scans with controlled artefacts as a data augmentation technique. This will avoid the labour-intensive collection and labelling process of scans with rare artefacts. Second, we propose a large pool of abstract and engineered image features developed to identify 9 different artefacts for structural MRI. Finally, we use an artefact-based feature selection block that, for each class of artefacts, finds the set of features that provide the best classification performance. We performed validation experiments on a large data set of scans with artificially-generated artefacts, and in a multiple sclerosis clinical trial where real artefacts were identified by experts, showing that the proposed pipeline outperforms traditional methods. In particular, our data augmentation increases performance by up to 12.5 percentage points on the accuracy, F1, F2, precision and recall. At the same time, the computation cost of our pipeline remains low - less than a second to process a single scan - with the potential for real-time deployment. Our artefact simulators obtained using adversarial learning enable the training of a quality control system for brain MRI that otherwise would have required a much larger number of scans in both supervised and unsupervised settings. We believe that systems for quality control will enable a wide range of high-throughput clinical applications based on the use of automatic image-processing pipelines.

Bootstrapping of repeated diffusion-weighted image datasets enables nonparametric quantification of the uncertainty in the inferred fiber orientation. The wild bootstrap and the residual bootstrap are model-based residual resampling methods which use a single dataset. Previously, the wild bootstrap method has been presented as an alternative to conventional bootstrapping for diffusion tensor imaging. Here we present a study of an implementation of model-based residual bootstrapping using q -ball analysis and compare the outputs with conventional bootstrapping. We show that model-based residual bootstrap q-ball generates results that closely match the output of the conventional bootstrap. Both the residual and conventional bootstrap of multifiber methods can be used to estimate the probability of different numbers of fiber populations existing in different brain tissues. Also, we have shown that these methods can be used to provide input for probabilistic tractography, avoiding existing limitations associated with data calibration and model selection.

This work shows that multi-fibre reconstruction techniques, such as Persistent Angular Structure (PAS) MRI or QBall Imaging, provide much more information than just discrete fibre orientations, which is all that previous tractography algorithms exploit from them. We show that the shapes of the peaks of the functions output by multiple-fibre reconstruction algorithms reflect the underlying distribution of fibres. Furthermore, we show how to exploit this extra information to improve Probabilistic Index of Connectivity (PICo) tractography. The method uses the Bingham distribution to model the uncertainty in fibre-orientation estimates obtained from peaks in the PAS or QBall Orientation Distribution Function (ODF). The Bingham model captures anisotropy in the uncertainty, allowing the method to track through fanning and bending structures, which previous methods do not recover reliably. We devise a new calibration procedure to construct a mapping from peak shape to Bingham parameters. We test the accuracy of the calibration using a bootstrap experiment. Finally, we show that exploiting the peak shape in this way can provide improved PICo tractography results.

A method is presented for determining paths of anatomical connection between regions of the brain using magnetic resonance diffusion tensor information. Level set theory, applied using fast marching methods, is used to generate three-dimensional time of arrival maps, from which connection paths between brain regions may be identified. The method is demonstrated in the normal brain and it is shown that major white matter tracts may be elucidated and that multiple connections and tract branching are allowed. Maps of connectivity between brain regions are also determined. Four options are described for estimating the degree of connectivity between regions.