パンダのデータフレームにpyprojを使用して投影を変換する

Question

これは間違いなく "木の木を見ることができません"という瞬間です。私は1時間このコードを見てきて、私が間違ったことを見ることはできません。私は顔で私を見つめていることを知っているが、私はそれを見ることができない！

私は、Pythonを使用して2つの地理座標系を変換しようとしています。

Iは、経度（x軸）と緯度（y軸）の値を有し、単一の点についてOSGB 1936に変換する、私は次の操作を実行できます

import numpy as np 
import pandas as pd 
import shapefile 
import pyproj 

inProj = pyproj.Proj(init='epsg:4326') 
outProj = pyproj.Proj(init='epsg:27700') 

x1,y1 = (-2.772048, 53.364265) 

x2,y2 = pyproj.transform(inProj,outProj,x1,y1) 

print(x1,y1) 
print(x2,y2) 

これは以下を生成します：合理的なようで、-2.772048の経度が348721.0103978の座標に変換されることを示唆している

-2.772048 53.364265 
348721.01039783185 385543.95241055806 

。

実際、私はこれをPandasデータフレームで行いたいと考えています。データフレームには経度と緯度を含む列が含まれており、変換された座標を含む2つの列（newLongとnewLat）を追加したいとします。

模範データフレームは次のようになります。

latitude longitude 
0 53.364265 -2.772048 
1 53.632481 -2.816242 
2 53.644596 -2.970592 

そして、私が書いたコードは次のとおりです。

生成

import numpy as np 
import pandas as pd 
import shapefile 
import pyproj 

inProj = pyproj.Proj(init='epsg:4326') 
outProj = pyproj.Proj(init='epsg:27700') 

df = pd.DataFrame({'longitude':[-2.772048,-2.816242,-2.970592],'latitude':[53.364265,53.632481,53.644596]}) 

def convertCoords(row): 
    x2,y2 = pyproj.transform(inProj,outProj,row['longitude'],row['latitude']) 
    return pd.Series({'newLong':x2,'newLat':y2}) 

df[['newLong','newLat']] = df.apply(convertCoords,axis=1) 

print(df) 

：

latitude longitude  newLong   newLat 
0 53.364265 -2.772048 385543.952411 348721.010398 
1 53.632481 -2.816242 415416.003113 346121.990302 
2 53.644596 -2.970592 416892.024217 335933.971216 

をしかし、今、それはニューロングようですnewLatの値が混在しています（上記のシングルポイント変換の結果と比較して）。

この結果を得るために私のワイヤーを交差させた場所はありますか？ （完全に明白な場合はお詫び申し上げます）

Answer 1

df[['newLong','newLat']] = df.apply(convertCoords,axis=1)を実行すると、df.apply出力の列のインデックスが作成されます。ただし、系列が辞書（本来は順序付けられていない）を使用して定義されているため、列の順序は任意です。

あなたは固定列の順序でシリーズを返すことを選ぶことができます

：あなたがあなたの代わりに結果を組み合わせること.joinを使用することができ、convertCoords出力標識保持したい場合は、

return pd.Series([x2, y2]) 

代わり：

return pd.Series({'newLong':x2,'newLat':y2}) 
... 
df = df.join(df.apply(convertCoords, axis=1))