R ggplot2结合geom_tile和coord_map（“moll”）

我使用标准投影制作了一张地图。我猜这是“mercator”。下面是代码，R ggplot2结合geom_tile和coord_map（“moll”）

p=ggplot(output) 

p=p+geom_tile(aes(x=lon,y=lat,fill=dcm)) 

p=p+scale_fill_gradientn("Depth of DCM (m)",colours=rgb(rgb[1:100,1],rgb[1:100,2],rgb[1:100,3]), 
limits=c(15,200),labels=c(25,50,75,100,125,150,175,200),breaks=c(25,50,75,100,125,150,175,200)) 

p = p+guides(fill = guide_colorbar(barwidth = 0.5, barheight = 9)) 

p=p+layer(data=coastline.world,geom="polygon",mapping=aes(x=longitude,y=latitude)) 



p=p+theme(text=element_text(family="Times",size=9)) 

p=p+theme(legend.title = element_text(face = 'plain')) 

p=p+guides(colour = guide_legend(title.hjust = 0.5))

当我尝试使用

p=p+coord_map("moll")

程序开始运行，但从未停止投影更改为 “mollweide”

。

你有什么想法吗？

感谢亚历克斯

的data.frame coastline.world来自输出的

library(oce) 
data(coastlineWorld) 
coastline.world=data.frame(longitude=coastlineWorld[["longitude"]],latitude=coastlineWorld[["latitude"]])

第一行：

structure(list(lon = c(-180, -179.5, -179, -178.5, -178, -177.5 
), lat = c(-59.5, -59.5, -59.5, -59.5, -59.5, -59.5), dcm = c(NA, 
41.4461206739867, 45.6921865291417, 48.135154847963, 48.4013604947836, 
46.9140989480546)), .Names = c("lon", "lat", "dcm"), row.names = c(NA, 
6L), class = "data.frame")

来源

2013-07-25 user2619599

欢迎来到Stack Overflow！如果你给了我们一个可重现的数据例子（即'output'和'coastline.world'），它会让我们重现问题并找出错误。你可以通过发布'dput（output）'和'dput（coastline.world）'的结果来给出一个可重复的例子，但是如果它们是大数据帧，你可以做'dput（head（coastline.world））'和' dput（head（output））'而不是（只要你确认问题仍然存在） –

任何更新？我正在努力解决类似的问题，我刚开始认为ggplot还没有准备好做出这样的情节。 –

这可能只是数据量太大。 coastline.world有超过400K行。我也改变了你的layer，其中有geom_polygon到geom_points。

所以我尝试了一些统一的抽样（选择每隔100行，占总数据的1％）。

newcw <- coastline.world[seq(1, nrow(coastline.world), by = 100),]

现在，当我试图

p <- p+layer(data=newcw,geom="point",mapping=aes(x=longitude,y=latitude)) 
p <- p+geom_point(data=newcw, mapping=aes(x=longitude,y=latitude))

然后是coord_map，呈现好。

p <- p+coord_map("moll")

这是我得到：

enter image description here

你可能有颜色修补得到你想要的东西。

希望能帮助你前进。

来源

2013-07-25 23:17:42

感谢您的帮助。 coastline.world的大小不是问题，因为我可以使用geom_polygon和coord_map（“moll”）绘制地图。当我尝试用geom_tile绘制输出时，问题就出现了。 – user2619599

R ggplot2结合geom_tile和coord_map（“moll”）

回答

相关问题