From f4827e0a3b5c5b07ccf82fc6cd3334079e75e02c Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Thu, 2 Apr 2020 22:36:37 +0200
Subject: Make sure that importing a Projected 3D CRS from WKT:2019 keeps the
 base geographic CRS as 3D (fixes #2122)

---
 src/iso19111/io.cpp | 6 ++++++
 1 file changed, 6 insertions(+)

(limited to 'src')

diff --git a/src/iso19111/io.cpp b/src/iso19111/io.cpp
index 6a1d7e32..8bfb4395 100644
--- a/src/iso19111/io.cpp
+++ b/src/iso19111/io.cpp
@@ -3839,6 +3839,12 @@ WKTParser::Private::buildProjectedCRS(const WKTNodeNNPtr &node) {
         ThrowNotExpectedCSType("Cartesian");
     }
 
+    if (cartesianCS->axisList().size() == 3 &&
+        baseGeodCRS->coordinateSystem()->axisList().size() == 2) {
+        baseGeodCRS = NN_NO_CHECK(util::nn_dynamic_pointer_cast<GeodeticCRS>(
+            baseGeodCRS->promoteTo3D(std::string(), dbContext_)));
+    }
+
     addExtensionProj4ToProp(nodeP, props);
 
     return ProjectedCRS::create(props, baseGeodCRS, conversion,
-- 
cgit v1.2.3