R/utils.R

In hanel/cordex: What the package does (short line)

hasDate = function(x){
  if (!inherits(x, 'Raster') ) stop(parse(text = x), ' is not a raster object.')
  if (is.null(x@z)) stop(parse(text = x), 'does not have z dimension.')
}

month = function(...){
  UseMethod('month')
}

month.default = data.table::month

month.RasterBrick = function(x) {
  hasDate(x)
  month(x@z[[1]])
}

month.RasterStack = month.RasterBrick

year = function(...){
  UseMethod('year')
}

year.default = data.table::year

year.RasterBrick = function(x) {
  hasDate(x)
  year(x@z[[1]])
}

monmean = function(x){
  
  if (is.null(x)) return(NULL)
  
  hasDate(x)
  mon = month(x)
  
  R = list()
  for (i in 1:12){
    rr = x[[which(mon == i)]]
    R[[i]] = mean(rr)
  }
  R = brick(R)
  R@z = list(1:12)
  names(R) = month.abb
  R
}

selper = function(x, per){
  
  y = year(x)
  if (length(which(y %in% per)) == 0) return(NULL)
  px = x[[which(y %in% per)]]
  px@z = list(Date = x@z[[1]][(y %in% per)])
 # names(px) = as.character(y[(y %in% per)])
  px
  
}

wgs2krov=function(wgs){
  
  wgs@proj4string=CRS('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs')
  spTransform(wgs,CRS('+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=0 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs'))

}

krov2wgs=function(krov){
  
  require(rgdal)
  krov@proj4string=CRS('+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=0 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs')
  spTransform(krov,CRS('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'))
  
}