Я пытаюсь определить высоту для заданной долготы.
library(raster)
prj_dd <- "+proj=utm +zone=43 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
elev <- getData('alt', country = 'IND')
crs(elev) <- prj_dd
latlon <- data.frame(lon = 77.125, lat = 8.125)
latlon_sp <- SpatialPoints(latlon, proj4string = CRS(prj_dd))
elev_sub <- raster::extract(elev, latlon_sp)
elev_sub
Я получаю значение NA, которое, как я предполагаю, связано с тем, что долгота лежит на краю / в океане. В таких случаях, как я могу назначить ближайшую сетку с данными о высоте этой долготе. Я пробовал использовать функцию буфера:
elev_sub <- raster::extract(elev, latlon_sp, buffer = 100)
class(elev_sub)
'list'
Она возвращает список, в котором есть NA. Может ли кто-нибудь помочь мне, как go об этом?