Hello list, I'm trying to rasterize some 3d-lines programtically from python. Using gdal 1.11.2 from the osgeo4w64 environment. Seems to work fine, except that I get a z-offset of 255 (one byte) in my output raster. Works fine if I use gdal_rasterize with the -3d option.
The following python code should reproduce the problem: from osgeo import ogr,gdal print gdal.__version__ georef=[0,1.0,0,10.5,0,-1.0] m_drv_ogr=ogr.GetDriverByName("Memory") line_ds = m_drv_ogr.CreateDataSource( "dummy") layer = line_ds.CreateLayer( "lines", None, ogr.wkbLineString25D) layerdefn=layer.GetLayerDefn() line=ogr.Geometry(ogr.wkbLineString25D) line.AddPoint(0,5.5,0) line.AddPoint(10,5.5,10) feature=ogr.Feature(layerdefn) feature.SetGeometry(line) res=layer.CreateFeature(feature) layer.ResetReading() m_drv_gdal=gdal.GetDriverByName("MEM") raster_ds=m_drv_gdal.Create("dummy",10,10,1,gdal.GDT_Float64) raster_ds.SetGeoTransform(georef) ok=gdal.RasterizeLayer(raster_ds,[1], layer, options=["BURN_VALUE_FROM=Z"]) A=raster_ds.ReadAsArray() print A Will output: [[ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 255. 256. 257. 258. 259. 260. 261. 262. 263. 264.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.] [ 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]] Cheers, Simon Kokkendorff, Dansish Geodata Agency
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev