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

我尝试以下SQL命令：PostGIS的边框查询返回奇怪的结果

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，但是“食物棒2”和“食物棒3”具有47.32和47.42的ymax值。哪里有问题？

来源

2014-11-04 JustWonder

使用ST_Intersects或ST_Contains来代替&&。 – 2014-11-04 21:54:37

Here is your geography envelope：

url_mapper

这与您的查询点。信封被夷为平地，笛卡尔空间ST_Segmentize：

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

point_in_bbox

那么，你是正确的，因为点应该是地理范围内的，但您使用的&&边框运营商，而忽略了几何形状。对于地理信封上的边框看起来是这样的：

point_in_bbox_w_env

这表明在边框中的所有点。

修复查询像这样的东西：

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

来源

2014-11-04 22:03:09

谢谢迈克。这是一个很好的解释。 – JustWonder 2014-11-04 22:21:25

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

回答

相关问题