<div dir="ltr">Hi,<div><br></div><div>I'm not sure that there's enough information here for me to help you.  You'll have to be specific about how the counts are different than you expect, how you're creating your tiles and reassembling them.  Keep in mind that points outside the provided bounds aren't considered when making the raster.  This typically means that some points that would be considered as contributing to a cell will get discarded when that cell is on the edge of the raster.  You're probably better off making the bounds one cell larger than you need (depending on the radius argument), and then trimming those edge cells away with GDAL or some such.  This is behavior that could be changed, but I believe that when I wrote the code, I was trying to keep it consistent with points2grid, which was the software that PDAL originally used to create rasters.</div><div><br></div><div>Hope that helps,</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 13, 2018 at 8:15 AM, Lars <span dir="ltr"><<a href="mailto:laasunde@hotmail.com" target="_blank">laasunde@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
Hello,
<p align="LEFT" style="font-size:12pt;margin-top:0px;margin-bottom:0px" dir="LTR">
<br>
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"></p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR">We have created a pipeline using pgpointcloud reader and gdal writer (see attached file).</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"></p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"><br>
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR">Using the pipeline to create a full raster for an area works great. However if we tile the same area we are seeing some 'artifacts' along the edge of the raster (band 5, Count). The count
 is different along the edges of the tiled raster compared with a full raster (see attached screenshot taken from QGis showing band 5 of rasters).</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"><br>
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"></p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR">We use the 'where' option in the reader.pgpointcloud to provide points that is writer.gdal.bounds + margin.
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"></p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"><br>
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR">We would like to know why these 'artifacts' occur and how to prevent them from happening?</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"></p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR"><br>
</p>
<p align="LEFT" style="margin-top:0px;margin-bottom:0px" dir="LTR">kind regards, Lars</p>
<br>
</div>
</div>

<br>______________________________<wbr>_________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/pdal</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</div>