如何扩展现有Shapefile？求R/ArcMap操作步骤

一、R语言实现（优先方案） #

1. 环境准备与数据读取 #

先安装并加载空间数据处理所需的包，读取你的空间数据：

# 首次运行安装依赖包
install.packages(c("sf", "dplyr", "mapview", "rgeos"))

# 加载包
library(sf)
library(dplyr)
library(mapview)
library(rgeos)

# 读取你的RDS格式空间数据（若为Shapefile，替换为 st_read("你的文件路径.shp")）
myshape <- readRDS('bay.RDS')

# 查看原始数据范围并可视化
print(st_bbox(myshape))
mapview(myshape)  # 弹出交互式地图查看原始区域

2. 坐标投影转换（关键） #

如果你的数据是WGS84地理坐标（EPSG:4326），需先转换为投影坐标（如UTM），否则距离计算会出错。请根据你的区域选择对应EPSG代码（例如中国东部用32650）：

# 转换为投影坐标系
myshape_proj <- st_transform(myshape, crs = 32650)

3. 分方向扩展边界 #

通过方向缓冲区分别向内陆（沿海岸向上）和海域（沿海岸向下）扩展，合并后得到扩大的范围：

# 转为sp对象以使用方向缓冲功能
myshape_sp <- as(myshape_proj, "Spatial")

# 提取海岸线边界线
boundary_sp <- gBoundary(myshape_sp)

# 向内陆扩展（side=1，宽度10000米=10公里，可自行调整）
buffer_land <- gBuffer(boundary_sp, width = 10000, side = 1)
# 向海域扩展（side=2，宽度20000米=20公里，可自行调整）
buffer_sea <- gBuffer(boundary_sp, width = 20000, side = 2)

# 转回sf对象并合并所有区域
buffer_land_sf <- st_as_sf(buffer_land)
buffer_sea_sf <- st_as_sf(buffer_sea)
new_shape <- st_union(myshape_proj, buffer_land_sf, buffer_sea_sf) %>% st_cast("POLYGON")

# 可视化对比原始与扩展后的区域
mapview(new_shape, color = "blue", alpha = 0.5) + mapview(myshape_proj, color = "red", alpha = 0.3)

注：side=1和side=2对应不同的扩展方向，若方向错误，互换两者即可。

4. 输出结果（可选） #

将扩展后的区域保存为Shapefile：

st_write(new_shape, "扩展后的边界.shp", delete_layer = TRUE)

二、ArcMap/ArcGIS Pro操作步骤 #

1. 数据加载与投影转换 #

• 打开软件，添加你的Shapefile到地图。
• 若数据为地理坐标，右键图层 → 【数据】→ 【投影】，选择合适的投影坐标系（如UTM）。

2. 提取海岸线边界 #

• 打开【ArcToolbox】→ 【数据管理工具】→ 【要素】→ 【要素转线】，输入原始Shapefile，输出海岸线边界线。

3. 分方向生成缓冲区 #

• 打开【ArcToolbox】→ 【分析工具】→ 【邻域分析】→ 【缓冲区】：
  • 输入海岸线边界线，设置缓冲距离（如10公里），【侧类型】选【左侧】（对应内陆方向，测试后调整），输出内陆缓冲区。
  • 重复操作，设置海域方向的缓冲距离（如20公里），【侧类型】选【右侧】，输出海域缓冲区。

4. 合并区域 #

• 打开【ArcToolbox】→ 【数据管理工具】→ 【要素】→ 【合并】，将原始Shapefile、内陆缓冲区、海域缓冲区作为输入，输出合并后的新Shapefile。

5. 验证调整 #

在地图中查看结果，若方向或距离不符合需求，调整缓冲参数后重新操作。

内容的提问来源于stack exchange，提问作者Salvador