使用Python在数值网格上插值/填充间隙

Question

我使用的是MODIS陆地表面温度数据，它带有纬度和经度信息，在目前的格式下，我发现没有什么帮助。在

此信息存储整个图像的每个5x5像素窗口的中心坐标，而不是每个像素的坐标。很难解释，但本质上，这一信息比完整的MODIS图像小5倍，因为在垂直和水平方向上，每5个像素只有一个经纬度值。在

可以在这里找到一些示例数据：http://e4ftl01.cr.usgs.gov//MODIS_Dailies_F/MOLT/MOD11_L2.041/2015.08.15/MOD11_L2.A2015227.0215.041.2015228093331.hdf

我的问题如下：

如何插入/填充此信息之间的间隙以获得每个像素的纬度/经度，而不是5x5窗口？

如果MODIS数据中的纬度变量是这样的（当然简化了）：

[[1.0  2.0  3.0  4.0  5.0]
[ 1.1  2.1  2.1  4.1  5.1]
[ 1.2  2.2  3.2  4.2  5.2]
[ 1.3  2.3  3.3  4.3  5.3]
[ 1.4  2.4  3.4  4.4  5.4]]

如何使这个网格大5倍（25 x 25），以便在水平和垂直方向的值之间进行插值？在

Answer 1

我遇到了完全相同的问题，遇到了python-geotiepoints，它正好满足您的需要，唯一的依赖项是numpy和scipy。你只要做到：

from geotiepoints import modis5kmto1km

lons_1km, lats_1km = modis5kmto1km(lon_5km, lat_5km)

它会根据需要进行插值。在