在香草 Postgres 查询中使用 PostGIS 中的边界框答案

【问题标题】：Use bounding box from PostGIS in vanilla Postgres query在香草 Postgres 查询中使用 PostGIS 中的边界框
【发布时间】：2020-12-30 10:30:04
【问题描述】：

我有一个普通 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;

现在，我想生成一个给定点的边界框，并在位置查询中使用边界框结果。我正在使用以下 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 查询。关于如何将两者合并的任何想法？最好保留索引。

【问题讨论】：

标签： sql postgresql gis postgis

【解决方案1】：

如果您想将 bbox 坐标作为分隔值获取，您可能需要查看 ST_XMax、ST_YMax、ST_XMin、ST_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

旁注：将几何值存储为数字可能看起来很简单，但这几乎不是更好的选择 - 特别是在处理多边形时！所以我真的建议你将这些值存储为geometry 或geography，乍一看可能看起来很复杂，但从长远来看肯定会有所回报。

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

【讨论】：