File di forma per la conversione raster in R?

Question

Ho uno shapefile scaricato da worldwildlife.org per le ecoregioni terrestri del mondo. Il file può essere caricato qui: http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Viene come un file di forma standard e vorrei fare due cose con esso. Primo: prendere la shapefile dalla mia directory locale e agganciarlo in misura di orientale del Nord America (ext = misura (-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) 

Sono sicuro che devo usare la funzione rasterize nel pacchetto "raster", ma non riesco ancora a farlo funzionare correttamente. Gradirei qualche suggerimento su come farlo.

Answer 1

Avete ragione di pensare che dovreste utilizzare raster (anziché le classi spaziali raster sp) per dati raster spaziali. È inoltre necessario utilizzare rgdal (anziché maptools) per leggere, scrivere e manipolare in altro modo i dati vettoriali spaziali.

Questo dovrebbe iniziare:

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)