关于yf2003测绘方位角的计算公式

mycad

yfy2003和gyl的 测绘方位角计算很方便，首先感谢一下他的无私奉献精神。但是公式不完整；下面是本人使用C#计算的代码；公式引用 yfy2003 和gyl的。感谢他们。

              // yfy2003

                //Azimuth = Math.Acos(D_21x / Dis) * sgn(D_21y)*180/Math.PI  + 360;
                //if (Azimuth>360)
                 // Azimuth=Azimuth-360;

//gyl

Azimuth = 180 - (90 * sgn(D_21y)) - Math.Atan(D_21x / D_21y)*180/Math.PI  ;
            this.textBox5.Text = Azimuth.ToString();

        }

        private int sgn(double n)
        {
            if (n > 0) return 1;
            else if (n < 0) return -1;
            return 0;  

             //Sign(y)   (y>=0?1:-1)  

        }

经过测试结果均正确