GEO - Functions

[wurzelzwerg]

Zitat von jfheins:

Das ist der modulo-Operator. "mod" unter Delphi.

Ja, klar. Nur in Delphi geht mod nur mit Ganzzahlen.

Ich hab das jetzt so gemacht:

Delphi-Quellcode:

procedure AddDistance(var Lat, Lon: Extended; direction, Distance: Extended);
var
  Lat1, Lon1, dist, winkel, Lat2, Lon2: Extended;
begin
  Lat1:= (Lat * PI) / 180;
  Lon1:= (Lon * PI) / 180;
  dist:= Distance / 6378388;
  winkel:= (direction * PI) / 180;
  Lat2:= arcsin(Sin(Lat1) * Cos(dist) + Cos(Lat1) * Sin(dist) * Cos(-winkel));
  if (Cos(Lat1) = 0) then
    Lon2:= Lon1
  else
  begin
     lon2:= lon1 - arcsin(sin(-winkel) * sin(dist) / cos(lat1)) + PI;
     lon2:= (lon2 - (Trunc(lon2 / (2*PI)) * (2*PI))) - PI;
  end;
  Lat:= (Lat2 * 180.0) / PI;
  Lon:= (Lon2 * 180.0) / PI;
end;

[wurzelzwerg]

DecimalLon, DecimalLat haben einen Fehler.
Man muss bei W und S das Ergebnis negieren. So wie das jetzt ist wird nur der Anfang(Grad) negiert.

8°28''58.3"W ergibt jetzt -7,48286151885986. Das sollte aber -8,48286151885986 sein.

Delphi-Quellcode:

 
//if (Pos('°W', Longitude) > 0) or (Longitude[length(Longitude)] = 'W') then
//   Longitude := '-' + Longitude;
....
// erst am Ende prüfen
if (Pos('°W', Longitude) > 0) or (Longitude[length(Longitude)] = 'W') then
  Result := -Result;

[QuickAndDirty]

@WurzelZwerg
Bist du dir sicher?
Weil Longitude hat doch am Ende der Funktion gar keine 'W' '°' Zeichen mehr oder?

[Luckie]

Ich finde deine Parameterreihenfolge etwas seltsam: (Lat1, Lat2, Lon1, Lon2: Extended) . Man muss erst die Breite der ersten Koordinate angeben, dann die Breite der zweiten usw. Warum nicht: (Lat1, Lon1, Lat2, Lon2: Extended) , also paarweise, wie man sie auf dem Papier in der Regel auch vorliegen hat?

[Rollo62]

Zitat von wurzelzwerg:

Zitat von jfheins:

Das ist der modulo-Operator. "mod" unter Delphi.

Ja, klar. Nur in Delphi geht mod nur mit Ganzzahlen.

Wo ein Wille ist ...