From d25e791bb8d21a8491b97206cf5f58ea090c7e16 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Mon, 22 Jun 2015 16:01:25 +0200 Subject: [PATCH] Debugging. --- pylot/core/pick/Picker.py | 14 +++++++++++--- 1 file changed, 11 insertions(+), 3 deletions(-) diff --git a/pylot/core/pick/Picker.py b/pylot/core/pick/Picker.py index 3ca67ef5..d0c7a08f 100644 --- a/pylot/core/pick/Picker.py +++ b/pylot/core/pick/Picker.py @@ -292,6 +292,7 @@ class PragPicker(AutoPicking): self.Pick = None self.SNR = None self.slope = None + pickflag = 0 #smooth CF ismooth = int(round(self.Tsmooth / self.dt)) cfsmooth = np.zeros(len(self.cf)) @@ -354,15 +355,22 @@ class PragPicker(AutoPicking): #now decide which pick: left or right? if flagpick_l > 0 and flagpick_r > 0 and cfpick_l <= cfpick_r: self.Pick = pick_l + pickflag = 1 elif flagpick_l > 0 and flagpick_r > 0 and cfpick_l >= cfpick_r: self.Pick = pick_r + pickflag = 1 + else: + print 'PragPicker: Could not find reliable onset!' + self.Pick = None + pickflag = 0 if self.getiplot() > 1: p = plt.figure(self.getiplot()) p1, = plt.plot(Tcfpick,cfipick, 'k') p2, = plt.plot(Tcfpick,cfsmoothipick, 'r') - p3, = plt.plot([self.Pick, self.Pick], [min(cfipick), max(cfipick)], 'b', linewidth=2) - plt.legend([p1, p2, p3], ['CF', 'Smoothed CF', 'Pick']) + if pickflag > 0: + p3, = plt.plot([self.Pick, self.Pick], [min(cfipick), max(cfipick)], 'b', linewidth=2) + plt.legend([p1, p2, p3], ['CF', 'Smoothed CF', 'Pick']) plt.xlabel('Time [s] since %s' % self.Data[0].stats.starttime) plt.yticks([]) plt.title(self.Data[0].stats.station) @@ -371,6 +379,6 @@ class PragPicker(AutoPicking): plt.close(p) else: - self.Pick = None print 'PragPicker: No initial onset time given! Check input!' + self.Pick = None return