[pdal] Z value change on reprojection in pdal 1.6

[Doug Newcomb]

Ok,
Finally diving into pdal more deeply.  I installed the docker  instance of
pdal 1.6 and tried a bit of command line reprojection:

I downloaded some data from the NOAA website and  re-projected from Decimal
Degrees to Albers Meters Conus ( EPSG:5070)
pdal translate /data/20100503_LAS_S23_1836_10.laz
/data/albers_20100503_LAS_S23_1836_10.laz -f filters.reprojection
 --filters.reprojection.out_srs="EPSG:5070"

 For the input data , lasinfo gives :


  Minimum and Maximum Attributes (min,max)
---------------------------------------------------------
  Min X, Y, Z: -76.2571767, 37.9028129, 0.024
  Max X, Y, Z: -76.2527221, 37.9068695, 29.966
  Bounding Box: -76.2571767, 37.9028129, -76.2527221, 37.9068695

for the output laz file, lasinfo gives:

  Minimum and Maximum Attributes (min,max)
---------------------------------------------------------
  Min X, Y, Z: 1707096.44, 1829911.69, -35.72
  Max X, Y, Z: 1707557.20, 1830350.12, -5.77
  Bounding Box: 1707096.44, 1829911.69, 1707557.20, 1830350.12


pdal info for the input file gives:


{
  "filename": "\/data\/20100503_LAS_S23_1836_10.laz",
  "pdal_version": "1.6.0 (git-version: 94ad83)",
  "stats":
  {
    "bbox":
    {
      "EPSG:4326":
      {
        "bbox":
        {
          "maxx": -76.2527221,
          "maxy": 37.9068695,
          "maxz": -5.777103023,
          "minx": -76.2571767,
          "miny": 37.9028129,
          "minz": -35.71281856
        },
        "boundary": {
   "coordinates" : [
      [
         [ -76.257176700000002, 37.902812900000001 ],
         [ -76.257176700000002, 37.906869499999999 ],
         [ -76.2527221, 37.906869499999999 ],
         [ -76.2527221, 37.902812900000001 ],
         [ -76.257176700000002, 37.902812900000001 ]
      ]
   ],
   "type" : "Polygon"
}

      },
      "native":
      {
        "bbox":
        {
          "maxx": -76.2527221,
          "maxy": 37.9068695,
          "maxz": 29.966,
          "minx": -76.2571767,
          "miny": 37.9028129,
          "minz": 0.024
        },
        "boundary": {

pdal info on the output laz file gives:
{
  "filename": "\/data\/albers_20100503_LAS_S23_1836_10.laz",
  "pdal_version": "1.6.0 (git-version: 94ad83)",
  "stats":
  {
    "bbox":
    {
      "EPSG:4326":
      {
        "bbox":
        {
          "maxx": -76.25185156,
          "maxy": 37.90686951,
          "maxz": -5.77,
          "minx": -76.25806657,
          "miny": 37.9021939,
          "minz": -35.72
        },
        "boundary": {
   "coordinates" : [
      [
         [ -76.258066569999997, 37.903041880000004 ],
         [ -76.257028300000002, 37.90686951 ],
         [ -76.251851560000006, 37.90602148 ],
         [ -76.252890100000002, 37.9021939 ],
         [ -76.258066569999997, 37.903041880000004 ]
      ]
   ],
   "type" : "Polygon"
}

      },
      "native":
      {
        "bbox":
        {
          "maxx": 1707557.2,
          "maxy": 1830350.12,
          "maxz": -5.77,
          "minx": 1707096.44,
          "miny": 1829911.69,
          "minz": -35.72
        },
        "boundary": {


Am I missing a switch in the translate command?


Doug
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20171113/387353dc/attachment.html>