Task 2 - Read the Postal Codes into sf objects
nyc_zipcode_sf <- st_read("data/nyc/ZIP_CODE_040114.shp")
## Reading layer `ZIP_CODE_040114' from data source
## `C:\Users\elish\OneDrive - Hunter - CUNY\GTECH38520\R-spatial\data\nyc\ZIP_CODE_040114.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 263 features and 12 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: 913129 ymin: 120020.9 xmax: 1067494 ymax: 272710.9
## Projected CRS: NAD83 / New York Long Island (ftUS)
Task 3 - Read NYS health facilities data into sf objects
nys_healthfacilities <- read_csv("NYS_Health_Facility.csv",show_col_types = FALSE,
lazy = FALSE,)
# Remove missing values
nys_healthfacilities <- nys_healthfacilities %>%
filter(!is.na(`Facility Latitude`) & !is.na(`Facility Longitude`))
# Geographic coordinates
nys_healthfacilities_sf <- st_as_sf(
nys_healthfacilities,
coords = c("Facility Longitude", "Facility Latitude"),
crs = 4326
)
str(nys_healthfacilities_sf)
## sf [3,848 × 35] (S3: sf/tbl_df/tbl/data.frame)
## $ Facility ID : num [1:3848] 204 620 1156 2589 3455 ...
## $ Facility Name : chr [1:3848] "Hospice at Lourdes" "Charles T Sitrin Health Care Center Inc" "East Side Nursing Home" "Wellsville Manor Care Center" ...
## $ Short Description : chr [1:3848] "HSPC" "NH" "NH" "NH" ...
## $ Description : chr [1:3848] "Hospice" "Residential Health Care Facility - SNF" "Residential Health Care Facility - SNF" "Residential Health Care Facility - SNF" ...
## $ Facility Open Date : chr [1:3848] "06/01/1985" "02/01/1989" "08/01/1979" "02/01/1989" ...
## $ Facility Address 1 : chr [1:3848] "4102 Old Vestal Road" "2050 Tilden Avenue" "62 Prospect St" "4192A Bolivar Road" ...
## $ Facility Address 2 : chr [1:3848] NA NA NA NA ...
## $ Facility City : chr [1:3848] "Vestal" "New Hartford" "Warsaw" "Wellsville" ...
## $ Facility State : chr [1:3848] "New York" "New York" "New York" "New York" ...
## $ Facility Zip Code : chr [1:3848] "13850" "13413" "14569" "14895" ...
## $ Facility Phone Number : num [1:3848] 6.08e+09 3.16e+09 5.86e+09 5.86e+09 7.17e+09 ...
## $ Facility Fax Number : num [1:3848] NA NA NA NA NA ...
## $ Facility Website : chr [1:3848] NA NA NA NA ...
## $ Facility County Code : num [1:3848] 3 32 60 2 14 ...
## $ Facility County : chr [1:3848] "Broome" "Oneida" "Wyoming" "Allegany" ...
## $ Regional Office ID : num [1:3848] 3 3 1 1 1 7 1 7 5 7 ...
## $ Regional Office : chr [1:3848] "Central New York Regional Office" "Central New York Regional Office" "Western Regional Office - Buffalo" "Western Regional Office - Buffalo" ...
## $ Main Site Name : chr [1:3848] NA NA NA NA ...
## $ Main Site Facility ID : num [1:3848] NA NA NA NA NA ...
## $ Operating Certificate Number: chr [1:3848] "0301501F" "3227304N" "6027303N" "0228305N" ...
## $ Operator Name : chr [1:3848] "Our Lady of Lourdes Memorial Hospital Inc" "Charles T Sitrin Health Care Center, Inc" "East Side Nursing Home Inc" "Wellsville Manor LLC" ...
## $ Operator Address 1 : chr [1:3848] "169 Riverside Drive" "Box 1000 Tilden Avenue" "62 Prospect Street" "4192a Bolivar Road" ...
## $ Operator Address 2 : chr [1:3848] NA NA NA NA ...
## $ Operator City : chr [1:3848] "Binghamton" "New Hartford" "Warsaw" "Wellsville" ...
## $ Operator State : chr [1:3848] "New York" "New York" "New York" "New York" ...
## $ Operator Zip Code : chr [1:3848] "13905" "13413" "14569" "14897" ...
## $ Cooperator Name : chr [1:3848] NA NA NA NA ...
## $ Cooperator Address : chr [1:3848] NA NA NA NA ...
## $ Cooperator Address 2 : chr [1:3848] NA NA NA NA ...
## $ Cooperator City : chr [1:3848] NA NA NA NA ...
## $ Cooperator State : chr [1:3848] "New York" "New York" "New York" "New York" ...
## $ Cooperator Zip Code : chr [1:3848] NA NA NA NA ...
## $ Ownership Type : chr [1:3848] "Not for Profit Corporation" "Not for Profit Corporation" "Business Corporation" "LLC" ...
## $ Facility Location : chr [1:3848] "(42.097095, -75.975243)" "(43.05497, -75.228828)" "(42.738979, -78.12867)" "(42.126461, -77.967834)" ...
## $ geometry :sfc_POINT of length 3848; first list element: 'XY' num [1:2] -76 42.1
## - attr(*, "sf_column")= chr "geometry"
## - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "names")= chr [1:34] "Facility ID" "Facility Name" "Short Description" "Description" ...
Task 4 - Read NYS retail food store data into sf objects
nys_retailfood <- read_csv("nys_retail_food_store_xy.csv",show_col_types = FALSE,
lazy = FALSE,)
# Remove missing values
nys_retailfood <- nys_retailfood %>%
filter(!is.na(X) & !is.na(Y))
# sf object
nys_retailfood_sf <- st_as_sf(nys_retailfood,coords=c("X","Y"),crs = 4326)
Task 5 - Verify in mapview
mapview(nyc_zipcode_sf)
mapview(nys_healthfacilities_sf)
mapview(nys_retailfood_sf, popup = NULL)
Task 6 - Save into RData file
save(nyc_zipcode_sf, nys_healthfacilities_sf, nys_retailfood_sf, file = "spatial_data.RData")