Przeanalizowałem zestaw danych punktów GPS za pomocą density.ppp, aby wygenerować rodzaj mapy cieplnej intensywności punktów, jak pokazano poniżej:

Chciałbym jednak, aby obraz był ograniczony do granic pliku kształtu, podobnie jak poniżej:

Pierwszy obraz nazywa się jako

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

gdzie przypadki @ coords są współrzędnymi GPS każdego interesującego punktu, a x to shapefile, który podaje zarys jednostki geograficznej.

Drugi obraz jest wywoływany przy użyciu tego kodu:

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

Czy ktoś wie, jak można to zrobić? Być może nie jest to możliwe z działką () i potrzebuję innego pakietu.

Na marginesie, następnym krokiem, jaki planuję, będzie nałożenie tego obrazu na mapę zaimportowaną z GoogleEarth. Nie jestem jeszcze pewien, jak to zrobić, ale opublikuję odpowiedź, jeśli i kiedy to rozpracuję

Wielkie dzięki