Postgis Raster: Загрузка Tif-Subdataset

Question

Я хотел бы загрузить подканал из файла TIF в мою базу данных PostGIS.

gdalinfo дает мне следующие Infos о TIF-файл:

Driver: GTiff/GeoTIFF 
Files: pathto/imagename.tif 
Size is 611, 1 
Coordinate System is `' 
Metadata: 
    TIFFTAG_IMAGEDESCRIPTION=CXmlTiff type 
    TIFFTAG_RESOLUSTIONUNIT=2 (pixels/inch) 
    TIFFTAG_XRESOLUTION=100 
    TIFFTAG_YRESOLUTION=100 
ImageStructure Metadata: 
    INTERLEAVE=BAND 
Subdatasets: 
    SUBDATASET_1_NAME=GTIFF_DIR:1:pathto/imagename.tif 
    SUBDATASET_1_DESC=Page 1 (611P x 1L x 1B) 
    SUBDATASET_2_NAME=GTIFF_DIR:2:pathto/imagename.tif 
    SUBDATASET_2_DESC=Page 1 (392P x 514352L x 1B) 
Corner Coordinates: 
Upper Left ( 0.0, 0.0) 
Lower Left ( 0.0, 1.0) 
Upper Right (611.0, 0.0) 
Lower Right (611.0, 1.0) 
Center  (305.5, 0.5) 
Band 1 Block=611x1 Type=Byte, ColorInterp=Gray 

Я хотел бы, чтобы загрузить вторую subdataset (392x514352).

Прямо сейчас, используя raster2pgsql и имя файла в качестве входных данных, загружается только первый поднабор, второй игнорируется. Существует только одна полоса, поэтому использование -b-переключателя не работает. Любой совет?

Большое спасибо

Томас

Answer 1

Subdatasets не распространены в растров, и я никогда не видел их в TIFFs (обычно, вы видите полосы, не subdatasets). Он выглядит как a similar bug was fixed for HDF5 files, но это, вероятно, не поможет вашему GTiff. Я бы посоветовал вам создать биг-код для этой функции.

В то же время gdal_translate может использоваться с опцией -sds, чтобы скопировать все поднаборные файлы вашего файла в отдельные выходные файлы, которые вы могли бы использовать с raster2pgsql.