require(raster)
require(sp)
r <- raster(ncol=1000, nrow=1000)
r[] <- rpois(ncell(r), lambda=1)
r <- calc(r, fun=function(x) { x[x >= 1] <- 1; return(x) } )
xy <- SpatialPointsDataFrame(SpatialPoints(cbind(-50, seq(-80, 80, by=20))),
data.frame(ID=seq(1,9,1)))
land.metrics(x=xy, y=r, bw=0.72, bkgd = 0, Trace=FALSE)Run the code above in your browser using DataLab