HOMOR: Higher Order Model Outlier Rejection for high b-value MR diffusion data

Diffusion MR images are prone to artefacts caused by head movement and cardiac pulsation. Previous techniques for the automated voxel-wise detection of signal intensity outliers have relied on the fit of the diffusion tensor to the data (RESTORE). However, the diffusion tensor cannot appropriately model more than a single fibre population, which may lead to inaccuracies when identifying outlier voxels in crossing fibre regions, particularly when high b-values are used to obtain increased angular contrast. HOMOR (higher order model outlier rejection) was developed to overcome this limitation and is introduced in this study. HOMOR is closely related to RESTORE, but employs a higher order model capable of resolving multiple fibre populations within a voxel. Using high b-value (b=3000 s/mm2) diffusion data from a population of 90 healthy participants, as well as simulations, HOMOR was found to identify a decreased number of outlier voxels compared to RESTORE primarily within areas of crossing, bending and fanning fibres. At lower b-values, however, RESTORE and HOMOR give similar results, which is demonstrated using diffusion data acquired at b=1000 s/mm2 in a mixed cohort. This study demonstrates that, although RESTORE is suitable for low b-value data, HOMOR is better suited for high b-value data.