GEO - Functions

[wurzelzwerg]

Danke,
habe auch diesen Java-Code gefunden:

Code:

double lat1 = (52.523405d * Math.PI) / 180.0d;
double lon1 = (13.411400d * Math.PI) / 180.0d;
      
double distance = (150000.0d) / 6378388.0d;
double winkel = (208.78d * Math.PI) / 180.0d;
      
double lat2 = Math.asin(Math.sin(lat1) * Math.cos(distance) + Math.cos(lat1) * Math.sin(distance) * Math.cos(-winkel));
double lon2;
if(Math.cos(lat1) == 0){
   lon2 = lon1;
}else{
   lon2 = ((lon1 - Math.asin(Math.sin(-winkel) * Math.sin(distance) / Math.cos(lat1)) + Math.PI) % (2*Math.PI)) - Math.PI;
}
lat2 = (lat2 * 180.0d) / Math.PI;
lon2 = (lon2 * 180.0d) / Math.PI;
      
log.debug("Lat: " + lat2 + " Lon: " + lon2);

Soweit ok, nur wie setzt man das hier in Dephi um:

Code:

lon2 = ((lon1 - Math.asin(Math.sin(-winkel) * Math.sin(distance) / Math.cos(lat1)) + Math.PI) % (2*Math.PI)) - Math.PI;

Das hier funktioniert nicht:

Delphi-Quellcode:

lon2:= lon1 - arcsin(sin(-winkel) * sin(distance) / cos(lat1)) + PI;
lon2:= (lon2 - (Trunc(lon2 / (2*PI)) * (2*PI))) - PI;  // falsch

Was macht "%" in Java?

[jfheins]

Hi,

Zitat von wurzelzwerg:

Was macht "%" in Java?

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

[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 ...