从多个数据集R中提取坐标答案

【问题标题】：Extract coordinates from multiple datasets R从多个数据集R中提取坐标
【发布时间】：2016-10-04 22:09:53
【问题描述】：

我发现使用从 shapefile 获得的特定地理区域从具有经度和纬度坐标的多个数据帧中提取栅格值时遇到问题。所有地理区域的范围不同，都有对应的坐标数据出现数据框。我尝试使用光栅包中的掩码函数：

首先 - 将 shapefile、栅格图层和数据框加载为列表

shp<-list.files(pattern = "shp")
shape<-lapply(shp,readShapeSpatial,proj4string=CRS("+proj=longlat +ellps=WGS84"))

rastFiles<-list.files(pattern="bil") #List of bioclimatic variables 
layers<-stack(rastFiles)

oco<-list.files(pattern = ".csv")
oco<-lapply(oco,read.table,header=T,sep = ",")
oco<-lapply(oco,"[", c(6,5))
y<-lapply(oco, na.omit)

第二 - 操作 shapefile，裁剪栅格图层并创建我需要的具有地理扩展名的蒙版。

p<- lapply(shape, function(i) {
    i[i@data$ORIGIN==1,] #specifying the region of shapefile
    p1<-SpatialPolygons(i@polygons,proj4string=i@proj4string) #transform specific region in shapefile
    nr<-crop(layers,extent(p1)) #crop raster stack with polygons extent
    m<-mask(nr,p1) #raster mask using poligon as model
    })

因此，对象“p”是一个包含 RasterBrick 对象的列表。为了从每个数据帧中提取单元格值坐标，我尝试了

value<-lapply(p,extract,y)

控制台返回给我：

Error in round(y) : non-numeric argument to mathematical function

检查“p”的结构后，我发现栅格值不可用：

.. ..@ data    :Formal class '.MultipleRasterData' [package "raster"] with 14 slots
.. .. .. ..@ values    : int [1:99660, 1:20] NA NA NA NA NA NA NA NA NA NA ...

但是，我可以使用列表运算符在栅格堆栈中一一提取数据帧值：

extract(nat[[1]],y[[1]])

通过这种方式，我进行了搜索以了解为什么我的值的 lapply 函数不起作用。

感谢您的帮助！

【问题讨论】：

标签： r lapply shapefile r-raster

【解决方案1】：

看看你的最后一行代码，我假设你想在nat和y的第一个元素上应用相同的函数（extract），然后是两者的第二个元素，等等。这是mapply 做了什么。因此，您的代码将如下所示：

mapply(extract, nat, y)

【讨论】：