PostGIS中的拓扑关系及函数（一）

PostGIS提供了若干拓扑关系函数，本文介绍如下函数：

1. ST_3DIntersects      #相交
2. ST_Contains          #包含
3. ST_ContainsProperly  #完全包含
4. ST_Covers            #覆盖
5. ST_CoveredBy         #被覆盖
6. ST_Crosses           #空间相交

ST_3DIntersects-相交

• 定义
  boolean ST_3DIntersects( geometry geomA , geometry geomB );
• 说明
  重叠或者接触，或者在空间上交叉，都返回True，否则返回false。
• 注意
• 此函数支持3d，不会删除z索引。
• 此函数支持多面体曲面。
• 此函数支持三角形和不规则三角网曲面（TIN）。
• 此方法也由sfgal后端提供。
• 此方法实现SQL/MM规范。
• 示例

SELECT ST_3DIntersects(pt, line), ST_Intersects(pt, line)
	FROM (SELECT 'POINT(0 0 2)'::geometry As pt, 'LINESTRING (0 0 1, 0 2 3)'::geometry As line) As foo;
 st_3dintersects | st_intersects
-----------------+---------------
 f               | t
(1 row)

ST_Contains-包含

• 定义
  boolean ST_Contains(geometry geomA, geometry geomB);
• 说明
  几何A包含几何B当且仅当B的点不在A的外部，且B的内部至少有一点位于A的内部。此定义的一个重要微妙之处在于A不包含其边界，但A包含自身。与之形成对比的是，在ST_Contains中，几何体A本身并不正确包含。
  如果几何图形B完全位于几何图形A内，则返回TRUE。要使此函数有意义，源几何图形必须具有相同的坐标投影，并且具有相同的SRID。ST_Contains是ST_Inside的反面。所以ST_Contains（A，B）意味着ST_Within（B，A）内，除非是无效的几何图形，不管是否定义，结果总是错误的。
• 注意
• 该方法实现了sql1.1的OpenGIS简单特性实现规范。
• 示例

-- 一个圆在另一个圆中
SELECT ST_Contains(smallc, bigc) As smallcontainsbig,
	   ST_Contains(bigc,smallc) As bigcontainssmall,
	   ST_Contains(bigc, ST_Union(smallc, bigc)) as bigcontainsunion,
	   ST_Equals(bigc, ST_Union(smallc, bigc)) as bigisunion,
	   ST_Covers(bigc, ST_ExteriorRing(bigc)) As bigcoversexterior,
	   ST_Contains(bigc, ST_ExteriorRing(bigc)) As bigcontainsexterior
FROM (SELECT ST_Buffer(ST_GeomFromText('POINT(1 2)'), 10) As smallc,
			 ST_Buffer(ST_GeomFromText('POINT(1 2)'), 20) As bigc) As foo;

-- 结果
  smallcontainsbig | bigcontainssmall | bigcontainsunion | bigisunion | bigcoversexterior | bigcontainsexterior
------------------+------------------+------------------+------------+-------------------+---------------------
 f                | t                | t                | t          | t        | f

-- 示例演示包含与完全包含之间的差异
SELECT ST_GeometryType(geomA) As geomtype, ST_Contains(geomA,geomA) AS acontainsa, ST_ContainsProperly(geomA, geomA) AS acontainspropa,
   ST_Contains(geomA, ST_Boundary(geomA)) As acontainsba, ST_ContainsProperly(geomA, ST_Boundary(geomA)) As acontainspropba
FROM (VALUES ( ST_Buffer(ST_Point(1,1), 5,1) ),
			 ( ST_MakeLine(ST_Point(1,1), ST_Point(-1,-1) ) ),
			 ( ST_Point(1,1) )
	  ) As foo(geomA);

  geomtype    | acontainsa | acontainspropa | acontainsba | acontainspropba
--------------+------------+----------------+-------------+-----------------
ST_Polygon    | t          | f              | f           | f
ST_LineString | t          | f              | f           | f
ST_Point      | t          | t              | f           | f

ST_ContainsProperly-完全包含

• 定义
  boolean ST_ContainsProperly(geometry geomA, geometry geomB);
• 说明
  如果B与A的内部相交，但不与边界（或外部）相交，则返回true。
  A不能正确地包含自身，但确实包含自身。
  其他几何体的每个点都是该几何体内部的一个点。
