用pymap3d计算熊猫列的有效方法？

Question

我有一个带有地理坐标列的 pandas 数据框，我正在使用 pymap3d 将位置转换为其他坐标系。我实现的一个典型功能是：

def append_enu(df, observerlla):
e = []
n = []
u = []
for index, row in df.iterrows():
    tmpe, tmpn, tmpu = pm.geodetic2enu( row["lat_deg"], row["lon_deg"], row["alt_m"], *observerlla )
    e.append(tmpe)
    n.append(tmpn)
    u.append(tmpu)
df["enu_e_m"] = e
df["enu_n_m"] = n
df["enu_u_m"] = u
return df    

这可行，但我发现它非常慢。 （我的表有超过 700000 行，我正在为 3 个不同的坐标系添加转换。）是否有一种“更 Pythonic”的方式来正确利用 Pandas 允许数据框操作的方式？

Answer 1

这是我想出的解决方案，它修改了原始 DataFrame 而不是从函数返回一个新的：

(df['enu_e_m'], df['enu_n_m'], df['enu_u_m']) = pm.geodetic2enu(df['lat_deg'], df['lon_deg'], df['alt_m'], *observerlla)

所以，是的，“它确实有效。”