🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

5034閲覧

geopandas.sjoinの仕様を知りたい

xail2222

総合スコア1508

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2021/03/11 14:43

編集2021/03/12 23:03

前提・実現したいこと

geopandas.sjoinのパラメータに設定できる内容を知りたい。
特にopパラメータと許容誤差は指定できるのかについて

ネットでGithubを見てみたら
opは{None, 'intersects', 'within', 'contains', 'overlaps', 'crosses', 'touches'}ぽい?

sjoinは空間結合ということでいいと思うのだけどその説明は以下を参考に
リンク

該当のソースコード

sjoinの使用例

python

1from shapely.geometry import LineString 2from shapely.geometry import Point 3import geopandas as gpd 4 5lines= [LineString([(0, 0), (5, 0), (5, 5), (10, 5)]), 6 LineString([(5, 0), (30, 0), (30, 5), (0, 5)])] 7 8lines2= [LineString([(0, 0), (6, 0.1), (5, 5), (10, 5)]), 9 LineString([(5, 0), (30, 0), (30, 5), (0, 5)])] 10 11pnts = [Point(6, 0.1),Point(30,2.5)] 12 13gdfl=gpd.GeoDataFrame(geometry=lines) 14gdfl2=gpd.GeoDataFrame(geometry=lines2) 15gdfp=gpd.GeoDataFrame(geometry=pnts) 16 17gdfsj=gpd.sjoin(gdfp,gdfl) 18 19for index,row in gdfsj.iterrows(): 20 print(str(index)+ ':' +str(row['index_right'])) 21 print(row['geometry']) 22 23print('**') 24 25gdfsj2=gpd.sjoin(gdfp,gdfl2) 26 27for index,row in gdfsj2.iterrows(): 28 print(str(index)+ ':' +str(row['index_right'])) 29 print(row['geometry']) 30

上記の出力結果は
1:1
POINT (30 2.5)
**
0:0
POINT (6 0.1)
1:1
POINT (30 2.5)

対処案

許容誤差を指定できないならバッファを取って結合する?

python

1def sjoin2(gdfl,gdfr,distance=0): 2 if distance==0: 3 return gpd.sjoin(gdfl,gdfr) 4 else: 5 dfl=gdfl.drop(columns='geometry') 6 gdfl2=gpd.GeoDataFrame(dfl,geometry=gdfl.buffer(distance)) 7 gdfsj=gpd.sjoin(gdfl2,gdfr).drop(columns='geometry') 8 return gdfsj.merge(gpd.GeoDataFrame(geometry=gdfl['geometry']),left_index=True,right_index=True,sort=False)

…数倍遅い感じですかね…
許容誤差を加味して結合するのはどのようにするのが良いのでしょうか。

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

Python 3.9.2
geopandas 0.8.2

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

自己解決

Sjoin Arguments

書いてあった。リンクしてあったのでたどれば見つかる感じでした…

その仕様が

The Shapely User Manual-Binary Predicates

に書いてあった。

以下であってるっぽい。
intersects
contains
within
touches
crosses
overlaps

また、許容誤差については、やはりバッファを作って処理すると書いてありました。

以上です。

投稿2021/05/26 23:16

xail2222

総合スコア1508

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

xail2222

2021/12/19 17:09

geopandas.GeoDataFrame.sjoin_nearestがバージョン0.10.?から追加になったみたい。許容誤差を指定して結合できそうな感じです
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問