veg_map <- function(fsoils,fmap,fwr) {
veg_f <- raster(fsoils)
f_soilslevs <- levels(fsoils)[[1]]$SoilCat%>%as.character()
veg_f[fsoils==which(f_soilslevs=="ALLUVIUM")] <- 1 # 1=Alluvial Savanna
veg_f[fsoils==which(f_soilslevs=="CALCSAND_DEEP")] <- 2 # 2=Dune Fynbos-Thicket Mosaic
veg_f[fsoils==which(f_soilslevs=="CALCSAND_SHALLOW")] <- 3 # 3=Limestone Fynbos
veg_f[fsoils==which(f_soilslevs=="FERTILE LOAMS AND CLAYS")] <- 4 # 4=Renosterveld
veg_f[fsoils==which(f_soilslevs=="ACID SANDS_DEEP")] <- 5 # 3=Sand Fynbos
veg_f[fsoils==which(f_soilslevs=="ACID SANDS_SHALLOW")] <- 6 # 6=Sandstone Fynbos
veg_f[veg_f==4&(fwr<0.5)&fmap>=350]<-7 # Grassland
veg_f[veg_f==4&(fwr<0.5)&fmap<350]<-11 # Grassland
veg_f[veg_f==4&(fmap>600)] <- 8 # Version.2: Shale Fynbos
# Line above changed from that below on 23-11-2018
#veg_f[veg_f==4&(fwr>=0.5)&(fmap>650)] <- 8 # Version.1: Shale Fynbos
# Line above changed from that below on 19-11-2018
#veg_f[veg_f==4&(fwr>=0.66)&(fmap>650)] <- 8 # Shale Fynbos
veg_f[!is.na(forest)]<-9 # Forest
veg_f[!is.na(thicket)]<-10 # Albany Thicket

veg_name <- c("Alluvial Savanna",#1
         "Dune Fynbos-Thicket mosaic",
         "Limestone Fynbos",
         "Renosterveld",
         "Sand Fynbos",
         "Sandstone Fynbos",
         "Grassland",
         "Shale Fynbos",
         "Forest",
         "Albany Thicket",
         "Nama-Karoo")
  veg_col <-
  c("#fdbf6f",  #"Alluvial Savanna",
    "#a6cee3",     #"Dune Fynbos-Thicket mosaic",
    "#1f78b4",    #"Limestone Fynbos",
    "#cab2d6", #"Renosterveld",
     "#fb9a99",    #"Sand Fynbos",
     "#e31a1c",    #"Sandstone Fynbos",
     "#ffffb3",    #"Grassland",
    "#6a3d9a", #"Shale Fynbos",
    "#33a02c", #"Forest",
    "#b2df8a",#"Albany Thicket"
    "brown")# Nama-Karoo )
  in.map <- values(veg_f)%>%unique(na.rm=T)%>%sort()
  levels(veg_f) <- data.frame(ID=1:length(veg_name),
                              veg_name,
                              veg_col,
                              stringsAsFactors = F)#[in.map,]
  return(veg_f)
}

veg_map_plotter <- function(fveg) {

  fcols <- levels(fveg)[[1]]$veg_col
  names(fcols)<-levels(fveg)[[1]]$ID%>%as.character()
  fname <- levels(fveg)[[1]]$veg_name
  names(fname)<-levels(fveg)[[1]]$ID%>%as.character()
  fveg%>%
  gplot(.) + 
  geom_raster(aes(fill = factor(value)),na.rm=T) +
  coord_equal(xlim=c(19.89288,23.65788),ylim=c(-36.03007,-33.70507))+
  #scale_fill_manual(values=levels(fveg)[[1]]$veg_col,
                  #   "Vegetation\nType",
                  # labels=levels(fveg)[[1]]$veg_name,drop=FALSE)+
  scale_fill_manual(values=fcols,
                    "Vegetation\nType",
                  labels=fname,drop=FALSE)+
  
    ylab("")+xlab("")+
  theme_bw()+
  theme(legend.position = "bottom")+
    #theme(legend.position = "right")+
#  theme(legend.position='none')+
  geom_polygon(data=provs,aes(x=long,y=lat,group=group),
               inherit.aes = F,  
               fill=NA,color=rgb(0,0,0,1),na.rm=T)->
    figure
  return(figure)
}