使用R的writeOGR向MySQL导入多面边界失败，求解决方案

我明白你在用GDAL的writeOGR导入空间多边形数据到MySQL时碰到了图层创建失败的问题，咱们先试着排查下当前方法的问题，再看看其他可行的导入路径。

先排查当前writeOGR失败的可能原因 #

在尝试替代方法之前，先确认是不是基础环节出了问题：

• 验证MySQL连接权限：先确保你的数据库账号有足够权限（创建表、写入数据），可以用RMySQL先测试连接是否正常：

library(RMySQL)
con <- dbConnect(MySQL(), dbname = "testdb", username = "testuser", password = "****", host = "127.0.0.1")
# 能列出表说明连接没问题
dbListTables(con)
dbDisconnect(con)

• 简化驱动参数：先去掉layer_options='ENGINE=Aria'试试，有时候特定存储引擎的兼容性会引发问题，用默认引擎先测试：

writeOGR(data_2013, dsn=DSN, driver='MySQL', layer='testinglayer', overwrite_layer = TRUE)

• 查看详细错误日志：GDAL的R接口有时候不会输出完整错误，你可以用命令行工具ogrinfo或ogr2ogr测试，能得到更具体的报错信息：

# 先测试连接数据库
ogrinfo MySQL:testdb,user=testuser,host=127.0.0.1,password=**** -so -al
# 尝试用命令行导入
ogr2ogr -f MySQL MySQL:testdb,user=testuser,host=127.0.0.1,password=**** /Path/to/shapefile/testinglayer.shp -nln testinglayer -overwrite

• 检查GDAL依赖：确保你的GDAL安装了完整的MySQL驱动依赖，比如Linux下需要libmysqlclient-dev，Windows下需要对应的MySQL客户端库。

替代导入方法 #

如果上述排查没解决问题，试试下面几种更可靠的方法：

方法1：用sf包直接写入（推荐） #

sf现在是R空间数据处理的主流工具，对MySQL的空间支持更完善：

1. 把你的SpatialPolygonsDataFrame转换成sf对象：

library(sf)
data_2013_sf <- st_as_sf(data_2013)

2. 连接MySQL并写入数据：

library(RMySQL)
con <- dbConnect(MySQL(), dbname = "testdb", username = "testuser", password = "****", host = "127.0.0.1")
# 直接写入，支持覆盖现有表
st_write(data_2013_sf, con, "testinglayer", overwrite = TRUE)
# 如果需要指定存储引擎，添加layer_options参数
# st_write(data_2013_sf, con, "testinglayer", overwrite = TRUE, layer_options = "ENGINE=Aria")
dbDisconnect(con)

方法2：先导出GeoJSON，再用MySQL原生函数导入 #

这种方法适合需要更精细控制表结构的场景：

1. 把R中的空间数据导出为GeoJSON：

st_write(data_2013_sf, "/Path/to/your_data.geojson")

2. 在MySQL中手动创建带空间字段的表（根据你的数据属性调整字段）：

CREATE TABLE testinglayer (
  id INT AUTO_INCREMENT PRIMARY KEY,
  -- 替换成你数据里的实际属性列
  name VARCHAR(255),
  area DOUBLE,
  geom GEOMETRY
);

3. 用MySQL的JSON函数批量导入数据：

LOAD DATA INFILE '/Path/to/your_data.geojson'
INTO TABLE testinglayer
FIELDS TERMINATED BY '\n'
(@json)
SET 
  name = JSON_UNQUOTE(JSON_EXTRACT(@json, '$.properties.name')),
  area = JSON_EXTRACT(@json, '$.properties.area'),
  geom = ST_GeomFromGeoJSON(JSON_EXTRACT(@json, '$.geometry'));

注意：要确保MySQL的secure_file_priv配置允许读取该GeoJSON文件的路径，或者把文件放到MySQL允许的目录下。

方法3：用dbplyr结合sf导入 #

如果你习惯用dplyr风格的代码，可以用dbplyr来完成导入：

library(dbplyr)
library(sf)
library(RMySQL)

data_2013_sf <- st_as_sf(data_2013)
con <- dbConnect(MySQL(), dbname = "testdb", username = "testuser", password = "****", host = "127.0.0.1")

# 用copy_to把sf对象写入数据库
copy_to(con, data_2013_sf, "testinglayer", overwrite = TRUE)
dbDisconnect(con)

copy_to会自动处理空间字段的类型映射，操作起来比较简洁。

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