GEO - Functions

[wurzelzwerg]

Super!

Wie würde man einer Koordinate x Meter addieren/subtrahieren?
(also nur Lat +- x Meter oder Lon +- x Meter. Nicht diagonal)
Bin leider eine absolute Mathe-Null.

Gruss
Udo

[jfheins]

Hi,
also du kannst die Erdoberfläche lokal als Ebene annähern.
Das geht aber immer nur in der Umgebung von einem Bestimmten Punkt!
Da alle Längengrade gleich lang sind, entspricht ein Grad in Nord-Süd-Richtung ca. 111,7km.
Bei dem Breitengrad ist das etwas schwieriger:
Nehmen wir z.B. Hamburg: 53,55° N, 10° O.
An dieser Stelle entspricht ein Grad in Ost-West-Richtung ungefähr 67,2km.
In Freiburg im Breisgau sind es schon 74,7 km.

Das errechnet man wie folgt:

Zitat:

Bogenlänge von einem Grad geographischer Länge: 2*PI*6400km/360° = 111,7 km/°
Bogenlänge eines Grads geographischer Breite: 2*PI*6400km * cos(Breitengrad)/360°

Ich hoffe das war jetzt verständlich. Ich weiß leider nicht, was du schon für Mathe-Vorkenntnisse hast

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

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