在普通 Postgres 查询中使用来自 PostGIS 的边界框
Use bounding box from PostGIS in vanilla Postgres query
我有一个 vanilla 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 查询边界框的结果用于 vanilla lat/lng Postgres 查询。关于如何合并两者的任何想法?最好保留索引。
如果您想将 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
,乍一看这可能看起来很复杂,但在 运行.
上绝对是值得的
这个答案可能会阐明 distance/containment 涉及多边形的查询:
我有一个 vanilla 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 查询边界框的结果用于 vanilla lat/lng Postgres 查询。关于如何合并两者的任何想法?最好保留索引。
如果您想将 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
,乍一看这可能看起来很复杂,但在 运行.
这个答案可能会阐明 distance/containment 涉及多边形的查询: