The paper focuses on a new method for the inference of a parametric random spheroid from the observations of its 2D orthogonal projections. Such a stereological problem is well-known from the literature when the projections come from only one deterministic spheroid. Nevertheless, when the spheroid is random itself, the estimation of its distribution is not straightforward. From a theoretical viewpoint, it is shown that the semi-axes of the spheroid and the ones of the projected ellipses are linked through a random polynomial of degree two which admits two real random positive roots. The likelihood can be formulated in terms of the coefficients of the random polynomial, but is not analytically tractable. Assuming that the random spheroid is parameterized by a set of parameters θ real , an approximation of the maximum a posteriori is used to estimate θ real. The estimator is based on the so-called approximate Bayesian computation method and a kernel density technique. As an illustration, the case of a spheroids population, whose semi-major axis follows a gamma distribution and the flattening coefficient a truncated normal distribution, is studied. The numerical results demonstrate that the bias of the estimator is very low, with a reasonable variance, both for the first and the second order moments of the semi-axes. The proposed method enables to recover some 3D morphological characteristics of a population of independent and