[gdal-dev] 3D Shapefiles

[Opel Fahrer]

Hi,

I stumbled across an unexpected problem when I tried to write my own shapefile following the tutorial found here:

http://www.gdal.org/ogr/ogr_apitut.html

I want to store 3D coordinates, which should be possible I guess. However, it seems that doing something like

//========
OGRLayer *poLayer = poDS->CreateLayer( "points", NULL, wkbPoint, NULL );
OGRFeature *poFeature = OGRFeature::CreateFeature( poLayer->GetLayerDefn() );

OGRPoint pt(1, 2, 3);
poFeature->SetGeometry( &pt );
poLayer->CreateFeature( poFeature );
//========

does NOT create a 3D feature, but flattens to 2D. Whatever I tried, I could not figure out how to get around this. 

Any help is appreciated. I built my gdal libraries from the latest stable sources 1.6.2.

Regards