小编RMF*_*ish的帖子

将 shapefile 转换为光栅

我在栅格化 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本身的问题吗?

请和谢谢!

r shapefile rasterize

3
推荐指数
1
解决办法
1万
查看次数

标签 统计

r ×1

rasterize ×1

shapefile ×1