Dear gdal-developers,
I noticed a strange behaviour while reading Pleaides Jpeg2000 4 bands EO
products using overviews. I use GDAL 1.10.1 with Erdas ECW4.3 to support
jpeg2000 files. The following code
GDALDataset * dataset = static_cast<GDALDataset *>(GDALOpen(inputFilename,
GA_ReadOnly));
< Image has : 4 bands, and WxH : 10296, 6573>
< Image DataType : GDT_UInt16>
<define parameters : offset={0,0}, desiredNbOfSamples=167424,
data.width=500, data.height=400>
// Just get and show twice the same data :
for (int index = 0; index < 2; index++)
{
int pixelSizeInBytes = nbBands * 8;
int rowSizeInBytes = pixelSizeInBytes * data.width;
for (int band=1;band<=nbBands;band++)
{
GDALRasterBand* rasterBand = dataset->GetRasterBand(band)->
GetRasterSampleOverview(desiredNbOfSamples);
rasterBand->RasterIO( GF_Read, offset[0], offset[1], data.width, data.
height,
data.buffer + 8*(band - 1), data.width, data.
height, GDT_Float64, pixelSizeInBytes, rowSizeInBytes );
}
showData(data);
}
GDALClose(dataset);
The result of the code is intended to show twice the same data, however the
result is that the first shown image is correct and the second image is
horizontally distorted (compressed) as 1/4 of the good width.
Could anyone, please explain me whether I missed something ?
Cheers,
Victor
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
