R/makeRaster.R

In KNMI-DataLab/SkyViewFactor: Calculate the Sky view factor from las files

#'Create a raster
#'@title Make a raster
#'@param spatialDF a spatialpointsdataframe object
#'@param xres resolution in the x direction
#'@param yres resolution in the y-direction
#'@param pro the projection of the raster
#'@export
#'
makeRaster<-function(spatialDF, xres, yres, pro){
  dummyRaster<-raster(nrow=10,ncol=10,crs=pro) #dummy raster with projection
  extent(dummyRaster)<-extent(spatialDF)
  res(dummyRaster)<-c(xres,yres) # set the resolution
  r<-rasterize(spatialDF,dummyRaster,field="Z") #rasterizing the spatial points to a 1x1 grid
  rm(dummyRaster)
  gc()
  return(r)
}