<html>

  <head>

    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">

    <style type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>

  </head>

  <body bidimailui-charset-is-forced="true" text="#000000"

    bgcolor="#FFFFFF">

    <br>

    <div class="moz-cite-prefix">On 10/21/19 1:33 PM, Tais wrote:<br>

    </div>

    <blockquote type="cite"

      cite="mid:d581c296-ec9e-b2a1-7733-45722647afa2@ulb.ac.be">Hi all,

      <br>

      <br>

      I would like to compute a raster density map ("heat map") from

      vector points using v.kernel. The points represent house

      addresses. However, I would need to weight the points with a value

      stored in the attribute table (number of inhabitants). I saw on

      the module manual page that this functionality is a "known issue".

      <br>

      <br>

      I imagined a solution to by-passe the limitation: duplicate each

      point as many time as needed (according to the value stored in the

      attribute table) and use this new layer directly in v.kernel.

      However, I think it will definitely not be the most efficient way

      to do the trick. Do you have an alternative in mind ?

      <br>

      <br>

      Of course the best solution would be that someone in the

      development community would have the time and the kindness to

      implement this function directly in the module (I'm not skilled in

      C and thus cannot try myself).

      <br>

      <br>

    </blockquote>

    <p><br>

    </p>

    <p>Maybe not the answer you are looking for, but you can create a

      weighted kernel density in R. I can send an R script that is

      intended to be run from within a GRASS session. It accepts a point

      vector and attrib column, and creates a kernel density raster,

      weighted by the attribute values (using the R defaults for

      bandwidth and Gaussian kernel)</p>

    <p>Requires, of course, that R is installed, with the packages:

      spatstat, raster, rgrass7.<br>

    </p>

    <p><br>

    </p>

    <blockquote type="cite"

      cite="mid:d581c296-ec9e-b2a1-7733-45722647afa2@ulb.ac.be">NB:

      Sorry if I should have posted this on the developer mailing list.

      I hesitated but decided to post here finally.

      <br>

      <br>

      Best,

      <br>

      <br>

    </blockquote>

    <pre class="moz-signature" cols="72">-- 

Micha Silver

Ben Gurion Univ.

Sde Boker, Remote Sensing Lab

cell: +972-523-665918</pre>

  </body>

</html>