Question

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

Answer 1

使用两个GDAL命令行工具：

gdal_merge.py将乐队拼接成一个栅格
gdalwarp重新投影栅格

Answer 2

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

gdal_merge.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

如何使用Python GDAL拼接MODIS hdf文件？

2 个答案: