[GRASS-user] r.basin: Completes on Some Basins, Fails on Others

[Margherita Di Leo]

Hi Micha,

The next addition might be how feed these crosspoints X and Y straight into
> the easting and northing of r.basins...
>
>
you mean something like (untested):

#!/usr/bin/env python

fin = open("crosspoints","r")
fout = open("commands","w")
linesout = (line.rstrip().split() for line in fin)
cmd = 'r.basin map=dem prefix=bas%s easting=%s northing=%s
threshold=10000\n'
fout.writelines(cmd % (n, easting, northing) for n,(easting,northing) in
enumerate(linesout))

?

ciao madi


-- 
Dr. Margherita Di Leo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20120607/f2864ca7/attachment-0001.html>