Как создать мозаику из файлов MODIS hdf с помощью Python GDAL?

Я загрузил большое количество HDF-файлов MODIS и извлек необходимые полосы из каждого файла. Как совместить все hdf-файлы вместе с помощью GDAL lib в Python и перепроецировать их? Спасибо!


person oops    schedule 22.01.2013    source источник
comment
Попробуйте gdal_merge.py, а затем gdalwarp, чтобы перепроецировать его.   -  person Mike T    schedule 22.01.2013
comment
Хорошо, я мог бы также оставить это в области ответов ...   -  person Mike T    schedule 23.01.2013


Ответы (2)


Используйте два инструмента командной строки GDAL:

  1. gdal_merge.py, чтобы объединить каналы в один растр.
  2. gdalwarp для перепроецирования растра
person Mike T    schedule 23.01.2013

Вы также можете попробовать 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

person lucadelu    schedule 26.04.2014
comment
Если gdal не может открыть файл, возможно, вы указываете неправильный указатель на набор данных или вам не хватает драйвера, необходимого для открытия типа файла. - person Rutger Kassies; 28.04.2014