In a maximum a posteriori tomosynthesis method to reconstruct a three-dimensional image from two-dimensional x-ray images, a Geman prior function to reduce the noise is used in which the edges and boundaries of tissue structures remain visible. The method is parameterized by the estimated noise value of the attenuation coefficients and the estimated average tissue attenuation value.