Question

我正在尝试学习如何用ggplot2绘制shapefile，我发现了一个wiki：

https://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefiles

我下载并解压缩了我桌面上的新目录（http://www.epa.gov/nheerl/arm/documents/design_doc/ecoregion_design.zip）。

我开始阅读这篇文章，一行一行，但我遇到了fortify的问题。

require("rgdal")
require("maptools")
require("ggplot2")
require("plyr")

> my_dsn <- "directory here"
> 
>   utah = readOGR(dsn=my_dsn, layer="eco_l3_ut")
OGR data source with driver: ESRI Shapefile 
Source: "directory here", layer: "eco_l3_ut"
with 10 features and 7 fields
Feature type: wkbPolygon with 2 dimensions
>   utah@data$id = rownames(utah@data)
>   utah.points = fortify(utah, region="id")
Error: isTRUE(gpclibPermitStatus()) is not TRUE

我查看了?fortify帮助页面并获得了如何使用它的清单。我认为它被正确使用。我还在SO上找到了另一个带有相同错误的问题，

Getting error while using fortify function in R (ggplot)

但我无法弄清楚如何应用我的解决方案所说的内容。如果我现在可以访问Dropbox，我可能会帮助，但我没有。

Answer 1

简短的回答是安装rgeos包并将其包含在顶部。

可以在此处（以及其他地方）找到更长的答案：https://stat.ethz.ch/pipermail/r-sig-geo/2012-July/015685.html

通过在ggplot wiki中绘制shapefile的错误

1 个答案: