Я хочу создать карту countour из x, y и z и обрезать данные с помощью шейп-файла, используя ggplot

У меня есть набор данных с Eastings, Northings и Age. Я хочу построить возраст для моего набора данных. Набор данных загружается по следующему адресу: https://www.dropbox.com/s/zmtxluhwuupb9iv/age%20at%202100.csv

Код, который я использовал для создания простого графика и контура, выглядит следующим образом:

age2100 <- read.csv("age at 2100.csv",header=TRUE, sep=",")
p <- ggplot(age2100, aes(x=Eastings, y=Northings, colour=("Age")))+geom_point()
p+geom_density2d()

Вывод, который я получил, выглядит следующим образом:

Я не смог загрузить изображение, поэтому найдите его в папке Dropbox по следующему адресу:

https://www.dropbox.com/s/5k5gefqi7wz6v36/age.png

Теперь я хочу добавить шейп-файл на фон и обрезать контурные линии в шейп-файле.

Желаемый результат выглядит следующим образом:

https://www.dropbox.com/s/tqhx2t8cqqozwbt/7Q10%20flow%20dye%20from%20P%20and%20S.JPG

Как перейти от рисунка 1 к желаемому результату? Если бы кто-нибудь мог дать мне предложения, это было бы здорово.

Большое спасибо.

Дждбаба