<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Thanks Stefan<br>
<br>
It seems first of all v.out.ascii is slow. Following your proposaI,
the tmp-file has only 5300 observations after 9min15s. As I have
650.000 points this export will need approx. 18h! Can I speed this
up somehow? (Database-driver is dbf) <br>
<br>
Patrick<br>
<br>
P.S. g.region -up & v.info are given below.<br>
<br>
<br>
##g.region -up<br>
projection: 99 (Swiss. Obl. Mercator)<br>
zone: 0<br>
datum: ch1903<br>
ellipsoid: bessel<br>
north: 1297000<br>
south: 1074200<br>
west: 2484400<br>
east: 2834800<br>
nsres: 25<br>
ewres: 25<br>
rows: 8912<br>
cols: 14016<br>
cells: 124910592<br>
<br>
##v.info<br>
Number of points: 647957 Number of centroids:
0 |<br>
| Number of lines: 0 Number of boundaries:
0 |<br>
| Number of areas: 0 Number of islands:
0 |<br>
|
|<br>
| Map is 3D:
No |<br>
| Number of dblinks:
1 |<br>
|
|<br>
| Projection: Swiss. Obl.
Mercator |<br>
|
|<br>
| N: 1294822 S:
1075554 |<br>
| E: 2831184 W: 2486245 <br>
<br>
<br>
<div class="moz-cite-prefix">On 30.07.2015 13:59, Blumentrath,
Stefan wrote:<br>
</div>
<blockquote
cite="mid:c970dc642e40488bb5331efa96649b98@NINSRV23.nina.no"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;
color:black;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";
color:black;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;
color:black;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">Hi Patrick,<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">From my experience r.in.xyz is very fast.<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">How many points are you processing (v.info pt)
and what are your region settings (g.region –up)?<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">In order to find out where time is spend yo
could split upt the process into two steps:<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">v.out.ascii input=pt output=./tmp column=VAL<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">r.in.xyz input=./tmp z=4 output=pt method=sum<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">Cheers<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US">Stefan<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1
1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:windowtext"
lang="EN-US">From:</span></b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:windowtext"
lang="EN-US"> <a class="moz-txt-link-abbreviated" href="mailto:grass-user-bounces@lists.osgeo.org">grass-user-bounces@lists.osgeo.org</a>
[<a class="moz-txt-link-freetext" href="mailto:grass-user-bounces@lists.osgeo.org">mailto:grass-user-bounces@lists.osgeo.org</a>] <b>On
Behalf Of </b>patrick s.<br>
<b>Sent:</b> 30. juli 2015 10:18<br>
<b>To:</b> Markus Neteler<br>
<b>Cc:</b> GRASS user list<br>
<b>Subject:</b> Re: [GRASS-user] v.to.rast<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span
lang="EN-US">Markus<br>
<br>
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
<i>v.to.rast</i> takes a few minutes for conversion, the
combination "<i>v.out.ascii input=pt output=- column=VAL |
r.in.xyz input=- z=4 output=pt method=sum</i>" has been
running for multiple hours and is still in progress.
<br>
<br>
Is there any alternative to run this conversion? I need it
to run neighborhood analyses as kernel density with
population field, i.e.
<i>r.neighbors in=pt out=pt_dens -c size=300 meth=sum</i>.
Maybe there is an approach to directly use the vector data
and avoid conversion- something as
<i>v.neighbors meth=sum</i>? <br>
<br>
Thanks for you help,<br>
Patrick<br>
<br>
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span lang="EN-US">On 24.07.2015 04:13,
Markus Neteler wrote:<o:p></o:p></span></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<pre><span lang="EN-US">On Thu, Jul 23, 2015 at 11:00 AM, patrick s. </span><a moz-do-not-send="true" href="mailto:patrick_gis@gmx.net"><span lang="EN-US"><patrick_gis@gmx.net></span></a><span lang="EN-US"> wrote:<o:p></o:p></span></pre>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<pre><span lang="EN-US">Dear all<o:p></o:p></span></pre>
<pre><span lang="EN-US"><o:p> </o:p></span></pre>
<pre><span lang="EN-US">I am puzzled on the behavior of v.to.rast. When several points fall into one<o:p></o:p></span></pre>
<pre><span lang="EN-US">gridcell, the raster seems to get one value but not the sum of these. Is<o:p></o:p></span></pre>
<pre><span lang="EN-US">there a way to sum these up instead?<o:p></o:p></span></pre>
</blockquote>
<pre><span lang="EN-US"><o:p> </o:p></span></pre>
<pre><span lang="EN-US">Yes. I have added a related example here:<o:p></o:p></span></pre>
<pre><span lang="EN-US"><o:p> </o:p></span></pre>
<pre><a moz-do-not-send="true" href="http://grass.osgeo.org/grass70/manuals/v.to.rast.html#convert-vector-points-to-raster-with-raster-cell-binning"><span lang="EN-US">http://grass.osgeo.org/grass70/manuals/v.to.rast.html#convert-vector-points-to-raster-with-raster-cell-binning</span></a><span lang="EN-US"><o:p></o:p></span></pre>
<pre><span lang="EN-US"><o:p> </o:p></span></pre>
<pre><span lang="EN-US">(while it does not really fit to that manual page it is expected<o:p></o:p></span></pre>
<pre><span lang="EN-US">there. Perhaps we need to really enhance v.to.rast to do such a job<o:p></o:p></span></pre>
<pre>right away).<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre>HTH<o:p></o:p></pre>
<pre>Markus<o:p></o:p></pre>
</blockquote>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</blockquote>
<br>
</body>
</html>