9

我从 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 函数，但我仍然无法让它正常工作。我将不胜感激有关如何执行此操作的任何建议。

4

1 回答 1

14

您认为应该将raster（而不是sp栅格空间类）用于空间栅格数据是正确的。您还应该使用rgdal（而不是maptools）来读取、写入和以其他方式操作空间矢量数据。

这应该让你开始：

library(rgdal)
library(raster)

## Read in the ecoregion shapefile (located in R's current working directory)
teow <- readOGR(dsn = "official_teow/official", layer = "wwf_terr_ecos")

## Set up a raster "template" to use in rasterize()
ext <-  extent (-95, -50, 24, 63)
xy <- abs(apply(as.matrix(bbox(ext)), 1, diff))
n <- 5
r <- raster(ext, ncol=xy[1]*n, nrow=xy[2]*n)

## Rasterize the shapefile
rr <-rasterize(teow, r)

## A couple of outputs
writeRaster(rr, "teow.asc")
plot(rr)

在此处输入图像描述

于 2013-02-21T00:25:20.910 回答