marcel/pylot
3
0
Fork 0
Code Issues 7 Pull Requests 1 Releases Wiki Activity

Stabilized AICPicker by introducing 1st derivative of CF to find global maximum.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-03-10 16:18:32 +01:00
parent 380cccdf16
commit ea68b38f7e
1 changed files with 16 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
pylot/core/pick/Picker.py
View File

@ -52,7 +52,7 @@ class AutoPicking(object):
:type: float :type: float
:param: Pick1, initial (prelimenary) onset time, starting point for PragPicker and :param: Pick1, initial (prelimenary) onset time, starting point for PragPicker and
EarlLatePick EarlLatePicker
:type: float :type: float
''' '''
@ -189,8 +189,16 @@ class AICPicker(AutoPicking):
#remove offset #remove offset
offset = abs(min(aic) - min(aicsmooth)) offset = abs(min(aic) - min(aicsmooth))
aicsmooth = aicsmooth - offset aicsmooth = aicsmooth - offset
#get maximum of CF as starting point #get maximum of 1st derivative of CF (more stable!) as starting point
icfmax = np.argmax(aic) diffcf = np.diff(aicsmooth)
#find NaN's
nn = np.isnan(diffcf)
if len(nn) > 1:
diffcf[nn] = 0
#taper CF to get rid off side maxima
tap = np.hanning(len(diffcf))
diffcf = tap * diffcf * max(abs(aicsmooth))
icfmax = np.argmax(diffcf)
#find minimum in front of maximum #find minimum in front of maximum
lpickwindow = int(round(self.PickWindow / self.dt)) lpickwindow = int(round(self.PickWindow / self.dt))
@ -222,6 +230,11 @@ class AICPicker(AutoPicking):
#find maximum within slope determination window #find maximum within slope determination window
#'cause slope should be calculated up to first local minimum only! #'cause slope should be calculated up to first local minimum only!
imax = np.argmax(self.Data[0].data[islope]) imax = np.argmax(self.Data[0].data[islope])
if imax == 0:
print 'AICPicker: Maximum for slope determination right at the beginning of the window!'
print 'Choose longer slope determination window!'
pdb.set_trace()
return
islope = islope[0][0 :imax] islope = islope[0][0 :imax]
dataslope = self.Data[0].data[islope] dataslope = self.Data[0].data[islope]
#calculate slope as polynomal fit of order 1 #calculate slope as polynomal fit of order 1