R/density.R

In n3ssuno/RKS: Regional Knowledge Space

#' @export

# - improve documentation

local_density <- function(occt, FUN = proximity) {
    #-- Preliminary steps and checks
    info <- data_info(occt)
    occt <- remove_zeros(occt)

    #-- Main function
    if (info$n_dims == 3) {
        w <- apply(occt, 3, local_density)
    } else {
        Phi <- FUN(occt)
        den <- Matrix::rowSums(Phi, na.rm = TRUE)
        Mcp <- rta(occt, binary = TRUE)
        w <- Matrix::tcrossprod(Mcp, Phi / den)
    }

    return(w)
}