將點 (Lat-Lon) 添加到 Marmap 圖中

Question

我想繪制一系列海豹軌跡的經緯度點，每個點根據屬性着色，繪制到顯示水深（100m 等高線）和海岸線的地圖上。 我學會了如何使用 marmap 和 ggplot2 創建地圖來顯示測深+海岸線。 代碼在這里：

dat <- getNOAA.bathy(-58,-62.5,43,46.0,res=0, keep=TRUE)
plot(dat,image=TRUE,bpal = list(c(min(dat), 0, "darkblue", "blue","lightblue"), c(0, max(dat), "gray90","gray10")),drawlabels=TRUE,deep=c(-500,200,0),shallow=c(-500,100,0),step=c(500,100,0),lwd=c(1,1,1),lty=c(1,1,1),land=TRUE)+
  scaleBathy(dat, deg=1.232, x="bottomleft", inset=5) #100km

這創建了一個有用的地圖。 但是，我對如何在這張地圖上添加海豹軌跡感到猶豫。

我可以在 ggmap 中做到這一點（使用下面的代碼），但我更喜歡 marmap 地圖

myLocation <- c(-62.5,43,-58,46)
seal_map2<-get_map(location=myLocation,maptype="watercolor",source="stamen",zoom=10)
ggmap(seal_map2)+
geom_point(data=sealtrack,aes(color=category),size=0.5)+
scale_color_gradientn(colours=rainbow(6), breaks=seq(1,6,by=1)) 

任何指導將不勝感激

Answer 1

在“強化”它之后，您應該能夠將來自 marmap 的測深信息作為等高線圖層添加到您的圖上。 沒有您的數據，很難確保它正常工作（並且 NOAA 服務器現在對我來說已關閉）：

library(ggplot2)
library(marmap)
dat <- getNOAA.bathy(-58,-62.5,43,46.0,res=0, keep=TRUE)
dat <- fortify(dat)
ggmap(seal_map2) +
   geom_contour(dat, aes(x = x, y = y, z = z)) +
   geom_point(data=sealtrack,aes(color=category),size=0.5) +
   scale_color_gradientn(colours=rainbow(6), breaks=seq(1,6,by=1))