python - 如何使用 Python GDAL 镶嵌 MODIS hdf 文件？

Question

我已经下载了大量 MODIS hdf 文件并从每个文件中提取了所需的波段。如何使用 Python 中的 GDAL lib 将所有 hdf 文件拼接在一起并重新投影？谢谢！

Answer 1

使用两个 GDAL 命令行工具：

1. gdal_merge.py将波段镶嵌成一个栅格
2. gdalwarp重新投影光栅

Answer 2

您也可以尝试 pyModis [0]，它在开发版本中使用 Python GDAL（从今天开始，我需要对其进行更多测试）。

gdal_merge.py 不适用于 MODIS HDF（可能是一个错误），这是返回的错误

Traceback (most recent call last):
  File "/usr/local/bin/gdal_merge.py", line 510, in <module>
    sys.exit(main())
  File "/usr/local/bin/gdal_merge.py", line 390, in main
    file_infos = names_to_fileinfos( names )
  File "/usr/local/bin/gdal_merge.py", line 122, in names_to_fileinfos
    if fi.init_from_name( name ) == 1:
  File "/usr/local/bin/gdal_merge.py", line 147, in init_from_name
    self.band_type = fh.GetRasterBand(1).DataType
AttributeError: 'NoneType' object has no attribute 'DataType'

[0] https://github.com/lucadelu/pyModis