Calculates ambient occlusion using sky visibility algorithm. Based on Duguet & Girardeau-Montaut (2004).
compute_pcv(las, radius = 1, num_directions = 60, ncpu = 4)Numeric vector of PCV values (sky visibility) for each point
LAS object
Numeric radius for neighborhood search
Integer number of directional rays to cast
Integer number of CPUs to use for parallel processing
The PCV algorithm computes the visible portion of the sky from each point by:
Casting rays in multiple directions around each point
Computing the maximum elevation angle (horizon) in each direction
Calculating the average sky visibility across all directions
This function uses an optimized C++ implementation with OpenMP parallelization for improved performance on large point clouds.