Aktualisierung...

pvl.py

@@ -5,10 +5,12 @@ import sys              #importieren saemtlicher Module
 import os
 import random as rdm
 import re as re         #regular express fuer Bibo regulaerer ausdruecke
-import math as sqrt       #modul fuer mathematische Ausdruecke
+
+from math import sin, cos, sqrt, atan2, radians     #Teilelemente aus math fuer mathematische Ausdruecke
 
 
 #Begruessung
+os.system("clear")                                      #loeschen des Bildschirms
 print ("-")*50
 print ("Luftlinienberechner")
 print ("Erstellt von Niko Hessler und Martin Mittrenga")
@@ -19,80 +21,110 @@ print ("-")*50
 #========================================================================
 
 
-#Datei einlesen mit print als Test
-os.system("clear")
-f = open("/data/share/nhessler/DE.tab", "r")
-for line in f:
-    print line.rstrip()                 #in liste schreiben und regularen ausdruck aus liste lesen!!!!!!!!!
-
-f.close()
+#Datei einlesen
+datenbank = open("/data/share/nhessler/pvl/DE.tab", "r").readlines() #Einlesen der Datei
 
 
 #========================================================================
 
 
+#Startpunkt Auswahl
+startpunkt = (raw_input("Geben sie Ihren Startpunkt an: "))
+print ("-")*50
+
+
+idee = []
+for line in datenbank:
+    if re.search(startpunkt, line):
+
+        idee.append(line.split())          #Schreiben als liste
+
+for a in range(len(idee)):
+	print "%i. idee:%s"%(a+1,idee[a])
+print ("-")*50
+
+choice = int(raw_input("Welche Option aus der Liste ist gewuenscht? "))-1
+print idee[choice]
+latitudestart = idee[choice][4]
+longitudestart = idee[choice][5]
+
+print ("-")*50
 
 
-#Startpunkt Auswahl
-startpunkt=(raw_input("Geben sie Ihren Startpunkt an: "))
 
-if re.search(startpunkt, line):
-        print ("Gefunden, super.")
-else:
-        print ("Nicht gefunden, schade.")
 
 
 
 
 #Endpunkt Auswahl
-endpunkt=(raw_input("Geben sie Ihren Endpunkt an: "))
+endpunkt = (raw_input("Geben sie Ihren Endpunkt an: "))
+print ("-")*50
 
-if re.search(endpunkt, line):
-        print ("Gefunden, super.")
-else:
-        print ("Nicht gefunden, schade.")
 
+idee = []                           #leeren der Liste
+for line in datenbank:
+    if re.search(endpunkt, line):
 
+        idee.append(line.split())          #Schreiben als liste
 
+for a in range(len(idee)):
+	print "%i. idee:%s"%(a+1,idee[a])              #%s string eingefuegt    #%i dezimal/int
+print ("-")*50
 
+choice = int(raw_input("Welche Option aus der Liste ist gewuenscht? "))-1
 
-#Berechnung allgemein
+latitudeende = idee[choice][4]
+longitudeende = idee[choice][5]
 
-lonstart=0    #Wie bekommen wir die Ergebnisse aus der DB?
-lonend=0
-latstart=0
-latend=0
 
 
+#========================================================================
 
 
 
-pos_x_lon = 111,3*cos(lat)*(pos_x_lon-pos_y_lon)
-pos_y_lon = 111,3*(lat1-lat2)
+print "-"*50
 
+print "Zwischenausgabe der vier uebergebenen Variablen zur Berechnung:"
+print latitudestart
+print longitudestart
+print latitudeende
+print longitudeende
 
+print "-"*50
 
-#Berechnung fuer Longitude
-c_square = pos_x_lon**2 + pos_y_lon**2
-c = int(sqrt(c_square))
-if ((c_square - c**2) == 0):
-        print c
-        print "Danke fuer die Eingabe."
-        print "-"*50
 
+#Berechnung
+R = 6373.0 #Radius der Erde in Kilometer
 
+#Startpunkt
+lat1 = radians(float(latitudestart))
+lon1 = radians(float(longitudestart))
 
+#Endpunkt
+lat2 = radians(float(latitudeende))
+lon2 = radians(float(longitudeende))
+
+#Differenz
+dlon = lon2-lon1
+dlat = lat2-lat1
+
+
+a = sin(dlat / 2)**2 + cos(lat1) * cos(lat2) * sin(dlon / 2)**2
+c = 2*atan2(sqrt(a), sqrt(1-a))
+
+distance = R * c
 
 print "-"*50
-print ("Ihre Entfernung betraegt XXX Kilometer. Gute Fahrt!")
+print ("Ihre Entfernung betraegt:"), distance, ("Kilometer. Gute Fahrt!")
 print "-"*50
 
 
+
 #======================================================================
 
 
 
-while True:											#Sicherheitsabfrage ob wirklich beendet werden soll
+while True:		#Sicherheitsabfrage ob wirklich beendet werden soll
 
         eingabe=raw_input("Wollen Sie Wirklich beenden? Bitte 'J' oder 'N' eingeben: ")
         if eingabe=="N":