python pandas用函数算距离时报错

在网上copy来的一个算经纬度的方法,但在实际操作运行到 lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2]) 的时候到报错:”TypeError: cannot convert the series to <class 'float'》,显示series不能转换为float,请问下各位大佬这个问题怎么解决 需要怎么修改脚本。

import numpy as np
import pandas as pd 

test_data = pd.DataFrame(pd.read_excel("test.xlsx"))

def haversine(lon1, lat1, lon2, lat2):   # 经度1,纬度1,经度2,纬度2 
  """函数作用:计算两个经纬度之间的距离"""
  # 将十进制度数转化为弧度
  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 * 1000

test_data["距离(米)"] = test_data.eval('@haversine(经度,纬度,Longitude,Latitude)')


"""
表内容举例:
name    经度  纬度  Longitude   Latitude
A   113.545761  23.106667   113.09461   23.39367
B   113.545761  23.106667   113.173833  23.336054
C   113.545761  23.106667   113.1799    23.43946
D   113.545761  23.106667   113.1975278 23.38591667
E   113.545761  23.106667   113.1975278 23.38591667
F   113.545761  23.106667   113.219664  23.100377
G   113.545761  23.106667   113.219664  23.100377
H   113.545761  23.106667   113.22253   23.18483
I   113.545761  23.106667   113.22253   23.18483

"""

程序运行到

3个回答

题主你的报错原因是map函数的第二个参数应该是列表,而你传入的四列数据相当于是一个二维的dataframe,这样每列就相当于一个Series,而函数接受的参数应该是四个数字组成的列表,也就是一行数据,就是说你一次性传入所有的数据是不行的。必须一行一行的计算,每行数据计算得到一个距离,针对你的问题,我们可以使用pandas的apply函数每次传一行数据进入函数,具体可以将代码改成如下(题主你的sin cos sqrt函数都没有加np.这个我不知道是为什么,应该是要加的):
import numpy as np
import pandas as pd

def haversine(df): # 经度1,纬度1,经度2,纬度2
"""函数作用:计算两个经纬度之间的距离"""
# 将十进制度数转化为弧度
# df.tolist()作用是将传入的一行数据转化为列表
lon1, lat1, lon2, lat2 = map(np.radians, df.tolist())
# haversine公式
dlon = lon2 - lon1
dlat = lat2 - lat1
a = np.sin(dlat/2)**2 + np.cos(lat1) * np.cos(lat2) * np.sin(dlon/2)**2
c = 2 * a * np.sin(np.sqrt(a))
r = 6371 # 地球平均半径,单位为公里
return c * r * 1000

test_data = pd.read_excel("test.xlsx")

这里的apply函数作用是将dataframe中每行的四个值传入函数haversine,然后计算距离,这样每行数据能返回一个距离

test_data["距离(米)"] = test_data[['经度', '纬度', 'Longitude', 'Latitude']].apply(haversine, axis=1)
print(test_data)

weixin_42653818
Tyche_BO 回复道恒渐远: 哦哦 多谢了~
3 个月之前 回复
qq_43615577
道恒渐远 回复Tyche_BO: 我查了一下math.asin()和np.arcsin()应该是一样的,都是反正弦函数
3 个月之前 回复
weixin_42653818
Tyche_BO 但我查了下numpy没有asin的函数只有一个np.arcsin,网上查了下asin和np.arcsin好像差不多,但又不确认。。。纠结啊 大佬你知道吗、?
3 个月之前 回复
weixin_42653818
Tyche_BO 今天试了下你这个方法貌似可行。但c = 2 * a * np.sin(np.sqrt(a))这里面应该是asin不是sin,
3 个月之前 回复
weixin_42653818
Tyche_BO 昨晚我尝试用apply代替了map,也把所有三角函数改为np. 但最后输出结果报错显示series是可变的,应该是顺序不能再对应原来表格了
3 个月之前 回复

看下表是不是没有去掉表头,或者里面有空值等,无法转换为float

在[lon1, lat1, lon2, lat2]前面加上
print([lon1, lat1, lon2, lat2])
看看输出什么

weixin_42653818
Tyche_BO [0 113.545761 1 113.545761 2 113.545761 3 113.545761 4 113.545761 5 113.545761 6 113.545761 7 113.545761 8 113.545761 Name: 经度, dtype: float64, 0 23.106667-------复制一列出来一共4列就没截了但分布都显示那4列经纬度。我怀疑 map(radians, [lon1, lat1, lon2, lat2])这里根本不是把这里面的列值转化为弧度(float)二手把这一列series对象本身尝试转换为float所以失败
3 个月之前 回复

首先你的radians函数的调用,应该是np.radians 吧
第二,输出一下你读入的数据看看具体格式是什么样的,判断一下格式对不对 test_data.head()

Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!