我想要一张地图,以一定比例显示整个区域的价值分布。这是我的数据框:

head(trip)
   TIP      LON     LAT
1 1.318878 -73.9932 40.7223
2 1.370667 -73.9932 40.7222
3 0.933232 -73.9772 40.7559
4 1.268699 -73.9932 40.7628
5 1.265304 -73.9932 40.7429
6 1.193437 -73.9852 40.7447


我使用以下代码生成了地图:

map <- ggmap::get_map("new york", zoom = 11)
nyc_map <- ggmap::ggmap(map, legend="topleft")


该地图已正确显示。

我还用数据表示生成了一个图层:

ggplot(aes(LON,LAT, fill = TIP), data=as.data.frame(trip)) +
geom_tile() +
scale_fill_continuous(low="white", high="blue") +
coord_equal()


并且这也会生成并显示。

问题是当我想这样做时:

nyc_map +
ggplot(aes(LON,LAT, fill = TIP), data=as.data.frame(trip)) +
geom_tile() +
scale_fill_continuous(low="white", high="blue") +
coord_equal()


我收到以下错误:

Error in p + o : non-numeric argument to binary operator
In addition: Warning message:
Incompatible methods ("+.gg", "Ops.data.frame") for "+"


如果您能帮助我加入这两个对象,将不胜感激。

最佳答案

查看您的数据,最好使用geom_point代替geom_tile。其原因是这种数据点在地图上更清晰可见。

如何将ggmapggplot组合为一个绘图的示例:

library(ggplot2)
library(ggmap)

nyc_map <- get_map("new york", zoom = 12, maptype = "hybrid")

ggmap(nyc_map) +
  geom_point(data=trip, aes(x=LON, y=LAT, fill=TIP), size=6, shape=21, position="jitter") +
  scale_fill_continuous(low="white", high="blue")


这段示例代码给出了以下映射:



在上面的代码中,我使用了position="jitter",因为其中两个点是重叠的。

使用的数据:

trip <- read.table(text="TIP      LON     LAT
1.318878 -73.9932 40.7223
1.370667 -73.9932 40.7222
0.933232 -73.9772 40.7559
1.268699 -73.9932 40.7628
1.265304 -73.9932 40.7429
1.193437 -73.9852 40.7447", header=TRUE)

07-24 13:58