আমি ডাব্লুজিএস ৮৮ লেট / দীর্ঘের পয়েন্ট পেয়েছি এবং আমি তাদের মধ্যে "ছোট" (5 কিলোমিটার কম বলার চেয়ে কম) পরিমাপ করতে চাই।
আমি http://www.movable-type.co.uk/scriptts/latlong.html থেকে হ'রসিন সূত্রটি ব্যবহার করতে পারি এবং এটি খুব ভালভাবে কাজ করে।
যদিও আমি পাইথন শেপলি লাইব্রেরিগুলি ব্যবহার করতে চাই, যাতে আমি কেবল দূরত্বের চেয়ে আরও বেশি অপারেশন করতে পারি এবং যে পরিমাণে আমি কাজ করছি তার জন্য সমতল পৃথিবী যথেষ্ট পরিমাণে আনুমানিক। ভৌগলিক সমন্বয়গুলি কার্টেসিয়ান কোর্ডে নির্ভরযোগ্যভাবে প্রজেক্ট করতে, আমি পাইথনের ব্যবহার করছি proj4তবে আমার চেয়ে বড় ত্রুটি পাওয়া যাবে বলে মনে হচ্ছে।
যদি আমি স্থানীয় ইউটিএম জোন ব্যবহার করি তবে আমি কয়েক মিটার হ্যাওয়ারিনের মধ্যে পার্থক্য পেয়েছি, যা ভাল। তবে আমি ইউটিএম জোনটি বের করতে চাই না (পয়েন্টগুলি বিশ্বব্যাপী হতে পারে), তাই আমি "গোলাকৃতির মার্কেটেটার" দিয়ে চেষ্টা করেছি তবে এখন হেরসাইন এবং প্রজেক্টেড দূরত্বের মধ্যে পার্থক্য 100% এরও বেশি ভাল are এটি কি গোলাকৃতির মার্কেটারের পক্ষে সত্যই সঠিক? আমি যা চাই তা হ'ল বিশ্বের যে কোনও স্থানে একে অপরের 5 কিলোমিটারের মধ্যে দুটি পয়েন্টের জন্য একটি কার্যক্ষম কার্টেসিয়ান প্রক্ষেপণ।
from shapely.geometry import Point
from pyproj import Proj
proj = Proj(proj='utm',zone=27,ellps='WGS84')
#proj = Proj(init="epsg:3785") # spherical mercator, should work anywhere...
point1_geo = (-21.9309694, 64.1455718)
point2_geo = (-21.9372481, 64.1478206)
point1 = proj(point1_geo[0], point1_geo[1])
point2 = proj(point2_geo[0], point2_geo[1])
point1_cart = Point(point1)
point2_cart = Point(point2)
print "p1-p2 (haversine)", hdistance(point1_geo, point2_geo)
print "p1-p2 (cartesian)", point1_cart.distance(point2_cart)
এই মুহুর্তে, তাদের মধ্যবর্তী স্থানে হরসিন দূরত্ব 394 মি এবং ইউটিএম জোন 27, 395 মি ব্যবহার করছে। তবে আমি যদি গোলাকার মারকেটর ব্যবহার করি তবে কার্টেসিয়ান দূরত্ব 904 মিটার, যা বন্ধ।