R/expandedge.R

In forestSAS: Forest Spatial Structure Analysis Systems

expandedge <-
function(X,xwidth,ywidth,id=1:X$n){
  stopifnot(spatstat::is.ppp(X))
  x<-y<-NULL
  df.X<-as.data.frame(X)
  df.X.ppp<-spatstat::ppp(x=df.X$x,y=df.X$y,
                          xrange=X$window$xrange,yrange=X$window$yrange)
  if(is.null(X$marks)){
    df.X.ppp<-df.X.ppp
  }else{
    spatstat::marks(df.X.ppp)<-spatstat::marks(df.X.ppp)
  }
  if(xwidth>df.X.ppp$window$xrange[2]-df.X.ppp$window$xrange[1])
    stop("xwidth beyond the xrange")
  if(ywidth>df.X.ppp$window$yrange[2]-df.X.ppp$window$yrange[1])
    stop("ywidth beyond the yrange")
  totalX<-as.data.frame(cbind(x=df.X.ppp$x,y=df.X.ppp$y,df.X.ppp$marks,id))
  s.totalX<-totalX
  s.totalX$y<-s.totalX$y+ywidth
  uparea<-subset(totalX,y<=df.X.ppp$window$yrange[1]+ywidth)
  uparea$y<-uparea$y+(df.X.ppp$window$yrange[2]-df.X.ppp$window$yrange[1])
  downarea<-subset(totalX,y>=df.X.ppp$window$yrange[2]-ywidth)
  downarea$y<-downarea$y-(df.X.ppp$window$yrange[2]-df.X.ppp$window$yrange[1])
  updown_area<-rbind(uparea,totalX,downarea)
  s.updown_area<-updown_area
  rightarea<-subset(updown_area,x<=df.X.ppp$window$xrange[1]+xwidth)
  rightarea$x<-rightarea$x+(df.X.ppp$window$xrange[2]-df.X.ppp$window$xrange[1])
  leftarea<-subset(updown_area,x>=df.X.ppp$window$xrange[2]-xwidth)
  leftarea$x<-leftarea$x-(df.X.ppp$window$xrange[2]-df.X.ppp$window$xrange[1])
  leftright_area<-data.frame(rbind(leftarea,s.updown_area,rightarea))
  total_area<-as.data.frame(cbind(id=1:(nrow(leftarea)+nrow(rightarea)+nrow(updown_area)),old.id=leftright_area$id,leftright_area[-which(colnames(leftright_area)=="id")],row.names=NULL))
}