std_devs = np.diag(np.sqrt(cov))
kern_nx, kern_ny = np.round(scotts_factor * 2 * np.pi * std_devs)
Implementing a 2D, FFT-based Kernel Density Estimator in python, and comparing it to the SciPy implimentation