[GRASS-user] v.what.rast (large file)
Massimo Di Stefano
massimodisasha at yahoo.it
Mon Apr 9 04:06:56 EDT 2007
hi ,
i've a large point vector file (ascii 160 mb)
i've loaded it into grass using v.in.ascii : cat file.xyz |
v.in.ascii -r out=rst_k fs=" " col="x double, y double, z double"
then i added a column to the vector table using : v.db.addcol
map=rst_k layer=1 'columns=kval integer'
v.info:
Type of Map: Vector (level: 2)
Number of points: 1460445
Map is 3D: 0
Number of dblinks: 1
Projection: Transverse Mercator (zone
0) |
| N: 4507802.498 S:
4494627.516 |
| E: 2531867.500 W: 2523932.502
r.info :
Data Type: CELL
| Rows: 2635
| Columns: 1587
| Total Cells: 4181745
| Projection: Transverse Mercator (zone 0)
| N: 4507802.5 S: 4494627.5 Res: 5
| E: 2531867.5 W: 2523932.5 Res: 5
| Range of data: min = 1 max = 6
my targhet is to add the raster value to the kval column of my vector
points,
to do so i'm tring to use :
v.what.rast vector=rst_k raster=k6 layer=1 column=kval
i'm tring it on this region :
g.region:
projection: 99 (Transverse Mercator)
zone: 0
datum: rome40
ellipsoid: international
north: 4507802.498
south: 4494627.498
west: 2523932.502
east: 2531867.502
nsres: 5
ewres: 5
rows: 2635
cols: 1587
cells: 4181745
...after 24 hours at 100% cpu ... nothing is done ,
i force to kill the process :-(
but if i try it on a smoll area it works fine.
have you any suggestion about?
can i divide my area in many subregion
and then patch the results ?
thanks
Massimo
More information about the grass-user
mailing list