ヤマレコBLOG

GPSデータからグラフ作り：３）経度・緯度で指定した2地点間の距離を求める

PHPでGPXファイルを取得することができたので、その取得した経度・緯度座標のリストから、距離を求める必要があります。

Webで検索したところ、正式っぽい数式が国土地理院のページにありました。
2点の平面直角座標x,yから測地線長および方向角を求める計算

・・・が、さっぱり理解できない。
仕方がないので他のページも調べてみると、Javascriptのプログラムまでついているページがありました。球面三角法という計算手法を使うようです。
緯度と経度による距離計算

内容を理解し切れてないのですが、ソースコードを参考（そのままコピー・・・）して使わせていただきました。
PHP向けに作り直したコードを乗せておきます。
function cal_len($na, $ta, $nb, $tb){ $na=pi()*$na/180; $nb=pi()*$nb/180; $ta=pi()*$ta/180; $tb=pi()*$tb/180; $na=$na-((11.55/60)*pi()/180)*sin(2*$na); $nb=$nb-((11.55/60)*pi()/180)*sin(2*$nb); $c=cos($na)*cos($nb)*cos($ta-$tb)+sin($na)*sin($nb); $s=sqrt(1-$c*$c); $t=$s/$c; $z=atan($t); $z=6378.137*$z; // 地球の半径を6378.137kmで計算 return $z; }
これを使って、やっとルートのグラフを生成するプログラムができました！
↓グラフはこんな感じです。


あとは、GPSに標高データが入っていないときでもグラフが作れるように、座標から標高データを取得する方法です。できるのかな・・・。