我从 worldwildlife.org 下载了世界陆地生态区的 shapefile。该文件可以在这里加载:http ://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world 。
它是一个标准的形状文件,我想用它做两件事。首先:从我的本地目录中获取 shapefile 并将其剪辑到北美东部的范围内(ext= extent (-95, -50, 24, 63))
# Read shapefile using package "maptools"
eco_shp <- readShapeLines("F:/01_2013/Ecoregions/Global/wwf_terr_ecos.shp",
proj4string=CRS("+proj=utm +zone=33 +datum=WGS84"))
# Set the desired extent for the final raster using package "raster"
ext <- extent(-95, -50, 24, 63)
我确信我必须使用“raster”包中的 rasterize 函数,但我仍然无法让它正常工作。我将不胜感激有关如何执行此操作的任何建议。