geotools工具

# # v0.0 : DGR, le 07 fevrier 2004 # v0.1 : DGR, le 08 fevrier 2004 mise a jour des controles (4326, 2154) # v0.2 : DGR, le 09 fevrier 2004 mise a jour des controles (7412 --correction boggue geotools2--) # # Installation de la base EPSG (http://www.epsg.org/) sous postgreSQL Hypothese : 1.- OS Linux (installe, en fonctionnement) 2.- installation Mandrake de postgreSQL standard. 3.- postGIS installe dans contrib de pgsql (/usr/lib/pgsql). 4.- scripts SQL fournis par l'EPSG : 4.1.- EPSG_v65.mdb_DDL.sql (dictionnaire de la base) ; 4.2.- EPSG_v65.mdb_DML.sql (donnees). 5.- fonctions en jeu : 5.1.- administrateur systeme (AS) ; 5.2.- administrateur postgreSQL : c'est lui qui lance le serveur postgreSQL (il est normalement non loggable) ; 5.3.- un utilisateur Unix, administrateur de la base a creer (US) . But : a partir des scripts fournis par l'EPSG, re-creer la base EPSG des systemes de references de coordonnees. Note : les scripts SQL ont ete debarasses des retour-chariots Windows, ainsi que des espaces en fin de ligne. Note : EPSG_v65.mdb_DDL.sql : #302, manque une virgule entre parameter_code et coord_op_method_code un begin; et un commit; ont ete rajoute (debut et fin) les contraintes (ALTER TABLE) ont ete deplacees dans EPSG_v65.mdb_DDL2.sql Note : EPSG_v65.mdb_DML.sql : l'ordre des tables n'est pas conforme aux contraintes d'integrite, un bon ordre est : epsg_namingsystem epsg_alias epsg_area epsg_change epsg_deprecation epsg_coordinateaxisname epsg_coordinatesystem epsg_unitofmeasure epsg_coordinateaxis epsg_ellipsoid epsg_primemeridian epsg_datum epsg_coordinatereferencesystem epsg_coordoperationmethod epsg_coordoperation epsg_coordoperationpath epsg_coordoperationparam epsg_coordoperationparamusage epsg_coordoperationparamvalue epsg_versionhistory Des begin; ont ete rajoutes avant les commit; Note : u script SQL EPSG_v65.mdb_DDL3.sql a ete cree pour detruire ces tables ! Methode : 1.- Creation des acces a la future base (AS) : Modification du fichier de configuration des autorisations d'acces /var/lib/pgsql/data/hb_pga.conf, (AS) ajoute : ---- ... host epsg 127.0.0.1 255.255.255.255 password web host epsg 192.168.0.0 255.255.255.0 password web ---- Elles indiquent que la base 'epsg' est accessible localement ou par la machine d'adresse IP 192.168.0.0, sous reserve de fournir un mot de passe pour un utilisation dont le nom existe dans le fichier /var/lib/pgsql/data/web : ---- ... epsg_reader ---- (AS) relance alors le serveur postgreSQL : $ service postgresql restart Arr�t de postgresql [ OK ] D�marrage de postgresql [ OK ] 2.- Creation de l'administrateur de la future base (AS) : (AS) cree l'utilisateur gestionnaire de la future base si l'utilisateur Unix postgres n'est pas loggable (sinon cette operation s'effectue en tant que postgres) : $ su - postgres -c "createuser --createdb --adduser --pwprompt --echo epsg_reader" Enter password for user "epsg_reader": Enter it again: CREATE USER "epsg_reader" WITH PASSWORD '#epsg' CREATEDB CREATEUSER CREATE USER L'administrateur de la future base postgreSQL est cree. 2.- Creation de la base (US) : L'utilisateur (administrateur) commence par creer la base postgreSQL : $ createdb -U epsg_reader -h localhost -E LATIN9 -e epsg "EPSG sous postgreSQL" Password: CREATE DATABASE "epsg" WITH ENCODING = 'LATIN9' CREATE DATABASE Password: COMMENT ON DATABASE "epsg" IS 'EPSG sous postgreSQL' COMMENT L'administrateur permet le langage PLpgSQL pour postGIS : $ createlang -U epsg_reader -h localhost plpgsql epsg Password: Password: Password: Password: Puis, il importe les tables associees a postGIS : $ psql -U epsg_reader -h localhost -f /usr/lib/pgsql/contrib/postgis/postgis.sql epsg Password: BEGIN ... COMMIT $ psql -U epsg_reader -h localhost -f /usr/lib/pgsql/contrib/postgis/spatial_ref_sys.sql epsg Password: BEGIN ... COMMIT Enfin, il se donne les acces a ces tables : $ psql -U epsg_reader -h localhost -c "grant select on geometry_columns, spatial_ref_sys to epsg_reader;" epsg Password: GRANT Note : pourquoi creer une base geo-localisee pour gerer des systemes de coordonnees ? L'idee est de pouvoir "verifier" les contraintes spatiales associees a un systeme et, eventuellement, de graphiquement les representer !-) 3.- Insertion des donnes de l'EPSG (US) : L'administrateur cree d'abord le dictionnaire EPSG : $ psql -U epsg_reader -h localhost -f /path/2/EPSG_v65.mdb_DDL.sql epsg Password: ... Puis, il insere les donnes : $ psql -U epsg_reader -h localhost -f /path/2/EPSG_v65.mdb_DML.sql epsg ... Puis, il lance les contraintes : $ psql -U epsg_reader -h localhost -f /path/2/EPSG_v65.mdb_DDL2.sql epsg ... 4.- Controle de la base (US) : $ psql -U epsg_reader -h localhost epsg Password: Welcome to psql, the PostgreSQL interactive terminal. Type: \copyright for distribution terms \h for help with SQL commands \? for help on internal slash commands \g or terminate with semicolon to execute query \q to quit epsg=# \dt List of relations Name | Type | Owner --------------------------------+-------+------------- epsg_alias | table | epsg_reader epsg_area | table | epsg_reader epsg_change | table | epsg_reader epsg_coordinateaxis | table | epsg_reader epsg_coordinateaxisname | table | epsg_reader epsg_coordinatereferencesystem | table | epsg_reader epsg_coordinatesystem | table | epsg_reader epsg_coordoperation | table | epsg_reader epsg_coordoperationmethod | table | epsg_reader epsg_coordoperationparam | table | epsg_reader epsg_coordoperationparamusage | table | epsg_reader epsg_coordoperationparamvalue | table | epsg_reader epsg_coordoperationpath | table | epsg_reader epsg_datum | table | epsg_reader epsg_deprecation | table | epsg_reader epsg_ellipsoid | table | epsg_reader epsg_namingsystem | table | epsg_reader epsg_primemeridian | table | epsg_reader epsg_unitofmeasure | table | epsg_reader epsg_versionhistory | table | epsg_reader geometry_columns | table | epsg_reader spatial_ref_sys | table | epsg_reader (22 rows) epsg=# \q Si la bibliotheque geotools2 est installe, un controle plus fin sera : $ java -cp /path/2/opengis-2.0-B2-DEV.jar:\ /path/2/resources-2.0-B2-DEV.jar:\ /path/2/cts-coordtrans-2.0-B2-DEV.jar:\ /path/2/postgresql-3.0.jar:$CLASSPATH \ org.geotools.cs.CoordinateSystemEPSGFactory \ -driver=org.postgresql.Driver \ -connection="jdbc:postgresql://localhost:5432/epsg?user=epsg_reader&password=#epsg" \ -implementation=org.geotools.cs.CoordinateSystemModifiedEPSGFactory \ 4326 2154 7412 <=== EPSG 4326 ===> GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"] ], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["�0,017 rad",0.017453292519943278], AXIS["Geodetic latitude",NORTH], AXIS["Geodetic longitude",EAST], AUTHORITY["EPSG","4326"] ] <=== EPSG 2154 ===> P