Как преобразовать HDF5 в geotiff с помощью gdal_translate или gdalwarp?

У меня есть HDF5 без географической привязки 10° на 10°, который я хочу преобразовать в geotiff. Я пробовал 2 подхода (1) gdalwarp и (2) gdal_translate, но не смог найти, как преобразовать файлы.

1) gdalwarp

gdalwarp -t_srs '+proj=longlat +datum=WGS84 +no_defs ' -te 40 30 50 40 -tr 0.0089285714 -0.0089285714 -multi HDF5:"/PathToTheHDF/FileName.HDF5"://SUBDATASET /PathToOutput/FileName_Subdataset.tif

У меня есть следующая ошибка:

ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates for ...
There is no affine transformation and no GCPs.
256

2) gdal_translate

gdal_translate -of GTiff -a_srs '+proj=longlat +datum=WGS84 +no_defs ' -a_ullr 40 40 50 30 HDF5:"/PathToTheHDF/FileName.HDF5"://SUBDATASET /PathToOutput/FileName_Subdataset.tif

и у меня есть следующая ошибка:

Input file size is 1120, 1120
0HDF5-DIAG: Error detected in HDF5 (1.8.11) thread 140288255100800:
  #000: ../../../src/H5Dio.c line 182 in H5Dread(): can't read data
    major: Dataset
    minor: Read failed
  #001: ../../../src/H5Dio.c line 550 in H5D__read(): can't read data
    major: Dataset
    minor: Read failed
  #002: ../../../src/H5Dchunk.c line 1837 in H5D__chunk_read(): unable to read raw data chunk
    major: Low-level I/O
    minor: Read failed
  #003: ../../../src/H5Dchunk.c line 2868 in H5D__chunk_lock(): data pipeline read failed
    major: Data filters
    minor: Filter operation failed
  #004: ../../../src/H5Z.c line 1150 in H5Z_pipeline(): required filter 'szip' is not registered
    major: Data filters
    minor: Read failed
  #005: ../../../src/H5PL.c line 293 in H5PL_load(): search in paths failed
    major: Plugin for dynamically loaded library
    minor: Can't get value
  #006: ../../../src/H5PL.c line 397 in H5PL__find(): can't open directory
    major: Plugin for dynamically loaded library
    minor: Can't open directory or file
ERROR 1: H5Dread() failed for block.
ERROR 1: HDF5:/export/GIO/gww-probav-test/modis/2014/01/21/MC10GWW_X22Y09_20140121_1KM.HDF5://HUE, band 1: IReadBlock failed at X offset 0, Y offset 0
ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0

Любая идея о том, что не так или как мне преобразовать эти данные?


hdf5 gdal geotiff
person Wraf    schedule 08.04.2014    source источник

Ответы (3)


arrow_upward
3
arrow_downward

Использование R и gdal_translate:

library(gdalUtils)

band <- 1

file_path <- "/myfolder/myfile.hdf"

system(paste0('gdal_translate ', get_subdatasets(file_path)[band],' myfile_', band','.tif'))
person DJack    schedule 25.08.2015

arrow_upward
10
arrow_downward 

gdal_translate your_HDF_file -sds output.tif

Параметр -sds разделит ваш многоканальный файл HDF5 на несколько одноканальных tif-файлов.

person Generic Wevers    schedule 16.09.2015

arrow_upward
0
arrow_downward

Попробуй это:

gdalwarp your_HDF_file_sds C:\Save\your\file\sample.tif
person doo    schedule 24.07.2014