帮我写一个python函数,函数中要用到地球赤道半径,能够返回两组经纬坐标的地球球面距离
时间: 2023-09-04 18:04:28 浏览: 143
### 回答1:
好的,以下是 Python 函数,可以计算两组经纬坐标的地球球面距离:
```python
import math
def calc_distance(lat1, lon1, lat2, lon2):
# 地球赤道半径,单位为米
earth_radius = 6378137.
# 将经纬度转换为弧度
lat1_rad = math.radians(lat1)
lon1_rad = math.radians(lon1)
lat2_rad = math.radians(lat2)
lon2_rad = math.radians(lon2)
# 计算经纬度差值
delta_lon = lon2_rad - lon1_rad
delta_lat = lat2_rad - lat1_rad
# 使用 Haversine 公式计算球面距离
a = math.sin(delta_lat / 2) ** 2 + math.cos(lat1_rad) * math.cos(lat2_rad) * math.sin(delta_lon / 2) ** 2
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))
distance = earth_radius * c
return distance
```
其中,`lat1`、`lon1` 表示第一组经纬度,`lat2`、`lon2` 表示第二组经纬度,函数返回的是两组经纬度之间的球面距离,单位为米。
### 回答2:
可以使用haversine公式来计算两组经纬度坐标在地球球面上的距离。该公式基于地球是一个近似于规则球体的假设。下面是一个使用Python编写的函数来计算地球球面距离的例子:
```python
import math
def haversine(lat1, lon1, lat2, lon2):
# 地球赤道半径(单位:千米)
R = 6371.0
# 将角度转换为弧度
lat1 = math.radians(lat1)
lon1 = math.radians(lon1)
lat2 = math.radians(lat2)
lon2 = math.radians(lon2)
# 计算经纬度的差值
dlon = lon2 - lon1
dlat = lat2 - lat1
# 使用haversine公式计算地球球面上的距离
a = math.sin(dlat/2)**2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon/2)**2
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
distance = R * c
return distance
# 示例调用
lat1 = 23.4567
lon1 = 45.6789
lat2 = 12.3456
lon2 = 34.5678
result = haversine(lat1, lon1, lat2, lon2)
print("地球球面距离:", result, "千米")
```
在这个示例中,`haversine()`函数接受四个参数,分别是第一组经纬度坐标的纬度(`lat1`)、经度(`lon1`),以及第二组经纬度坐标的纬度(`lat2`)、经度(`lon2`)。函数返回这两组坐标之间的地球球面距离(单位:千米)。在示例调用中,我们假设第一组经纬度坐标是(23.4567, 45.6789),第二组经纬度坐标是(12.3456, 34.5678),我们将结果打印出来。
### 回答3:
下面是一个Python函数,用于计算地球球面上两组经纬坐标的距离。
```python
import math
def calculate_distance(lat1, lon1, lat2, lon2):
# 地球赤道半径(单位:米)
R = 6371000
# 将经纬度转换为弧度
lat1_rad = math.radians(lat1)
lon1_rad = math.radians(lon1)
lat2_rad = math.radians(lat2)
lon2_rad = math.radians(lon2)
# 计算经纬度差值
dlat = lat2_rad - lat1_rad
dlon = lon2_rad - lon1_rad
# 使用Haversine公式计算球面距离
a = math.sin(dlat/2) * math.sin(dlat/2) + math.cos(lat1_rad) * math.cos(lat2_rad) * math.sin(dlon/2) * math.sin(dlon/2)
c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
distance = R * c
return distance
```
使用示例:
```python
# 经度和纬度(示例数据)
lat1 = 39.9042
lon1 = 116.4074
lat2 = 31.2304
lon2 = 121.4737
# 计算距离
distance = calculate_distance(lat1, lon1, lat2, lon2)
print("两组坐标之间的地球球面距离为:", distance, "米")
```
注意:在使用该函数时,经度应为正数表示东经,为负数表示西经;纬度应为正数表示北纬,为负数表示南纬。
阅读全文