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

Debuged and cleaned source code

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-02-25 09:59:59 +01:00
parent 1966a2b612
commit 4a48874f88
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
pylot/core/pick/CharFuns.py
View File

@ -17,8 +17,6 @@ autoregressive prediction: application ot local and regional distances, Geophys.
""" """
import numpy as np import numpy as np
from obspy.core import Stream from obspy.core import Stream
import pdb
import matplotlib.pyplot as plt
class CharacteristicFunction(object): class CharacteristicFunction(object):
''' '''
@ -235,7 +233,7 @@ class AICcf(CharacteristicFunction):
cf[ff] = 0 cf[ff] = 0
self.cf = cf - np.mean(cf) self.cf = cf - np.mean(cf)
self.xcf = xnp self.xcf = x
class HOScf(CharacteristicFunction): class HOScf(CharacteristicFunction):
@ -287,7 +285,7 @@ class HOScf(CharacteristicFunction):
if len(nn) > 1: if len(nn) > 1:
LTA[nn] = 0 LTA[nn] = 0
self.cf = LTA self.cf = LTA
self.xcf = xnp self.xcf = x
class ARZcf(CharacteristicFunction): class ARZcf(CharacteristicFunction):
@ -326,7 +324,7 @@ class ARZcf(CharacteristicFunction):
if len(nn) > 1: if len(nn) > 1:
cf[nn] = 0 cf[nn] = 0
self.cf = cf self.cf = cf
self.xcf = xnp self.xcf = x
def arDetZ(self, data, order, rind, ldet): def arDetZ(self, data, order, rind, ldet):
''' '''
@ -431,7 +429,7 @@ class ARHcf(CharacteristicFunction):
ldet = int(round(self.getTime1() / self.getIncrement())) #length of AR-determination window [samples] ldet = int(round(self.getTime1() / self.getIncrement())) #length of AR-determination window [samples]
lpred = int(np.ceil(self.getTime2() / self.getIncrement())) #length of AR-prediction window [samples] lpred = int(np.ceil(self.getTime2() / self.getIncrement())) #length of AR-prediction window [samples]
cf = np.zeros(tend - lpred + 1) cf = np.zeros(len(xenoise))
loopstep = self.getARdetStep() loopstep = self.getARdetStep()
arcalci = ldet + self.getOrder() - 1 #AR-calculation index arcalci = ldet + self.getOrder() - 1 #AR-calculation index
for i in range(ldet + self.getOrder() - 1, tend - 2 * lpred + 1): for i in range(ldet + self.getOrder() - 1, tend - 2 * lpred + 1):
@ -563,7 +561,7 @@ class AR3Ccf(CharacteristicFunction):
ldet = int(round(self.getTime1() / self.getIncrement())) #length of AR-determination window [samples] ldet = int(round(self.getTime1() / self.getIncrement())) #length of AR-determination window [samples]
lpred = int(np.ceil(self.getTime2() / self.getIncrement())) #length of AR-prediction window [samples] lpred = int(np.ceil(self.getTime2() / self.getIncrement())) #length of AR-prediction window [samples]
cf = np.zeros(tend - lpred + 1) cf = np.zeros(len(xenoise))
loopstep = self.getARdetStep() loopstep = self.getARdetStep()
arcalci = ldet + self.getOrder() - 1 #AR-calculation index arcalci = ldet + self.getOrder() - 1 #AR-calculation index
for i in range(ldet + self.getOrder() - 1, tend - 2 * lpred + 1): for i in range(ldet + self.getOrder() - 1, tend - 2 * lpred + 1):