[GRASS-user] v.to.rast

[patrick s.]

Markus

I tested the approach you proposed, but seems to be very slow for large 
datasets: I am processing ~1000.000 points of observation on a 25m-Grid 
across Switzerland, where only few cells have multiple points inside. 
While /v.to.rast/ takes a few minutes for conversion, the combination 
"/v.out.ascii input=pt output=- column=VAL | r.in.xyz input=- z=4 
output=pt method=sum/" has been running for multiple hours and is still 
in progress.

Is there any alternative to run this conversion? I need it to run 
neighborhood analyses as kernel density with population field, i.e. 
/r.neighbors in=pt out=pt_dens -c size=300 meth=sum/. Maybe there is an 
approach to directly use the vector data and avoid conversion- something 
as /v.neighbors meth=sum/?

Thanks for you help,
Patrick


On 24.07.2015 04:13, Markus Neteler wrote:
> On Thu, Jul 23, 2015 at 11:00 AM, patrick s. <patrick_gis at gmx.net> wrote:
>> Dear all
>>
>> I am puzzled on the behavior of v.to.rast. When several points fall into one
>> gridcell, the raster seems to get one value but not the sum of these. Is
>> there a way to sum these up instead?
> Yes. I have added a related example here:
>
> http://grass.osgeo.org/grass70/manuals/v.to.rast.html#convert-vector-points-to-raster-with-raster-cell-binning
>
> (while it does not really fit to that manual page it is expected
> there. Perhaps we need to really enhance v.to.rast to do such a job
> right away).
>
> HTH
> Markus

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20150730/e58b8048/attachment.html>