轮廓函数估计公式

Question

我计算了一些空间数据的等高线图（x = lon，y = lat，z）。

install.packages("akima")
library(akima)
lon <- runif(n=200,min=-180,max=180)
lat <- runif(n=200,min=-90,max=90)
z <- rnorm(n=200,mean=20,sd=1)
xo <- (-180:180)
yo <- (-90:90)
int <- interp(lon,lat,z,duplicate="mean")
contour(int,nlevels=50)

现在我正在寻找获得这些轮廓背后的数学函数（公式）的方法。换句话说，我想得到z = f（lat，lon）的函数来计算任何lat和lon值的z估计值。

有什么想法吗？提前谢谢！