shapelyのワーニング：ShapelyDeprecationWarningの対処がわからない。

Question

前提・実現したいこと

pandasのdataframeからGeoSeriesを作成する際にShapelyDeprecationWarningとワーニングが出るのだけど対処がわからない。

発生している問題・エラーメッセージ

C:\Program Files\Python310\lib\site-packages\pandas\core\dtypes\cast.py:118: ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.

該当のソースコード

import pandas as pd
from shapely.geometry import Point
dfb = pd.DataFrame(
  [
    [35.6585840, 139.7454316], 
    [35.7100069, 139.8108103], 
    [35.6842961, 139.8647716], 
  ], columns = ['緯度', '経度']
)
series_geom_b=dfb.apply(lambda x:Point(x['経度'],x['緯度']), axis=1)

もしかして。

以下であればワーニングは出ないですが、これでGeoDataFrame作成しろということでしょうか。

list_geom=[]
for index, row in dfb.iterrows():
    list_geom.append(Point(row['経度'],row['緯度']))

補足情報（FW/ツールのバージョンなど）

Shapely 1.8.0

Answer 1

調べた結果、ワーニングが出なくなった他の方法等を自己回答しておきます。

# ワーニング出ない
series_geom_b=dfb[['経度','緯度']].apply(lambda x:(x[0],x[1]),axis=1).apply(Point)

# 以下ではワーニングが出る
list_geom=[]
for index, row in dfb.iterrows():
    list_geom.append(Point(row['経度'],row['緯度']))

df=pd.DataFrame({'geom':list_geom})

# こちらでもワーニングが出る
sr=pd.Series(list_geom)

#　中を追っていくと以下のような感じでワーニングだった
import numpy as np
result = np.empty(len(list_geom), dtype="object")
result[:] = list_geom

# こちらはエラーは出なかった。
import geopandas as gpd
gdf=gpd.GeoDataFrame(geometry=list_geom)

よくわかってません…

Answer 2

ワーニングメッセージには以下の様に表示されています。

ShapelyDeprecationWarning: The array interface is deprecated and will no longer work in Shapely 2.0. Convert the '.coords' to a numpy array instead.
arr = construct_1d_object_array_from_listlike(values)

そこで、Convert 〜 の通りにしますとワーニングは表示されなくなるのですが、これは求める結果ではないと思われます。

>>> pd.__version__
'1.4.0'
>>> series_geom_b=dfb.apply(lambda x:np.array(Point(x['経度'],x['緯度']).coords), axis=1)
>>> series_geom_b
0     [[139.7454316, 35.658584]]
1    [[139.8108103, 35.7100069]]
2    [[139.8647716, 35.6842961]]
dtype: object
>>> type(series_geom_b)
<class 'pandas.core.series.Series'>

まぁ、以下の様にしてもよいかと思います(何だか奇妙ですが)。

>>> series_geom_b = dfb[['経度', '緯度']].apply(lambda x: [Point(x)], axis=1).str[0]
>>> series_geom_b
0     POINT (139.7454316 35.658584)
1    POINT (139.8108103 35.7100069)
2    POINT (139.8647716 35.6842961)
dtype: object
>>> type(series_geom_b)
<class 'pandas.core.series.Series'>

Comments

[xail2222]

何故ワーニングが出なくなったのかわかりませんが、ありがとうございました。