R/addmark.ppp.R

In forestSAS: Forest Spatial Structure Analysis Systems

addmark.ppp <-
function(X,add.mark,add.name="storey"){
  stopifnot(spatstat::is.ppp(X))
  df.X<-as.data.frame(X)
  df.addX<-data.frame(cbind(df.X,add.mark))
  names(df.addX)<-c(names(df.X),add.name)
  newX<-spatstat::ppp(X$x,X$y,window=X$window,
     marks=df.addX[,-match(c("x","y"),colnames(df.addX))])
  return(newX)
}