1、Delaunay三角剖分
from scipy.spatial.qhull import Delaunay
from scipy.interpolate import LinearNDInterpolator, griddata
tri = Delaunay(np.asarray([orig_lon.ravel(), orig_lat.ravel()]).T)
interpolator = LinearNDInterpolator(tri, orig_val.ravel())
2、scipy.interpolate.LinearNDInterpolator線性插值
https://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.LinearNDInterpolator.html