From fda68193f0ad0b02fa82135d2ac638d77a0101fa Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Fri, 19 Mar 2021 16:00:26 +0100 Subject: CRS::promoteTo3D(): propagate the extent from the 2D CRS (fixes #2587) This will help getting more consistent results between the 2D and 3D cases, as identified in https://github.com/OSGeo/PROJ/issues/2587#issue-836061171 --- src/iso19111/crs.cpp | 18 ++++++++++++++++++ 1 file changed, 18 insertions(+) (limited to 'src/iso19111/crs.cpp') diff --git a/src/iso19111/crs.cpp b/src/iso19111/crs.cpp index 9c4c04c3..20cf54c9 100644 --- a/src/iso19111/crs.cpp +++ b/src/iso19111/crs.cpp @@ -1000,6 +1000,24 @@ CRSNNPtr CRS::promoteTo3D(const std::string &newName, auto props = util::PropertyMap().set(common::IdentifiedObject::NAME_KEY, !newName.empty() ? newName : nameStr()); + const auto &l_domains = domains(); + if (!l_domains.empty()) { + auto array = util::ArrayOfBaseObject::create(); + for (const auto &domain : l_domains) { + auto extent = domain->domainOfValidity(); + if (extent) { + // Propagate only the extent, not the scope, as it might + // imply more that we can guarantee with the promotion to + // 3D. + auto newDomain = common::ObjectDomain::create( + util::optional(), extent); + array->add(newDomain); + } + } + if (!array->empty()) { + props.set(common::ObjectUsage::OBJECT_DOMAIN_KEY, array); + } + } const auto &l_identifiers = identifiers(); if (l_identifiers.size() == 1) { std::string remarks("Promoted to 3D from "); -- cgit v1.2.3