Spatial Mapping in R

Load all packages and libraries

install.packages('rasterVis')
install.packages('tmap')
install.packages('tigris')
install.packages('rgdal', type = "source", configure.args=c('--with-proj-include=/Library/Frameworks/PROJ.framework/Headers', '--with-proj-lib=/Library/Frameworks/PROJ.framework/unix/lib'))
install.packages("tmaptools")
install.packages("sf")
install.packages("leaflet")
library(sp)
library(tmap)
library(raster)
library(RColorBrewer)
library(classInt)
library(rasterVis)
library(tigris)
library(rgdal)
library("tmaptools")
library("sf")
library("leaflet")

Viz1: Create Ventura and LA County baselayer maps. I used library(tigris) to download LA and Ventura county census tract boundaries.Then I downloaded census level median home value data from American Factfinder. The last step was to use the merge function to join the census level data to the empty county boundaries. After that is was just all about fine-tuning the map. The interactive mapping mode fixes the problem of overcrowding the map.

Map1

legend.postion is used for plot mode. Use view.legend.position in tm_view to set the legend position in view mode.
legend.postion is used for plot mode. Use view.legend.position in tm_view to set the legend position in view mode.

Viz2: Add Hill and Woolsey Fire extent. For Viz2, I thought it would be interesting to add the extents of the Woosley fire to see the type of homes that were being affected by the fire. This could also be helpful in channeling recovery efforts since it would allow aid to be distributed to lower value housing areas first (assuming that’s the goal). For anyone interested: I had to scour the far reaches of the internet to find a downloadable fire extent data, not sure why it’s so hard. But I eventually found this google map: https://www.google.com/maps/d/u/0/viewer?mid=1IkW-WhQoRcP4WGIIP-h_v6WK8zb0jG3M&ll=34.126688966793836%2C-118.7345541136649&z=11. Then I downloaded a KML version of the extents then converted it using an online converter: https://mygeodata.cloud/conversion. This was all done outside of R but I figured this method was more realistic of a workflow I’d have outside this class.

Map1+
tm_shape(fire) +
  tm_fill(col = "Name", alpha = 0.8, palette = "OrRd", title = "Name of Fire") +
  tm_borders(alpha=0.5, lwd = 2, col = "red")+
  tm_scale_bar()

legend.postion is used for plot mode. Use view.legend.position in tm_view to set the legend position in view mode.
legend.postion is used for plot mode. Use view.legend.position in tm_view to set the legend position in view mode.

Viz3: Recreating GEOG7 Map. I thought it would be a good exercise to try to recreate one of the earlier maps I did for GEOG7 on R since it’ll be good practice for manipulating map elements plus it was kinda fun and nostalgic for me.

*It’s not quite exact as the reference: https://imgur.com/Mk6rozp but I think I got pretty close.