Я проанализировал набор данных точек GPS, используя плотность .ppp, чтобы сгенерировать тепловую карту интенсивности точек, как показано ниже:

Однако я бы хотел, чтобы изображение было ограничено границами шейп-файла, как показано ниже:

Первое изображение называется

x <- readShapePoly("dk.shp")
xlim<-c(min(912),max(920))
ylim<-c(min(8023),max(8030))
a<-ppp(cases@coords[,1], cases@coords[,2], xlim, ylim, unitname=c("km"))
plot(density.ppp(a, 0.1), col=COLORS)
plot(x, add=T, border="white")

где случаях @ координаты - это координаты GPS каждой интересующей точки, а х - шейп-файл, который обеспечивает схему географической единицы.

Второе изображение вызывается с помощью этого кода:

plot(x, axes=T, col=COLORS, border="White")

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

Кроме того, следующим шагом, который я планирую сделать, будет наложение этого изображения на карту, импортированную из GoogleEarth. Я еще не уверен, как это сделать, но опубликую ответ, если и когда я его решу

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