djay's blog

Sortie de PostGIS 1.4.1

Thu, 17 Dec 2009 14:22:07 +0100

Ce matin, la nouvelle version de PostGIS 1.4.1 a été annoncée par Paul Ramsey, le lien pour télécharger cette version est ici.

La version 1.4.1 de PostGIS est maintenant disponible. C'est une version mineure qui permet de résoudre certain problèmes rencontrés depuis la version 1.4.0.

#241, problème avec ST_LineCrossingDirection
#210, erreur de segmentation dans ST_Union pour les valeurs NULLs
#277, problème lors de l'utilisation de grandes valeurs numériques en entrée
#179, ST_MakeLine crash coté serveur avec des valeurs NULL
#253, resultat inconsistents lors de l'utilisation de l'opérateur ~=
#276, ST_AsGML produisant des chaînes GML non compatibles
#239, problème mémoire avec les fonctions ST_As*
#272, ST_LineCrossingDirection
#316, problème de récupération des géométries invalides avec une plus grande dimension

La liste des présentations FOSS4G 2009

Mon, 20 Jul 2009 13:50:04 +0200

La liste des présentations acceptées par le comité de sélection du prochain FOSS4G qui se tiendra à Sydney a été publiée aujourd'hui. Une quantité impressionnante de présentations a été proposée cette année, plus de 183 présentations pour seulement 85 sélectionnées, sans compter les 52 articles scientifiques proposés dont seulement 12 ont été retenus.

Au passage, cous constaterez sans doute que bon nombre de présentations mentionnent l'utilisation de PostGIS dans les solutions développées et qui seront présentées lors du FOSS4G 2009.

Autre point qui me semble important : le retour du grand Paul Ramsey à cette conférence pour la présentation intitulée The State of PostGIS, il sera aussi responsable de workshop mais ce n'est pas le sujet ici...

Problème avec les couches PostGIS depuis MapServer avec PostgreSQL 8.3

Wed, 03 Jun 2009 10:52:22 +0200

Lors de migration de vos serveur de bases de données PostgreSQL vers une version 8.3, un problème survient avec les versions de MapServer mises à disposition dans MS4W. Pour une couche de type PostGIS, vous obtenez en effet une message d'erreur du type :

Warning: [MapServer Error]: msDrawMap(): Failed to draw layer named 'layer1'. in D:\MapJaxDist\ms4w\apps\mapjax\interface\gestion\update\index.php on line 297

Warning: [MapServer Error]: prepare_database(): Error declaring cursor: ERROR: current transaction is aborted, commands ignored until end of transaction block With query string: DECLARE mycursor BINARY CURSOR FOR SELECT asbinary(force_collection(force_2d(wkb_geometry)),'NDR'),OID::text from septolis_departements_l93 WHERE wkb_geometry && setSRID('BOX3D(-109363.162923348 5996601,1450962.16291935 7163567)'::BOX3D, find_srid('','septolis_departements_l93','wkb_geometry') ) in D:\MapJaxDist\ms4w\apps\mapjax\interface\gestion\update\index.php on line 297

