Salve a tutti
oggi mi sono cimentato per la prima volta con uno script in python per
automatizzare l'operazione di esportazione di molti raster dal mapset,
sottoforma di geotiff.
Per la precisione vorrei esportare tutti i raster che iniziano con la
lettera maiuscola S.
Sono partito da un pezzo di codice trovato in rete e l'ho modificato.
#!/usr/bin/env python
import grass.script as gscript
import sys
def main():
gscript.run_command('g.region', flags='p')
rastlist = gscript.read_command('g.list',type='rast')
rastlist1 = rastlist.split(':', 1)[1]
rastlist2 = rastlist1.split('-', 1)[0]
RastListClean=rastlist2.split()
folderOut = '/media/giacomo/Volume/Livorno/outraster/raster'
for rasterOut in RastListClean:
if rasterOut.startswith("S"):
folderFileOut = folderOut + "/" + rasterOut + '.tiff'
grass.run_command('r.out.gdal', input=rasterOut, output=folderFileOut,
format='GTiff')
if __name__ == '__main__':
main()
Ho ottenuto questo messaggio di errore:
File "/home/giacomo/Desktop/prova.py", line 9, in main
rastlist1 = rastlist.split(':', 1)[1]
IndexError: list index out of range
Non mi aspettavo certo che tutto funzionasse alla prima, non sono ne un
genio di python tantomeno di GRASS (purtroppo).
La cosa che mi ha lasciato perplesso è che non riesco a capire come
debuggare il mio codice.
Non trovo una shell pyhon interattiva all'interno di GRASS grazie alla
quale è possibile testare parti di codice (QGIS ad esempio ce l'ha).
Magari c'è (sicuramente) e io non l'ho ancora trovata.
Sapete darmi qualche dritta?
W GRASS!
Grazie
_______________________________________________
GRASS-Italia mailing list
GRASS-Italia@listserv.unipr.it
http://listserv.unipr.it/mailman/listinfo/grass-italia
