<div dir="ltr">Hi All,<div><br></div><div>I am trying to implement a warp from a latitude, longitude, altitude (LLA) image into a custom coordinate system using the GDAL Warp API.  The transform between LLA and my coordinate requires scaling and shifting the elevation.  My coordinate system is  a local tangent plane on the Earth's surface where x is north, y is east, and z is pointed into the Earth.</div><div><br></div><div>When I follow the instructions in the <a href="http://www.gdal.org/warptut.html">GDAL Warp API tutorial</a>,  my custom transform function supplied to "<span style="font-size:13px;color:rgb(0,0,0);font-family:monospace,fixed;line-height:13px;white-space:pre-wrap;background-color:rgb(251,252,253)">psWarpOptions-></span><a class="" href="http://www.gdal.org/structGDALWarpOptions.html#a35647845a2629da876f379f82d58a2cd" style="font-size:13px;color:rgb(70,101,162);text-decoration:none;font-family:monospace,fixed;line-height:13px;white-space:pre-wrap;background-color:rgb(251,252,253)">pTransformerArg</a>" only receives input z=0 values.  From the source code, it looks like the internal GDALWarpKernel only passes z=0 values to to the pTransformerArg.  The Z values from the source elevation map are never used and therefore can never be re-projected.<br></div><div><br></div><div>How can GDAL properly transform the source elevation data if the source elevation is never passed to the transform function? </div><div><br></div><div> How are transforms between ellipsoids or elevation units performed if not through the <a class="" href="http://www.gdal.org/classGDALWarpOperation.html" style="font-size:13px;color:rgb(70,101,162);text-decoration:none;font-family:monospace,fixed;line-height:13px;white-space:pre-wrap;background-color:rgb(251,252,253)">GDALWarpOperation</a> API?<br></div><div><br></div><div><br></div><div>Thank you for any insights.</div><div><br></div><div>-Adam Stambler</div><div><br></div></div>