更改绘制geom_polygon，ggplot的顺序

Question

我正在使用R / ggplot创建一个世界地震区的全球地图。我在这里使用修改后的shapefile：http://gmo.gfz-potsdam.de/。

图中有9层，超过4400组。但是，有些图层应该有孔，但似乎没有。因此，为了解决这个问题，我需要首先渲染这些多边形，使它们出现在后面。使用geom_polygon时如何更改多边形的渲染顺序？

require("rgdal") # requires sp, will use proj.4 if installed
require("maptools")
require("ggplot2")
gpclibPermit() # required for fortify method
require("plyr")

fill_colours <- c("#dddddd", "#8FC463FF", "#46A742FF", "#FFFC00FF", "#F2B212FF", "#E9726AFF", "#D92C50FF", "#FF1B23FF", "#B1681FFF")

eq.shp = readOGR(dsn=".", "PGA_PolygonsPY v2_region")

# Create a discrete scale, there are probably more efficient ways of doing this
eq.shp$MAXVALUE[eq.shp$MAXVALUE==10.0] = 9
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 4.8] = 8
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 4.0] = 7
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 3.2] = 6
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 2.4] = 5
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 1.6] = 4
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.8] = 3
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.4] = 2
eq.shp$MAXVALUE[eq.shp$MAXVALUE== 0.2] = 1
eq.shp$MAXVALUE  <- sprintf("%d", eq.shp$MAXVALUE)

eq.shp@data$id = rownames(eq.shp@data)
eq.points = fortify(eq.shp, region='id')

# Join the data frames back together
eq.df = join(eq.points, eq.shp@data, by='id')

p <- ggplot()
p <- p + geom_polygon(data = eq.df, aes(long, lat, group=group, fill=MAXVALUE))
p <- p + scale_fill_manual(values=fill_colours) 
p <- p + coord_equal(ratio = 1, xlim=c(65, 85), ylim=c(35, 45))
p    

此代码生成以下图表

如果将上述图表与原始数据数据（使用GIS软件绘制）进行比较，则可以看到以下工件：   - 原始图中点（Long = 77，Lat = 42）周围的区域不存在。这是因为棕色多边形最后呈现。   - 同样对于周围的大绿色区域（长= 80，纬度= 38），它全部为黄色。这是因为黄色多边形在绿色多边形之后呈现。   - 我不太关心漏洞，例如附近区域（长= 74，纬度= 39）。

Answer 1

好的，我已经能够解决这个问题了。虽然这是一种解决方法，并且渲染时间增加了，但它确实符合我的要求。

我没有尝试使用ggplot对多边形的渲染进行排序/重新排序，而是通过使用for循环强制执行此操作，并根据需要渲染图层。

解决方案：

替换以下代码：

p <- p + geom_polygon(data = eq.df, aes(long, lat, group=group, fill=MAXVALUE))

以下内容：

for (i in seq(9, 1, -1)) {
    p <- p + geom_polygon(data = eq.df[eq.df$MAXVALUE==sprintf("%d", i),], aes(long, lat, group=group, fill=MAXVALUE))
}

这会导致多边形按照我希望的顺序渲染。然后我得到了与原版非常相似的东西：

但是，我可以看到我引入了新的问题，有一些没有孔的多边形，所以渲染它们会导致它们覆盖其他一些层。因此，我需要在多边形中创建孔以正确渲染它们。

感谢所有发表评论的人。

<强>更新 上面的解决方案可以重新排序多边形的渲染，但它并不能完全解决潜在的问题。您似乎需要更复杂的算法才能这样做。我也尝试按多边形区域排序，然后按多边形的加权区域排序，但都没有达到理想的解决方案。最终解决方案需要一些人工干预。