SpatiaLite Cookbook

	Recette #18 Railway et tampons
Février 2011

Table des matières

Le problème

Ceci est une sorte de conclusion visuelle de l'exercice précédent. le problème est le suivant:

créer une couche représentant les zones A, B, C, D et E précédemment définies.

CREATE TABLE railway_zones (
  id INTEGER NOT NULL PRIMARY KEY AUTOINCREMENT,
  railway_name TEXT NOT NULL,
  zone_name TEXT NOT NULL);

SELECT AddGeometryColumn('railway_zones', 'geometry',
23032, 'MULTIPOLYGON', 'XY');

Commençons par créer une nouvelle table:

Comme d'habitude, on définit d'abord la table, sans considérer la colonne géométrique
Puis nous créons la colonne géométrique, en utilisant AddGeometryColumn().
Spécifions le SRID 23032 [ED50 UTM zone 32] pour cette table (la table railway ayant le même SRID)
Spécifions le type géométrique MULTYPOLYGON : nous verrons plus tard les raisons de ce choix

INSERT INTO railway_zones
  (id, railway_name, zone_name, geometry)
SELECT NULL, name, 'A class [< 1Km]',
  CastToMultiPolygon(
    ST_Buffer(geometry, 1000.0))
FROM railways;

Rien de bien intéressant dans cette clause INSERT INTO ... SELECT ...
a part ces quelques points:

on utilise ST_Buffer() afin de créer un tampon de type POLYGON correspondant à la zone de classe A ( 1 km ).
note: il n'est pas toujours évident de deviner le type de géométrie créé par un tampon ( ST_Buffer );
parfois elle va créer des POLYGON, et d'autres des MULTYPOLYGON
(ceci dépend de la forme de la géométrie source, ainsi que du rayon choisi pour le tampon).
ainsi, c'est afin d'éviter toute incohérence, que nous avons spécifié une géométrie de type MULTIPOLYGON pour cette table.
on force maintenant ce type géométrique via la fonction CastToMultiPolygon()(cette fonction va transformer la géométrie résultante d'un tampon en multipolygon)

INSERT INTO railway_zones
  (id, railway_name, zone_name, geometry)
SELECT NULL, name, 'B class [< 2.5Km]',
  CastToMultiPolygon(
    ST_Difference(
      ST_Buffer(geometry, 2500.0),
      ST_Buffer(geometry, 1000.0)))
FROM railways;

INSERT INTO railway_zones
  (id, railway_name, zone_name, geometry)
SELECT NULL, name, 'C class [< 5Km]',
  CastToMultiPolygon(
    ST_Difference(
      ST_Buffer(geometry, 5000.0),
      ST_Buffer(geometry, 2500.0)))
FROM railways;

INSERT INTO railway_zones
  (id, railway_name, zone_name, geometry)
SELECT NULL, name, 'D class [< 10Km]',
  CastToMultiPolygon(
    ST_Difference(
      ST_Buffer(geometry, 10000.0),
      ST_Buffer(geometry, 5000.0)))
FROM railways;

INSERT INTO railway_zones
  (id, railway_name, zone_name, geometry)
SELECT NULL, name, 'E class [< 20Km]',
  CastToMultiPolygon(
    ST_Difference(
      ST_Buffer(geometry, 20000.0),
      ST_Buffer(geometry, 10000.0)))
FROM railways;

Creer les autres zones n'est pas plus compliqué.

won utilise simplement ST_Difference() afin d'obtenir la représentation désirée:
en d'autres termes, on doit créer des trous dans chaque zone afin d'exclure les autres zones déjà créées.

Vous pouvez dès à présent visualiser le résultat avec QGIS. Et c'est tout.

Table des matières

	Author: Alessandro Furieri a.furieri@lqt.it Traduced from English By RIVIERE Romain
	This work is licensed under the Attribution-ShareAlike 3.0 Unported (CC BY-SA 3.0) license.

	Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.3 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts.