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

[closes 195] in principle comparing automatic and manual picks works

Browse Source
This commit is contained in:
Sebastian Wehling-Benatelli 2016-04-12 04:55:04 +02:00
parent 9359338e81
commit 4a9c02e1d1
2 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
pylot/core/pick/compare.py
View File

@ -134,7 +134,12 @@ class Comparison(object):
pdf_dict = self.comparison pdf_dict = self.comparison
exp_array = list() exp_array = list()
for station, phases in pdf_dict.items(): for station, phases in pdf_dict.items():
try:
exp_array.append(phases[phase].expectation()) exp_array.append(phases[phase].expectation())
except KeyError as e:
print('{err_msg}; station = {station}, phase = {phase}'.format(
err_msg=str(e), station=station, phase=phase))
continue
return exp_array return exp_array
def get_std_array(self, phase): def get_std_array(self, phase):
@ -242,7 +247,8 @@ class PDFDictionary(object):
return pdf_picks return pdf_picks
comp_obj = Comparison(manual='/home/sebastianp/Data/Insheim/e0057.026.13.xml', auto='/home/sebastianp/Data/Insheim/e0057.026.13.xml') #comp_obj = Comparison(manual='/home/sebastianp/Data/Insheim/e0019.048.13.xml',
comp_obj.plot() # auto='/data/Geothermie/Insheim/EVENT_DATA/LOCAL/RefDB/e0019.048.13/autoPyLoT.xml')
comp_obj.hist_expectation() #comp_obj.plot()
comp_obj.hist_standard_deviation() #comp_obj.hist_expectation()
#comp_obj.hist_standard_deviation()

2
pylot/core/util/pdf.py
View File

@ -181,7 +181,7 @@ class ProbabilityDensityFunction(object):
self._x = np.array(x) self._x = np.array(x)
@classmethod @classmethod
def fromPick(self, lbound, barycentre, rbound, incr=0.0001, decfact=0.01, def fromPick(self, lbound, barycentre, rbound, incr=0.001, decfact=0.01,
type='gauss'): type='gauss'):
''' '''
Initialize a new ProbabilityDensityFunction object. Initialize a new ProbabilityDensityFunction object.