Ok ho risolto con grass 7 su windows, sono riuscito a caricare correttamente il file senza creare valori nulli utilizzando r3.in.ascii.
utilizzando r3.in.xyz invece non risce a caricare correttamente tutti i layer: r3.in.xyz.py --overwrite input=C:\Documents and Settings\ospite\Desktop\file_input_grass\r3.in.xyz\punti_non_interpolati_12880_CSV1.txt output=grid_40x23_xyz separator=; value_column=4 Processing horizontal slice 1 of 14 [0,0.999999999999999) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 2760 punti nella regione. Processing horizontal slice 2 of 14 [1,2) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 3 of 14 [2,3) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 1840 punti nella regione. Processing horizontal slice 4 of 14 [3,4) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 5 of 14 [4,5) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 6 of 14 [5,6) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 7 of 14 [6,7) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 8 of 14 [7,8) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 9 of 14 [8,9) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 0 punti nella regione. Processing horizontal slice 10 of 14 [9,10) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 11 of 14 [10,11) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 0 punti nella regione. Processing horizontal slice 12 of 14 [11,12) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. Processing horizontal slice 13 of 14 [12,13) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 0 punti nella regione. Processing horizontal slice 14 of 14 [13,14) ... Lettura dei dati in corso ... Scrittura della mappa in corso ... r.in.xyz completo. Tovati 920 punti nella regione. raster map(s) available in mapset <PERMANENT>: Creazione della mappa raster 3D Done. 3D raster map <grid_40x23_xyz> created. Credo sia dovuto al fatto che la profondità dei layer aumenta in maniera non lineare. infatti nel primo layer ha individuato 2760 punti (920x3) invece di 920. credo abbia raggruppato i primi tre layer del mio dataset di input in un unico layer. allo stesso modo il terzo layer ha 1840 punti (929x2), e di conseguenza i layer 9, 11, 13 hanno 0 punti. comunque con r3.in.ascii questo non si verifica. ciao grazie ancora Il 28 novembre 2012 13:45, Markus Neteler <nete...@osgeo.org> ha scritto: > Ciao Matteo, > > 2012/11/27 Matteo Paolini <matteo.paolini...@gmail.com>: >> Ciao Markus >> >> dunque; con grass 6.4.2 (stand alone)su windows, > > per i voxel (e velocità ecc) consiglio di usare GRASS 7: > > http://wingrass.fsv.cvut.cz/grass70/ > > Anche se non ancora rilasciato ufficialmente ci sono tanti > vantaggi. > >> tramite r3.in.ascii >> carico un file di testo costiutito da 12880 punti con griglia 40 x 23 >> x 14. Con r3.stats ho notato che nel file ci sono dei valori nulli >> circa 2000. Con r3.out.vtk visualizzo il volume con paraview ma come >> si vede nello screenshot in basso c'è una zona blu con limiti definiti >> con valori nulli (-10000). perchè si sono creati questi valori nulli? > > Se il problema persiste in GRASS 7 possiamo guardare, così non > so quale sia il problema. Su GRASS 7 i voxel sono stati migliorati > ma probabilmente e un problema della formattazione dei dati originali. > Meglio r3.in.xyz allora. > >> Come faccio con grass a sostiutire i valori nulli con dei valori che >> siano interpolati correttamente a partire dai valori vicini? ho visto >> che c'è r3.null ma posso inserire solo un valore manualmente > > Con interpolazione (oppure r3.mapcalc): > > http://grass.osgeo.org/grass70/manuals/topic_voxel.html > -> v.vol.rst: Interpolates point data to a 3D raster map using > regularized spline with tension (RST) algorithm. > > ciao > Markus > > _______________________________________________ > 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