We tackle the reconstruction problem of density images from indirect measurements with a novel variational approach: By implementing an appropriate modelling of the mass-conserving density transformation in the reconstruction process we obtain the first building block of our variational method.Suitable regularization for images with edges(total variation)and for reasonable deformations(hyperelastic)without self folding completes the functional.We focus on obtaining analytical results and conclude the talk with applications to cardiac PET.