根据两点经纬度求直线与正北方向的夹角

/// <summary>
        /// 根据两点经纬度求直线与正北方向的夹角
        /// </summary>
        /// <param name="A"></param>
        /// <param name="B"></param>
        /// <returns>AB连线与正北方向的角度（顺时针0~360）</returns>
        public  static double getAngle(MyLatLng A,MyLatLng B)
        {  
            double dx=(B.m_RadLo-A.m_RadLo)*A.Ed;
            double dy=(B.m_RadLa-A.m_RadLa)*A.Ec;
            double angle=0.0;
            angle=Math.Atan(Math.Abs(dx/dy))*180/Math.PI;
            double dLo=B.m_Longitude-A.m_Longitude;
            double dLa=B.m_Latitude-A.m_Latitude;
            if(dLo>0&&dLa<=0)
            {  
                angle=(90-angle)+90;
            }  
            else if(dLo<=0&&dLa<0)
            {  
                angle=angle+180;
            }else if(dLo<0&&dLa>=0)
            {  
                angle= (90-angle)+270;
            }  
            return angle;
        }
        public class MyLatLng 
        { 
            public static double Rc=6378137;  
            public static double Rj=6356725;
            public double m_LoDeg, m_LoMin, m_LoSec;
            public double m_LaDeg, m_LaMin, m_LaSec;
            public double m_Longitude, m_Latitude;  
            public double m_RadLo,m_RadLa;
            public double Ec;
            public double Ed;  
            public MyLatLng(double longitude,double latitude)
            {  
                m_LoDeg=(int)longitude;
                m_LoMin=(int)((longitude-m_LoDeg)*60);
                m_LoSec=(longitude-m_LoDeg-m_LoMin/60)*3600;
              
                m_LaDeg=(int)latitude;
                m_LaMin=(int)((latitude-m_LaDeg)*60);
                m_LaSec=(latitude-m_LaDeg-m_LaMin/60)*3600;
              
                m_Longitude=longitude;
                m_Latitude=latitude;
                m_RadLo=longitude*Math.PI/180;
                m_RadLa=latitude*Math.PI/180;
                Ec=Rj+(Rc-Rj)*(90-m_Latitude)/90;
                Ed=Ec*Math.Cos(m_RadLa);
            }
        }