Fit of a 2D histogram by 2D gaussians

I have a data distribution in a 2D histogram.
I need this histogram to be fitted by a sum of 2D gaussians.