安装必要的包（如果尚未安装）
if (!require("sf")) install.packages("sf")
if (!require("rnaturalearth")) install.packages("rnaturalearth")
if (!require("rnaturalearthdata")) install.packages("rnaturalearthdata")
加载包
library(sf)
library(rnaturalearth)
library(rnaturalearthdata)
设置工作路径（替换为你自己的路径）
setwd("F:/XAKS2025110601") # Windows示例
setwd("/path/to/your/directory") # Linux/Mac示例
从rnaturalearth获取全球海岸线数据
coastline <- ne_coastline(scale = "medium", returnclass = "sf")
查看数据基本信息
print("海岸线数据基本信息：")
print(coastline)
查看坐标系
st_crs(coastline)
可视化检查（可选）
plot(st_geometry(coastline), main = "全球海岸线地图")
导出为shapefile
output_file <- "global_coastline.shp"
st_write(coastline, output_file, delete_dsn = TRUE)
检查是否导出成功
if (file.exists(output_file)) {
print(paste("海岸线地图已成功导出到:", getwd(), "/", output_file))
} else {
print("导出失败，请检查错误信息")
}
同时导出相关的投影信息文件（.prj）
print("导出的文件包括：")
list.files(pattern = "global_coastline")
made:yafei 2025