تحويل إحداثيات إسقاط ESRI إلى LAT-LNG

Question

لدي مجموعات بيانات كبيرة من X، Y إحداثيات في "NAD 1983 StatePlane Michigan South Fips 2113 قدم" (ويعرف أيضا باسم ESRI 102690.). أود تحويلها إلى نقاط LNG للغاز الطبيعي المسال.

من الناحية النظرية، هذا شيء بجن تم تصميمه للتعامل معه، لكن الوثائق لم تمنح لي فكرة - يبدو أنها تصف حالات أكثر تعقيدا.

لقد حاولت استخدام واجهة بيثون, ، مثل ذلك:

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

ولكن هذا يعطي إخراج غير صحيح بعنف.

هناك مكتبة جافا سكريبت، ولكن لدي ~ 50،000 نقطة للتعامل معها، لذلك لا يبدو من المناسب.

---

ما كان يعمل بالنسبة لي:

لقد أنشأت ملف يسمى ptest مع كل زوج على خط خاص به، فإن x و y الإحداثيات مفصولة عن طريق الفضاء، مثل ذلك:

13304147.06410000000 288651.94040000000
...

ثم قمت بإنقاذ هذا الملف في الأمر وأنقل النتائج إلى ملف الإخراج:

$>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

إذا كنت بحاجة فقط إلى تحديث ويمكن أن تفعل بعض التعدين البيانات على الملفات النصية التي تستخدم ما تريد واستخدامها http://spatialreference.org/ref/esri/102690/ كمرجع.

على سبيل المثال، استخدم ProJ4 وتخزينه في ملف Shell / CMD واتصل بملف الإدخال الخاص بك باستخدام ProJ4 (إصدار Linux / Windows المتوفر) لا توجد مشكلة في حجم DataSet الخاص بك.

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}

لذلك في حالتك شيء مثل هذا:

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

إذا كان لديك إحداثيات في ملف TXT أو CSV أو XLS، فيمكنك القيام به CTRL + C وإدراجها http://cs2cs.mygeodata.eu. حيث يمكنك ضبط نظام الإدخال المناسب ونظام تنسيق الإخراج المطلوب. من الممكن إدراج الآلاف من الإحداثيات في تنسيقات مختلفة ...