Question

我发现使用从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
    })

结果，对象＆＃34; p＆＃34;是一个包含RasterBrick对象的列表。为了从每个数据帧中提取单元格值坐标，我尝试了

value<-lapply(p,extract,y)

控制台告诉我：

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

检查＆＃34; p＆＃34;的结构后我发现栅格值不可用：

.. ..@ 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 ...

但是，我可以使用列表运算符在rasterstack中逐个提取数据帧值：

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

通过这种方式，我搜索了解为什么我的价值观的功能不起作用。

感谢您的帮助！

Answer 1

只是查看最后一行代码，我假设您要在extract和nat的第一个元素上应用相同的函数（y），然后是第二行两者的元素等等。这是mapply的作用。所以，你的代码看起来像这样：

mapply(extract, nat, y)

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

1 个答案: