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

improved SNR plots

Browse Source
This commit is contained in:
Marcel Paffrath 2016-01-04 15:47:46 +01:00
parent 3c1be950b9
commit acd9e942f9
1 changed files with 11 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
pylot/core/active/surveyUtils.py
View File

@ -16,11 +16,13 @@ def setArtificialPick(shot_dict, traceID, pick):
def fitSNR4dist(shot_dict, shiftdist = 30, shiftSNR = 100): def fitSNR4dist(shot_dict, shiftdist = 30, shiftSNR = 100):
import numpy as np import numpy as np
import matplotlib.pyplot as plt
dists = [] dists = []
picks = [] picks = []
snrs = [] snrs = []
snr_sqrt_inv = [] snr_sqrt_inv = []
snrthresholds = [] snrthresholds = []
snrBestFit = []
for shot in shot_dict.values(): for shot in shot_dict.values():
for traceID in shot.getTraceIDlist(): for traceID in shot.getTraceIDlist():
if shot.getSNR(traceID)[0] >= 1: if shot.getSNR(traceID)[0] >= 1:
@ -31,18 +33,23 @@ def fitSNR4dist(shot_dict, shiftdist = 30, shiftSNR = 100):
fit = np.polyfit(dists, snr_sqrt_inv, 1) fit = np.polyfit(dists, snr_sqrt_inv, 1)
fit_fn = np.poly1d(fit) fit_fn = np.poly1d(fit)
for dist in dists: for dist in dists:
snrBestFit.append((1/(fit_fn(dist)**2)))
dist += shiftdist dist += shiftdist
snrthresholds.append((1/(fit_fn(dist)**2)) - shiftSNR * np.exp(-0.05 * dist)) snrthresholds.append((1/(fit_fn(dist)**2)) - shiftSNR * np.exp(-0.05 * dist))
plotFittedSNR(dists, snrthresholds, snrs) plotFittedSNR(dists, snrthresholds, snrs, snrBestFit)
return fit_fn #### ZU VERBESSERN, sollte fertige funktion wiedergeben return fit_fn #### ZU VERBESSERN, sollte fertige funktion wiedergeben
def plotFittedSNR(dists, snrthresholds, snrs): def plotFittedSNR(dists, snrthresholds, snrs, snrBestFit):
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
plt.interactive(True) plt.interactive(True)
fig = plt.figure() fig = plt.figure()
plt.plot(dists, snrs, '.', markersize = 1.0, label = 'SNR values') plt.plot(dists, snrs, 'b.', markersize = 2.0, label = 'SNR values')
plt.plot(dists, snrthresholds, 'r.', markersize = 1, label = 'Fitted threshold') dists.sort()
snrthresholds.sort(reverse = True)
snrBestFit.sort(reverse = True)
plt.plot(dists, snrthresholds, 'r', markersize = 1, label = 'Fitted threshold')
plt.plot(dists, snrBestFit, 'k', markersize = 1, label = 'Best fitted curve')
plt.xlabel('Distance[m]') plt.xlabel('Distance[m]')
plt.ylabel('SNR') plt.ylabel('SNR')
plt.legend() plt.legend()