Hi all,
I'm tring to read a raster map from my grassData.
I'm using this code
import os
import grass.script.setup as gsetup
from grass.pygrass import raster
gisbase = os.environ['GISBASE']
gisdbase = os.path.abspath("/home/myuser/grassData")
location = 'mylocation'
mapset = 'mymapset'
gsetup.init(gisbase, gisdbase, location, mapset)
elev=raster.RasterRow('elevation')
elev.open()
for row in elev:
print row
But I get all values with nan, where is the error?
If I do `elev.info()` I get this result
elevation@
rows: 4502
cols: 5554
north: 45.0743663458 south: 45.0500262069 nsres:5.40651685671e-06
east: 7.67457811722 west: 7.6445503225 ewres:5.40651687473e-06
range: 233.050003052, 310.049987793
proj: 3
Tnx
Lorenzo
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
