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 ΔBz,c 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 ΔBz,c for the reconstruction of the current flow and ohmic conductivity. Here, we demonstrate reliable ΔBz,c 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 ΔBz,c with respect to the current strength. For the more efficient SSFP-FID measurements, we demonstrate a strong influence of magnetic stray fields on the ΔBz,c 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 ΔBz,c fields as weak as 1 nT, caused by currents of 1 mA strength. Comparison of the ΔBz,c measurements with simulated ΔBz,c 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 ΔBz,c data. Reconstructions from non-corrected ΔBz,c measurements systematically overestimate the current densities. Comparing the current densities reconstructed from corrected ΔBz,c measurements and from simulated ΔBz,c images reveals an average coefficient of determination R2 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.