0

I read in my las catalog,

then chunked it using following code -

dir.create(glue('{2023}/tile'), showWarnings = FALSE)
opt_chunk_buffer(ctg) <- 10
opt_chunk_size(ctg) <- 100
opt_output_files(ctg) <- paste0(2023, "/tile/{XLEFT}_{YBOTTOM}")
opt_laz_compression(ctg) <- T
plot(ctg, chunk = TRUE)
ctg_tiled <- catalog_retile(ctg)

Then I tried making the DTM with the TIN method on my initial ctg file and also the retiled file -

dtm_tin <- rasterize_terrain(ctg_tiled, res = 0.25, algorithm = tin())
dtm_tin <- rasterize_terrain(ctg, res = 0.25, algorithm = tin())

BOTH gives same error (with all gray chunks, no green ones) - Error in any_list[[1]] : subscript out of bounds

How to fix this? Did I do something wrong with the retiling or is there a step I missed? What does the error even mean?

0 Answers0