Intro Thoughts
Status Quo
library(tidyverse)
nc <- sf::st_read(system.file("shape/nc.shp", package="sf"))
## Reading layer `nc' from data source
## `/Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library/sf/shape/nc.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 100 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -84.32385 ymin: 33.88199 xmax: -75.45698 ymax: 36.58965
## Geodetic CRS: NAD27
nc_ref <- nc |>
select(county_name = NAME, fips = FIPS)
options(sf2stat.ref_data = nc_ref)
getOption("sf2stat.ref_data", nc)
## Simple feature collection with 100 features and 2 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -84.32385 ymin: 33.88199 xmax: -75.45698 ymax: 36.58965
## Geodetic CRS: NAD27
## First 10 features:
## county_name fips geometry
## 1 Ashe 37009 MULTIPOLYGON (((-81.47276 3...
## 2 Alleghany 37005 MULTIPOLYGON (((-81.23989 3...
## 3 Surry 37171 MULTIPOLYGON (((-80.45634 3...
## 4 Currituck 37053 MULTIPOLYGON (((-76.00897 3...
## 5 Northampton 37131 MULTIPOLYGON (((-77.21767 3...
## 6 Hertford 37091 MULTIPOLYGON (((-76.74506 3...
## 7 Camden 37029 MULTIPOLYGON (((-76.00897 3...
## 8 Gates 37073 MULTIPOLYGON (((-76.56251 3...
## 9 Warren 37185 MULTIPOLYGON (((-78.30876 3...
## 10 Stokes 37169 MULTIPOLYGON (((-80.02567 3...
Experiment