<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=EN-US link="#0563C1" vlink="#954F72" style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal>I have a file with a projection shown, using REPORT_COMPD_CS=TRUE of:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>COMPOUNDCRS["WGS 84 / UTM zone 17N + unknown",<o:p></o:p></p><p class=MsoNormal> PROJCRS["WGS 84 / UTM zone 17N",<o:p></o:p></p><p class=MsoNormal> BASEGEOGCRS["WGS 84",<o:p></o:p></p><p class=MsoNormal> DATUM["World Geodetic System 1984",<o:p></o:p></p><p class=MsoNormal> ELLIPSOID["WGS 84",6378137,298.257223563,<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1]]],<o:p></o:p></p><p class=MsoNormal> PRIMEM["Greenwich",0,<o:p></o:p></p><p class=MsoNormal> ANGLEUNIT["degree",0.0174532925199433]],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",4326]],<o:p></o:p></p><p class=MsoNormal> CONVERSION["UTM zone 17N",<o:p></o:p></p><p class=MsoNormal> METHOD["Transverse Mercator",<o:p></o:p></p><p class=MsoNormal> ID["EPSG",9807]],<o:p></o:p></p><p class=MsoNormal> PARAMETER["Latitude of natural origin",0,<o:p></o:p></p><p class=MsoNormal> ANGLEUNIT["degree",0.0174532925199433],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",8801]],<o:p></o:p></p><p class=MsoNormal> PARAMETER["Longitude of natural origin",-81,<o:p></o:p></p><p class=MsoNormal> ANGLEUNIT["degree",0.0174532925199433],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",8802]],<o:p></o:p></p><p class=MsoNormal> PARAMETER["Scale factor at natural origin",0.9996,<o:p></o:p></p><p class=MsoNormal> SCALEUNIT["unity",1],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",8805]],<o:p></o:p></p><p class=MsoNormal> PARAMETER["False easting",500000,<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",8806]],<o:p></o:p></p><p class=MsoNormal> PARAMETER["False northing",0,<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",8807]]],<o:p></o:p></p><p class=MsoNormal> CS[Cartesian,2],<o:p></o:p></p><p class=MsoNormal> AXIS["easting",east,<o:p></o:p></p><p class=MsoNormal> ORDER[1],<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1]],<o:p></o:p></p><p class=MsoNormal> AXIS["northing",north,<o:p></o:p></p><p class=MsoNormal> ORDER[2],<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1]],<o:p></o:p></p><p class=MsoNormal> ID["EPSG",32617]],<o:p></o:p></p><p class=MsoNormal> VERTCRS["unknown",<o:p></o:p></p><p class=MsoNormal> VDATUM["unknown"],<o:p></o:p></p><p class=MsoNormal> CS[vertical,1],<o:p></o:p></p><p class=MsoNormal> AXIS["up",up,<o:p></o:p></p><p class=MsoNormal> LENGTHUNIT["metre",1,<o:p></o:p></p><p class=MsoNormal> ID["EPSG",9001]]]]]<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>When trying to gdalwarp to EGM2008 using gdalwarp -t_srs epsg:32617+3855, vertical datum reprojection does not occur unless I overwride the source srs with -s_srs epsg:32617. Using GDAL 3.6.2 and proj 9.1.1. With proj_debug, I can see that the PROJ_TRACE: vgridshift: is not occurring unless the srs is overridden. I’m fairly sure that this has worked in the past but I can’t say exactly which version. This is an older geotiff that was written about 7 years ago. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>-- <o:p></o:p></p></div><div><p class=MsoNormal>Michael Smith<o:p></o:p></p><p class=MsoNormal>US Army Corps / Remote Sensing GIS Center<o:p></o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>