• Willkommen im Geoclub - dem größten deutschsprachigen Geocaching-Forum. Registriere dich kostenlos, um alle Inhalte zu sehen und neue Beiträge zu erstellen.

Kreis und sekant schnittpunkt berechnen mit latitude longitu

ts.tsogoo

Geonewbie
Hallo Leute,

Ich versuche seit 3 Tagen das Problem zu Loesen. Ohne Erfolg. Hat jemand irgend ne Ahnung wie ich das Loesen kann.
Es geht um ein Dienswagen Verwaltungs System dabei ein klein Modul fuer durchschnittliche Geschwindigkeit Berechner. Dafuer habe ich das Ort in Kreis Zonen geteilt.
Man faehrt mit unterschiedlicher Geschwindigkeit fuer jeder Zone. Das heisst:
Ich fahre A bis D mit Geschwindigkeit v1 und D bis B mit v2. Die v1 und v2 sind mir bekannt.
Die Entfernung zwischen zwei Koordinaten kann ich berechnen. Das heisst Laenge von AB, AC, CB sind mir auch bekannt.

Die Koordinaten von Google Map macht das ganze fuer mich schwerer. Ansonsten es ist ein Kreis, Sekant Schnittpunkt Berechnung.

Ich bin fuer jeden Tipp dankbar.

Danke voraus!
howTo.png

R = 1 km
C = {48.137024,11.575249}
A = {48.139115,11.578081}
B = {48.146303,11.593102}
D = {?,?}
 

ColleIsarco

Geowizard
Moin moin,
ts.tsogoo schrieb:
Hallo Leute,

Ich versuche seit 3 Tagen das Problem zu Loesen. Ohne Erfolg. Hat jemand irgend ne Ahnung wie ich das Loesen kann.
[...]
Die Koordinaten von Google Map macht das ganze fuer mich schwerer. Ansonsten es ist ein Kreis, Sekant Schnittpunkt Berechnung.
Bei derart kleinen Entfernungen würde ich einfach die Korrdinaten in den Euklidischen Raum übertragen, sprich UTM- oder GK-Koordinaten verwenden und dann einfach den Schnittpunkt berechnen. Also: Kreis- und Geradengleichung aufstellen und losrechnen.

Gruß
ColleIsarco
 

Kalleson

Geowizard
ColleIsarco schrieb:
Moin moin,
ts.tsogoo schrieb:
Hallo Leute,

Ich versuche seit 3 Tagen das Problem zu Loesen. Ohne Erfolg. Hat jemand irgend ne Ahnung wie ich das Loesen kann.
[...]
Die Koordinaten von Google Map macht das ganze fuer mich schwerer. Ansonsten es ist ein Kreis, Sekant Schnittpunkt Berechnung.
Bei derart kleinen Entfernungen würde ich einfach die Korrdinaten in den Euklidischen Raum übertragen, sprich UTM- oder GK-Koordinaten verwenden und dann einfach den Schnittpunkt berechnen. Also: Kreis- und Geradengleichung aufstellen und losrechnen.

Gruß
ColleIsarco

Wie? Was? Die Erde ist keine Scheibe?

... was ich eigentlich sagen wollte: Collesarco hat recht. Die Abweichung bei den gegebenen Daten würfte vernachlässigbar sein.
 

moenk

Administrator
Teammitglied
Ist das ein Mystery? Wer weiß ob der Owner das richtig gerechnet hat.
Achja:
INiem.png
 
Intervall-Schachtelung mit Mopsos.

Code:
program SekanteKreis;

VAR
   LatA, LonA: real;
   LatB, LonB: real;
   LatC, LonC: real;
   LatD, LonD: real;

   W, R, Rd: real;
   DistAB: real;
   Inside: integer;
   F: real;
   Teiler: real;
   Fehler: real;
begin
   R := 1000;
   LatC := 48.137024;
   LonC := 11.575249;
   LatA := 48.139115;
   LonA := 11.578081;
   LatB := 48.146303;
   LonB := 11.593102;

   if CalcProjection(LatA, LonA, LatB, LonB, DistAB, W) then begin
      Writeln(DistAB:10:2);
      Teiler:=2;
      F := 1 / Teiler;
      repeat
         MakeProjection(LatA, LonA, DistAB * F, W, LatD, LonD);
         Rd := Distance(latC, lonC, latD, lonD);

         Fehler := RD - R;
         write(Fehler:10:2, Teiler:10:0);

         Teiler:=Teiler*2;
         if Fehler < 0 then begin
            F := F + 1 / Teiler
         end
         else begin
            F := F - 1 / Teiler
         end;

         writeln(' ',WGS84(LatD, LonD))
      until Abs(Fehler) < 0.1
   end;

end.

Der macht es dir auf dem Ellipsoid.
 

toc-rox

Geomaster
Identische Crosspostings in mehrere Foren finde ich nicht so gut ...
Siehe auch hier: http://forum.openstreetmap.org/viewtopic.php?id=18742

Gruß Klaus
 

S-Man42

Geomaster
Ich CrossPoste mal die Lösung:

Hi, wenn du nur die Lösung brauchst:
D müsste bei 48.14273015/11.58563466 liegen.

Ich habe hierfür GCC für Android verwendet, die Koordinatenberechnungen der Form ausführen kann.
 
Ah, wenn schon Lösungen gepostet werden werde ich mich mal anschließen:

Code:
   1374.19
     -3.99         2 N 48° 8.5626 E 11° 35.1355 48.14270925 11.58559098
    338.99         4 N 48° 8.6704 E 11° 35.3608 48.14450618 11.58934636
    167.46         8 N 48° 8.6165 E 11° 35.2481 48.14360773 11.58746863
     81.72        16 N 48° 8.5895 E 11° 35.1918 48.14315849 11.58652980
     38.86        32 N 48° 8.5760 E 11° 35.1636 48.14293387 11.58606038
     17.44        64 N 48° 8.5693 E 11° 35.1495 48.14282156 11.58582568
      6.72       128 N 48° 8.5659 E 11° 35.1425 48.14276540 11.58570833
      1.37       256 N 48° 8.5642 E 11° 35.1390 48.14273732 11.58564965
     -1.31       512 N 48° 8.5634 E 11° 35.1372 48.14272329 11.58562031
      0.03      1024 N 48° 8.5638 E 11° 35.1381 48.14273031 11.58563498
 
Oben