用C++求地球两点间距离

问题描述：

用C++求地球两点间距离
已知地球的平均半径为6371.393千米,假设在地球的同一纬度上,有两个处于不同经度的点A、B,编程序求出这两点之间的直线距离（即弦长）、弦长所对球心角的大小（单位为度）以及球面距离（即劣弧的长度）.
要求：
⑴计算结果要有尽可能高的精确度.
⑵ 要考虑到横跨东、西两个半球的情况.（提示：东经为正,西经为负；北纬为正,南纬为负.）
⑶ 要求输入数据之间以空格分隔.例如,北纬30度上东经10度到东经50度的距离,输入格式应为：30 10 50

答

根据你的题目写的

#include <iostream>
#include <string>
using namespace std;
#define EARTH_RADIUS 6371.393
#define PI 3.14159
int main(int argc,char ** argv)
{
\x05while(1)
\x05{
\x05\x05char szInput[128];
\x05\x05string strInput;
\x05\x05string strTmp;
\x05\x05char szTmp[128];
\x05\x05string::size_type Index;
\x05\x05float LatLongitude[3]={0};\x05//经纬度
\x05\x05float Radius;
\x05\x05float ArcLen;\x05\x05//弧长
\x05\x05float ChordLen;\x05\x05//弦长
\x05\x05float LonDiff;\x05//圆心角
\x05\x05int i=0;
\x05\x05cout<<"-----------------------------------------------"<<endl;
\x05\x05cout<<"请输入经纬度"<<endl;
\x05\x05cin.clear();
\x05\x05cin.sync();
\x05\x05cin.get(szInput,sizeof(szInput));
\x05\x05strInput=szInput;
\x05\x05while((Index=strInput.find(" ")) != string::npos)
\x05\x05{
\x05\x05\x05strTmp=strInput.substr(0,Index);
\x05\x05\x05_snprintf(szTmp,sizeof(szTmp),"%s",strTmp.c_str());
\x05\x05\x05sscanf(szTmp,"%f",&(LatLongitude[i]));
\x05\x05\x05strInput=strInput.substr(Index+1);
\x05\x05\x05i++;
\x05\x05}
\x05\x05_snprintf(szTmp,sizeof(szTmp),"%s",strInput.c_str());
\x05\x05sscanf(szTmp,"%f",&(LatLongitude[i]));
\x05\x05if(LatLongitude[0]>90 || LatLongitude[0]<-90)
\x05\x05{
\x05\x05\x05cout<<"输入纬度有误"<<endl;
\x05\x05\x05continue;
\x05\x05}
\x05\x05if(LatLongitude[1]<-180 || LatLongitude[1]>180 || LatLongitude[2]<-180 || LatLongitude[2]>180)
\x05\x05{
\x05\x05\x05cout<<"输入经度有误"<<endl;
\x05\x05\x05continue;
\x05\x05}
\x05\x05Radius=EARTH_RADIUS*sin((90-LatLongitude[0])/180*PI);
\x05\x05if(LatLongitude[2]>LatLongitude[1])
\x05\x05{
\x05\x05\x05LonDiff=LatLongitude[2]-LatLongitude[1];

\x05\x05}
\x05\x05else
\x05\x05{
\x05\x05\x05LonDiff=LatLongitude[1]-LatLongitude[2];
\x05\x05}\x05
\x05\x05if(LonDiff>180)
\x05\x05{
\x05\x05\x05LonDiff=360-LonDiff;
\x05\x05}
\x05\x05ChordLen=Radius*sin(LonDiff/2/180*PI)*2;
\x05\x05ArcLen=2*Radius*PI*(LonDiff/360);
\x05\x05cout<<"两点之间的弦长:"<<ChordLen<<"千米"<<endl;
\x05\x05cout<<"弦长所对应的圆心角:"<<LonDiff<<"度"<<endl;
\x05\x05cout<<"两点之间的劣弧长:"<<ArcLen<<"千米"<<endl;

\x05}
\x05
\x05return 0;
}

运行

相关推荐