Magnetic resonance current density imaging (MRCDI) and MR electrical impedance tomography (MREIT) are two emerging modalities, which combine weak time-varying currents injected via surface electrodes with magnetic resonance imaging (MRI) to acquire information about the current flow and ohmic conductivity distribution at high spatial resolution. The injected current flow creates a magnetic field in the head, and the component of the induced magnetic field ΔB parallel to the main scanner field causes small shifts in the precession frequency of the magnetization. The measured MRI signal is modulated by these shifts, allowing to determine ΔB for the reconstruction of the current flow and ohmic conductivity. Here, we demonstrate reliable ΔB measurements in-vivo in the human brain based on multi-echo spin echo (MESE) and steady-state free precession free induction decay (SSFP-FID) sequences. In a series of experiments, we optimize their robustness for in-vivo measurements while maintaining a good sensitivity to the current-induced fields. We validate both methods by assessing the linearity of the measured ΔB with respect to the current strength. For the more efficient SSFP-FID measurements, we demonstrate a strong influence of magnetic stray fields on the ΔB images, caused by non-ideal paths of the electrode cables, and validate a correction method. Finally, we perform measurements with two different current injection profiles in five subjects. We demonstrate reliable recordings of ΔB fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔB measurements with simulated ΔB images based on FEM calculations and individualized head models reveals significant linear correlations in all subjects, but only for the stray field-corrected data. As final step, we reconstruct current density distributions from the measured and simulated ΔB data. Reconstructions from non-corrected ΔB measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔB measurements and from simulated ΔB images reveals an average coefficient of determination R of 71%. In addition, it shows that the simulations underestimated the current strength on average by 24%. Our results open up the possibility of using MRI to systematically validate and optimize numerical field simulations that play an important role in several neuroscience applications, such as transcranial brain stimulation, and electro- and magnetoencephalography.