1

我正在尝试使用 ggplot 根据 @data 插槽中的值映射区域和标签。我的错误与此处、此处和此处的问题相同，但我无法将这些答案翻译为我的问题。我可以从 raster 包中的 getData 绘制美国地图中的州，但我从那里遇到了麻烦。例如，这些是我尝试过的数据集（我列出了包，因为我认为它们可能与问题有关）......

library(rgdal)
library(raster)
library(ggplot2)
library(maps)
library(maptools)
library(ggmap) # used for theme_nothing() later    
us <- getData("GADM",country="USA",level=1)
PA <- us[us$NAME_1 == 'Pennsylvania',]
AK <- us[us$NAME_1 == 'Alaska',]

当我尝试这个时，我得到了一张时髦的地图（PA 是一个多边形）。

ggplot(data=PA, aes(long, lat)) +
 geom_polygon(aes(group=group), color='black', fill=NA) +
 geom_text(data=PA, aes(long, lat, label = PA$NAME_1))

在此处输入图像描述

当我尝试这个时，我得到了美学错误（AK 由 2 个多边形组成）。

ggplot(data=AK, aes(long, lat)) +
 geom_polygon(aes(group=group), color='black', fill=NA) +
 geom_text(data=AK, aes(long, lat, label = AK$NAME_1))

为每个多边形定义的区域为每个多边形定义的区域错误：美学长度必须为 1 或与数据相同 (816071)：x、y、标签

我不知道为什么我必须将 $ 放在 geom_text 中并且已经存在数据调用，但是没有它它无法识别变量。我还尝试按照此代码中的 OBJECTID 为 2 个 AK 多边形着色，但得到相同的美学错误。

ggplot(data=AK, aes(long, lat)) +
 geom_polygon(aes(group=group), color='black', fill=AK$OBJECTID)

感谢您的任何见解。我真的很想更频繁地使用 R 进行映射。

4

3 回答 3

1

正如其他答案和评论中提到的，修改 geom_text 是要进行的第一个更改。我的主要补充是使用coord_map()，它是 ggplot2 的一部分。coord_map()为典型的经纬度映射正确缩放轴。

ggplot(data=AK, aes(long, lat)) +
  geom_polygon(aes(group=group), 
               color = "black", 
               fill = NA, 
               alpha = 0.2) + 
  coord_map()

这导致非常长（因此未显示）但准确的映射。这是由于默认比例对阿拉斯加的极端坐标进行了一些不正确的解释。

通过scale_x_continuous在 -190 和 -120 处对轴进行少量裁剪，可以获得更好的裁剪。

ggplot(data=AK, aes(long, lat)) +
  geom_polygon(aes(group=group), 
               color = "black", 
               fill = "darkgreen", 
               alpha = 0.2) + 
  coord_map()+
  scale_x_continuous(limits = c(-190, -120)) +
  annotate(geom = "text", 
           x = coordinates(AK)[2, 1], # note: [2,] refers to mainland label 
           y = coordinates(AK)[2, 2], 
           label = AK$NAME_1[1],
           size = 10)

阿拉斯加情节 请注意，阿拉斯加的标签坐标需要与宾夕法尼亚选择稍有不同，因为阿拉斯加的 data.frame 包含两个对象。coordinates(AK)[1,]提供最大岛屿的标签坐标，而coordinates(AK)[2,]是大陆的标签坐标。

为了进一步阅读，我正在使用人口普查形状文件绘制人口普查数据，Kevin Johnson 的关于在 R 中制作地图的博客文章非常有帮助。人口普查形状文件具有不同的分辨率级别，因此如果您不知道阿拉斯加轮廓中存在的许多坐标，那么 5m 或 20m 的人口普查形状文件版本可能更易于管理。

于 2016-01-09T06:05:49.440 回答
0

“Pennsylvania”文本条目过多是由于将所有long's 和lat' 值提供给 geom_text 函数造成的。您需要切换到仅采用单个坐标的函数。制作 PA 对象后，我可以做到这一点，诚然这看起来像是一个 hack。S4 对象类应该具有适当的提取器方法，但我在?SpatialPolygonsDataFrame帮助页面和链接中进行了搜索，但没有找到 labpt 插槽的提取器。（然后@aosmith 填写了那个空白，我查找了如何将纵横比设置为 1.0）：

plot(PA)
text(x= PA@polygons[[1]]@Polygons[[1]]@labpt[1],   
     y=PA@polygons[[1]]@Polygons[[1]]@labpt[2], 
     label=PA$NAME_1)

我尝试绘制阿拉斯加对象，但我的 R 会话挂断了。AK 是一个包含所有这些岛屿的相当复杂的空间对象。当您查看结构时，您会看到超过 2 个多边形：

..@ polygons   :List of 2
  .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
  .. .. .. ..@ Polygons :List of 4147

您可能需要遍历 4147 列表以找到面积最大的多边形或输入手动编码的坐标。这个绘图工作的 ggplot 版本PA可能从这个开始（尽管需要固定纵横比）：

ggplot(data=PA, aes(long, lat)) +
 geom_polygon(aes(group=group), color='black', fill=NA) + 
 annotate(geom="text", x=coordinates(PA)[1], y=coordinates(PA)[2], label = PA$NAME_1) + 
 coord_fixed(ratio=1)

在此处输入图像描述

于 2015-12-30T21:10:03.743 回答
0

使用 ggplot 解决我的问题仍然没有多大运气。我现在已经切换到 plot 和 spplot 并且已经找到了如何标记中心的多边形和 Labpt 插槽，如下所示......invisible(text(getSpPPolygonsLabptSlots(PA), labels=as.character(PA$NAME_1), cex=0.5))认为这可能会对某人有所帮助。

于 2016-01-04T20:34:02.177 回答