update gridDecimation filter.

alavenant · Jun 5, 2024 · 3dd71eb · 3dd71eb
1 parent 02ff3e2
commit 3dd71eb
Showing 1 changed file with 45 additions and 3 deletions.
diff --git a/test/unit/filters/GridDecimationFilterTest.cpp b/test/unit/filters/GridDecimationFilterTest.cpp
@@ -12,6 +12,8 @@
 #include <io/FauxReader.hpp>
 #include <filters/GridDecimationFilter.hpp>
 
+#include <filters/CropFilter.hpp>
+
 #include "Support.hpp"
 
 using namespace pdal;
@@ -32,9 +34,11 @@ TEST(DecimationFilterTest, GridDecimationFilterTest_test1)
     Stage& r = *(factory.createStage("readers.las"));
     r.setOptions(ro);
 
+    double resolution = 10.;
+
     Options gdOps;
     gdOps.add("output_type", "max");
-    gdOps.add("resolution", 10.);
+    gdOps.add("resolution", resolution);
     gdOps.add("value", "Classification=5");
 
     GridDecimationFilter filter;
@@ -49,14 +53,52 @@ TEST(DecimationFilterTest, GridDecimationFilterTest_test1)
     EXPECT_EQ(viewSet.size(), 1u);
 
     PointViewPtr view = *viewSet.begin();
-
+    
     int nbThreadPts (0);
     for (PointId i = 0; i < view->size(); ++i)
     {
         PointRef point = view->point(i);
         uint8_t classification = point.getFieldAs<uint8_t>(Dimension::Id::Classification);
         if (classification==5) nbThreadPts++;
     }
-
+    
     EXPECT_EQ(nbThreadPts, 400);
+
+    double xmin = 1639600.0;
+    double ymin = 1454500.02;
+
+    for (size_t l(0); l<20; l++){
+        for (size_t c(0); l<20; l++){
+
+            BOX2D dstBounds(xmin + c*resolution, ymin + l*resolution, xmin + (c+1)*resolution, ymin + (l+1)*resolution);
+
+            Options cropOpts;
+            cropOpts.add("bounds", dstBounds);
+            CropFilter cropfilter;
+            cropfilter.setOptions(cropOpts);
+
+            PointTable cropTable;
+            cropfilter.prepare(cropTable);
+            cropfilter.setInput(filter);
+            viewSet = cropfilter.execute(table);
+            view = *viewSet.begin();
+
+            int nbThreadPtsCrop (0);
+            double Zmax, ZmaxGrid;
+            for (PointId i = 0; i < view->size(); ++i)
+            {
+                PointRef point = view->point(i);
+                double z_pt = point.getFieldAs<double>(Dimension::Id::Z);
+                if (i==0) Zmax = z_pt;
+                else if (z_pt > Zmax) Zmax = z_pt;
+                uint8_t classification = point.getFieldAs<uint8_t>(Dimension::Id::Classification);
+                if (classification==5) {nbThreadPtsCrop++; ZmaxGrid=z_pt;}
+            }
+
+            EXPECT_EQ(nbThreadPtsCrop, 1);
+            EXPECT_EQ(Zmax, ZmaxGrid);
+        }
+    }
+
+
 }