createOperations(): fix vertical to geographic when synthetizing an operation that involves a vertical axis reversal

1 files changed, 12 insertions, 0 deletions

diff --git a/test/unit/test_operation.cpp b/test/unit/test_operation.cpp
index 9f6a808f..9cde8822 100644
--- a/test/unit/test_operation.cpp
+++ b/test/unit/test_operation.cpp
@@ -4675,6 +4675,18 @@ TEST(operation, vertCRS_to_geogCRS_context) {
             "+proj=pipeline +step +inv +proj=vgridshift +grids=egm08_25.gtx "
             "+multiplier=1");
     }
+    {
+        auto ctxt =
+            CoordinateOperationContext::create(authFactory, nullptr, 0.0);
+        auto list = CoordinateOperationFactory::create()->createOperations(
+            // NGVD29 depth (ftUS)
+            authFactory->createCoordinateReferenceSystem("6359"),
+            authFactory->createCoordinateReferenceSystem("4326"), ctxt);
+        ASSERT_EQ(list.size(), 1U);
+        EXPECT_EQ(
+            list[0]->exportToPROJString(PROJStringFormatter::create().get()),
+            "+proj=affine +s33=-0.304800609601219");
+    }
 }
 
 // ---------------------------------------------------------------------------