Abstract

Purpose To create a realistic in-silico head phantom for the second QSM Reconstruction Challenge and for future evaluations of processing algorithms for Quantitative Susceptibility Mapping (QSM). Methods We created a whole-head tissue property model by segmenting and post-processing high-resolution, multi-parametric MRI data acquired from a healthy volunteer. We simulated the steady-state magnetization using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based post-processing. We demonstrated some of the phantom’s properties, including the possibility of generating phase data that do not evolve linearly with echo time due to partial volume effects or complex distributions of frequency shifts within the voxel. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the model, such as the inclusion of pathologies, as well as the simulation of a wide range of acquisition protocols. Results The brain-part of the phantom features realistic morphology combined with realistic spatial variations in relaxation and susceptibility values. Simulation code allows adjusting the following parameters and effects: repetition time and echo time, voxel size, background fields, and RF phase biases. Additionally, diffusion weighted imaging data of the phantom is provided allowing future investigations of tissue microstructure effects in phase and QSM algorithms. Conclusion The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative MRI reconstruction algorithms.