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

Here I will try to show two examples of error with PostGIS and POLYGON.

1. Self-Intersection

The follwing POLYGON is not valid.

testgis=# SELECT IsValid('POLYGON((402 253,401 434,675 435,675 314,487 314,559 314,559 201,404 253,402 253))'::geometry);
INFO:  Self-intersection
INFO:  Self-intersection
 isvalid
---------
 f

He intersects hilmself. Have a look on the following picture!

2. Ring Self-Intersection

Suppose that you have a the following POLYGON

POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),(120 300, 300 280, 320 200, 
160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))

Let's try Summary() function.

testgis=# SELECT summary('POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),
(120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))'::geometry);
         summary
-------------------------

 Polygon[B] with 2 rings
    ring 0 has 5 points
    ring 1 has 10 points

(1 ligne)

With

• ring 0 = [380 340, 40 340, 40 20, 380 20, 380 340] --> My POLYGON
• ring 1 = [120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300] --> the hole of my polygon

If I try IsValid() function with PostGIS, I get

testgis=# SELECT IsValid('POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),
(120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))'::geometry);

INFO:  Ring Self-intersection
INFO:  Ring Self-intersection
 isvalid
---------
 f
(1 ligne)

That's not correct according to OGC Specification. A Polygon with a hole (ring 1) that self-intersects (i.e. POINT (320 200) )! Remember that it is just a case! This example has been taken from Geos SVN repository at ./tests/xmltester/TestValid.xml

RESUME

C'est sà»rement l'interface de programmation en C++ de PostgreSQL, la plus innovante. Elle constitue par rapport aux habituelles, l'une de la nouvelle génération. Son intérêt réside surtout avec son adéquation avec la bibliothèque STL (Standard Template Library). Comme le souligne son auteur, elle a été pensée «STL-compatible» . Le site de développement est http://thaiopensource.org/development/libpqxx/.

Autre point à  signaler, libpqxx est basée sur le mode transactionnel de PostgreSQL. Elle sort donc du cadre purement basique de la programmation «Requêtes - Réponses» pour l'interfaçage client-serveur. La richesse de ses fonctionnalités est vraiment impressionnante, déroutante parfois au premier abord. Elle permet un code lisible, bien stylé, appréciable.

INSTALLATION SOUS UBUNTU DAPPER

En pré-requis, il faut avoir un serveur PostgreSQL sous la main. Libpqxx est basée sur libpq. Sous Windows, j'ai pu la compilé avec MinGW/Msys sans souci. Je livre ici mes notes pour une compilation sous Ubuntu Dapper (ou Debian). Je ne l'ai pas testé sous Visual C++.

Son installation se fait en faisant

wget http://thaiopensource.org/download/software/libpqxx/libpqxx-2.6.9.tar.gz
tar xvzf libpqxx-2.6.9.tar.gz
cd libpqxx-2.6.9

Chez moi PostgreSQL est installé à  /opt/pgsql/, je dois effectuer une modification du fichier libpqxx.pc.in avant de continuer. Je fais celà  surtout pour pouvoir compiler mes programmes par la suite. On ouvre donc le fichier et on effectue le remplacement suivant

prefix=@prefix@
exec_prefix=@exec_prefix@
libdir=@libdir@
includedir=@includedir@

Name: libpqxx
Description: a C++ API to the PostgreSQL database management system.
Version: @VERSION@
Libs: -L${libdir} -L@with_postgres_lib@ -lpq -lpqxx
Cflags: -I${includedir} -I@with_postgres_include@

On peut maintenant compiler.

export PATH=/opt/pgsql/bin/:$PATH
./configure 
make
make install

En essayant d'utiliser pqxx-config, j'ai eu l'erreur suivante

#pqxx-config --libs
pqxx-config is deprecated; please use "/usr/bin/pkg-config --libs libpqxx"

En fait c'est surtout sur la commande pkg-config --libs libpqxx ou pkg-config --cflags libpqxx qu'il faut se baser. Par exemple, pour l'une des deux j'obtiens

# pkg-config --libs libpqxx
-L/usr/local/lib -L/opt/pgsql/lib -lpq -lpqxx

TEST DE REGRESSION

Toujours dans les sources, il faut se rendre dans le sous-répertoire test les sources et faire

cd test
chmod 777 -R .
su postgres
make check

UN PREMIER PROGRAMME

Pour la compilation de mon programme, je vais me doter du Makefile suivant

