He analizado un conjunto de datos de puntos GPS usando density.ppp para generar una especie de mapa de calor de intensidad de los puntos, como se muestra a continuación:

Sin embargo, me gustaría que la imagen se limite a los bordes del shapefile, similar a la siguiente:

La primera imagen se llama como

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")

donde los casos @coordos son las coordenadas GPS de cada punto de interés, y x es un shapefile que proporciona el esquema de la unidad geográfica.

La segunda imagen se llama utilizando este código:

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

¿Alguien sabe cómo se puede hacer esto? Quizás no sea posible con plot () y necesitaré otro paquete.

Además, el próximo paso que planeo hacer es superponer esta imagen sobre un mapa importado de GoogleEarth. Tampoco estoy seguro de cómo hacerlo, pero publicaré la respuesta cuando lo resuelva

muchas gracias