我在栅格化 shapefile 以在 0.5*0.5 网格上生成点时遇到问题。shapefile 表示全球珊瑚礁对综合威胁的风险级别(低 0、中 100、高 1000、极高 1500)的分类。
我从另一个工作正常的示例中提取了代码,但是当我为我的数据尝试它时,我从 plot 函数中什么也得不到。有关 shapefile 和我的代码的链接,请参见下文:
# Read shapefile into R
library(rgdal)
library(raster)
int.threat.2030 <- readOGR(dsn = "Global_Threats/Integrated_Future",
layer = "rf_int_2030_poly")
## Set up a raster "template" for a 0.5 degree grid
ext <- extent(-110, -50, 0, 35)
gridsize <- 0.5
r <- raster(ext, res=gridsize)
## Rasterize the shapefile
rr <- rasterize(int.threat.2030, r)
## Plot raster
plot(rr)
Run Code Online (Sandbox Code Playgroud)
我可能会出错的任何想法?是shapefile本身的问题吗?
请和谢谢!