এক্স, ওয়াই কে পিওপ্রজ এবং প্রজ 4 ব্যবহার করে ল্যাট / লম্বায় স্থানাঙ্কিত করে ভুল স্থানাঙ্ক ফেরত দেয়

Question

আমি একটি পাইথন স্ক্রিপ্ট লিখছি যা এক্স এবং ওয়াই সমন্বয়যুক্ত একাধিক এক্সএমএল ফাইল পড়ে এবং সেগুলি একক সিএসভি ফাইলের সাথে সংযুক্ত করে। অক্ষাংশ এবং দ্রাঘিমাংশ সিএসভিতে প্রয়োজনীয় ক্ষেত্র, তবে ওহিও উত্তর স্টেট প্লেন ইউএসএফটি-কে ডাব্লুজিএস 84 এ রূপান্তর করতে আমার এক্স, ওয়াই স্থানাঙ্কে অসুবিধা হচ্ছে।

>>> p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters
>>> p(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)
>>> p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code
>>> p1(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)

উপরের দুটি পদ্ধতিই একই ফলস্বরূপ ফিরে আসে, তবে এই দীর্ঘতম হডসন উপসাগরের কোথাও। আমি যখন আর্কম্যাপে স্থানাঙ্কগুলি প্লট করি তখন সঠিক লম্বা লম্বা হয়: -81.142311,41.688205।

* লক্ষ্য করুন সমস্ত ল্যাট লম্বা দীর্ঘ প্রদান করা হয়, ল্যাট যেহেতু প্রজ এটি ব্যবহার করে

কেউ কি জানেন যে আমি কেন প্রজ 4 এবং পাইপ্রোজ থেকে ভুল স্থানাঙ্ক পাব?

Answer 1

আমি যখন চালাচ্ছি তখন gdaltransform@geographichika এবং প্রোজ ৪.৪ সরঞ্জামের মতো একই ফলাফল পেয়েছি cs2cs:

$ gdaltransform -s_srs EPSG:3734 -t_srs EPSG:4326
739400.9 2339327.3             
-87.3195485720169 45.9860670658218 0

cs2cs +proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs +to +proj=lonlat +datum=WGS84
739400.9 2339327.3
87d19'10.375"W 45d59'9.841"N 0.000

আপনার পয়েন্টের x এবং y স্থানাঙ্কগুলির বিপরীত করা ফলস্বরূপ দেয় যা আপনি আর্কম্যাপে দেখছেন:

gdaltransform -s_srs EPSG:3734 -t_srs EPSG:4326
2339327.3 739400.9
-81.1423086719059 41.6882035384526 0

সুতরাং আপনার এক্স এবং y সঠিক রাউন্ডে স্থানাঙ্ক রয়েছে তা নিশ্চিত করার জন্য আপনাকে একটি ভিজ্যুয়াল চেক করতে হবে। এটি এমন একটি সমস্যা যা আমি অতীতে ছিলাম যখন আপনি দুটি ফলাফল পান যা যথেষ্ট পরিমাণে সমান হয় আপনি এটিকে গোলটি ত্রুটি বা কোনও কিছুতে রেখে দেন।

Answer 2

পাইপ্রোজ অনুমান করে যে আপনার স্থানাঙ্কগুলি মিটারে রয়েছে। আমি অনুমান করতে পারি যে ফুট / মিটার সম্পর্কিত কোনও বিষয়ই সমস্যার কারণ।

দীর্ঘস্থায়ী আর্গুমেন্ট সহ প্রজ শ্রেণীর উদাহরণকে কল করা ল্যাট ল্যান্ট / ল্যাট (ডিগ্রিতে) এক্স / ওয়াই নেটিভ ম্যাপ প্রজেকশন স্থানাঙ্কগুলিতে (মিটারে) রূপান্তর করবে

Preচ্ছিক কীওয়ার্ড 'preserv_units' সত্য হলে মানচিত্রের অভিক্ষেপ স্থানাঙ্কগুলির ইউনিটগুলি মিটার হতে বাধ্য হয় না।

http://pyproj.googlecode.com/svn/trunk/docs/pyproj.Proj-class.html

আপনার প্রাথমিক স্থানাঙ্কগুলি কি পায়ে রয়েছে? আপনি যখন আর্কম্যাপে ডেটা লোড করবেন তখন মানচিত্রটি কোন ইউনিট ব্যবহার করে?

এটি স্থানাঙ্কগুলি আরও কাছাকাছি পায়:

p1 = Proj(init="epsg:3734")
#1 foot = 0.3048 meters
conv = 0.3048
print p1(739400.91 * conv,2339327.3 * conv,inverse=True)
(-87.3195533069909, 45.98605408134072)

একটি অনুরূপ ইস্যু এখানে পাওয়া যাবে ।

Comments

আপনাকে অনেক ধন্যবাদ. সংরক্ষণ করুন_আনাইটস যুক্তিটি নিশ্চয়ই কৌতুকটি করেছে তবে স্থানাঙ্কগুলি এখনও ভুল। @ মের্সিভিকিং এই উত্তরটি আমাকে সঠিক স্থানাঙ্ক দিয়েছে। আমি উভয় উত্তরকে উত্তর হিসাবে চিহ্নিত করতে পারি কারণ তারা উভয়েই সহায়তা করেছিল।

— ব্রায়ান

ঠিক আছে যদি লোকেরা আমার চেয়ে জিওগ্রাফিকের উত্তরটিকে বেশি উত্সাহিত করে, তবে তা সমস্তই বেরিয়ে আসবে :) খুশী সব কিছু এইভাবে কাজ করেছে।

— MerseyViking 15

যেহেতু লিঙ্কটি অসম্পূর্ণ, এটা যে আপনি লিখতে পারেন দেখানোর জন্য সহায়ক হতে পারে:p1 = Proj( init="epsg:3734", preserve_units=True )

— BenjaminGolder

Answer 3

আমি আসলে ওএইচ দক্ষিণ রাজ্যের বিমানের গ্রিড বাদে একই জিনিসটি করার চেষ্টা করছিলাম এবং আমি আপনার প্রশ্ন জুড়ে এসেছি। আমি 3735 এর সাথে ভুল ফলাফল পাচ্ছিলাম, এখন আমি 3729 দিয়ে সঠিক ফলাফল পেয়েছি I আমি আশা করি আপনি যদি 3734 থেকে 3728 এ পরিবর্তন করেন তবে আপনি সঠিক ফলাফল পাবেন।

EPSG: 3728: NAD83 (NSRS2007) / ওহিও উত্তর (ftUS) EPSG: 3729: NAD83 (NSRS2007) / ওহিও দক্ষিণ (ftUS) EPSG: 3734: NAD83 / ওহিও উত্তর (ftUS) EPSG: 3735: NAD83 / ওহিও দক্ষিণ (ftUS)

আমি আপনার সরবরাহিত ল্যাট, লম্বা এবং এক ফুটেরও কম দুরে আছি।

পি 2 = পাইপ্রোজ.প্রজ (আরম্ভ = "ইপসিজি: 3728", সংরক্ষণ_ইনিতস = সত্য)

P2 (-81.142311,41.688205)

(2339326.6558868014, 739401.4226131936)

বনাম 2339327.3, 739400.91