Find the roots of a bivariate Spline on a rectangular grid using scipy.optimize

Question

I have data points on a two-dimensional grid representing some function g(x, y) that I do not know the analytical form of. I want to do a spline interpolation and find the p