<div dir="ltr">Micheal,<div>The problem is the z values of all of the point values are have been changed. I had created a raster from the mean value of the ground points from the pdal reprojected data. I then went back to NOAA's site and selected a small subset to compare and had their engine deliver the points in Albers Meters Conus NAVD88. I imported the points to create a raster in the same way and then subtraced the Z value of the pdal reprojected points from the NOAA reprojected points. The height values vary by 32-36 m ( see below) .</div><div><br></div><div>That was why I was asking if I had missed a step that ensured that the Z values stayed at the Geoid12A/12B/NAVD88 values.</div><div><br></div><div><div>Tue Nov 14 09:45:55 2017) </div><div><a href="http://r.info">r.info</a> map=noaa_pdal_vert_diff_m@dnewcomb </div><div> +----------------------------------------------------------------------------+</div><div> | Map: noaa_pdal_vert_diff_m@dnewcom Date: Tue Nov 14 09:36:19 2017 |</div><div> | Mapset: dnewcomb Login of Creator: dnewcomb |</div><div> | Location: seva_albers |</div><div> | DataBase: /data2/grass70_data |</div><div> | Title: noaa_pdal_vert_diff_m |</div><div> | Timestamp: none |</div><div> |----------------------------------------------------------------------------|</div><div> | |</div><div> | Type of Map: raster Number of Categories: 0 |</div><div> | Data Type: FCELL |</div><div> | Rows: 12465 |</div><div> | Columns: 20411 |</div><div> | Total Cells: 254423115 |</div><div> | Projection: NAD83 / Conus Albers |</div><div> | N: 1687158.47979 S: 1668162.01936 Res: 1.52398399 |</div><div> | E: 1653024.11274 W: 1621917.96408 Res: 1.52398945 |</div><div> | Range of data: min = 32.17395 max = 36.06071 |</div><div> | |</div><div> | Data Description: |</div><div> | generated by r.mapcalc |</div><div> | |</div><div> | Comments: |</div><div> | 2635_mean_ground_1_524_test_from_noaa_meters@dnewcomb - |</div><div> | 2635_mean_ground_1_524@PERMANENT |</div><div> | |</div><div> +----------------------------------------------------------------------------+</div></div><div><br></div><div><br></div><div><br></div><div>Doug</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 13, 2017 at 8:20 PM, Michael Smith <span dir="ltr"><<a href="mailto:michael.smith.erdc@gmail.com" target="_blank">michael.smith.erdc@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="white" lang="EN-US" link="#0563C1" vlink="#954F72"><div class="m_-6789161839749526003WordSection1"><p class="MsoNormal">Doug,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">It looks like perhaps the header metadata is not correct, which is what I think you are seeing from lasinfo. When you inspected the Decimal Degrees file it shows<u></u><u></u></p><span class=""><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"> "bbox":<u></u><u></u></p><p class="MsoNormal"> {<u></u><u></u></p><p class="MsoNormal"> "maxx": -76.2527221,<u></u><u></u></p><p class="MsoNormal"> "maxy": 37.9068695,<u></u><u></u></p><p class="MsoNormal"> <span style="background:yellow">"maxz": -5.777103023,</span><u></u><u></u></p><p class="MsoNormal"> "minx": -76.2571767,<u></u><u></u></p><p class="MsoNormal"> "miny": 37.9028129,<u></u><u></u></p><p class="MsoNormal"> "<span style="background:yellow">minz": -35.71281856</span><u></u><u></u></p><p class="MsoNormal"> },<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p></span><p class="MsoNormal">which match your albers output. <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><span style="font-size:10.5pt;color:black">-- <u></u><u></u></span></p><div><p class="MsoNormal"><span style="font-size:10.5pt;color:black">Michael Smith<u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="font-size:10.5pt;color:black">Remote Sensing/GIS Center<u></u><u></u></span></p></div><p class="MsoNormal"><span style="font-size:10.5pt;color:black">US Army Corps of Engineers</span><u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in"><p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">pdal <<a href="mailto:pdal-bounces@lists.osgeo.org" target="_blank">pdal-bounces@lists.osgeo.org</a>> on behalf of Doug Newcomb <<a href="mailto:newcomb.ttg@gmail.com" target="_blank">newcomb.ttg@gmail.com</a>><br><b>Date: </b>Monday, November 13, 2017 at 8:14 PM<br><b>To: </b>pdal <<a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a>><br><b>Subject: </b>[pdal] Z value change on reprojection in pdal 1.6<u></u><u></u></span></p></div><div><div class="h5"><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Ok,<u></u><u></u></p><div><p class="MsoNormal">Finally diving into pdal more deeply. I installed the docker instance of pdal 1.6 and tried a bit of command line reprojection:<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">I downloaded some data from the NOAA website and re-projected from Decimal Degrees to Albers Meters Conus ( EPSG:5070) <u></u><u></u></p></div><div><p class="MsoNormal">pdal translate /data/20100503_LAS_S23_1836_<wbr>10.laz /data/albers_20100503_LAS_S23_<wbr>1836_10.laz -f filters.reprojection --filters.reprojection.out_<wbr>srs="EPSG:5070"<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> For the input data , lasinfo gives :<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> Minimum and Maximum Attributes (min,max)<u></u><u></u></p></div><div><p class="MsoNormal">------------------------------<wbr>---------------------------<u></u><u></u></p></div><div><p class="MsoNormal"> Min X, Y, Z: <span class="m_-6789161839749526003gmail-apple-tab-span"> </span>-76.2571767, 37.9028129, 0.024<u></u><u></u></p></div><div><p class="MsoNormal"> Max X, Y, Z: <span class="m_-6789161839749526003gmail-apple-tab-span"> </span>-76.2527221, 37.9068695, 29.966<u></u><u></u></p></div><div><p class="MsoNormal"> Bounding Box:<span class="m_-6789161839749526003gmail-apple-tab-span"> </span>-76.2571767, 37.9028129, -76.2527221, 37.9068695<u></u><u></u></p></div></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">for the output laz file, lasinfo gives:<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><div><p class="MsoNormal"> Minimum and Maximum Attributes (min,max)<u></u><u></u></p></div><div><p class="MsoNormal">------------------------------<wbr>---------------------------<u></u><u></u></p></div><div><p class="MsoNormal"> Min X, Y, Z: <span class="m_-6789161839749526003gmail-apple-tab-span"> </span>1707096.44, 1829911.69, -35.72<u></u><u></u></p></div><div><p class="MsoNormal"> Max X, Y, Z: <span class="m_-6789161839749526003gmail-apple-tab-span"> </span>1707557.20, 1830350.12, -5.77<u></u><u></u></p></div><div><p class="MsoNormal"> Bounding Box:<span class="m_-6789161839749526003gmail-apple-tab-span"> </span>1707096.44, 1829911.69, 1707557.20, 1830350.12<u></u><u></u></p></div></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">pdal info for the input file gives:<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><div><p class="MsoNormal">{<u></u><u></u></p></div><div><p class="MsoNormal"> "filename": "\/data\/20100503_LAS_S23_<wbr>1836_10.laz",<u></u><u></u></p></div><div><p class="MsoNormal"> "pdal_version": "1.6.0 (git-version: 94ad83)",<u></u><u></u></p></div><div><p class="MsoNormal"> "stats":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "EPSG:4326":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "maxx": -76.2527221,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxy": 37.9068695,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxz": -5.777103023,<u></u><u></u></p></div><div><p class="MsoNormal"> "minx": -76.2571767,<u></u><u></u></p></div><div><p class="MsoNormal"> "miny": 37.9028129,<u></u><u></u></p></div><div><p class="MsoNormal"> "minz": -35.71281856<u></u><u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "boundary": {<u></u><u></u></p></div><div><p class="MsoNormal"> "coordinates" : [<u></u><u></u></p></div><div><p class="MsoNormal"> [<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.257176700000002, 37.902812900000001 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.257176700000002, 37.906869499999999 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.2527221, 37.906869499999999 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.2527221, 37.902812900000001 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.257176700000002, 37.902812900000001 ]<u></u><u></u></p></div><div><p class="MsoNormal"> ]<u></u><u></u></p></div><div><p class="MsoNormal"> ],<u></u><u></u></p></div><div><p class="MsoNormal"> "type" : "Polygon"<u></u><u></u></p></div><div><p class="MsoNormal">}<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "native":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "maxx": -76.2527221,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxy": 37.9068695,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxz": 29.966,<u></u><u></u></p></div><div><p class="MsoNormal"> "minx": -76.2571767,<u></u><u></u></p></div><div><p class="MsoNormal"> "miny": 37.9028129,<u></u><u></u></p></div><div><p class="MsoNormal"> "minz": 0.024<u></u><u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "boundary": {<u></u><u></u></p></div></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">pdal info on the output laz file gives:<u></u><u></u></p></div><div><div><p class="MsoNormal">{<u></u><u></u></p></div><div><p class="MsoNormal"> "filename": "\/data\/albers_20100503_LAS_<wbr>S23_1836_10.laz",<u></u><u></u></p></div><div><p class="MsoNormal"> "pdal_version": "1.6.0 (git-version: 94ad83)",<u></u><u></u></p></div><div><p class="MsoNormal"> "stats":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "EPSG:4326":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "maxx": -76.25185156,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxy": 37.90686951,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxz": -5.77,<u></u><u></u></p></div><div><p class="MsoNormal"> "minx": -76.25806657,<u></u><u></u></p></div><div><p class="MsoNormal"> "miny": 37.9021939,<u></u><u></u></p></div><div><p class="MsoNormal"> "minz": -35.72<u></u><u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "boundary": {<u></u><u></u></p></div><div><p class="MsoNormal"> "coordinates" : [<u></u><u></u></p></div><div><p class="MsoNormal"> [<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.258066569999997, 37.903041880000004 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.257028300000002, 37.90686951 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.251851560000006, 37.90602148 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.252890100000002, 37.9021939 ],<u></u><u></u></p></div><div><p class="MsoNormal"> [ -76.258066569999997, 37.903041880000004 ]<u></u><u></u></p></div><div><p class="MsoNormal"> ]<u></u><u></u></p></div><div><p class="MsoNormal"> ],<u></u><u></u></p></div><div><p class="MsoNormal"> "type" : "Polygon"<u></u><u></u></p></div><div><p class="MsoNormal">}<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "native":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "bbox":<u></u><u></u></p></div><div><p class="MsoNormal"> {<u></u><u></u></p></div><div><p class="MsoNormal"> "maxx": 1707557.2,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxy": 1830350.12,<u></u><u></u></p></div><div><p class="MsoNormal"> "maxz": -5.77,<u></u><u></u></p></div><div><p class="MsoNormal"> "minx": 1707096.44,<u></u><u></u></p></div><div><p class="MsoNormal"> "miny": 1829911.69,<u></u><u></u></p></div><div><p class="MsoNormal"> "minz": -35.72<u></u><u></u></p></div><div><p class="MsoNormal"> },<u></u><u></u></p></div><div><p class="MsoNormal"> "boundary": {<u></u><u></u></p></div></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Am I missing a switch in the translate command?<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Doug<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div><p class="MsoNormal">______________________________<wbr>_________________ pdal mailing list <a href="mailto:pdal@lists.osgeo.org" target="_blank">pdal@lists.osgeo.org</a> <a href="https://lists.osgeo.org/mailman/listinfo/pdal" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/pdal</a><u></u><u></u></p></div></div>
</blockquote></div><br></div>