<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body>
<br>
<div class="moz-forward-container">
<table class="moz-email-headers-table" cellspacing="0"
cellpadding="0" border="0">
<tbody>
<tr>
<th valign="BASELINE" nowrap="nowrap" align="RIGHT">Subject:
</th>
<td>Re: [pdal] z-units conversion</td>
</tr>
<tr>
<th valign="BASELINE" nowrap="nowrap" align="RIGHT">Date: </th>
<td>Thu, 19 Mar 2020 11:09:42 -0600</td>
</tr>
<tr>
<th valign="BASELINE" nowrap="nowrap" align="RIGHT">From: </th>
<td>Paul White <a class="moz-txt-link-rfc2396E" href="mailto:paul@dgi.com"><paul@dgi.com></a></td>
</tr>
<tr>
<th valign="BASELINE" nowrap="nowrap" align="RIGHT">To: </th>
<td>Andrew Bell <a class="moz-txt-link-rfc2396E" href="mailto:andrew.bell.ia@gmail.com"><andrew.bell.ia@gmail.com></a></td>
</tr>
</tbody>
</table>
<br>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<div class="moz-cite-prefix">Hi Andrew,<br>
<br>
Thanks for your reply.<br>
<br>
On 3/19/2020 10:31 AM, Andrew Bell wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CACJ51z13CZcXT+Hiec9ix+0ZhuCui2B+Cgr7O-jYJQ_fXrRGcQ@mail.gmail.com">
<meta http-equiv="content-type" content="text/html;
charset=UTF-8">
<div dir="ltr">Without seeing your pipeline it's hard to
answer. </div>
</blockquote>
Here's my .json file:<br>
<br>
<font size="+1"><tt>[</tt><tt><br>
</tt><tt> "Job520660_36076_92_01.laz",</tt><tt><br>
</tt><tt> {</tt><tt><br>
</tt><tt> "type":"filters.colorization",</tt><tt><br>
</tt><tt> "raster":"google_wms_satellite_1478.tif"</tt><tt><br>
</tt><tt> },</tt><tt><br>
</tt><tt> "pdal_colorization_Job520660_36076_92_01.laz"</tt><tt><br>
</tt><tt>]</tt></font><br>
<br>
and I execute it like this in a Anaconda3 Prompt window:<br>
<br>
<font size="+1"><tt>pdal pipeline colorize.json</tt></font><br>
<br>
<br>
<blockquote type="cite"
cite="mid:CACJ51z13CZcXT+Hiec9ix+0ZhuCui2B+Cgr7O-jYJQ_fXrRGcQ@mail.gmail.com">
<div dir="ltr">You can assign any SRS you like at output with
the a_srs option, </div>
</blockquote>
Yes, I'm aware of that, using 'spatialreference',
'filters.reprojection', 'in_srs', 'out_srs', etc. <br>
In this use case, the SRS of the input .laz is "EPSG: 32618". In
the case of EPSG:32618,<br>
it is a UTM projection and by definition, the x/y units are in
meters. All well and good...<br>
<br>
However, SRS definitions do not assign a unit value to Z. <br>
<br>
My input .laz file (shown above) has the z-units in feet
(actually, they are identified<br>
in the .laz as "Linear_Foot_US_Survey"). Here is the result of
running lasinfo on the<br>
original .laz:<br>
<br>
<font size="+1"><tt> Keyed_Information:</tt><tt><br>
</tt><tt> GTModelTypeGeoKey (Short,1): ModelTypeProjected</tt><tt><br>
</tt><tt> ProjectedCSTypeGeoKey (Short,1):
PCS_WGS84_UTM_zone_18N</tt><tt><br>
</tt><tt> ProjLinearUnitsGeoKey (Short,1): Linear_Meter</tt><tt><br>
</tt><tt> VerticalCSTypeGeoKey (Short,1): Unknown-3855</tt><tt><br>
</tt><tt> <b>VerticalUnitsGeoKey (Short,1):
Linear_Foot_US_Survey</b></tt><tt><br>
</tt><tt> End_Of_Keys.</tt></font><br>
<br>
After the 'filters.colorization' step, the z-units are
changed/renamed to "Linear_Meter":<br>
<br>
<font size="+1"><tt> Keyed_Information:</tt><tt><br>
</tt><tt> GTModelTypeGeoKey (Short,1): ModelTypeProjected</tt><tt><br>
</tt><tt> GTRasterTypeGeoKey (Short,1): RasterPixelIsArea</tt><tt><br>
</tt><tt> GTCitationGeoKey (Ascii,8): "unknown"</tt><tt><br>
</tt><tt> GeogCitationGeoKey (Ascii,7): "WGS 84"</tt><tt><br>
</tt><tt> GeogAngularUnitsGeoKey (Short,1): Angular_Degree</tt><tt><br>
</tt><tt> ProjectedCSTypeGeoKey (Short,1):
PCS_WGS84_UTM_zone_18N</tt><tt><br>
</tt><tt> ProjLinearUnitsGeoKey (Short,1): Linear_Meter</tt><tt><br>
</tt><tt> VerticalCSTypeGeoKey (Short,1): Unknown-3855</tt><tt><br>
</tt><tt> VerticalCitationGeoKey (Ascii,15): "EGM2008
height"</tt><tt><br>
</tt><tt> VerticalDatumGeoKey (Short,1): Unknown-1027</tt><tt><br>
</tt><tt> <b>VerticalUnitsGeoKey (Short,1): Linear_Meter</b></tt><tt><br>
</tt><tt> End_Of_Keys.</tt></font><br>
<br>
<br>
As you can see, the 'VerticalUnitsGeoKey' has been changed to
meters.<br>
<br>
This renders the 'colorized' .laz incompatible with all of my
other data sources<br>
(DEM's, DSM's, shapefiles and other 3D data types, etc.) in a
3D/4D co-visualized<br>
application environment.<br>
<br>
I'm looking for a PDAL filter to either leave the
'VerticalUnitsGeoKey' alone (output <br>
.laz same as input .laz) or set/name it to 'Linear_Foot_US_Survey'
or even just plain 'Feet'.<br>
<br>
Any suggestions would be most helpful.<br>
<br>
Thanks for your time,<br>
Paul<br>
<br>
<br>
<br>
<blockquote type="cite"
cite="mid:CACJ51z13CZcXT+Hiec9ix+0ZhuCui2B+Cgr7O-jYJQ_fXrRGcQ@mail.gmail.com">
<div dir="ltr">but I'm not sure what else might be going on.</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Wed, Mar 18, 2020 at 6:08
PM Paul White <<a href="mailto:paul@dgi.com"
moz-do-not-send="true">paul@dgi.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">I have an input .las
which is UTM 18 N-- x/y units are in meters and <br>
z-units are in 'Linear_Foot_US_Survey'.<br>
When process the input .las through any PDAL filters, the
output z-units <br>
are renamed to 'Linear_Meter'.<br>
How can I designate the z-units to remain in feet?<br>
<br>
_______________________________________________<br>
pdal mailing list<br>
<a href="mailto:pdal@lists.osgeo.org" target="_blank"
moz-do-not-send="true">pdal@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/pdal"
rel="noreferrer" target="_blank" moz-do-not-send="true">https://lists.osgeo.org/mailman/listinfo/pdal</a></blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr" class="gmail_signature">Andrew Bell<br>
<a href="mailto:andrew.bell.ia@gmail.com" target="_blank"
moz-do-not-send="true">andrew.bell.ia@gmail.com</a></div>
<span class="BEGIN-ANTISPAM-VOTING-LINKS"></span>
<div class="BEGIN-ANTISPAM-VOTING-LINKS" style="display: inline">
<!-- BEGIN-ANTISPAM-VOTING-LINKS -->
<div><span style="font-size: medium; color: black; font-style:
normal; font-weight: normal; background-color: white;
display: inline">
<hr> NOTE: This message was trained as non-spam. If this
is wrong, please correct the training as soon as possible.
<br>
<a rel="nofollow" target="canit_note"
href="https://rio.dgi.com:8443/b.php?c=s&i=012ggw7io&m=d965d3c93268"
moz-do-not-send="true">Spam</a><br>
<a rel="nofollow" target="canit_note"
href="https://rio.dgi.com:8443/b.php?c=n&i=012ggw7io&m=d965d3c93268"
moz-do-not-send="true">Not spam</a><br>
<a rel="nofollow" target="canit_note"
href="https://rio.dgi.com:8443/b.php?c=f&i=012ggw7io&m=d965d3c93268"
moz-do-not-send="true">Forget previous vote</a><br>
</span></div>
<!-- END-ANTISPAM-VOTING-LINKS --> </div>
<span class="END-ANTISPAM-VOTING-LINKS"></span> </blockquote>
<br>
</div>
</body>
</html>