使用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

我的问题如下：

如何插入/填充此信息之间的间隙以获得每个像素的Lat / Lon，而不是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)

它将根据需要进行插值。