【原】卫星方位角计算

设观测点经纬度为(Φ2，λ2)，卫星经纬度为(Φ1，λ1)

　　X=Φ2-Φ1；Y=λ2-λ1

　　az=atan[tan(abs(X))/sin(abs(Y))]，az的取值范围为-90到90度

分析：正北0度，东方90度，南方180度，西方270度

　　Y>0，观测地为北方，卫星在南方，az在90-270度；Y<0，观测地为南方，卫星在北方，az在-90到90度

　　X>0，观测地为东方，卫星在西方；X<0，观测地为西方，卫星在东方

　　假设X>0，Y>0；卫星位于西南（180-270度），az=pi+az

　　假设X<0，Y>0；卫星位于东南（90-180度），az=pi-az;

　　假设X>0，Y<0；卫星位于西北（270-360度），az=2*pi-az

　　假设X<0，Y<0；卫星位于东北（0-90度），az=az

结论

　　如果Y>0，那么公式，az=180+atan[tan(X)/sin(Y)]

　　如果Y<0，那么公式，az=atan[tan(X)/sin(Y)]，设西方270-360度为-90到0度