如何在使用ggplot2创建的等值区域地图中添加标签?

时间:2015-03-10 11:24:02

标签: r ggplot2 gis

我正在尝试并且当前未能将文本注释添加到我在ggplot2中创建的等值区域地图中。我试图用它的名字标记每个多边形(一个地方政府区域)。

在我继续之前,我知道在(非常好的)教程here中已经提出了类似的问题。但是,我尝试了几种不成功的方法,并认为我可能偶然发现了一个不同的问题。我怀疑我的代码失败的原因之一是我试图注释geom_polygon(),而我见过的其他方法详细说明了如何注释geom_map对象。也就是说,我无法想到为什么用geom_polygon不可能做到这一点。

我在下面提供了我的代码。 You can download my data from here。数据框包括我加入强化shapefile的数据。我尝试追加的标签位于" LGA_NAME11"列中。

## LOAD PACKAGES
require(ggplot2)
require(rgdal)
require(dplyr)

## SET GGPLOT THEME
theme_clean <- function(base_size = 12) {
    require(grid)
    theme_grey(base_size) %+replace%
            theme(
                    axis.title = element_blank(),
                    axis.text = element_blank(),
                    panel.background = element_blank(),
                    panel.grid = element_blank(),
                    axis.ticks.length = unit(0,"cm"), 
                    axis.ticks.margin = unit(0,"cm"),
                    panel.margin = unit(0,"lines"),
                    plot.margin = unit(c(0, 0, 0, 0), "lines"),
                    complete = TRUE
            )}

## SET COLOUR PALETTES
palette1 <- c("#f2f0f7", "#dadaeb", "#bcbddc", "#9e9ac8", "#756bb1", "#54278f")

## SET LABEL NAMES
lgaNamesSydney <- aggregate(cbind(long, lat) ~ LGA_NAME11, data=sydneyMapData, FUN = function(x) mean(range(x)))
lgaNamesSydney <- lgaNamesSydney %>% rename(lga = LGA_NAME11)
lgaNamesSydney$angle <- 0

## ATTEMPT TO PLOT MAP WITH LABELS
ggplot(sydneyMapData) +
    aes(long, lat, group=group, fill=Factor1) +
    geom_polygon() +
    geom_text(data=lgaNamesSydney, aes(long, lat, label = LGA_NAME11, angle=angle, map_id =NULL), size=2.5) +
    scale_fill_manual(values = palette1) +
    labs(fill="Drop Bears \nper 1000 population") +
    coord_map(projection = "mercator") +
    theme_clean()

如果有人有任何建议,如果有人能说出我错的地方,我将非常感激。提前致谢。

顺便说一句,我认识到标记等级可以减损视觉体验,但是老板特别要求他们!

Here is my plot sans labels

1 个答案:

答案 0 :(得分:5)

你可以添加像这样的标签

# http://stackoverflow.com/questions/9441778/improve-centering-county-names-ggplot-maps
centroids <- setNames(do.call("rbind.data.frame", by(sydneyMapData, sydneyMapData$group, function(x) {Polygon(x[c('long', 'lat')])@labpt})), c('long', 'lat')) 
centroids$label <- sydneyMapData$LGA_NAME11[match(rownames(centroids), sydneyMapData$group)]

ggplot(sydneyMapData, aes(long, lat, group=group, fill=Factor1)) +
  geom_polygon(colour = "white") +
  with(centroids, annotate(geom="text", x = long, y=lat, label = label, size = 2.5)) 

enter image description here