如何在大型曲线网格上插入2D函数？

Question

我有一些表面数据为x_mesh, y_mesh, z_mesh。

x_mesh和y_mesh是根据地理位置上的经度和纬度逐步生成的（因此，它们的形状相同且规则，但x和y的步长不是）。

z_mesh是地图中心从切线平面到地球椭球的高度。

我可以使用matplotlib.pyplot.pcolormesh(x_mesh, y_mesh, z_mesh)轻松绘制曲面。它有效。

现在我想通过鼠标设置一条线，并以某种方式在此线下取一个Z轮廓。我需要一些内插器来制作xy - ＆gt; z，但不知道要使用哪一个。

我试着这样做：

scipy.interpolate.interp2d(x_mesh, y_mesh, z_mesh)

但它给了我一个错误：OverflowError: Too many data points to interpolate

您对如何插入此类数据有任何想法吗？

P.S。 geotiff不是很大，它是6K x 6K像素。我看到pcolormesh以某种方式插入像素之间的颜色值，如果我放大。