sort based on GPStime for COPC output (#94)

hobuinc · Nov 29, 2021 · 6421a7f · abellgithub · Jan 6, 2022 · hobu
1 parent c28d537
commit 6421a7f
Show file tree

Hide file tree

Showing 2 changed files with 19 additions and 1 deletion.
diff --git a/bu/Processor.cpp b/bu/Processor.cpp
@@ -450,7 +450,7 @@ void Processor::writeEptFile(const std::string& filename, pdal::PointTableRef ta
     pdal::Options wopts;
     wopts.add("extra_dims", "all");
     wopts.add("software_id", "Entwine 1.0");
-    wopts.add("compression", "laszip");
+    wopts.add("compression", "lazperf");
     wopts.add("filename", filename);
     wopts.add("offset_x", m_b.offset[0]);
     wopts.add("offset_y", m_b.offset[1]);
@@ -470,6 +470,19 @@ void Processor::writeEptFile(const std::string& filename, pdal::PointTableRef ta
     w->execute(table);
 }
 
+void Processor::sortChunk(const VoxelKey& key, pdal::PointViewPtr view)
+{
+
+    using namespace pdal;
+    auto cmp = [](const PointRef& p1, const PointRef& p2)
+    {
+        bool result = p1.compare(pdal::Dimension::Id::GpsTime, p2);
+        return result;
+    };
+
+    std::stable_sort(view->begin(), view->end(), cmp);
+
+}
 void Processor::createChunk(const VoxelKey& key, pdal::PointViewPtr view)
 {
     using namespace pdal;
@@ -480,7 +493,11 @@ void Processor::createChunk(const VoxelKey& key, pdal::PointViewPtr view)
         return;
     }
 
+    if (view->layout()->hasDim(Dimension::Id::GpsTime))
+        sortChunk(key, view);
+
     PointLayoutPtr layout = view->layout();
+
     int ebCount {0};
     for (DimType dim : m_extraDims)
         ebCount += layout->dimSize(dim.m_id);

diff --git a/bu/Processor.hpp b/bu/Processor.hpp
@@ -59,6 +59,7 @@ class Processor
     void writeEptFile(const std::string& filename, pdal::PointTableRef table,
         pdal::PointViewPtr view);
     void createChunk(const VoxelKey& key, pdal::PointViewPtr view);
+    void sortChunk(const VoxelKey& key, pdal::PointViewPtr view);
     void fillPointBuf(pdal::PointRef& point, std::vector<char>& buf);
 
     VoxelInfo m_vi;