From f1e3720a0e77c86b52836b5321ad29dcd0bb69f4 Mon Sep 17 00:00:00 2001
From: Marcel Paffrath <marcel.paffrath@rub.de>
Date: Mon, 10 Apr 2017 17:06:39 +0200
Subject: [PATCH] gathering some ideas for map projection

---
 map_projection.py     | 56 +++++++++++++++++++++++++++++++++++++++++++
 pylot/RELEASE-VERSION |  2 +-
 2 files changed, 57 insertions(+), 1 deletion(-)
 create mode 100644 map_projection.py

diff --git a/map_projection.py b/map_projection.py
new file mode 100644
index 00000000..cef4e4b1
--- /dev/null
+++ b/map_projection.py
@@ -0,0 +1,56 @@
+from mpl_toolkits.basemap import Basemap
+import matplotlib.pyplot as plt
+import numpy as np
+import obspy
+
+#import QtPyLoT
+from pylot.core.util.dataprocessing import read_metadata
+
+#pf=QtPyLoT.main()
+
+def onpick(event):
+    ind = event.ind
+    print(ind)
+
+
+metadata=read_metadata('/data/Geothermie/Insheim/STAT_INFO/MAGS2_net.dless')
+parser=metadata[1]
+
+lat=[]
+lon=[]
+stcl=[]
+picks=[]
+picks_rel=[]
+
+for station in parser.stations:
+    stcl.append(station[0].station_call_letters)
+    lat.append(station[0].latitude)
+    lon.append(station[0].longitude)
+
+for station in stcl:
+    try:
+        picks.append(pf.autopicks[station]['P']['mpp'])
+    except:
+        picks.append(np.nan)
+
+minp = min(picks)
+
+for pick in picks:
+    if type(pick) is obspy.core.utcdatetime.UTCDateTime:
+        pick -= minp
+    picks_rel.append(pick)
+
+
+m = Basemap(projection='hammer',lon_0=0,lat_0=0)
+m.drawmapboundary(fill_color='aqua')
+m.drawcountries()
+m.drawstates()
+m.fillcontinents(color='grey', lake_color='aqua')
+m.drawcoastlines()
+m.scatter(lon, lat, s=50, c=picks_rel, latlon=True, zorder=10, picker=True)
+
+ax = plt.gca() # IMPROVE!!!!
+
+ax.figure.canvas.mpl_connect('pick_event', onpick)
+
+plt.show()
diff --git a/pylot/RELEASE-VERSION b/pylot/RELEASE-VERSION
index 3a0d6f3c..9d801014 100644
--- a/pylot/RELEASE-VERSION
+++ b/pylot/RELEASE-VERSION
@@ -1 +1 @@
-0656-dirty
+d370-dirty