Jose Gomez-Dans wrote:
You are absolutely right. It's just that you have to call external program such as gdal_rasterize.
Ah -- I see what you mean by cumbersome now.
Is the rasterization stuff supposed to be available to the Python GDAL bindings?
That would be nice -- anyone know? I suppose looking at the gdal_rasterize code might help you figure it out.
Another option is to rasterize to an in-memory bitmap with another lib, like PIL or pyCairo, or mapnik, or maybe mapserver, though I don't know if you could keep it in memory with mapserver.
In the end, I ended up rasterising rather than testing each geometry. Blindingly fast ;)
I"m glad that worked out. -Chris -- Christopher Barker, Ph.D. Oceanographer Emergency Response Division NOAA/NOS/OR&R (206) 526-6959 voice 7600 Sand Point Way NE (206) 526-6329 fax Seattle, WA 98115 (206) 526-6317 main reception [EMAIL PROTECTED] _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
