aka non-parametric blobology.

opit got together with a mathematician colleague and wrote down how the 2D Gaussian KDE works and how it is implemented in scipy. A .pdf with the equations and comments is here.

[illustration from Wikipedia.]