Le portail francophone de la géomatique

Je crois même qu'il est dans l'annuaire de liens...qui reste une solution d'avenir en complément des moteurs de recherches.

Merci pour ces premières réponses.

J'ai mal formulé mon titre, car je ne cherche pas la liste des projections du monde : je cherche à trouver des projections à partir d'un set de coordonnées dans une projection inconnue. Les deux sites (que je connais bien) sont des aides efficaces, mais aucun n'égale l'efficacité de Dogeo pour savoir rapidement dans quelle projection on se situe.

Le site TWCC dispose d'un nombre plus important de projections mais il faut manuellement faire défiler les projections pour voir leurs coordonnées quand on clique sur la carte, et il n'y a pas de moyen simple d'afficher toutes les possibilités pour un set de données à projection inconnue.

Dans l'idée, si un pays a un nombre réduit de projections, je peux m'en sortir, mais s'il y en a 10, 20, 30, je ne veux pas perdre du temps à toutes les tester. L'outil dogeo pour ça est très efficace pour retrouver ses petits mais j'aimerais avoir le même genre des résultats... avec plus de points sur la carte pour plus de projections possibles.

Salut,

Voici un exemple sommaire avec OGR en Python. C'est pas l'API la plus élégante mais il a de bonnes ressources ici et le code doit pouvoir se transposer facilement dans une autre API (Arcpy, pyQGIS, Fiona ...).

# -*- coding:utf-8 -*-
from osgeo import osr
from osgeo import ogr

#unknow coords (ici Toulouse en Lambert 2 carto 27572)
x = 528000
y = 1845000

#init wgs84 proj
wgs84 = osr.SpatialReference()
wgs84.ImportFromEPSG(4326)

#projection potentielles à tester (ici projections Fr NTF et RGF)
potentialEPSG = [
3942, 3943, 3944, 3945, 3946, 3947, 3948, 3949, 3950, 2154,
27571, 27572, 27573, 27574, 27561, 27562, 27563, 27564
]

#create shapefile in wgs84
driver = ogr.GetDriverByName("ESRI Shapefile")
data_source = driver.CreateDataSource("coordsTest.shp")
layer = data_source.CreateLayer("coordsTest", wgs84, ogr.wkbPoint)
layer.CreateField(ogr.FieldDefn("EPSG", ogr.OFTInteger))


for epsg in potentialEPSG:
    #reproj from potential epsg to wgs84
    prj = osr.SpatialReference()
    prj.ImportFromEPSG(epsg)
    transfo = osr.CoordinateTransformation(prj,wgs84)
    result = transfo.TransformPoint(x, y)

    #Add reproj point and data to wgs84 shp
    feature = ogr.Feature(layer.GetLayerDefn())
    point = ogr.Geometry(ogr.wkbPoint)
    point.AddPoint(result[0], result[1])
    feature.SetGeometry(point)
    feature.SetField("EPSG", epsg)
    layer.CreateFeature(feature)
    feature.Destroy()

data_source.Destroy()

L'approche me parait plus efficace qu’enchaîner des commandes ogr2ogr. Le problème est de savoir quels sont les systèmes de projection que l'on veut tester. Si à chaque fois on connaît le pays dans lequel se trouve le point il suffirait de faire une requête sur la base EPSG (http://www.epsg.org/) pour extraire au préalable toutes les projections concernant ce pays puis les tester.

Retour promis : ca marche du tonnerre, je suis ravi
Merci !!!!

Pour ceux qui souhaitent l'utiliser (exemple) :
1. Ouvrir QGIS
2. Ouvrir la console Python (Menu Extension/Console Python)
3. Cliquer sur le bouton "Editeur" : une feuille de papier avec un crayon
4. Coller le code / enregistrer sous... quelque part
5. Changer x = 528000 et y = 1845000 pour de nouvelles coordonnées si besoin
6. Changer les chiffres dans potentialEPSG = [3942, ...27564] si besoin
7. Changer l'emplacement du fichier créé, par exemple : data_source = driver.CreateDataSource("C:/tmp/coordsTest.shp")
8. Cliquer sur le bouton Exécuter le script (triangle vers la droite). Agrandir la fenêtre vers le haut si vous ne le voyez pas.
9. Admirez la maîtrise de Dominique

Tu veux dire en adaptant le code ? Why not.  A suivre !