aboutsummaryrefslogtreecommitdiff
path: root/test
diff options
context:
space:
mode:
Diffstat (limited to 'test')
-rw-r--r--test/unit/CMakeLists.txt8
-rw-r--r--test/unit/Makefile.am12
-rw-r--r--test/unit/basic_test.cpp44
-rw-r--r--test/unit/pj_transform_test.cpp583
4 files changed, 593 insertions, 54 deletions
diff --git a/test/unit/CMakeLists.txt b/test/unit/CMakeLists.txt
index 8a25d7a3..68a6e60f 100644
--- a/test/unit/CMakeLists.txt
+++ b/test/unit/CMakeLists.txt
@@ -58,13 +58,13 @@ endif()
include_directories(${CMAKE_SOURCE_DIR}/include)
include_directories(${SQLITE3_INCLUDE_DIR})
-add_executable(proj_test_unit
+add_executable(proj_pj_transform_test
main.cpp
- basic_test.cpp)
-target_link_libraries(proj_test_unit
+ pj_transform_test.cpp)
+target_link_libraries(proj_pj_transform_test
gtest
${PROJ_LIBRARIES})
-add_test(NAME proj_test_unit COMMAND proj_test_unit)
+add_test(NAME proj_pj_transform_test COMMAND proj_pj_transform_test)
add_executable(proj_errno_string_test
main.cpp
diff --git a/test/unit/Makefile.am b/test/unit/Makefile.am
index 913a3f29..236ad444 100644
--- a/test/unit/Makefile.am
+++ b/test/unit/Makefile.am
@@ -9,17 +9,17 @@ AM_CXXFLAGS = @CXX_WFLAGS@ @NO_ZERO_AS_NULL_POINTER_CONSTANT_FLAG@
PROJ_LIB ?= ../../data
-noinst_PROGRAMS = basic_test
+noinst_PROGRAMS = pj_transform_test
noinst_PROGRAMS += pj_phi2_test
noinst_PROGRAMS += proj_errno_string_test
noinst_PROGRAMS += test_cpp_api
noinst_PROGRAMS += gie_self_tests
-basic_test_SOURCES = basic_test.cpp main.cpp
-basic_test_LDADD = ../../src/libproj.la ../../test/googletest/libgtest.la
+pj_transform_test_SOURCES = pj_transform_test.cpp main.cpp
+pj_transform_test_LDADD = ../../src/libproj.la ../../test/googletest/libgtest.la
-basic_test-check: basic_test
- ./basic_test
+pj_transform_test-check: pj_transform_test
+ PROJ_LIB=$(PROJ_LIB) ./pj_transform_test
pj_phi2_test_SOURCES = pj_phi2_test.cpp main.cpp
pj_phi2_test_LDADD = ../../src/libproj.la ../../test/googletest/libgtest.la
@@ -45,4 +45,4 @@ gie_self_tests_LDADD = ../../src/libproj.la ../../test/googletest/libgtest.la @S
gie_self_tests-check: gie_self_tests
PROJ_LIB=$(PROJ_LIB) ./gie_self_tests
-check-local: basic_test-check pj_phi2_test-check proj_errno_string_test-check test_cpp_api-check gie_self_tests-check
+check-local: pj_transform_test-check pj_phi2_test-check proj_errno_string_test-check test_cpp_api-check gie_self_tests-check
diff --git a/test/unit/basic_test.cpp b/test/unit/basic_test.cpp
deleted file mode 100644
index d565f839..00000000
--- a/test/unit/basic_test.cpp
+++ /dev/null
@@ -1,44 +0,0 @@
-/******************************************************************************
- *
- * Project: PROJ
- * Purpose: Test
- * Author: Even Rouault <even dot rouault at spatialys dot com>
- *
- ******************************************************************************
- * Copyright (c) 2018, Even Rouault <even dot rouault at spatialys dot com>
- *
- * Permission is hereby granted, free of charge, to any person obtaining a
- * copy of this software and associated documentation files (the "Software"),
- * to deal in the Software without restriction, including without limitation
- * the rights to use, copy, modify, merge, publish, distribute, sublicense,
- * and/or sell copies of the Software, and to permit persons to whom the
- * Software is furnished to do so, subject to the following conditions:
- *
- * The above copyright notice and this permission notice shall be included
- * in all copies or substantial portions of the Software.
- *
- * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
- * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
- * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
- * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
- * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
- * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
- * DEALINGS IN THE SOFTWARE.
- ****************************************************************************/
-
-#include "gtest_include.h"
-#include <proj.h>
-#include <memory>
-
-namespace {
-
-// ---------------------------------------------------------------------------
-
-TEST( dumy, dummy ) {
- std::unique_ptr<PJ, decltype(&::proj_destroy)> pj(
- ::proj_create(PJ_DEFAULT_CTX, "+proj=merc +ellps=clrk66 +lat_ts=33"),
- &::proj_destroy);
- ASSERT_NE(pj, nullptr);
-}
-
-} // namespace
diff --git a/test/unit/pj_transform_test.cpp b/test/unit/pj_transform_test.cpp
new file mode 100644
index 00000000..95329bdb
--- /dev/null
+++ b/test/unit/pj_transform_test.cpp
@@ -0,0 +1,583 @@
+/******************************************************************************
+ *
+ * Project: PROJ
+ * Purpose: Test pj_transform() legacy interface
+ * Author: Even Rouault <even dot rouault at spatialys dot com>
+ *
+ ******************************************************************************
+ * Copyright (c) 2018, Even Rouault <even dot rouault at spatialys dot com>
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a
+ * copy of this software and associated documentation files (the "Software"),
+ * to deal in the Software without restriction, including without limitation
+ * the rights to use, copy, modify, merge, publish, distribute, sublicense,
+ * and/or sell copies of the Software, and to permit persons to whom the
+ * Software is furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included
+ * in all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
+ * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
+ * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+ * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+ * DEALINGS IN THE SOFTWARE.
+ ****************************************************************************/
+
+#define ACCEPT_USE_OF_DEPRECATED_PROJ_API_H
+
+#include "gtest_include.h"
+#include <memory>
+
+// PROJ include order is sensitive
+// clang-format off
+#include <projects.h>
+#include <proj_api.h>
+// clang-format on
+
+namespace {
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_longlat) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_EQ(x, 2 * DEG_TO_RAD);
+ EXPECT_EQ(y, 49 * DEG_TO_RAD);
+
+ x = 182 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_EQ(x, 182 * DEG_TO_RAD);
+ EXPECT_EQ(y, 49 * DEG_TO_RAD);
+
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_proj) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=utm +zone=31 +datum=WGS84");
+ double x = 3 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 500000, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_proj_tometer) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=utm +zone=31 +datum=WGS84 +to_meter=1000");
+ double x = 3 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 500, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, proj_to_longlat) {
+ auto src = pj_init_plus("+proj=utm +zone=31 +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 500000;
+ double y = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 3 * DEG_TO_RAD, 1e-12);
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, proj_to_proj) {
+ auto src = pj_init_plus("+proj=utm +zone=31 +datum=WGS84");
+ auto dst = pj_init_plus("+proj=utm +zone=31 +datum=WGS84");
+ double x = 500000;
+ double y = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 500000, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_geocent) {
+ auto src = pj_init_plus("+proj=longlat +R=2");
+ auto dst = pj_init_plus("+proj=geocent +R=2");
+ double x = 0;
+ double y = 0;
+ double z = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ EXPECT_NEAR(z, 0, 1e-8);
+
+ x = 90 * DEG_TO_RAD;
+ y = 0;
+ z = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 0, 1e-8);
+ EXPECT_NEAR(y, 2, 1e-8);
+ EXPECT_NEAR(z, 0, 1e-8);
+
+ x = 0;
+ y = 90 * DEG_TO_RAD;
+ z = 0.1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 0, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ EXPECT_NEAR(z, 2 + 0.1, 1e-8);
+
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_geocent_to_meter) {
+ auto src = pj_init_plus("+proj=longlat +R=2");
+ auto dst = pj_init_plus("+proj=geocent +R=2 +to_meter=1000");
+ double x = 0;
+ double y = 0;
+ double z = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2e-3, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ EXPECT_NEAR(z, 0, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, geocent_to_longlat) {
+ auto src = pj_init_plus("+proj=geocent +R=2");
+ auto dst = pj_init_plus("+proj=longlat +R=2");
+ double x = 0;
+ double y = 2;
+ double z = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 0, 1e-12);
+
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, geocent_to_meter_to_longlat) {
+ auto src = pj_init_plus("+proj=geocent +to_meter=1000 +R=2");
+ auto dst = pj_init_plus("+proj=longlat +R=2");
+ double x = 0;
+ double y = 2e-3;
+ double z = 0;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 0, 1e-12);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, pm) {
+ auto src = pj_init_plus("+proj=longlat +pm=3 +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +pm=1 +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, (2 + 3 - 1) * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_EQ(y, 49 * DEG_TO_RAD) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_geoc_to_longlat) {
+ auto src = pj_init_plus("+proj=longlat +geoc +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 48.809360314691766 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_to_longlat_geoc) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +geoc +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 48.809360314691766 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, ellps_to_ellps_noop) {
+ auto src = pj_init_plus("+proj=longlat +ellps=clrk66");
+ auto dst = pj_init_plus("+proj=longlat +ellps=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, towgs84_3param_noop) {
+ auto src = pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=1,2,3");
+ auto dst = pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=1,2,3");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 10;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, towgs84_7param_noop) {
+ auto src =
+ pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=1,2,3,4,5,6,7");
+ auto dst =
+ pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=1,2,3,4,5,6,7");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 10;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_towgs84_3param_to_datum) {
+ auto src = pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=0,1,0");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 90 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ double z = 10;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10 + 1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_towgs84_3param_to_datum_no_z) {
+ auto src = pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=0,1,0");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 90 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_towgs84_7param_to_datum) {
+ auto src =
+ pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=0,1,0,0,0,0,0.5");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 90 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ double z = 10;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 14.189073500223458, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, datum_to_longlat_towgs84_3param) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=0,1,0");
+ double x = 90 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ double z = 10 + 1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, datum_to_longlat_towgs84_7param) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst =
+ pj_init_plus("+proj=longlat +ellps=WGS84 +towgs84=0,1,0,0,0,0,0.5");
+ double x = 90 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ double z = 14.189073500223458;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 90 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, ellps_grs80_towgs84_to_datum_wgs84) {
+ auto src = pj_init_plus("+proj=longlat +ellps=GRS80 +towgs84=0,0,0");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 10;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-15) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-15) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, longlat_nadgrids_to_datum) {
+ auto src = pj_init_plus("+proj=longlat +ellps=clrk66 +nadgrids=conus");
+ auto dst = pj_init_plus("+proj=longlat +datum=NAD83");
+ double x = -100 * DEG_TO_RAD;
+ double y = 40 * DEG_TO_RAD;
+ double z = 10;
+ int ret = pj_transform(src, dst, 1, 0, &x, &y, &z);
+ EXPECT_TRUE(ret == 0 || ret == PJD_ERR_FAILED_TO_LOAD_GRID);
+ if (ret == 0) {
+ EXPECT_NEAR(x, -100.00040583667015 * DEG_TO_RAD, 1e-12)
+ << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 40.000005895651363 * DEG_TO_RAD, 1e-12)
+ << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10.000043224543333, 1e-8);
+ }
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, nadgrids_noop) {
+ auto src = pj_init_plus("+proj=longlat +ellps=clrk66 +nadgrids=conus");
+ auto dst = pj_init_plus("+proj=longlat +ellps=clrk66 +nadgrids=conus");
+ double x = -100 * DEG_TO_RAD;
+ double y = 40 * DEG_TO_RAD;
+ double z = 10;
+ int ret = pj_transform(src, dst, 1, 0, &x, &y, &z);
+ EXPECT_TRUE(ret == 0);
+ if (ret == 0) {
+ EXPECT_NEAR(x, -100 * DEG_TO_RAD, 1e-15) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 40 * DEG_TO_RAD, 1e-15) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ }
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, datum_to_longlat_nadgrids) {
+ auto src = pj_init_plus("+proj=longlat +datum=NAD83");
+ auto dst = pj_init_plus("+proj=longlat +ellps=clrk66 +nadgrids=conus");
+ double x = -100.00040583667015 * DEG_TO_RAD;
+ double y = 40.000005895651363 * DEG_TO_RAD;
+ double z = 10.000043224543333;
+ int ret = pj_transform(src, dst, 1, 0, &x, &y, &z);
+ EXPECT_TRUE(ret == 0 || ret == PJD_ERR_FAILED_TO_LOAD_GRID);
+ if (ret == 0) {
+ EXPECT_NEAR(x, -100 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 40 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 10, 1e-8);
+ }
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, long_wrap) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84 +lon_wrap=180");
+ double x = -1 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 359 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 0 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, src_vto_meter) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84 +vto_meter=1000");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 1000, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, dest_vto_meter) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84 +vto_meter=1000");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 1000;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, src_axis_neu_with_z) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84 +axis=neu");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 49 * DEG_TO_RAD;
+ double y = 2 * DEG_TO_RAD;
+ double z = 1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, src_axis_neu_without_z) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84 +axis=neu");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 49 * DEG_TO_RAD;
+ double y = 2 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, src_axis_swd) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84 +axis=swd");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84");
+ double x = 49 * DEG_TO_RAD;
+ double y = 2 * DEG_TO_RAD;
+ double z = -1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, -2 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, -49 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, dst_axis_neu) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84 +axis=neu");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, 49 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, 2 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, 1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, dst_axis_swd) {
+ auto src = pj_init_plus("+proj=longlat +datum=WGS84");
+ auto dst = pj_init_plus("+proj=longlat +datum=WGS84 +axis=swd");
+ double x = 2 * DEG_TO_RAD;
+ double y = 49 * DEG_TO_RAD;
+ double z = 1;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, &z), 0);
+ EXPECT_NEAR(x, -49 * DEG_TO_RAD, 1e-12) << x / DEG_TO_RAD;
+ EXPECT_NEAR(y, -2 * DEG_TO_RAD, 1e-12) << y / DEG_TO_RAD;
+ EXPECT_NEAR(z, -1, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+// ---------------------------------------------------------------------------
+
+TEST(pj_transform_test, init_epsg) {
+ auto src = pj_init_plus("+init=epsg:4326");
+ ASSERT_TRUE(src != nullptr);
+ auto dst = pj_init_plus("+init=epsg:32631");
+ double x = 3 * DEG_TO_RAD;
+ double y = 0 * DEG_TO_RAD;
+ EXPECT_EQ(pj_transform(src, dst, 1, 0, &x, &y, nullptr), 0);
+ EXPECT_NEAR(x, 500000, 1e-8);
+ EXPECT_NEAR(y, 0, 1e-8);
+ pj_free(src);
+ pj_free(dst);
+}
+
+} // namespace
generated by cgit v1.2.3 (git 2.25.1) at 2025-12-29 02:16:09 +0000