From 3d75435aa978588cf929a57286a2ae13d8ea128c Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Tue, 20 Apr 2021 11:20:16 +0200 Subject: Merge pull request #2678 from rouault/revert_ebe3425bf Revert "proj_create_crs_to_crs_from_pj(): do not use PROJ_SPATIAL_CRITERION_PARTIAL_INTERSECTION if area is specified" --- test/unit/gie_self_tests.cpp | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) (limited to 'test') diff --git a/test/unit/gie_self_tests.cpp b/test/unit/gie_self_tests.cpp index 92f1c758..59b78497 100644 --- a/test/unit/gie_self_tests.cpp +++ b/test/unit/gie_self_tests.cpp @@ -804,6 +804,30 @@ TEST(gie, proj_create_crs_to_crs_outside_area_of_use) { // --------------------------------------------------------------------------- +TEST(gie, proj_create_crs_to_crs_with_area_large) { + + // Test bugfix for https://github.com/OSGeo/gdal/issues/3695 + auto area = proj_area_create(); + proj_area_set_bbox(area, -14.1324, 49.5614, 3.76488, 62.1463); + auto P = proj_create_crs_to_crs(PJ_DEFAULT_CTX, "EPSG:4277", "EPSG:4326", + area); + proj_area_destroy(area); + ASSERT_TRUE(P != nullptr); + PJ_COORD c; + + c.xyzt.x = 50; // Lat in deg + c.xyzt.y = -2; // Long in deg + c.xyzt.z = 0; + c.xyzt.t = HUGE_VAL; + c = proj_trans(P, PJ_FWD, c); + EXPECT_NEAR(c.xy.x, 50.00065628, 1e-8); + EXPECT_NEAR(c.xy.y, -2.00133989, 1e-8); + + proj_destroy(P); +} + +// --------------------------------------------------------------------------- + TEST(gie, proj_trans_generic) { // GDA2020 to WGS84 (G1762) auto P = proj_create( -- cgit v1.2.3