10

我正在尝试使用 R 中的函数 fortify 来绘制地图。但是我在处理数据时总是出错。我使用以下代码读取了形状文件：

shape <- readShapeSpatial("Lond_City/lond_city.shp")

可以在以下保管箱位置找到 shapefile： https ://www.dropbox.com/sh/d4w6saailydtu1r/5oIa56xV6S

然后我尝试使用以下命令：

shape1 <- fortify(shape,region="data")

在上面的代码中，我应该用什么来代替“数据”。这是一个看起来令人困惑的地方。我可以在 R 中打开数据或 shapefile。但是，当我运行 fortify 行时，我收到以下错误：

'[.data.frame'(attr, , region) 中的错误：选择了未定义的列

谁能告诉我如何通过示例正确强化形状文件？我想使用 fortify 的原因是我想将形状文件与数据点结合起来。

太感谢了。

贾巴巴

4

2 回答 2

12

您的问题是您需要提供fortify一个实际存在的列。如果您str(lon.shape)从下面的示例代码中键入，您将看到没有region列：

> str(lon.shape)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
  ..@ data       :'data.frame': 1 obs. of  7 variables:
  .. ..$ ons_label: chr "00AA"
  .. ..$ name     : chr "City of London"
  .. ..$ label    : chr "02AA"
  .. ..$ X_max    : num 533843
  .. ..$ y_max    : num 182198
  .. ..$ X_min    : num 0
  .. ..$ y_min    : num 0
  ..@ polygons   :List of 1
  .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
  .. .. .. ..@ Polygons :List of 1
  .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
  .. .. .. .. .. .. ..@ labpt  : num [1:2] 532464 181220
  .. .. .. .. .. .. ..@ area   : num 3151465
  .. .. .. .. .. .. ..@ hole   : logi FALSE
  .. .. .. .. .. .. ..@ ringDir: int 1
  .. .. .. .. .. .. ..@ coords : num [1:771, 1:2] 531027 531029 531036 531074 531107 ...
  .. .. .. ..@ plotOrder: int 1
  .. .. .. ..@ labpt    : num [1:2] 532464 181220
  .. .. .. ..@ ID       : chr "0"
  .. .. .. ..@ area     : num 3151465
  ..@ plotOrder  : int 1
  ..@ bbox       : num [1:2, 1:2] 530967 180404 533843 182198
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "x" "y"
  .. .. ..$ : chr [1:2] "min" "max"
  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
  .. .. ..@ projargs: chr "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs"

相反，请尝试使用name，如下所示：

lon.df <- fortify(lon.shape, region = "name")

产生这个数据框：

> head(lon.df)
      long      lat order  hole piece            group             id
1 531026.9 181611.1     1 FALSE     1 City of London.1 City of London
2 531028.5 181611.2     2 FALSE     1 City of London.1 City of London
3 531036.1 181611.5     3 FALSE     1 City of London.1 City of London
4 531074.0 181610.3     4 FALSE     1 City of London.1 City of London
5 531107.0 181609.3     5 FALSE     1 City of London.1 City of London
6 531117.1 181608.9     6 FALSE     1 City of London.1 City of London

这是一种解决方法，从开始到结束：

library(rgdal)
library(ggplot2)
library(rgeos)

shape.dir <- "c:/test/london" # use your directory name here

lon.shape <- readOGR(shape.dir, layer = "lond_city")
lon.df <- fortify(lon.shape, region = "name")

ggplot(lon.df, aes(x = long, y = lat, group = group)) +
    geom_polygon(colour = "black", fill = "grey80", size = 1) +
    theme()

...导致以下结果：

截屏

于 2013-01-17T03:01:25.543 回答
2

你应该label把data.

require(maptools); require(ggplot2)
shape <- readShapeSpatial("Lond_City/lond_city.shp")
shape1 <- fortify(shape,region="label")

为了找出这一点，我检查了以下数据元素shape：

> shape@data
  ons_label           name label    X_max    y_max X_min y_min
0      00AA City of London  02AA 533842.7 182198.4     0     0

这表明ons_label或name也会起作用，并且可能更适合您的目的。

请注意，您的形状文件有点不寻常，因为似乎只有一组。

于 2013-01-17T03:00:32.480 回答