PostGIS与GEOS缓冲折线结果差异问题排查

核心差异在于计算模型和距离单位的本质不同：

• GEOS是纯平面几何库，它默认将输入坐标的单位当作距离单位。如果你的折线使用WGS84（EPSG:4326）经纬度坐标，GEOS会把「度」作为距离单位计算缓冲——比如你传入100作为缓冲距离，GEOS会理解为100度，这和PostGIS中geography类型的100米缓冲完全不是一个量级。
• PostGIS的geography类型基于椭球面/球面计算，ST_Buffer默认接受的距离单位是米，内部会自动处理球面坐标的曲率，生成符合真实地理空间的缓冲区域。

解决方法 #

方法1：投影转换（推荐，精度高） #

将WGS84经纬度坐标转换为局部米制投影坐标系（比如UTM分带、地方投影），再用GEOS执行缓冲操作，最后转换回WGS84：

1. 确定你的几何所在的UTM分带或合适的米制投影（例如EPSG:32650对应UTM 50N，单位为米）
2. 用坐标转换库（如Perl的Geo::Proj4）将折线从EPSG:4326转换到目标投影
3. 用GEOS对转换后的平面坐标执行缓冲（此时输入的距离单位为米）
4. 将缓冲结果转换回EPSG:4326

方法2：近似单位转换（仅适合赤道附近小范围场景） #

如果不想做投影转换，可采用赤道附近的近似换算：1米≈8.98315×10⁻⁶度。但该方法在高纬度地区误差极大，仅适合小范围、低纬度的简单场景。

验证步骤 #

1. 检查Perl代码中缓冲距离的单位：确认你传入GEOS的数值是否是米（而非直接用PostGIS的米数值，因为GEOS在WGS84下会把它当作度）
2. 确认Geo::Geos是否处理了坐标投影：若库本身不支持投影，需先通过其他工具完成坐标转换再传入GEOS

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