library(sf)
## Linking to GEOS 3.12.1, GDAL 3.8.4, PROJ 9.3.1; sf_use_s2() is TRUE
library(tmap)
## Breaking News: tmap 3.x is retiring. Please test v4, e.g. with
## remotes::install_github('r-tmap/tmap')
library(tigris)
## To enable caching of data, set `options(tigris_use_cache = TRUE)`
## in your R script or .Rprofile.
library(tidycensus)
#Q1
var <- c(poptotal='B03002_001E',
hispanic='B03002_012E')
st <-"WA"
ct <-"King"
df <- get_acs(geography = "tract", variables = var, count=ct,
state = st, output ="wide", year = 2021, geometry = TRUE)
## Getting data from the 2017-2021 5-year ACS
## Downloading feature geometry from the Census website. To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.
## Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 15 kB Downloading: 31 kB Downloading: 31 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 35 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 52 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 65 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 83 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 99 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 110 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 130 kB Downloading: 150 kB Downloading: 150 kB Downloading: 150 kB Downloading: 150 kB Downloading: 150 kB Downloading: 150 kB Downloading: 160 kB Downloading: 160 kB Downloading: 160 kB Downloading: 160 kB Downloading: 160 kB Downloading: 160 kB Downloading: 180 kB Downloading: 180 kB Downloading: 190 kB Downloading: 190 kB Downloading: 190 kB Downloading: 190 kB Downloading: 190 kB Downloading: 190 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 200 kB Downloading: 220 kB Downloading: 220 kB Downloading: 220 kB Downloading: 220 kB Downloading: 220 kB Downloading: 220 kB Downloading: 230 kB Downloading: 230 kB Downloading: 240 kB Downloading: 240 kB Downloading: 240 kB Downloading: 240 kB Downloading: 250 kB Downloading: 250 kB Downloading: 250 kB Downloading: 250 kB Downloading: 250 kB Downloading: 250 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 270 kB Downloading: 280 kB Downloading: 280 kB Downloading: 290 kB Downloading: 290 kB Downloading: 300 kB Downloading: 300 kB Downloading: 300 kB Downloading: 300 kB Downloading: 320 kB Downloading: 320 kB Downloading: 320 kB Downloading: 320 kB Downloading: 320 kB Downloading: 320 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 330 kB Downloading: 350 kB Downloading: 350 kB Downloading: 350 kB Downloading: 350 kB Downloading: 370 kB Downloading: 370 kB Downloading: 380 kB Downloading: 380 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 400 kB Downloading: 420 kB Downloading: 420 kB Downloading: 420 kB Downloading: 420 kB Downloading: 430 kB Downloading: 430 kB Downloading: 430 kB Downloading: 430 kB Downloading: 450 kB Downloading: 450 kB Downloading: 460 kB Downloading: 460 kB Downloading: 460 kB Downloading: 460 kB Downloading: 470 kB Downloading: 470 kB Downloading: 490 kB Downloading: 490 kB Downloading: 500 kB Downloading: 500 kB Downloading: 500 kB Downloading: 500 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 520 kB Downloading: 530 kB Downloading: 530 kB Downloading: 530 kB Downloading: 530 kB Downloading: 550 kB Downloading: 550 kB Downloading: 550 kB Downloading: 550 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 590 kB Downloading: 610 kB Downloading: 610 kB Downloading: 610 kB Downloading: 610 kB Downloading: 640 kB Downloading: 640 kB Downloading: 650 kB Downloading: 650 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 660 kB Downloading: 670 kB Downloading: 670 kB Downloading: 670 kB Downloading: 670 kB Downloading: 680 kB Downloading: 680 kB Downloading: 700 kB Downloading: 700 kB Downloading: 700 kB Downloading: 700 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 710 kB Downloading: 720 kB Downloading: 720 kB Downloading: 740 kB Downloading: 740 kB Downloading: 760 kB Downloading: 760 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 780 kB Downloading: 800 kB Downloading: 800 kB Downloading: 820 kB Downloading: 820 kB Downloading: 830 kB Downloading: 830 kB Downloading: 830 kB Downloading: 830 kB Downloading: 970 kB Downloading: 970 kB Downloading: 970 kB Downloading: 970 kB Downloading: 970 kB Downloading: 970 kB Downloading: 990 kB Downloading: 990 kB Downloading: 990 kB Downloading: 990 kB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.1 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.2 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.3 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB Downloading: 1.4 MB
df$hispanic_pct <- df$hispanic/df$poptotal
df <- na.omit(df)
tm_shape(df)+
tm_polygons(col = "hispanic_pct", title = "Hispanic %")+
tm_compass(type = "8star", position = c("right", "top"), size = 1)+
tm_scale_bar(position = c("left", "bottom"))+
tm_layout(legend.position = c("center", "bottom"),
legend.outside = TRUE,
legend.outside.position = "right")
#Q2
ADA_ramp<-st_read(dsn = "C:/Users/dejmo/OneDrive/Documents/URP Cert/URP 5363- Planning Methods I/Project 2/ADA_(American_with_Disabilities_Act)_Curb_Ramps___ada_ramp_point.shp")
## Reading layer `ADA_(American_with_Disabilities_Act)_Curb_Ramps___ada_ramp_point' from data source `C:\Users\dejmo\OneDrive\Documents\URP Cert\URP 5363- Planning Methods I\Project 2\ADA_(American_with_Disabilities_Act)_Curb_Ramps___ada_ramp_point.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 4676 features and 26 fields
## Geometry type: POINT
## Dimension: XY
## Bounding box: xmin: 1225426 ymin: 70955.09 xmax: 1411831 ymax: 287486.6
## Projected CRS: NAD83(HARN) / Washington North (ftUS)
king_tracts<- tracts(state = "WA", county = "King", cb=T)
## Retrieving data for the year 2022
## | | | 0% | |== | 2% | |=== | 4% | |===== | 7% | |======= | 10% | |========= | 13% | |=========== | 16% | |================================== | 48% | |==================================== | 51% | |====================================== | 54% | |======================================= | 55% | |========================================== | 59% | |============================================ | 62% | |============================================== | 65% | |================================================ | 68% | |=================================================== | 72% | |===================================================== | 75% | |======================================================= | 78% | |========================================================= | 82% | |============================================================ | 85% | |============================================================= | 88% | |================================================================ | 91% | |================================================================== | 95% | |======================================================================| 100%
tm_shape(king_tracts)+
tm_borders()+
tm_shape(ADA_ramp)+
tm_dots(col="red", title = "ADA Ramp Location") +
tm_layout(frame = FALSE)
#Q3
metro_transit_routes<-st_read(dsn = "C:/Users/dejmo/OneDrive/Documents/URP Cert/URP 5363- Planning Methods I/Project 2/Transit_Routes_for_King_County_Metro___transitroute_line.shp")
## Reading layer `Transit_Routes_for_King_County_Metro___transitroute_line' from data source `C:\Users\dejmo\OneDrive\Documents\URP Cert\URP 5363- Planning Methods I\Project 2\Transit_Routes_for_King_County_Metro___transitroute_line.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 146 features and 9 fields
## Geometry type: MULTILINESTRING
## Dimension: XY
## Bounding box: xmin: 1225211 ymin: 70794.9 xmax: 1405246 ymax: 320456.3
## Projected CRS: NAD83(HARN) / Washington North (ftUS)
tm_shape(king_tracts)+
tm_borders()+
tm_shape(metro_transit_routes)+
tm_lines(lwd=2, col = "blue") +
tm_layout(frame = FALSE)
#Q4
parks<-st_read(dsn = "C:/Users/dejmo/OneDrive/Documents/URP Cert/URP 5363- Planning Methods I/Project 2/Parks_in_King_County___park_area.shp")
## Reading layer `Parks_in_King_County___park_area' from data source
## `C:\Users\dejmo\OneDrive\Documents\URP Cert\URP 5363- Planning Methods I\Project 2\Parks_in_King_County___park_area.shp'
## using driver `ESRI Shapefile'
## Warning in CPL_read_ogr(dsn, layer, query, as.character(options), quiet, : GDAL
## Message 1: organizePolygons() received an unexpected geometry. Either a
## polygon with interior rings, or a polygon with less than 4 points, or a
## non-Polygon geometry. Return arguments as a collection.
## Warning in CPL_read_ogr(dsn, layer, query, as.character(options), quiet, : GDAL
## Message 1: Geometry of polygon of fid 1124 cannot be translated to Simple
## Geometry. All polygons will be contained in a multipolygon.
## Warning in CPL_read_ogr(dsn, layer, query, as.character(options), quiet, : GDAL
## Message 1: organizePolygons() received an unexpected geometry. Either a
## polygon with interior rings, or a polygon with less than 4 points, or a
## non-Polygon geometry. Return arguments as a collection.
## Warning in CPL_read_ogr(dsn, layer, query, as.character(options), quiet, : GDAL
## Message 1: Geometry of polygon of fid 1441 cannot be translated to Simple
## Geometry. All polygons will be contained in a multipolygon.
## Simple feature collection with 1467 features and 13 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 1220335 ymin: 53899.29 xmax: 1583156 ymax: 287209.5
## Projected CRS: NAD83(HARN) / Washington North (ftUS)
#set tmap to automatically check and fix geometries
tmap_options(check.and.fix = TRUE)
sf::st_is_valid(parks)
## [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [13] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [25] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [37] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [49] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [61] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [73] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [85] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [97] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [109] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [121] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [133] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [145] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [157] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [169] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [181] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [193] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [205] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [217] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [229] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [241] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [253] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [265] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [277] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [289] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [301] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [313] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [325] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [337] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [349] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [361] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [373] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [385] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [397] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [409] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [421] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [433] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [445] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [457] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [469] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [481] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [493] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [505] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [517] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [529] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [541] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [553] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [565] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [577] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [589] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [601] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [613] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [625] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [637] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [649] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [661] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [673] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [685] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [697] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [709] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [721] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [733] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [745] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [757] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [769] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [781] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [793] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [805] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [817] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [829] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [841] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [853] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [865] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [877] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [889] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [901] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [913] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [925] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [937] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [949] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [961] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [973] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [985] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [997] TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1009] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1021] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1033] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1045] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1057] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1069] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1081] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1093] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1105] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1117] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE FALSE TRUE TRUE TRUE
## [1129] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1141] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1153] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1165] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1177] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1189] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1201] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1213] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1225] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1237] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1249] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1261] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1273] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1285] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1297] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1309] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1321] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1333] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1345] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1357] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1369] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1381] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1393] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1405] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1417] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1429] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1441] TRUE FALSE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1453] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
## [1465] TRUE TRUE TRUE
sf::st_make_valid(parks)
## Simple feature collection with 1467 features and 13 fields
## Geometry type: GEOMETRY
## Dimension: XY
## Bounding box: xmin: 1220335 ymin: 53899.29 xmax: 1583156 ymax: 287209.5
## Projected CRS: NAD83(HARN) / Washington North (ftUS)
## First 10 features:
## OBJECTID KC_FAC_FID KCPARKFID SITENAME SITETYPE
## 1 1 2662 6223001 Jenkins Creek Park Park Site
## 2 2 2805 5250702 Street End Park Park Site
## 3 3 5069 4261003 Blue Heron Park Park Site
## 4 4 5080 7242602 Azalea Park Park Site
## 5 5 2814 5250503 Tot Lot Park Park Site
## 6 6 7821 4231802 Linde Hill Park Park Site
## 7 7 8227 4261004 Whispering Willow Park Park Site
## 8 8 7217 4230203 Thayer Place Park Site
## 9 9 7133 3252406 Coe Play Park Park Site
## 10 10 2529 4232501 Bicentennial Park - Tukwila Park Site
## OWNER OWNERTYPE MANAGER MANAGETYPE
## 1 City of Covington City City of Covington City
## 2 City of Kirkland City City of Kirkland City
## 3 City of Lake Forest Park City City of Lake Forest Park City
## 4 City of Snoqualmie City City of Snoqualmie City
## 5 City of Kirkland City City of Kirkland City
## 6 City of Burien City City of Burien City
## 7 City of Lake Forest Park City City of Lake Forest Park City
## 8 City of Seattle City City of Seattle City
## 9 City of Seattle City City of Seattle City
## 10 City of Tukwila City City of Tukwila City
## MAINTD_BY MAINTTYPE SHAPE_Leng SHAPE_Area
## 1 City of Covington City 5632.3248 965530.821
## 2 City of Kirkland City 1489.6380 41080.857
## 3 City of Lake Forest Park City 746.2705 32094.864
## 4 City of Snoqualmie City 1225.8320 87437.482
## 5 City of Kirkland City 620.0069 22800.505
## 6 City of Burien City 846.9323 41232.508
## 7 City of Lake Forest Park City 729.5020 26398.692
## 8 City of Seattle City 239.6966 3341.999
## 9 City of Seattle City 341.9664 7012.632
## 10 City of Tukwila City 1324.8975 49923.600
## geometry
## 1 POLYGON ((1327665 135120.9,...
## 2 POLYGON ((1302344 248195.4,...
## 3 POLYGON ((1284219 277618.1,...
## 4 POLYGON ((1382829 199424.9,...
## 5 POLYGON ((1302069 251672.8,...
## 6 POLYGON ((1266982 177656.5,...
## 7 POLYGON ((1284264 278048.6,...
## 8 POLYGON ((1290335 189511.4,...
## 9 POLYGON ((1262672 237215.8,...
## 10 POLYGON ((1290668 169833.3,...
tm_shape(king_tracts)+
tm_borders()+
tm_shape(parks)+
tm_polygons(col = "green") +
tm_layout(frame = FALSE)
## Warning: The shape parks is invalid (after reprojection). See sf::st_is_valid
#Warning message: The shape parks is invalid (after reprojection). See sf::st_is_valid
#Q5
metro_transit_routes <- st_make_valid(metro_transit_routes)
parks <- st_make_valid(parks)
metro_transit_routes<-st_transform(metro_transit_routes, crs = 4326)
parks<-st_transform(parks, crs = 4326)
clipped_routes <- st_intersection(metro_transit_routes, parks)
## Warning: attribute variables are assumed to be spatially constant throughout
## all geometries
#Q7
library(raster)
## Loading required package: sp
DEM_SW <- raster("C:/Users/dejmo/OneDrive/Documents/URP Cert/URP 5363- Planning Methods I/Project 2/USGS_1_n48w123_20240327.tif")
DEM_SW <- setMinMax(DEM_SW)
cellStats(DEM_SW, max)
## [1] 1346.843
# view raster attributes
DEM_SW
## class : RasterLayer
## dimensions : 3612, 3612, 13046544 (nrow, ncol, ncell)
## resolution : 0.0002777778, 0.0002777778 (x, y)
## extent : -123.0017, -121.9983, 46.99833, 48.00167 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=NAD83 +no_defs
## source : USGS_1_n48w123_20240327.tif
## names : Layer_1
## values : -13.81634, 1346.843 (min, max)
#Get min and max cell values from raster
#NOTE: this code may fail if the raster is too large
cellStats(DEM_SW, min)
## [1] -13.81634
cellStats(DEM_SW, max)
## [1] 1346.843
#Q8
plot(DEM_SW < 350, main = "Elevation Criteria", col = c("#ffffff", "#0000ff"))
#Q9
library(ggplot2)
library(rasterVis)
## Warning: package 'rasterVis' was built under R version 4.4.2
## Loading required package: lattice
plot(terrain(DEM_SW, opt = "slope", unit = "degrees") < 3,
main = "Slope Criteria",
col = c("#ffffff", "#ff9900"))
#Q10
plot(terrain(DEM_SW, opt = "aspect", unit = "degrees") >22.5 &
terrain(DEM_SW, opt = "aspect", unit = "degrees") <157.5,
main = "Aspect Criteria")
#Q11 - Writeup
#Very few of King County, WA's census tracts show a concentration of Hispanic population living in poverty above 50%
#The initial focus during Project 1 was King County, WA's disabled population, specifically those with mobility and vision impairments. I continued this focus during Project 2 and obtained GIS data pertaining to ADA ramps, metro transit routes and their relationship with existing parks as all individuals should have access to recreational areas. Available data shows King County has both an extensive transit system plus numerous parks throughout but ADA ramp presence is concentrated in clusters in specific areas.
#Based on available DEM data for King County, WA, the min. elevation is -13.81634 while the max elevation is 1346.843.
knitr::opts_chunk$set(echo = TRUE)
This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
summary(cars)
## speed dist
## Min. : 4.0 Min. : 2.00
## 1st Qu.:12.0 1st Qu.: 26.00
## Median :15.0 Median : 36.00
## Mean :15.4 Mean : 42.98
## 3rd Qu.:19.0 3rd Qu.: 56.00
## Max. :25.0 Max. :120.00
You can also embed plots, for example:
Note that the echo = FALSE parameter was added to the
code chunk to prevent printing of the R code that generated the
plot.