azimuth.cpp

来自「C 程序 球面方位角算法： azimuth.cpp CRC校验快速算法： C」· C++ 代码 · 共 15 行

//－－已知Ａ(Ja,Wa)、Ｂ(Jd,Wd)两点，求Ａ测Ｂ的球面方位角[0,360)－－－－－－－－－－－－－－－－－

double trua(double Ja, double Wa, double Jd, double Wd)
{
     int m;
     double ta, Lx;
     Lx=atan(cos(Jd-Ja)*tan(Wa));
     m=Lx>Wd?1:0;
     ta=m*M_PI+pow(-1,m)*asin(cos(Wd)*sin(Jd-Ja)/sin(acos(sin(Wa)*sin(Wd)
        +cos(Wa)*cos(Wd)*cos(Jd-Ja))));
     ta= ta>=0?ta:ta+2*M_PI;
     return(ta);
}