PostGIS常用函数有哪些?ST_Buffer怎么写?

别再手动画缓冲区了！PostGIS函数实战指南

你是不是还在用桌面GIS软件一个个点选要素、设置半径、点击“确定”来生成缓冲区？在处理成千上万条道路或地块时，这种操作不仅耗时，还容易出错。更惨的是——当你老板临时要求把缓冲区从50米改成100米时，你得重头再来一遍。别慌，今天Dr. Gis带你用PostGIS的ST_Buffer函数，一行SQL解决战斗。

PostGIS常用函数速览：你的空间分析瑞士军刀

PostGIS不是魔法，但它确实像一把多功能瑞士军刀——小巧但功能强大。以下是我最常使用的几类函数，覆盖90%的日常需求：

• 几何构造函数：ST_GeomFromText, ST_Point, ST_MakeLine —— 把文本或坐标转成几何对象。
• 空间关系判断：ST_Intersects, ST_Contains, ST_DWithin —— 判断两个要素是否相交、包含或在一定距离内。
• 几何操作函数：ST_Buffer, ST_Union, ST_Difference —— 生成缓冲区、合并、裁剪几何体。
• 测量与属性：ST_Area, ST_Length, ST_Centroid —— 计算面积、长度、质心。

我在参与某市国土空间规划项目时，曾用ST_DWithin快速筛选出距离地铁站500米内的所有待开发地块，原本需要半天的手动操作，现在30秒搞定。这就是空间数据库的力量。

ST_Buffer详解：给你的数据“穿件外套”

想象一下，你要给一条河流加个“保护带”，或者给一个工厂画个“影响范围圈”。这其实就是缓冲区（Buffer）——以某个要素为中心，向外扩展指定距离形成的区域。在PostGIS里，ST_Buffer就是干这个的。

它的基本语法如下：

ST_Buffer(geometry geomA, float radius);

举个实际例子：假设你有一张表roads，里面存着城市道路的几何数据（字段名为geom），你想为每条路生成50米的缓冲区：

SELECT road_name, ST_Buffer(geom, 50) AS buffer_geom FROM roads;

就这么简单。但别急，这里有几个坑我踩过，你得注意：

1. 单位问题：如果你的数据是WGS84（经纬度），那50代表的是“50度”，不是50米！这时候你需要先投影到平面坐标系，比如UTM，或者使用ST_Transform转换后再缓冲。
2. 负值缓冲：对多边形使用负半径，可以实现“向内收缩”，常用于生成核心区或排除边缘干扰。
3. 端点样式：默认是圆头（endcap=round），你也可以指定'endcap=square'或'endcap=flat'来改变缓冲区两端形状。

进阶写法示例（带投影和样式控制）：

SELECT 
  road_name,
  ST_Buffer(
    ST_Transform(geom, 32650), -- 转换到UTM Zone 50N
    50,                        -- 50米
    'endcap=flat join=mitre'   -- 平头+斜角连接
  ) AS buffer_geom
FROM roads;

实战案例：找出学校周边200米内的餐饮店

我们来串烧一下前面学的函数。假设有两张表：schools（学校）和restaurants（餐馆），都包含geom字段。

目标：找出所有距离学校200米内的餐馆，并标注是哪所学校附近的。

SELECT 
  r.name AS restaurant_name,
  s.name AS school_name,
  ST_Distance(r.geom, s.geom) AS distance_meters
FROM restaurants r, schools s
WHERE ST_DWithin(r.geom, s.geom, 200)
ORDER BY s.name, distance_meters;

这里用了ST_DWithin做空间筛选（比先缓冲再相交更快），再用ST_Distance计算精确距离。效率高，语义清晰。

总结与行动号召

今天我们拆解了PostGIS中最实用的几类函数，重点攻克了ST_Buffer的写法和常见陷阱。记住：空间分析不是点鼠标，而是写逻辑。掌握这些函数，你就掌握了自动化处理海量空间数据的钥匙。

你在项目中用过哪些PostGIS神函数？有没有被ST_Buffer的单位坑过？欢迎在评论区分享你的血泪史或神操作，我们一起升级打怪！