Geografic-Information-System/swig/python/examples/rasteraccess.py

#!/usr/bin/python

# run within GRASS Spearfish session

import os, sys
import swig.grass as grasslib
import swig.raster as grassrast

if not os.environ.has_key("GISBASE"):
    print "You must be in GRASS GIS to run this program."
    sys.exit(1)

if len(sys.argv)==2:
  input = sys.argv[1]
else:
  input = raw_input("Raster Map Name? ")

# initialize
grasslib.G_gisinit('')

# find map in search path
mapset = grasslib.G_find_cell2(input,'')

# determine the inputmap type (CELL/FCELL/DCELL) */
data_type = grasslib.G_raster_map_type(input, mapset)

infd = grasslib.G_open_cell_old(input, mapset)
inrast = grasslib.G_allocate_raster_buf(data_type)

rown=0
while 1:
    myrow = grasslib.G_get_raster_row(infd, inrast, rown, data_type)
    print rown,myrow[0:10]
    rown = rown+1
    if rown==476:break

grasslib.G_close_cell(inrast)
grasslib.G_free(cell)