R. ggmap。使用具有不同数据集的多边形覆盖几何密度图

时间:2017-07-23 19:27:22

标签: r ggmap

全部,

我是编程和R的新手,但在这个问题上已经花了很长时间。

我正在尝试根据两个不同的数据集叠加两个地图。第一个是密度图,第二个是几何多边形。我可以分别绘制两个绘图,但是当我将它们绘制在一起时,我得到以下消息:"错误:提供给连续刻度的离散值"。我尝试添加inherit.aes= FALSE因为我怀疑这可能是基于早期帖子的问题,但它没有帮助

US_map=get_map(location="united states", zoom=4)
ggmap(US_map, extent="device")+
  geom_polygon(inherit.aes=FALSE, 
data=fortify(coords),aes(x=coords[,1],y=coords[,2], fill="orange", 
colour="black", alpha=1))+
  annotate("text",x=long_average,y=35,label=share_in_frame_perc,size=4)+
geom_density_2d(inherit.aes=FALSE, data=Mill_dots_T,aes(x=Mill_dots_T[,1], 
y=Mill_dots_T[,2]), size=0.3)+
  stat_density2d(data=Mill_dots_T,
             aes(x=Mill_dots_T[,1], y=Mill_dots_T[,2], 
fill=..level..,alpha=..level..), size=0.01, bins=16, geom="polygon")+
  scale_fill_gradientn(colours = rev(rainbow(7)), breaks=c(0.006, 0.0125, 
0.025, 0.05,0,01,0.02,0.03))+
  scale_alpha(range= c(0.2,1), guide=FALSE)

我的数据看起来有点像这样:

Mill_dots_T
long_rect (first column)
-70
-70
-80
-80
-70
...

lat_rect (second column)
35
45
45
35
35
...

" Mill_dots_T"数据看起来有点像这样

Longitude (first column)
-86,8
-86,8
-86,8
-86,8
-86,3
-86,3
-86,3
-86,3
....
Latitude (second column)
33,5
33,5
33,5
33,5
33,5
32,3
32,3
32,3
32,3
32,3
32,3
32,3
32,3
32,3
....

0 个答案:

没有答案