我有以下SpatialPolygonsDataFrame。
require(raster)
usa <- getData('GADM', country='USA', level=2)
metro <- subset(usa, NAME_1=="Nebraska" & NAME_2 %in% c("Dodge","Douglas","Sarpy","Washington"))
plot(metro)
我希望能够复制以下地图边界（由颜色定义）：
有谁知道一个好的进攻计划吗？我意识到这是一个有点手动的过程。我已经下载了所有地理更详细的美国人口普查文件。我希望可以汇总更详细的地理水平来回答上述问题，但不幸的是，这些地区的排列不一样。
是否已经有一个 R 函数可以帮助协助这个手动过程？至少，我希望能够利用 4 个县区的周边。
【问题讨论】：
使用rgdal 包中的writeOGR 创建metro 对象的shapefile。然后安装 QGIS (http://www.qgis.org/)，一个免费的开源 GIS，并将 shapefile 作为新图层加载。
然后您可以编辑图层、添加新多边形、编辑线等，然后保存为 shapefile 以读回 R。
此外，您可以“地理参考”您的图像（通过识别图像上已知的经纬度点）并将其作为栅格图层加载到 QGIS 中。这使您的新领域数字化变得更加容易。你所需要的只是一些特定点的经纬度坐标，例如多边形的角或线的交点，然后 QGIS 有一个地理参考插件可以做到这一点。
我认为您不会找到任何像 QGIS 一样适合在图像上数字化新几何图形的 R 代码。
半小时后（和二十年的经验，不是你需要的全部）我得到了这个：
不过，我并没有将您的新边界精确地数字化，只是为了速度。 QGIS 屏幕截图显示了四个都市区下的五个彩色区域。
第一步是地理配准。此屏幕截图显示了如何对 PNG 进行地理参考 - 红线是在通过匹配控制点将 PNG 转换为 GeoTIFF 后，在 PNG 上以透明方式绘制的都市区 shapefile。
第二步是使用 QGIS 编辑工具来分割、连接和创建新的多边形。然后我只是给它们上色并添加标签来美化它。
我可以将这些文件全部打包给你整理，但这真的不需要那么长时间，你会从中学到很多东西。另外，这可能是一个 gis.stackexchange.com 问题...
【讨论】：