AnalyticDB:ジオメトリモデル

ジオメトリの測定

ST_Length および ST_Area を使用して、距離や面積を計算します。これらの関数は、ジオメトリの座標参照系に基づく単位（SRID 3857 の場合はメートル）で結果を返します。

-- 道路の総延長（キロメートル単位）
SELECT SUM(ST_Length(the_geom)) / 1000 AS km_roads FROM bc_roads;

     km_roads
------------------
 70842.1243039643
(1 row)

-- 特定の市町村の面積（ヘクタール単位）
SELECT ST_Area(the_geom) / 10000 AS hectares
FROM bc_municipality
WHERE name = 'PRINCE GEORGE';

      hectares
------------------
 32657.9103824927
(1 row)

空間関係のテスト

GanosBase Geometry は PostGIS のすべての空間関係関数をサポートしています。以下の表に、最も一般的な関数を示します。

例 — ST_Contains を使用して市町村ごとの道路延長を算出：

SELECT m.name, SUM(ST_Length(r.the_geom)) / 1000 AS roads_km
FROM bc_roads AS r, bc_municipality AS m
WHERE ST_Contains(m.the_geom, r.the_geom)
GROUP BY m.name
ORDER BY roads_km;

例 — ST_Covers を使用して 2 つの円のカバー範囲を比較：

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)

ST_Covers と ST_Contains は、境界点の扱い方が異なります。ST_Covers は A が B の境界を覆う場合に true を返しますが、ST_Contains はそうしません。

例 — ST_Disjoint を使用してポイントとラインの非共有点をテスト：

-- ライン上にないポイント：非共有点
SELECT ST_Disjoint('POINT(0 0)'::geometry, 'LINESTRING ( 2 0, 0 2 )'::geometry);
-- 結果：t

-- ライン上のポイント：非共有点ではない
SELECT ST_Disjoint('POINT(0 0)'::geometry, 'LINESTRING ( 0 0, 0 2 )'::geometry);
-- 結果：f

例 — ST_Overlaps を使用して重複と交差をテスト：

SELECT ST_Overlaps(a, b) AS a_overlap_b,
       ST_Crosses(a, b)  AS a_crosses_b,
       ST_Intersects(a, b) AS a_intersects_b,
       ST_Contains(b, a) AS b_contains_a
FROM (
  SELECT ST_GeomFromText('POINT(1 0.5)') AS a,
         ST_GeomFromText('LINESTRING(1 0, 1 1, 3 5)') AS b
) AS foo;

 a_overlap_b | a_crosses_b | a_intersects_b | b_contains_a
-------------+-------------+----------------+--------------
 f           | f           | t              | t

例 — ST_Relate を使用して空間関係をテスト：

SELECT ST_Relate(ST_GeometryFromText('POINT(1 2)'), ST_Buffer(ST_GeometryFromText('POINT(1 2)'), 2), '0FFFFF212');

 st_relate
-----------
 t

例 — ST_Touches を使用してポイントがラインに接しているかをテスト：

-- 内部ポイント：接していない（ST_Touches は境界接触のみを対象とする）
SELECT ST_Touches('LINESTRING(0 0, 1 1, 0 2)'::geometry, 'POINT(1 1)'::geometry);
-- 結果：f

-- ラインの端点：接している
SELECT ST_Touches('LINESTRING(0 0, 1 1, 0 2)'::geometry, 'POINT(0 2)'::geometry);
-- 結果：t

例 — ST_Within を使用してネストされた円の包含を確認：

SELECT
  ST_Within(smallc, smallc) AS 小円_内_小円,
  ST_Within(smallc, bigc)   AS 小円_内_大円,
  ST_Within(bigc, smallc)   AS 大円_内_小円,
  ST_Within(ST_Union(smallc, bigc), bigc) AS 和集合_内_大円,
  ST_Within(bigc, ST_Union(smallc, bigc)) AS 大円_内_和集合,
  ST_Equals(bigc, ST_Union(smallc, bigc)) AS 大円_等_和集合
FROM (
  SELECT
    ST_Buffer(ST_GeomFromText('POINT(50 50)'), 20) AS smallc,
    ST_Buffer(ST_GeomFromText('POINT(50 50)'), 40) AS bigc
) AS foo;

 smallinsmall | smallinbig | biginsmall | unioninbig | biginunion | bigisunion
--------------+------------+------------+------------+------------+------------
 t            | t          | f          | t          | t          | t
(1 row)

ジオメトリの有効性チェック

ST_IsSimple を使用して、ジオメトリが自己交差していないかを確認し、ST_NRings を使用してポリゴン内のリング数をカウントします。

-- 有効な閉じたリング（シンプルポリゴン）
SELECT ST_IsSimple(ST_GeomFromText('POLYGON((1 2, 3 4, 5 6, 1 2))'));
-- 結果：t

-- 自己交差するラインストリング（シンプルではない）
SELECT ST_IsSimple(ST_GeomFromText('LINESTRING(1 1,2 2,2 3.5,1 3,1 2,2 1)'));
-- 結果：f

-- 複数のリングを持つ最大面積の市町村を検索
SELECT gid, name, ST_Area(the_geom) AS area
FROM bc_municipality
WHERE ST_NRings(the_geom) > 1
ORDER BY area DESC
LIMIT 1;

 gid | name   | area
-----+--------+------------------
  12 | Anning | 257374619.430216
(1 row)