PostGIS中ST_Within函数失效：多边形内坐标却返回未找到

问题描述 #

传入的经纬度坐标确认处于数据库存储的多边形范围内，但使用ST_Within的SQL查询始终返回Location not found，而非预期的Location found，相关代码如下：

@app.get('/polygons/<latitude>/<longitude>')
def verify_polygon(latitude, longitude):
    try:
        conn = connect_db()
        cur = conn.cursor()

        cur.execute(f'SELECT id_0 FROM public."polygons-c3" WHERE ST_Within(ST_SetSRID(ST_MakePoint({longitude}, {latitude}), 4326), geom)')
        result = cur.fetchone()
        cur.close()
        conn.close()

        if result:
            return jsonify({'status': 'Location found', 'lote': result[0]}), 200
        else:
            return jsonify({'status': 'Location not found'}), 404
    except Exception as e:
        return jsonify({'error': str(e)}), 500

排查与解决步骤 #

1. 确认坐标系一致性 #

数据库中geom字段的SRID必须与传入点的SRID（代码中为4326）一致，否则空间判断会失效：

• 执行SQL查询geom的SRID：

  SELECT ST_SRID(geom) FROM public."polygons-c3" LIMIT 1;
  

• 如果返回值不是4326，需将传入点转换为对应坐标系，例如geom为3857（Web墨卡托）时，修改SQL为：

  SELECT id_0 FROM public."polygons-c3" 
  WHERE ST_Within(
    ST_Transform(ST_SetSRID(ST_MakePoint({longitude}, {latitude}), 4326), 3857), 
    geom
  );
  

2. 手动验证点坐标的查询结果 #

排除代码参数传递的问题，直接在数据库中用目标坐标执行查询：

SELECT id_0 FROM public."polygons-c3" 
WHERE ST_Within(ST_SetSRID(ST_MakePoint(你的实际经度值, 你的实际纬度值), 4326), geom);

如果手动查询能返回结果，说明代码中存在参数解析或拼接的问题。

3. 检查多边形拓扑有效性 #

多边形存在自相交、闭合错误等拓扑问题时，ST_Within的判断逻辑会异常：

• 查询无效多边形及错误原因：

  SELECT id_0, ST_IsValidReason(geom) FROM public."polygons-c3" WHERE NOT ST_IsValid(geom);
  

• 修复无效多边形：

  UPDATE public."polygons-c3" SET geom = ST_MakeValid(geom) WHERE NOT ST_IsValid(geom);
  

4. 处理边界点的判断逻辑 #

ST_Within仅判断点是否严格在多边形内部（不包含边界），若点刚好在多边形边界上，会返回false：

• 改用ST_Intersects包含边界判断：

  SELECT id_0 FROM public."polygons-c3" 
  WHERE ST_Intersects(ST_SetSRID(ST_MakePoint({longitude}, {latitude}), 4326), geom);
  

5. 修复SQL注入风险与参数化查询 #

当前代码用字符串格式化拼接SQL，存在注入风险，同时可能因坐标格式问题导致查询异常，改用参数化查询：

cur.execute(
    'SELECT id_0 FROM public."polygons-c3" WHERE ST_Within(ST_SetSRID(ST_MakePoint(%s, %s), 4326), geom)',
    (longitude, latitude)
)

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