Abstract: In the past decade, low-field NMR relaxation and diffusion measurements in grossly inhomogeneous fields have been used to characterize pore size distribution of porous media. Estimation of these distributions from the measured magnetization data plays a central role in the inference of insitu petro-physical and fluid properties such as porosity, permeability, and hydrocarbon viscosity. In general, inversion of the relaxation and/or diffusion distribution from NMR data is a non-unique and ill-conditioned problem. It is often solved in the literature by finding the smoothest relaxation distribution that fits the measured data by use of regularization. In this paper, estimation of these distributions is further constrained by linear functionals of the measurement that can be directly estimated from the measured data. These linear functionals include Mellin, Fourier–Mellin, and exponential Haar transforms that provide moments, porosity, and tapered areas of the distribution, respectively. The addition of these linear constraints provides more accurate estimates of the distribution in terms of a reduction in bias and variance in the estimates. The resulting distribution is also more stable in that it is less sensitive to regularization. Benchmarking of this algorithm on simulated data sets shows a reduction of artefacts often seen in the distributions and, in some cases, there is an increase of resolution in the features of the T 2 distribution. This algorithm can be applied to data obtained from a variety of pulse sequences including CPMG, inversion and saturation recovery and diffusion editing, as well as pulse sequences often deployed down-hole. [Copyright &y& Elsevier]