У меня есть шейп-файл, загруженный с worldwildlife.org для наземных экорегионов мира. Файл можно загрузить здесь:http://worldwildlife.org/publications/terrestrial-ecoregions-of-the-world.

Он поставляется в виде стандартного файла формы, и я хотел бы сделать с ним две вещи. Во-первых: возьмите шейп-файл из моего локального каталога и обрежьте его в пределах восточной части Северной Америки (доб = экстент (-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)

Я уверен, что мне нужно использовать функцию растеризации в пакете «растр», но я все еще не могу заставить ее работать правильно. Буду признателен за любые предложения о том, как это сделать.