Partie 10 : Les relations spatiales

Jusqu’à présent, nous avons utilisé uniquement des fonctions qui permettent de mesurer (ST_Area, ST_Length), de sérialiser (ST_GeomFromText) ou désérialiser (ST_AsGML) des géométries. Ces fonctions sont toutes utilisées sur une géométrie à la fois.

Les base de données spatiales sont puissantes car elle ne se contentent pas de stocker les géométries, elle peuvent aussi vérifier les relations entre les géométries.

Pour les questions comme “Quel est le plus proche garage à vélo près du parc ?” ou “Ou est l’intersection du métro avec telle rue ?”, nous devrons comparer les géométries représentant les garages à vélo, les rues et les lignes de métro.

Le standard de l’OGC définit l’ensemble de fonctions suivantes pour comparer les géométries.

ST_Equals

ST_Equals(geometry A, geometry B) teste l’égalité spatiale de deux géométries.

ST_Equals retourne TRUE si les deux géométries sont du même type et ont des coordonnées x.y identiques.

Premièrement, essayons de récupérer la représentation d’un point de notre table nyc_subway_stations. Nous ne prendrons que l’entrée : ‘Broad St’.

SELECT name, the_geom, ST_AsText(the_geom)
FROM nyc_subway_stations
WHERE name = 'Broad St';

   name   |                      the_geom                      |      st_astext
----------+----------------------------------------------------+-----------------------
 Broad St | 0101000020266900000EEBD4CF27CF2141BC17D69516315141 | POINT(583571 4506714)

Maintenant, copiez / collez la valeur affichée pour tester la fonction ST_Equals:

SELECT name
FROM nyc_subway_stations
WHERE ST_Equals(the_geom, '0101000020266900000EEBD4CF27CF2141BC17D69516315141');

Broad St

Note

La représentation du point n’est pas vraiment compréhensible (0101000020266900000EEBD4CF27CF2141BC17D69516315141) mais c’est exactement la représentation des coordonnées. Pour tester l’égalité, l’utilisation de ce format est nécessaire.

ST_Intersects, ST_Disjoint, ST_Crosses et ST_Overlaps

ST_Intersects, ST_Crosses, et ST_Overlaps teste si l’intérieur des géométries s’intersecte, se croise ou se chevauche.

ST_Intersects(geometry A, geometry B) retourne t (TRUE) si l’intersection ne renvoie pas un ensemble vide de résultats. Intersects retourne le résultat exactement inverse de la fonction disjoint.

L’opposé de ST_Intersects est ST_Disjoint(geometry A , geometry B). Si deux géométries sont disjointes, elle ne s’intersectent pas et vice-versa. En fait, il est souvent plus efficace de tester si deux géométries ne s’intersectent pas que de tester si elles sont disjointes du fait que le test d’intersection peut être spatialement indexé alors que le test disjoint ne le peut pas.

Pour les comparaisons de couples de types multipoint/polygon, multipoint/linestring, linestring/linestring, linestring/polygon, et linestring/multipolygon, ST_Crosses(geometry A, geometry B) retourne t (TRUE) si les résultats de l’intersection sont à l’intérieur des deux géométries.

ST_Overlaps(geometry A, geometry B) compare deux géométries de même dimension et retourne TRUE si leur intersection est une géométrie différente des deux fournies mais de même dimension.

Essayons de prendre la station de métro de Broad Street et de déterminer son voisinage en utilisant la fonction ST_Intersects :

SELECT name, boroname
FROM nyc_neighborhoods
WHERE ST_Intersects(the_geom, '0101000020266900000EEBD4CF27CF2141BC17D69516315141');

        name        | boroname
--------------------+-----------
 Financial District | Manhattan

ST_Touches

ST_Touches teste si deux géométries se touchent en leur contours extérieurs, mais leur contours intérieurs ne s’intersectent pas

ST_Touches(geometry A, geometry B) retourne TRUE soit si les contours des géométries s’intersectent ou si l’un des contours intérieurs de l’une intersecte le contour extérieur de l’autre.

ST_Within et ST_Contains

ST_Within et ST_Contains teste si une géométrie est totalement incluse dans l’autre.

ST_Within(geometry A , geometry B) retourne TRUE si la première géométrie est complètement contenue dans l’autre. ST_Within teste l’exact opposé au résultat de ST_Contains.

ST_Contains(geometry A, geometry B) retourne TRUE si la seconde géométrie est complètement contenue dans la première géométrie.

ST_Distance et ST_DWithin

Une question fréquente dans le domaine du SIG est “trouver tous les éléments qui se trouvent à une distance X de cet autre élément”.

La fonction ST_Distance(geometry A, geometry B) calcule la plus courte distance entre deux géométries. Cela est pratique pour récupérer la distance entre les objets.

SELECT ST_Distance(
  ST_GeometryFromText('POINT(0 5)'),
  ST_GeometryFromText('LINESTRING(-2 2, 2 2)'));

3

Pour tester si deux objets sont à la même distance d’un autre, la fonction ST_DWithin fournit un test tirant profit des index. Cela est très utile pour répondre a une question telle que: “Combien d’arbres se situent dans un buffer de 500 mètres autour de cette route ?”. Vous n’avez pas à calculer le buffer, vous avez simplement besoin de tester la distance entre les géométries.

En utilisant de nouveau notre station de métro Broad Street, nous pouvons trouver les rues voisines (à 10 mètres de) de la station :

SELECT name
FROM nyc_streets
WHERE ST_DWithin(
        the_geom,
        '0101000020266900000EEBD4CF27CF2141BC17D69516315141',
        10
      );

     name
--------------
   Wall St
   Broad St
   Nassau St

Nous pouvons vérifier la réponse sur une carte. La station Broad St est actuellement à l’intersection des rues Wall, Broad et Nassau.

Liste des fonctions

ST_Contains(geometry A, geometry B) : retourne TRUE si aucun des points de B n’est à l’extérieur de A, et au moins un point de l’intérieur de B est à l’intérieur de A.

ST_Crosses(geometry A, geometry B) : retourne TRUE si la géométrie A a certains, mais pas la totalité, de ses points à l’intérieur de B.

ST_Disjoint(geometry A , geometry B) : retourne TRUE si les géométries ne s’intersectent pas - elles n’ont aucun point en commun.

ST_Distance(geometry A, geometry B) : retourne la distance cartésienne en 2 dimensions minimum entre deux géométries dans l’unité de la projection.

ST_DWithin(geometry A, geometry B, radius) : retourne TRUE si les géométries sont distante (radius) l’une de l’autre.

ST_Equals(geometry A, geometry B) : retourne TRUE si les géométries fournies représentent la même géométrie. L’ordre des entités n’est pas pris en compte.

ST_Intersects(geometry A, geometry B) : retourne TRUE si les géométries s’intersectent - (ont un espace en commun) et FALSE si elles n’en ont pas (elles sont disjointes).

ST_Overlaps(geometry A, geometry B) : retourne TRUE si les géométries ont un espace en commun, sont de la même dimension, mais ne sont pas complètement contenues l’une dans l’autre.

ST_Touches(geometry A, geometry B) : retourne TRUE si les géométries ont au moins un point en commun, mais leur intérieurs ne s’intersectent pas.

ST_Within(geometry A , geometry B) : retourne TRUE si la géométrie A est complètement à l’intérieur de B