Punto medio geografico tra due coordinate

Question

Ho utilizzato il sito web Moveable-Type per aiutarmi in alcuni calcoli geocoordinati ed è stato molto utile, tuttavia, ho un bug nel mio calcolo del punto centrale tra due coordinate. Il mio risultato è vicino al previsto, ma non abbastanza vicino:

posA = {47.64570362, -122.14073746} 
posB = {47.64316917, -122.14032175} 

risultato atteso (preso dalla calcolatrice caratteri mobili) = 47 ° 38'40 "N, 122 ° 08'26" W = {47.644444, -122.140556} mio risultato : {49.6054801645915, -122.14052959995759}

Ecco il mio codice:

private Geocoordinate MidPoint(Geocoordinate posA, Geocoordinate posB) 
{ 
    Geocoordinate midPoint = new Geocoordinate(); 

    double dLon = DegreesToRadians(posB.Longitude - posA.Longitude); 
    double Bx = Math.Cos(DegreesToRadians(posB.Latitude)) * Math.Cos(dLon); 
    double By = Math.Cos(DegreesToRadians(posB.Latitude)) * Math.Sin(dLon); 

    midPoint.Latitude = RadiansToDegrees(Math.Atan2(Math.Sin(DegreesToRadians(posA.Latitude)) + Math.Sin(DegreesToRadians(posB.Latitude)), 
       Math.Sqrt((Math.Cos(DegreesToRadians(posA.Latitude)) + Bx) * (Math.Cos(DegreesToRadians(posA.Latitude))) + Bx) + By * By)); 

    midPoint.Longitude = posA.Longitude + RadiansToDegrees(Math.Atan2(By, Math.Cos(DegreesToRadians(posA.Latitude)) + Bx)); 

    return midPoint; 
} 

ho un paio di metodi privati ​​per fare la conversione tra gradi e radianti e schiena. E.g.

private double DegreeToRadian(double angle) 
{ 
    return Math.PI * angle/180.0; 
} 

Non riesco a capire perché i miei risultati sono spenti di un paio di gradi sul valore Lat. Qualche idea?

Grazie

Answer 1

Hai posizionato alcune parentesi errate. Ho segnato il posto nel codice.

private Geocoordinate MidPoint(Geocoordinate posA, Geocoordinate posB) 
{ 
    Geocoordinate midPoint = new Geocoordinate(); 

    double dLon = DegreesToRadians(posB.Longitude - posA.Longitude); 
    double Bx = Math.Cos(DegreesToRadians(posB.Latitude)) * Math.Cos(dLon); 
    double By = Math.Cos(DegreesToRadians(posB.Latitude)) * Math.Sin(dLon); 

    midPoint.Latitude = RadiansToDegrees(Math.Atan2(
       Math.Sin(DegreesToRadians(posA.Latitude)) + Math.Sin(DegreesToRadians(posB.Latitude)), 
       Math.Sqrt(
        (Math.Cos(DegreesToRadians(posA.Latitude)) + Bx) * 
        (Math.Cos(DegreesToRadians(posA.Latitude)) + Bx) + By * By))); 
       // (Math.Cos(DegreesToRadians(posA.Latitude))) + Bx) + By * By)); // Your Code 

    midPoint.Longitude = posA.Longitude + RadiansToDegrees(Math.Atan2(By, Math.Cos(DegreesToRadians(posA.Latitude)) + Bx)); 

    return midPoint; 
} 

Answer 2

Alcuni metodi (convenzioni ad esempio WGS84) includono schiacciamento della terra. (at least that's what it says here)