Warning: [MapServer Error]: msPOSTGISLayerRetrievePGVersion(): Error executing POSTGIS statement (msPOSTGISLayerRetrievePGVersion():select substring(version() from 12 for (position('on' in version()) - 13)) in D:\MapJaxDist\ms4w\apps\mapjax\interface\gestion\update\index.php on line 297

Fatal error: Call to a member function saveImage() on a non-object in D:\MapJaxDist\ms4w\apps\mapjax\interface\gestion\update\index.php on line 298

La partie importante de ces messages d'erreur mentionne l'erreur MapServer dans la fonction msPOSTGISLayerRetrievePGVersion. La requête utilisée pour récupérer le numéro de version de PostgreSQL ne fonctionne pas. Si on regarde de plus prêt, rien de très compliqué ici, on utilise la chaine retournée par la fontion version() de PostgreSQL.

Dans notre cas, la fonction version() renvoit : PostgreSQL 8.3.7, compiled by Visual C++ build 1400. Ainsi il n'est pas possible d'utiliser la position de la chaîne 'on' dans ce résultat étant donnée qu'elle n'y est pas.

La solution la plus simple est d'utiliser le script SQL ci-dessous qui se contente de renommer la fonction version() intiallement présente dans PostgreSQL en initial_version_func puis de déclarer une fonction version qui renvoie une chaîne correcte pour ce cas d'utilisation :

; Renommage de la fonction version
ALTER FUNCTION version() RENAME TO initial_version_func;

; Définition de la nouvelle fonction version
CREATE OR REPLACE FUNCTION version() RETURNS VARCHAR AS
$$
DECLARE
   vers text;
   rec record;
BEGIN
   FOR rec in EXECUTE $q$ SELECT 'PostgreSQL 8.3.7 on win32, compiled Visual C++ build 1400' as version$q$
   LOOP
     RETURN rec.version;
   END LOOP;
   RETURN '';
END
$$ LANGUAGE plpgsql;

Le problème est ainsi corrigé.

Il est à noté que, grâce aux corrections de Paul Ramsey, cette erreur ne devrait plus être présente dans les versions de MapServer sorties depuis le 11 avril 08 (cf. changeset 8003).

MapServer 5.4.0 affichage de symboles aux extrémintés d'une ligne

Tue, 28 Apr 2009 16:41:27 +0200

Avec MapServer 5.4.0, l'affichage des extrémités des couches linéaires est maintenant facilité et possible à l'aide d'une seule couche.

Plus besoin donc de pré-traitement pour récupérer le début et la fin de la ligne requêtes suivantes, qui se faisait néanmoins très facilement avec des couches PostGIS. En utilisant, par exemple, les définitions de couches (1 et 2) ci-dessous ( la définition de la couche ligne de base est omise), respectivement pour le point de départ et le point final (on suppose que la table est <table>):


Définition de couche 1 :


 LAYER

     CONNECTION ""

     CONNECTIONTYPE POSTGIS
    DATA "wkb_geometry1 from (SELECT startPoint( <table>.wkb_geometry) as wkb_geometry1,<table>.ogc_fid, CASE WHEN plus(pow((y(e)-y(s)),2),pow((x(e)-x(s)),2)) = 0 THEN 1 ELSE (180*(acos((x(e)-x(s))/sqrt(plus(pow((y(e)-y(s)),2),pow((x(e)-x(s)),2))))/pi()))-90 END as angle FROM (select ogc_fid, pointN(wkb_geometry,2) as e, startPoint(wkb_geometry) as s from <table>) as foo,  <table> WHERE foo.ogc_fid=<table>.ogc_fid)  as foo1 using unique ogc_fid using SRID=27572"
     NAME "espt_start"

     PROCESSING "CLOSE_CONNECTION=DEFER"

     STATUS OFF

     TOLERANCE 0

     TYPE POINT

     UNITS METERS

     CLASS

         NAME "espt_start1"

         STYLE

             ANGLE [angle]

             COLOR 250 8 255

             SIZE 10

             SYMBOL "triangle_isocele"

             WIDTH 1

         END

     END

  END 

Définition de couche 2 :


 LAYER

     CONNECTION ""

     CONNECTIONTYPE POSTGIS
    DATA "wkb_geometry1 from (SELECT endPoint( 	<table>.wkb_geometry) as wkb_geometry1,<table>.ogc_fid, CASE WHEN plus(pow((y(e)-y(s)),2),pow((x(e)-x(s)),2)) = 0 THEN 1 ELSE (180*(acos((x(e)-x(s))/sqrt(plus(pow((y(e)-y(s)),2),pow((x(e)-x(s)),2))))/pi()))-270 END as angle FROM (select ogc_fid, pointN(wkb_geometry,numPoints(wkb_geometry)-1) as s, endPoint(wkb_geometry) as e from <table>) as foo,  poisson.total_especes2 WHERE foo.ogc_fid=<table>.ogc_fid)  as foo1 using unique ogc_fid using SRID=27572"
     NAME "espt_end"

     PROCESSING "CLOSE_CONNECTION=DEFER"

     STATUS OFF

     TOLERANCE 0

     TYPE POINT

     UNITS METERS

     CLASS

         NAME "espt_end1"

         STYLE

             ANGLE [angle]

             COLOR 250 8 255

             SIZE 10

             SYMBOL "triangle_isocele"

             WIDTH 1

     END

         END

     END



L'astuce consistait donc simplement à ajouter deux couches supplémentaires, l'une pour le point de départ et l'autre pour le point final (si votre symbol est orienté) en plus de la couche de la ligne elle-même. Avec la nouvelle version de MapServer 5.4.0, on peut utiliser l'option GEOMTRANSFORM d'une balise STYLE qui, comme son nom l'indique, permet de réaliser une opération sur une géométrie (end, start, buffer ...). Par exemple, si vous souhaitez afficher une barre perpendiculaire aux extrémités d'une ligne, voilà la définition d'une couche MapServer qui utilise ce nouveau paramètre :

 LAYER

      CONNECTION ""

     CONNECTIONTYPE POSTGIS

     DATA "wkb_geometry FROM <table>"

     NAME "espt"

     PROCESSING "CLOSE_CONNECTION=DEFER"

     STATUS OFF

     TOLERANCE 0

     TYPE LINE

     UNITS METERS

     CLASS

         NAME "1"

         STYLE

         ANGLE 360

         COLOR 250 8 255

         SYMBOL 0

         WIDTH 2

     END

     STYLE

         ANGLE 360

         ANTIALIAS TRUE

         COLOR 0 0 255

         SIZE 5

         SYMBOL "l_oant_01_line"

         WIDTH 2

     END

     STYLE

         ANGLE AUTO

         ANTIALIAS TRUE

         COLOR 250 8 255

         MAXSIZE 100

         SIZE 10

         SYMBOL "openoffice_137"

         WIDTH 1

         GEOMTRANSFORM "end([wkb_geometry])"

     END

     STYLE

         ANGLE AUTO

         ANTIALIAS TRUE

         COLOR 250 8 255

         MAXSIZE 100

         SIZE 10

         SYMBOL "openoffice_137"

         WIDTH 1

         GEOMTRANSFORM "start([wkb_geometry])"

     END

        END

  END 


Il est donc maintenant possible de facilement afficher un symbole quelconque au début, à la fin ou pour les plus gourmands les deux 
Attention, lors d'une interrogation de données, vous constaterez sans doute le même résultat que l'image ci-dessous, le trait horizontal est répété tout au long de la ligne :



Lors de la sauvegarde des mapfiles à l'aide de l'API mapscript avec cette nouvelle version, j'ai rencontré un problème pour les couches ayant un STYLE utilisant ANGLE AUTO. Cette partie de la définition de la couche était  remplacée par ANGLE 360. Du coup, j'ai créé le "patch" qui suit (applicable à l'aide de la commande patch -P0 mapserver.patch depuis le répertoire des sources de MapServer) pour corriger ce problème :

 

diff -ru ./mapfile.c ./mapfile.c

--- ./mapfile.c    2009-04-07 22:30:23.000000000 +0200

+++ ./mapfile.c 2009-04-28 23:52:43.000000000 +0200

@@ -1911,7 +1911,7 @@

   style->minscaledenom=style->maxscaledenom = -1.0;

   style->offsetx = style->offsety = 0; /* no offset */

   style->antialias = MS_FALSE;

-  style->angle = 360;

+  style->angle = 0;

   style->autoangle= MS_FALSE;
   style->opacity = 100; /* fully opaque */

   style->_geomtransformexpression = NULL;

@@ -2127,7 +2127,9 @@

   fprintf(stream, "      STYLE\n");

   if(style->numbindings > 0 && style->bindings[MS_STYLE_BINDING_ANGLE].item)

      fprintf(stream, "        ANGLE [%s]\n", style->bindings[MS_STYLE_BINDING_ANGLE].item);

-  else if(style->angle != 0) fprintf(stream, "        ANGLE %g\n", style->angle);

+  else if(style->angle != 0 && style->autoangle!=MS_TRUE)

+     fprintf(stream, "        ANGLE %g\n", style->angle);

+  else if(style->autoangle==MS_TRUE) fprintf(stream, "        ANGLE AUTO \n");
   if(style->antialias) fprintf(stream, "        ANTIALIAS TRUE\n");

   writeColor(&(style->backgroundcolor), stream, "BACKGROUNDCOLOR", "        ");

On utilise donc la propriété autoangle de style pour vérifier si l'angle doit être ajouté ou non. Au passage, vous noterez que l'on n'ajoute plus un angle de 360 par défaut si aucun angle n'est spécifié pour la couche en question (ref. #2991).
Un autre détail, en passant, concernant cette nouvelle version 5.4.0 de MapServer : la nécessité d'avoir le mot clef SYMBOLSET au début du fichier symbols.sym (la balise END n'est pas requise).



Une nouvelle version de PROJ.4 (4.6.1)
Fri, 05 Sep 2008 11:57:42 +0200
Le 1er septembre 2008, Frank Warmerdam annonçait qu'une nouvelle version de PROJ été sortie : PROJ-4.6.1.
Au programme : 

Mise à jour de nad/epsg à la version EPSG 6.17. Correction du problème de précision introduit dans la version précédente.
Ajout de la logique pour la projection tmerc afin qu'elle échoue plutôt que renvoiyer des résultats éronés dans des zones éloigné de plus de 90° du méridien centre (#5). Cette modification est temporaire et restera le temps qu'une meilleure solution soit trouvée.
Correction de gestion de texte supplémentaire dans cs2cs
Renommage de INSTALL en INSTALL.txt
La projection egc a été généralisée pour inclure la latitude de l'origine.
Ajout de la projection glabsgm (Gauss Laborde / Sphere Geometric Mean), et sa variation gstmerc.
Création du catalogue d'initialisation nad/IGNF.
Ajout du fichier ntf_r93.gsb
Ajout de l'option de compilation /Op dans nmake.opt pour éviter le bug d'optimisation de VC7  (#12)
Correction de plusieurs scripts de test qui fonctionneront lorsque les fichiers grid seront disponibles.



Fonction PL/PGSQL permettant la génération de GPolyline
Sat, 10 May 2008 17:39:09 +0200
Vous trouverez ici le code d'une fonction PL/PGSQL permettant de produire des fichiers XML utilisable depuis l'API Google Map pour produire des GPolyline.


Mapserver 5.0.0 sous Gentoo
Fri, 21 Sep 2007 10:24:59 +0200
Le FOSS4G 2007 s'approchant, vous aurez sans doute déjà remarqué que la version de mapserver 5.0.0 avec le support de la librairie agg est sortie. Cette version est dorénavent disponible pour la ditribution Gentoo Linux depuis ce matin.
Nous en avons profiter pour mettre à jour note système de génération de document statistiques utilisant un serveur OpenOffice sur le système geolive. Nous avons veillé à conserver deux versions, l'une utilisant le format de sortie AGG et l'autre utlisant GD.
Je vais me rendre au FOSS4G 2007 au Canada et je compte bien, de là-bas, en profiter pour mettre à jour de nombreux ebuilds pour la distribution Gentoo Linux et plus particulièrement ceux de PostGIS.
En espérant croiser certains d'entre vous à Victoria.


Mapserver 4.10.3 sous gentoo
Thu, 06 Sep 2007 10:35:37 +0200
     
La dernière version stable de mapserver de la branche 4.10, à savoir la version 4.10.3, corrigeant entre autre une faille de sécurité est maintenant disponible pour la distribution Gentoo Linux.
Le patch permettant de compiler à la fois le support pour php4 et php5 a été mis à jour pour suivre les modifications apportées au fichier configure.in du projet. Il est disponible en libre téléchargement, pour de plus amples informations sur ce patch vous pouvez vous reporter à la page décrivant l'ancienne version de ce patch.
Je vous invite donc à tester ce nouvel ebuild et le cas échéant à remplir une demande de correction de bogues à l'aide de cette interface.


La version 1.3.1 de PostGIS
Tue, 14 Aug 2007 09:55:41 +0200
La version 1.3.1 de PostGIS est maintenant disponible. Cette version introduit de nouvelles fonctionnalités et apporte diverses améliorations concernant les performances.

Début de la migration des noms de fonction pour respecter la convention de nommage de la  SQL-MM, utilisant le type (SP) pour prefix.
Création de racourcis pour "point-in-polygon" mis en cache et indexés pour les fonctions ST_Contains, ST_Intersects, ST_Within, ST_Dijoint.
Ajout de l'utilisation automatique des indexes pour les fonctions relationnelles (cf. cette note en bas de page).
Extension du support du type géométrique curve : accesseurs et certaines fonctions de traitement.
Ajout du dialecte "Hibernate"

Version anglaise de la nouvelle.


Sur la route vers mapserver 5.0
Tue, 07 Aug 2007 10:59:17 +0200
La nouvelle version 5.0 de MapServer, qui devrait sortir le 12 septembre, un peu avant le FOSS4G2007 donc, offrira de nouvelles fonctionnalités comme par exemple la création de graphiques statistiques dynamiques sur les carte. Comme le dit très bien Howard Butler, il en sera enfin fini des phrases du genre "ha, ça c'est une image générée par mapserver !". En effet, comme vous le savez sans doute, les versions actuelles de MapServer utilisent la librairie gd pour rendre les images. Cette librairie n'est plus maintenue par les initiateurs du projet depuis maintenant un petit moment, bien que des efforts aient été menés afin de palier à ce problème. La nouvelle version de MapServer utilisera donc une autre librairie : AGG, qui semble bien plus efficace et rend des images d'une meilleure qualité. Regardez par exemple une sortie utilisant gd et la sorite équivalente utilisant agg. Il faut d'ailleurs savoir que l'utilisation de la librairie agg en lieu et place de la librairie gd est déjà disponible sur la version svn de MapServer.


Connecteur PostGIS pour FDO
Mon, 30 Apr 2007 11:02:25 +0200
Il existe une version alpha du connecteur PostGIS pour FDO (mise à disposition afin d'effectuer des tests uniquement).
Actuellement, ce connecteur peut être utilisé avec la version 8.2.X de PostgreSQL et PostGIS version 1.1.6 à 1.2.1.
Vous pouvez télécharger les binaires (fdopostgis-20070428-r390.zip) ici


Postgis 1.1.6 et 1.2.0 sous gentoo
Wed, 10 Jan 2007 13:31:43 +0100
  
Les ebuilds pour les versions 1.1.6 et 1.2.0 ont été ajoutés à l'arborescence officielle de portage le 6 janvier. Je tiens d'ailleurs à remercier le dévelopeur gentoo Tiziano Müller (aka dev-zero) pour ses remarques constructive lors de la réalisation de ces deux ebuilds.


Problème de "mise à jour légère" pour passer à la version 1.2.0
Tue, 19 Dec 2006 20:12:57 +0100
Si vous avez installé la version 1.2.0 de PostGIS et que vous possédiez déjà des bases de données spatiales utilisant une version antérieur, vous avez sans doute rencontré des problèmes lors de la première tentative de "mise à jour légère" de votre base.


Ebuild de mapserver.
Mon, 27 Nov 2006 11:28:42 +0100
 
Très tôt dans la matinée du samedi 25 novembre, je me suis enfin décider à poster l'ebuild de mapserver dans l'arborescence officielle du portage de gentoo.
L'ensemble des supports disponibles, à l'exception de celui pour csharp, sont activables via les drapeaux USE dont vous avez l'habitude. Le support pour SOS a de plus été ajouté pour cette nouvelle version de l'ebuild et est activé à l'aide des deux drapeaux proj et xml. Comme c'était déjà le cas dans les versions précédantes, l'ebuild permet de compiler le support phpmapscript pour les diverses versions de php que vous avez installé, vous pouvez en effet installer le support à la fois pour php4 et php5 (vous trouverez de plus amples informations en anglais ici concernant la configuration d'apache dans ce cas de figure).


Un nouvel ebuild pour proj 4.5.0
Sat, 25 Nov 2006 00:02:29 +0100
 
La toute denière version de proj 4.5.0 est maintenant disponible pour les heureux utilisateurs de Gentoo. Les anciennes versions (la 4.7.7-r1 et la 4.7.8 ) qui tranaient encore dans l'arborescence officelle de gentoo ont été par la même occasion supprimées. Merci à Tom Payne pour ses conseils et remarques  que vous pouvez consulter ici.
Je suis de plus devenu le responsable offciel de ce paquet pour le projet Gentoo Linux, comme vous pouvez le constatez en consultant les métadonnées du paquet.