用python将卫星沿轨数据插值到等经纬度网格上

两米阳光

我想将卫星沿轨数据（有经纬度信息，但是不规则，不是持续递增或递减）插值到等经纬度网格上，查阅了很久，尝试了多种方法，都失败了，但是有interpolate.griddata和interpolate.NearestNDInterpolator还算勉强出结果，但依然是不对的，原来没有数据的地方插值出了数据，想请教各位，这是什么原因，感谢！主要代码如下：
spd, lat, lon = read_nc(filename)

lon = lon.ravel()
lat = lat.ravel()
spd = spd.ravel()lat0 = np.arange(-90.0, 90, 0.5)lon0 = np.arange(0.0, 360, 0.5)lon_new, lat_new = np.meshgrid(lon0, lat0)spd_new = interpolate.griddata((lon, lat), spd, (lon_new, lat_new), method='nearest').reshape(360, 720)

两米阳光

  主要代码行：
   spd, lat, lon = read_nc(filename)

    lon = lon.ravel()
    lat = lat.ravel()
    spd = spd.ravel()
    lat0 = np.arange(-90.0, 90, 0.5)
    lon0 = np.arange(0.0, 360, 0.5)
    lon_new, lat_new = np.meshgrid(lon0, lat0)
    spd_new = interpolate.griddata((lon, lat), spd, (lon_new, lat_new), method='nearest').reshape(360, 720)

myj0219

6666666666666666666

MissLH88

这个插值我老师之前和我说他们一般用双线性插值，我当时就自己写的函数，算的很慢，但是结果看着还挺正常