PostGIS中如何获取与指定LineString相交的多边形数组

嘿，刚接触PostGIS不用慌，这个需求其实挺容易实现的，我来给你一步步拆解清楚～

核心思路 #

咱们要用到PostGIS里专门判断空间相交关系的函数ST_Intersects()，它会返回true只要两个几何对象有任何空间交集——不管是边界接触、部分重叠还是完全穿过，都能检测到。

具体查询语句 #

假设你要传入的LineString是WKT格式（比如'LINESTRING(0 0, 10 10)'这种文本格式），直接用下面的SQL就能搞定：

-- 替换成你自己的表名、几何列名、目标LineString和空间参考系
SELECT *
FROM polygon
WHERE ST_Intersects(polygon.geom, ST_GeomFromText('LINESTRING(0 0, 10 10)', 4326));

关键细节说明： #

• polygon.geom：把这里换成你polygon表中存储多边形几何数据的列名（常见的命名是geom或者geometry）
• ST_GeomFromText()：负责把WKT格式的文本转换成PostGIS能识别的几何对象，第二个参数4326是空间参考系的EPSG编码（也就是常用的WGS84经纬度坐标系，如果你用的是其他坐标系，替换成对应的编码即可）
• 如果你的传入数据是WKB二进制格式，可以用ST_GeomFromWKB()代替ST_GeomFromText()

优化查询性能（可选但强烈推荐） #

如果你的polygon表数据量不小，为了让查询跑得更快，一定要给几何列创建空间索引：

-- 同样替换成你的表名和几何列名
CREATE INDEX idx_polygon_geom ON polygon USING GIST(geom);

空间索引能大幅提升空间关系查询的效率，数据量越大效果越明显。

额外注意事项 #

• 如果你不确定自己的几何列用的是什么空间参考系，可以用这个语句查询：

  SELECT ST_SRID(geom) FROM polygon LIMIT 1;
  

• 必须保证传入的LineString和polygon表几何列的空间参考系一致，否则可能出现错误或者不准确的结果。如果不一致，可以用ST_Transform()来转换坐标系：

  -- 例子：把EPSG3857的LineString转换成EPSG4326，和polygon列匹配
  SELECT *
  FROM polygon
  WHERE ST_Intersects(
    polygon.geom,
    ST_Transform(ST_GeomFromText('LINESTRING(0 0, 10 10)', 3857), 4326)
  );
  

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