我正在尝试提取给定纬度的海拔。
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值,我认为这是因为lon lat位于边缘/海洋中 在这种情况下,如何将带有高程数据的最近网格分配给该坐标。我尝试使用缓冲功能:
elev_sub <- raster::extract(elev, latlon_sp, buffer = 100)
class(elev_sub)
'list'
它返回具有NA的列表。谁能帮我解决这个问题?