library(raster)
library(rgbif)
library(rgdal)

# Download data from gadm.org
adm <- getData("GADM", country = "FI", level = 0)
# Project to YKJ
adm.YKJ <- spTransform(adm, CRS("+init=epsg:2393"))

# Get some example data
key <- name_backbone(name = "Boloria selene", kingdom = "animals")$speciesKey
# Limit to 200 data points
selene.occ <- occ_search(taxonKey = key, return = "data", country = "FI", limit = 200)
# Crate a SpatialPointsDataFrame for WGS84 -> YKJ transformation
sp.selene <- SpatialPointsDataFrame(coords = selene.occ[, c("longitude", "latitude")], 
    data = data.frame(spp = rep("Boloria selene", 200)), proj4string = CRS("+proj=longlat +datum=WGS84"))
sp.selene.YKJ <- spTransform(sp.selene, CRS("+init=epsg:2393"))

# Fill
plot(adm.YKJ, lwd = 1, border = "grey", col = "grey")
points(sp.selene.YKJ, pch = 21, col = "black", bg = "red")
title(main = "Boloria selene")

plot of chunk plot


# Or just the outlines
plot(adm.YKJ, lwd = 1, border = "darkgrey")
points(sp.selene.YKJ, pch = 21, col = "black", bg = "red")
title(main = "Boloria selene")

plot of chunk plot