R에서 쉐이프 파일 - 래스터 변환?

Question

전 세계의 지구 생태 지역에 대한 worldwildlife.org에서 다운로드 한 shapefile이 있습니다. 파일은 여기에로드 할 수 있습니다 : http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

표준 모양 파일로 제공되며 두 가지를하고 싶습니다. 첫째 : 내 로컬 디렉토리에서 Shape 파일을 가지고 내가 사용해야하는 확신 동부 북미 (내선 = 범위 (-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) 

의 정도에 클립 rasterize 함수는 "래스터"패키지에 있지만 여전히 제대로 작동하지 않습니다. 이 작업을 수행하는 방법에 대한 제안은 감사하겠습니다.

Answer 1

공간 래스터 데이터의 경우 raster (sp 래스터 Spatial 클래스가 아닌)을 사용해야한다고 생각하는 것이 좋습니다. 또한 공간 벡터 데이터를 읽거나 쓰고 조작하는 데 maptools이 아닌 rgdal을 사용해야합니다.

이것은 당신이 시작할 수 있어야합니다

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)