R Markdown Tutorial for Maps
Highlighting portions of a Country Map
library(ggplot2)
library(rnaturalearth)
library(rnaturalearthdata)
##
## Attaching package: 'rnaturalearthdata'
## The following object is masked from 'package:rnaturalearth':
##
## countries110
# Get Philippines map
ph_map <- ne_countries(country = "Philippines", scale = "medium", returnclass = "sf")
# Plot map
ggplot(data = ph_map) +
geom_sf(fill = "lightblue", color = "black") +
theme_minimal() +
labs(title = "Map of the Philippines")
# Load libraries
library(ggplot2)
library(sf)
## Linking to GEOS 3.13.0, GDAL 3.10.1, PROJ 9.5.1; sf_use_s2() is TRUE
library(rnaturalearth)
library(rnaturalearthdata)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
# Load Philippines country boundary
ph_map <- ne_countries(country = "Philippines", scale = "medium", returnclass = "sf")
# Load Philippines provincial shapefile (GADM Level 1)
ph_prov <- st_read("C:/Users/User/Desktop/Spatial Analysis/gadm41_PHL_1.json/gadm41_PHL_1.json", layer = "gadm41_PHL_1")
## Reading layer `gadm41_PHL_1' from data source
## `C:\Users\User\Desktop\Spatial Analysis\gadm41_PHL_1.json\gadm41_PHL_1.json'
## using driver `GeoJSON'
## Simple feature collection with 81 features and 11 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 116.9283 ymin: 4.5869 xmax: 126.6052 ymax: 21.0699
## Geodetic CRS: WGS 84
# Check available province names
unique(ph_prov$NAME_1) # Verify province names in GADM
## [1] "Abra" "AgusandelNorte" "AgusandelSur"
## [4] "Aklan" "Albay" "Antique"
## [7] "Apayao" "Aurora" "Basilan"
## [10] "Bataan" "Batanes" "Batangas"
## [13] "Benguet" "Biliran" "Bohol"
## [16] "Bukidnon" "Bulacan" "Cagayan"
## [19] "CamarinesNorte" "CamarinesSur" "Camiguin"
## [22] "Capiz" "Catanduanes" "Cavite"
## [25] "Cebu" "CompostelaValley" "DavaodelNorte"
## [28] "DavaodelSur" "DavaoOriental" "DinagatIslands"
## [31] "EasternSamar" "Guimaras" "Ifugao"
## [34] "IlocosNorte" "IlocosSur" "Iloilo"
## [37] "Isabela" "Kalinga" "LaUnion"
## [40] "Laguna" "LanaodelNorte" "LanaodelSur"
## [43] "Leyte" "Maguindanao" "Marinduque"
## [46] "Masbate" "MetropolitanManila" "MisamisOccidental"
## [49] "MisamisOriental" "MountainProvince" "NegrosOccidental"
## [52] "NegrosOriental" "NorthCotabato" "NorthernSamar"
## [55] "NuevaEcija" "NuevaVizcaya" "OccidentalMindoro"
## [58] "OrientalMindoro" "Palawan" "Pampanga"
## [61] "Pangasinan" "Quezon" "Quirino"
## [64] "Rizal" "Romblon" "Samar"
## [67] "Sarangani" "Siquijor" "Sorsogon"
## [70] "SouthCotabato" "SouthernLeyte" "SultanKudarat"
## [73] "Sulu" "SurigaodelNorte" "SurigaodelSur"
## [76] "Tarlac" "Tawi-Tawi" "Zambales"
## [79] "ZamboangadelNorte" "ZamboangadelSur" "ZamboangaSibugay"
# Correct province names for Davao Region (Region XI)
davao_region <- ph_prov %>% filter(NAME_1 %in% c("DavaodelNorte", "DavaodelSur", "DavaoOriental","CompostelaValley", "DavaoOccidental"))
# Plot the Philippines with Davao Region highlighted
ggplot() +
geom_sf(data = ph_prov, fill = "lightgray", color = "black") + # Other provinces in gray
geom_sf(data = davao_region, fill = "red", color = "black") + # Davao Region in red
theme_minimal() +
labs(title = "Highlighting Davao Region in the Philippines")
Baranggay Level Maps within the City
# Load required libraries
library(sf)
library(ggplot2)
library(dplyr)
# Load the GADM Level 3 shapefile (Update the file path)
ph_l3<- st_read("C:/Users/User/Desktop/Spatial Analysis/gadm41_PHL_3.json/gadm41_PHL_3.json", layer = "gadm41_PHL_3")
## Reading layer `gadm41_PHL_3' from data source
## `C:\Users\User\Desktop\Spatial Analysis\gadm41_PHL_3.json\gadm41_PHL_3.json'
## using driver `GeoJSON'
## Simple feature collection with 41948 features and 16 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 116.9283 ymin: 4.5869 xmax: 126.6053 ymax: 21.0699
## Geodetic CRS: WGS 84
# Check available columns
colnames(ph_l3)
## [1] "GID_3" "GID_0" "COUNTRY" "GID_1" "NAME_1" "NL_NAME_1"
## [7] "GID_2" "NAME_2" "NL_NAME_2" "NAME_3" "VARNAME_3" "NL_NAME_3"
## [13] "TYPE_3" "ENGTYPE_3" "CC_3" "HASC_3" "geometry"
# Check unique city/municipality names
unique(ph_l3$NAME_2)
## [1] "Bangued" "Boliney"
## [3] "Bucay" "Bucloc"
## [5] "Daguioman" "Danglas"
## [7] "Dolores" "LaPaz"
## [9] "Lacub" "Lagangilang"
## [11] "Lagayan" "Langiden"
## [13] "Licuan-Baay" "Luba"
## [15] "Malibcong" "Manabo"
## [17] "Peñarrubia" "Pidigan"
## [19] "Pilar" "Sallapadan"
## [21] "SanIsidro" "SanJuan"
## [23] "SanQuintin" "Tayum"
## [25] "Tineg" "Tubo"
## [27] "Villaviciosa" "Buenavista"
## [29] "ButuanCity" "CabadbaranCity"
## [31] "Carmen" "Jabonga"
## [33] "Kitcharao" "LasNieves"
## [35] "Magallanes" "MainitLake"
## [37] "Nasipit" "RemediosT.Romualdez"
## [39] "Santiago" "Tubay"
## [41] "BayuganCity" "Bunawan"
## [43] "Esperanza" "Loreto"
## [45] "Prosperidad" "Rosario"
## [47] "SanFrancisco" "SanLuis"
## [49] "SantaJosefa" "Sibagat"
## [51] "Talacogon" "Trento"
## [53] "Veruela" "Altavas"
## [55] "Balete" "Banga"
## [57] "Batan" "Buruanga"
## [59] "Ibajay" "Kalibo"
## [61] "Lezo" "Libacao"
## [63] "Madalag" "Makato"
## [65] "Malay" "Malinao"
## [67] "Nabas" "NewWashington"
## [69] "Numancia" "Tangalan"
## [71] "Bacacay" "BatoLake"
## [73] "Camalig" "Daraga"
## [75] "Guinobatan" "Jovellar"
## [77] "LegazpiCity" "Libon"
## [79] "LigaoCity" "Malilipot"
## [81] "Manito" "Oas"
## [83] "PioDuran" "Polangui"
## [85] "Rapu-Rapu" "SantoDomingo"
## [87] "TabacoCity" "Tiwi"
## [89] "Anini-Y" "Barbaza"
## [91] "Belison" "Bugasong"
## [93] "Caluya" "Culasi"
## [95] "Hamtic" "Laua-An"
## [97] "Libertad" "Pandan"
## [99] "Patnongon" "SanJose"
## [101] "SanRemigio" "Sebaste"
## [103] "Sibalom" "Tibiao"
## [105] "TobiasFornier" "Valderrama"
## [107] "Calanasan" "Conner"
## [109] "Flora" "Kabugao"
## [111] "Luna" "Pudtol"
## [113] "SantaMarcela" "Baler"
## [115] "Casiguran" "Dilasag"
## [117] "Dinalungan" "Dingalan"
## [119] "Dipaculao" "MariaAurora"
## [121] "Akbar" "Al-Barka"
## [123] "HadjiMohammadAjul" "IsabelaCity"
## [125] "LamitanCity" "Lantawan"
## [127] "Maluso" "Sumisip"
## [129] "Tipo-Tipo" "Tuburan"
## [131] "UngkayaPukan" "Abucay"
## [133] "Bagac" "BalangaCity"
## [135] "Dinalupihan" "Hermosa"
## [137] "Limay" "Mariveles"
## [139] "Morong" "Orani"
## [141] "Orion" "Samal"
## [143] "Basco" "Itbayat"
## [145] "Ivana" "Mahatao"
## [147] "Sabtang" "Uyugan"
## [149] "Agoncillo" "Alitagtag"
## [151] "Balayan" "BatangasCity"
## [153] "Bauan" "Calaca"
## [155] "Calatagan" "Cuenca"
## [157] "Ibaan" "Laurel"
## [159] "Lemery" "Lian"
## [161] "LipaCity" "Lobo"
## [163] "Mabini" "Malvar"
## [165] "MataasNaKahoy" "Nasugbu"
## [167] "PadreGarcia" "SanNicolas"
## [169] "SanPascual" "SantaTeresita"
## [171] "SantoTomas" "Taal"
## [173] "Taallake" "Talisay"
## [175] "TanauanCity" "Taysan"
## [177] "Tingloy" "Tuy"
## [179] "Atok" "BaguioCity"
## [181] "Bakun" "Bokod"
## [183] "Buguias" "Itogon"
## [185] "Kabayan" "Kapangan"
## [187] "Kibungan" "LaTrinidad"
## [189] "Mankayan" "Sablan"
## [191] "Tuba" "Tublay"
## [193] "Almeria" "Biliran"
## [195] "Cabucgayan" "Caibiran"
## [197] "Culaba" "Kawayan"
## [199] "Maripipi" "Naval"
## [201] "Albuquerque" "Alicia"
## [203] "Anda" "Antequera"
## [205] "Baclayon" "Balilihan"
## [207] "Batuan" "BienUnido"
## [209] "Bilar" "Calape"
## [211] "Candijay" "Catigbian"
## [213] "Clarin" "Corella"
## [215] "Cortes" "Dagohoy"
## [217] "Danao" "Dauis"
## [219] "Dimiao" "Duero"
## [221] "GarciaHernandez" "Guindulman"
## [223] "Inabanga" "Jagna"
## [225] "Jetafe" "Lila"
## [227] "Loay" "Loboc"
## [229] "Loon" "Maribojoc"
## [231] "Panglao" "Pres.CarlosP.Garcia"
## [233] "Sagbayan" "SanMiguel"
## [235] "Sevilla" "SierraBullones"
## [237] "Sikatuna" "TagbilaranCity"
## [239] "Talibon" "Trinidad"
## [241] "Tubigon" "Ubay"
## [243] "Valencia" "Baungon"
## [245] "Cabanglasan" "Damulog"
## [247] "Dangcagan" "DonCarlos"
## [249] "Impasug-Ong" "Kadingilan"
## [251] "Kalilangan" "Kibawe"
## [253] "Kitaotao" "Lantapan"
## [255] "Libona" "MalaybalayCity"
## [257] "Malitbog" "ManoloFortich"
## [259] "Maramag" "Pangantucan"
## [261] "Quezon" "SanFernando"
## [263] "Sumilao" "Talakag"
## [265] "ValenciaCity" "Angat"
## [267] "Balagtas" "Baliuag"
## [269] "Bocaue" "Bulacan"
## [271] "Bustos" "Calumpit"
## [273] "DoñaRemediosTrinidad" "Guiguinto"
## [275] "Hagonoy" "MalolosCity"
## [277] "Marilao" "MeycauayanCity"
## [279] "Norzagaray" "Obando"
## [281] "Pandi" "Paombong"
## [283] "Plaridel" "Pulilan"
## [285] "SanIldefonso" "SanJosedelMonteCity"
## [287] "SanRafael" "SantaMaria"
## [289] "Abulug" "Alcala"
## [291] "Allacapan" "Amulung"
## [293] "Aparri" "Baggao"
## [295] "Ballesteros" "Buguey"
## [297] "Calayan" "Camalaniugan"
## [299] "Claveria" "Enrile"
## [301] "Gattaran" "Gonzaga"
## [303] "Iguig" "Lal-Lo"
## [305] "Lasam" "Pamplona"
## [307] "Peñablanca" "Piat"
## [309] "Rizal" "Sanchez-Mira"
## [311] "SantaAna" "SantaPraxedes"
## [313] "SantoNiño" "Solana"
## [315] "Tuao" "TuguegaraoCity"
## [317] "Basud" "Capalonga"
## [319] "Daet" "JosePanganiban"
## [321] "Labo" "Mercedes"
## [323] "Paracale" "SanLorenzoRuiz"
## [325] "SanVicente" "SantaElena"
## [327] "Vinzons" "Baao"
## [329] "Balatan" "Bato"
## [331] "Bombon" "Buhi"
## [333] "BuhiLake" "Bula"
## [335] "Cabusao" "Calabanga"
## [337] "Camaligan" "Canaman"
## [339] "Caramoan" "DelGallego"
## [341] "Gainza" "Garchitorena"
## [343] "Goa" "IrigaCity"
## [345] "Lagonoy" "Libmanan"
## [347] "Lupi" "Magarao"
## [349] "Milaor" "Minalabac"
## [351] "Nabua" "NagaCity"
## [353] "Ocampo" "Pasacao"
## [355] "Pili" "Presentacion"
## [357] "Ragay" "Sagnay"
## [359] "Sipocot" "Siruma"
## [361] "Tigaon" "Tinambac"
## [363] "Catarman" "Guinsiliban"
## [365] "Mahinog" "Mambajao"
## [367] "Sagay" "Cuartero"
## [369] "Dao" "Dumalag"
## [371] "Dumarao" "Ivisan"
## [373] "Jamindan" "Ma-Ayon"
## [375] "Mambusao" "Panay"
## [377] "Panitan" "Pontevedra"
## [379] "PresidentRoxas" "RoxasCity"
## [381] "Sapi-An" "Sigma"
## [383] "Tapaz" "Bagamanoc"
## [385] "Baras" "Caramoran"
## [387] "Gigmoto" "Panganiban"
## [389] "SanAndres" "Viga"
## [391] "Virac" "Alfonso"
## [393] "Amadeo" "Bacoor"
## [395] "Carmona" "CaviteCity"
## [397] "Dasmariñas" "GeneralEmilioAguinaldo"
## [399] "GeneralMarianoAlvarez" "GeneralTrias"
## [401] "Imus" "Indang"
## [403] "Kawit" "Maragondon"
## [405] "Mendez" "Naic"
## [407] "Noveleta" "Silang"
## [409] "TagaytayCity" "Tanza"
## [411] "Ternate" "TreceMartiresCity"
## [413] "Alcantara" "Alcoy"
## [415] "Alegria" "Aloguinsan"
## [417] "Argao" "Asturias"
## [419] "Badian" "Balamban"
## [421] "Bantayan" "Barili"
## [423] "BogoCity" "Boljoon"
## [425] "Borbon" "Carcar"
## [427] "Catmon" "CebuCity"
## [429] "Compostela" "Consolacion"
## [431] "Cordoba" "Daanbantayan"
## [433] "Dalaguete" "DanaoCity"
## [435] "DanaoLake" "Dumanjug"
## [437] "Ginatilan" "Lapu-LapuCity"
## [439] "Liloan" "Madridejos"
## [441] "Malabuyoc" "MandaueCity"
## [443] "Medellin" "Minglanilla"
## [445] "Moalboal" "Oslob"
## [447] "Pinamungahan" "Poro"
## [449] "Ronda" "Samboan"
## [451] "SantaFe" "Santander"
## [453] "Sibonga" "Sogod"
## [455] "Tabogon" "Tabuelan"
## [457] "TalisayCity" "ToledoCity"
## [459] "Tudela" "Laak"
## [461] "Maco" "Maragusan"
## [463] "Mawab" "Monkayo"
## [465] "Montevista" "Nabunturan"
## [467] "NewBataan" "Pantukan"
## [469] "Asuncion" "BraulioE.Dujali"
## [471] "Kapalong" "NewCorella"
## [473] "PanaboCity" "SamalCity"
## [475] "TagumCity" "Talaingod"
## [477] "Bansalan" "DavaoCity"
## [479] "DigosCity" "DonMarcelino"
## [481] "JoseAbadSantos" "Kiblawan"
## [483] "Magsaysay" "Malalag"
## [485] "Malita" "Matanao"
## [487] "Padada" "SantaCruz"
## [489] "Sarangani" "Sulop"
## [491] "Baganga" "Banaybanay"
## [493] "Boston" "Caraga"
## [495] "Cateel" "GovernorGeneroso"
## [497] "Lupon" "Manay"
## [499] "MatiCity" "Tarragona"
## [501] "Basilisa" "Cagdianao"
## [503] "Dinagat" "Libjo"
## [505] "Tubajon" "Arteche"
## [507] "Balangiga" "Balangkayan"
## [509] "BoronganCity" "Can-Avid"
## [511] "GeneralMacarthur" "Giporlos"
## [513] "Guiuan" "Hernani"
## [515] "Jipapad" "Lawaan"
## [517] "Llorente" "Maslog"
## [519] "Maydolong" "Oras"
## [521] "Quinapondan" "Salcedo"
## [523] "SanJulian" "SanPolicarpo"
## [525] "Sulat" "Taft"
## [527] "Jordan" "NuevaValencia"
## [529] "SanLorenzo" "Sibunag"
## [531] "Aguinaldo" "AlfonsoLista"
## [533] "Asipulo" "Banaue"
## [535] "Hingyon" "Hungduan"
## [537] "Kiangan" "Lagawe"
## [539] "Lamut" "Mayoyao"
## [541] "Tinoc" "Adams"
## [543] "Bacarra" "Badoc"
## [545] "Bangui" "Banna"
## [547] "BatacCity" "Burgos"
## [549] "Carasi" "Currimao"
## [551] "Dingras" "Dumalneg"
## [553] "LaoagCity" "Marcos"
## [555] "NuevaEra" "Pagudpud"
## [557] "Paoay" "PaoayLake"
## [559] "Pasuquin" "Piddig"
## [561] "Pinili" "Sarrat"
## [563] "Solsona" "Vintar"
## [565] "Alilem" "Banayoyo"
## [567] "Bantay" "Cabugao"
## [569] "CandonCity" "Caoayan"
## [571] "Cervantes" "Galimuyod"
## [573] "GregorioDelPilar" "Lidlidda"
## [575] "Magsingal" "Nagbukel"
## [577] "Narvacan" "Quirino"
## [579] "SanEmilio" "SanEsteban"
## [581] "Santa" "SantaCatalina"
## [583] "SantaLucia" "Sigay"
## [585] "Sinait" "Sugpon"
## [587] "Suyo" "Tagudin"
## [589] "ViganCity" "Ajuy"
## [591] "Alimodian" "Anilao"
## [593] "Badiangan" "Balasan"
## [595] "Banate" "BarotacNuevo"
## [597] "BarotacViejo" "Batad"
## [599] "Bingawan" "Cabatuan"
## [601] "Calinog" "Carles"
## [603] "Concepcion" "Dingle"
## [605] "Duenas" "Dumangas"
## [607] "Estancia" "Guimbal"
## [609] "Igbaras" "IloiloCity"
## [611] "Janiuay" "Lambunao"
## [613] "Leganes" "Leon"
## [615] "Maasin" "Miagao"
## [617] "Mina" "NewLucena"
## [619] "Oton" "PassiCity"
## [621] "Pavia" "Pototan"
## [623] "SanDionisio" "SanEnrique"
## [625] "SanJoaquin" "SantaBarbara"
## [627] "Sara" "Tigbauan"
## [629] "Tubungan" "Zarraga"
## [631] "Angadanan" "Aurora"
## [633] "BenitoSoliven" "Cabagan"
## [635] "CauayanCity" "Cordon"
## [637] "DelfinAlbano" "Dinapigue"
## [639] "Divilacan" "Echague"
## [641] "Gamu" "Ilagan"
## [643] "Jones" "Maconacon"
## [645] "Mallig" "Naguilian"
## [647] "Palanan" "Ramon"
## [649] "ReinaMercedes" "Roxas"
## [651] "SanAgustin" "SanGuillermo"
## [653] "SanManuel" "SanMariano"
## [655] "SanMateo" "SanPablo"
## [657] "SantiagoCity" "Tumauini"
## [659] "Balbalan" "Lubuagan"
## [661] "Pasil" "Pinukpuk"
## [663] "TabukCity" "Tanudan"
## [665] "Tinglayan" "Agoo"
## [667] "Aringay" "Bacnotan"
## [669] "Bagulin" "Balaoan"
## [671] "Bangar" "Bauang"
## [673] "Caba" "Pugo"
## [675] "SanFernandoCity" "SanGabriel"
## [677] "Santol" "Sudipen"
## [679] "Tubao" "Alaminos"
## [681] "Bay" "Biñan"
## [683] "Cabuyao" "CalambaCity"
## [685] "Calauan" "Cavinti"
## [687] "Famy" "Kalayaan"
## [689] "KalibatoLake" "Lagunalake"
## [691] "Liliw" "LosBaños"
## [693] "Luisiana" "Lumban"
## [695] "Mabitac" "Magdalena"
## [697] "Majayjay" "Nagcarlan"
## [699] "Paete" "Pagsanjan"
## [701] "Pakil" "PalakpakinLake"
## [703] "Pangil" "Pila"
## [705] "SampalocLake" "SanPabloCity"
## [707] "SanPedro" "SantaRosaCity"
## [709] "Siniloan" "Victoria"
## [711] "Waterbody" "Bacolod"
## [713] "Baloi" "Baroy"
## [715] "IliganCity" "Kapatagan"
## [717] "Kauswagan" "Kolambugan"
## [719] "Lala" "Linamon"
## [721] "Maigo" "Matungao"
## [723] "Munai" "Nunungan"
## [725] "PantaoRagat" "Pantar"
## [727] "PoonaPiagapo" "Salvador"
## [729] "Sapad" "SultanNagaDimaporo"
## [731] "Tagoloan" "Tangcal"
## [733] "Tubod" "BacolodKalawi"
## [735] "Balabagan" "Balindong"
## [737] "Bayang" "Binidayan"
## [739] "Buadiposo-Buntong" "Bubong"
## [741] "Bumbaran" "Butig"
## [743] "Calanogas" "DapaoLake"
## [745] "Ditsaan-Ramain" "Ganassi"
## [747] "Kapai" "LanaoLake"
## [749] "Lumba-Bayabao" "LumbacaUnayan"
## [751] "Lumbatan" "Lumbayanague"
## [753] "Madalum" "Madamba"
## [755] "Maguing" "Malabang"
## [757] "Marantao" "MarawiCity"
## [759] "Marogong" "Masiu"
## [761] "Mulondo" "Pagayawan"
## [763] "Piagapo" "Picong"
## [765] "PoonaBayabao" "Pualas"
## [767] "Saguiaran" "SultanDumalondong"
## [769] "TagoloanII" "Tamparan"
## [771] "Taraka" "Tubaran"
## [773] "Tugaya" "Wao"
## [775] "Abuyog" "Alangalang"
## [777] "Albuera" "Babatngon"
## [779] "Barugo" "BaybayCity"
## [781] "Burauen" "Calubian"
## [783] "Capoocan" "Carigara"
## [785] "Dagami" "Dulag"
## [787] "Hilongos" "Hindang"
## [789] "Inopacan" "Isabel"
## [791] "Jaro" "Javier"
## [793] "Julita" "Kananga"
## [795] "Leyte" "Macarthur"
## [797] "Mahaplag" "Matag-Ob"
## [799] "Matalom" "Mayorga"
## [801] "Merida" "OrmocCity"
## [803] "Palo" "Palompon"
## [805] "Pastrana" "Tabango"
## [807] "Tabontabon" "TaclobanCity"
## [809] "Tanauan" "Tolosa"
## [811] "Tunga" "Villaba"
## [813] "Ampatuan" "Barira"
## [815] "Buldon" "Buluan"
## [817] "BuluanLake" "CotabatoCity"
## [819] "DatuAbdullahSanki" "DatuAnggalMidtimbang"
## [821] "DatuBlahT.Sinsuat" "DatuOdinSinsuat"
## [823] "DatuPaglas" "DatuPiang"
## [825] "DatuSaudi-Ampatuan" "DatuUnsay"
## [827] "Gen.S.K.Pendatun" "Guindulungan"
## [829] "Kabuntalan" "Mamasapano"
## [831] "Mangudadatu" "Matanog"
## [833] "NorthernKabuntalan" "Pagagawan"
## [835] "Pagalungan" "Paglat"
## [837] "Pandag" "Parang"
## [839] "RajahBuayan" "ShariffAguak"
## [841] "SouthUpi" "SultanKudarat"
## [843] "SultanMastura" "SultanSaBarongis"
## [845] "Talayan" "Talitay"
## [847] "Upi" "Boac"
## [849] "Gasan" "Mogpog"
## [851] "Torrijos" "Aroroy"
## [853] "Baleno" "Balud"
## [855] "Cataingan" "Cawayan"
## [857] "Dimasalang" "Mandaon"
## [859] "MasbateCity" "Milagros"
## [861] "Mobo" "Monreal"
## [863] "Palanas" "PioV.Corpuz"
## [865] "Placer" "SanJacinto"
## [867] "Uson" "KalookanCity"
## [869] "LasPiñas" "MakatiCity"
## [871] "Malabon" "Mandaluyong"
## [873] "Manila" "Marikina"
## [875] "Muntinlupa" "Navotas"
## [877] "Parañaque" "PasayCity"
## [879] "PasigCity" "Pateros"
## [881] "QuezonCity" "Taguig"
## [883] "Valenzuela" "Aloran"
## [885] "Baliangao" "Bonifacio"
## [887] "Calamba" "DonVictorianoChiongbian"
## [889] "Jimenez" "LopezJaena"
## [891] "OroquietaCity" "OzamisCity"
## [893] "Panaon" "SapangDalaga"
## [895] "Sinacaban" "TangubCity"
## [897] "Alubijid" "Balingasag"
## [899] "Balingoan" "Binuangan"
## [901] "CagayandeOroCity" "ElSalvadorCity"
## [903] "GingoogCity" "Gitagum"
## [905] "Initao" "Jasaan"
## [907] "Kinoguitan" "Lagonglong"
## [909] "Laguindingan" "Lugait"
## [911] "Manticao" "Medina"
## [913] "Naawan" "Opol"
## [915] "Salay" "Sugbongcogon"
## [917] "Talisayan" "Villanueva"
## [919] "Barlig" "Bauko"
## [921] "Besao" "Bontoc"
## [923] "Natonin" "Paracelis"
## [925] "Sabangan" "Sadanga"
## [927] "Sagada" "Tadian"
## [929] "BacolodCity" "BagoCity"
## [931] "Binalbagan" "CadizCity"
## [933] "Calatrava" "Candoni"
## [935] "Cauayan" "EnriqueB.Magalona"
## [937] "EscalanteCity" "HimamaylanCity"
## [939] "Hinigaran" "Hinoba-An"
## [941] "Ilog" "Isabela"
## [943] "KabankalanCity" "LaCarlotaCity"
## [945] "LaCastellana" "Manapla"
## [947] "MoisesPadilla" "Murcia"
## [949] "Pulupandan" "SagayCity"
## [951] "SalvadorBenedicto" "SanCarlosCity"
## [953] "SilayCity" "SipalayCity"
## [955] "Toboso" "Valladolid"
## [957] "VictoriasCity" "Amlan"
## [959] "Ayungon" "Bacong"
## [961] "BaisCity" "Basay"
## [963] "BayawanCity" "Bindoy"
## [965] "CanlaonCity" "Dauin"
## [967] "DumagueteCity" "GuihulnganCity"
## [969] "Jimalalud" "LaLibertad"
## [971] "Mabinay" "Manjuyod"
## [973] "Siaton" "Sibulan"
## [975] "TanjayCity" "Tayasan"
## [977] "Vallehermoso" "Zamboanguita"
## [979] "Alamada" "Aleosan"
## [981] "Antipas" "Arakan"
## [983] "Banisilan" "Kabacan"
## [985] "KidapawanCity" "Libungan"
## [987] "M'Lang" "Magpet"
## [989] "Makilala" "Matalam"
## [991] "Midsayap" "Pigkawayan"
## [993] "Pikit" "Tulunan"
## [995] "Allen" "Biri"
## [997] "Bobon" "Capul"
## [999] "Catubig" "Gamay"
## [1001] "Laoang" "Lapinig"
## [1003] "LasNavas" "Lavezares"
## [1005] "LopedeVega" "Mapanas"
## [1007] "Mondragon" "Palapag"
## [1009] "Pambujan" "SanAntonio"
## [1011] "SanRoque" "SilvinoLobos"
## [1013] "Aliaga" "Bongabon"
## [1015] "CabanatuanCity" "Cabiao"
## [1017] "Carranglan" "Cuyapo"
## [1019] "Gabaldon" "GapanCity"
## [1021] "GeneralMamertoNatividad" "GeneralTinio"
## [1023] "Guimba" "Jaen"
## [1025] "Laur" "Licab"
## [1027] "Llanera" "Lupao"
## [1029] "MuñozCity" "Nampicuan"
## [1031] "PalayanCity" "Pantabangan"
## [1033] "Peñaranda" "SanJoseCity"
## [1035] "SanLeonardo" "SantaRosa"
## [1037] "Talavera" "Talugtug"
## [1039] "Zaragoza" "AlfonsoCastaneda"
## [1041] "Ambaguio" "Aritao"
## [1043] "Bagabag" "Bambang"
## [1045] "Bayombong" "Diadi"
## [1047] "DupaxDelNorte" "DupaxDelSur"
## [1049] "Kasibu" "Kayapa"
## [1051] "Solano" "Villaverde"
## [1053] "AbradeIlog" "Calintaan"
## [1055] "Looc" "Lubang"
## [1057] "Mamburao" "Paluan"
## [1059] "Sablayan" "Baco"
## [1061] "Bansud" "Bongabong"
## [1063] "Bulalacao" "CalapanCity"
## [1065] "Gloria" "Mansalay"
## [1067] "Naujan" "NaujanLake"
## [1069] "Pinamalayan" "Pola"
## [1071] "PuertoGalera" "SanTeodoro"
## [1073] "Socorro" "Aborlan"
## [1075] "Agutaya" "Araceli"
## [1077] "Balabac" "Bataraza"
## [1079] "Brooke'sPoint" "Busuanga"
## [1081] "Cagayancillo" "Coron"
## [1083] "Culion" "Cuyo"
## [1085] "Dumaran" "ElNido"
## [1087] "Linapacan" "Narra"
## [1089] "PuertoPrincesaCity" "SofronioEspanola"
## [1091] "Taytay" "AngelesCity"
## [1093] "Apalit" "Arayat"
## [1095] "Bacolor" "Candaba"
## [1097] "Floridablanca" "Guagua"
## [1099] "Lubao" "Mabalacat"
## [1101] "Macabebe" "Magalang"
## [1103] "Masantol" "Mexico"
## [1105] "Minalin" "Porac"
## [1107] "SanSimon" "SantaRita"
## [1109] "Sasmuan" "Agno"
## [1111] "Aguilar" "AlaminosCity"
## [1113] "Asingan" "Balungao"
## [1115] "Bani" "Basista"
## [1117] "Bautista" "Bayambang"
## [1119] "Binalonan" "Binmaley"
## [1121] "Bolinao" "Bugallon"
## [1123] "Calasiao" "DagupanCity"
## [1125] "Dasol" "Infanta"
## [1127] "Labrador" "Laoac"
## [1129] "Lingayen" "Malasiqui"
## [1131] "Manaoag" "Mangaldan"
## [1133] "Mangatarem" "Mapandan"
## [1135] "Natividad" "Pozzorubio"
## [1137] "Rosales" "SanFabian"
## [1139] "Sison" "Sual"
## [1141] "Tayug" "Umingan"
## [1143] "Urbiztondo" "UrdanetaCity"
## [1145] "Villasis" "Agdangan"
## [1147] "Alabat" "Atimonan"
## [1149] "Burdeos" "Calauag"
## [1151] "Candelaria" "Catanauan"
## [1153] "GeneralLuna" "GeneralNakar"
## [1155] "Guinayangan" "Gumaca"
## [1157] "Hinunangan" "Jomalig"
## [1159] "Lopez" "Lucban"
## [1161] "LucenaCity" "Macalelon"
## [1163] "Mauban" "Mulanay"
## [1165] "PadreBurgos" "Pagbilao"
## [1167] "Panukulan" "Patnanungan"
## [1169] "Perez" "Pitogo"
## [1171] "Polillo" "Real"
## [1173] "Sampaloc" "SanNarciso"
## [1175] "Sariaya" "Tagkawayan"
## [1177] "TayabasCity" "Tiaong"
## [1179] "Unisan" "Aglipay"
## [1181] "Cabarroguis" "Diffun"
## [1183] "Maddela" "Nagtipunan"
## [1185] "Saguday" "Angono"
## [1187] "AntipoloCity" "Binangonan"
## [1189] "Cainta" "Cardona"
## [1191] "Jala-Jala" "Pililla"
## [1193] "Rodriguez" "Tanay"
## [1195] "Teresa" "Banton"
## [1197] "Cajidiocan" "Corcuera"
## [1199] "Ferrol" "Magdiwang"
## [1201] "Odiongan" "Romblon"
## [1203] "Almagro" "Basey"
## [1205] "CalbayogCity" "Calbiga"
## [1207] "CatbaloganCity" "Daram"
## [1209] "Gandara" "Hinabangan"
## [1211] "Jiabong" "Marabut"
## [1213] "Matuguinao" "Motiong"
## [1215] "Pagsanghan" "Paranas"
## [1217] "Pinabacdao" "SanJorge"
## [1219] "SanJosedeBuan" "SanSebastian"
## [1221] "SantaMargarita" "SantoNino"
## [1223] "Tagapul-An" "Talalora"
## [1225] "Tarangnan" "Villareal"
## [1227] "Zumarraga" "Alabel"
## [1229] "Glan" "Kiamba"
## [1231] "Maasim" "Maitum"
## [1233] "Malapatan" "Malungon"
## [1235] "EnriqueVillanueva" "Larena"
## [1237] "Lazi" "Maria"
## [1239] "Siquijor" "Barcelona"
## [1241] "Bulan" "Bulusan"
## [1243] "Castilla" "Donsol"
## [1245] "Gubat" "Irosin"
## [1247] "Juban" "Matnog"
## [1249] "PrietoDiaz" "SantaMagdalena"
## [1251] "SorsogonCity" "GeneralSantosCity"
## [1253] "KoronadalCity" "LakeSebu"
## [1255] "Norala" "Polomolok"
## [1257] "Surallah" "T'Boli"
## [1259] "Tampakan" "Tantangan"
## [1261] "Tupi" "Anahawan"
## [1263] "Hinundayan" "Libagon"
## [1265] "Limasawa" "MaasinCity"
## [1267] "Macrohon" "Pintuyan"
## [1269] "SaintBernard" "SanRicardo"
## [1271] "Silago" "TomasOppus"
## [1273] "Bagumbayan" "Columbio"
## [1275] "Isulan" "Kalamansig"
## [1277] "Lambayong" "Lebak"
## [1279] "Lutayan" "Palimbang"
## [1281] "PresidentQuirino" "Sen.NinoyAquino"
## [1283] "TacurongCity" "HadjiPanglimaTahil"
## [1285] "Indanan" "Jolo"
## [1287] "KalingalanCaluang" "Lugus"
## [1289] "Luuk" "Maimbung"
## [1291] "OldPanamao" "Pandami"
## [1293] "PanglimaEstino" "Pangutaran"
## [1295] "Pata" "Patikul"
## [1297] "Siasi" "Talipao"
## [1299] "Tapul" "Tongkil"
## [1301] "Bacuag" "Claver"
## [1303] "Dapa" "DelCarmen"
## [1305] "Gigaquit" "Mainit"
## [1307] "Malimono" "SanBenito"
## [1309] "SantaMonica" "SurigaoCity"
## [1311] "Tagana-An" "Barobo"
## [1313] "Bayabas" "BisligCity"
## [1315] "Cagwait" "Cantilan"
## [1317] "Carrascal" "Hinatuan"
## [1319] "Lanuza" "Lianga"
## [1321] "Lingig" "Madrid"
## [1323] "Marihatag" "Tagbina"
## [1325] "Tago" "TandagCity"
## [1327] "Anao" "Bamban"
## [1329] "Camiling" "Capas"
## [1331] "Gerona" "Mayantoc"
## [1333] "Moncada" "Paniqui"
## [1335] "Pura" "Ramos"
## [1337] "SanClemente" "SantaIgnacia"
## [1339] "TarlacCity" "Bongao"
## [1341] "Languyan" "Mapun"
## [1343] "PanglimaSugala" "Sapa-Sapa"
## [1345] "Sibutu" "Simunul"
## [1347] "Sitangkai" "SouthUbian"
## [1349] "Tandubas" "TurtleIslands"
## [1351] "Botolan" "Cabangan"
## [1353] "Castillejos" "Iba"
## [1355] "Masinloc" "OlongapoCity"
## [1357] "Palauig" "SanFelipe"
## [1359] "SanMarcelino" "Subic"
## [1361] "Bacungan" "Baliguian"
## [1363] "DapitanCity" "DipologCity"
## [1365] "Godod" "Gutalac"
## [1367] "JoseDalman" "Kalawit"
## [1369] "Katipunan" "Labason"
## [1371] "Liloy" "Manukan"
## [1373] "Mutia" "Pinan"
## [1375] "Polanco" "Pres.ManuelA.Roxas"
## [1377] "Salug" "SergioOsmenaSr."
## [1379] "Siayan" "Sibuco"
## [1381] "Sibutad" "Sindangan"
## [1383] "Siocon" "Sirawai"
## [1385] "Tampilisan" "Bayog"
## [1387] "Dimataling" "Dinas"
## [1389] "Dumalinao" "Dumingag"
## [1391] "Guipos" "Josefina"
## [1393] "Kumalarang" "Labangan"
## [1395] "Lakewood" "LakewoodLake"
## [1397] "Lapuyan" "Mahayag"
## [1399] "Margosatubig" "Midsalip"
## [1401] "Molave" "PagadianCity"
## [1403] "RamonMagsaysay" "Sominot"
## [1405] "Tabina" "Tambulig"
## [1407] "Tigbao" "Tukuran"
## [1409] "VincenzoA.Sagun" "ZamboangaCity"
## [1411] "Buug" "Diplahan"
## [1413] "Imelda" "Ipil"
## [1415] "Kabasalan" "Mabuhay"
## [1417] "Malangas" "Naga"
## [1419] "Olutanga" "Payao"
## [1421] "RosellerLim" "Siay"
## [1423] "Talusan" "Titay"
## [1425] "Tungawan"
# Filter for Davao City (all barangays)
davao_barangays <- ph_l3 %>% filter(NAME_2 == "DavaoCity")
# Check unique barangay names
unique(davao_barangays$NAME_3)
## [1] "Acacia" "Agdao" "Alambre"
## [4] "AlejandraNavarro" "AlfonsoAngliongtoS" "Angalan"
## [7] "Atan-Awe" "Baganihan" "BagoAplaya"
## [10] "BagoGallera" "BagoOshiro" "Baguio"
## [13] "Balengaeng" "Baliok" "BangkasHeights"
## [16] "Bantol" "Baracatan" "Barangay1-A"
## [19] "Barangay10-A" "Barangay11-B" "Barangay12-B"
## [22] "Barangay13-B" "Barangay14-B" "Barangay15-B"
## [25] "Barangay16-B" "Barangay17-B" "Barangay18-B"
## [28] "Barangay19-B" "Barangay2-A" "Barangay20-B"
## [31] "Barangay21-C" "Barangay22-C" "Barangay23-C"
## [34] "Barangay24-C" "Barangay25-C" "Barangay26-C"
## [37] "Barangay27-C" "Barangay28-C" "Barangay29-C"
## [40] "Barangay3-A" "Barangay30-C" "Barangay31-D"
## [43] "Barangay32-D" "Barangay33-D" "Barangay34-D"
## [46] "Barangay35-D" "Barangay36-D" "Barangay37-D"
## [49] "Barangay38-D" "Barangay39-D" "Barangay4-A"
## [52] "Barangay40-D" "Barangay5-A" "Barangay6-A"
## [55] "Barangay7-A" "Barangay8-A" "Barangay9-A"
## [58] "Bato" "Bayabas" "BiaoEscuela"
## [61] "BiaoGuianga" "BiaoJoaquin" "Binugao"
## [64] "Bucana" "Buda" "Buhangin"
## [67] "Bunawan" "Cabantian" "Cadalian"
## [70] "Calinan" "Callawa" "Camansi"
## [73] "Carmen" "CatalunanGrande" "CatalunanPequeño"
## [76] "Catigan" "Cawayan" "Centro"
## [79] "Colosas" "Communal" "CrossingBayabas"
## [82] "Dacudao" "Dalag" "Dalagdag"
## [85] "Daliao" "DaliaonPlantation" "DatuSalumay"
## [88] "Dominga" "Dumoy" "Eden"
## [91] "Fatima" "Gatungan" "Gov.PacianoBangoy"
## [94] "Gov.VicenteDuterte" "Gumalang" "Gumitan"
## [97] "Ilang" "Inayangan" "Indangan"
## [100] "Kap.TomasMonteverd" "Kilate" "Lacson"
## [103] "Lamanan" "Lampianao" "Langub"
## [106] "Lapu-Lapu" "LeonGarcia,Sr." "Lizada"
## [109] "LosAmigos" "Lubogan" "Lumiad"
## [112] "Ma-A" "Mabuhay" "Magsaysay"
## [115] "Magtuod" "Mahayag" "Malabog"
## [118] "Malagos" "Malamba" "Manambulan"
## [121] "Mandug" "ManuelGuianga" "Mapula"
## [124] "Marapangi" "Marilog" "MatinaAplaya"
## [127] "MatinaBiao" "MatinaCrossing" "MatinaPangi"
## [130] "Megkawayan" "Mintal" "Mudiang"
## [133] "Mulig" "NewCarmen" "NewValencia"
## [136] "Pampanga" "Panacan" "Panalum"
## [139] "Pandaitan" "Pangyan" "Paquibato"
## [142] "ParadiseEmbak" "RafaelCastillo" "Riverside"
## [145] "Salapawan" "Salaysay" "Saloy"
## [148] "SanAntonio" "SanIsidro" "SantoNiño"
## [151] "Sasa" "Sibulan" "Sirawan"
## [154] "Sirib" "Suawan" "Subasta"
## [157] "Sumimao" "Tacunan" "Tagakpan"
## [160] "Tagluno" "Tagurano" "Talandang"
## [163] "Talomo" "TalomoRiver" "Tamayong"
## [166] "Tambobong" "Tamugan" "Tapak"
## [169] "Tawan-Tawan" "Tibuloy" "Tibungco"
## [172] "Tigatto" "Toril" "Tugbok"
## [175] "Tungakalan" "Ubalde" "Ula"
## [178] "VicenteHizonSr." "Waan" "Wangan"
## [181] "WilfredoAquino" "Wines"
ggplot() +
geom_sf(data = davao_barangays, fill = "lightgray", color = "black") +
theme_minimal() +
labs(title = "Barangay-Level Map of Davao City")
# Select barangays to highlight
selected_barangays <- davao_barangays %>% filter(NAME_3 %in% c("Buhangin", "Agdao"))
# Plot with highlights
ggplot() +
geom_sf(data = davao_barangays, fill = "lightgray", color = "black") + # Other barangays in gray
geom_sf(data = selected_barangays, fill = "red", color = "black") + # Highlighted barangays in red
theme_minimal() +
labs(title = "Highlighted Barangays in Davao City")
Plotting Population density in the City Map
pop_data <- read.csv("C:/Users/User/Desktop/Spatial Analysis/pop.csv")
pop_data <- pop_data %>%
mutate(pop_density = population / area_km2)
davao_map <- davao_barangays %>%
left_join(pop_data, by = "NAME_3")
ggplot(davao_map) +
geom_sf(aes(fill = pop_density), color = "black") +
scale_fill_viridis_c(option = "magma",direction=-1, name = "Pop. Density (per km²)") +
theme_minimal() +
labs(title = "Population Density of Davao City Barangays")
scale_fill_viridis_c(option = "plasma", name = "Pop. Density")
## <ScaleContinuous>
## Range:
## Limits: 0 -- 1
theme_minimal() +
theme(
plot.title = element_text(size = 16, face = "bold"),
legend.position = "bottom"
)
## List of 136
## $ line :List of 6
## ..$ colour : chr "black"
## ..$ linewidth : num 0.5
## ..$ linetype : num 1
## ..$ lineend : chr "butt"
## ..$ arrow : logi FALSE
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_line" "element"
## $ rect :List of 5
## ..$ fill : chr "white"
## ..$ colour : chr "black"
## ..$ linewidth : num 0.5
## ..$ linetype : num 1
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_rect" "element"
## $ text :List of 11
## ..$ family : chr ""
## ..$ face : chr "plain"
## ..$ colour : chr "black"
## ..$ size : num 11
## ..$ hjust : num 0.5
## ..$ vjust : num 0.5
## ..$ angle : num 0
## ..$ lineheight : num 0.9
## ..$ margin : 'margin' num [1:4] 0points 0points 0points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : logi FALSE
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ title : NULL
## $ aspect.ratio : NULL
## $ axis.title : NULL
## $ axis.title.x :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 1
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 2.75points 0points 0points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.title.x.top :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 0
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 0points 2.75points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.title.x.bottom : NULL
## $ axis.title.y :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 1
## ..$ angle : num 90
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 2.75points 0points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.title.y.left : NULL
## $ axis.title.y.right :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 1
## ..$ angle : num -90
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 0points 0points 2.75points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : chr "grey30"
## ..$ size : 'rel' num 0.8
## ..$ hjust : NULL
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : NULL
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text.x :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 1
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 2.2points 0points 0points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text.x.top :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : NULL
## ..$ vjust : num 0
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 0points 2.2points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text.x.bottom : NULL
## $ axis.text.y :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : num 1
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 2.2points 0points 0points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text.y.left : NULL
## $ axis.text.y.right :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : num 0
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 0points 0points 2.2points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.text.theta : NULL
## $ axis.text.r :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : num 0.5
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : 'margin' num [1:4] 0points 2.2points 0points 2.2points
## .. ..- attr(*, "unit")= int 8
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ axis.ticks : list()
## ..- attr(*, "class")= chr [1:2] "element_blank" "element"
## $ axis.ticks.x : NULL
## $ axis.ticks.x.top : NULL
## $ axis.ticks.x.bottom : NULL
## $ axis.ticks.y : NULL
## $ axis.ticks.y.left : NULL
## $ axis.ticks.y.right : NULL
## $ axis.ticks.theta : NULL
## $ axis.ticks.r : NULL
## $ axis.minor.ticks.x.top : NULL
## $ axis.minor.ticks.x.bottom : NULL
## $ axis.minor.ticks.y.left : NULL
## $ axis.minor.ticks.y.right : NULL
## $ axis.minor.ticks.theta : NULL
## $ axis.minor.ticks.r : NULL
## $ axis.ticks.length : 'simpleUnit' num 2.75points
## ..- attr(*, "unit")= int 8
## $ axis.ticks.length.x : NULL
## $ axis.ticks.length.x.top : NULL
## $ axis.ticks.length.x.bottom : NULL
## $ axis.ticks.length.y : NULL
## $ axis.ticks.length.y.left : NULL
## $ axis.ticks.length.y.right : NULL
## $ axis.ticks.length.theta : NULL
## $ axis.ticks.length.r : NULL
## $ axis.minor.ticks.length : 'rel' num 0.75
## $ axis.minor.ticks.length.x : NULL
## $ axis.minor.ticks.length.x.top : NULL
## $ axis.minor.ticks.length.x.bottom: NULL
## $ axis.minor.ticks.length.y : NULL
## $ axis.minor.ticks.length.y.left : NULL
## $ axis.minor.ticks.length.y.right : NULL
## $ axis.minor.ticks.length.theta : NULL
## $ axis.minor.ticks.length.r : NULL
## $ axis.line : list()
## ..- attr(*, "class")= chr [1:2] "element_blank" "element"
## $ axis.line.x : NULL
## $ axis.line.x.top : NULL
## $ axis.line.x.bottom : NULL
## $ axis.line.y : NULL
## $ axis.line.y.left : NULL
## $ axis.line.y.right : NULL
## $ axis.line.theta : NULL
## $ axis.line.r : NULL
## $ legend.background : list()
## ..- attr(*, "class")= chr [1:2] "element_blank" "element"
## $ legend.margin : 'margin' num [1:4] 5.5points 5.5points 5.5points 5.5points
## ..- attr(*, "unit")= int 8
## $ legend.spacing : 'simpleUnit' num 11points
## ..- attr(*, "unit")= int 8
## $ legend.spacing.x : NULL
## $ legend.spacing.y : NULL
## $ legend.key : list()
## ..- attr(*, "class")= chr [1:2] "element_blank" "element"
## $ legend.key.size : 'simpleUnit' num 1.2lines
## ..- attr(*, "unit")= int 3
## $ legend.key.height : NULL
## $ legend.key.width : NULL
## $ legend.key.spacing : 'simpleUnit' num 5.5points
## ..- attr(*, "unit")= int 8
## $ legend.key.spacing.x : NULL
## $ legend.key.spacing.y : NULL
## $ legend.frame : NULL
## $ legend.ticks : NULL
## $ legend.ticks.length : 'rel' num 0.2
## $ legend.axis.line : NULL
## $ legend.text :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : 'rel' num 0.8
## ..$ hjust : NULL
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : NULL
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ legend.text.position : NULL
## $ legend.title :List of 11
## ..$ family : NULL
## ..$ face : NULL
## ..$ colour : NULL
## ..$ size : NULL
## ..$ hjust : num 0
## ..$ vjust : NULL
## ..$ angle : NULL
## ..$ lineheight : NULL
## ..$ margin : NULL
## ..$ debug : NULL
## ..$ inherit.blank: logi TRUE
## ..- attr(*, "class")= chr [1:2] "element_text" "element"
## $ legend.title.position : NULL
## $ legend.position : chr "bottom"
## $ legend.position.inside : NULL
## $ legend.direction : NULL
## $ legend.byrow : NULL
## $ legend.justification : chr "center"
## $ legend.justification.top : NULL
## $ legend.justification.bottom : NULL
## $ legend.justification.left : NULL
## $ legend.justification.right : NULL
## $ legend.justification.inside : NULL
## $ legend.location : NULL
## $ legend.box : NULL
## $ legend.box.just : NULL
## $ legend.box.margin : 'margin' num [1:4] 0cm 0cm 0cm 0cm
## ..- attr(*, "unit")= int 1
## $ legend.box.background : list()
## ..- attr(*, "class")= chr [1:2] "element_blank" "element"
## $ legend.box.spacing : 'simpleUnit' num 11points
## ..- attr(*, "unit")= int 8
## [list output truncated]
## - attr(*, "class")= chr [1:2] "theme" "gg"
## - attr(*, "complete")= logi TRUE
## - attr(*, "validate")= logi TRUE
ggplot(davao_map) +
geom_sf(aes(fill = pop_density), color = "black") +
geom_text(data = davao_map, aes(label = NAME_3, geometry = geometry),
stat = "sf_coordinates", size = 3, check_overlap = TRUE) +
scale_fill_viridis_c(option = "magma", name = "Pop. Density") +
theme_minimal() +
labs(title = "Population Density of Davao City Barangays")
## Warning in st_point_on_surface.sfc(sf::st_zm(x)): st_point_on_surface may not
## give correct results for longitude/latitude data
