proj_create_vertical_crs_ex(): add a ACCURACY option to provide an explicit accuracy, or derive it from the grid name if it is known

1 files changed, 19 insertions, 0 deletions

diff --git a/src/iso19111/factory.cpp b/src/iso19111/factory.cpp
index a9d82268..a011f397 100644
--- a/src/iso19111/factory.cpp
+++ b/src/iso19111/factory.cpp
@@ -1256,6 +1256,25 @@ DatabaseContext::getNonDeprecated(const std::string &tableName,
     return res;
 }
 
+// ---------------------------------------------------------------------------
+
+std::vector<operation::CoordinateOperationNNPtr>
+DatabaseContext::getTransformationsForGridName(
+    const DatabaseContextNNPtr &databaseContext, const std::string &gridName) {
+    auto sqlRes = databaseContext->d->run(
+        "SELECT auth_name, code FROM grid_transformation "
+        "WHERE grid_name = ? OR grid_name = "
+        "(SELECT original_grid_name FROM grid_alternatives "
+        "WHERE proj_grid_name = ?)",
+        {gridName, gridName});
+    std::vector<operation::CoordinateOperationNNPtr> res;
+    for (const auto &row : sqlRes) {
+        res.emplace_back(AuthorityFactory::create(databaseContext, row[0])
+                             ->createCoordinateOperation(row[1], true));
+    }
+    return res;
+}
+
 //! @endcond
 
 // ---------------------------------------------------------------------------