输入必须是数组、列表、元组或标量 pyproj

Question

我有一个 DF，我正在尝试将 eastings/northings 转换为 long/lats。我的 df 看起来像这样：

import pandas as pd
import numpy as np
import pyproj

    Postcode    Eastings    Northings
0   AB101AB     394235      806529
1   AB101AF     394181      806429
2   AB101AG     394230      806469
3   AB101AH     394371      806359
4   AB101AL     394296      806581

我正在使用众所周知的代码块将东距和北距转换为 long/lats 并将这些 long/lats 作为新列添加到 df:

def proj_transform(df):
    bng = pyproj.Proj("+init=EPSG:27700")
    wgs84 = pyproj.Proj("+init=EPSG:4326")
    lats = pd.Series()
    lons = pd.Series()
    for idx, val in enumerate(df['Eastings']):
        lon, lat = pyproj.transform(bng, wgs84, df['Eastings'][idx], df['Northings'][idx])
        lats.set_value(idx, lat)
        lons.set_value(idx, lon)
    df['lat'] = lats
    df['lon'] = lons
    return df

df_transform = proj_transform(my_df)

但是，我不断收到以下错误，“输入必须是数组、列表、元组或标量”。有没有人知道我哪里出错了？

Answer 1

您可以将 DataFrame.apply 与 axis=1 一起使用并更改函数，例如：

def proj_transform(x):
    e = x['Eastings']
    n = x['Northings']
    bng = pyproj.Proj("+init=EPSG:27700")
    wgs84 = pyproj.Proj("+init=EPSG:4326")
    lon, lat = pyproj.transform(bng, wgs84, e, n)

    return pd.Series([lon, lat])

my_df[['lat','lon']] = my_df.apply(proj_transform, axis=1)

Answer 2

这是最快的方法：

https://gis.stackexchange.com/a/334307/144357

from pyproj import Transformer

trans = Transformer.from_crs(
    "EPSG:27700",
    "EPSG:4326",
    always_xy=True,
)
xx, yy = trans.transform(my_df["Eastings"].values, my_df["Northings"].values)
my_df["X"] = xx
my_df["Y"] = yy

也有参考价值：

输入必须是数组、列表、元组或标量 pyproj

input must be an array, list, tuple or scalar pyproj

numpy

pandas

geopandas

pyproj