কিছুক্ষণ আগে আমি হ'রসিন সূত্রের একটি উদাহরণ লিখেছিলাম এবং এটি আমার ওয়েবসাইটে প্রকাশ করেছি:
/**
* Calculates the great-circle distance between two points, with
* the Haversine formula.
* @param float $latitudeFrom Latitude of start point in [deg decimal]
* @param float $longitudeFrom Longitude of start point in [deg decimal]
* @param float $latitudeTo Latitude of target point in [deg decimal]
* @param float $longitudeTo Longitude of target point in [deg decimal]
* @param float $earthRadius Mean earth radius in [m]
* @return float Distance between points in [m] (same as earthRadius)
*/
function haversineGreatCircleDistance(
$latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo, $earthRadius = 6371000)
{
// convert from degrees to radians
$latFrom = deg2rad($latitudeFrom);
$lonFrom = deg2rad($longitudeFrom);
$latTo = deg2rad($latitudeTo);
$lonTo = deg2rad($longitudeTo);
$latDelta = $latTo - $latFrom;
$lonDelta = $lonTo - $lonFrom;
$angle = 2 * asin(sqrt(pow(sin($latDelta / 2), 2) +
cos($latFrom) * cos($latTo) * pow(sin($lonDelta / 2), 2)));
return $angle * $earthRadius;
}
➽ দ্রষ্টব্য যে আপনি প্যারামিটারের সাথে পাস করার সাথে একই ইউনিটে দূরত্ব ফিরে পেয়েছেন $earthRadius
। ডিফল্ট মান 71৩71১০০০০ মিটার তাই ফলটি [এম] এও হবে। মাইল থেকে ফলাফল পেতে, আপনি যেমন 3959 মাইল উত্তীর্ণ করতে পারেন $earthRadius
এবং ফলাফলটি হবে [মাইল]। আমার মতে এসআই ইউনিটগুলির সাথে লেগে থাকা ভাল অভ্যাস, অন্যথায় করার কোনও বিশেষ কারণ না থাকলে।
সম্পাদনা:
হিসাবে TreyA সঠিকভাবে নির্দিষ্ট, Haversine সূত্র সঙ্গে দুর্বলতা রয়েছে বিপরীত পৃথিবীপৃষ্ঠ-সংক্রান্ত পয়েন্ট কারণ ত্রুটি rounding (যদিও এটি হয় ছোট দূরত্বের জন্য স্থিতিশীল)। তাদের কাছাকাছি যেতে, আপনি পরিবর্তে ভিনসেন্টি সূত্রটি ব্যবহার করতে পারেন ।
/**
* Calculates the great-circle distance between two points, with
* the Vincenty formula.
* @param float $latitudeFrom Latitude of start point in [deg decimal]
* @param float $longitudeFrom Longitude of start point in [deg decimal]
* @param float $latitudeTo Latitude of target point in [deg decimal]
* @param float $longitudeTo Longitude of target point in [deg decimal]
* @param float $earthRadius Mean earth radius in [m]
* @return float Distance between points in [m] (same as earthRadius)
*/
public static function vincentyGreatCircleDistance(
$latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo, $earthRadius = 6371000)
{
// convert from degrees to radians
$latFrom = deg2rad($latitudeFrom);
$lonFrom = deg2rad($longitudeFrom);
$latTo = deg2rad($latitudeTo);
$lonTo = deg2rad($longitudeTo);
$lonDelta = $lonTo - $lonFrom;
$a = pow(cos($latTo) * sin($lonDelta), 2) +
pow(cos($latFrom) * sin($latTo) - sin($latFrom) * cos($latTo) * cos($lonDelta), 2);
$b = sin($latFrom) * sin($latTo) + cos($latFrom) * cos($latTo) * cos($lonDelta);
$angle = atan2(sqrt($a), $b);
return $angle * $earthRadius;
}