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

[change] docstrings in compare.py

Browse Source
This commit is contained in:
Darius Arnold 2017-10-16 21:02:04 +02:00
parent 4d37abc7d8
commit fab1f53a41
1 changed files with 33 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -109,10 +109,9 @@ class Comparison(object):
Comparison is carried out with the help of pdf representation of the picks Comparison is carried out with the help of pdf representation of the picks
and a probabilistic approach to the time difference of two onset and a probabilistic approach to the time difference of two onset
measurements. measurements.
:param a: filename for pickset A :param type: type of the returned `~pylot.core.util.pdf.ProbabilityDensityFunction` object.
:type a: str Possible values: 'exp' and 'gauss', representing the type of branches of the PDF
:param b: filename for pickset B :type type: str
:type b: str
:return: dictionary containing the resulting comparison pdfs for all picks :return: dictionary containing the resulting comparison pdfs for all picks
:rtype: dict :rtype: dict
""" """
@ -142,8 +141,7 @@ class Comparison(object):
istations = range(nstations) istations = range(nstations)
fig, axarr = plt.subplots(nstations, 2, sharex='col', sharey='row') fig, axarr = plt.subplots(nstations, 2, sharex='col', sharey='row')
for n in istations: for n, station in enumerate(stations):
station = stations[n]
if station not in self.comparison.keys(): if station not in self.comparison.keys():
continue continue
compare_pdf = self.comparison[station] compare_pdf = self.comparison[station]
@ -190,6 +188,20 @@ class Comparison(object):
return self.get_array(phase, 'standard_deviation') return self.get_array(phase, 'standard_deviation')
def hist_expectation(self, phases='all', bins=20, normed=False): def hist_expectation(self, phases='all', bins=20, normed=False):
"""
Plot a histogram of the expectation values of the PDFs.
Expectation represents the time difference between two most likely arrival times
:param phases: type of phases to compare
:type phases: str
:param bins: number of bins in histogram
:type bins: int
:param normed: Normalize histogram
:type normed: bool
:return: None
:rtype: None
"""
phases.strip() phases.strip()
if phases.find('all') is 0: if phases.find('all') is 0:
phases = 'ps' phases = 'ps'
@ -210,6 +222,20 @@ class Comparison(object):
plt.show() plt.show()
def hist_standard_deviation(self, phases='all', bins=20, normed=False): def hist_standard_deviation(self, phases='all', bins=20, normed=False):
"""
Plot a histogram of the compared standard deviation values of two arrivals.
Standard deviation of two compared picks represents the combined uncertainties/pick errors
(earliest possible pick, latest possible pick)
:param phases: type of phases to compare
:type phases: str
:param bins: number of bins in histogram
:type bins: int
:param normed: Normalize histogram
:type normed: bool
:return: None
:rtype: None
"""
phases.strip() phases.strip()
if phases.find('all') == 0: if phases.find('all') == 0:
phases = 'ps' phases = 'ps'
@ -370,7 +396,7 @@ class PDFDictionary(object):
class PDFstatistics(object): class PDFstatistics(object):
""" """
This object can be used to get various statistic values from probabillity density functions. This object can be used to get various statistic values from probability density functions.
Takes a path as argument. Takes a path as argument.
""" """