ensure WKTv1 writing for LAS 1.4

PDAL · Feb 21, 2020 · c4a673b · c4a673b
1 parent 366bd3b
commit c4a673b
Show file tree

Hide file tree

Showing 3 changed files with 25 additions and 0 deletions.
diff --git a/io/LasWriter.cpp b/io/LasWriter.cpp
@@ -507,6 +507,8 @@ bool LasWriter::addWktVlr()
     if (wkt.empty())
         return false;
 
+    // LAS 1.4 requires WKTv1
+    wkt = SpatialReference::convertToWKT1(wkt);
     std::vector<uint8_t> wktBytes(wkt.begin(), wkt.end());
     // This tacks a NULL to the end of the data, which is required by the spec.
     wktBytes.resize(wktBytes.size() + 1, 0);

diff --git a/pdal/SpatialReference.cpp b/pdal/SpatialReference.cpp
@@ -438,6 +438,28 @@ std::string SpatialReference::prettyWkt(const std::string& wkt)
     return outWkt;
 }
 
+std::string SpatialReference::convertToWKT1(const std::string& wkt)
+{
+    std::string outWkt(wkt);
+
+// WKT flavors only exist in GDAL 3.x+
+#if GDAL_VERSION_MAJOR >= 3
+    OGRScopedSpatialReference srs = ogrCreateSrs(wkt);
+    if (!srs)
+        return outWkt;
+
+    char *buf = nullptr;
+
+    char **papszOptions = nullptr;
+    papszOptions = CSLAddNameValue(papszOptions, "FORMAT", "WKT1_GDAL");
+    srs->exportToWkt(&buf, papszOptions);
+
+    outWkt = buf;
+    CPLFree(buf);
+#endif
+    return outWkt;
+}
+
 
 int SpatialReference::getUTMZone() const
 {

diff --git a/pdal/SpatialReference.hpp b/pdal/SpatialReference.hpp
@@ -164,6 +164,7 @@ class PDAL_DLL SpatialReference
     static int calculateZone(double lon, double lat);
     static bool isWKT(const std::string& wkt);
     static std::string prettyWkt(const std::string& wkt);
+    static std::string convertToWKT1(const std::string& wkt);
 
 private:
     std::string m_wkt;