From fa4c5c645fc5061a306edaef41a9a091189ce528 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Mon, 16 Aug 2021 13:49:03 +0200 Subject: Merge pull request #2807 from rouault/fix_2804 createOperations(): fix missing deg<-->rad conversion when transforming with a CRS that has a fallback-to-PROJ4-string behaviour and is a BoundCRS of a GeographicCRS (fixes #2804) --- test/unit/test_operationfactory.cpp | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) (limited to 'test') diff --git a/test/unit/test_operationfactory.cpp b/test/unit/test_operationfactory.cpp index 47c6636c..56d9f743 100644 --- a/test/unit/test_operationfactory.cpp +++ b/test/unit/test_operationfactory.cpp @@ -6108,6 +6108,30 @@ TEST(operation, createOperation_fallback_to_proj4_strings) { // --------------------------------------------------------------------------- +TEST(operation, createOperation_fallback_to_proj4_strings_bound_of_geog) { + auto objSrc = PROJStringParser().createFromPROJString( + "+proj=longlat +geoc +ellps=GRS80 +towgs84=0,0,0 +type=crs"); + auto src = nn_dynamic_pointer_cast(objSrc); + ASSERT_TRUE(src != nullptr); + + auto objDest = PROJStringParser().createFromPROJString( + "+proj=longlat +geoc +ellps=clrk66 +towgs84=0,0,0 +type=crs"); + auto dest = nn_dynamic_pointer_cast(objDest); + ASSERT_TRUE(dest != nullptr); + + auto op = CoordinateOperationFactory::create()->createOperation( + NN_CHECK_ASSERT(src), NN_CHECK_ASSERT(dest)); + ASSERT_TRUE(op != nullptr); + EXPECT_EQ(op->exportToPROJString(PROJStringFormatter::create().get()), + "+proj=pipeline " + "+step +proj=unitconvert +xy_in=deg +xy_out=rad " + "+step +inv +proj=longlat +geoc +ellps=GRS80 +towgs84=0,0,0 " + "+step +proj=longlat +geoc +ellps=clrk66 +towgs84=0,0,0 " + "+step +proj=unitconvert +xy_in=rad +xy_out=deg"); +} + +// --------------------------------------------------------------------------- + TEST( operation, createOperation_fallback_to_proj4_strings_regular_with_datum_to_projliteral) { -- cgit v1.2.3