我有一个包含推文坐标的数据框,例如:
library(ggplot2)
df <- data.frame(long = c(-58.1, -58.2, -58.3, -58.4, -58.5, -55),
lat = c(-34.1, -34.2, -34.3, -34.4, -34.5, -25))
我想绘制布宜诺斯艾利斯的大都会区,即AMBA。 它由以下区域定义: 纵向:( - 58,-59) 纬度:( - 34,-35)
我的数据框中有一行在AMBA区域之外。 (对问题(2)很重要)
con <- url("http://gadm.org/data/rda/ARG_adm2.RData")
print(load(con))
close(con)
ggmap <- fortify(gadm, region = "NAME_2")
lim <- data.frame(lon = c(-59, -58), lat = c(-35, -34))
ggplot(data=ggmap, aes(x=long, y=lat)) +
scale_x_continuous(limits = c(-59,-58)) +
scale_y_continuous(limits = c(-35,-34)) +
geom_polygon(data = ggmap, fill = "grey80", aes(group=group)) +
geom_path(color="white",aes(group=group)) +
geom_point(data = df, aes(x = lon, y = lat, colour = "red"), alpha = 30/100)
答案 0 :(得分:7)
尝试使用coord_map
代替scale_x_continuous
/ scale_y_continuous
。
设置坐标系的限制将缩放绘图(就像你用放大镜看它一样),并且不会像刻度上的设置限制那样改变基础数据。
ggplot(data=ggmap, aes(x=long, y=lat)) +
geom_polygon(data=ggmap, fill="grey80", aes(group=group)) +
geom_path(color="white",aes(group=group)) +
geom_point(data=df, aes(x=long, y=lat), colour="red", alpha=30/100) +
coord_map(xlim=-c(59, 58), ylim=-c(35,34))