使用PostgreSQL进行多边形点的SQL查询

时间:2009-06-18 12:59:09

标签: sql postgresql gis rdbms postgis

我有以下简单的表格:

CREATE TABLE tbl_test
(
  id serial NOT NULL,
  poly polygon NOT NULL
)
WITH (OIDS=FALSE);

然后我尝试插入一个带有多边形的行:

insert into tbl_test values(1, PolyFromText('POLYGON((0 0, 10 10, 10 0, 0 0))'))

并遇到这个错误:

  

列“poly”是多边形类型,但表达式是几何类型

哪个是蹩脚的。所以我的第一个问题是:

  1. 我真的要演员吗?
  2. 无论如何,在投射后它起作用。现在我正在尝试做一个简单的ST_Contains查询:

    select id, poly from tbl_test where ST_Contains(poly, Point(GeomFromText('POINT(9 2)')))
    

    出现错误:

    ERROR:  function st_contains(polygon, point) does not exist
    LINE 1: select id, poly from tbl_test where ST_Contains(poly, Point(...
                                                ^
    HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
    

    我应该做什么?

    以下作品:

    select st_contains(st_geomfromtext('POLYGON((0 0, 10 10, 10 0, 0 0))'), st_geomfromtext('POINT(0 0)'))
    

    但这可能是因为两个参数都是Geometry类型。对表数据的实际查询不起作用。

    答案:

    土井!问题是我创建的DB不是基于postgis模板DB(因此没有相关的函数和几何列表等)。 最后,我可以简单地说,PostGIS要求您向数据库添加数百个函数,行和几个表的方式只是为了让您获得GIS支持是完全蹩脚的。它使模式的备份更加复杂并且非常容易出错(如果您忽略调用AddGeometryColumn并且只是自己添加几何列,则天堂禁止)。

2 个答案:

答案 0 :(得分:9)

多边形是PostGIS基础的基础Postgres类型。使用PostGIS函数select AddGeometryColumn(...)启用几何列。否则你正在使用直的多边形:

=> create table gt (id int, space polygon);
=> insert into gt values (1, '((2,2),(3,4),(3,6),(1,1))');
INSERT 0 1
=> select point(space) from gt where id = 1;
    point    
-------------
 (2.25,3.25)
(1 row)

这是多边形的中心点

=> select circle(space) from gt where id = 1;
             circle             
--------------------------------
 <(2.25,3.25),1.93994028704315>
(1 row)

这是多边形的最小边界圆,表示为Postgres circle类型。记录所有几何运算符here: http://www.postgresql.org/docs/8.3/interactive/functions-geometry.html基本多边形没有任何投影数据,SRID等,因此如果它与PostGIS一起工作,它可能只是默认为预设并且变得幸运。但是,当然有很多情况下你只需要在地理空间范围内的几何形状。

答案 1 :(得分:6)

好的,很奇怪,我发现以下更简单的语法有效:

insert into tbl_test (poly) values ('(0,0),(0,10),(10, 10), (0, 0)')

select * from tbl_test where poly @> '(2, 8)'

但我正在努力弄清楚这些功能和运营商之间的区别。这种较短的语法(不符合OpenGIS标准)是否利用了相同的空间索引等?