Ciao Giulia. Non ti preoccupare, l'impatto con GRASS è sempre un po' doloroso inizialmente! ;) Provo a rispondere ad alcune domande, anche se essendo di corsa non potrò farlo in maniera approfondita.
1- HD esterno. Per GRASS è indifferente che tu utilizzi un HD interno o esterno, l'unica cosa è evitare l'utilizzo di nomi di directory e file con spazi, perché questo può dare dei problemi. Se non è il tuo caso, ci potrebbero essere problemi di permessi di scrittura sull'HD esterno o in specifiche directory. Dacci qualche elemento in più per aiutarti... 2 - Non so che tipo di elaborazioni dovrai fare, ma spero tanto non sia necessario avere l'intero mosaico di DEM in memoria! Mi spiego, definire una regione anche enorme non dà alcun problema, il problema eventualmente è elaborare un intero, gigantesco, raster. In questo caso servono accorgimenti specifici, che potremo indicarti se fosse questa la tua necessità. Per definire i parametri della regione, ti conviene importare una sola tile (come hai fatto) e poi usare il comando "g.region rast=il_nome_della_tua_tile", così anche la risoluzione viene impostata a quella del raster. Il DEM NASA è in proiezione WGS84, con risoluzione di 30' d'arco, mentre la risoluzione in metri non è costante su tutto il globo (mediamente 1' d'arco alle nostre latitudini sono una trentina di metri). Quindi, se imposti la location con una proiezione WGS84, lascia che g.region si vada a pescare da sé la risoluzione corretta... Per avere invece il bounding box di tutte le tue tile, potresti usare l'utility gdal "gdalbuildvrt" [1]. E' un'utility molto comoda, perché senza richiedere tanta memoria, ti crea una sorta di "mosaico virtuale" di tutte le tue tile, senza fisicamente crearne uno, e ti calcola anche il boundingbox completo. Es.: >gdalbuildvrt dem_mosaico.vrt dem_tiles/*.tif In uscita otterrai un file .vrt, che è un file xml. La cosa bella è che in GRASS (tramite r.in.gdal) puoi importarti direttamente questo file (dem_mosaico.vrt) e se usi l'opzione "-e" (extend location) la regione della tua location verrà estesa in modo da comprendere tutto il tuo mosaico. 3 - FileMaker non lo conosco. Hai dati con attributi geometrici puntuali? Puoi tranquillamente esportare in dbf, ma la parte geometrica come pensavi di gestirla? Anche qui, dacci qualche informazione in più. Scusa la frettolosità, magari in una prossima email cercherò di spiegare meglio... Giovanni [1] http://www.gdal.org/gdalbuildvrt.html Il giorno 10 novembre 2011 13:28, Giulia De Nobili <xeil...@gmail.com> ha scritto: > Salve a tutti, > sono una studentessa di archeologia e sto impostando il GIS per la mia > tesi di laurea. Ho sempre usato ArcGIS ma ho deciso di passare a > GRASS... in più non ho mai affrontato un progetto così grosso e sono > seriamente in difficoltà! Spero che qualcuno con molta pazienza possa > consigliarmi nelle fasi iniziali, poi dovrei riuscire a cavarmela! > > Quello che devo fare è un'analisi della distribuzione dei siti > archeologici nell'area della Trans-caucasia. E' un'area piuttosto > ampia e le nostre informazioni sono decisamente lacunose quindi ho > deciso di utilizzare i DEM della NASA fatti sulle immagini ASTER. Ne > ho scaricati un centinaio che dovrebbero coprire grosso modo tutta > l'area di mio interesse. > > Prima di tutto ho il problema che ero abituata a lavorare con dati > contenuti su un HD esterno perchè ho a disposizione solo il mio > portatile. Ma non riesco a creare la location sull'HD esterno, non so > se è un problema del mio HD o di GRASS che effettivamente non permette > di fare questa cosa... > > Ho creato la location sul computer e ho utilizzato una delle tiles del > DEM per impostare i parametri della regione. Però in questo modo non > ha settato la risoluzione che so essere un parametro importante per le > analisi successive. Io volevo impostare la risoluzione perchè > coincidesse con quella del DEM. Nel pdf allegato alle immagini che ho > scaricato c'è scritto che la risoluzione orizzontale è di 75m. Allora > ho provato a fare g.region res=75, ho provato a fare un g.region ma a > quel punto ottengo: > > projection: 3 (Latitude-Longitude) > zone: 0 > datum: wgs84 > ellipsoid: wgs84 > north: 44:22:40N > south: 28:23:55N > west: 33:15:40E > east: 51:20:42E > nsres: 15:58:45 > ewres: 18:05:02 > rows: 1 > cols: 1 > cells: 1 > > e mi suona male che celle, colonne e righe siano =1!!! Allora ho > provato a fare dei controlli e ho realizzato che l'unità di misura in > questo momento è in gradi e quindi suppongo che io abbia impostato la > risoluzione per equivalere a 75° e non 75m, è così???? ma a questo > punto non so che pesci pigliare... come faccio a cambiare l'unità di > misura da gradi a metri? mi sarebbe utile anche in seguito per > calcolare le distanze tra i siti... non vorrei dover ogni volta fare > le conversioni. > > Secondo voi mi conviene unire tutte le tiles del DEM in un unico file? > tanto avrò sempre bisogno di lavorare sull'intera superficie...ma ci > mette una vita a caricare tutto! > Tra l'altro io per andare meglio lavoro con l'interfaccia grafica (gis > manager) ma non mi lascia selezionare più di un file per volta, tipo: > se dovessi fare r.patch c'è un modo per selezionare in un'unica volta > tutte le tiles del DEM? > > Dopodichè io ho a disposizione un database fatto con FileMaker di > tutti i siti archeologici che conosciamo in zona, uno dei campo > contiene anche le coordinate geografiche. Pensavo di esportarlo > in .dbf, aggiungere una voce ID (non so se devo numerarli da 0 o da 1) > per ottenere lo shape file... l'avrò fatto mille volte in ArcGis ma > non sono totalmente sicura di quale sia la procedura con GRASS, mi > aiutate? > > grazie mille > > Giulia > > > > > > _______________________________________________ > GRASS-Italia mailing list > GRASS-Italia@listserv.unipr.it > http://listserv.unipr.it/mailman/listinfo/grass-italia >
_______________________________________________ GRASS-Italia mailing list GRASS-Italia@listserv.unipr.it http://listserv.unipr.it/mailman/listinfo/grass-italia
