AnalyticDB:Model geometri

Ukur geometri

Gunakan ST_Length dan ST_Area untuk menghitung jarak dan luas. Kedua fungsi ini mengembalikan hasil dalam satuan sistem referensi koordinat geometri tersebut (meter untuk SRID 3857).

-- Total panjang jalan dalam kilometer
SELECT SUM(ST_Length(the_geom)) / 1000 AS km_roads FROM bc_roads;

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

-- Luas wilayah administratif tertentu dalam hektar
SELECT ST_Area(the_geom) / 10000 AS hectares
FROM bc_municipality
WHERE name = 'PRINCE GEORGE';

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

Uji hubungan spasial

GanosBase Geometry mendukung seluruh fungsi hubungan spasial PostGIS. Tabel berikut mencantumkan fungsi-fungsi yang paling umum digunakan.

Contoh — cari panjang jalan per wilayah administratif menggunakan 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;

Contoh — bandingkan cakupan antara dua lingkaran menggunakan ST_Covers:

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 dan ST_Contains berbeda dalam cara menangani titik batas. ST_Covers mengembalikan true ketika A menutupi batas B; ST_Contains tidak demikian.

Contoh — uji titik dan garis yang saling lepas menggunakan ST_Disjoint:

-- Titik tidak berada pada garis: saling lepas
SELECT ST_Disjoint('POINT(0 0)'::geometry, 'LINESTRING ( 2 0, 0 2 )'::geometry);
-- Hasil: t

-- Titik berada pada garis: tidak saling lepas
SELECT ST_Disjoint('POINT(0 0)'::geometry, 'LINESTRING ( 0 0, 0 2 )'::geometry);
-- Hasil: f

Contoh — uji tumpang tindih dan interseksi menggunakan 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

Contoh — uji hubungan spasial menggunakan ST_Relate:

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

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

Contoh — uji apakah suatu titik menyentuh garis menggunakan ST_Touches:

-- Titik interior: tidak menyentuh (ST_Touches hanya berlaku untuk kontak batas)
SELECT ST_Touches('LINESTRING(0 0, 1 1, 0 2)'::geometry, 'POINT(1 1)'::geometry);
-- Hasil: f

-- Titik ujung garis: menyentuh
SELECT ST_Touches('LINESTRING(0 0, 1 1, 0 2)'::geometry, 'POINT(0 2)'::geometry);
-- Hasil: t

Contoh — periksa kandungan untuk lingkaran bersarang menggunakan ST_Within:

SELECT
  ST_Within(smallc, smallc) AS smallinsmall,
  ST_Within(smallc, bigc)   AS smallinbig,
  ST_Within(bigc, smallc)   AS biginsmall,
  ST_Within(ST_Union(smallc, bigc), bigc) AS unioninbig,
  ST_Within(bigc, ST_Union(smallc, bigc)) AS biginunion,
  ST_Equals(bigc, ST_Union(smallc, bigc)) AS bigisunion
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)

Periksa validitas geometri

Gunakan ST_IsSimple untuk memverifikasi bahwa geometri tidak memiliki self-intersection, dan ST_NRings untuk menghitung jumlah ring dalam sebuah poligon.

-- Ring tertutup yang valid (poligon simple)
SELECT ST_IsSimple(ST_GeomFromText('POLYGON((1 2, 3 4, 5 6, 1 2))'));
-- Hasil: t

-- Linestring yang saling berpotongan (tidak simple)
SELECT ST_IsSimple(ST_GeomFromText('LINESTRING(1 1,2 2,2 3.5,1 3,1 2,2 1)'));
-- Hasil: f

-- Temukan wilayah administratif dengan luas terbesar yang memiliki lebih dari satu ring
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)