Mam plik shapefile pobrany z worldwildlife.org dla ziemskich ekoregionów świata. Plik można załadować tutaj:http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Jest to standardowy plik kształtu i chciałbym zrobić z nim dwie rzeczy. Najpierw: weź plik shapefile z mojego lokalnego katalogu i przypnij go do wschodniej części Ameryki Północnej (ext = range (-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 muszę użyć funkcji rasteryzacji w pakiecie „raster”, ale wciąż nie mogę go poprawnie uruchomić. Byłbym wdzięczny za wszelkie sugestie, jak to zrobić.