Read example data
map <- read_sf(system.file("shapes/columbus.shp",
package = "spData"), quiet = TRUE)
map
## Simple feature collection with 49 features and 20 fields
## Geometry type: POLYGON
## Dimension: XY
## Bounding box: xmin: 5.874907 ymin: 10.78863 xmax: 11.28742 ymax: 14.74245
## CRS: NA
## # A tibble: 49 × 21
## AREA PERIMETER COLUMBUS_ COLUMBUS_I POLYID NEIG HOVAL INC CRIME OPEN
## <dbl> <dbl> <dbl> <dbl> <dbl> <int> <dbl> <dbl> <dbl> <dbl>
## 1 0.309 2.44 2 5 1 5 80.5 19.5 15.7 2.85
## 2 0.259 2.24 3 1 2 1 44.6 21.2 18.8 5.30
## 3 0.192 2.19 4 6 3 6 26.4 16.0 30.6 4.53
## 4 0.0838 1.43 5 2 4 2 33.2 4.48 32.4 0.394
## 5 0.489 3.00 6 7 5 7 23.2 11.3 50.7 0.406
## 6 0.283 2.34 7 8 6 8 28.8 16.0 26.1 0.563
## 7 0.257 2.55 8 4 7 4 75 8.44 0.178 0
## 8 0.205 2.14 9 3 8 3 37.1 11.3 38.4 3.48
## 9 0.501 3.17 10 18 9 18 52.6 17.6 30.5 0.527
## 10 0.247 2.09 11 10 10 10 96.4 13.6 34.0 1.55
## # ℹ 39 more rows
## # ℹ 11 more variables: PLUMB <dbl>, DISCBD <dbl>, X <dbl>, Y <dbl>, NSA <dbl>,
## # NSB <dbl>, EW <dbl>, CP <dbl>, THOUS <dbl>, NEIGNO <dbl>,
## # geometry <POLYGON>