Skip to content

Instantly share code, notes, and snippets.

View stefan789's full-sized avatar

stefan789

  • Technical University München
View GitHub Profile
@stefan789
stefan789 / liveplot.py
Created June 23, 2017 12:20
Fluxgate readout KoreaPC
import matplotlib.pyplot as plt
import matplotlib.animation as animation
import time
import sys
fig = plt.figure()
ax1 = fig.add_subplot(1,1,1)
ax1.ticklabel_format(style="plain")
fil = sys.argv[1]
import numpy as np
import json
import cloudant
import datetime as dt
import pytz
berlin = pytz.timezone("Europe/Berlin")
# set up db and define functions
#acct = cloudant.Account(uri="http://10.155.59.88:5984")
@stefan789
stefan789 / allan.py
Last active September 10, 2015 12:30 — forked from mgmarino/allan.py
allan deviaton
import numpy as np
def allanVar(data, sf=1):
'''
calculates allan variance
according to eq 8.13a in
David W. Allan, John H. Shoaf and Donald Halford: Statistics of Time and Frequency Data Analysis,
NBS Monograph 140, pages 151–204, 1974
'''
import json
import cloudant
import datetime as dt
# set up db and define functions
#acct = cloudant.Account(uri="http://10.155.59.88:5984")
acct = cloudant.Account(uri="http://raid.nedm1:5984")
res = acct.login("nedm_user", "")
assert res.status_code == 200
@stefan789
stefan789 / allan.py
Created August 19, 2015 21:26
allan deviaton
import numpy as np
def allanVar(data, sf=1):
'''
calculates allan variance
according to eq 8.13a in
David W. Allan, John H. Shoaf and Donald Halford: Statistics of Time and Frequency Data Analysis,
NBS Monograph 140, pages 151–204, 1974
'''
leng = len(data)
@stefan789
stefan789 / plt_dates.py
Last active September 11, 2015 09:02
Plot db data with timezones
fig, ax = plt.subplots()
fig.set_size_inches(8,6)
fmt = mdates.DateFormatter("%H:%M", tz=berlin)
ax.xaxis.set_major_formatter(fmt)
ax.plot_date([i.astimezone(berlin) for i in data["Date"]], data["Bz1"], "-", label="y1")
fig.autofmt_xdate()
#ax.set_xlim([dt.datetime(2015,8,7,17,0,0), dt.datetime(2015,8,7,19,0,0)])
#ax.set_ylim([-30,40])
plt.legend(loc='center right',
import numpy as np
import json
import cloudant
import datetime as dt
import pytz
berlin = pytz.timezone("Europe/Berlin")
# set up db and define functions
#acct = cloudant.Account(uri="http://10.155.59.88:5984")
import ROOT
import numpy as np
import re
import argparse
ROOT.gSystem.Load("/home/stefan/Programs/TWaveform/lib/libWaveWaveBase.so")
p = argparse.ArgumentParser()
p.add_argument('filena')
p.add_argument('samplerate', type=float, help="sample rate")
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1 import AxesGrid
import argparse
parser = argparse.ArgumentParser()
parser.add_argument("comp", help = "b field component to be plotted, Bx, By, Bz or normB")
args = parser.parse_args()
import cloudant
import json
import codecs
# Authentication
acct = cloudant.Account(uri="http://raid.nedm1:5984")
res = acct.login("username", "password")
assert res.status_code == 200
# Grab the correct database