Bearing and Reverse bearing about projnet4geoapi HOT 3 OPEN

/// <summary>
		/// 
		/// </summary>
		public enum CompassDirection
        {
			N,
			NE,
			E,
			SE,
			S,
			SW,
			W,
			NW
        }

		/// <summary>
		/// 
		/// </summary>
		/// <param name="bearing"></param>
		/// <returns></returns>
		[System.Runtime.CompilerServices.MethodImpl(System.Runtime.CompilerServices.MethodImplOptions.AggressiveInlining)]
		public static CompassDirection GetCompassDirection(double bearing)
        {
			if (bearing >= 0 && bearing <= 22.5)
			{
				return CompassDirection.N;
			}
			else if (bearing >= 22.5 && bearing <= 67.5)
			{
				return CompassDirection.NE;
			}
			else if (bearing >= 67.5 && bearing <= 112.5)
			{
				return CompassDirection.E;
			}
			else if (bearing >= 112.5 && bearing <= 157.5)
			{
				return CompassDirection.SE;
			}
			else if (bearing >= 157.5 && bearing <= 205.5)
			{
				return CompassDirection.S;
			}
			else if (bearing >= 202.5 && bearing <= 247.5)
			{
				return CompassDirection.SW;
			}
			else if (bearing >= 247.5 && bearing <= 292.5)
			{
				return CompassDirection.W;
			}
			else if (bearing >= 292.5 && bearing <= 337.5)
			{
				return CompassDirection.NW;
			}
            else
            {
				return CompassDirection.N;
            }
		}
[System.Runtime.CompilerServices.MethodImpl(System.Runtime.CompilerServices.MethodImplOptions.AggressiveInlining)]
		public static double DegreeBearing(double lon1, double lat1, double lon2, double lat2)
		{
			var dLon = ToRadians(lon2 - lon1);
			var dPhi = Math.Log(Math.Tan(ToRadians(lat2) / 2 + Math.PI / 4) / Math.Tan(ToRadians(lat1) / 2 + Math.PI / 4));
			if (Math.Abs(dLon) > Math.PI) dLon = dLon > 0 ? -(2 * Math.PI - dLon) : (2 * Math.PI + dLon);
			return ToBearing(Math.Atan2(dLon, dPhi));
		}


[System.Runtime.CompilerServices.MethodImpl(System.Runtime.CompilerServices.MethodImplOptions.AggressiveInlining)]
		public static double ToBearing(double radians)
		{
			// convert radians to degrees (as bearing: 0...360)
			return (ToDegrees(radians) + 360) % 360;
		}