A computer-implemented method enables a proton density distribution to be obtained. In one embodiment, the method comprises acquiring nuclear magnetic resonance data from porous media; inverting the nuclear magnetic resonance data via a global optimization algorithm to determine a proton density distribution within the porous media; and outputting the determined proton density distribution.