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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
package cn.wisenergy.common.utils;
/**
* @Author zyj
* @Date 2020/5/6 18:18
* @Description 计算一个经纬度在另一个经纬度什么方向
*/
public class LngLatDirection {
public static void main(String[] args) {
LngLatDirection lngLatDirection = new LngLatDirection();
String str = getDirection(39.915322, 116.404341, 39.954713, 116.45982);
System.out.println(str);
}
/**
* @param lat1 纬度1
* @param lng1 经度1
* @param lat2 纬度2
* @param lng2 经度2
* @return 方向
*/
public static String getDirection(double lat1, double lng1, double lat2, double lng2) {
double jiaodu = getAngle(lat1, lng1, lat2, lng2);
if ((jiaodu <= 10) || (jiaodu > 350))
return "东";
if ((jiaodu > 10) && (jiaodu <= 80))
return "东北";
if ((jiaodu > 80) && (jiaodu <= 100))
return "北";
if ((jiaodu > 100) && (jiaodu <= 170))
return "西北";
if ((jiaodu > 170) && (jiaodu <= 190))
return "西";
if ((jiaodu > 190) && (jiaodu <= 260))
return "西南";
if ((jiaodu > 260) && (jiaodu <= 280))
return "南";
if ((jiaodu > 280) && (jiaodu <= 350))
return "东南";
return "";
}
private static double getAngle(double lat1, double lng1, double lat2, double lng2) {
double x1 = lng1;
double y1 = lat1;
double x2 = lng2;
double y2 = lat2;
double pi = Math.PI;
double w1 = y1 / 180 * pi;
double j1 = x1 / 180 * pi;
double w2 = y2 / 180 * pi;
double j2 = x2 / 180 * pi;
double ret;
if (j1 == j2) {
if (w1 > w2)
return 270; // 北半球的情况，南半球忽略
else if (w1 < w2)
return 90;
else
return -1;// 位置完全相同
}
ret = 4 * Math.pow(Math.sin((w1 - w2) / 2), 2) - Math.pow(
Math.sin((j1 - j2) / 2) * (Math.cos(w1) - Math.cos(w2)), 2);
ret = Math.sqrt(ret);
double temp = (Math.sin(Math.abs(j1 - j2) / 2) * (Math.cos(w1) + Math
.cos(w2)));
ret = ret / temp;
ret = Math.atan(ret) / pi * 180;
if (j1 > j2) { // 1为参考点坐标
if (w1 > w2)
ret += 180;
else
ret = 180 - ret;
} else if (w1 > w2)
ret = 360 - ret;
return ret;
}
}