The QSM reconstruction comprised the total field estimation using an adaptive, quadratic fit of the mGRE phase,48 (link) local field calculation with the projection onto dipole fields (PDF) method,47 (link) and susceptibility estimation using the morphology-enabled dipole inversion with automatic uniform CSF zero-reference algorithm.46 (link),52 (link)-54 (link) In stroke patients, the total field was estimated using a linearfit of the mGRE phase,55 (link) as 3D flow compensation was not available on the clinical scanner. The CBF maps (mL/100 g/min) were generated from the ASL data using FuncTool (GE Healthcare). All images were co-registered and interpolated to the resolution of the QSM maps using the FSL FLIRT algorithm.56 (link),57 (link)