Aggreagate with majority

majority_aggregate.R

library(raster)

# HUOM! Kaikissa hakemistopoluissa pitää Windows-kenoviiva korvata normikenoviivalla
# esim. C:\Kansio\tiedosto.txt -> C:/Kansio/tiedosto.txt

workspace <- "[KANSIO-JOSSA-LAHTORASTERIT]"
# Oletus: rasterit ovat yllä mainitussa kansiossa
org.raster.file <- "[ALKUPERAISEN-RASTERIN-NIMI]"
agg.raster.file <- "[AGGREGOIDUN-RASTERIN-NIMI]"

setwd(workspace)

custom.modal <- function(x, ...) {
  return(modal(x, ties="highest", ...))
}

org.raster <- raster(org.raster.file)
(aggregate(org.raster, fact=3, fun=custom.modal, 
           filename=agg.raster.file, option=c("COMPRESS=YES")))