pythonで地図上に円を書き、プロットがその円内にあるか判定したい

Question

前提・実現したいこと

pythonでfoliumのような地図モジュールを使って、ある中心地点から半径Xmの円を描写し、
その円の中にプロット（緯度、経度）の数を数えたいです。

イメージ図　この図だと円内のプロット数は4

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

円の描写はできるのですが、プロットが円の中に入っているかどうかの判定ロジックに悩んでいます。

該当のソースコード

参考コード
python

import folium

m = folium.Map(location=[35.681382, 139.76608399999998], zoom_start=12)

folium.Marker([35.658581, 139.745433], popup='Tokyo tower', icon=folium.Icon(color='blue')).add_to(m)
folium.Marker([35.710063, 139.8107], popup='Tokyo skytree', icon=folium.Icon(color='blue', icon='bookmark')).add_to(m)

folium.CircleMarker(
    location=[35.681382, 139.76608399999998],
    radius=2000,
    popup='Tokyo Station',
    color='#3186cc',
    fill_color='#3186cc'
).add_to(m)
m

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

python(3.8.3)

Answer 1

Leaflet の距離計算メソッドとかが、Wrap されていないようで geopy などで距離計算させる必要があるようです。

• PyLearn MAP 01: 地図プログラミングの基礎

※上のリンクが適切かどうかは精査できてません。

CircleMarker の中心座標とプロットの座標から、geopyのdistance で距離を計算し、CircleMarker の半径より小さいかどうかで判定するようなイメージになるでしょうか。

【以下、folim ではうまく使えませんでした。】

Circle の getBounds() メソッドで LatLngBounds を生成し、contains(<LatLng> latlng) でプロットの緯度経度を渡せば、含まれているかどうかの判定ができそうです。Leaflet 自体のドキュメントを見て話をしているので folium の方で呼べるか未検証ですが。

参考： Documentation - Leaflet - a JavaScript library for interactive maps

Comments

[shin33]

わざわざ検証までいただきありがとうございます。 getBounds()まではラッパーされてないのでしょうか。私の方でもドキュメント読んでみます。

[Yasumichi]

getBouds は get_bounds で呼び出せるのですが、結果はただのリストになってしまうので contains が呼び出せませんでした。

なんか、元々の Leaflet でも四角で計算してるようなので使えそうにありません。

https://github.com/substack/point-in-polygon

というプラグインがベースにしたアルゴリズムもリンクが切れているし。
とりあえずは、aokikenichi さんの方法で考えるしかないかもしれません。

[Yasumichi]

geopy で距離計算する例を見つけたので追記しました。

Answer 2

地球規模ならば地理座標系を意識する必要がありますが、このサイズならば平面近似で全く問題ないのでは？

円を描いているならば中心座標はあるはずで

中心座標と点xの距離　<　円の半径　ならば　円内
中心座標と点xの距離　=　円の半径　ならば　円周上
中心座標と点xの距離　>　円の半径　ならば　円の外