tobeycarman

#!/usr/bin/env python

# Tobey Carman 2014
# Experiments with Tempo DB Python API
# to extract data for Cultivar

import datetime
import os
import re
import json

tobeycarman

# from scipy.signal.detrend docs

from scipy import signal
randgen = np.random.RandomState(9)
npoints = 1e3
noise = randgen.randn(npoints)
trended_noise = 3 + 2*np.linspace(0, 1, npoints) + noise
print "detrended signal minus noise, max difference from original signal: ",
print (signal.detrend(trended_noise) - noise).max()

tobeycarman

import matplotlib.pyplot as plt
import numpy as np

PAR = [
  99.85782790688292,
  0.00029915877415765486,
  0.00010509505480959243,
  0.12320735471319781,
  1.8037705277837454e-05,
  1.634040105529734e-05,