PKGCONFIG=/usr/bin/pkg-config
CXXFLAGS = -g -Wall $(shell $(PKGCONFIG) --cflags libpqxx)
LDLIBS   = $(shell $(PKGCONFIG) --libs libpqxx)

exes = test1

all:	$(exes)

test1: test1.cc

clean:
	rm -f $(exes)

Si l'utilisation d'un Makefile peut s'avérer déroutant, on pourra essayer

g++ -g -Wall -I/usr/local/include -I/opt/pgsql/include      test1.cc  -L/usr/local/lib -L/opt/pgsql/lib -lpq -lpqxx   -o test1

Mon programme test1.cc sera le suivant, juste un petit programme de connexion

/* test1.cc */

#include <stdlib.h>	// Nécessaire pour exit(), EXIT_SUCCESS, et EXIT_FAILURE
#include <iostream> 	// requis pour cerr
#include <pqxx/pqxx>	// définitions pour libpqxx

using namespace pqxx;
using namespace std;

int main( int argc, const char * argv[] )
{

  try 
  {
    connection    myConnection( argc > 1 ? argv[1] : "" );

    myConnection.activate();
    if (myConnection.is_open())
    {
       string host = to_string( myConnection.server_version() );
       host.replace(1,1,"."); host.replace(3,1,".");
       cout<<" INFORMATIONS:"<<endl;
        cout<<" Paramètres de connexion host="<<myConnection.hostname();
       cout<<" dbname="<<myConnection.dbname()<<" user="<<myConnection.username()<<endl;
       cout<<" Serveur - Version "<<host<<endl;
    }
  }	
  catch( runtime_error & e )
  {
    cerr << "Echec à  la connexion avec erreur := " << e.what();
    exit( EXIT_FAILURE );
  }
  catch( exception & e )
  {
    cerr << e.what();
    exit( EXIT_FAILURE );
  }
  catch( ... )
  {
    cerr << "Erreur/exception inconnue !" << endl;
    exit( EXIT_FAILURE );
  }

  exit( EXIT_SUCCESS );

}

Ensuite on tape

make

Pour une utilisation standard, je vais passer à  test1 les paramètres de connexion

./test1 "host=localhost dbname=testgis user=postgres password=empr888"

qui me renverra

 INFORMATIONS:
 Paramètres de connexion host=localhost dbname=testgis user=postgres
 Serveur - Version 8.2.3

On peut aussi utiliser des variables d'environnement de PostgreSQL, par exemple

# PGHOST=localhost PGUSER=postgres ./test1 "dbname=testgis password=empr888"
 INFORMATIONS:
 Paramètres de connexion host=localhost dbname=testgis user=postgres
 Serveur - Version 8.2.3

Bien sà»r, en cas d'échec de connexion, j'aurais par exemple sur un hà´te distant

# PGHOST=192.168.0.5 PGUSER=postgres ./test1 "dbname=testgis password=empr888"
Echec à  la connexion avec erreur := FATAL:  authentification par mot de passe échouée pour l'utilisateur  «postgres»

Suite à  mon précédent billet, je me suis amusé à  tester les données NavTeq de ADCI avec le module PgRoutin de KaMap. Un exemple de démo avec les données de Orkney est visible à  kappa.gonfi.org/dev/gonfi-ka-map/htdocs/index_routing.html

J'ai réussi à  obtenir par exemple la screenshot suivante pour mes propres données

Une démo est disponible chez moi à  http://www.davidgis.fr/ka-map/index_routing.html

Pour commencer il faut suivre les indications de http://ka-map.ominiverdi.org/wiki/index.php/KaRouting. Malheureusement il s'avère que les indications fournies sont incomplètes . Reprenons ici les détails de l'installation

On s'assure par exemple d'avoir le jeu de tests de PgRouting provenant de http://www.postlbs.org/postlbs-cms/en/download à  savoir «pgRouting Sample Application». Chez moi, j'ai décompressé le répertoire dans mon "/var/www/". J'ai ensuite copié tout ce qui était associé à  streets.shp dans /var/www/routingj/data/

cp streets.* /var/www/routingj/data/

J'ai ensuite crée la mapfile suivante /var/www/routingj/maps/kamap-streets.map dont le contenu est

MAP
  IMAGETYPE   PNG
  EXTENT      -77.1166763305664 38.7924194335938 -76.909538269043 38.9952430725098
  #SIZE        500 489.585635359116
  SIZE 200 195.834254143646
  IMAGECOLOR  255 255 255
  SHAPEPATH   "../data"

  FONTSET     "../fonts/fonts.list"
  
