Data source esotici per GDAL
Da Wikipedia, l'enciclopedia libera.
Data source esotici per GDAL
Contrariamente a quanto riportato correntemente sui manuali di gdal_translate e di r.in.gdal, il raster map di input non è necessariamente un file, ma segue una sintassi che è strettamente legata al relativo driver che viene usato da gdal. Per formati particolarmente complessi che sono multi dataset e multibanda come HDF, si scopre che senza la corretta sintassi l'importazione nel DB di GRASS fallisce miseramente.
Mi sono imbattutto per la prima volta in questa feature di GDAL cercando di importare i file HDF di un sensore Hyperion (iperspettrale a 242 bande, distribuito in formato HDF4). Il prodotto in questione si trasforma in un GeoTiff con un comando del tipo seguente:
gdal_translate -sds HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":0 EO1H0220392005249110KM.tif
dove EO1H0220392005249110KM.L1R è il filename del dataset in questione. Lo stesso tipo di sintassi può essere utilizzata direttamente da GRASS, ad esempio
r.in.gdal -e in=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":0 out=EO1H0220392005249110KM
L'informazione sul subset da consultare lo si desume utilizzando gdalinfo che restituisce per il prodotto L1R di questo caso i subset:
Subdatasets: SUBDATASET_1_NAME=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":0 SUBDATASET_1_DESC=[7210x242x256] EO1H0220392005249110KM.L1R (16-bit integer) SUBDATASET_2_NAME=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":1 SUBDATASET_2_DESC=[242x256] Spectral Center Wavelengths (32-bit floating-point) SUBDATASET_3_NAME=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":2 SUBDATASET_3_DESC=[242x256] Spectral Bandwidths (32-bit floating-point) SUBDATASET_4_NAME=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":3 SUBDATASET_4_DESC=[242x256] Gain Coefficients (32-bit floating-point) SUBDATASET_5_NAME=HDF4_SDS:HYPERION_L1:"EO1H0220392005249110KM.L1R":4 SUBDATASET_5_DESC=[7210x242x256] Flag Mask (8-bit unsigned integer)
--FrancescoLovergine 00:12, Mar 22, 2007 (CET)
