<div dir="ltr"><div class="gmail_default" style="font-family:arial,sans-serif">Miguel,</div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif">I can't seem to find this in the docs, but I think what you want is "MERGE_ALG" set to ADD instead of the default REPLACE.</div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif"><a href="https://github.com/OSGeo/gdal/blob/master/gdal/alg/gdalrasterize.cpp#L1165">https://github.com/OSGeo/gdal/blob/master/gdal/alg/gdalrasterize.cpp#L1165</a><br></div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif">This doesn't address the buffering (I guess you might still want to buffer the whole set), but it should take care of the additive step.  </div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif">I actually implemented this as part of a programming test for my current job - to build heatmaps of satellite coverage. </div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,sans-serif">Best regards,</div><div class="gmail_default" style="font-family:arial,sans-serif">Frank</div><div class="gmail_default" style="font-family:arial,sans-serif"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, May 4, 2021 at 1:30 PM Miguel A. Manso <<a href="mailto:m.manso@upm.es">m.manso@upm.es</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all<br>
<br>
I have a multiline layer type (tracking) and I want to make a heat map <br>
(densities) on a buffer of them, so that in those places where several <br>
buffers overlap the density is added as many times as lines overlap.<br>
<br>
Now I'm doing it with a python script that reads from a shp the lines, <br>
generates a buffer and rasterize the polygon filling it with a given <br>
value with gda.RasterizeLayer(). I read the image as an array and with <br>
numpy I add it to another array (accumulating). When I finish I convert <br>
it to an image and save it in a geotiff.<br>
<br>
This procedure is relatively slow. Especially if you have to do it for <br>
many lines.<br>
<br>
My question is along the lines of whether the RasterizeLayer function <br>
itself could do this process in a more efficient way, making the burning <br>
of the objects not an OR operation but an arithmetic sum...<br>
<br>
Any suggestions would be welcome.<br>
<br>
Regards<br>
Miguel A. Manso<br>
<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><font face="monospace">---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>light and sound - activate the windows | +1 650-701-7823<br>and watch the world go round - Rush    | Geospatial Software Developer</font></div></div></div></div></div>