使用gdal / python重新投影和镶嵌MODIS level2

Question

我有.hdf格式的MODIS大气产品（MOD04）的数据集，我想重新投影并将它们保存为Geotiff。我知道可以使用gdalwarp来实现它，但问题是我需要在保存图像之前进行一些操作。基本上，我需要的是从hdf文件中读取不同的子数据集，使用numpy进行一些修改，最后将数组保存为Geotiff文件。

我不是使用gdal的专家，但我一般都知道它是如何工作的。我必须注意，当我用gdal打开MOD04 hdf文件时，似乎gdal无法读取其GeoProjection和GeoTransform数据！

>>> src = 'D:/Test/MYD04_L2.A2007231.0745.006.2014054160056.hdf'
>>> dataset = gdal.Open(src,gdal.GA_ReadOnly)
>>> print dataset.GetProjection()

>>> print dataset.GetGeoTransform()
(0.0, 1.0, 0.0, 0.0, 0.0, 1.0)

以上是阅读我的hdf文件的示例。当我直接打开子数据集时，它是相同的。我想知道是否有人知道如何解决这个问题？

Answer 1

我不认为您的数据已经投射过，因此使用gdalwarp进行重新投影是行不通的。首先需要将MODIS Swath数据重新采样到网格中。您可以使用MRTSwath或PyResample等工具执行此操作。

理论上gdalwarp应该能够将lat / lon网格作为地理定位阵列处理，但由于扫描传感器的蝴蝶结效应，我认为它不能正常工作。