[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