• 注意
• 该方法实现了sql1.1的OpenGIS简单特性实现规范。
• 示例

--一个圆在另一个圆中
	SELECT ST_ContainsProperly(smallc, bigc) As smallcontainspropbig,
	ST_ContainsProperly(bigc,smallc) As bigcontainspropsmall,
	ST_ContainsProperly(bigc, ST_Union(smallc, bigc)) as bigcontainspropunion,
	ST_Equals(bigc, ST_Union(smallc, bigc)) as bigisunion,
	ST_Covers(bigc, ST_ExteriorRing(bigc)) As bigcoversexterior,
	ST_ContainsProperly(bigc, ST_ExteriorRing(bigc)) As bigcontainsexterior
	FROM (SELECT ST_Buffer(ST_GeomFromText('POINT(1 2)'), 10) As smallc,
	ST_Buffer(ST_GeomFromText('POINT(1 2)'), 20) As bigc) As foo;
	--Result
smallcontainspropbig | bigcontainspropsmall | bigcontainspropunion | bigisunion | bigcoversexterior | bigcontainsexterior
------------------+------------------+------------------+------------+-------------------+---------------------
 f                     | t                    | f                    | t          | t                 | f

ST_Covers-覆盖

• 定义

boolean ST_Covers(geometry geomA, geometry geomB);
boolean ST_Covers(geography geogpolyA, geography geogpointB);

• 说明
  如果几何体/地理B中没有点在几何体/地理位置A之外，则返回1（真）。
• 示例

    --几何体对象示例
	--圆覆盖圆
SELECT ST_Covers(smallc,smallc) As smallinsmall,
	ST_Covers(smallc, bigc) As smallcoversbig,
	ST_Covers(bigc, ST_ExteriorRing(bigc)) As bigcoversexterior,
	ST_Contains(bigc, ST_ExteriorRing(bigc)) As bigcontainsexterior
FROM (SELECT ST_Buffer(ST_GeomFromText('POINT(1 2)'), 10) As smallc,
	ST_Buffer(ST_GeomFromText('POINT(1 2)'), 20) As bigc) As foo;
	--结果
 smallinsmall | smallcoversbig | bigcoversexterior | bigcontainsexterior
--------------+----------------+-------------------+---------------------
 t            | f              | t                 | f
(1 row)	

-- 地理位置对象示例
-- 一个具有300米缓冲区的点和另外一个点及和这个点的10米缓冲区的覆盖关系比较
SELECT ST_Covers(geog_poly, geog_pt) As poly_covers_pt,
	ST_Covers(ST_Buffer(geog_pt,10), geog_pt) As buff_10m_covers_cent
	FROM (SELECT ST_Buffer(ST_GeogFromText('SRID=4326;POINT(-99.327 31.4821)'), 300) As geog_poly,
				ST_GeogFromText('SRID=4326;POINT(-99.33 31.483)') As geog_pt ) As foo;

 poly_covers_pt | buff_10m_covers_cent
----------------+------------------
 f              | t

ST_CoveredBy-被覆盖

• 与ST_Coveres参数相反，不再赘述。

ST_Crosses-空间相交

• 定义
  boolean ST_3DIntersects( geometry geomA , geometry geomB );
• 说明
  ST_cross接受两个几何体对象，如果它们的交集“空间交叉”，即几何体具有一些但不是所有内部点，则返回TRUE。几何图形内部的交集不能是空集，并且其尺寸必须小于两个输入几何图形的最大尺寸。此外，两个几何图形的交集不得等于任一源几何图形。否则，返回FALSE。
• 示例
• 以下图示均返回TRUE。
  

-- 考虑一个用户有两个表的情况：一个是道路表和一个高速公路表。
CREATE TABLE roads (
  id serial NOT NULL,
  the_geom geometry,
  CONSTRAINT roads_pkey PRIMARY KEY (road_id)
);
CREATE TABLE highways (
  id serial NOT NULL,
  the_gem geometry,
  CONSTRAINT roads_pkey PRIMARY KEY (road_id)
);
--要确定穿过高速公路的道路列表，请使用类似的查询：
SELECT roads.id
FROM roads, highways
WHERE ST_Crosses(roads.the_geom, highways.the_geom);