使用GDAL C API编写Raster

Question

我正试图了解GDAL C API。 我有一些包含int，float或double值的数组，当我尝试使用GDALRasterIO将它们转储到tiff文件时，我会得到带有0值的垂直波段的栅格（我怀疑未初始化的值...

这是一个最小的例子：

#include "gdal.h"

void create_new_tiff(char * DstFilename, int nx, int ny);

int main()
{
char  filename[10] = "test.tif";
int nx, ny;

nx = 600;
ny = 500;

create_new_tiff(filename, nx, ny);

return 0;
}

create_new_tiff函数：

void create_new_tiff(char * DstFilename, int nx, int ny)
{

const char *pszFormat = "GTiff";
int test[nx*ny];
GDALDatasetH hDstDS;
GDALRasterBandH hBand;
char **papszOptions = NULL;

GDALAllRegister();
GDALDriverH hDriver = GDALGetDriverByName( pszFormat );

hDstDS = GDALCreate( hDriver, DstFilename, nx, ny, 1, GDT_Int16, papszOptions);
hBand = GDALGetRasterBand( hDstDS, 1 );

for(int i=0; i<nx*ny; i++) test[i] = 4;

CPLErr Err = GDALRasterIO( hBand, GF_Write, 0, 0, nx, ny, test, nx, ny, GDT_Int16, 0, 0);
GDALClose( hDstDS );
}

看起来我混淆了数据类型，但无法弄清楚出了什么问题......

由于

PS：以下是QGis中栅格的样子：

Answer 1

错误来自于我假设“int”实际上默认为2个字节，这是错误的。 int的大小取决于实现。见is int by default long int in C?

见kyle评论