GeotoolKit - Géoréférencement

Par défaut GeotoolKit tentera d'utiliser une base de données Derby/JavaDB à condition d'avoir les .jar : geotk-epsg et derby dans le projet.

Si la base de données est inexistante, celui-ci la créera automatiquement. Elle sera placée à l'endroit suivant :
Linux = /home/utilisateur/.geotoolkit.org/EPSG/version
Mac = /Users/utilisateur/Library/GeotoolKit.org/EPSG/version
Windows = C:\Documents and Settings\utilisateur\Application Data\Geotoolkit.org\EPSG\version
Cette base de données n'est accessible qu'à une seule application à la fois.
Donc si vous avez plus d'une application utilisant GeotoolKit en même temps il vous faudra utiliser une des solutions qui suivent.

De la même façon que pour Derby/JavaDB, GeotoolKit se chargera de créer automatiquement la base de données dans le même répertoire à condition d'avoir les .jar geotk-epsg et hsql dans le projet.
Attention, dans ce cas, à ne pas avoir le jar de derby présent auquel cas celui-ci sera utilisé en premier.

Il faut inclure les jars odbc et jdbc-odbc bridge dans votre projet.

1- Téléchargez la base Access EPSG depuis http://www.epsg.org et copiez-la où vous le souhaitez.
2- Ouvrez le panneau de contrôle ODBC de Windows.
3- Cliquez sur Ajouter et sélectionnez : Microsoft Access Driver (*.mdb).
4- Nommez la source : EPSG.
5- Cliquez sur sélectionner et indiquez votre base EPSG_version.mdb.
6- Créez ou modifiez le fichier : C:\Documents and Settings\utilisateur\Application Data\Geotoolkit.org\EPSG\DataSource.properties pour qu'il contienne :

#
# Connection parameters to the EPSG database.
#
URL = jdbc:odbc:EPSG

Cette solution est dédiée au serveur ou à tout réseau afin d'avoir une base commune ou quand de multiples connections simultanées sont nécessaires.
Il faut inclure les .jar postgresql et geotk-epsg dans le projet.

1- Créer une base de données postgresql vide. Les lignes suivantes créent un rôle « Geotoolkit » et une base de données « Referencing » en UTF-8 sur votre machine locale.
(À adapter selon vos besoins)

createuser --createdb Geotoolkit
createdb --username=Geotoolkit --host=localhost --encoding=UTF-8 Referencing "CRS definitions"
createlang plpgsql --dbname=Referencing --username=Geotoolkit

2- Lancer le jar geotk-epsg comme suit :

java -jar geotk-epsg.jar create jdbc:postgresql://localhost/Referencing Geotoolkit --schema=epsg

3- Nettoyage de la base :

vacuumdb --full --analyze --username Geotoolkit --host localhost Referencing
psql -c 'ALTER USER \"Geotoolkit\" WITH NOCREATEDB' Referencing

4- Créer ou modifier le fichier DataSource.properties situé : Linux = /home/utilisateur/.geotoolkit.org/EPSG/DataSource.properties
Mac = /Users/utilisateur/Library/GeotoolKit.org/EPSG/DataSource.properties
Windows = C:\Documents and Settings\utilisateur\Application Data\Geotoolkit.org\EPSG\DataSource.properties

#
# Connection parameters to the EPSG database.
#
URL  = jdbc:postgresql://localhost/Referencing
user = Geotoolkit
password = *****

Pour diverses raisons il peut être intéressant d'avoir une version embarquée plus compacte :
- pas d'accès disque ;
- que quelques projections nécessaires ;
- réduire le nombre de fichiers jar ;
- pour les applets java et javafx.

Cette fois, il est nécessaire de renseigner toutes les projections dont on aura besoin dans un fichier nommé epsg.properties.
Placez ensuite celui-ci dans votre projet au chemin suivant :
org\geotoolkit\referencing\factory\epsg\epsg.properties.
Chaque projection doit être sur une seule ligne :

4326=GEOGCS["WGS 84", ... AUTHORITY["EPSG","4326"]]
3395=PROJCS["WGS 84 / World Mercator", ... AUTHORITY["EPSG","3395"]]

La majeure partie du temps, la projection sera lue automatiquement avec les différents fichiers SIG. Mais il arrive de vouloir obtenir une projection directement, souvent à partir d'un code du type EPSG:XXXX.
Voici comment procéder :

CoordinateReferenceSystem crs = CRS.decode("EPSG:27582");

Bien souvent les données sont des projections différentes, il faut donc en reprojeter l'une ou l'autre pour faire notre travail sereinement. Voici comment faire pour obtenir cette transformation mathématique et l'utiliser :

CoordinateReferenceSystem inCRS = CRS.decode("EPSG:4326");
CoordinateReferenceSystem outCRS = CRS.decode("EPSG:27582");

MathTransform trs = CRS.findMathTransform(inCRS, outCRS);

DirectPosition d1 = new DirectPosition2D(inCRS, 45, 56);
System.out.println("FROM EPSG:4326 = " + d1);

DirectPosition d2 = trs.transform(d1, null);
System.out.println("TO EPSG:27582 = " + d2);

Cette opération demande à GeotoolKit de disposer d'énormément d'informations ainsi que dans certains cas de grilles de transformations. Deux types d'erreurs surviennent régulièrement.

Type : EPSG:XXX Not found…
Erreur classique quand vous avez mal configuré votre source de données pour les projections.
Type : Missing Bursa Wolf… Les paramètres Bursa-Wolf sont des valeurs permettant de définir la transformation vers un système neutre. Il est possible de forcer à trouver une transformation au prix d'une perte de précision plus ou moins conséquente, pour cela il faut exécuter le code suivant au démarrage de votre application :

Hints.putSystemDefault(Hints.LENIENT_DATUM_SHIFT, Boolean.TRUE);

Contrairement à la librairie Proj4, GeotoolKit va beaucoup plus loin dans la définition et l'exactitude de celle-ci. Mais surtout elle ne se limite pas à deux dimensions.

Vous pouvez construire des projections avec autant d'axes que vous voulez. Voici par exemple une projection à quatre dimensions :

CoordinateReferenceSystem crs2D = CRS.decode("EPSG:27582");
TemporalCRS temporalAxis = DefaultTemporalCRS.JULIAN;
VerticalCRS verticalAxis = DefaultVerticalCRS.GEOIDAL_HEIGHT;
CompoundCRS crs4D = new DefaultCompoundCRS("MyCRS4D", crs2D, temporalAxis, verticalAxis);

Ainsi que quelques méthodes utilitaires :

verticalAxis = CRS.getVerticalCRS(crs4D);
temporalAxis = CRS.getTemporalCRS(crs4D);

Je vous invite à explorer la javadoc pour ceux qui veulent aller plus loin encore et créer de toutes pièces les projections et formules mathématiques.