diff options
Diffstat (limited to 'test/unit')
| -rw-r--r-- | test/unit/test_c_api.cpp | 211 | ||||
| -rw-r--r-- | test/unit/test_factory.cpp | 658 |
2 files changed, 846 insertions, 23 deletions
diff --git a/test/unit/test_c_api.cpp b/test/unit/test_c_api.cpp index 20fb0583..b624e063 100644 --- a/test/unit/test_c_api.cpp +++ b/test/unit/test_c_api.cpp @@ -1776,29 +1776,58 @@ TEST_F(CApi, proj_context_set_database_path_null) { // --------------------------------------------------------------------------- -TEST_F(CApi, proj_context_set_database_path_main_memory_one_aux) { +TEST_F(CApi, proj_context_set_database_path_aux) { + + const std::string auxDbName( + "file:proj_test_aux.db?mode=memory&cache=shared"); + + sqlite3 *dbAux = nullptr; + sqlite3_open_v2( + auxDbName.c_str(), &dbAux, + SQLITE_OPEN_READWRITE | SQLITE_OPEN_CREATE | SQLITE_OPEN_URI, nullptr); + ASSERT_TRUE(dbAux != nullptr); + ASSERT_TRUE(sqlite3_exec(dbAux, "BEGIN", nullptr, nullptr, nullptr) == + SQLITE_OK); + { + auto ctxt = DatabaseContext::create(); + const auto dbStructure = ctxt->getDatabaseStructure(); + for (const auto &sql : dbStructure) { + ASSERT_TRUE(sqlite3_exec(dbAux, sql.c_str(), nullptr, nullptr, + nullptr) == SQLITE_OK); + } + } - auto c_path = proj_context_get_database_path(m_ctxt); - ASSERT_TRUE(c_path != nullptr); - std::string path(c_path); - const char *aux_db_list[] = {path.c_str(), nullptr}; - - // This is super exotic and a miracle that it works. :memory: as the - // main DB is empty. The real stuff is in the aux_db_list. No view - // is created in the ':memory:' internal DB, but as there's only one - // aux DB its tables and views can be directly queried... - // If that breaks at some point, that wouldn't be a big issue. - // Keeping that one as I had a hard time figuring out why it worked ! - // The real thing is tested by the C++ - // factory::attachExtraDatabases_auxiliary - EXPECT_TRUE(proj_context_set_database_path(m_ctxt, ":memory:", aux_db_list, - nullptr)); + ASSERT_TRUE(sqlite3_exec( + dbAux, + "INSERT INTO geodetic_crs VALUES('OTHER','OTHER_4326','WGS " + "84',NULL,'geographic 2D','EPSG','6422','EPSG','6326'," + "NULL,0);", + nullptr, nullptr, nullptr) == SQLITE_OK); + ASSERT_TRUE(sqlite3_exec(dbAux, "COMMIT", nullptr, nullptr, nullptr) == + SQLITE_OK); - auto source_crs = proj_create_from_database(m_ctxt, "EPSG", "4326", - PJ_CATEGORY_CRS, false, - nullptr); // WGS84 - ASSERT_NE(source_crs, nullptr); - ObjectKeeper keeper_source_crs(source_crs); + const char *const aux_db_list[] = {auxDbName.c_str(), nullptr}; + + EXPECT_TRUE( + proj_context_set_database_path(m_ctxt, nullptr, aux_db_list, nullptr)); + + sqlite3_close(dbAux); + + { + auto crs = proj_create_from_database(m_ctxt, "EPSG", "4326", + PJ_CATEGORY_CRS, false, + nullptr); // WGS84 + ASSERT_NE(crs, nullptr); + ObjectKeeper keeper_source_crs(crs); + } + + { + auto crs = proj_create_from_database(m_ctxt, "OTHER", "OTHER_4326", + PJ_CATEGORY_CRS, false, + nullptr); // WGS84 + ASSERT_NE(crs, nullptr); + ObjectKeeper keeper_source_crs(crs); + } } // --------------------------------------------------------------------------- @@ -2728,6 +2757,15 @@ TEST_F(CApi, proj_context_get_database_metadata) { // --------------------------------------------------------------------------- +TEST_F(CApi, proj_context_get_database_structure) { + auto list = proj_context_get_database_structure(m_ctxt, nullptr); + ASSERT_NE(list, nullptr); + ASSERT_NE(list[0], nullptr); + proj_string_list_destroy(list); +} + +// --------------------------------------------------------------------------- + TEST_F(CApi, proj_clone) { auto obj = proj_create(m_ctxt, "+proj=longlat"); ObjectKeeper keeper(obj); @@ -5265,4 +5303,135 @@ TEST_F(CApi, proj_crs_is_derived) { } } +// --------------------------------------------------------------------------- + +TEST_F(CApi, proj_get_insert_statements) { + { + auto session = proj_insert_object_session_create(nullptr); + EXPECT_NE(session, nullptr); + + EXPECT_EQ(proj_insert_object_session_create(nullptr), nullptr); + + proj_insert_object_session_destroy(nullptr, session); + } + + { proj_insert_object_session_destroy(nullptr, nullptr); } + { + auto wkt = "GEOGCRS[\"myGDA2020\",\n" + " DATUM[\"GDA2020\",\n" + " ELLIPSOID[\"GRS_1980\",6378137,298.257222101,\n" + " LENGTHUNIT[\"metre\",1]]],\n" + " PRIMEM[\"Greenwich\",0,\n" + " ANGLEUNIT[\"Degree\",0.0174532925199433]],\n" + " CS[ellipsoidal,2],\n" + " AXIS[\"geodetic latitude (Lat)\",north,\n" + " ORDER[1],\n" + " ANGLEUNIT[\"degree\",0.0174532925199433]],\n" + " AXIS[\"geodetic longitude (Lon)\",east,\n" + " ORDER[2],\n" + " ANGLEUNIT[\"degree\",0.0174532925199433]]]"; + auto crs = proj_create_from_wkt(m_ctxt, wkt, nullptr, nullptr, nullptr); + ObjectKeeper keeper_from_wkt(crs); + EXPECT_NE(crs, nullptr); + + { + char *code = + proj_suggests_code_for(m_ctxt, crs, "HOBU", false, nullptr); + ASSERT_NE(code, nullptr); + EXPECT_EQ(std::string(code), "MYGDA2020"); + proj_string_destroy(code); + } + + { + char *code = + proj_suggests_code_for(m_ctxt, crs, "HOBU", true, nullptr); + ASSERT_NE(code, nullptr); + EXPECT_EQ(std::string(code), "1"); + proj_string_destroy(code); + } + + const auto sizeOfStringList = [](const char *const *list) { + if (list == nullptr) + return -1; + int size = 0; + for (auto iter = list; *iter; ++iter) { + size += 1; + } + return size; + }; + + // No session specified: we use a temporary session + for (int i = 0; i < 2; i++) { + auto list = proj_get_insert_statements( + m_ctxt, nullptr, crs, "HOBU", "XXXX", false, nullptr, nullptr); + ASSERT_NE(list, nullptr); + ASSERT_NE(list[0], nullptr); + EXPECT_EQ(std::string(list[0]), + "INSERT INTO geodetic_datum VALUES('HOBU'," + "'GEODETIC_DATUM_XXXX','GDA2020','','EPSG','7019'," + "'EPSG','8901',NULL,NULL,NULL,0);"); + EXPECT_EQ(sizeOfStringList(list), 4); + proj_string_list_destroy(list); + } + + // Pass an empty list of allowed authorities + // We cannot reuse the EPSG ellipsoid and prime meridian + { + const char *const allowed_authorities[] = {nullptr}; + auto list = + proj_get_insert_statements(m_ctxt, nullptr, crs, "HOBU", "XXXX", + false, allowed_authorities, nullptr); + EXPECT_EQ(sizeOfStringList(list), 6); + proj_string_list_destroy(list); + } + + // Allow only PROJ + // We cannot reuse the EPSG ellipsoid and prime meridian + { + const char *const allowed_authorities[] = {"PROJ", nullptr}; + auto list = + proj_get_insert_statements(m_ctxt, nullptr, crs, "HOBU", "XXXX", + false, allowed_authorities, nullptr); + EXPECT_EQ(sizeOfStringList(list), 6); + proj_string_list_destroy(list); + } + + // Allow EPSG + { + const char *const allowed_authorities[] = {"EPSG", nullptr}; + auto list = + proj_get_insert_statements(m_ctxt, nullptr, crs, "HOBU", "XXXX", + false, allowed_authorities, nullptr); + EXPECT_EQ(sizeOfStringList(list), 4); + proj_string_list_destroy(list); + } + + auto session = proj_insert_object_session_create(m_ctxt); + EXPECT_NE(session, nullptr); + + { + auto list = proj_get_insert_statements( + m_ctxt, session, crs, "HOBU", "XXXX", false, nullptr, nullptr); + ASSERT_NE(list, nullptr); + ASSERT_NE(list[0], nullptr); + EXPECT_EQ(std::string(list[0]), + "INSERT INTO geodetic_datum VALUES('HOBU'," + "'GEODETIC_DATUM_XXXX','GDA2020','','EPSG','7019'," + "'EPSG','8901',NULL,NULL,NULL,0);"); + proj_string_list_destroy(list); + } + + // Object already inserted: return empty list + { + auto list = proj_get_insert_statements( + m_ctxt, session, crs, "HOBU", "XXXX", false, nullptr, nullptr); + ASSERT_NE(list, nullptr); + ASSERT_EQ(list[0], nullptr); + proj_string_list_destroy(list); + } + + proj_insert_object_session_destroy(m_ctxt, session); + } +} + } // namespace diff --git a/test/unit/test_factory.cpp b/test/unit/test_factory.cpp index 454e4f2e..e0616baa 100644 --- a/test/unit/test_factory.cpp +++ b/test/unit/test_factory.cpp @@ -2828,10 +2828,12 @@ TEST(factory, attachExtraDatabases_auxiliary) { ASSERT_TRUE(dbAux != nullptr); ASSERT_TRUE(sqlite3_exec(dbAux, "BEGIN", nullptr, nullptr, nullptr) == SQLITE_OK); + + std::vector<std::string> tableStructureBefore; { auto ctxt = DatabaseContext::create(); - const auto dbStructure = ctxt->getDatabaseStructure(); - for (const auto &sql : dbStructure) { + tableStructureBefore = ctxt->getDatabaseStructure(); + for (const auto &sql : tableStructureBefore) { if (sql.find("CREATE TRIGGER") == std::string::npos) { ASSERT_TRUE(sqlite3_exec(dbAux, sql.c_str(), nullptr, nullptr, nullptr) == SQLITE_OK); @@ -2864,6 +2866,9 @@ TEST(factory, attachExtraDatabases_auxiliary) { auto gcrs = nn_dynamic_pointer_cast<GeographicCRS>(crs); EXPECT_TRUE(gcrs != nullptr); } + + const auto dbStructure = ctxt->getDatabaseStructure(); + EXPECT_EQ(dbStructure, tableStructureBefore); } { @@ -3434,4 +3439,653 @@ TEST(factory, getUnitList) { } } +// --------------------------------------------------------------------------- + +TEST(factory, objectInsertion) { + + // Cannot nest startInsertStatementsSession + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + EXPECT_THROW(ctxt->startInsertStatementsSession(), FactoryException); + } + + { + auto ctxt = DatabaseContext::create(); + // Tolerated withtout explicit stop + ctxt->startInsertStatementsSession(); + } + + { + auto ctxt = DatabaseContext::create(); + // Tolerated + ctxt->stopInsertStatementsSession(); + } + + // getInsertStatementsFor() must be preceded with + // startInsertStatementsSession() + { + auto ctxt = DatabaseContext::create(); + EXPECT_THROW(ctxt->getInsertStatementsFor(GeographicCRS::EPSG_4326, + "EPSG", "4326", true), + FactoryException); + } + + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + // Nothing to do + EXPECT_TRUE(ctxt->getInsertStatementsFor(GeographicCRS::EPSG_4326, + "EPSG", "4326", true) + .empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + GeographicCRS::EPSG_4326->datum(), + GeographicCRS::EPSG_4326->datumEnsemble(), + GeographicCRS::EPSG_4326->coordinateSystem()); + + EXPECT_EQ(ctxt->suggestsCodeFor(crs, "HOBU", true), "1"); + EXPECT_EQ(ctxt->suggestsCodeFor(crs, "HOBU", false), "MY_EPSG_4326"); + + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "1234", true); + + EXPECT_EQ(ctxt->suggestsCodeFor(crs, "HOBU", true), "1235"); + + ASSERT_EQ(sql.size(), 2U); + EXPECT_EQ(sql[0], "INSERT INTO geodetic_crs VALUES('HOBU','1234','my " + "EPSG:4326','','geographic " + "2D','EPSG','6422','EPSG','6326',NULL,0);"); + EXPECT_EQ( + sql[1], + "INSERT INTO usage " + "VALUES('HOBU','USAGE_GEODETIC_CRS_1234','geodetic_crs','HOBU','" + "1234','PROJ','EXTENT_UNKNOWN','PROJ','SCOPE_UNKNOWN');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "1234", true).empty()); + ctxt->stopInsertStatementsSession(); + AuthorityFactory::create(ctxt, std::string("EPSG")) + ->createGeographicCRS("4326"); + EXPECT_THROW(AuthorityFactory::create(ctxt, std::string("HOBU")) + ->createGeographicCRS("1234"), + NoSuchAuthorityCodeException); + } + + // Geographic 3D CRS, with known usage + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto usages = AuthorityFactory::create(ctxt, std::string("EPSG")) + ->createGeographicCRS("4979") + ->domains(); + auto array(ArrayOfBaseObject::create()); + for (const auto &usage : usages) { + array->add(usage); + } + auto props = + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4979"); + props.set(ObjectUsage::OBJECT_DOMAIN_KEY, + nn_static_pointer_cast<BaseObject>(array)); + const auto crs = + GeographicCRS::create(props, GeographicCRS::EPSG_4979->datum(), + GeographicCRS::EPSG_4979->datumEnsemble(), + GeographicCRS::EPSG_4979->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "4979", false); + ASSERT_EQ(sql.size(), 2U); + EXPECT_EQ(sql[0], "INSERT INTO geodetic_crs VALUES('HOBU','4979','my " + "EPSG:4979','','geographic " + "3D','EPSG','6423','EPSG','6326',NULL,0);"); + EXPECT_EQ( + sql[1], + "INSERT INTO usage " + "VALUES('HOBU','USAGE_GEODETIC_CRS_4979','geodetic_crs','HOBU','" + "4979','EPSG','1262','EPSG','1176');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "4979", false).empty()); + ctxt->stopInsertStatementsSession(); + } + + // BoundCRS of Geocentric CRS, with new usage + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + auto props = + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4978"); + auto array(ArrayOfBaseObject::create()); + const auto extent = Extent::createFromBBOX(1, 2, 3, 4); + optional<std::string> scope; + scope = "my scope"; + array->add(ObjectDomain::create(scope, extent)); + props.set(ObjectUsage::OBJECT_DOMAIN_KEY, + nn_static_pointer_cast<BaseObject>(array)); + const auto crs = GeodeticCRS::create( + props, NN_NO_CHECK(GeodeticCRS::EPSG_4978->datum()), + NN_NO_CHECK(nn_dynamic_pointer_cast<CartesianCS>( + GeodeticCRS::EPSG_4978->coordinateSystem()))); + const auto boundCRS = BoundCRS::createFromTOWGS84( + crs, std::vector<double>{1, 2, 3, 4, 5, 6, 7}); + const auto sql = + ctxt->getInsertStatementsFor(boundCRS, "HOBU", "4978", false); + ASSERT_EQ(sql.size(), 4U); + EXPECT_EQ( + sql[0], + "INSERT INTO geodetic_crs VALUES('HOBU','4978','my " + "EPSG:4978','','geocentric','EPSG','6500','EPSG','6326',NULL,0);"); + EXPECT_EQ(sql[1], + "INSERT INTO scope VALUES('HOBU','SCOPE_geodetic_crs_4978'," + "'my scope',0);"); + EXPECT_EQ(sql[2], + "INSERT INTO extent VALUES('HOBU','EXTENT_geodetic_crs_4978'," + "'unknown','unknown',2,4,1,3,0);"); + EXPECT_EQ( + sql[3], + "INSERT INTO usage VALUES('HOBU','USAGE_GEODETIC_CRS_4978'," + "'geodetic_crs','HOBU','4978','HOBU'," + "'EXTENT_geodetic_crs_4978','HOBU','SCOPE_geodetic_crs_4978');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(boundCRS, "HOBU", "4978", false) + .empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown datum, numeric code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + Ellipsoid::WGS84, optional<std::string>(), + PrimeMeridian::GREENWICH); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true); + ASSERT_EQ(sql.size(), 4U); + EXPECT_EQ(sql[0], + "INSERT INTO geodetic_datum VALUES('HOBU','1','my " + "datum','','EPSG','7030','EPSG','8901',NULL,NULL,NULL,0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown datum, alpha code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + Ellipsoid::WGS84, optional<std::string>(), + PrimeMeridian::GREENWICH); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "MY_EPSG_4326", false); + + EXPECT_EQ(ctxt->suggestsCodeFor(crs, "HOBU", false), "MY_EPSG_4326_2"); + + ASSERT_EQ(sql.size(), 4U); + EXPECT_EQ(sql[0], + "INSERT INTO geodetic_datum " + "VALUES('HOBU','GEODETIC_DATUM_MY_EPSG_4326','my " + "datum','','EPSG','7030','EPSG','8901',NULL,NULL,NULL,0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "MY_EPSG_4326", false) + .empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown ellipsoid, numeric code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto ellipsoid = Ellipsoid::createFlattenedSphere( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my ellipsoid"), + Length(6378137), Scale(295)); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + ellipsoid, optional<std::string>(), PrimeMeridian::GREENWICH); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true); + ASSERT_EQ(sql.size(), 5U); + EXPECT_EQ( + sql[0], + "INSERT INTO ellipsoid VALUES('HOBU','1','my " + "ellipsoid','','PROJ','EARTH',6378137,'EPSG','9001',295,NULL,0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown ellipsoid, alpha code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto ellipsoid = Ellipsoid::createTwoAxis( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my ellipsoid"), + Length(6378137), Length(6378136)); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + ellipsoid, optional<std::string>(), PrimeMeridian::GREENWICH); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false); + ASSERT_EQ(sql.size(), 5U); + EXPECT_EQ(sql[0], "INSERT INTO ellipsoid " + "VALUES('HOBU','ELLPS_GEODETIC_DATUM_XXXX','my " + "ellipsoid','','PROJ','EARTH',6378137,'EPSG','9001'," + "NULL,6378136,0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown prime meridian, numeric code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto pm = PrimeMeridian::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "My meridian"), + Angle(10)); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + Ellipsoid::WGS84, optional<std::string>(), pm); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true); + ASSERT_EQ(sql.size(), 5U); + EXPECT_EQ(sql[0], "INSERT INTO prime_meridian VALUES('HOBU','1','My " + "meridian',10,'EPSG','9122',0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Geographic 2D CRS with unknown prime meridian, alpha code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto pm = PrimeMeridian::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "My meridian"), + Angle(10)); + const auto datum = GeodeticReferenceFrame::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my datum"), + Ellipsoid::WGS84, optional<std::string>(), pm); + const auto crs = GeographicCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my EPSG:4326"), + datum, GeographicCRS::EPSG_4326->coordinateSystem()); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false); + ASSERT_EQ(sql.size(), 5U); + EXPECT_EQ(sql[0], "INSERT INTO prime_meridian " + "VALUES('HOBU','PM_GEODETIC_DATUM_XXXX','My " + "meridian',10,'EPSG','9122',0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Projected CRS, numeric code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto crs = ProjectedCRS::create( + PropertyMap().set(IdentifiedObject::NAME_KEY, "my projected CRS"), + GeographicCRS::EPSG_4807, + Conversion::createUTM(PropertyMap(), 31, true), + CartesianCS::createEastingNorthing(UnitOfMeasure::METRE)); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true); + ASSERT_EQ(sql.size(), 4U); + EXPECT_EQ(sql[0], + "INSERT INTO conversion VALUES('HOBU','1'," + "'UTM zone 31N',''," + "'EPSG','9807','Transverse Mercator'," + "'EPSG','8801','Latitude of natural origin',0,'EPSG','9122'," + "'EPSG','8802','Longitude of natural origin',3,'EPSG','9122'," + "'EPSG','8805','Scale factor at natural origin',0.9996," + "'EPSG','9201'," + "'EPSG','8806','False easting',500000,'EPSG','9001'," + "'EPSG','8807','False northing',0,'EPSG','9001'," + "NULL,NULL,NULL,NULL,NULL,NULL," + "NULL,NULL,NULL,NULL,NULL,NULL,0);"); + EXPECT_EQ(sql[1], + "INSERT INTO usage " + "VALUES('HOBU','USAGE_CONVERSION_1','conversion','HOBU','1','" + "PROJ','EXTENT_UNKNOWN','PROJ','SCOPE_UNKNOWN');"); + EXPECT_EQ( + sql[2], + "INSERT INTO projected_crs VALUES('HOBU','XXXX','my projected " + "CRS','','EPSG','4400','EPSG','4807','HOBU','1',NULL,0);"); + EXPECT_EQ( + sql[3], + "INSERT INTO usage " + "VALUES('HOBU','USAGE_PROJECTED_CRS_XXXX','projected_crs','HOBU','" + "XXXX','PROJ','EXTENT_UNKNOWN','PROJ','SCOPE_UNKNOWN');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Vertical CRS, known vertical datum, numeric code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + PropertyMap propertiesVDatum; + propertiesVDatum.set(Identifier::CODESPACE_KEY, "EPSG") + .set(Identifier::CODE_KEY, 5101) + .set(IdentifiedObject::NAME_KEY, "Ordnance Datum Newlyn"); + auto vdatum = VerticalReferenceFrame::create(propertiesVDatum); + PropertyMap propertiesCRS; + propertiesCRS.set(IdentifiedObject::NAME_KEY, "my height"); + const auto uom = + UnitOfMeasure("my unit", 3.0, UnitOfMeasure::Type::LINEAR); + const auto crs = VerticalCRS::create( + propertiesCRS, vdatum, VerticalCS::createGravityRelatedHeight(uom)); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true); + ASSERT_EQ(sql.size(), 5U); + EXPECT_EQ(sql[0], "INSERT INTO coordinate_system VALUES" + "('HOBU','CS_VERTICAL_CRS_XXXX','vertical',1);"); + EXPECT_EQ(sql[1], "INSERT INTO unit_of_measure VALUES" + "('HOBU','MY_UNIT','my unit','length',3,NULL,0);"); + EXPECT_EQ(sql[2], + "INSERT INTO axis VALUES('HOBU'," + "'CS_VERTICAL_CRS_XXXX_AXIS_1','Gravity-related height','H'," + "'up','HOBU','CS_VERTICAL_CRS_XXXX',1,'HOBU','MY_UNIT');"); + EXPECT_EQ(sql[3], + "INSERT INTO vertical_crs VALUES('HOBU','XXXX','my height'," + "'','HOBU','CS_VERTICAL_CRS_XXXX','EPSG','5101',0);"); + EXPECT_EQ(sql[4], + "INSERT INTO usage VALUES('HOBU','USAGE_VERTICAL_CRS_XXXX'," + "'vertical_crs','HOBU','XXXX','PROJ','EXTENT_UNKNOWN'," + "'PROJ','SCOPE_UNKNOWN');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", true).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Vertical CRS, unknown vertical datum, alpha code + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + PropertyMap propertiesVDatum; + propertiesVDatum.set(IdentifiedObject::NAME_KEY, "my datum"); + auto vdatum = VerticalReferenceFrame::create(propertiesVDatum); + PropertyMap propertiesCRS; + propertiesCRS.set(IdentifiedObject::NAME_KEY, "my height"); + const auto crs = VerticalCRS::create( + propertiesCRS, vdatum, + VerticalCS::createGravityRelatedHeight(UnitOfMeasure::METRE)); + const auto sql = + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false); + ASSERT_EQ(sql.size(), 4U); + EXPECT_EQ(sql[0], + "INSERT INTO vertical_datum VALUES('HOBU'," + "'VERTICAL_DATUM_XXXX','my datum','',NULL,NULL,NULL,0);"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), IComparable::Criterion::EQUIVALENT)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE( + ctxt->getInsertStatementsFor(crs, "HOBU", "XXXX", false).empty()); + ctxt->stopInsertStatementsSession(); + } + + // Compound CRS + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto wkt = + "COMPD_CS[\"unknown\"," + "PROJCS[\"NAD_1983_2011_StatePlane_South_Carolina_FIPS_3900_USFT\"," + "GEOGCS[\"NAD83(2011)\"," + "DATUM[\"NAD83_National_Spatial_Reference_System_2011\"," + "SPHEROID[\"GRS 1980\",6378137,298.257222101004," + "AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"1116\"]]," + "PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433," + "AUTHORITY[\"EPSG\",\"9122\"]]]," + "PROJECTION[\"Lambert_Conformal_Conic_2SP\"]," + "PARAMETER[\"latitude_of_origin\",31.8333333333333]," + "PARAMETER[\"central_meridian\",-81]," + "PARAMETER[\"standard_parallel_1\",32.5]," + "PARAMETER[\"standard_parallel_2\",34.8333333333333]," + "PARAMETER[\"false_easting\",1999996]," + "PARAMETER[\"false_northing\",0]," + "UNIT[\"US survey foot\",0.304800609601219," + "AUTHORITY[\"EPSG\",\"9003\"]]," + "AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]]," + "VERT_CS[\"NAVD88 height (ftUS)\"," + "VERT_DATUM[\"North American Vertical Datum 1988\",2005," + "AUTHORITY[\"EPSG\",\"5103\"]]," + "UNIT[\"US survey foot\",0.304800609601219," + "AUTHORITY[\"EPSG\",\"9003\"]]," + "AXIS[\"Up\",UP],AUTHORITY[\"EPSG\",\"6360\"]]]"; + const auto crs = + nn_dynamic_pointer_cast<CRS>(WKTParser().createFromWKT(wkt)); + ASSERT_TRUE(crs != nullptr); + const auto sql = ctxt->getInsertStatementsFor(NN_NO_CHECK(crs), "HOBU", + "XXXX", false); + ASSERT_EQ(sql.size(), 6U); + EXPECT_EQ(sql[4], + "INSERT INTO compound_crs VALUES('HOBU','XXXX','unknown'," + "'','HOBU','COMPONENT_XXXX_1','EPSG','6360',0);"); + EXPECT_EQ(sql[5], + "INSERT INTO usage VALUES('HOBU','USAGE_COMPOUND_CRS_XXXX'," + "'compound_crs','HOBU','XXXX','PROJ','EXTENT_UNKNOWN'," + "'PROJ','SCOPE_UNKNOWN');"); + const auto identified = + crs->identify(AuthorityFactory::create(ctxt, std::string())); + ASSERT_EQ(identified.size(), 1U); + EXPECT_EQ( + *(identified.front().first->identifiers().front()->codeSpace()), + "HOBU"); + EXPECT_TRUE(identified.front().first->isEquivalentTo( + crs.get(), + IComparable::Criterion::EQUIVALENT_EXCEPT_AXIS_ORDER_GEOGCRS)); + EXPECT_EQ(identified.front().second, 100); + EXPECT_TRUE(ctxt->getInsertStatementsFor(NN_NO_CHECK(crs), "HOBU", + "XXXX", false) + .empty()); + ctxt->stopInsertStatementsSession(); + } + + // DynamicGeodeticReferenceFrame + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto datum = AuthorityFactory::create(ctxt, "EPSG") + ->createDatum("1165"); // ITRF2014 + const auto sql = ctxt->getInsertStatementsFor(datum, "HOBU", "XXXX", + false, {"HOBU"}); + const auto datumNew = + AuthorityFactory::create(ctxt, "HOBU")->createDatum("XXXX"); + EXPECT_TRUE(datumNew->isEquivalentTo( + datum.get(), IComparable::Criterion::EQUIVALENT)); + ctxt->stopInsertStatementsSession(); + } + + // DynamicVerticalReferenceFrame + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto datum = AuthorityFactory::create(ctxt, "EPSG") + ->createDatum("1096"); // Norway Normal Null 2000 + const auto sql = ctxt->getInsertStatementsFor(datum, "HOBU", "XXXX", + false, {"HOBU"}); + const auto datumNew = + AuthorityFactory::create(ctxt, "HOBU")->createDatum("XXXX"); + EXPECT_TRUE(datumNew->isEquivalentTo( + datum.get(), IComparable::Criterion::EQUIVALENT)); + ctxt->stopInsertStatementsSession(); + } + + // geodetic DatumEnsemble, and add members inline + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto ensemble = AuthorityFactory::create(ctxt, "EPSG") + ->createDatumEnsemble("6326"); // WGS84 + const auto sql = ctxt->getInsertStatementsFor(ensemble, "HOBU", "XXXX", + false, {"HOBU"}); + const auto ensembleNew = + AuthorityFactory::create(ctxt, "HOBU")->createDatumEnsemble("XXXX"); + EXPECT_TRUE(ensembleNew->isEquivalentTo( + ensemble.get(), IComparable::Criterion::EQUIVALENT)); + ctxt->stopInsertStatementsSession(); + } + + // geodetic DatumEnsemble, and reference members + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + const auto ensemble = AuthorityFactory::create(ctxt, "EPSG") + ->createDatumEnsemble("6326"); // WGS84 + const auto sql = + ctxt->getInsertStatementsFor(ensemble, "HOBU", "XXXX", false); + const auto ensembleNew = + AuthorityFactory::create(ctxt, "HOBU")->createDatumEnsemble("XXXX"); + EXPECT_TRUE(ensembleNew->isEquivalentTo( + ensemble.get(), IComparable::Criterion::EQUIVALENT)); + ctxt->stopInsertStatementsSession(); + } + + // vertical DatumEnsemble + { + auto ctxt = DatabaseContext::create(); + ctxt->startInsertStatementsSession(); + // British Isles height ensemble + const auto ensemble = + AuthorityFactory::create(ctxt, "EPSG")->createDatumEnsemble("1288"); + const auto sql = ctxt->getInsertStatementsFor(ensemble, "HOBU", "XXXX", + false, {"HOBU"}); + const auto ensembleNew = + AuthorityFactory::create(ctxt, "HOBU")->createDatumEnsemble("XXXX"); + EXPECT_TRUE(ensembleNew->isEquivalentTo( + ensemble.get(), IComparable::Criterion::EQUIVALENT)); + ctxt->stopInsertStatementsSession(); + } +} + } // namespace |