Tim Michelsen wrote: > Hello, > taking up on the post on accessing raster data from python [1] I'd like > to ask for your advice on writing arrays (numpy arrays) to python > raster/vector. > > I have a array of data of which I know the four boundary coordinates. > The array is a numpy array. > > Is there any example on how I can save this to a geoformat like Arc/Info > ASCII Grid or GeoTiff (including the world file)? > > How can I save a array of the form > xi yi zi > to a gdal raster? > > Are there also examples or tutorials for vector data input/output with > python? > I know that matplotlib basemap has a shapefile interface and can plot on > basemaps. > > Thanks and kind regards, > Tim > > [1]: http://permalink.gmane.org/gmane.comp.python.gis/15 >
Hi Tim, I don't work with rasters these days and am not familiar with the new Numpy interoperability features of GDAL and OGR. You'll need to search the GDAL site and mailing list for examples. Some Numpy arrays can be adapted to Shapely geometries. See http://gispython.org/shapely/manual.html#numpy-array-interface and then converted to WKT or WKB (and thereby directly into PostGIS if you want), but Shapely is not really a geodata package. Sean _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
