2

我有 3 个 shp 文件用于创建城市地图：

  • land.shp（绘制水面以上土地的多边形）
  • road.shp（绘制所有道路的多边形，注意其中一些是“圆形道路”，即中间有一个洞）
  • building.shp（多边形绘制所有建筑物）

我使用 QGIS 绘制了我想要的地图，然后我使用 ggplot 播放land.shp，然后使用 road.shp 和建筑 shp 再做一次。下面是谷歌地图的输出来说明我的问题：

在此处输入图像描述

你可以看到它们之间有 2 座桥梁和一些海（我没有海 shp，我只是将背景设置为蓝色），用蓝点标记。在 R 中，该区域应该是一个洞，但它全是灰色的。同样的问题是用红点标记的灰色区域，这是一块土地，另一个用绿点标记的灰色区域，是被道路包围的建筑物。

我将在 road.shp 的洞中拥有陆地/海洋/建筑物，我无法使用 R 显示它们。

谁能教我如何在 R 中显示 road.shp 层背后的东西？谢谢。

4

1 回答 1

11

绘制带孔多边形的一种常见约定是：

  • 具有逆时针方向前进的点的封闭多边形形成实心形状
  • 点顺时针前进的封闭多边形形成一个洞

所以，让我们构建一些数据并绘制：

library(ggplot2)

ids <- letters[1:2]

# IDs and values to use for fill colour
values <- data.frame(
  id = ids,
  value = c(4,5)
)

# Polygon position
positions <- data.frame(
  id = rep(ids, each = 10),
  #     shape      hole         shape        hole
  x = c(1,4,4,1,1, 2,2,3,3,2,   5,10,10,5,5, 6,6,9,9,6),
  y = c(1,1,4,4,1, 2,3,3,2,2,   5,5,10,10,5, 6,9,9,6,6)
)

# Merge positions and values
datapoly <- merge(values, positions, by=c("id"))

ggplot(datapoly, aes(x=x, y=y)) + 
  geom_polygon(aes(group=id, fill=factor(value))) +
  scale_fill_discrete("Key")

在此处输入图像描述

于 2012-08-20T07:40:00.270 回答