Re: [gdal-dev] Problem with file source projection doing a vertical datum reprojection

2023-04-12 Thread Even Rouault
Le 12/04/2023 à 11:45, Michael Smith a écrit : I guess I thought GTIFF_REPORT_COMPD_CS would only change the display of the output, not change the interpretation of the SRS? It is documented in https://gdal.org/drivers/raster/gtiff.html as: """GTIFF_REPORT_COMPD_CS : Can be set to TRUE to

Re: [gdal-dev] Problem with file source projection doing a vertical datum reprojection

2023-04-12 Thread Michael Smith
I guess I thought GTIFF_REPORT_COMPD_CS would only change the display of the output, not change the interpretation of the SRS? I can confirm that setting GTIFF_REPORT_COMPD_CS=FALSE does indeed allow the vertical datum reprojection to happen but I guess I have a poor understanding of just what

Re: [gdal-dev] Problem with file source projection doing a vertical datum reprojection

2023-04-12 Thread Even Rouault
Michael, it is well possible this might have worked in the past due to slightly different behavior of GDAL and/or PROJ, but the current behavior seems consistent to me. You're trying to transform between 2 compound CRS, the source one with a unknown vertical datum and the target one with

[gdal-dev] Problem with file source projection doing a vertical datum reprojection

2023-04-12 Thread Michael Smith
I have a file with a projection shown, using REPORT_COMPD_CS=TRUE of: COMPOUNDCRS["WGS 84 / UTM zone 17N + unknown",     PROJCRS["WGS 84 / UTM zone 17N",     BASEGEOGCRS["WGS 84",     DATUM["World Geodetic System 1984",     ELLIPSOID["WGS 84",6378137,298.257223563,