[Openev-discuss] writing a GDALDataset to file

[Carolyn J. <Car...@ve...>]

OK, I give, this is driving me nuts. I am trying to create a GDALDataset 
in Python, fit it out with metadata and raster data, and write it to a 
file. So simple!

I am checking to see that the geotiff driver supports Create() (it 
does), creating the dataset with Create(), and then running the 
following code, but nothing happens after I call Create(), at which time 
the file has size 0.

    outformat = 'GTiff'
    out_driver = gdal.GetDriverByName(outformat)
    metadata = out_driver.GetMetadata()
    print metadata
    if metadata.has_key(gdal.DCAP_CREATE) \
           and metadata[gdal.DCAP_CREATE] == 'YES':
        print 'Driver %s supports Create() method.' % outformat
    else:
        print 'Driver %s does not support Create() method.' % outformat
    dd_obj.demdiff = out_driver.Create(dd_obj.ddfilename, xsize, ysize, 
1, gdal.GDT_Float32)
    print "demdiff object: ", dd_obj.demdiff
    dd_obj.demdiff.SetDescription("DEM Difference produced by DEM Fusion 
Prototype")
    dd_obj.demdiff.SetGeoTransform(std_geomatrix)
    dd_obj.demdiff.SetProjection(std_srs.ExportToWkt())
    dd_obj.demdiff.GetRasterBand(1).WriteArray(ddchunk)

So what is wrong with this?
:) Carolyn

-- 
_____________________________________
Carolyn P Johnston
Sr. Research Engineer, Vexcel Corporation
1690 38th Street, Boulder, CO 80301
ph: 303-583-0273
fax: 303-583-0246
http://www.vexcel.com