Eu tenho um shapefile baixado do worldwildlife.org para as ecorregiões terrestres do mundo. O arquivo pode ser carregado aqui:http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Ele vem como um arquivo de formato padrão e eu gostaria de fazer duas coisas com ele. Primeiro: pegue o shapefile do meu diretório local e prenda-o ao leste da América do Norte (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)

Tenho certeza que tenho que usar a função rasterize no pacote "raster", mas ainda não consigo fazer funcionar corretamente. Eu gostaria de receber sugestões sobre como fazer isso.