Extent no longer dependent of the subcircle option #518

NEWS.md

@@ -109,6 +109,10 @@ New functions are mostly convenient features that simplify some workflow without
 
 15. New functions `header()`, `payload()`, `phb()`, `vlr()`, `evlr()` to get the corresponding data from a `LAS` object.
 
+### Enhancement
+
+1. The bounding box of the CHM not longer depends on the `subcircle` tweak. See [#518](https://github.com/r-lidar/lidR/issues/518).
+
 ### Documentation
 
 - Man pages of `classify_*`, `rasterize_*`, `*_metrics`, `segment_*` and `normalize_*` were grouped.

R/rasterize_canopy.R

@@ -16,13 +16,8 @@ rasterize_canopy.LAS = function(las, res = 1, algorithm = p2r(), ...)
   dots <- list(...)
   pkg <- if (is.null(dots$pkg)) getOption("lidR.raster.default") else dots$pkg
 
-  # Some algorithm have an extra option 'subscircle' that need to buffer the layout
-  # Must be rewritten because it is a hack !
-  subcircle <- as.list(environment(algorithm))$subcircle
-  subcircle <- if (is.null(subcircle)) 0 else subcircle
-
   # Compute the raster that encompass the point cloud
-  layout <- if (!is_a_number(res)) raster_template(res) else raster_layout(las, res, buffer = subcircle, format = "template")
+  layout <- if (!is_a_number(res)) raster_template(res) else raster_layout(las, res, format = "template")
   layout <- raster_materialize(layout, pkg = "stars")
 
   # Compute the elevation for each cells
@@ -32,7 +27,7 @@ rasterize_canopy.LAS = function(las, res = 1, algorithm = p2r(), ...)
   # Quantize
   z <- round(z, 3)
 
-  layout <- raster_layout(las, res, buffer = subcircle)
+  layout <- raster_layout(las, res)
   layout <- raster_materialize(layout, pkg = pkg)
   layout <- raster_set_values(layout, z)
   raster_names(layout) <- "Z"