[gdal-dev] Multidimensional raster support in GDAL

[Ari Jolma]

Even Rouault kirjoitti 21.11.2017 klo 16:06:
>
> gdal_translate (actually GDALDatasetCopyWholeRaster()) uses the 
> INTERLEAVE=PIXEL/BAND metadata item of the IMAGE_STRUCTURE domain as a 
> hint whether to read band by band or all bands together
>
> So if you add the following, it should request all bands :
>
> poDS->SetMetadataItem("INTERLEAVE", "PIXEL", "IMAGE_STRUCTURE")
>

I don't see that happening in the case of WCS driver:

PIXEL_INTERLEAVE=TRUE is my option to set the metadata item but it could 
be set with option INTERLEAVE=PIXEL, which GDALDatasetCopyWholeRaster() 
examines.

8<-------------------
gdal_translate -oo PIXEL_INTERLEAVE=TRUE -oo OriginAtBoundary -oo 
CACHE=wcs_cache -srcwin 0 0 2 2 
"WCS:http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?version=1.0.0&coverage=BGS_EMODNET_CentralMed-MCol" 
x
GNM: GNMRegisterAllInternal
GNM: RegisterGNMFile
GNM: RegisterGNMdatabase
Setting INTERLEAVE to PIXEL
GDAL: 
GDALOpen(WCS:http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?version=1.0.0&coverage=BGS_EMODNET_CentralMed-MCol, 
this=0x1256210) succeeds as WCS.
Input file size is 2977, 3883
0GTiff: Reopen with strip chop enabled
GDAL: GDAL_CACHEMAX = 585 MB
GDAL: GDALDatasetCopyWholeRaster(): 2*2 swaths, bInterleave=1
URL=http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?SERVICE=WCS&REQUEST=GetCoverage&VERSION=1.0.0&COVERAGE=BGS_EMODNET_CentralMed-MCol&FORMAT=GTiff&BBOX=9.83125,44.05624829,14.09792008,46.18958333&WIDTH=1024&HEIGHT=512&CRS=EPSG:4326&bands=1
HTTP: 
Fetch(http://194.66.252.155/cgi-bin/BGS_EMODnet_bathymetry/ows?SERVICE=WCS&REQUEST=GetCoverage&VERSION=1.0.0&COVERAGE=BGS_EMODNET_CentralMed-MCol&FORMAT=GTiff&BBOX=9.83125,44.05624829,14.09792008,46.18958333&WIDTH=1024&HEIGHT=512&CRS=EPSG:4326&bands=1)
HTTP: libcurl/7.47.0 OpenSSL/1.0.2g zlib/1.2.8 libidn/1.32 librtmp/2.3
WCS: GDALOpenResult() on content-type: image/tiff
GDAL: GDALOpen(/vsimem/wcs/0x1256210/wcsresult.dat, this=0x12a2610) 
succeeds as GTiff.
ERROR 1: Returned tile does not match expected band configuration.
Got 3 bands instead of one although the coverage has band range type.

GDAL: GDALClose(/vsimem/wcs/0x1256210/wcsresult.dat, this=0x12a2610)
ERROR 1: wcs_cache/znQnc.xml, band 1: IReadBlock failed at X offset 0, Y 
offset 0: Returned tile does not match expected band configuration.
Got 3 bands instead of one although the coverage has band range type.

GDAL: GDALClose(wcs_cache/znQnc.xml, this=0x1256210)
GDAL: In GDALDestroy - unloading GDAL shared library.
8<-------------------



Ari

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20171122/159191be/attachment.html>