if (FALSE) {
mba <- function(n = 1, m = 1, h = 8, extend = TRUE) {
f <- function(las, where) {
res <- MBA::mba.points(las@data, where, n, m , h, extend)
return(res$xyz.est[,3])
}
f <- plugin_dtm(f)
return(f)
}
LASfile <- system.file("extdata", "Topography.laz", package="lidR")
las = readLAS(LASfile)
dtm = rasterize_terrain(las, algorithm = mba())
}
Run the code above in your browser using DataLab