Add WKT function

fmidev · Jan 5, 2024 · 5e1e807 · 5e1e807
1 parent f0d38bc
commit 5e1e807
Show file tree

Hide file tree

Showing 2 changed files with 29 additions and 0 deletions.
diff --git a/himan-lib/include/grid.h b/himan-lib/include/grid.h
@@ -255,6 +255,7 @@ class regular_grid : public grid
 	virtual double Dj() const;
 
 	virtual std::string Proj4String() const;
+	virtual std::string WKT(const std::map<std::string, std::string>& opts = {}) const;
 	virtual std::unique_ptr<OGRPolygon> Geometry() const;
 	virtual earth_shape<double> EarthShape() const override;
 

diff --git a/himan-lib/source/grid.cpp b/himan-lib/source/grid.cpp
@@ -229,6 +229,34 @@ std::string regular_grid::Proj4String() const
 	return proj;
 }
 
+std::string regular_grid::WKT(const std::map<std::string, std::string>& opts) const
+{
+	if (itsSpatialReference == nullptr)
+	{
+		throw std::runtime_error("Spatial reference instance not initialized");
+	}
+
+	char** projopts = nullptr;
+
+	// https://gdal.org/doxygen/classOGRSpatialReference.html#ae986da88649783b5c194de55c46890a5
+	for (const auto& opt : opts)
+	{
+		projopts = CSLAddNameValue(projopts, opt.first.c_str(), opt.second.c_str());
+	}
+
+	char* wkt;
+	if (itsSpatialReference->exportToWkt(&wkt) != OGRERR_NONE)
+	{
+		throw std::runtime_error("Failed to get wkt");
+	}
+
+	std::string w(wkt);
+	CPLFree(wkt);
+	CSLDestroy(projopts);
+
+	return w;
+}
+
 point regular_grid::Projected(const point& latlon) const
 {
 	double projX = latlon.X(), projY = latlon.Y();