使用ggplot创建具有2维值的等值线(多边形)贴图

时间:2016-05-11 20:35:02

标签: r ggplot2 maps

我想创建一个可视化二维值的等值区域图(=颜色根据有序对(v1, v2)填充多边形,其中v1v2是有序因子)。

以下是结果应如何显示的示例: enter image description here

我认为具有两个维度的颜色矩阵图例清楚地说明了我想要实现的目标。

我想使用ggplot2::geom_polygon实现此目的。

最小例子:

ids <- factor(c("1.1", "2.1", "1.2", "2.2", "1.3", "2.3"))

values <- data.frame(
    id = ids,
    v1 = factor(c("Hi","Med","Med","Hi","Lo","Lo"), 
                levels=c("Lo", "Med", "Hi"), ordered=TRUE),
    v2 = factor(c("Hi","Lo","Lo","Med","Med","Hi"), 
                levels=c("Lo", "Med", "Hi"), ordered=TRUE))

positions <- data.frame(
  id = rep(ids, each = 4),
  x = c(2, 1, 1.1, 2.2, 1, 0, 0.3, 1.1, 2.2, 1.1, 1.2, 2.5, 1.1, 0.3,
        0.5, 1.2, 2.5, 1.2, 1.3, 2.7, 1.2, 0.5, 0.6, 1.3),
  y = c(-0.5, 0, 1, 0.5, 0, 0.5, 1.5, 1, 0.5, 1, 2.1, 1.7, 1, 1.5,
        2.2, 2.1, 1.7, 2.1, 3.2, 2.8, 2.1, 2.2, 3.3, 3.2))

datapoly <- merge(values, positions, by=c("id"))

我想在上面的例子后面的一个中结合以下两个地图。 polgons的颜色填充应该根据有序对(v1, v2),当然我需要一个颜色矩阵图例。

library("ggplot2")
ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=v1, group=id))
ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=v2, group=id))

enter image description here

2 个答案:

答案 0 :(得分:1)

在这里晚点参加聚会,但是,为了以后遇到这个问题的任何人的利益,https://github.com/wmurphyrd/colorplaner似乎确实可以满足您的要求。 usage部分中的第二个示例甚至提供了地图中着色的示例。

答案 1 :(得分:-1)

如何interaction(v1, v2)并在生成的互动中映射单一比例?

# using your inputs
values$v1_2 <- interaction(values$v1, values$v2)
...
datapoly <- merge(values, positions, by=c("id"))

library(ggplot2)
library(scales)
library(RColorBrewer)
ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=v1_2, group=id)) +
  scale_fill_brewer(palette= "Blues")

enter image description here

或者在上面的评论中使用链接到的@gregor解决方案。