Autre projection ESRI coordonnées de latitude-lng

Question

J'ai un grand ensemble de données de coordonnées x, y dans "NAD 1983 StatePlane Michigan FIPS Sud 2113 Feet" (alias 102690 Esri). Je voudrais les convertir en des points lat-lng.

En théorie, cela est quelque chose proj est conçu pour gérer, mais la documentation n'a pas donné moi un indice -. il semble décrire des cas beaucoup plus complexes

Je l'ai essayé d'utiliser un python Interface , comme suit:

from pyproj import Proj
p = Proj(init='esri:102690')
sx = 13304147.06410000000 #sample points
sy = 288651.94040000000
x2, y2 = p(sx, sy, inverse=True)

Mais cela donne une sortie incorrecte d'une manière extravagante.

Il y a une bibliothèque Javascript, mais je ~ 50.000 points à traiter, ce qui ne semble pas approprié.

---

Ce qui a fonctionné pour moi:

I créé un fichier appelé ptest avec chaque paire sur une ligne, coordonnées x et y séparées par un espace, comme suit:

13304147.06410000000 288651.94040000000
...

Ensuite, je nourris ce fichier dans la commande et les résultats passepoilées à un fichier de sortie:

$>cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 
+lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 
+datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to 
+proj=latlon ptest > out.txt

Answer 1

Si vous ne devez reprojeter et peut faire des données sur l'exploitation de vos fichiers texte utiliser tout ce que vous aimez et utilisez http://spatialreference.org/ref/esri/102690/ comme référence.

Par exemple, utilisez le Proj4 et la stocker dans un fichier shell / cmd et appelez votre fichier d'entrée avec proj4 (version linux / windows disponible) aucun problème avec la taille de votre ensemble de données.

cs2cs +proj=latlong +datum=NAD83 +to +proj=utm +zone=10  +datum=NAD27 -r <<EOF
cs2cs -f %.16f +proj=utm +zone=20N +to +proj=latlon - | awk '{print $1 " " $2}

dans votre cas quelque chose comme ceci:

cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 +lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to +proj=latlon

http://trac.osgeo.org/proj/wiki/man_cs2cs

http://trac.osgeo.org/proj/

Answer 2

Si vous avez coordonnées dans TXT, CSV ou un fichier XLS, vous pouvez faire CTRL + C et les insérer http: // cs2cs.mygeodata.eu où vous pouvez définir l'entrée et la sortie désirée appropriée du système de coordonnées. Il est possible d'insérer des milliers de coordonnées dans différents formats ...