REFERENCE
   IMAGE "/var/www/keymap-streets.png"
  EXTENT -77.1166763305664 38.7924194335938 -76.909538269043 38.9952430725098
  STATUS ON
  COLOR -1 -1 -1
  OUTLINECOLOR 255 0 0
  SIZE 200 164
END

  WEB
    IMAGEPATH "/var/www/tmp/"
    IMAGEURL  "/tmp/"
  END
		
  PROJECTION
   "init=epsg:4326" # JGD2000
  END

  UNITS DD

  SYMBOL
    NAME "hatch"
    TYPE hatch
  END
  
  LAYER 
    NAME "roads"
    DATA streets
    STATUS DEFAULT
    TYPE LINE

    CLASS
      COLOR 155 155 155
    END

  END 

END 

Bien sà»r pour calculer l'extent et les dimensions adéquates de l'image pour width=500, je me suis servi des requêtes suivantes

routing=# select extent(the_geom) from streets_edges;
                                  extent
---------------------------------------------------------------------------
 BOX(-77.1166763305664 38.7924194335938,-76.909538269043 38.9952430725098)
(1 ligne)

et de

routing=# SELECT
               500*(abs(Ymax(foo.extent)-Ymin(foo.extent))/abs(Xmax(foo.extent)-Xmin(foo.extent)))
               AS hauteur
         FROM (SELECT Extent(the_geom) FROM streets_edges) AS foo;
     hauteur
------------------
 489.585635359116
(1 ligne)

Générer la keymap pour la légende

Pour connaà®tre les dimensions pour width=200

routing=# SELECT
               200*(abs(Ymax(foo.extent)-Ymin(foo.extent))/abs(Xmax(foo.extent)-Xmin(foo.extent)))
               AS hauteur
         FROM (SELECT Extent(the_geom) FROM streets_edges) AS foo;
     hauteur
------------------
 195.834254143646
(1 ligne)

Le script php suivant a amplement suffit pour générer la keymap /var/www/keymap-streets.png

<html>
<body>
<?php
    $sw_MapFile = "/var/www/routingj/maps/kamap-streets.map";
    
    $map = ms_newMapObj( $sw_MapFile );
    
    $image = $map->draw();
//$image_url = $image->saveWebImage(MS_PNG,1,1,0);
$image->saveImage("/var/www/keymap-streets.png");
echo "<IMG   
BORDER=0 
SRC='".$image_url."'
width='".$map->width."' height='".$map->height."'/>
<BR>"; 
?>
</body>
</html>

Configuration de PgRouting avec Kamap

On peut télécharger la CVS disponible à  http://ka-map.ominiverdi.org/CVSimg/ka-map-cvs-image.tar.gz

cd /var/www/
wget http://ka-map.ominiverdi.org/CVSimg/ka-map-cvs-image.tar.gz
tar xvzf ka-map-cvs-image.tar.gz

Ce qui donnera le répertoire ka-map.

Comme précisé dans le wiki

i
cd ka-map
cp include/config.routing.dist.php include/config.php

Par rapport au contenu proposé voici le mien avec les différences notables en gras

 // $aszGMap = array (
 //        'title' => 'GMap 75',
 //        'path' => '/var/www/gmap/htdocs/gmap75.map',
 //        'scales' => array( 40000000, 25000000, 12000000, 7500000, 3000000, 1000000 ),
 //        'format' =>'PNG'
 //        /* Sample authorized_users entry. See auth.php for more details:
 //         * ,'authorized_users' => array('popplace' => array('user1', 'user2'),
 //         *                              'park'     => array('user1')
 //         */
 //);

 $aszKaRoute = array (
         'title' => 'Routing',
         'path' => '/var/www/routingj/maps/kamap-streets.map',
         'scales' => array(100000,50000,10000 ),
         'format' =>'PNG'
         /* Sample authorized_users entry. See auth.php for more details:
          * ,'authorized_users' => array('popplace' => array('user1', 'user2'),
          *                              'park'     => array('user1')
          */
 );
 $aszMapFiles = array( 'routeMap'=>$aszKaRoute);


/******************************************************************************
 * figure out which map file to use and set up the necessary variables for
 * the rest of the code to use.  This does need to be done on every page load
 * unfortunately.
 *
 * szMap should be set to the default map file to use but can change if
 * this script is called with map=.
 */
$szMap = 'routeMap';

/******************************************************************************
 * Configure kaRouting settings
 * the array set the second mapfile you need to call and let kaRouting.php 
 * override and the postgres connection string
 */
