R使用缓冲区从栅格中提取值

时间:2020-05-28 19:18:53

标签: r raster

我正在尝试提取给定纬度的海拔。

  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的列表。谁能帮我解决这个问题?

0 个答案:

没有答案