<html>
  <head>

    <meta http-equiv="content-type" content="text/html; charset=UTF-8">
  </head>
  <body>
    <p> </p>
    <div class="moz-forward-container">
      <p>Greetings,<br>
        <br>
        My goal is to reproject/align a third party GRS80 point cloud
        data source to an existing UTM10 WGS84 (EPSG 32610)
        georeferenced point cloud file/"ground truth".<br>
        <br>
        Running Lasinfo on the third party data source (to be
        corrected/reprojected/aligned to Ground Truth) shows the
        following:<br>
        <br>
        <br>
        <i>// 'OGC Transformation Record'</i><i><br>
        </i><i>    WKT OGC COORDINATE SYSTEM:</i><i><br>
        </i><i>    COMPD_CS["NAD83(2011) / Conus Albers + NAVD88
          height",PROJCS["NAD83(2011) / Conus
          Albers",GEOGCS["NAD83(2011)",DATUM["NAD83 (National Spatial
          Reference System 2011)",SPHEROID["GRS
1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","1116"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","6318"]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",29.5],PARAMETER["standard_parallel_2",45.5],PARAMETER["latitude_of_center",23],PARAMETER["longitude_of_center",-96],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["meter",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","6350"]],VERT_CS["NAVD88
          height",VERT_DATUM["North American Vertical Datum
1988",2005,AUTHORITY["EPSG","5103"]],UNIT["meter",1,AUTHORITY["EPSG","9001"]],AXIS["Up",UP],AUTHORITY["EPSG","5703"]]]
          //</i><br>
        <br>
        <br>
        I am able to successfully reproject the file in PDAL so that the
        relative 2D alignment in Easting/Northing is corrected - using a
        command like this:<br>
        <br>
        <i><span style="color: rgb(29, 28, 29); font-family: Slack-Lato,
            Slack-Fractions, appleLogo, sans-serif; font-size: 15px;
            font-variant-ligatures: common-ligatures; font-variant-caps:
            normal; font-weight: 400; letter-spacing: normal;
            text-align: left; text-indent: 0px; text-transform: none;
            white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(255,
            255, 255); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;">                        
            //pdal translate filename.laz filename_utm.laz reprojection
            --filters.reprojection.out_srs="EPSG:32610"//</span></i><br>
        <br>
        <br>
        <b>Reviewing the output, this shows proper 2D Grid alignment
          between ground truth and third party data, but a delta of ~30M
          in elevation between the ground planes. </b><b>This needs to
          be rectified.<br>
        </b></p>
      <b>I've recently found this site describing a process to reproject
        a .las file in respect to various vertical datums.</b><br>
      <p> <br>
        <a moz-do-not-send="true"
          href="http://www.gadom.ski/2017/05/31/vdatum-with-pdal.html"
          class="moz-txt-link-freetext">http://www.gadom.ski/2017/05/31/vdatum-with-pdal.html</a><br>
        <br>
        <i><span style="color: rgb(17, 17, 17); font-family:
            -apple-system, BlinkMacSystemFont, "Segoe UI",
            Roboto, Helvetica, Arial, sans-serif, "Apple Color
            Emoji", "Segoe UI Emoji", "Segoe UI
            Symbol"; font-size: 16px; font-variant-ligatures:
            normal; font-variant-caps: normal; font-weight: 400;
            letter-spacing: normal; text-align: start; text-indent: 0px;
            text-transform: none; white-space: normal; word-spacing:
            0px; -webkit-text-stroke-width: 0px; background-color:
            rgb(253, 253, 253); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;">//In
            order to use<span> </span></span></i><i><a
            href="https://www.pdal.io/" style="color: rgb(42, 122, 226);
            text-decoration: none; font-family: -apple-system,
            BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
            Arial, sans-serif, "Apple Color Emoji",
            "Segoe UI Emoji", "Segoe UI Symbol";
            font-size: 16px; font-variant-ligatures: normal;
            font-variant-caps: normal; font-weight: 400; letter-spacing:
            normal; text-align: start; text-indent: 0px; text-transform:
            none; white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(253,
            253, 253);" moz-do-not-send="true">PDAL</a></i><i><span
            style="color: rgb(17, 17, 17); font-family: -apple-system,
            BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
            Arial, sans-serif, "Apple Color Emoji",
            "Segoe UI Emoji", "Segoe UI Symbol";
            font-size: 16px; font-variant-ligatures: normal;
            font-variant-caps: normal; font-weight: 400; letter-spacing:
            normal; text-align: start; text-indent: 0px; text-transform:
            none; white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(253,
            253, 253); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;"><span> </span>to
            convert vertical coordinates to a new vertical datum, we’ll
            need to fetch the appropriate geoid. PDAL uses<span> </span></span></i><i><a
            href="http://proj4.org/" style="color: rgb(42, 122, 226);
            text-decoration: none; font-family: -apple-system,
            BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
            Arial, sans-serif, "Apple Color Emoji",
            "Segoe UI Emoji", "Segoe UI Symbol";
            font-size: 16px; font-variant-ligatures: normal;
            font-variant-caps: normal; font-weight: 400; letter-spacing:
            normal; text-align: start; text-indent: 0px; text-transform:
            none; white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(253,
            253, 253);" moz-do-not-send="true">proj.4</a></i><i><span
            style="color: rgb(17, 17, 17); font-family: -apple-system,
            BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
            Arial, sans-serif, "Apple Color Emoji",
            "Segoe UI Emoji", "Segoe UI Symbol";
            font-size: 16px; font-variant-ligatures: normal;
            font-variant-caps: normal; font-weight: 400; letter-spacing:
            normal; text-align: start; text-indent: 0px; text-transform:
            none; white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(253,
            253, 253); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;"><span> </span>under
            the hood, and proj.4 accepts geoid’s in NOAA VDatum’s<span> </span></span></i><i><code
            class="language-plaintext highlighter-rouge"
            style="font-size: 15px; border: 1px solid rgb(232, 232,
            232); border-radius: 3px; background-color: rgb(238, 238,
            255); padding: 1px 5px; color: rgb(17, 17, 17);
            font-variant-ligatures: normal; font-variant-caps: normal;
            font-weight: 400; letter-spacing: normal; text-align: start;
            text-indent: 0px; text-transform: none; white-space: normal;
            word-spacing: 0px; -webkit-text-stroke-width: 0px;
            text-decoration-thickness: initial; text-decoration-style:
            initial; text-decoration-color: initial;">gtx</code></i><i><span
            style="color: rgb(17, 17, 17); font-family: -apple-system,
            BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
            Arial, sans-serif, "Apple Color Emoji",
            "Segoe UI Emoji", "Segoe UI Symbol";
            font-size: 16px; font-variant-ligatures: normal;
            font-variant-caps: normal; font-weight: 400; letter-spacing:
            normal; text-align: start; text-indent: 0px; text-transform:
            none; white-space: normal; word-spacing: 0px;
            -webkit-text-stroke-width: 0px; background-color: rgb(253,
            253, 253); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;"><span> </span>file
            format. We can fetch the 12B geoid from the<span> </span></span></i><i><a
            href="https://vdatum.noaa.gov/download.php" style="color:
            rgb(42, 122, 226); text-decoration: none; font-family:
            -apple-system, BlinkMacSystemFont, "Segoe UI",
            Roboto, Helvetica, Arial, sans-serif, "Apple Color
            Emoji", "Segoe UI Emoji", "Segoe UI
            Symbol"; font-size: 16px; font-variant-ligatures:
            normal; font-variant-caps: normal; font-weight: 400;
            letter-spacing: normal; text-align: start; text-indent: 0px;
            text-transform: none; white-space: normal; word-spacing:
            0px; -webkit-text-stroke-width: 0px; background-color:
            rgb(253, 253, 253);" moz-do-not-send="true">VDatum download
            website</a></i><i><span style="color: rgb(17, 17, 17);
            font-family: -apple-system, BlinkMacSystemFont, "Segoe
            UI", Roboto, Helvetica, Arial, sans-serif, "Apple
            Color Emoji", "Segoe UI Emoji", "Segoe
            UI Symbol"; font-size: 16px; font-variant-ligatures:
            normal; font-variant-caps: normal; font-weight: 400;
            letter-spacing: normal; text-align: start; text-indent: 0px;
            text-transform: none; white-space: normal; word-spacing:
            0px; -webkit-text-stroke-width: 0px; background-color:
            rgb(253, 253, 253); text-decoration-thickness: initial;
            text-decoration-style: initial; text-decoration-color:
            initial; display: inline !important; float: none;">. Inside
            the 12B zipfile there is a file called<span> </span></span></i><i><code
            class="language-plaintext highlighter-rouge"
            style="font-size: 15px; border: 1px solid rgb(232, 232,
            232); border-radius: 3px; background-color: rgb(238, 238,
            255); padding: 1px 5px; color: rgb(17, 17, 17);
            font-variant-ligatures: normal; font-variant-caps: normal;
            font-weight: 400; letter-spacing: normal; text-align: start;
            text-indent: 0px; text-transform: none; white-space: normal;
            word-spacing: 0px; -webkit-text-stroke-width: 0px;
            text-decoration-thickness: initial; text-decoration-style:
            initial; text-decoration-color: initial;">g2012b_conus.gtx</code></i><span
          style="color: rgb(17, 17, 17); font-family: -apple-system,
          BlinkMacSystemFont, "Segoe UI", Roboto, Helvetica,
          Arial, sans-serif, "Apple Color Emoji", "Segoe
          UI Emoji", "Segoe UI Symbol"; font-size: 16px;
          font-style: normal; font-variant-ligatures: normal;
          font-variant-caps: normal; font-weight: 400; letter-spacing:
          normal; orphans: 2; text-align: start; text-indent: 0px;
          text-transform: none; white-space: normal; widows: 2;
          word-spacing: 0px; -webkit-text-stroke-width: 0px;
          background-color: rgb(253, 253, 253);
          text-decoration-thickness: initial; text-decoration-style:
          initial; text-decoration-color: initial; display: inline
          !important; float: none;"><i>; this is the geoid for the
            Continental United States (CONUS). We’ll use this file to
            convert our data.//</i><br>
          <br>
          <br>
        </span>I've downloaded the VDatum collection as stated - That
        said, I've been unable to get the process to work properly.<br>
        <br>
        <i>//pdal translate filename.laz filename_utm.laz reprojection
          --filters.reprojection.out_srs="+init=EPSG:32610
          +geoidgrids=~/Downloads/vdatum/g2012b_conus.gtx"
          --writers.las.a_srs=EPSG:32610+4326</i><i><br>
        </i><i>(pdal translate Error) GDAL failure (1) SetCompoundCS()
          fails, vertical component is not VERT_CS.</i><i><br>
        </i><i>(pdal translate writers.las Error) GDAL failure (1) PROJ:
          proj_crs_get_coordinate_system: Object is not a SingleCRS//</i><br>
        <br>
        <br>
        Does anyone have any Tips/Pointers as to how to resolve/debug
        this issue? <br>
        <br>
        Thanks.</p>
    </div>
  </body>
</html>