Tengo un shapefile descargado de worldwildlife.org para las ecorregiones terrestres del mundo. El archivo se puede cargar aquí:http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Viene como un archivo de forma estándar y me gustaría hacer dos cosas con él. Primero: tome el shapefile de mi directorio local y recórtelo en una extensión del este de América del Norte (ext = extensión (-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)

Estoy seguro de que tengo que usar la función rasterizar en el paquete "raster" pero todavía no puedo hacer que funcione correctamente. Agradecería cualquier sugerencia sobre cómo hacer esto.