Wie gesagt, ein winziger Umrechnungsfaktor, Anpassung der Ellipse, sollte noch mit rein.
Kann man sich bestimmt für die Mitte seines Bereichs genau ausrechnen,
aber so für Deutschland kannst etwa 2/3 nehmen.
Das sind ja fast Quadrate, aber in der Breite 3° und in der Höhe 2°.
https://www.mapsofworld.com/lat_long...-lat-long.html
Bielefeld-Magdeburg (3° Lon) müsste fast so weit sein, wie Bielefeld-Mainz (2° Lat).
d := sqrt(pow(lat1-lat2, 2) + pow((lon1-lon2) * 0.66, 2)); // ich hoffe es ist an der richtigen Stelle drin