1

我有一个列表（从 Stata dta 文件导入），其中包含为地图定义空间多边形的条目，以提供一个想法：

> typeof(aux3)
[1] "list"
> mode(aux3)
[1] "list"
> head(aux3)
  _ID        _X       _Y
1   1        NA       NA
2   1 -22.23933 64.56315
3   1 -22.25667 64.56053
4   1 -22.25026 64.56653
5   1 -22.27167 64.57319
6   1 -22.30311 64.55409
> tail(head(aux3,20000))
      _ID       _X       _Y
19995   2 21.21593 60.24696
19996   2 21.21650 60.24337
19997   2 21.23972 60.24913
19998   2 21.22203 60.23304
19999   2 21.21618 60.23332
20000   2 21.22092 60.23930
> #etc.

我想把它转换成一个数据结构，我可以用它在 R 中创建一个地图（没有更大的难度，也不需要广泛的专业知识，以前从未做过）。我推断 sp 包的 SpacialPolygons 类型是最简单的选择。此外，从定义来看，似乎 SpacialPolygons 方法（在此包中定义，请参阅包的文档，第 79 页）似乎是从列表转换为此数据类型的正确方法。

不幸的是，该方法并不那么简单，我需要一些帮助。我的（天真）尝试产生了一个我不理解的错误，并且在谷歌搜索中没有出现任何有趣的结果：

> library(maptools)
Loading required package: sp
Checking rgeos availability: FALSE
    Note: when rgeos is not available, polygon geometry     computations in maptools depend on gpclib,
    which has a restricted licence. It is disabled by default;
    to enable gpclib, type gpclibPermit()
> SP<-SpatialPolygons(aux3)
Error in SpatialPolygons(aux3) : 
  cannot get a slot ("area") from an object of type "integer"

上面的列表可以转换为 SpacialPolygon 吗？如果是这样，怎么做？如果不是，我应该选择什么格式？谢谢。

4

1 回答 1

1

如果你只想绘制地图，你可以使用 ggplot。我无法测试这个，因为我没有你的数据，但应该可以。

library(ggplot2)
aux_df <- data.frame(aux3)
# as per comment: the list-to-dataframe conversion prepends the column names with X
ggplot(aux_df, aes(x=X_X, y=X_Y, group=X_ID) + 
  geom_polygon(color='black', fill=NA) +
  coord_map() +
  theme_classic()

以上假设列表在 和 中aux3具有经度和纬度，_X以及_Y由 定义的多边形_ID。它默认应用墨卡托投影，但您可以使用不同的参数将其更改为coord_map().

这是一个使用俄亥俄县地图的工作示例（任意选择）。

library(ggplot2)
ohio_poly <- map_data('county', region='ohio')
ggplot(ohio_poly, aes(x=long, y=lat, group=group)) +
  geom_polygon(color='black', fill=NA) +
  coord_map() +
  theme_classic()

在此处输入图像描述

于 2015-08-10T15:38:46.757 回答