R/possess.R

In Schuch666/inventory: A R-package for modeling emission inventories from geographic data

#' Read a georeference image and extract values for using in griding
#'
#' @description extract values from a image to be used on each region
#'
#' @param geoemiss a sf from geoemiss
#' @param filename image filename
#' @param plots plot individual regions
#' @param verbose display additional information
#'
#' @return a sf with the area shapes and the total emission and a raster image for each area
#'
#' @import sf
#' @import raster
#'
#' @export
#'
#' @examples
#' states <- sf::st_read(paste(system.file("extdata", package = "inventory"),"/states.shp",sep=""),
#'                       quiet = TRUE)
#' Nox    <- geoemiss(geom = states,variable = "Nox",names = c("sp","rj"),values = c(1000,20))
#' image  <- paste(system.file("extdata", package = "inventory"),"/tiny.tif",sep="")
#' ras    <- possess(Nox,image,plots = FALSE)
#' test   <- griding(ras,variable = "Nox", res = 0.1, type = "local", plot = TRUE)


possess <- function(geoemiss, filename = NA, plots = T,verbose = T){
  if(is.na(filename)){
    cat("choose a file:\n")   # nocov
    filename <- file.choose() # nocov
    cat(paste(filename,"\n")) # nocov
  }

  r  <- raster::crop(raster::raster(filename),sf::as_Spatial(geoemiss$geometry))
  s  <- sf::as_Spatial(geoemiss$geometry)
  l  <- list()
  geoemiss$image <- NULL

  for(j in 1:nrow(geoemiss)){
    cro <- raster::mask(r,s[j])
    cro <- raster::crop(cro,s[j])
    if(plots)
      plot(cro,main =  paste(geoemiss$region[j])) # nocov
    l[[j]] <- cro
  }

  geoemiss$image <- l
  return(geoemiss)
}