SHP文件与PostGIS数据库Geom字段：WKT和EWKB相互转换SQL

如果你曾经和 GIS 数据打过交道，一定遇到过两个“老熟人”：一个是磁盘上常见的 SHP 文件，另一个是数据库里神通广大的 PostGIS。SHP 文件是 ESRI 公司的矢量数据格式，而 PostGIS 则是 PostgreSQL 的空间扩展，它让我们的关系型数据库也能处理空间数据。今天，我们就来聊聊它们之间的一座关键桥梁——WKT 和 EWKB 这两种格式的相互转换。### 为什么需要转换？SHP 文件存储几何信息时，使用的是二进制格式（其实内部就是某种编码的二进制），而 PostGIS 在数据库中存储几何字段时，则默认使用一种叫 EWKB（Extended Well-Known Binary）的二进制格式。但我们在查看数据、调试问题、或者写 SQL 查询时，通常希望看到的是人可读的文本格式，比如 WKT（Well-Known Text）。WKT 长这样：POINT(116.39 39.9)，一看就懂。所以，转换的必要性在于：数据库内部用 EWKB 高效存储，而我们人（或者脚本）需要用 WKT 来理解和操作。另外，SHP 文件导入到 PostGIS 时，工具（比如 shp2pgsql）会自动完成格式转换，但如果我们想手动控制，或者从数据库导出数据到 SHP，就需要掌握这些 SQL 转换函数了。### 核心函数：ST_AsEWKT 和 ST_GeomFromEWKT在 PostGIS 中，最常用的两个转换函数就是：- ST_AsEWKT(geometry)：把 EWKB 二进制几何对象转换为 EWKT 文本（带 SRID 的 WKT）。- ST_GeomFromEWKT(text)：把 EWKT 文本解析回几何对象（EWKB）。小提示：WKT 和 EWKT 的区别在于，EWKT 多了一个 SRID=xxxx; 前缀，例如 SRID=4326;POINT(116.39 39.9)。SRID 是空间参考标识符，它告诉数据库这个坐标属于哪个坐标系（比如 4326 就是 WGS84 经纬度坐标系）。#### 代码示例 1：从 PostGIS 查询并输出 WKT假设我们有一张表 land_parcels，里面有个 geom 字段（类型是 geometry）。如果我们想直接查看这条记录的几何信息，可以这样写 SQL：sql-- 查询地块信息，并把 geom 字段转换为可读的 EWKT 文本SELECT parcel_id, ST_AsEWKT(geom) AS geom_wktFROM land_parcelsLIMIT 5;结果可能会是这样：parcel_id | geom_wkt----------+-----------------------------------1 | SRID=4326;POLYGON((116.39 39.9, ...))2 | SRID=4326;POINT(116.41 39.92)你看，这样我们就能直接看到几何形状了。如果只需要纯 WKT（不带 SRID），可以使用 ST_AsText(geom)，它返回的是不带 SRID 的 WKT。### 从 WKT 到 EWKB：插入或更新几何字段反过来，当我们有了一段 WKT 文本（比如从外部系统获取的），想把数据写入 PostGIS 的几何字段时，就需要用到 ST_GeomFromEWKT 或者 ST_GeomFromText（不带 SRID 版本）。注意，如果你用的是 ST_GeomFromText，还需要额外指定 SRID，否则 SRID 会变成 0。#### 代码示例 2：插入一条带 WKT 的数据sql-- 向地块表插入一条新记录，geom 字段使用 EWKT 文本转换而来INSERT INTO land_parcels (parcel_id, geom)VALUES ( 1001, ST_GeomFromEWKT('SRID=4326;POLYGON((116.39 39.9, 116.40 39.9, 116.40 39.91, 116.39 39.91, 116.39 39.9))'));如果你只有普通的 WKT（没 SRID），可以这样：sql-- 使用 ST_GeomFromText 并显式指定 SRID 为 4326INSERT INTO land_parcels (parcel_id, geom)VALUES ( 1002, ST_SetSRID(ST_GeomFromText('POINT(116.42 39.93)'), 4326));这里我用了 ST_SetSRID 来把 SRID 设置为 4326，这样数据才能正确叠加到地图上。### 实战：结合 SHP 文件处理现在我们来做一个更贴近实际的应用场景。假设你有一个 SHP 文件，但你不想用命令行工具导入，而是想自己写一段 Python 脚本读取 SHP 文件，然后生成 SQL 插入到 PostGIS。你可以使用 geopandas 读取 SHP，然后利用 to_wkt() 方法获取 WKT，再拼接 SQL 语句。当然，更推荐直接使用 geopandas.to_postgis()，但为了理解转换原理，我们手动演示一下。pythonimport geopandas as gpd# 读取 SHP 文件gdf = gpd.read_file('my_shapefile.shp')# 确保坐标系是 EPSG:4326（WGS84）gdf = gdf.to_crs(epsg=4326)# 遍历每一行，生成 SQL 插入语句for idx, row in gdf.iterrows(): wkt = row.geometry.wkt # 获取 WKT 文本 # 注意：这里要转义单引号，避免 SQL 注入或语法错误 safe_wkt = wkt.replace("'", "''") sql = f"INSERT INTO target_table (geom) VALUES (ST_GeomFromEWKT('SRID=4326;{safe_wkt}'));" print(sql) # 执行 sql，可以用 psycopg2 连接数据库执行这段 Python 代码展示了如何从 SHP 文件中提取 WKT，并生成对应的插入 SQL。核心就是 ST_GeomFromEWKT 函数。### 进阶：EWKB 的十六进制表示有时候，你会看到 EWKB 的十六进制字符串，例如 0101000020E6100000... 这种。这其实就是 EWKB 的十六进制编码。PostGIS 也支持直接转换：- ST_AsEWKB(geometry) 返回 EWKB 的 bytea 类型。- encode(ST_AsEWKB(geom), 'hex') 可以得到十六进制字符串。- ST_GeomFromEWKB(decode(hex_string, 'hex')) 则从十六进制还原。不过日常使用中，我们更常用 WKT，因为可读性高。EWKB 主要用于程序间的传输。### 总结今天我们聊了 SHP 文件和 PostGIS 之间的核心转换逻辑：WKT 是人类友好的文本，EWKB 是机器高效的二进制。通过 PostGIS 提供的 ST_AsEWKT、ST_GeomFromEWKT 等函数，我们可以轻松地在两者之间切换。在实际项目中，无论是导入 SHP 文件，还是导出数据到其他系统，掌握这些 SQL 转换技巧，能让你的 GIS 开发更加灵活高效。记住：数据库里存的是 EWKB，但你的 SQL 里写的往往是 WKT，二者之间的那座桥，就是这些看似简单却极其强大的函数。希望这篇文章能帮你少走一些弯路！