Objectives. (1) Develop a computational pipeline for three-dimensional
fast neural magnetic detection electrical impedance tomography (MDEIT), (2)
determine whether constant current or constant voltage is preferable for MDEIT,
(3) perform reconstructions of simulated neural activity in a human head
model with realistic noise and compare MDEIT to EIT and (4) perform a
two-dimensional study in a saline tank for MDEIT with optically pumped
magnetometers (OPMs) and compare reconstruction algorithms. Approach.
Forward modelling and image reconstruction were performed with a realistic
model of a human head in three dimensions and at three noise levels for
four perturbations representing neural activity. Images were compared using
the error in the position and size of the reconstructed perturbations. Two-
dimensional MDEIT was performed in a saline tank with a resistive perturbation
and one OPM. Six reconstruction algorithms were compared using the error
in the position and size of the reconstructed perturbations. Main Results.
A computational pipeline was developed in COMSOL Multiphysics, reducing
the Jacobian calculation time from months to days. MDEIT reconstructed
images with a lower reconstruction error than EIT with a mean difference of
7.0%, 5.5% and 11% for three noise cases representing current noise, reduced
current source noise and reduced current source and magnetometer noise. A
rank analysis concluded that the MDEIT Jacobian was less rank-deficient than
the EIT Jacobian. Reconstructions of a phantom in a saline tank had a best
reconstruction error of 13%, achieved using 0th order Tikhonov regularisation
with simulated noise-based correction. Significance. This study demonstrated
that three-dimensional MDEIT for neural imaging is feasible and that MDEIT
reconstructed superior images to EIT, which can be explained by the lesser rank
deficiency of the MDEIT Jacobian. Reconstructions of a perturbation in a saline
tank demonstrated a proof of principle for two-dimensional MDEIT with OPMs
and identified the best reconstruction algorithm.