R 将具有等高线的二维密度图覆盖到地图上

所以我现在正在研究如何将等高线叠加到地图上，到目前为止我已经有了这个。任何帮助都将不胜感激。提前谢谢

worldmap = getMap(resolution = "high")
NrthAm = worldmap[which(worldmap$REGION =="North America"),]
plot(NrthAm, xlim=c(-123.35,-122.65),
     ylim=c(49,49.35), main = "Pokemon in Vancouver")

newnew = "CREATE VIEW location AS SELECT latitude, longitude FROM Vanpoke"
query1 = "SELECT latitude, longitude from location"
point = dbGetQuery(dbcon,query1)

kde = kde2d(point$longitude,point$latitude)
image(kde)
contour(kde)

编辑1：我添加了一些示例数据

   latitude longitude
1 49.33597 -123.0458
2 49.27268 -122.8834
3 49.25722 -122.9652
4 49.30699 -123.0204
5 49.19668 -122.9262
6 49.33281 -123.0636

---

Edit2：我在地图上添加了一个密度图，但我不确定它是否正确。

除非您想让我们访问您的数据库（强烈反对），否则您需要向我们提供地图上显示的

query1

的样本。请告知您使用的非基本软件包，可能存在歧义。这个代码是不是很好，这个问题还是很好的，但是考虑一下如何使它具有可重复性，这样就更容易让人拿起它并尝试它（没有访问你的私人数据和数据库）。这包括示例代码（特别是在这里，列出非基本R包）、不在包中的示例明确数据（例如，

dput（head（x））

或

data.frame（x=…，y=…）

）和预期输出。参考文献：，还有。谢谢你的回复，我对这一切都有点陌生。所以我使用RSQLite、rworldmap和rworldetra来解决这个问题。对于查询1，它是一个带有经度和纬度的表。我希望它像笛卡尔坐标系一样显示在地图上，x=经度，y=纬度。请理解“这是一个带有经度和纬度的表”是一回事，但我们不能将该描述用于猜测以外的任何东西（是

runif（2,30,35）

和

runif（2100102）

足够好吗？如果是，请使用它）。如果你认为我们应该考虑一些数据来帮助你解决这个问题，你的责任就是给我们一些东西。我不在乎你使用的是

RSQLite

，你可以简化你的问题，只给出示例数据。另外，您似乎正在使用

质量

，对吗？（我不得不追查这类信息，这让我有点泄气。）我对信息的缺乏感到非常抱歉，因为我还在想办法。我在帖子中添加了一些示例数据和一个绘图。感谢您的时间。除非您想让我们访问您的数据库（强烈反对），否则您需要从地图上显示的

query1

中向我们提供您想要的样本。请告知您使用的非基本软件包，可能存在歧义。这个代码是不是很好，这个问题还是很好的，但是考虑一下如何使它具有可重复性，这样就更容易让人拿起它并尝试它（没有访问你的私人数据和数据库）。这包括示例代码（特别是在这里，列出非基本R包）、不在包中的示例明确数据（例如，

dput（head（x））

或

data.frame（x=…，y=…）

）和预期输出。参考文献：，还有。谢谢你的回复，我对这一切都有点陌生。所以我使用RSQLite、rworldmap和rworldetra来解决这个问题。对于查询1，它是一个带有经度和纬度的表。我希望它像笛卡尔坐标系一样显示在地图上，x=经度，y=纬度。请理解“这是一个带有经度和纬度的表”是一回事，但我们不能将该描述用于猜测以外的任何东西（是

runif（2,30,35）

和

runif（2100102）

足够好吗？如果是，请使用它）。如果你认为我们应该考虑一些数据来帮助你解决这个问题，你的责任就是给我们一些东西。我不在乎你使用的是

RSQLite

，你可以简化你的问题，只给出示例数据。另外，您似乎正在使用

质量

，对吗？（我不得不追查这类信息，这让我有点泄气。）我对信息的缺乏感到非常抱歉，因为我还在想办法。我在帖子中添加了一些示例数据和一个绘图。谢谢你抽出时间。