I'm working in R and I am new to maps and GIS and I have no idea how it works. But I am trying to return which part of the city specific coordinates are located.
I got the shapefile from here: 

> 
> <https://opendata.arcgis.com/datasets/ac4e22c0d8444b8b81a0ee9af3273f40_2.zip>
> 
> 
> 

I read the file and transform it:

```
library(rgdal)
library(mapview)    
nhood <- readOGR(dsn = "hverfahlutar.shp")
nhood2 = spTransform(nhood, CRS("+proj=longlat +datum=WGS84"))

```

the dataframe I am working with is this:

```
schools <- structure(list(school = skolar <- c("Vesturbæjarskóli", "Landakotsskóli", 
                      "Barnaskóli Hjallastefnunnar í Reykjavík ", "Tjarnarskóli", "Austurbæjarskóli", 
                      "Langholtsskóli", "Vogaskóli", "Suðurhlíðarskóli", "Brúarskóli", 
                      "Hlíðaskóli"), long = c(64.1497778, 64.1478234, 64.1287022, 64.1406399, 
                                              64.1417735, 64.1380849, 64.131696, 64.1216781, 64.1256253, 64.130929
                      ), lat = c(-21.959836, -21.9517121, -21.9311303, -21.9365311, 
                                 -21.9237807, -21.8623288, -21.8606524, -21.9040609, -21.9130514, 
                                 -21.9122735)), .Names = c("school", "long", "lat"), row.names = c(NA, 
                                                                                                   -10L), class = c("tbl_df", "tbl", "data.frame"))

```

I can plot it with points so the data seems to be okay:

```
plot(nhood2)
points(schools$long ~ schools$lat, col = "red", cex = 1)

```

What to do next I do not know. I think I need to use the over function but it has not been working for me.

What I need to do is, for each coordinates in each row of school. I need to return the corresponding value in nhood2$HEITI (the name of the neighbourhood.)

I think I need to prepare the data in the schools dataframe for this, but I don't understand how... or why.