使用PostGIS识别几何重叠求助：市政多边形与缓冲圆匹配失败

核心问题点及修复步骤： #

1. 经纬度顺序错误
   PostGIS的ST_MakePoint要求参数顺序是经度(longitude)在前，纬度(latitude)在后，你写反了，这会导致生成的点完全偏离目标位置，缓冲圆自然找不到匹配的多边形。

2. 缓冲单位与坐标系不匹配
   如果你的数据是WGS84（EPSG:4326，经纬度格式），直接调用ST_Buffer的话，单位是度，1000度的缓冲范围远远超出实际需求。解决方式二选一：

• 转投影坐标系缓冲：先把经纬度点转成以米为单位的投影（比如UTM），再做缓冲，示例：

  ST_Buffer(ST_Transform(ST_MakePoint(points.longitude, points.latitude), 32633), 1000)
  

  注：32633是UTM 33N的EPSG码，需根据你的数据所在区域选择对应投影。
• 球面地理缓冲：直接基于地理坐标计算米级缓冲，无需转投影，更简便：

  ST_Buffer(ST_SetSRID(ST_MakePoint(points.longitude, points.latitude), 4326)::geography, 1000)::geometry
  

3. 空间算子选择不当
   ST_OVERLAPS只检测两个几何互相重叠且彼此不包含的场景，无法覆盖“圆包含多边形”或“多边形与圆边缘触碰”的情况。改用ST_Intersects，它能捕捉所有空间相交的情形（包含、触碰、重叠都算）。

4. 表连接逻辑优化
   用显式JOIN代替笛卡尔积，逻辑更清晰，也方便后续添加索引优化性能。

修正后的完整查询（推荐球面缓冲方案） #

SELECT 
    points.common_name, 
    municipalities.id, 
    municipalities.name
FROM points
JOIN municipalities ON ST_Intersects(
    -- 生成1000米半径的缓冲圆
    ST_Buffer(ST_SetSRID(ST_MakePoint(points.longitude, points.latitude), 4326)::geography, 1000)::geometry,
    -- 将文本多边形转为带SRID的geometry
    ST_SetSRID(ST_GeomFromText(municipalities.geoPolygon), 4326)
);

额外验证步骤： #

• 检查municipalities.geoPolygon的WKT格式是否合法，比如是否是标准的POLYGON((x1 y1, x2 y2, ..., x1 y1))结构，有没有语法错误。
• 单独测试缓冲圆是否正确：执行SELECT ST_AsText(ST_Buffer(ST_SetSRID(ST_MakePoint(你的目标经度, 目标纬度), 4326)::geography, 1000)::geometry);，查看输出的WKT是否符合预期范围。

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