# Load and prepare data ----------

# Read in GCNP boundaries

dsn = unzip("D:/birasafab overall/Downloads/rwa_adm_2006_nisr_wgs1984_20181002_shp.zip")
data <- st_read(dsn[30])
## Reading layer `rwa_adm3_2006_NISR_WGS1984_20181002' from data source 
##   `C:\Users\fabrice.birasa\Documents\My files\Learning\R learning\Creating 3D maps\rwa_adm3_2006_NISR_WGS1984_20181002.shp' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 416 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 28.86173 ymin: -2.84023 xmax: 30.89975 ymax: -1.047167
## Geodetic CRS:  WGS 84

Creating the first rwandan map with sf package

#rwanda_map = 

  ggplot(data)+
  geom_sf(
    aes(fill =ADM1_EN)
    )+
  ggthemes::theme_map()+
  labs(title = "This is the map of Rwanda subdivided into 416 sectors \n",
       caption = "Created By Birasafab | Email: birasafab@gmail.com | data source: Rwandan Shapefiles",
       fill = "Province/CoK:")+
  theme(legend.position = "top",
        plot.title =element_text(size = 16))