PostGIS导入Shapefile?常用命令是什么?

为什么你用 shp2pgsql 导入 Shapefile 总报错？

“表不存在”、“编码不对”、“几何类型不匹配”——这三大高频报错，我在国土空间规划项目里帮实习生 debug 过不下二十次。别慌，这不是你操作失误，而是 PostGIS 的“脾气”没摸透。今天我就手把手带你把 Shapefile 安稳送进 PostgreSQL 的怀抱。

Shapefile 导入的本质：不是复制，是翻译

想象一下，你要把一本繁体竖排的古籍（Shapefile）录入到现代数据库（PostGIS）里。你不能直接扫描粘贴，得先“转码”——这就是 shp2pgsql 干的事。它负责三件事：

1. 读取 .shp、.dbf、.prj 文件；
2. 按目标数据库结构生成 SQL 插入语句；
3. 执行导入或输出 SQL 脚本供你审核。

Dr. Gis 小贴士：我第一次给某市做智慧交通底图时，就因忽略 .prj 文件导致全市路网偏移了300米——投影没对上，等于地图穿了错袜子。

实战命令详解：从新手到老鸟都该收藏

最基础命令长这样：

shp2pgsql -s 4326 roads.shp public.roads | psql -U postgres -d gisdb

拆解来看：

• -s 4326：指定源数据坐标系（WGS84），不写这个，90%概率出错！
• roads.shp：你的 Shapefile 文件名（无需带路径时默认当前目录）。
• public.roads：目标模式.表名，不写模式默认 public。
• | psql ...：管道符直连数据库执行，省去中间文件。

但真实项目哪有这么理想？下面这些“进阶开关”才是保命符：

避坑指南：三个你绝对想不到的雷区

雷区一：字段名大小写敏感
Shapefile 的字段在 .dbf 里可能是大写，但 PostgreSQL 默认转小写。如果你代码里写 SELECT RoadName，实际表里是 roadname ——直接报错。解决方案：用双引号包裹或统一转小写。

雷区二：多部件几何(MultiPart)被拆散
一个湖泊包含多个岛屿（MultiPolygon），导入后可能变成多个独立 Polygon。加参数 -S 可强制保持单部件，但更推荐用 ST_Collect 后处理聚合。

雷区三：空表不报错，你以为成功了
如果 .shp 文件损坏或路径错误，命令可能静默执行且返回“成功”，但表里一条记录都没有。务必加 -v 参数开启 verbose 模式看详细日志。

终极技巧：先生成 SQL 再人工审核

生产环境别直接管道导入！先输出 SQL 文件人工检查：

shp2pgsql -s 4326 -I -W UTF-8 roads.shp roads > import.sql
psql -U postgres -d gisdb -f import.sql

这样你能打开 import.sql 确认建表语句、坐标系、字段映射是否正确——相当于给数据进口装了个“安检门”。

总结：导入 Shapefile 的黄金 checklist

下次操作前，默念这五步：

1. ✅ 确认 .prj 存在且坐标系明确（用 -s 指定）
2. ✅ 中文字段？加 -W UTF-8
3. ✅ 首次导入用 -c，追加用 -a
4. ✅ 大数据量必加 -I 建索引
5. ✅ 生产环境先输出 SQL 审核，再执行

掌握这些，你导入的成功率能从 60% 飙到 99%。你在导入时踩过什么奇葩的坑？或者有更骚的操作？评论区等你来Battle！