bk_program_map_wo_Kosovo.R

  library(ggplot2)
  library(ggmap)
  library(rgdal)

## Loading required package: sp

## rgdal: version: 1.2-13, (SVN revision 686)
##  Geospatial Data Abstraction Library extensions to R successfully loaded
##  Loaded GDAL runtime: GDAL 2.2.0, released 2017/04/28
##  Path to GDAL shared files: C:/Users/IRCUser/Documents/R/R-3.4.2/library/rgdal/gdal
##  Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
##  Path to PROJ.4 shared files: C:/Users/IRCUser/Documents/R/R-3.4.2/library/rgdal/proj
##  Linking to sp version: 1.2-5

  library(ggrepel)
  library(cowplot)

## 
## Attaching package: 'cowplot'

## The following object is masked from 'package:ggmap':
## 
##     theme_nothing

## The following object is masked from 'package:ggplot2':
## 
##     ggsave

  library(rgdal)
  
  ### Importing shape Files 
  
  serbiamap <- readOGR("C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer="SRB_adm0")

## OGR data source with driver: ESRI Shapefile 
## Source: "C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer: "SRB_adm0"
## with 1 features
## It has 70 fields
## Integer64 fields read as strings:  ID_0 OBJECTID_1

  bulgariamap <- readOGR("C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer="BGR_adm0")

## OGR data source with driver: ESRI Shapefile 
## Source: "C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer: "BGR_adm0"
## with 1 features
## It has 70 fields
## Integer64 fields read as strings:  ID_0 OBJECTID_1

  hungarymap <- readOGR("C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer="HUN_adm0")

## OGR data source with driver: ESRI Shapefile 
## Source: "C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer: "HUN_adm0"
## with 1 features
## It has 70 fields
## Integer64 fields read as strings:  ID_0 OBJECTID_1

  kosovomap <- readOGR("C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer="KO__adm0")

## OGR data source with driver: ESRI Shapefile 
## Source: "C:/Users/IRCUser/Box Sync/R Analysis/Shapefile/Serbia_shapefile", layer: "KO__adm0"
## with 1 features
## It has 70 fields

  plot(serbiamap, 
       col = "#FDC82F", 
       border="#000000")
  
  plot(bulgariamap, 
       add  = TRUE, 
       border="000000",
       col = "#FDC82F")
  
  plot(hungarymap, 
       add  = TRUE, 
       border="000000",
       col = "#FDC82F")
  
  plot(kosovomap, 
       add  = TRUE, 
       border="000000",
       col = "#FDC82F")

  ### Importing Data Set
  
  tr <- read.csv("C:/Users/IRCUser/Box Sync/R Analysis/Protection Monitoring Mapping/bkcentersdata.csv", comment.char="#")
  df <- NULL
  df$NAME <- tr$NAME
  df$lat<- tr$lat
  df$lon<- tr$lon
  df$PARTNER<- tr$PARTNER

  
   df<-as.data.frame(df) 
  area.points <- fortify(serbiamap)

## Regions defined for each Polygons

  area.points2 <- fortify(bulgariamap)

## Regions defined for each Polygons

  area.points3 <- fortify(hungarymap)

## Regions defined for each Polygons

 ### Plotting Map 
  ircmap <- get_map(location = c(lon = mean(df$lon), lat = mean(df$lat)), zoom = 5,
                        source = "google" , maptype = "hybrid", color="bw", scale = 1)

## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=44.201329,21.195751&zoom=5&size=640x640&scale=1&maptype=hybrid&language=en-EN&sensor=false

  ircmap <- ggmap(ircmap, darken = c(0.7, "black")) +
    geom_polygon(aes(x = long,
                     y = lat,
                     group = group),
                 data = area.points,
                 alpha = 0.99,size = 0.95,color="#cc6600",fill="#FDC82F")+
    geom_polygon(aes(x = long,
                     y = lat,
                     group = group),
                 data = area.points2,
                 alpha = 0.99,size = 0.95,color="#cc6600",fill="#FDC82F")+
    geom_polygon(aes(x = long,
                     y = lat,
                     group = group),
                 data = area.points3,
                 alpha = 0.99,size = 0.95,color="#cc6600",fill="#FDC82F")+

    
    geom_point(data = df, aes(x = lon, y = lat, color= PARTNER), stroke = 2, fill="white",shape = 21, alpha=0.9)+
    
    geom_label_repel (data = df, aes(x =lon, y = lat,label = NAME), fontface = 'bold', color = 'black',
                      box.padding = 0.35, point.padding = 0.5,force = 3, size = 2)

 
   plot(ircmap)