cov = np.cov(xyi)
scotts_factor = np.power(n, -1.0 / 6)
std_devs = np.diag(np.sqrt(cov))
kern_nx, kern_ny = np.round(scotts_factor * 2 * np.pi * std_devs)
inv_cov = np.linalg.inv(cov * scotts_factor ** 2)
Внедрение 2D-оценщика плотности ядра на основе БПФ в python и сравнение его с реализацией SciPy