[Gfoss] [GRASS 7] Grass 7 e python script

[junkhead@funbox.it]

Ciao ho installato da snv i nuovi sorgenti, effetivamente lo
script di test è cambiato adesso è il seguente (per
chiarezza):

--- ini script 
#!/usr/bin/env python

"""
Run within GRASS session
Run this before starting python to append module search
path:

@code
export
PYTHONPATH=$PYTHONPATH:/usr/local/grass-7.0.svn/etc/python/
@endcode

Check with "import sys; sys.path"
or:

@code
sys.path.append("/usr/local/grass-7.0.svn/etc/python")
@endcode
"""

import os, sys
from grass.lib import grass
from grass.lib import 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
grass.G_gisinit('')

# find map in search path
mapset = grass.G_find_raster2(input, '')

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

infd = grassrast.Rast_open_old(input, mapset)
inrast = grassrast.Rast_allocate_buf(data_type)

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

grassrast.Rast_close(inrast)

--- fine script 

anche io ho il medesimo errore "TypeError: 'int' object is
unsubscriptable"

ma questo penso sia dovuto al fatto che la funzione
Rast_get_row restituiscae  1, 0 o -1 e non un array
(http://download.osgeo.org/grass/grass7_progman/raster_2get__row_8c.html#ec254ec7e611f98dd42974e8b7051475)
risolto questo segnalo un successivo errore nella righa
grassrast.Rast_close(inrast) che si risolve chiudendo la
giusta risorsa cioè grassrast.Rast_close(infd).

Premetto che sono un ex utente arcgis passato con estrema
gioia a grass perciò sono ancora molto inesperto su questo
nuovo universo.
Dopo un discreto successo in termini di produttivita lato
user desktop vorrei imparare a fare qualche script... nulla
di faraonico...ma qualcosa del tipo:

apri un raster 
crea un nuovo raster
leggi il valore della cella x,y e mettilo in a
leggi il valore della cella x-1,y e mettilo in b
se a>b scrivi 1 nella cella x,y dei un nuovo raster 

Mi trovo solo un pò perplesso da non aver trovato troppe
risorse con esempi per quanto rigurada la programmazione di
script qualcuno sa suggerirmi qualche risorsa?

grazie mille maurizio



----- Original Message -----
Da : Markus Neteler <neteler a osgeo.org>
A : "junkhead a funbox.it" <maurizio.melani a gislab.it>
Cc: gfoss a faunalia.it
Oggetto : Re: [Gfoss] [GRASS 7] Grass 7 e python script
Data : Sat, 10 Oct 2009 18:25:04 +0200

> 2009/10/10 junkhead a funbox.it <maurizio.melani a gislab.it>:
> > Salve ho installato grass 7 dal snv e devo dire che
> > funziona tutto all
> > 100%... a parte gli script python...
> >
> > Ad esempio se lancio lo script rasteraccess.py (fornito
> > nella
> > distribuzione) ho in uscita l'errore :
> >
> > File "./rasteraccess.py", line 40, in <module>
> >    data_type = grass.G_raster_map_type(input,
mapset)
> > AttributeError: 'module' object has no attribute
> > 'G_raster_map_type'
> >
> > Credo di aver fatto la giusta installazione visto che
> > alcune righe
> > funzionano... (tipo grass.G_find_cell2(input, '')..)
> 
> si, in GRASS 7 la rasterlib e la gislib sono finalmente
> state separate. La funzione si chiama ora
> G_find_raster2() e tutte G_raster_*() sono
> Rast_*().
> 
> > confrontando le reference di grass 6
> >
> (http://download.osgeo.org/grass/grass6_progman/opencell_8
> > c.html) con le
> > nuove di grass 7
> >
> (http://download.osgeo.org/grass/grass7_progman/opencell_8
> > c.html) ho
> > notato che la proprietà ha cambiato nome
> > G_raster_map_type in grass 7 si
> > invoca con Rast_map_type .... però provando ancora non
> > funziona...
> 
> Ho appena fatto (in SVN) gli aggiornamenti per riportarlo
> allo stato di GRASS 7.
> Mi da ancora un errore di "TypeError: 'int' object is
> unsubscriptable" ma spero che sia una banalità che non
ho
> visto.
> 
> Potresti aggiornare dal SVN e riprovare?
> Magari stasera Martin Landa da un'occhiata.
> 
> ciao
> Markus