R:从生成的图像中剪出一个shapefile边界

Question

我制作了一张普通克里格预言的图像。我有一个边界线的形状文件，我想以该形状文件的形式裁剪普通的克里格预测。

这是我用来生成图像的代码:映像(OK.pred、loc=grid、axes=F、useRaster=TRUE)。我只想从图像中剪出一个物体--当我绘制它们时，它们完美地叠加在一起。

它几乎与这里的问题相同，https://gis.stackexchange.com/questions/167170/is-it-possible-to-clip-a-shapefile-to-an-image-in-r，但是我对R还比较陌生，并且完全搞不懂netcdf文件部分。

我找到了一堆关于如何剪辑光栅的代码，但我只是想不出如何将图像保存到变量中，更不用说将其转换为光栅来剪辑它了。任何帮助都将不胜感激！

 OK.pred<-krige.conv(gambling.geo,coords = gambling.geo$coords, data=gambling.geo$data, locations=grid,krige=krige.control(obj.model=gambling.vario.wls))

ordinarykrig = image(OK.pred,loc=grid,axes=F,useRaster=TRUE)

Macau <- readOGR("MAC_adm0.shp")

x <- crop(?...)

Answer 1

摘自http://leg.ufpr.br/geoR/tutorials/kc2sp.R：您需要将克里格输出转换为空间对象，然后才能将其传递给掩码()。下列各方应这样做：

OK.pred<-krige.conv(gambling.geo,coords = gambling.geo$coords, data=gambling.geo$data, locations=grid,krige=krige.control(obj.model=gambling.vario.wls))
GT.s <- points2grid(SpatialPoints(as.matrix(grid)))
reorder <- as.vector(matrix(1:nrow(grid), nc=slot(GT.s, "cells.dim")[2])[,slot(GT.s, "cells.dim")[2]:1])
SGDF.s <- SpatialGridDataFrame(grid=GT.s, data=as.data.frame(OK.pred[1:2])[reorder,])
r<-raster(SGDF.s)
x<-mask(r, Macau)