在原始Postgres查询中使用PostGIS中的边界框

时间:2020-09-11 22:04:09

标签: sql postgresql gis postgis

我的数据库结构是普通Postgres:

CREATE TABLE IF NOT EXISTS locations (
  name text NOT NULL,
  lat double precision NOT NULL,
  lng double precision NOT NULL,
);
CREATE INDEX ON locations(lat,lng);

当我要计算边界框中左下角和右上角的所有位置时,请使用以下查询:

SELECT * FROM locations 
WHERE lat >= min_lat AND
WHERE lat <= max_lat AND
WHERE lng >= min_lng AND
WHERE lng <= max_lng; 

现在,我想生成一个给定点的边界框,并在locations查询中使用边界框结果。我正在使用以下PostGIS查询来生成边界框:

SELECT 
  ST_Extent(
   ST_Envelope(
    ST_Rotate(
     ST_Buffer(
      ST_GeomFromText('POINT (-87.6297982 41.8781136)',4326)::GEOGRAPHY,160934)::GEOMETRY,0)));

结果:BOX(-89.568160053866 40.4285062983089,-85.6903925527536 43.3273499289221)

但是,我不确定一次调用如何将PostGIS查询边界框中的结果用于原始lat / lng Postgres查询中。关于如何将两者合并的任何想法?最好这样以保留索引。

1 个答案:

答案 0 :(得分:3)

如果要获取bbox坐标作为单独的值,则可以查看ST_XMaxST_YMaxST_XMinST_YMin。嵌入您的查询的以下CTE应该给您一个提示:

WITH j (geom) AS (
 SELECT 
  ST_Extent(ST_Envelope(
   ST_Rotate(ST_Buffer(
    ST_GeomFromText('POINT(-87.6297982 41.8781136)',4326)::GEOGRAPHY,160934)::GEOMETRY,0)))
)
SELECT
    ST_XMax(geom),ST_YMax(geom),
    ST_XMin(geom),ST_YMin(geom)
FROM j


      st_xmax      |     st_ymax     |      st_xmin      |     st_ymin      
-------------------+-----------------+-------------------+------------------
 -85.6903925527536 | 43.327349928921 | -89.5681600538661 | 40.4285062983098

附带说明:将几何值存储为数字可能看起来很简单,但这绝不是更好的选择-特别是在处理多边形时!因此,我真的建议您将这些值存储为geometrygeography,乍一看似乎很复杂,但从长远来看肯定会有所收获。

此答案可能有助于您了解涉及多边形的距离/包含性查询:Getting all Buildings in range of 5 miles from specified coordinates