如何使用pythongdal拼接MODIS hdf文件？

Question

我下载了大量的MODIS hdf文件，并从每个文件中提取了所需的波段档案。怎么办我可以使用Python中的GDAL lib将所有hdf文件拼接在一起并重新投影吗？谢谢！在

Answer 1

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

格达尔_合并.py不适用于MODIS HDF（可能是一个bug），这是返回的错误

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

Answer 2

使用两个GDAL命令行工具：

1. gdal_merge.py拼接成一个光栅
2. gdalwarp重新投影光栅