Bienvenue sur PostGIS.fr

Bienvenue sur PostGIS.fr , le site de la communauté des utilisateurs francophones de PostGIS.

PostGIS ajoute le support d'objets géographique à la base de données PostgreSQL. En effet, PostGIS "spatialise" le serverur PostgreSQL, ce qui permet de l'utiliser comme une base de données SIG.

Maintenu à jour, en fonction de nos disponibilités et des diverses sorties des outils que nous testons, nous vous proposons l'ensemble de nos travaux publiés en langue française.

source: trunk/workshop-foss4g/geometry_returning.rst @ 36

Revision 36, 9.1 KB checked in by djay, 13 years ago (diff)

fix typo

Partie 18 : fonction de construction de géométries

System Message: WARNING/2 (<string>, line 4)

Title underline too short.

Partie 18 : fonction de construction de géométries
==================================================

Toute les fonctions que nous avons vu jusqu'à présent traitent les géométries "comme elles sont" est retournent

Les "fonctions de construction de géométries" prennent des géométries en entrée et retourne de nousvelles formes.

ST_Centroid / ST_PointOnSurface

Un besoin commun lors de la création de requêtes spatiales est de remplacer une entité polygonale par un point représentant cette entité. Cela esr utile pour les jointures spatiales (comme indiqué ici : :ref:`polypolyjoins`) car utiliser :command:`ST_Intersects(geometry,geometry)` avec deux polygones impliquera un double comptage : un polygone pour le contour externe intersectera dans les deux senses; le replacer par un point le forcera à être dans un seul sens, pas les deux.

System Message: ERROR/3 (<string>, line 19); backlink

Unknown interpreted text role "ref".

System Message: ERROR/3 (<string>, line 19); backlink

Unknown interpreted text role "command".
  • :command:`ST_Centroid(geometry)` retourne le point qui est approximativement au centre de la masse de la géométrie passé en paramÚtre. C'est un calcul simle et rapide, mais parfois non proftable, car le point retourné peut se trouver à l'extérieur de l'entité elle-même. Si l'entité fournie est convee (imaginez la lettre 'C') le centroïde renvoyé pourrait ne pas être à l'intérieur du polygone.

    System Message: ERROR/3 (<string>, line 21); backlink

    Unknown interpreted text role "command".

  • :command:`ST_PointOnSurface(geometry)` retourne un point qui est obligatoirement dans l'entité passée en paramÚtre. Cette fonction coûte plus cher en ressource que le calcul du centroïd.

    System Message: ERROR/3 (<string>, line 22); backlink

    Unknown interpreted text role "command".

./geometry_returning/centroid.jpg

ST_Buffer

L'opération de zone tampon est souvent disponible dans les outils SIG, il est aussi disponible dans PostGIS. La fonction :command:`ST_Buffer(geometry,distance)` prend en paramÚtre une géométrie et une distance et retourne une zone tampon dont le contour est à une distance données de la géométrie d'origine.

System Message: ERROR/3 (<string>, line 30); backlink

Unknown interpreted text role "command".
./geometry_returning/st_buffer.png

Par exemple, si les services des parks américains souhaitaient renforcer la zone du traffique maritime autour de l'île 'Liberty', ils pourraient construire une zone tampon de 500 mÚtres autour de l'île. L'île de 'Liiberty' est représenté par un seul bloque dans notre table ``nyc_census_blocks`, nous pouvons donc facilement réaliser ce calcul.

System Message: WARNING/2 (<string>, line 34); backlink

Inline literal start-string without end-string.
-- Création d'une nouvelle table avec une zone tampon de 500 m autour de 'Liberty Island'
CREATE TABLE libery_island_zone AS
SELECT ST_Buffer(the_geom,500) AS the_geom
FROM nyc_census_blocks
WHERE blkid = '360610001009000';
-- Mise à jour de la table geometry_columns
SELECT Populate_Geometry_Columns();
./geometry_returning/liberty_positive.jpg

La fonction :command:`ST_Buffer` permet aussi d'utiliser des valeur négative pour le paramÚtre distance et construit un polygone inclue dans celui passé en paramÚtre. Pour les points et les lignes vous obtiendrez simplement un résultat vide.

System Message: ERROR/3 (<string>, line 49); backlink

Unknown interpreted text role "command".
./geometry_returning/liberty_negative.jpg

ST_Intersection

Une autre opération classique présente dans les SIGS - le chevauchement - crée une nouvelle entité en calculant la zone correpondant à l'intersection de deux polygones supperposés. Le résultat à la propriété de permettre de reconstruire les entité de base à l'aide de ce résultat.

La fonction :command:`ST_Intersection(geometry A, geometry B)` retourne la zone géographique (ou une ligne, ou un point) que les deux géométries on en commun. Si les géométries sont disjointes, la fontion retourne une géométrie vide.

System Message: ERROR/3 (<string>, line 59); backlink

Unknown interpreted text role "command".
-- Quelle est l'aire que ces deux cercles ont en commun ?
-- Utilisons la fonction ST_Buffer pour créer ces cercles !

