Dear Grass community,
I am trying to use r.walk from python and the output file seems to be created
correctly at the end of the python run.
However:
- while python is running I get the following error (but code keeps on running)
: ERROR: Unable to create file
<C:\A_GRID\PM\Grass\GrassDB/Location_L3/PERMANENT/.tmp/1824.0>
- All values of the output raster seem to be equal to 0.0 when I try to access
them through a Grass garray. But I open the raster in Grass GUI all, I can see
raster values were correctly populated for the whole raster. I get a nice cost
surface.
Here is the core of my code : (coord is an array containing coordinates of
point in the reference system of the current region)
a = garray.array()
grass.run_command("g.remove", rast=cost_output)
grass.run_command("r.walk", elevation="mnt_l3_toulouse@PERMANENT",
friction="mnt_l3_toulouse@PERMANENT", output=cost_output,
coordinate=[p1[0],p1[1]],_lambda=0)
a.read(cost_output)
for p2 in coord:
x = int((p2[0]-w)/ewres)
y = int((n-p2[1])/nsres)
print a[y,x]
Thanks in advance!
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
