Stabilized slope determination.

This commit is contained in:
Ludger Küperkoch 2017-05-31 11:41:24 +02:00
parent 60c882fe04
commit d7abd3bef9

View File

@ -240,7 +240,15 @@ class AICPicker(AutoPicker):
ax.set_yticks([])
ax.set_title(self.Data[0].stats.station)
return
islope = islope[0][0:imax]
if len(islope) < 2:
# get local maximum of AIC function to calculate properly slope
islope = np.where((self.Tcf <= min([self.Pick + tslope, len(self.Data[0].data)])) \
& (self.Tcf >= self.Pick))
imax = np.argmax(aicsmooth[islope])
islope = islope[0][0:imax]
dataslope = self.Data[0].data[islope]
# calculate slope as polynomal fit of order 1
xslope = np.arange(0, len(dataslope), 1)