Extract values from multiple rasters using lat / long and save as one CSV file

Question

I would like to extract values from multiple rasters with a single lat / long or Shapefile. I have tried every tutorial and stackoverflow post and am getting nowhere. Maybe someone could please help. Here is what I have tried so far:

# Imports
library(raster)
library(sp)

# Raster List
tmp <- lapply(list.files("D:/rasters", pattern = ".tif$", full.names = TRUE), raster)

# Lat and Lon
coords <- data.frame("lat" = c(29.940668), "lon" = c(48.363455))

# Create spatial points
pts <- SpatialPoints(coords = coords, 
              proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0  "))

# Extract the points
as.data.frame(cbind(coords, do.call("rbind", lapply(tmp, extract, pts))))

The code runs but the return values are all NA and I checked and I should be getting back some values.

I have also tried to read the values in as a shapefile using readOGR as in Extracting values from raster according to Lat and long of values?, but I get an error

Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv, : Cannot open layer

I'm at a loss why nothing is working for me.

Answer 1

Thanks all! I finally got it working per the code below.

# Imports
library(raster)
library(sp)

# Raster List
tmp <- lapply(list.files("D:/rasters", pattern = ".tif$", full.names = TRUE), raster)

# Lon and lat without quotes per phloem's comment above  
coords <- data.frame(lon = 47.667, lat = 29.317)

# Create spatial points
pts <- SpatialPoints(coords = coords, 
              proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0  "))

# Extract the points
result <- as.data.frame(cbind(coords, do.call("rbind", lapply(tmp, extract, pts))))

# Write to csv file
write.csv(result, file="filename.csv")

I was also able to finally read in the shapefile too but abandoned that method since the first one finally worked. Nonetheless, here is the code that finally worked to read in the shapefile:

point_location <- readOGR("foldername","shapefile name without extension")

Thanks again!