Many thanks for your suggestions. I will give it try and let you know which one works.
Cheers Sab > I used StarSpan to do something similar (but with less vector feature: > 6000). It worked pretty well. Maybee worth a try? > > http://starspan.projects.atlas.ca.gov/doku/doku.php > > Daniel > > On Fri, Nov 26, 2010 at 10:41 AM, Markus Metz > <[email protected]> wrote: >> On Fri, Nov 26, 2010 at 12:32 PM, <[email protected]> wrote: >>> Dear list >>> >>> I am trying to run v.rast.stats for a large vector file. It contains >>> around 700000 buildings representing buildings. >>> >>> here is my g.region info >>> g.region -p >>> projection: 99 (Transverse Mercator) >>> zone: 0 >>> datum: osgb36 >>> ellipsoid: airy >>> north: 183010.45829976 >>> south: 153516.19024768 >>> west: 513122.49174945 >>> east: 555953.95182415 >>> nsres: 4.99987592 >>> ewres: 5.00017045 >>> rows: 5899 >>> cols: 8566 >>> cells: 50530834 >>> >>> THe problem I am having is the process is extremely slow. It has done >>> 150 >>> features out of 700000 areas after 2 hours. >>> I need just the average raster value within each building (feature). >>> My questions are: >>> - Is there something I can do to speed up the process? >>> - Is there a way just to get the average value rather than the whole >>> set >>> of statistical parameters (I am not even using -e option) >>> >> You could try v.rast.stats2 from grass-addons, it should complete in >> less than 1 hour. When using v.rast.stats2 and grass 6.4, you also >> need r,univar.zonal from grass-addons (not needed for grass 6.5). >> >> Markus M >> _______________________________________________ >> grass-user mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/grass-user >> > _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