$kaRoute = Array(
			'routeMap' => '/var/www/routingj/maps/mapfile.map',//the second mapfile -should be changed with an SLD file
			'pgConn' => 'dbname=routing host=localhost port=5432 user=postgres password=empr888',
			'myGraphTable' => 'streets_edges',
			'costColumn' => 'length',
			'reverseCostColumn' => 'reverse_cost',
			'directed' => 'false',//check pgRouting docs for this value
			'has_reverse_cost' => 'true',//check pgRouting docs for this value
			'buffer' => '0.01'//should be related to your map unit
		);

La pour la suite, il faut créer la mapfile /var/www/routingj/maps/mapfile.map dont le contenu sera - conformément au wiki -

MAP
  IMAGETYPE   png
  EXTENT      139.602 35.413 139.678 35.489
  SIZE        550 450 

   IMAGECOLOR  255 255 255
  SHAPEPATH   "data" #there are no shapes 

   FONTSET      "../fonts/fonts.list"
  WEB
    IMAGEPATH "../tmp/"
    IMAGEURL  "/routingj/tmp/" 
  END 

  PROJECTION
   "init=epsg:4612" # JGD2000. the kanagawa one
  END 

  UNITS DD 

  SYMBOL
    NAME "hatch"
    TYPE hatch
  END 

LAYER
      NAME "path"
      CONNECTION "user=postgres password=empr888 dbname=routing host=localhost port=5432"
      CONNECTIONTYPE postgis
      STATUS OFF 
      TYPE LINE 

      CLASS
        NAME "path"
        STYLE
          COLOR 155 0 0
          WIDTH 2
       END
      END
 END


END

Rappel: Création de la table streets_edges

Par la suite, on s'assurera pour la création de la table streets_egdes de suivre les indications de http://www.davidgis.fr/documentation/win32/html/apa.html aux sections «A.8.1. Importation d'un jeu de données NavTeq» et «A.8.2. Noeuds du réseau et direction pour le routage».

Recommencdations supplémentaires du wiki

Pour être conforme au wiki, il faut ensuite modifier son fichier de configuration d'Apache

Alias /ka-map/ "/var/www/ka-map/htdocs/"

       Options Indexes
       AllowOverride None
       Order allow,deny
       Allow from all

Et ensuite de créer la table suivante dont le contenu sera

CREATE TABLE kroute (
    gid serial NOT NULL,
    sessionid text,
    the_geom geometry NOT NULL,
    date timestamp without time zone,
    CONSTRAINT enforce_dims_the_geom CHECK ((ndims(the_geom) = 2)),
    CONSTRAINT enforce_geotype_the_geom CHECK (((geometrytype(the_geom) = 'MULTILINESTRING'::text) OR (the_geom IS NULL))),
    CONSTRAINT enforce_srid_the_geom CHECK ((srid(the_geom) = -1))
);

Pour lancer l'application, il suffira de pointer vers http://localhost/ka-map/index_routing.html

Hier je suis allé le voir au Gaumont Multiplexe de Montpellier. Pour la séance de 16h30 - première semaine -, y'avait un de ces mondes ! J'suis arrivé, j'ai fait la queue. Les panneaux interactifs qui fournissaient le nombre de places restantes mentionnait 94 places restantes. Même pas deux minutes plus tard, il en restait 84.

Du coup dans la salle, j'ai été aux premièrs rangs. Mais celà  ne m'a pas emêché de matter le film . Que dire si ce n'est que c'était un bon film d'action. Moi qui a baigné dans l'univers de Marvel quand zétais petit: Iron-Man, X-Men. Comme troisième opus, c'était vraiment pas mal. C'est surtout Venom qui m'a fait craqué aussi bien fait que dans la BD avec ses nombreuses dents asserrées. Pour le reste, c'était plutà´t pas mal que d'avoir 3 super-vilains que d'un seul pour une fois.

Par contre avant le film, j'ai vu la bande-annonce de Transformers. De la balle, ça promet en juillet cet été. Un de mes dessins animés préférés porté à  l'écran. Je ne manquerais pas le RDV. Ce qui est dommage c'est que dans la bande-annonce, on ne devine pas les robots. C'est uniquement à  la fin en voyant le symbole suivant

Y'a aussi Pirates des Caraà¯bes 3 qui devrait bientà´t sortir. On aura droit à  de bons films cette année. Au moins j'aurais l'occasion de partir plus souvent au cinéma cet été.