This paper introduces a new approach, which combines empirical Bayes modeling with recent advances in Markov chain Monte Carlo filters for hidden Markov models, to address long-standing challenging problems in reconstruction of 3D images, with uncertainty quantification, from noisy 2D pixels in cryogenic electron microscopy and other applications such as brain network development in infants.