Thanks a lot, Roger, it worked like a charm, just fiddling a little with parameters and than run it in python. It's really magic, I struggled with this question for two days before writing to the list. I think that instead of making another script you could add it as an option to gdal_contour.py and use it alternatively with -l (default) or -p (output to polygons). By the way, another enhancement for gdal_contour could be adding an options to cut the output with a line or polygon layer (say the shoreline or other boundary). c
c > > On Wed, Jul 3, 2013 at 8:04 PM, Roger Veciana i Rovira <[email protected] > > wrote: > >> Hi, I was just working on that. >> I think that you need a script like gdal_contour, but creatig polygons >> instead of lines. The algorithm used is Marching Squares. >> You can see the code, without explanation, here: >> https://github.com/rveciana/geoexamples/tree/master/python/raster_isobands >> In some weeks I'll add some more info in a README file and transform the >> file into a working script, not just a function. But you can use it just >> changing the file name and the intervals. >> >> Roger >> >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
