[gdal-dev] Writing raster - 0 in every even column

[Goo Creations]

Hi all

I have a strange problem when writing to a GDAL dataset. I'm currently doing
the following:

*GDALDataset *mInputDataset = (GDALDataset*) GDALOpen("in.tif",
GA_ReadOnly);
GDALDataType type = mInputDataset->GetRasterBand(1)->GetRasterDataType();
int mWidth = mInputDataset->GetRasterXSize();
int mHeight = mInputDataset->GetRasterYSize();

GDALDataset *mOutputDataset = mInputDataset->GetDriver()->Create("out.tif",
mWidth, mHeight, mBands, type, 0);

int *data = <<get data for entire band>>
mOutputDataset->GetRasterBand(1)->RasterIO(GF_Write, 0, 0, aligner.width(),
aligner.height(), data, aligner.width(), aligner.height(), type, 0, 0);
mOutputDataset->FlushCache();
*

The input image provided has a datatype of UINT16. When I hardcode all
"type" variables in the above example to GDT_Int32, then the image is
correctly created, but when I leave the example as is (hence with GDT_UInt16
assigned to "type"), then I get a strange problem that every second column
has a value of 0. So every odd column has the correct values in, but every
even column has values of 0. I've checked the data array, and all the values
in there are correct.

Does anyone know why this is happening?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110603/b0627a9e/attachment-0001.html