localMax: localMax

In ISAAKiel/pathARR: Pathes in Archaeology

Description

Local Centres of Infrastructure from Monument Location

Usage

localMax(df, r = 5000, sw = 10, pd = 500)

Arguments

df	data frame, containing coordinates of path associated features
r	numeric, sets the size of moving window and sigma in the used density function in meter, default = 5000
sw	numeric, value used for defining threshold for defining a local maximum. A higher value means more local maxima are stored, default = 10. It's suggested not to change this parameter and instead work with the size of the moving window!
pd	numeric, dimension of pixels for intermediate raster image of the used density function in meter, default = 500

Details

In absence of known contemporary nodes, like settlements, the density of monument locations can be used to reconstruct nodes of infrastructure. This function will use the Kernel Density Approach combined with a moving window, in order to detect nodes in areas with less monuments (from a global perspective).

Value

SpatialPointsDataFrame with coordinates of local maxima

Author(s)

Franziska Faupel <[e-mail]>

Oliver Nakoinz <[e-mail]>

Hendrik Raese <[e-mail]>

Examples

set.seed(123)

# Creating random test data 
testmatrix <- data.frame(
 x = abs(runif(100, 1000, 1500)), 
 y = abs(runif(100, 1000, 1500)))

pd=25
sw=10
r=100
maxima <- localMax(df=testmatrix, r=r, sw=sw, pd=pd)

# Plot the result
maxima <- data.frame(maxima)

ggplot2::ggplot() +
 ggplot2::geom_point(data= testmatrix, ggplot2::aes(x=x, y=y)) +
 ggplot2::geom_point(data=maxima, 
                     ggplot2::aes(x=mx, y=my, color="red"), 
                     shape=15, 
                     size=3, 
                     show.legend = FALSE)

ISAAKiel/pathARR documentation built on Sept. 10, 2021, 7:54 a.m.