我想要一张地图,以一定比例显示整个区域的价值分布。这是我的数据框:
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
。其原因是这种数据点在地图上更清晰可见。
如何将ggmap
和ggplot
组合为一个绘图的示例:
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)