Fix computation of voxel indices (#1901)

PDAL · Apr 6, 2018 · ed68fb5 · ed68fb5
1 parent a7046ff
commit ed68fb5
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 6 deletions.
diff --git a/filters/VoxelCenterNearestNeighborFilter.cpp b/filters/VoxelCenterNearestNeighborFilter.cpp
@@ -79,9 +79,9 @@ PointViewSet VoxelCenterNearestNeighborFilter::run(PointViewPtr view)
         double y = view->getFieldAs<double>(Dimension::Id::Y, id);
         double x = view->getFieldAs<double>(Dimension::Id::X, id);
         double z = view->getFieldAs<double>(Dimension::Id::Z, id);
-        size_t r = static_cast<size_t>(floor(y - bounds.miny) / m_cell);
-        size_t c = static_cast<size_t>(floor(x - bounds.minx) / m_cell);
-        size_t d = static_cast<size_t>(floor(z - bounds.minz) / m_cell);
+        size_t r = static_cast<size_t>((y - bounds.miny) / m_cell);
+        size_t c = static_cast<size_t>((x - bounds.minx) / m_cell);
+        size_t d = static_cast<size_t>((z - bounds.minz) / m_cell);
         populated_voxels.emplace(std::make_tuple(r, c, d));
     }
 

diff --git a/filters/VoxelCentroidNearestNeighborFilter.cpp b/filters/VoxelCentroidNearestNeighborFilter.cpp
@@ -83,9 +83,9 @@ PointViewSet VoxelCentroidNearestNeighborFilter::run(PointViewPtr view)
         double y = view->getFieldAs<double>(Dimension::Id::Y, id);
         double x = view->getFieldAs<double>(Dimension::Id::X, id);
         double z = view->getFieldAs<double>(Dimension::Id::Z, id);
-        size_t r = static_cast<size_t>(floor(y - bounds.miny) / m_cell);
-        size_t c = static_cast<size_t>(floor(x - bounds.minx) / m_cell);
-        size_t d = static_cast<size_t>(floor(z - bounds.minz) / m_cell);
+        size_t r = static_cast<size_t>((y - bounds.miny) / m_cell);
+        size_t c = static_cast<size_t>((x - bounds.minx) / m_cell);
+        size_t d = static_cast<size_t>((z - bounds.minz) / m_cell);
         populated_voxel_ids[std::make_tuple(r, c, d)].push_back(id);
     }