はじめに以前『火山灰追跡モデルPUFF』についてnoteを書きました。原論文を見ながらゼロからPythonで実装するという、なかなかハードで良い経験でした。 このとき火山灰の動きを計算する方法として、Pythonのpyprojライブラリを活用しました。pyprojは地球が楕円体であることを考慮して、地理的な位置情報を使った計算ができるライブラリです。 pyprojの使い方をweb検索したところ、ある2地点の緯度経度から距離を計算する方法について書かれたweb記事等はたくさん見つかりました。 例えば大阪管区気象台と東京管区気象台の距離は、以下のコードで計算できます。 import pyproj grs80 = pyproj.Geod(ellps='GRS80') # 東京管区気象台 tokyo_lat = 35.69 tokyo_lon = 139.75 # 大阪管区気象台 osaka_la
![Pythonで出発点と進む距離・方向から到着点の緯度経度を計算する方法|加藤芳樹 (Yoshiki Kato)](https://cdn-ak-scissors.b.st-hatena.com/image/square/f4132360b16a6da8b9f729ca168ee96bc917ba02/height=288;version=1;width=512/https%3A%2F%2Fassets.st-note.com%2Fproduction%2Fuploads%2Fimages%2F88277064%2Frectangle_large_type_2_3ff4370440371a1f83acf30a443fb5de.png%3Ffit%3Dbounds%26quality%3D85%26width%3D1280)