You could use the freely available Level 2 GADM shapefile data for Greenland, which comes in EPSG:4326 (Lat/Lon projection) by default, and import it into R using the GDAL functionalities offered by the rgdal
package. Then, define the geographic extent for Washington Land and crop the previously imported shapefile. Make sure you have all the required dependencies installed (e.g. rgeos
)
# Required packages
library(rgdal)
library(raster)
# Import Greenland shapefile data
shp.greenland <- readOGR(dsn = "/path/to/shapefile/folder",
layer = "GRL_adm2")
# Define cropping extent for Washington Land
ext.washington.land <- extent(-67, -63, 78, 82)
# Subset Greenland shapefile
shp.washington.land <- crop(shp.greenland, ext.washington.land)
Update:
Now, if you'd like to include glaciated areas, you could download e.g. the corresponding 10m shapefile from Natural Earth, import and crop it according to the predefined extent, and add it to your plot.
# Import glaciated areas
shp.glaciers <- readOGR(dsn = "/path/to/shapefile/folder",
layer = "ne_10m_glaciated_areas")
# Subset (and plot) glaciated areas by extent of Washington Land
shp.glaciers <- crop(shp.glaciers, ext.washington.land)
plot(shp.washington.land, col = "darkolivegreen3")
plot(shp.glaciers, col = "azure", add = TRUE)
