source:
trunk/workshop-foss4g/spatial_relationships.rst
@
31
Revision 22, 9.1 KB checked in by djay, 13 years ago (diff) |
---|
Partie 10 : relations spatiales
Jusqu'à maintenant nous n'avons utilisé que des fonctions qui permettent de mesurer (:command:`ST_Area`, :command:`ST_Length`), de serialiser (:command:`ST_GeomFromText`) ou désérialiser (:command:`ST_AsGML`) des géométries. Ces fonctions ont en commun de fonctionner uniquement sur une géométrie à la fois.
Les base de données spatiales sont puissantes car elle ne font pas que stoquer les géométries, elle ont aussi la faculté de vérifier les relations entre les géométries.
Pour les questions comme "Quel est le plus proche garage à vélo prêt du park ?" ou "Ou est l'intersection du métros avec telle rue ?" nous devrons comparer les géométries représentant les garage à vélo, les rues et les lignes de métros.
Le standard de l'OGC définit l'ensemble suivant de fonctions pour comparer les géométries.
ST_Equals
:command:`ST_Equals(geometry A, geometry B)` test l'égalité spatiale de deux géométries.
ST_Equals retourne TRUE si les deux géométries sont du même type, 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 afficher pour tester la fonction :command:`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é, utiliser ce format est nécessaire.
ST_Intersects, ST_Disjoint, ST_Crosses et ST_Overlaps
:command:`ST_Intersects`, :command:`ST_Crosses`, et :command:`ST_Overlaps` test si l'intérieur des géométries s'intersect, se croise ou se chevauche.
:command:`ST_Intersects(geometry A, geometry B)` retourne t (TRUE) si l'intersection ne rénvoit pas un ensemble vide de résultats. Intersects retourne le résultat exactement inverse de la fonction disjoint.
L'opposé de ST_Intersects est :command:`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 éfficace de tester si deux géométries ne s'intersectent pas que de tester si elles sont dijointes 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, :command:`ST_Crosses(geometry A, geometry B)` retourne t (TRUE) si les résultats de l'intersection
For multipoint/polygon, multipoint/linestring, linestring/linestring, linestring/polygon, and linestring/multipolygon comparisons, :command:`ST_Crosses(geometry A, geometry B)` returns t (TRUE) if the intersection results in a geometry whose dimension is one less than the maximum dimension of the two source geometries and the intersection set is interior to both source geometries.
:command:`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 sont voisinage en utilisant la fonction :command:`ST_Intersects` :
SELECT name, boroname FROM nyc_neighborhoods WHERE ST_Intersects(the_geom, '0101000020266900000EEBD4CF27CF2141BC17D69516315141');
name | boroname --------------------+----------- Financial District | Manhattan
ST_Touches
:command:`ST_Touches` test si deux géométries se touchent en leur contour extérieur, mais leur contours intérieur ne s'intersectent pas
:command:`ST_Touches(geometry A, geometry B)` retourn TRUE soit si les contours des géométries s'intersectent ou si l'un des contours intérieur de l'une intersecte le contour extérieur de l'autre.
ST_Within et ST_Contains
:command:`ST_Within` et :command:`ST_Contains` test si une géométrie est totalement incluse dans l'autre.
:command:`ST_Within(geometry A , geometry B)` retourne TRUE si la premiÚre géométries est complÚtement contenue dans l'autre. ST_Within test l'exact opposé au résultat de ST_Contains.
:command:`ST_Contains(geometry A, geometry B)` retourne TRUE si la seconde géométries est complÚtement contenue dans la premiÚre géométrie.
ST_Distance et ST_DWithin
Une question qui arrive fréquemment dans le domaine du SIG est "trouver tout les trucs qui se trouve à une distance X de cet autre truc".
La fonction :command:`ST_Distance(geometry A, geometry B)` calcule la plus courte distance entre deux géoémé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 obets sont à la même distance d'un autre, la fonction :command:`ST_DWithin` fournit une test tirant proffit des indexes. Cela est trÚs utile pour répondre au questions du genre : "Combien d'arbre se situe 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étros 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émétries nes 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) : retourn TRUE si les géométries fournis représentent la même géométrie. L'ordre des entités n'est pas prit 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 contenu 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