计算球面上两点间最短距离多语言版本

最新推荐文章于 2024-05-11 17:17:43 发布

Robotec

最新推荐文章于 2024-05-11 17:17:43 发布

阅读量2.4k

点赞数

分类专栏： GIS 文章标签： GIS 球面距离

本文链接：https://blog.csdn.net/ifreely0204/article/details/53300602

版权

本文介绍了使用Haversine球面半正矢公式来计算球面上两点间最短距离的方法，并提供了多种编程语言的实现示例。

摘要由CSDN通过智能技术生成

Haversine球面半正矢公式，用来计算球面上连点之间的最短距离，
这里有几种语言版本的实现方式，就摘了几种，还有其他例如Smalltalk, Tcl, NCL的就不贴了。

//C#
using System;
namespace HaversineFormula
{
    /// <summary>
    /// 距离类型，英里/公里
    /// </summary>
    public enum DistanceType { Miles, Kilometers };
    /// <summary>
    /// Specifies a Latitude / Longitude point.
    /// </summary>
    public struct Position
    {
        public double Latitude;
        public double Longitude;
    }
    class Haversine
    {
        /// <summary>
        /// 返回英里或者公里的距离
        /// latitude / longitude points.
        /// </summary>
        /// <param name=”pos1″></param>
        /// <param name=”pos2″></param>
        /// <param name=”type”></param>
        /// <returns></returns>
        public double Distance(Position pos1, Position pos2, DistanceType type)
        {
            double R = (type == DistanceType.Miles)&nbsp;? 3960&nbsp;: 6371;
            double dLat = this.toRadian(pos2.Latitude - pos1.Latitude);
            double dLon = this.toRadian(pos2.Longitude - pos1.Longitude);
            double a = Math.Sin(dLat / 2) * Math.Sin(dLat / 2) +
                Math.Cos(this.toRadian(pos1.Latitude)) * Math.Cos(this.toRadian(pos2.Latitude)) *
                Math.Sin(dLon / 2) * Math.Sin(dLon / 2);
            double c = 2 * Math.Asin(Math.Min(1, Mat