A system and related methods for gamma emission breast tomosynthesis in which a set of two-dimensional images of a breast taken at different angular views are reconstructed into a three-dimensional map of the breast. The system applies an expectation maximization technique having integrated regularization, resolution recovery and attenuation correction to improve the clarity of the resulting three-dimensional map.