<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Nov 22, 2017 at 1:02 PM, Farwell, Michael <span dir="ltr"><<a href="mailto:MFARWE01@harris.com" target="_blank">MFARWE01@harris.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">





<div lang="EN-US">
<div class="gmail-m_-261676721297304404WordSection1">
<p class="MsoNormal">Hey All,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Using PDAL 1.5, I’m trying to decompress a LAZ file and then re-project it into UTM using the following pipeline:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">{<u></u><u></u></p>
<p class="MsoNormal">  "pipeline":[<u></u><u></u></p>
<p class="MsoNormal">    {<u></u><u></u></p>
<p class="MsoNormal">      "type":"readers.las",<u></u><u></u></p>
<p class="MsoNormal">      "compression":"laszip",<u></u><u></u></p>
<p class="MsoNormal">      "filename": "/tmp/in.laz"<u></u><u></u></p>
<p class="MsoNormal">    },<u></u><u></u></p>
<p class="MsoNormal">    {<u></u><u></u></p>
<p class="MsoNormal">      "type":"filters.reprojection",<u></u><u></u></p>
<p class="MsoNormal">      "out_srs":"EPSG:32619"<u></u><u></u></p>
<p class="MsoNormal">    },<u></u><u></u></p>
<p class="MsoNormal">    {<u></u><u></u></p>
<p class="MsoNormal">      "type":"writers.las",<u></u><u></u></p>
<p class="MsoNormal">      "filename": "/tmp/out.las",<u></u><u></u></p>
<p class="MsoNormal">     "minor_version": 2<u></u><u></u></p>
<p class="MsoNormal">    }<u></u><u></u></p>
<p class="MsoNormal">  ]<u></u><u></u></p>
<p class="MsoNormal">}<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Since the minor_version is “2”, PDAL is using the GeoKeyDirectoryTag  to store the georeferencing information.  Unfortunately, PDAL is writing out an empty GeoDoubleParamsTag.<u></u><u></u></p>
<p class="MsoNormal"><u></u></p></div></div></blockquote><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">...</span></div><div><span style="font-size:12.8px"><br></span></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span style="font-size:12.8px">Is there any way to prevent PDAL from writing out this VLR if it’s empty?</span></blockquote><div><br></div><div>Short answer: no.  I checked the standard and I don't see anything suggesting this behavior is invalid, so if you have other software that doesn't accept this input, perhaps it should be fixed.<br></div><div><br></div><div>That said, the change to not write empty GeoDoubles and GeoAscii VLRs is trivial. Perhaps others have thoughts on this?</div><div><br></div></div>-- <br><div class="gmail_signature">Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</div></div>