Dans cette section nous allons creuser un peu plus profondément comment les géométries simples sont construites en SQL. En réalité, vous utiliserez probablement un SIG comme QGIS pour créer des géométries complexes en utilisant leurs outils de numérisation ; cependant, comprendre comment ils sont formulés peut être pratique pour écrire des requêtes et comprendre comment la base de données est assemblée.
Objectif de cette leçon : Mieux comprendre comment créer des entités spatiales directement dans PostgreSQL/PostGIS.
Retournons à notre base de données address et faisons en sorte que la table des rues corresponde à celles précédemment étudiées, c’est-à-dire qu’il existe une contrainte géométrique, un index et en entrée dans la table de métadonnées geometry_columns.
Modifiez la table streets de sorte à ce qu’il y ait une colonne géométrie de type ST_PolyLigne.
N’oubliez pas de faire la mise à jour d’accompagnement à la table des colonnes géométriques !
Ajoutez également une contrainte pour empêcher toutes géométries ajoutées qui ne sont pas des POLYLIGNES ou nulles.
Créez un index spatial sur la nouvelle colonne géométrique
Insérons maintenant une polyligne dans notre table rues. Dans ce cas nous mettrons à jour un enregistrement de rue existant:
update streets
set the_geom = 'SRID=4326;LINESTRING(20 -33, 21 -34, 24 -33)'
where streets.id=2;
Jetez un coup d’œil aux résultats dans QGIS. (Vous pouvez devoir faire un clic droit sur la couches streets dans la ‘Légende’ des couches, et choisissez ‘Zoomer sur l’emprise de la couche’.)
Créez maintenant plus d’entrées streets - certaines dans QGIS et d’autres depuis la ligne de commande.
Créer des polygones est tout aussi facile. Une chose à se souvenir est que par définition les polygones ont au moins quatre sommets, avec le dernier et le premier co-localisés:
insert into cities (name, the_geom)
values ('Tokyo', 'SRID=4326;POLYGON((10 -10, 5 -32, 30 -27, 10 -10))');
Note
Un polygone requière des doubles parenthèses autour de sa liste de coordonnées ; c’est pour vous permettre d’ajouter des polygones complexes avec de multiples aires non connectées. Par exemple
insert into cities (name, the_geom)
values ('Tokyo Outer Wards',
'SRID=4326;POLYGON((20 10, 20 20, 35 20, 20 10),
(-10 -30, -5 0, -15 -15, -10 -30))'
);
Si vous suivez cette étape, vous pouvez vérifier ce qui a été fait en chargeant le jeu de données villes dans QGIS, en ouvrant sa table attributaire, et en sélectionnant la nouvelle entrée. Notez comment les deux nouveaux polygones se comportent comme un polygone.
Notre schéma devrait désormais ressembler à ça :
Créez des contours de ville en calculant l’enveloppe convexe minimum de toutes les adresses pour cette ville et en calculant un tampon autour de cette zone.
Avec les fonctions du Modèle-SFS, vous avez une large variété d’options pour accéder aux sous-objets des Géométries SFS. Lorsque vous voulez sélectionner le premier point du sommet de chaque géométrie de polygone dans la table myPolygonTable, vous devez faire cela de cette manière :
Transformez les contours du polygone en une polyligne:
select st_boundary(geometry) from myPolygonTable;
Sélectionnez le premier point du sommet de la polyligne résultante :
select st_startpoint(myGeometry)
from (
select st_boundary(geometry) as myGeometry
from myPolygonTable) as foo;
PostGIS supporte toutes les fonctions conformes à la norme OGC SFS/MM. Toutes ces fonctions commencent avec ST_.
Pour découper une sous-partie de vos données, vous pouvez utiliser la fonction ST_INTERSECT() Pour éviter les géométries vides, utilisez:
where not st_isempty(st_intersection(a.the_geom, b.the_geom))
select st_intersection(a.the_geom, b.the_geom), b.*
from clip as a, road_lines as b
where not st_isempty(st_intersection(st_setsrid(a.the_geom,32734),
b.the_geom));
À partir d’une table d’un point donné, vous voulez générer une polyligne. L’ordre des points est défini par leur id. Une autre méthode d’ordre pourrait être un horodatage, comme celui que vous obtenez lorsque vous capturez des points de navigation avec un récepteur GPS.
Pour créer une polyligne à partir d’une nouvelle couche de points appelée ‘points’, vous pouvez exécuter la commande suivante:
select ST_LineFromMultiPoint(st_collect(the_geom)), 1 as id
from (
select the_geom
from points
order by id
) as foo;
Pour voir comme cela fonctionne sans créer une nouvelle couche, vous pourriez aussi lancer cette commande sur la couche ‘people’, même si faire cela dans le monde réel aurait peu de sens.
Pour déceler la différence entre deux tables ayant la même structure, vous pouvez utiliser le mot-clé de PostgreSQL``EXCEPT``:
select * from table_a
except
select * from table_b;
Comme résultat, vous obtiendrez tous les enregistrements de la table_a qui ne sont pas stockés dans la table_b.
Vous pouvez définir où Postgre doit stocker ses données sur le disque en créant des espaces de table:
CREATE TABLESPACE homespace LOCATION '/home/pg';
Lorsque vous créez une base de données, vous pouvez indiquer quel espace de table utiliser, ex:
createdb --tablespace=homespace t4a
Vous avez appris comment créer des géométries complexes en utilisant les déclarations PostGIS. Gardez à l’esprit que vous avez amélioré votre connaissance globale sur les bases de données spatiales, à travers un outil SIG. Vous n’avez généralement pas besoin d’utiliser ces déclarations manuellement mais disposer de compétences sur leur utilisation et leur structure vous aidera à utiliser un SIG, particulièrement si vous rencontrez des erreurs qui pourraient apparaître comme incompréhensibles autrement.