Integrate 2D kernel density estimate

Question

I have a x,y distribution of points for which I obtain the KDE through scipy.stats.gaussian_kde. This is my code and how the output looks (the

Answer 1

Currently, it is available

kernel.integrate_box([-np.inf,-np.inf], [2.5,1.5])