Shapefile do konwersji rastrowej w R?

Question

Mam plik kształtu pobrany ze strony worldwildlife.org dla naziemnych ekoregionów świata. Plik można załadować tutaj: http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Jest dostarczany jako standardowy plik kształtu i chciałbym zrobić z nim dwie rzeczy. pierwsze: wziąć shapefile z mojego lokalnego katalogu i przypiąć go do stopnia wschodniej Ameryce Północnej (EXT = stopniu (-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) 

Jestem pewien, że musiał użyć funkcja rasteryzacji w pakiecie "raster", ale nadal nie jestem w stanie sprawić, by działała poprawnie. Byłbym wdzięczny za wszelkie sugestie, jak to zrobić.

Answer 1

Masz rację, myśląc, że powinieneś używać raster (zamiast sp rastrowych klas przestrzennych) dla przestrzennych danych rastrowych. Powinieneś także użyć rgdal (zamiast maptools) do czytania, pisania i manipulowania przestrzennymi danymi wektorowymi.

To powinno Ci zacząć:

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)