My original catalog:
has a mean density of 2.2 points/m2, but some areas have a higher density. :
> plot(lidR::grid_density(ctg, 1), main = "Density grid 1x1 m ", breaks = c(0,2,3,5,16), col = topo.colors(4))
I want to homogenize the point cloud to the maximum target density of 3 points/m2. To do so I use the function
decimate_points with the algorithm
ctg_thinned = decimate_points(ctg_clip, homogenize(density = 3, res = 25, use_pulse = FALSE)) > plot(lidR::grid_density(ctg_thinned , 1), breaks = c(0,2,3,5,16), col = topo.colors(4))
How is it possible that I get higher densities than 3 and even more, higher maximum densities than before?
Maybe is not the correct algortihm I am using. If not, how should I do I do it?