Вычисление координаты по углу и расстоянию от координаты

Question

Я написал этот метод расширения для IPoint.

public static IPoint Offset(this IPoint point, double angle, double distanceInMeters) 
{ 
    var radians = Math.PI * angle/180; 
    var distanceX = distanceInMeters * Math.Cos(radians); 
    var distanceY = distanceInMeters * Math.Sin(radians); 
    var earthRadius = 6371000; 

    var y = point.Y + ((distanceY/earthRadius) * 180/Math.PI); 
    var x = point.X + ((distanceX/(earthRadius * Math.Cos(y * 180/Math.PI))) * 180/Math.PI); 


    return new Point(x, y); 
} 

Он отлично работает, когда я ставлю в угол 0, 90, 180 и 270, то возвращают координаты в заданном расстоянии от начальной точки. Но когда я начинаю идти под углом, который не указывает точно Север, Восток и т. Д. Я получаю неправильные расстояния.

Где я могу пойти не так? Альтернатива Есть некоторые библиотеки для использования?

Answer 1

Попробуйте эту формулу. Поскольку я считаю, что вы должны перевести широту и долготу в радианы, а затем вернуться к градусам.

public static Point Offset(this Point point, double angle, double distanceInMeters) 
    { 
     double rad = Math.PI * angle/180; 

     double xRad = Math.PI * point.X/180; // convert to radians 
     double yRad = Math.PI * point.Y/180; 

     double R = 6378100; //Radius of the Earth in meters 
     double x = Math.Asin(Math.Sin(xRad) * Math.Cos(distanceInMeters/ R) 
           + Math.Cos(xRad) * Math.Sin(distanceInMeters/ R) * Math.Cos(rad)); 

     double y = yRad + Math.Atan2(Math.Sin(rad) * Math.Sin(distanceInMeters/ R) * Math.Cos(xRad), Math.Cos(distanceInMeters/ R) - Math.Sin(xRad) * Math.Sin(x)); 

     x = x * 180/Math.PI; // convert back to degrees 
     y = y * 180/Math.PI; 

     return new Point(x, y); 
    }