axmedmaxamuud/geosom

Somalia Geospatial Database # geosom geosom is providing ready-to-use geographic boundary shapes and other country datasets in Somalia for easy spatial analysis and mapping workflows in R. ## Installation You can install the development version from Github: ``` r # install.packages("devtools") remotes::install_github("axmedmaxamuud/geosom") ``` ## Available Layers ``` r library(geosom) library(ggplot2) library(dplyr) list_geosom_layers() #> layer description #> 1 admin0 Somalia national boundary #> 2 admin1 Administrative level 1 boundaries #> 3 admin2 Administrative level 2 boundaries #> 4 states Somalia state boundaries ``` ## Load administrative boundaries ``` r admin0 admin1name admin1pcod admin0name #> 1 Awdal SO11 Somalia #> 2 Bakool SO25 Somalia #> 3 Banadir SO22 Somalia #> 4 Bari SO16 Somalia #> 5 Bay SO24 Somalia #> 6 Galgaduud SO19 Somalia #> 7 Gedo SO26 Somalia #> 8 Hiraan SO20 Somalia #> 9 Middle Juba SO27 Somalia #> 10 Lower Juba SO28 Somalia ``` ``` r head(get_admin_names("admin2")[, 1:3], 10) #> admin2name admin2pcod admin1name #> 1 Banadir SO2201 Banadir #> 2 Borama SO1101 Awdal #> 3 Baki SO1102 Awdal #> 4 Lughaye SO1103 Awdal #> 5 Zeylac SO1104 Awdal #> 6 Hargeysa SO1201 Woqooyi Galbeed #> 7 Berbera SO1202 Woqooyi Galbeed #> 8 Gebiley SO1203 Woqooyi Galbeed #> 9 Burco SO1301 Togdheer #> 10 Buuhoodle SO1302 Togdheer ``` ## Basic map ``` r admin1 ## Population data ``` r admin2 % left_join( population_admin2_2026, by = c("admin2pcod" = "admin2_pcode") ) head(admin2_pop[, 2:10], 10) #> Simple feature collection with 10 features and 9 fields #> Geometry type: MULTIPOLYGON #> Dimension: XY #> Bounding box: xmin: 42.6 …