a07_visual_maps.R
In ALUES: Agricultural Land Use Evaluation System

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(echo = TRUE)

## -----------------------------------------------------------------------------
library(ALUES)
y <- MarinduqueLT
banana_suit <- suit("banana", terrain=y)
banana_ovsuit <- overall_suit(banana_suit[["soil"]], method="average")

## ---- eval=FALSE--------------------------------------------------------------
#  library(ggmap)
#  library(raster)
#  library(reshape2)
#  map_lvl0 <- getData("GADM", country = "PHL", level = 0)
#  map_lvl2 <- getData("GADM", country = "PHL", level = 2)
#  
#  prov <- map_lvl2[map_lvl2$NAME_1 == as.character("Marinduque"),]
#  munic_coord <- coordinates(prov)
#  munic_coord <- data.frame(munic_coord)
#  munic_coord$label <- prov@data$NAME_2
#  
#  val <- banana_suit[["soil"]][[2]]
#  val["Overall Suitability"] <- banana_ovsuit[,1]
#  d_map <- melt(as.matrix(val))
#  d_map$Lon <- rep(y$Lon, ncol(val)); d_map$Lat <- rep(y$Lat, ncol(val))
#  
#  fill <- "#FFF7BC"; shadow <- "#9ECAE1"; ncol <- 3; size <- 3; alpha <- 1
#  text_opts <- list(alpha = 1, angle = 0, colour = "black", family = "sans", fontface = 1, lineheight = 1, size = 3)
#  labels <- list(title = "", xlab = "", ylab = "")
#  
#  p1 <- ggplot() + geom_polygon(data = prov, aes(long + 0.008, lat - 0.005, group = group), fill = shadow) +
#    geom_polygon(data = prov, aes(long, lat, group = group), colour = "grey50", fill = fill) +
#    geom_tile(aes(x = Lat, y = Lon, fill = value), data = d_map, size = size, alpha = alpha) +
#    facet_wrap(~ Var2, ncol = ncol) +
#    geom_polygon(data = prov, aes(long, lat, group = group), colour = "#4E4E4C", alpha = 0) +
#    geom_label(data = munic_coord, aes(x = X1, y = X2, label = label), alpha = 0.5,
#              angle = text_opts$angle, colour = "white", fill = "black", family = text_opts$family,
#              fontface = text_opts$fontface,
#              lineheight = text_opts$lineheight, size = text_opts$size) +
#    coord_equal() + ggtitle(as.character(labels$title)) + xlab(as.character(labels$xlab)) + ylab(as.character(labels$ylab)) +
#    scale_fill_gradientn(name = "Score\n", colors = c("red", "#FFDF00")) +
#    scale_x_continuous(breaks = round(seq(min(d_map$Lat) + 0.05, max(d_map$Lat), len = 3), 2)) +
#    theme(panel.background = element_rect(fill = '#F7E7CE'),
#          strip.background = element_rect(fill = "#D4BF96"),
#          strip.text.x = element_text(size = 12),
#          axis.text.x = element_text(size=12),
#          legend.text=element_text(size=12),
#          legend.title=element_text(size=12),
#          axis.text.y = element_text(size=12), legend.position = c(0.85, 0.25))
#  p1

## ---- out.width="100%", echo=FALSE--------------------------------------------
knitr::include_graphics("../vignettes/img/scores1.jpg")

## ---- eval=FALSE--------------------------------------------------------------
#  val <- banana_suit[["soil"]][[3]]
#  val["Overall Suitability"] <- banana_ovsuit[,2]
#  d_map <- melt(as.matrix(val))
#  d_map$Lon <- rep(y$Lon, ncol(val)); d_map$Lat <- rep(y$Lat, ncol(val))
#  
#  d_map$Class <- factor(d_map$value, levels=c("N", "S3", "S2", "S1"))
#  
#  p1 <- ggplot() + geom_polygon(data = prov, aes(long + 0.008, lat - 0.005, group = group), fill = shadow) +
#    geom_polygon(data = prov, aes(long, lat, group = group), colour = "grey50", fill = fill) +
#    geom_tile(aes(x = Lat, y = Lon, fill = Class), data = d_map, size = size, alpha = alpha) +
#    facet_wrap(~ Var2, ncol = ncol) +
#    geom_polygon(data = prov, aes(long, lat, group = group), colour = "#4E4E4C", alpha = 0) +
#    geom_label(data = munic_coord, aes(x = X1, y = X2, label = label), alpha = 0.5,
#               angle = text_opts$angle, colour = "white", fill = "black", family = text_opts$family,
#               fontface = text_opts$fontface,
#               lineheight = text_opts$lineheight, size = text_opts$size) +
#    coord_equal() + ggtitle(as.character(labels$title)) + xlab(as.character(labels$xlab)) + ylab(as.character(labels$ylab)) +
#    scale_colour_discrete(name = "Class\n", breaks=c("N", "S3", "S2", "S1"), labels=c("N", "S3", "S2", "S1")) +
#    scale_x_continuous(breaks = round(seq(min(d_map$Lat) + 0.05, max(d_map$Lat), len = 3), 2)) +
#    theme(panel.background = element_rect(fill = '#F7E7CE'),
#          strip.background = element_rect(fill = "#D4BF96"),
#          strip.text.x = element_text(size = 12),
#          axis.text.x = element_text(size=12),
#          legend.text=element_text(size=12),
#          legend.title=element_text(size=12),
#          axis.text.y = element_text(size=12), legend.position = c(0.85, 0.25))
#  p1

## ---- out.width="100%", echo=FALSE--------------------------------------------
knitr::include_graphics("../vignettes/img/classes1.jpg")