【C++】 緯度経度より距離を求める

はじめに

• 緯度経度より距離を求める方法・備忘録 ｜ 地球は丸い！

上記のサイトを参考というか、このサイトではPHPで書かれているのをC++に書き換えてみました。では、どうぞ↓

ソースコード

distanceMeasurement.cpp

#include <iostream>
#include <math.h>

#define rad2deg(a) ((a)/M_PI * 180.0) /* rad を deg に換算するマクロ関数 */
#define deg2rad(a) ((a)/180.0 * M_PI) /* deg を rad に換算するマクロ関数 */

using namespace std;

typedef struct{
	char positionName[120];
	double latitude;
	double longitude;
}coordinate;

int main(){
	double earth_r = 6378.137;
	double laRe,loRe,NSD,EWD,distance;
	coordinate A;
	coordinate B;
	cout<<"１つ目の場所についての情報を入れてください。"<<endl;
	cout<<"地名：";
	cin>>A.positionName;
	cout<<"緯度：";
	cin>>A.latitude;
	cout<<"経度：";
	cin>>A.longitude;
	cout<<"２つ目の場所についての情報を入れてください。"<<endl;
	cout<<"地名：";
	cin>>B.positionName;
	cout<<"緯度：";
	cin>>B.latitude;
	cout<<"経度：";
	cin>>B.longitude;
	laRe = deg2rad(B.latitude - A.latitude);
	loRe = deg2rad(B.longitude - A.longitude);
	NSD = earth_r*laRe;
	EWD = cos(deg2rad(A.latitude))*earth_r*loRe;
	distance = sqrt(pow(NSD,2)+pow(EWD,2));
	cout<<A.positionName<<"と"<<B.positionName<<"の距離は"<<distance<<"km"<<endl;
}

実行結果

１つ目の場所についての情報を入れてください。
地名：大阪駅
緯度：34.702113
経度：135.494807
２つ目の場所についての情報を入れてください。
地名：東京駅
緯度：35.681541
経度：139.767103
大阪駅と東京駅の距離は405.91km

補足

座標はGoogleMapを利用して取得しました。
今回取得した距離はおおよその直線距離で実際の距離とは異なります。
精度を高めるにはもっと複雑な計算が必要ですが今回は簡易な方法で行いました。