makestatkde: makestatkde

In ISAAKiel/pathAAR: Pathes in Archaeology

Description

Calculating Static KDE

Usage

makestatkde(pppm, f_sd1 = 4, sgdf, df, x = 1, y = 2, num)

Arguments

pppm	PointPattern object of Spatial Data
f_sd1	factor defining size of the first kernel, which generate the structure of dynamic kernel
sgdf	grid object to store results
df	data frame, containing coordinates of monuments
x	numeric, indicating column number of x coordinates
y	numeric, indicating column number of y coordinates
num	numeric, number of grid cells

Details

description of function

Value

SpatialGrid Object

Author(s)

Franziska Faupel <[e-mail]>

Oliver Nakoinz <[e-mail]>

Hendrik Raese <[e-mail]>

Examples

# Crating Test Data Randomly
testmatrix <- data.frame(x = abs(rnorm(100)*50), y = abs(rnorm(100)*50))
# Setting geographical frame
xmin    <- 0
xmax    <- max(testmatrix$x)
ymin    <- 0
ymax    <- max(testmatrix$y)
ext_ai <- raster::extent(xmin, xmax, ymin, ymax)

sv <- (xmax-xmin)/(ymax-ymin)
rw     <- 10   # width of raster defined in m

rows  <- round((ymax-ymin)/rw, 0) + 1                                    
colums <- round((xmax-xmin)/rw, 0) + 1                                      
v <- cbind(1:(colums*rows))                                              
df <- data.frame(v)                                                         
gt      <- sp::GridTopology(c(xmin, ymin), c(rw, rw), c(colums, rows))
sgdf    <- sp::SpatialGridDataFrame(gt, df)

pppm <- spatstat::ppp(testmatrix[,1], testmatrix[,2], 
                       window = spatstat::owin(
                       xrange=c(sgdf@bbox[1,1],sgdf@bbox[1,2]),
                       yrange=c(sgdf@bbox[2,1],sgdf@bbox[2,2]), 
                       unitname="m"))
                       
base_kde <- makestatkde(pppm, sgdf=sgdf, df=testmatrix, x=1, y=2, num=length(df[,1]))

ISAAKiel/pathAAR documentation built on Sept. 7, 2021, noon