Suite

Passer de la zone d'étendue ArcGIS à la paire lat long


J'ai une copie de certains fichiers journaux ArcGIS avec des informations sur l'étendue de la zone de délimitation. Cet exemple de point doit correspondre au Montana.

Étendue :-12711495.260857,6062474.786336,-12711267.136963,66062864.358207

Je sais que le format est xMin, yMin, xMax, yMax pour les points. J'ai pris les points pour faire un point central.

xCentre = (xMax + xMin)/2,0 yCentre = (yMax + yMin)/2,0

Il semble que ces points soient en mètres et j'ai du mal à trouver une formule pour les traduire en paires de points lat long. Existe-t-il une formule directe que je peux utiliser dans mon script pour le traitement ? Je n'ai pas de copie d'arcpy à utiliser, ce qui semble être le moyen le plus suggéré pour résoudre le problème.


J'ai trouvé de l'aide dans ce post gis : Comment convertir les coordonnées projetées en latitude/longitude à l'aide de Python ?

et fait une fonction python pour résoudre mon problème

def convertEPSGtoWGS(pointX, pointY): # Système de référence spatiale inputEPSG = 3857 outputEPSG = 4326 # crée une géométrie à partir des coordonnées point = ogr.Geometry(ogr.wkbPoint) point.AddPoint(pointX, pointY) # crée une transformation de coordonnées inSpatialRef = osr. SpatialReference() inSpatialRef.ImportFromEPSG(inputEPSG) outSpatialRef = osr.SpatialReference() outSpatialRef.ImportFromEPSG(outputEPSG) coordTransform = osr.CoordinateTransformation(inSpatialRef, outSpatialRef) #26 transform print point.Transform #(coord print point inSpatialRef, outSpatialRef) #26 transform print point.Transform #(coord print point in GetX(), point.GetY()