相关疑难解决方法(0)

我们每天制作地图,显示我们地区30个不同区域的温度计算水平,每个区域根据水平填充不同的颜色.这张地图看起来像

现在我想将地图生成切换为R.我已经下载了省级和市级边界(您可以找到整个西班牙的边界或这里是我所在地区的子集)并设法根据Hadley的示例使用ggplot2绘制它们.

我还可以生成一个包含两列的ascii文件:identifier(CODINE)和每日级别.你可以在这里下载.

这是我的第一个尝试使用R和ggplot2绘制shapefile的脚本,因此可能存在错误并且肯定可以改进,欢迎提出建议.以下代码(基于Hadley之前提到的)对我有用:

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

# Reading municipal boundaries

esp = readOGR(dsn=".", layer="lineas_limite_municipales_etrs89")

muni=subset(esp, esp$PROV1 == "46" | esp$PROV1 == "12" | esp$PROV1 == "3")
muni@data$id = rownames(muni@data)
muni.points = fortify(muni, region="id")
muni.df = join(muni.points, muni@data, by="id")

# Reading province boundaries

prov = readOGR(dsn=".", layer="poligonos_provincia_etrs89")

pr=subset(prov, prov$CODINE == "46" | prov$CODINE == "12" | prov$CODINE == "03" )
pr@data$id = rownames(pr@data)
pr.points = …