R/matrix2raster.r

In oceanmap: A Plotting Toolbox for 2D Oceanographic Data

matrix2raster <- function(z,x,y,layer,proj="+proj=longlat"){
#   #inst.pkg('raster')

  if(!is.matrix(z)){
    cat('\nconverting array to RasterStack object')
    
    if(missing(layer)) layer <- 1:dim(z)[3]
    h <- raster()
    for(ts in layer){
      h.add <- raster(t(z[,,ts])[ncol(z):1,])
      h <- raster::addLayer(h,h.add)
    }
  }else{
    cat('\nconverting matrix to RasterLayer')    
    h <- raster(t(z)[ncol(z):1,])
    
  }
  if(!missing(x) & !missing(y)){
    
    raster::extent(h) <- raster::extent(c(range(x),range(y)))
  }else{
    warning('\nconverting without geographical coordinates!') 
  }
  raster::projection(h) <- proj
  return(h)
}