• Gps坐标距离计算C#实现


    园子里找到两钟实现方式,做一记录,回头再认真学习,先拿来就用吧:

    1、@旋风描述的算法:

    场景:已知两个GPS点的经纬度坐标信息。计算两点的距离。

    1. 距离/纬度关系

      GPS: 22.514519,113.380301

      GPS: 22.511962,113.380301

    距离: 284.6439379583341  

    jl_wd=284.6439379583341/(22.51451-22.511962)

    = 111712.69150641055729984301412873 (米/度)

    2. 距离/经度关系

    GPS: 22.514866, 113.388444

    GPS: 22.514866, 113.379378

    距离: 932.2997762326453

    jl_jd =932.2997762326453/(113.388444-113.379378)

    =102834.74258026089786013677476285(米/度)

    3.勾股定理 

    勾股定理:a2+b=c2

     

    看上图中已知道A,B的坐标后,就可以求得A点与B点的纬度差值的绝对值

    wd_c=|(N1,-N2)|

    A点与B点的经度差值的绝对值

    jd_c=|(E2-E1)|

    知道了经度、纬度和距离的关系,又知道了A点与B点的经度差和纬度差,那么就可以求出b的长度和a的长度,a和b求出来后就可以用直三角形的勾股定理求出c边的长度也就是A点也B点的距离。

    b=wd_c*jl_wd

    a=jd_c*jl_jd

    c=√(a+b2)

    计算两点GPS坐标距离 
    
    Code highlighting produced by Actipro CodeHighlighter (freeware)http://www.CodeHighlighter.com/-->   /// <summary>
         ///计算两点GPS坐标的距离
         /// </summary>
         /// <param name="n1">第一点的纬度坐标</param>
         /// <param name="e1">第一点的经度坐标</param>
         /// <param name="n2">第二点的纬度坐标</param>
         /// <param name="e2">第二点的经度坐标</param>
         /// <returns></returns>
            public static double Distance(double n1, double e1, double n2, double e2)
            {
                double jl_jd = 102834.74258026089786013677476285;
                double jl_wd = 111712.69150641055729984301412873;
                double b = Math.Abs((e1 - e2) * jl_jd);
                double a = Math.Abs((n1 - n2) * jl_wd);
                return Math.Sqrt((a * a + b * b));
    
            }
    View Code

    2、据@国之将亡五毛成灾 方案,要考虑地球是球体的,考虑到地球的弧度,经度每个地方长度并不是平均的,经度之间的距离随纬度影响,参考算法

    private const double EARTH_RADIUS = 6378.137;
    private static double rad(double d)
    {
    return d * Math.PI / 180.0;
    }
    
    public static double GetDistance(double lat1, double lng1, double lat2, double lng2)
    {
    double radLat1 = rad(lat1);
    double radLat2 = rad(lat2);
    double a = radLat1 - radLat2;
    double b = rad(lng1) - rad(lng2);
    double s = 2 * Math.Asin(Math.Sqrt(Math.Pow(Math.Sin(a/2),2) +
    Math.Cos(radLat1)*Math.Cos(radLat2)*Math.Pow(Math.Sin(b/2),2)));
    s = s * EARTH_RADIUS;
    s = Math.Round(s * 10000) / 10000;
    return s;
    }
  • 相关阅读:
    浅谈折半搜索
    NOIP PJ/CSP-J 题目选做
    SDN第五次上机作业
    SDN第四次上机作业
    SDN第三次上机作业
    SDN第二次上机作业
    必看
    关于tensor
    permute与transpose
    1823. 找出游戏的获胜者
  • 原文地址:https://www.cnblogs.com/kelite/p/3527731.html
Copyright © 2020-2023  润新知