Le blog de Jean David TECHER, un Réunionnais à Saint-Priest/Lyon

C'est en mattant le JT de France2 grâce à télévidilo que j'ai pu avoir cette image de Fillon:

Fillon - tout jeune d'avant - a sûrement trouvé le bon Filon en bossant avec NS . Filons à Fillon les dossiers qui défilent pour voir s'il se défile!

...Y'a pas à dire avec sa coiffe qui fait bien années 80!.........Je ne vais pas ici me mettre à critiquer l'équipe gouvernementale qui va être mis en place! Chacun aura son mot à dire! Le quinquennat est un exercice nouveau, qui change bien du septennat. Et même si Ségolène avait gagné, le chantier annoncé aurait été aussi énorme à entreprendre. Donc chaque jour qui passe et un jour important. Deux hommes différents avec une personnalité différente...On verra bien où tout celà nous mènera!

Mais mettre Bernard Kouchner au ministère des affaires étrangères - du moins s'il a accepté -, personnellement j'approuve ! Il y a beaucoup de dossiers qu'il connaît bien. Même si les deux hommes - lui et NS - ne sont pas d'accord sur bien des points, sur la politique étrangère à au moins 70%, il partage les mêmes idées et connaissent bien les dossiers épineux.

J'ai posté un message suite à question posée par un utilisateur de postgis dans la mailing-list à http://postgis.refractions.net/pipermail/postgis-users/2007-May/015680.html. Je n'ai pas essayé de la coder en C bien que celà soit possible mais en PL/PGSQL celà reste quand même aussi faisable. La facilité vient surtout du fait du type e géométrie employé : MULTILINESTRING. Pour cette dernière, PostGIS dispose de fonctionnalités intéressantes pour écrire une fonction en PL/PGSQL: pointn(), etc...

LA FONCTION

Il faut commencer par écrire un type en PostgreSQL qui aura la structure suivante

CREATE TYPE PointSetFromLinetype AS (idx_geom int4,idx_point  int4,point_x double precision, point_y double precision);

• idx_geom: l'indice de la nième LINESTRING;
• idx_point: qui sera le nième point dans la LINESTRING en cours;
• point_x, point_y: les coordonnées X et Y du point.

Quant à la fonction, on peut par exemple imaginer la fonction suivante

CREATE OR REPLACE FUNCTION points_from_linetype(geometry) RETURNS SETOF PointSetFromLinetype AS $$  
DECLARE
geom ALIAS FOR $1;
Idx_Geom int4 :=1;
Num_Geom int4 :=0;
Nb_Points_In_Ring int4;
Idx_Point int4 := 1;
j PointSetFromLinetype;
BEGIN 
/*
    Get the number of points
*/
SELECT INTO Num_Geom NumGeometries(multi(geom))+1 ;
/* Parsing the Geometry */
WHILE Idx_Geom < Num_Geom
     LOOP
         SELECT INTO Nb_Points_In_Ring npoints(GeometryN(multi(geom),Idx_Geom))+1;
          WHILE Idx_Point < Nb_Points_In_Ring
          LOOP
         	SELECT INTO j Idx_Geom,Idx_Point,x(pointn(GeometryN(multi(geom),Idx_Geom),Idx_Point)),y(pointn(GeometryN(multi(geom),Idx_Geom),Idx_Point));
         	RETURN NEXT  j;
      	        Idx_Point := Idx_Point + 1;
          END LOOP; -- End - FOR 1 
          Idx_Point := 1;
          Idx_Geom := Idx_Geom + 1;
     END  LOOP; -- End - while 
END;
$$ LANGUAGE plpgsql stable;

EXEMPLES D'UTILISATION

Pour le cas d'une simple LINESTRING, on aura

testgis=# SELECT * FROM points_from_linetype('LINESTRING(2 3,3 2,5 4)'::geometry);
 idx_geom | idx_point | point_x | point_y
----------+-----------+---------+---------
        1 |         1 |       2 |       3
        1 |         2 |       3 |       2
        1 |         3 |       5 |       4
(3 lignes)

Dans le cas d'un MULTILINESTRING, ce sera

testgis=# SELECT * FROM points_from_linetype('MULTILINESTRING((0 0,1 1,1 2),(2 3,3 2,5 4))'::geometry);
 idx_geom | idx_point | point_x | point_y
----------+-----------+---------+---------
        1 |         1 |       0 |       0
        1 |         2 |       1 |       1
        1 |         3 |       1 |       2
        2 |         1 |       2 |       3
        2 |         2 |       3 |       2
        2 |         3 |       5 |       4
(6 lignes)

Je ne fais que relayer un billet de postgresqlfr.org paru hier à http://www.postgresqlfr.org/?q=node/1233. Comme expliqué, il faudra se préparer à migrer bientôt pour la 8.3 pour les développeurs sous Windows.

Ayé j'ai fait mon premier fichier kml personnel - hors contexte professionnel bien sûr -. Il est disponible à http://www.davidgis.fr/download/jeandavidtecher.kml

Le bâtiment en question représente le bâtiment où je vis à savoir la Résidence Prado Lido. Le polygone qui est dessiné a pour représentation dans le système de projection WGS 84 pour une hauteur de 12 mètres -un peu ficitf je l'avoue , je ne connais pas vraiment la hauteur de mon bâtiment -:

SELECT GeomFromEWKT('SRID=4326;POLYGON((3.896026111319043 43.63077293064538 12,3.896125412083016 43.63129020374376 12,
3.896017722403089 43.63130336362247 12,3.896023316938252 43.6313529652266 12,
3.895547802871119 43.63140560406845 12,3.895538013117296 43.63135701466003 12,
3.895430322945158 43.63136612505183 12,3.895325433631279 43.63084986305638 12,
3.895428927275387 43.6308407529833 12,3.895543607566332 43.63132462175059 12,
3.896003736551513 43.63126793394599 12,3.895911429755532 43.63078305367532 12,
3.896026111319043 43.63077293064538 12))');