Create a tessellation from a subset of sampled data.

Sample_Tessellation.R

# create a tessalation
sampleVoronoi <- function(x, sampled,rw=NULL) {

require(raster)

  vor <- tile.list( deldir(x[sampled,1], x[sampled,2],rw=rw) )

  vorPoly <- lapply( 1:length(vor) ,
    function(i) {
      tmp <- cbind( vor[[i]]$x, vor[[i]]$y )
      tmp <- rbind( tmp, tmp[1,])
      colnames(tmp) <- c('x','y')
      return( Polygons(list(Polygon(tmp)), ID=i) )
    }
  )

  vorDat <- x[sampled,]
  vorDat <- as.data.frame(vorDat)

  rownames(vorDat) <- sapply(slot(SpatialPolygons(vorPoly), 'polygons'), slot,'ID')

  return( SpatialPolygonsDataFrame(SpatialPolygons(vorPoly),data=vorDat) )
}