<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif">Thanks for all your help and suggestions.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">After a short break, I relooked at the issue today and found it to be a case of mistaken spatial reference system.</div><div class="gmail_default" style="font-family:tahoma,sans-serif">The grid and the las data had different srs which meant there were points falling outside of the<span style="font-family:arial,sans-serif"> moturiki 1953 region.</span><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><span style="font-family:arial,sans-serif">In effect creating infinte values.</span></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><span style="font-family:arial,sans-serif"><br></span></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 22, 2019 at 11:48 PM Andrew Bell <<a href="mailto:andrew.bell.ia@gmail.com">andrew.bell.ia@gmail.com</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"><div dir="ltr">This is a good thought, Thejus, but that issue only comes up when using the --json option with PDAL translate.  Also, even if options were getting dropped/replaced, changing the scale won't ever get you Inf values.  They almost always come from strange reprojections.<div><br></div><div>If anyone is interested, this is the issue that Thejus was referring to: <a href="https://github.com/PDAL/PDAL/issues/2538" target="_blank">https://github.com/PDAL/PDAL/issues/2538</a> .  I haven't had a chance to fix it yet.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 22, 2019 at 7:33 AM thejus kambi <<a href="mailto:thejus.kambi@gmail.com" target="_blank">thejus.kambi@gmail.com</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"><div dir="ltr"><div>Hello Andrew,</div><div><br></div><div>For one of my queries, you had mentioned that the readers and writers can/will be replaced by command line parameters. Do you think something like that is happening?</div><div>The query was for "translate" command and not pipeline.<br></div><div><br></div><div>Warm regards,</div><div>Thejus<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 22, 2019 at 6:02 AM Andrew Bell <<a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</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"><div dir="ltr"><div dir="ltr"><div><br></div><div>Hi Rebecca,</div><div><br></div><div>I'm assuming that you can write the file without the reprojection filter.  If that is indeed the case, then Adam is right and some point is being reprojected to an infinite Z value (points that are unable to be transformed are dropped).  I'm not sure what you would like in the end, but you can write the output to some format you can examine (text, for example), to determine the points that are problematic.  There could also be some problem with the proj specification you've provided, but I'm not familiar enough with those details to be of much help.</div><div><br></div><div>Best,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 22, 2019 at 1:36 AM adam steer <<a href="mailto:adam.d.steer@gmail.com" target="_blank">adam.d.steer@gmail.com</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"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="arial, sans-serif">Hi Rebecca</font></div><div class="gmail_quote"><div><font face="arial, sans-serif"> </font></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><font face="arial, sans-serif">I ran lasinfo on the las file and get the following result:<br></font></div><div><font face="arial, sans-serif">Z                3696      58127<br></font></div></div></blockquote><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Ah those are unscaled values - so scaled to metres it would be 36.96 and 581.27?</font></div><div><font face="arial, sans-serif"><br></font></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><font face="arial, sans-serif">What would be the best way to determine if the geoid is introducing any divide-by-zeros?</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I loaded the gtx into postgis and these are the results from using the ST_SummaryStats:</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">'count' 'sum' 'mean'  'stddev'        'min'   'max'<br>14364 | 4215.86899040639 | 0.293502435979281 | 0.0288004424645927 | 0.17399999499321 | 0.49099999666214<br></font></div></div></blockquote><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">That’s a cool way to do it! I wondered if there were null or nodata values in the gtx file for some of your las points, but it appears not…</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">My last easy check idea is to make sure all the las points sit within the moturiki 1953 region.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I don’t feel like I’m being super helpful, sorry - just walking down a list of things that might produce an ‘-Infinity'….</font></div><div><br></div><div><font face="arial, sans-serif">Regards,</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Adam</font></div><div><font face="arial, sans-serif"> </font></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div style="font-family:tahoma,sans-serif"></div><div style="font-family:tahoma,sans-serif"><br></div><div style="font-family:tahoma,sans-serif">Cheers</div><div style="font-family:tahoma,sans-serif"><br></div><div style="font-family:tahoma,sans-serif">Rebecca</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 22, 2019 at 3:50 PM adam steer <<a href="mailto:adam.d.steer@gmail.com" target="_blank">adam.d.steer@gmail.com</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"><div dir="ltr">Hi Rebecca<div><br></div><div>What are the min and max Z values for your laz file? does it have a sneaky infinite value in there?</div><div><br></div><div>…and is your geoid grid introducing any divide-by-zeros to create un-handleable values?</div><div><br></div><div>Cheers</div><div><br></div><div>Adam</div></div></blockquote></div></blockquote></div><div dir="ltr" class="gmail-m_135176191207528863gmail-m_6620279518001276780gmail-m_-5718771113507098226gmail-m_-4092078182299936422gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_135176191207528863gmail-m_6620279518001276780gmail-m_-5718771113507098226gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div></div>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_135176191207528863gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</blockquote></div>