%pylab inline

prefix = 'http://service.iris.edu/irisws/timeseries/1/query?'

SCNL_parameters = 'net=IU&sta=ANMO&loc=00&cha=BHZ&'

times = 'starttime=2005-01-01T00:00:00&endtime=2005-01-02T00:00:00&'

output = 'output=ascii'

import urllib2
f = urllib2.urlopen(prefix + SCNL_parameters + times + output)
timeseries = f.read()

type(timeseries)

data = timeseries.split('\n')
len(data)

data[0]

data[1]

data[-1]

data[-2]

dates = []
values = []

for line in data[1:-1]:
    date, val = line.split()
    dates.append(date)
    values.append(float(val))

plot(values)

len(values)/20.0

_/(60*60.0)

len(values)/44100.0

plot(dates, values)

dates[0]

dates = array(dates, dtype='datetime64')

dates[0]

plot(dates, values)

plot(dates, values)
positions = linspace(dates[0].astype(int64), dates[-1].astype(int64), 24)
xticks(positions, arange(24));

from scipy.io import wavfile

wavfile.write('earthquake.wav', 44100, array(values))

!aplay earthquake.wav

array(values).dtype

audio_data = (array(values)).astype(int16)

wavfile.write('earthquake.wav', 44100, audio_data)

!aplay earthquake.wav

audio_data.shape

audio_data.size/44100.0

from scipy.signal import resample

audio_data2 = resample(audio_data, 5* 44100)

wavfile.write('earthquake2.wav', 44100, audio_data2.astype(int16))

!aplay earthquake2.wav

plot(audio_data2)

audio_data3 = resample(audio_data[400000:1000000], 30* 44100)

wavfile.write('earthquake3.wav', 44100, audio_data3.astype(int16))

!aplay earthquake3.wav

plot(audio_data3)

prefix = 'http://service.iris.edu/irisws/timeseries/1/query?'

SCNL_parameters = 'net=IU&sta=ANMO&loc=00&cha=BHZ&'

times = 'starttime=2014-03-28T00:00:00&endtime=2014-03-29T00:00:00&'

output = 'output=ascii'

import urllib2
f = urllib2.urlopen(prefix + SCNL_parameters + times + output)
timeseries2 = f.read()

data2 = timeseries2.split('\n')



dates2 = []
values2 = []
for line in data2[1:-1]:
    date, val = line.split()
    dates2.append(date)
    values2.append(float(val))
    
plot(values2)

audio_data = (array(values2)+ 50000).astype(int16)
wavfile.write('earthquake4.wav', 44100, audio_data.astype(int16))

!aplay earthquake4.wav

from keys import oer_key

url = 'https://openexchangerates.org/api/latest.json?app_id=%s&base=%s'%(oer_key, 'USD')
f = urllib2.urlopen(url)
timeseries = f.read()

len(timeseries.split('\n'))

import json

s = json.loads(timeseries)

s

values = []
for rate in s['rates'].values():
    values.append(rate)

plot(values)

loop = values*100
wavfile.write('currencies.wav', 44100, array(loop).astype(int16))

!aplay currencies.wav