postgresql – PostGIS边界框查询返回奇怪的结果

我尝试了以下SQL命令:

CREATE TABLE places(
    lat_lng geography(Point,4326),place_name varchar(50)
);

CREATE INDEX places_lat_lng_idx ON places USING gist(lat_lng);

INSERT INTO places values ('POINT(-126.4 45.32)','Food Bar1');
INSERT INTO places values ('POINT(-126.4 47.32)','Food Bar2');
INSERT INTO places values ('POINT(-125.4 47.42)','Food Bar3');

SELECT place_name,ST_AsText(lat_lng) as point FROM places WHERE places.lat_lng && ST_MakeEnvelope(-130.0,44.0,-100.0,46.7,4326);

结果是:

place_name |        point        
------------+---------------------
 Food Bar1  | POINT(-126.4 45.32)
 Food Bar2  | POINT(-126.4 47.32)
 Food Bar3  | POINT(-125.4 47.42)

这对我来说不合适,因为ymax是46.7,但“Food Bar2”和“Food Bar3”的ymax值分别为47.32和47.42.问题出在哪儿?

解决方法

Here is your geography envelope

在这里,您的查询点.使用ST_Segmentize将信封展平为笛卡尔空间:

SELECT ST_Segmentize(ST_MakeEnvelope(-130.0,4326)::geography,50000);

所以你是正确的,这些点应该在地理范围内,但你使用了&&边界框操作符,忽略几何形状.地理信封的边界框如下所示:

它显示了边界框中的所有点.

使用以下内容修复查询:

SELECT place_name,ST_AsText(lat_lng) as point
FROM places
WHERE ST_Intersects(ST_MakeEnvelope(-130.0,places.lat_lng)

相关文章

文章浏览阅读601次。Oracle的数据导入导出是一项基本的技能,...
文章浏览阅读553次。开头还是介绍一下群,如果感兴趣polardb...
文章浏览阅读3.5k次,点赞3次,收藏7次。折腾了两个小时多才...
文章浏览阅读2.7k次。JSON 代表 JavaScript Object Notation...
文章浏览阅读2.9k次,点赞2次,收藏6次。navicat 连接postgr...
文章浏览阅读1.4k次。postgre进阶sql,包含分组排序、JSON解...