update freetwon thematic analysis report #113

Author: Saif Shabou

baseline-indicators/biodiversity/scripts/biodiverisy_data_exploration.R

@@ -318,14 +318,17 @@ birds_in_geo_adm2_stat %>%
          yaxis = list(title = 'Number of observations in 2020'))
 
 
-# protected areas -----
+# protected areas map -----
 
 protected_areas_sanjose = st_read("./data/biodiversity/data/biodiversity_map_layers_protected_areas_sanjose.geojson",
                                 quiet = TRUE)
 
 protected_areas_costarica = st_read("./data/biodiversity/data/biodiversity_map_layers_protected_areas_costarica.geojson",
                                   quiet = TRUE)
 
+protected_areas_costarica = st_read("https://cities-urbanshift.s3.eu-west-3.amazonaws.com/baseline-indicators/biodiversity/data/biodiversity_map_layers_protected_areas_sanjose.geojson",
+                                    quiet = TRUE)
+
 leaflet(data = boundary, height = 500, width = "100%") %>% 
   addTiles() %>%
   addProviderTiles(providers$Stamen.TonerLite, group = "Toner Lite") %>%
@@ -335,7 +338,7 @@ leaflet(data = boundary, height = 500, width = "100%") %>%
   addPolygons(data = boundary,
               group = "Administrative boundaries",
               stroke = TRUE, color = col_Black, weight = 3,dashArray = "3",
-              smoothFactor = 0.3, fill = TRUE, fillOpacity = 0.2,
+              smoothFactor = 0.3, fill = FALSE, fillOpacity = 0.2,
               highlight = highlightOptions(
                 weight = 5,
                 color = "#666",
@@ -372,6 +375,27 @@ leaflet(data = boundary, height = 500, width = "100%") %>%
   ) %>% 
   hideGroup("Municipality boundaries")
 
-# percent of protected areas by municipality
+# percent of protected areas by municipality ----
+
+protectedarea_percent_by_municipality = read.csv("./data/biodiversity/data/protectedarea_percent_by_municipality.csv")
 
+protectedarea_percent_by_municipality = read.csv("https://cities-urbanshift.s3.eu-west-3.amazonaws.com/baseline-indicators/biodiversity/data/protectedarea_percent_by_municipality.csv")
 
+# rename columns
+protectedarea_percent_by_municipality = protectedarea_percent_by_municipality %>% 
+  rename_with(.cols = 1, ~'municpality_name') %>% 
+  rename_with(.cols = 2, ~'percent_protected_area')
+
+# plot chart
+protectedarea_percent_by_municipality %>% 
+  arrange(desc(percent_protected_area)) %>%
+  plot_ly(height = 500, width = 900) %>% 
+  add_trace(x = ~factor(municpality_name), 
+            y = ~percent_protected_area, 
+            marker = list(color = col_BoldRiverBlue),
+            type = "bar",
+            orientation = "v")  %>% 
+  layout(title = "Percent of protected area by municipality (2020)",
+         xaxis = list(title = '', categoryorder = "array",categoryarray = ~percent_protected_area),
+         yaxis = list(title = 'Percent of protected area (%)'))
+  

baseline-indicators/greenspace/scripts/merge_outputs.R

@@ -21,7 +21,7 @@ path = paste(getwd(),"github/cities-urbanshift/baseline-indicators/greenspace",s
 
 # Define city
 
-selected_city = "RWA-Kigali"
+selected_city = "CRI-San_Jose" # "SLE-Freetown_core" #"RWA-Kigali"
 
 
 # read outputs indicators

thematic-analysis/Freetown/Freetown-thematic-analysis.Rmd

@@ -970,7 +970,11 @@ These analyses provide a proof of concept of using global remote sensing open da
 
 The trained model presents a high accuracy (95%) based on the validation set extracted from the labeled data. The application of this model enabled us to explore the spatial distribution of `urban` and `non urban` areas in the whole Freetown region at different periods. The outputs show an increase in urban expansion encroaching on the Western Area Forest Reserve. 
 
-We proposed an urban expansion indicator based on the percentage of urban area between between 2016 and 2021. Aggregated at the ward level, the proposed indicator highlighted the cities that present the highest urban expansion, such as `Thnderhill` and `Mamba Ridge II` (more than 20% increase of urban areas in the last 6 years).
+- **Wards level:**  We proposed an urban expansion indicator based on the percentage of urban area between between 2016 and 2021. Aggregated at the ward level, the proposed indicator highlighted the cities that present the highest urban expansion, such as `Thnderhill` and `Mamba Ridge II` (more than 20% increase of urban areas in the last 6 years).
+
+- **Planning areas**: FCC provided us with the boundaries of planning areas used for urban planning policies. The aggregation of urban expansion indicator at the planning areas level highlighted the importance of urban expansion in two planning areas: `zone 8` and `zone 9` with respectively `19.3%` and `12.9%` increase in urban area.
+
+- **Urban expansion by land use class:** FCC provided us with the land use shapefile dividing Freetown municipality into 13 distinct land classes. The aggregation of urban expansion indicator by land use class shows that the urban expansion is mainly located in `waterbody`, `woodland` and `urban agriculture` areas.
 
 The results of this analysis may be used to better plan prevent urban expansion in the identified areas and localities as part of the development of Freetown master plan.
 
@@ -1001,6 +1005,23 @@ While the results of the proposed classification model provide insights on urban
   - [Labeled urban areas, geojson file](https://storage.googleapis.com/urbanshift/freetown/labeled_data/urban.geojson)
   - [Labeled non-urban areas, geojson file](https://storage.googleapis.com/urbanshift/freetown/labeled_data/non_urban.geojson)
 
-- **Classification maps**:
+- **Classification maps (rasters)**:
   - [Freetown classification map 2016, GoeTiff file](https://storage.googleapis.com/urbanshift/freetown/classified_maps/classified_map_2016.tif)
   - [Freetown classification map 2021,  GoeTiff file](https://storage.googleapis.com/urbanshift/freetown/classified_maps/classified_map_2021.tif)
+
+## Download output maps
+
+
+- [Land classification 2016, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/land-classification-2016.png)
+
+- [Land classification 2021, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/land-classification-2021.png)
+
+- [Urban expansion - 2016-2021 - region, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/urban-expansion-2016-2021.png)
+
+- [Urban expansion - 2016-2021 - wards, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/urban-expansion-2016-2021-wards.png)
+
+- [Urban expansion - 2016-2021 - planning areas, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/urban-expansion-2016-2021-planning-areas.png)
+
+- [Land use map, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/land-use-map.png)
+
+- [Urban expansion - 2016-2021 - land use class, png file](https://cities-urbanshift.s3.eu-west-3.amazonaws.com/thematic-analysis/freetown/maps/urban-expansion-2016-2021-by-land-use-classes.png)