R/p1Calc.R

In BioCalc: Bioclim Variables Calculator

# Author: Julian Ramirez, dawnpatrolmustaine@gmail.com
# Date :  December 2009
# Version 0.1
# Licence GPL v3

p1Calc <- function(rlist, outfile, format='') {
	
	if (!is.list(rlist)) {
		stop('First argument should be a list or rasters (tmean)')
	}
	
	if (!file.exists(outfile)) {
			cat("", "\n", "Annual Mean Temperature (P1)", "\n")
			
			tmpstack <- stack(rlist)
			rs <- raster(tmpstack, 1)
			
			if (!canProcessInMemory(rs, 3)) {
				b1fun <- function(x) {round(mean(x))}
				tmpstack <- stack(rlist)
				p1 <- calc(tmpstack, b1fun, progress="text", filename=outfile, format=format, overwrite=TRUE)
				rm(tmpstack)
			} else {
				
				pb <- pbCreate(12, type='text', style=3)
				for (i in 1:12) {
					if (i == 1) { 
						p1 <- readAll(raster(tmpstack, i))
					} else {
						p1 <- p1 + readAll(raster(tmpstack, i))
					}
					pbStep(pb, i)
				}
				pbClose(pb)
				
				p1 <- round(p1 / 12)
				p1 <- writeRaster(p1, outfile, format=format, overwrite=TRUE)
				
			}
		} else {
			cat("", "\n", "File bio_1 already exists, skipping calculation, but loading", "\n")
			p1 <- raster(outfile)
		}
		return(p1)
}