Это похоже на C #.
Сначала нужно определить toRadians
и toDegrees
:
double toRadians(double degrees) {
double sign = Math.Sign(degrees);
while(Math.Abs(degrees) > 360) {
degrees -= sign * 360;
}
return Math.PI * degrees / 180;
}
double toDegrees(double radians) {
double sign = Math.Sign(radians);
while(Math.Abs(radians) > 2 * Math.PI) {
radians -= sign * 2 * Math.PI;
}
return 180 * radians / Math.PI;
}
Затем для использования тригонометрических функций необходимо использовать Math.Sin
, Math.Cos
и т. Д.
double dist = Math.Sin(lat1) * Math.Sin(lat2)
+ Math.Cos(lat1) * Math.Cos(lat2) * Math.Cos(theta);
и
dist = toDegrees(Math.Acos(dist)) * 60 * 1.1515 * 1.609344 * 1000;
Комментарии:
public static double distance (double lat1, double lon1, double lat2, double lon2) {
double lat1 = Convert.ToDouble(latitude);
double lon1 = Convert.ToDouble(longitude);
double lat2 = Convert.ToDouble(destlat);
double lon2 = Convert.ToDouble(destlon);
Что это? Где определены latitude
, longitude
, destlat
и destlon
? Кроме того, кажется, что у вас есть lat1
, lon1
lat2
и lon2
в качестве параметров для этого метода, так что вы не можете определить здесь локальных пользователей с тем же именем.
double theta = toRadians(lon1-lon2);
lat1 = toRadians(lat1);
lon1 = toRadians(lon1);
lat2 = toRadians(lat2);
lon2 = toRadians(lon2);
Это плохой стиль. Если lat1
представляет широту в градусах, гораздо лучше вычислить эквивалентное радианам значение lat1
следующим образом:
double lat1Radians = toRadians(lat1);
Замените вышеприведенное на:
double theta = toRadians(lon1-lon2);
double lat1Radians = toRadians(lat1);
double lon1Radians = toRadians(lon1);
double lat2Radians = toRadians(lat2);
double lon2Radians = toRadians(lon2);
И наконец:
double dist = sin(lat1) * sin(lat2)
+ cos(lat1) * cos(lat2) * cos(theta);
dist = toDegrees(acos(dist)) * 60 * 1.1515 * 1.609344 * 1000;
Это тоже плохой стиль. Первая формула и вторая формула не могут одновременно представлять расстояние, которое вы пытаетесь вычислить. Вы должны присвоить результат первой формулы переменной с более значимым именем. В худшем случае, по крайней мере, сделайте следующее:
double temp = Math.Sin(lat1) * Math.Sin(lat2)
+ Math.Cos(lat1) * Math.Cos(lat2) * Math.Cos(theta);
double dist = toDegrees(Math.Acos(dist)) * 60 * 1.1515 * 1.609344 * 1000;
return dist;