grazie Paolo ho scaricato lo script...poi faccio un po' di prove "Non sicuro che creare un DTM in formato shape point sia una cosa saggia, ma immagino tu sappia cosa stai facendo."
sul fatto che sappia cosa sto facendo...di solito è raro cmq i file raster sono dei dati da satellite di precision farming con risuluzione 1m...ci sono diversi dati tra cui anche DTM e aspect...non sono io che devo usarli ma chi mi ha chiesto di "capire cos'erano questi .adf" aveva bisogno di una tabella da usare con SPSS. I raster riguardano una regione molto piccola e se io faccio tutte le trasformazioni con risoluzioni 1m dovrei ottenere il dato del pixel nel punto del file vettoriale shape. ...non so se questo vuol dire sapere cosa sto facendo...però mi è sembrata una buona soluzione... certo potrei dirgli di usare i moduli di statistica raster che ci sono in grass ma non so ne esattamente cosa deve fare... grazie ancora francesco > Message: 1 > Date: Sun, 27 Apr 2008 23:20:36 +0200 > From: Paolo Zatelli <[EMAIL PROTECTED]> > Subject: Re: [GRASS-Italia] file aux e annessi - script > To: [EMAIL PROTECTED], Mailing list utenti italiani di > GRASS <[EMAIL PROTECTED]> > Message-ID: <[EMAIL PROTECTED]> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > francesco sdraiati ha scritto: > >> Grazie Paolo per l'aiuto >> vado quindi a spiegarmi....credo sia abbastanza semplice ma il problema >> è che mi manca la conoscenza sulla sintassi >> >> Ho creato una location in WGS84 32N con cod EPSG 32632 >> >> LOCATION_NAME='br_pedo'; >> MAPSET='aux'; >> GISDBASE='/home/francesco/grass'; >> >> i passaggi che faccio sono i seguenti >> >> r.in.gdal -o 'input=/media/............/w001001.adf' output=aspect_prova >> >> r.to.vect input=aspect_prova output=aspect_prova feature=point >> >> v.out.ogr input=aspect type=point dsn=aspect olayer=aspect layer=1 >> format=ESRI_Shapefile >> >> uno dei problemi è che per ogni file che devo importare esiste una >> cartella omonima la quale contiene diversi file ma l'unico che mi >> interessa è sempre per tutti il file w001001.adf >> Per esempio: se voglio importare l'aspect devo andare nella cartella >> aspect e selezionare il file w001001.adf >> se voglio importare il file DTM devo andare nella cartella DTM e >> selezionare il file w001001.adf >> ....ecc >> >> Quindi importo tanti file w001001.adf (raster) e voglio costruire un >> file .shp aspect, uno DTM (vector chiaramente)....ecc >> >> spero di essermi espresso con chiarezza >> grazie >> Francesco >> >> > > Ciao Francesco, > credo di avere capito cosa ti serve. Ho creato lo script raster2shape.sh > che ho inserito sulla nostra pagina del software > http://www.ing.unitn.it/~grass/software.html > alla voce "GRASS scripts". > > Questo script importa tutti i file che trova nelle sottodirectory (di > qualunque profondita') a partire dalla directory da cui si lancia che > corrispondono ad un pattern dato (ad es. *.adf oppure w001001.adf) e > crea uno shape corrispondente (di tipo point) con il nome della > directory in cui e' stato trovato il file. > Il trucco e' tutto qui: > map_name=`dirname "$FILE"| rev | cut -d '/' -f1 | rev` > > > Dove sono ora non ho file adf per fare la prova ma dovrebbe andare, la > selezione dei nomi funziona correttamente, mi fido di avere usato i > comandi di GRASS con la sintassi giusta. > > *In ogni caso fai dei test!* > > Lo script imposta la regione sulla mappa importata (sia come estensione > che come risoluzione), se non e' quello che vuoi devi cambiare la riga > che contiene g.region secondo le tue necessita'. > > Non sicuro che creare un DTM in formato shape point sia una cosa saggia, > ma immagino tu sappia cosa stai facendo. > > Il messaggio e' arrivato in lista da alcuni giorni, ti rispondo solo ora > perche', contrariamente a quanto credono i miei studenti :) , non sono > sempre collegato in rete, soprattutto se ho finito il lavoro quando il > tempo regge e' probabile sia su qualche sentiero con la mia fida scale. > > Fammi sapere se funziona e soprattutto se risolve il problema. > > ciao Paolo > _______________________________________________ GRASS-Italia mailing list [email protected] http://listserv.unipr.it/mailman/listinfo/grass-italia
