2

I've read the Overlapping LAS tiles post, but it doesn't pertain to the following situation.

I'm trying to process a las catalog where I clip the catalog at specific plot locations with a 30 m radius. enter image description here

As a result of the 30 m radius, the clipped tiles overlap and I see edge artifacts in the output (as indicated might happen in the readLAScatalog warning).
enter image description here

Is it possible to process these tiles as a catalog without edge artifacts? Or will I need to read each file individually using readLAS? I've tried processing the catalog sequentially, but the output still has edge effects.

1 Answer 1

2

Your case is covered in the book in this section. The idea here is that your files are no longer a wall-to-wall coverage of a territory. The LAScatalog engine is designed to process las files that encompass a territory where each files is spatially related to its neighbourhood.

In your case the files are independent. You basically want to loop through each file. In this case you can put this information in the catalog

opt_independent_files(ctg) <- TRUE

Each file will be processed independently, without buffer. This is equivalent to a regular loop (but with parallel capabilities and real time monitoring).

In addition depending of what you are doing you may want to avoid final merging

opt_merge(ctg) = FALSE
2
  • That's great and thanks for link!
    – derelict
    Feb 11, 2021 at 18:40
  • 1
    I edited my answer to show you how to do not merge the output
    – JRR
    Feb 11, 2021 at 19:21

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.