<div dir="ltr">Hi Raife, <div><br></div><div>I think what is missing is the explicit vertical datum associated with that code.  I've come across this myself when converting a DTM to ellipsoidal heights.  Here's an example from my bash history that converts the heights with gdalwarp:</div><div><br></div><div>gdalwarp -s_srs "+proj=longlat +datum=WGS84 +no_defs +geoidgrids=egm96_15.gtx" -t_srs "EPSG:4326" dem-geoid.tif dem-hae.tif<br></div><div><br></div><div>I think the trick is just figuring out the right proj strings/finding the right EPSG compound code (there's an example at the bottom of this page <a href="https://gdal.org/programs/gdalwarp.html#gdalwarp">https://gdal.org/programs/gdalwarp.html#gdalwarp</a>).  I think under the hood it's using this fellow, so at the end of the day you just have to find the right grid:</div><div><br></div><div><a href="https://gdal.org/api/gdal_alg.html?highlight=applyvert#_CPPv426GDALApplyVerticalShiftGrid12GDALDatasetH12GDALDatasetHiddPPCKc">https://gdal.org/api/gdal_alg.html?highlight=applyvert#_CPPv426GDALApplyVerticalShiftGrid12GDALDatasetH12GDALDatasetHiddPPCKc</a><br></div><div><br></div><div>P</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Apr 18, 2021 at 9:17 AM Raife Crookes <<a href="mailto:raife@tp-regen.co.uk">raife@tp-regen.co.uk</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">





<div lang="EN-GB" style="overflow-wrap: break-word;">
<div class="gmail-m_8613101082346526930WordSection1">
<p class="MsoNormal">Good evening,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am hoping that someone could assist me with the completion of a reprojection of a geotiff lidar data set to WGS84.
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">The transformation should me making an elevation change of roughly 50m for known points however after running GDAL warp the elevations are staying the same when viewed through QGIS.
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">For some reason it seems to be defaulting to the WGS84 geoid height not the ellipsoid height I require for import onto our drone. (My knowledge is that EPSG:4326 should equate to the wgs lat/long and ellipsoid height.)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">GDALWARP instructions:<u></u><u></u></p>
<p class="MsoNormal">gdalwarp  -s_srs EPSG:27700 -t_srs EPSG:4326<u></u><u></u></p>
<p class="MsoNormal">gdalwarp  -s_srs -t_srs EPSG:4326<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Proj info return:<u></u><u></u></p>
<p class="MsoNormal">C:\>projinfo -s EPSG:27700 -t EPSG:4326 -o PROJ<u></u><u></u></p>
<p class="MsoNormal">Candidate operations found: 1<u></u><u></u></p>
<p class="MsoNormal">Note: using '--spatial-test intersects' would bring more results (7)<u></u><u></u></p>
<p class="MsoNormal">-------------------------------------<u></u><u></u></p>
<p class="MsoNormal">Operation No. 1:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">unknown id, Inverse of British National Grid + OSGB 1936 to WGS 84 (9), 1 m, UK - Britain and UKCS 49<span style="font-family:Arial,sans-serif">┬░</span>46'N to 61<span style="font-family:Arial,sans-serif">┬░</span>01'N, 7<span style="font-family:Arial,sans-serif">┬░</span>33'W
 to 3<span style="font-family:Arial,sans-serif">┬░</span>33'E<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">PROJ string:<u></u><u></u></p>
<p class="MsoNormal">+proj=pipeline +step +inv +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +step +proj=hgridshift +grids=OSTN15_NTv2_OSGBtoETRS.gsb +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap
 +order=2,1<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I can send copies of the original file if requested. <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Hoping that its something stupid like im using the wrong EPSG code but have a terrible feeling that a catch all to prevent users accidentally using the ellipsoid height is catching me out.
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">(Ground height OSGB – 30m  in WGS84 = 82m)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Really look forward to your help/<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Kind regards<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Raife<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>

_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>