请帮我写一段代码,要求是在一定经纬度范围内随机生成5个点,设点的移动速度为70m每分钟,让这五个点在该经纬度范围内随机运动1h,其运动范围不能超过先前设置的经纬度范围,求每个点每分钟的经纬度坐标,将结果输入excel表,其中一列为节点编号,一列为时间,一列为纬度,一列为经度。
时间: 2024-05-08 15:18:44 浏览: 116
好的,我可以帮您写这段代码。以下是代码实现:
```python
import random
import math
import pandas as pd
# 设置经纬度范围
lon_range = [114.0, 114.5] # 经度范围
lat_range = [22.5, 23.0] # 纬度范围
# 获取两个经纬度之间的距离(单位:米)
def get_distance(lat1, lon1, lat2, lon2):
rad_lat1 = math.radians(lat1)
rad_lat2 = math.radians(lat2)
a = rad_lat1 - rad_lat2
b = math.radians(lon1) - math.radians(lon2)
s = 2 * math.asin(math.sqrt(math.pow(math.sin(a / 2), 2) + math.cos(rad_lat1) * math.cos(rad_lat2) * math.pow(math.sin(b / 2), 2)))
s *= 6371.393 * 1000
return s
# 随机生成5个点
points = []
for i in range(5):
point = {}
point['id'] = i + 1
point['lat'] = random.uniform(lat_range[0], lat_range[1])
point['lon'] = random.uniform(lon_range[0], lon_range[1])
points.append(point)
# 计算每个点每分钟的经纬度坐标
results = []
for i in range(60):
for point in points:
# 随机生成角度和速度
angle = random.uniform(0, 360)
speed = 70
# 计算移动距离
distance = (speed / 60) * 1000
# 计算纬度和经度的偏移量
lat_offset = distance * math.cos(math.radians(angle)) / 111319.9
lon_offset = distance * math.sin(math.radians(angle)) / (111319.9 * math.cos(math.radians(point['lat'])))
# 计算新的经纬度坐标
new_lat = point['lat'] + lat_offset
new_lon = point['lon'] + lon_offset
# 判断新坐标是否在规定范围内
if lat_range[0] <= new_lat <= lat_range[1] and lon_range[0] <= new_lon <= lon_range[1]:
# 保存结果
result = {}
result['id'] = point['id']
result['time'] = i + 1
result['lat'] = new_lat
result['lon'] = new_lon
results.append(result)
# 将结果输出到excel表
results_df = pd.DataFrame(results, columns=['id', 'time', 'lat', 'lon'])
results_df.to_excel('result.xlsx', index=False)
```
该代码会随机生成5个点,然后计算每个点每分钟的经纬度坐标,并将结果保存到excel表中。该代码中使用了三个函数:
1. `get_distance()`:计算两个经纬度坐标之间的距离。
2. `random.uniform()`:生成一个指定范围内的随机数。
3. `math.radians()` 和 `math.cos()`:将角度转化为弧度,并计算余弦值。
请注意,该代码没有解决随机生成的5个点之间会不会重合的问题。如果您需要保证这5个点不重合,请自行在生成点的循环中添加相应的判断逻辑。
请问还有什么我可以为您效劳的吗?
阅读全文