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>