如何使用rgdal绘制Shapefile时仅显示轮廓而不显示内部边界？

当然可以实现！要去掉行政区域的内部边界，核心就是把同属一个大区域的子多边形合并成单一多边形，这样绘图时就只会渲染外部边界了。下面给你两种实用的解决方案：

方法1：基于属性字段合并（适合有明确层级的行政数据） #

如果你的Shapefile里有用来区分大区域的属性字段（比如省级ID、城市名称这类，同一大区域下的子区域对应相同值），可以用sp包的unionSpatialPolygons来合并：

library(rgdal)
library(sp)

# 读取Shapefile，替换成你的文件路径和图层名
shp <- readOGR(dsn = "./shapefile_folder", layer = "admin_boundaries")

# 假设用来分组的字段是"province_code"，按这个字段合并子多边形
merged_shp <- unionSpatialPolygons(shp, IDs = shp@data$province_code)

# 绘图，此时内部边界就消失了
plot(merged_shp, col = "lightblue", border = "darkgray")

方法2：自动合并所有相邻/重叠多边形（无属性字段也能用） #

如果没有明确的分组字段，或者想直接合并所有相连的多边形，可以用rgeos包的gUnaryUnion工具，它会自动识别并合并相邻的多边形：

library(rgdal)
library(rgeos)

shp <- readOGR(dsn = "./shapefile_folder", layer = "admin_boundaries")

# 合并所有多边形（独立的大区域会被保留为单独的多边形）
merged_shp <- gUnaryUnion(shp)

# 绘图
plot(merged_shp, col = "lightgreen", border = "black")

小提示 #

如果合并时遇到拓扑错误（比如多边形重叠、缝隙导致合并失败），可以先修复拓扑问题再合并：

# 修复拓扑
shp_fixed <- gBuffer(shp, byid = TRUE, width = 0)
# 再执行合并操作
merged_shp <- gUnaryUnion(shp_fixed)

内容的提问来源于stack exchange，提问作者René Nyffenegger