使用gdal/python重新投影和拼接MODIS level2

Question

我有一个.hdf格式的MODIS大气产品（MOD04）数据集，我想重新投影并保存为Geotiff。我知道使用gdalwalp是可能的，但是问题是在保存图像之前我需要做一些操作。基本上，我需要从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条带数据重新采样到栅格中。例如，您可以使用MRTSwath或{a1}之类的工具来实现这一点。在

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