2021-11-16-01

JLChen

2021-11-16 b8d7b4fa382f7c779fb9e81a9f049d94a9083fa3

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
using System; 
using System.Collections.Generic; 
using System.Linq; 
using System.Text; 
using System.Threading.Tasks; 
 
namespace Infrastructure.Service.Helper 
{ 
    public class CalculatedDistance 
    { 
        private static double HaverSin(double theta) 
        { 
            var v = Math.Sin(theta / 2); 
            return v * v; 
        } 
 
        static double EARTH_RADIUS = 6371.0;//km 地球半径 平均值，千米 
 
        /// <summary> 
        /// 给定的经度1，纬度1；经度2，纬度2. 计算2个经纬度之间的距离。 
        /// </summary> 
        /// <param name="lat1">经度1</param> 
        /// <param name="lon1">纬度1</param> 
        /// <param name="lat2">经度2</param> 
        /// <param name="lon2">纬度2</param> 
        /// <returns>距离（公里、千米）</returns> 
        public static int Distance(double lat1, double lon1, double lat2, double lon2) 
        { 
            //用haversine公式计算球面两点间的距离。 
            //经纬度转换成弧度 
            lat1 = ConvertDegreesToRadians(lat1); 
            lon1 = ConvertDegreesToRadians(lon1); 
            lat2 = ConvertDegreesToRadians(lat2); 
            lon2 = ConvertDegreesToRadians(lon2); 
 
            //差值 
            var vLon = Math.Abs(lon1 - lon2); 
            var vLat = Math.Abs(lat1 - lat2); 
 
            //h is the great circle distance in radians, great circle就是一个球体上的切面，它的圆心即是球心的一个周长最大的圆。 
            var h = HaverSin(vLat) + Math.Cos(lat1) * Math.Cos(lat2) * HaverSin(vLon); 
 
            var distance = (int)((2 * EARTH_RADIUS * Math.Asin(Math.Sqrt(h))) * 1000); 
 
            return distance; 
        } 
 
        /// <summary> 
        /// 将角度换算为弧度。 
        /// </summary> 
        /// <param name="degrees">角度</param> 
        /// <returns>弧度</returns> 
        private static double ConvertDegreesToRadians(double degrees) 
        { 
            return degrees * Math.PI / 180; 
        } 
 
        private static double ConvertRadiansToDegrees(double radian) 
        { 
            return radian * 180.0 / Math.PI; 
        } 
    } 
}