已知2点的经纬度如何计算距离
时间: 2024-05-06 13:15:27 浏览: 181
可以使用Haversine公式来计算2点之间的距离,公式如下:
a = sin²(Δlat/2) + cos(lat1) * cos(lat2) * sin²(Δlong/2)
c = 2 * atan2( √a, √(1−a) )
d = R * c
其中,Δlat和Δlong分别为两点的纬度差和经度差,lat1和lat2分别为两点的纬度,R为地球半径(一般取平均值6371km),d为两点间距离。
下面是一个Python实现示例:
```python
from math import radians, sin, cos, sqrt, atan2
def distance(lat1, long1, lat2, long2):
R = 6371 # 地球半径,单位为公里
lat1, long1, lat2, long2 = map(radians, [lat1, long1, lat2, long2])
dlat = lat2 - lat1
dlong = long2 - long1
a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlong/2)**2
c = 2 * atan2(sqrt(a), sqrt(1-a))
distance = R * c
return distance
```
其中,lat1和long1为第一个点的纬度和经度,lat2和long2为第二个点的纬度和经度,函数返回的是两点之间的距离,单位为公里。
相关问题
已知2点的经纬度计算距离java实现
可以使用 Haversine formula 来计算两点之间的距离,以下是一个 Java 实现的示例代码:
```java
public class DistanceCalculator {
private static final double EARTH_RADIUS = 6371; // 地球半径,单位为千米
public static double calculateDistance(double lat1, double lon1, double lat2, double lon2) {
double radLat1 = Math.toRadians(lat1);
double radLat2 = Math.toRadians(lat2);
double radLon1 = Math.toRadians(lon1);
double radLon2 = Math.toRadians(lon2);
double a = Math.sin((radLat1 - radLat2) / 2) * Math.sin((radLat1 - radLat2) / 2)
+ Math.cos(radLat1) * Math.cos(radLat2) * Math.sin((radLon1 - radLon2) / 2) * Math.sin((radLon1 - radLon2) / 2);
double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
return EARTH_RADIUS * c;
}
}
```
其中,`lat1` 和 `lon1` 表示第一个点的纬度和经度,`lat2` 和 `lon2` 表示第二个点的纬度和经度。函数返回值为两点之间的距离,单位为千米。
matlab已知两点经纬度求距离
可以使用 MATLAB 自带的 geodetic2ecef 和 distance 函数来计算两点之间的距离。
假设有两个点的经纬度分别为 (lat1, lon1) 和 (lat2, lon2),则可以按照以下步骤计算它们之间的距离:
1. 将经纬度转换为地心坐标系下的坐标。可以使用 geodetic2ecef 函数来实现:
```
[x1, y1, z1] = geodetic2ecef(lat1, lon1, 0);
[x2, y2, z2] = geodetic2ecef(lat2, lon2, 0);
```
2. 计算两点之间的距离。可以使用 distance 函数来计算:
```
dist = distance('gc', lat1, lon1, lat2, lon2);
```
其中,'gc' 表示计算大圆航线距离。注意,distance 函数的输入参数是经纬度,而不是地心坐标系下的坐标。
完整的 MATLAB 代码如下:
```
lat1 = 39.9; % 第一个点的纬度
lon1 = 116.4; % 第一个点的经度
lat2 = 31.2; % 第二个点的纬度
lon2 = 121.5; % 第二个点的经度
[x1, y1, z1] = geodetic2ecef(lat1, lon1, 0);
[x2, y2, z2] = geodetic2ecef(lat2, lon2, 0);
dist = distance('gc', lat1, lon1, lat2, lon2);
fprintf('两点之间的距离为 %.2f 千米\n', dist);
```
注意,这只是一个简单的示例,实际应用中需要考虑更多的因素,例如地球的椭球形状等。
阅读全文