PostGIS距离计算不准?投影转换如何做?

为什么你算出来的距离“差了十万八千里”？

上周一位在测绘院工作的朋友给我发来一段SQL，一脸困惑：“Dr. Gis，我用ST_Distance算两个点的距离，结果是0.8，单位是度？可实际明明是80多公里啊！”——这几乎是每个刚接触PostGIS的人都会踩的坑。别慌，这不是软件bug，而是你忘了“投影转换”这个关键动作。

地理坐标系（如WGS84）下的“度”，不是长度单位！它代表的是球面上的角度。直接拿经纬度算距离，就像用“橘子皮上的弧度”去量“剥开后平铺的直线距离”，当然不准。

投影的本质：把地球“剥皮摊平”的艺术

想象一下，你要把一个橙子完整地剥开，然后平铺在桌面上——无论你怎么努力，总会有些地方被拉伸或压缩。地图投影就是干这件事的：把三维球面“投影”到二维平面，同时尽量保持形状、面积、方向或距离的准确性。而我们计算距离时，需要的是“等距投影”或“局部保距投影”。

我在参与某市国土空间规划项目时，就吃过这个亏。当时用WGS84直接计算地块边界长度，结果误差高达7%，差点导致方案返工。后来统一转成UTM投影，误差立刻降到0.1%以内。

三步搞定PostGIS投影转换

解决之道非常简单，核心就两步：1）确认当前数据的SRID；2）转换到合适的投影坐标系。下面我手把手教你操作：

第一步：查清你的数据“身份证”——SRID

执行以下SQL，看看你的geometry字段用的是什么坐标系：

SELECT ST_SRID(geom) FROM your_table LIMIT 1;

如果返回4326，说明是WGS84地理坐标系——这正是问题根源。

第二步：选择合适的投影坐标系

中国地区推荐使用：EPSG:3857（Web墨卡托，适合可视化）、EPSG:4547（CGCS2000 / 3-degree Gauss-Kruger zone 37，适合华北地区精确计算）。全球通用可选UTM分区，比如北京在EPSG:32650（北半球第50带）。

第三步：动手转换并计算

使用 ST_Transform 函数进行坐标系转换，再调用 ST_Distance：

-- 假设原数据为WGS84 (SRID=4326)，目标投影为UTM 50N (SRID=32650)
SELECT 
  ST_Distance(
    ST_Transform(point_a, 32650), 
    ST_Transform(point_b, 32650)
  ) AS distance_meters
FROM your_table;

现在，结果单位是米，精度误差通常小于1米，完全满足工程需求。

进阶技巧：创建投影索引提升性能

如果你频繁做距离查询，建议预先创建投影后的几何列并建立空间索引：

-- 添加新列存储投影后几何
ALTER TABLE your_table ADD COLUMN geom_utm geometry(Point, 32650);

-- 批量转换
UPDATE your_table SET geom_utm = ST_Transform(geom, 32650);

-- 创建空间索引
CREATE INDEX idx_your_table_geom_utm ON your_table USING GIST (geom_utm);

这样后续所有距离计算都基于投影坐标，速度和精度双丰收。

总结：距离不准？先问投影！

PostGIS本身没有错，错的是我们忘了“坐标系思维”。记住这个黄金法则：地理坐标用于存储和交换，投影坐标用于计算和分析。下次遇到距离、面积不准，第一反应就该是——检查并转换投影！

你在项目中还遇到过哪些“投影坑”？或者对UTM/高斯克吕格如何选择有疑问？欢迎在评论区留言，我会一一回复。也别忘了分享给身边正在被距离计算折磨的小伙伴！