# -*- coding: utf-8 -*-
"""
Created on Thu Jun 14 18:44:06 2018
@author: dell
"""
import re
import pandas as pd
import numpy as np
from math import radians, cos, sin, asin, sqrt
#经纬度计算
def haversine(lon1, lat1, lon2, lat2): # 经度1,纬度1,经度2,纬度2 (十进制度数)
"""
Calculate the great circle distance between two points
on the earth (specified in decimal degrees)
"""
# 将十进制度数转化为弧度
lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])
# haversine公式
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
c = 2 * asin(sqrt(a))
r = 6371 # 地球平均半径,单位为公里
return c * r
#速度计算函数
def speedcal(lng1,lat1,lng2,lat2,time3):
distance=haversine(lng1,lat1,lng2,lat2)
speed=distance/time3*3600
return speed
计算道路等级
最新推荐文章于 2023-06-20 10:09:03 发布
该博客介绍了一个计算道路等级的程序,通过经纬度计算距离,速度,然后根据不同的速度范围给道路评级(A, B, C, D)。程序读取数据,提取小时信息,按路段和小时进行数据切片,计算每条道路在不同时间段的服务水平,并保存结果到CSV文件。"
84839193,1407364,React热更新实践与状态保存,"['前端开发', 'React', 'Webpack', '开发工具']
摘要由CSDN通过智能技术生成