Binary logical operations (11)

counties_with_point <- nc[unlist(inter), ]
counties_with_point$NAME

##  [1] "Cabarrus"    "Beaufort"    "Davie"       "Clay"        "Northampton"
##  [6] "Pender"      "Warren"      "Halifax"     "Robeson"     "Surry"

cbind(nc_points, areaname = counties_with_point$NAME)

## Simple feature collection with 10 features and 1 field
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -83.74344 ymin: 34.62668 xmax: -76.97503 ymax: 36.5447
## Geodetic CRS:  NAD27
##       areaname                          x
## 1     Cabarrus POINT (-80.69259 35.47426)
## 2     Beaufort POINT (-76.97503 35.30125)
## 3        Davie  POINT (-80.61781 35.8667)
## 4         Clay POINT (-83.74344 35.01048)
## 5  Northampton  POINT (-77.42125 36.5447)
## 6       Pender POINT (-77.70961 34.68814)
## 7       Warren  POINT (-78.1314 36.28169)
## 8      Halifax POINT (-77.67933 36.36195)
## 9      Robeson POINT (-79.41414 34.62668)
## 10       Surry POINT (-80.67994 36.40195)