From d87573cf62b339f22b24595fd05714f1b4e3bda3 Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Wed, 28 Nov 2018 20:13:42 +0100
Subject: exportToWKT WKT1_GDAL: export axis by default for GeocentricCRS

---
 src/crs.cpp | 9 +++++++++
 1 file changed, 9 insertions(+)

(limited to 'src/crs.cpp')

diff --git a/src/crs.cpp b/src/crs.cpp
index eec7a926..aa9209cd 100644
--- a/src/crs.cpp
+++ b/src/crs.cpp
@@ -861,7 +861,16 @@ void GeodeticCRS::_exportToWKT(io::WKTFormatter *formatter) const {
     if (!isWKT2) {
         unit._exportToWKT(formatter);
     }
+
+    const auto oldAxisOutputRule = formatter->outputAxis();
+    if (oldAxisOutputRule ==
+            io::WKTFormatter::OutputAxisRule::WKT1_GDAL_EPSG_STYLE &&
+        isGeocentric()) {
+        formatter->setOutputAxis(io::WKTFormatter::OutputAxisRule::YES);
+    }
     cs->_exportToWKT(formatter);
+    formatter->setOutputAxis(oldAxisOutputRule);
+
     ObjectUsage::baseExportToWKT(formatter);
 
     if (!isWKT2 && !formatter->useESRIDialect()) {
-- 
cgit v1.2.3