From 590449950057566d39e2997404422b00844343b6 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Wed, 3 Jun 2020 14:09:01 +0200 Subject: normalizeForVisualization(): make it switch axis for EPSG:5482 (RSRGD2000 / RSPS2000) Fixes #2254 --- test/unit/test_crs.cpp | 32 ++++++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) (limited to 'test') diff --git a/test/unit/test_crs.cpp b/test/unit/test_crs.cpp index 8fdd8a01..c3530a53 100644 --- a/test/unit/test_crs.cpp +++ b/test/unit/test_crs.cpp @@ -1057,6 +1057,38 @@ TEST(crs, EPSG_5042_projected_south_pole_east_north) { // --------------------------------------------------------------------------- +TEST(crs, EPSG_5482_projected_south_pole_south_west) { + auto dbContext = DatabaseContext::create(); + auto factory = AuthorityFactory::create(dbContext, "EPSG"); + auto crs = factory->createCoordinateReferenceSystem("5482"); + auto proj_crs = nn_dynamic_pointer_cast(crs); + ASSERT_TRUE(proj_crs != nullptr); + auto op = CoordinateOperationFactory::create()->createOperation( + factory->createCoordinateReferenceSystem("4764"), + NN_NO_CHECK(proj_crs)); + ASSERT_TRUE(op != nullptr); + auto proj_string = "+proj=pipeline " + "+step +proj=axisswap +order=2,1 " + "+step +proj=unitconvert +xy_in=deg +xy_out=rad " + "+step +proj=stere +lat_0=-90 +lon_0=180 +k=0.994 " + "+x_0=5000000 +y_0=1000000 +ellps=GRS80 " + "+step +proj=axisswap +order=2,1"; + EXPECT_EQ(op->exportToPROJString(PROJStringFormatter::create().get()), + proj_string); + + auto opNormalized = op->normalizeForVisualization(); + auto proj_string_normalized = + "+proj=pipeline " + "+step +proj=unitconvert +xy_in=deg +xy_out=rad " + "+step +proj=stere +lat_0=-90 +lon_0=180 +k=0.994 " + "+x_0=5000000 +y_0=1000000 +ellps=GRS80"; + EXPECT_EQ( + opNormalized->exportToPROJString(PROJStringFormatter::create().get()), + proj_string_normalized); +} + +// --------------------------------------------------------------------------- + TEST(crs, geodetic_crs_both_datum_datum_ensemble_null) { EXPECT_THROW(GeodeticCRS::create( PropertyMap(), nullptr, nullptr, -- cgit v1.2.3