SELECT ST_AsText(ST_Intersection(
  ST_Buffer('POINT(0 0)', 2),
  ST_Buffer('POINT(3 0)', 2)
));
./geometry_returning/intersection.jpg

ST_Union

Dans l'exemple précédent nous intersections des géométries, créant une nouvelle géométrie unique à partir de deux entités. La commade :command:`ST_Union` fait l'inverse, elle prend en paramÚtre des géométries et supprime les parties communes. Il y a deux versions possibles de la fonction :command:`ST_Union` :

System Message: ERROR/3 (<string>, line 78); backlink

Unknown interpreted text role "command".

System Message: ERROR/3 (<string>, line 78); backlink

Unknown interpreted text role "command".
  • :command:`ST_Union(geometry, geometry)`: une version avec deux paramÚtres qui prend les géométries et rentourne l'union des deux. Par exemple, nos deux cercles ressemble à ce qui suit si nous utilisons l'opération union plutÃŽt que l'intersection.

    System Message: ERROR/3 (<string>, line 80); backlink

    Unknown interpreted text role "command".

    -- Quelle est l'aire totale des ces deux cercles ?
    -- Utilisons ST_Buffer pour créer les cercles !
    
    SELECT ST_AsText(ST_Union(
      ST_Buffer('POINT(0 0)', 2),
      ST_Buffer('POINT(3 0)', 2)
    ));
    
    ./geometry_returning/union.jpg
  • :command:`ST_Union([geometry])`: une version agrégée qui prendre un ensemble de géométries et retourne une géométrie contenant l'ensemble des géométries rassemblées. La fonction égrégée ST_Union peut être utilisé grâce au SQL GROUP BY our créer un ensemble rassemblant des sous-ensembles de géométries basiques. Cela est trÚs puissant,

    System Message: ERROR/3 (<string>, line 95); backlink

    Unknown interpreted text role "command".

Par exemple la fonction d'agrégation :command:`ST_Union`, considÚrons notre table nyc_census_blocks. As an example of :command:`ST_Union` aggregation, consider our nyc_census_blocks table. Census geography is carefully constructed so that larger geographies can be built up from smaller ones. So, we can create a census tracts map by merging the blocks that form each tract (as we do later in :ref:`creatingtractstable`). Or, we can create a county map by merging blocks that fall within each county.

System Message: ERROR/3 (<string>, line 97); backlink

Unknown interpreted text role "command".

System Message: ERROR/3 (<string>, line 97); backlink

Unknown interpreted text role "command".

System Message: ERROR/3 (<string>, line 97); backlink

Unknown interpreted text role "ref".

To carry out the merge, note that the unique key blkid actually embeds information about the higher level geographies. Here are the parts of the key for Liberty Island we used earlier:

360610001009000 = 36 061 00100 9000
36     = State of New York
061    = New York County (Manhattan)
000100 = Census Tract
9      = Census Block Group
000    = Census Block

So, we can create a county map by merging all geometries that share the same first 5 digits of their blkid.

-- Création d'une table nyc_census_counties en regroupant les bloques
CREATE TABLE nyc_census_counties AS
SELECT
  ST_Union(the_geom) AS the_geom,
  SubStr(blkid,1,5) AS countyid
FROM nyc_census_blocks
GROUP BY countyid;
-- Mise à jour de la table geometry_columns
SELECT Populate_Geometry_Columns();
./geometry_returning/union_counties.png

An area test can confirm that our union operation did not lose any geometry. First, we calculate the area of each individual census block, and sum those areas grouping by census county id.

SELECT SubStr(blkid,1,5) AS countyid, Sum(ST_Area(the_geom)) AS area
FROM nyc_census_blocks
GROUP BY countyid;
 countyid |       area
----------+------------------
 36005    | 109807439.720947
 36047    | 184906575.839355
 36061    | 58973521.6225586
 36081    | 283764734.207275
 36085    | 149806077.958252

Ensuite nous calculons l'aire de chaque zone de nos nouveaux polygones de régions de la table count :

SELECT countyid, ST_Area(the_geom) AS area
FROM nyc_census_counties;
 countyid |       area
----------+------------------
 36005    | 109807439.720947
 36047    | 184906575.839355
 36061    | 58973521.6225586
 36081    | 283764734.207275
 36085    | 149806077.958252

La même réponse ! Nous avons construit avec succÚs une table des régions de NYC à partir de nos données initiales.

Liste des fonctions

ST_AsText(text): retourne la représentation Well-Known Text (WKT) de la geometry/geography sans métadonnée SRID.

ST_Buffer(geometry, distance): For geometry: Returns a geometry that represents all points whose distance from this Geometry is less than or equal to distance. Calculations are in the Spatial Reference System of this Geometry. For geography: Uses a planar transform wrapper.

ST_Intersection(geometry A, geometry B): Returns a geometry that represents the shared portion of geomA and geomB. The geography implementation does a transform to geometry to do the intersection and then transform back to WGS84.

ST_Union(): Returns a geometry that represents the point set union of the Geometries.

substring(string [from int] [for int]): PostgreSQL string function to extract substring matching SQL regular expression.

sum(expression): PostgreSQL aggregate function that returns the sum of records in a set of records.

Note: See TracBrowser for help on using the repository browser.