In [1]:
%matplotlib inline
import requests
import pandas as pd
In [2]:
import sys
path = '/home/usgs/github/cgsn-parsers/'
sys.path.append(path)
from cgsn_parsers.parsers.parse_mopak import Parser
In [3]:
#Coastal Endurance - WA Shelf Surface Piercing Profiler Mooring
baseurl = "https://rawdata.oceanobservatories.org/files/CE07SHSM/D00003/cg_data/dcl11/mopak/"
fname = "20160505_220005.mopak.log"
In [4]:
# Coastal Pioneer - Central Surface Mooring
baseurl = "https://rawdata.oceanobservatories.org/files/CP01CNSM/D00006/cg_data/dcl11/mopak/"
fname = "20170116_150009.mopak.log"
In [5]:
# initialize the Parser object for METBK
mopak = Parser(baseurl + fname)
In [6]:
r = requests.get(mopak.infile, verify=True) # use verify=False for expired certificate
mopak.raw = r.content
In [7]:
len(mopak.raw), mopak.raw[-5:]
Out[7]:
In [8]:
mopak.parse_data()
In [9]:
mopak.data.keys()
Out[9]:
In [10]:
df = pd.DataFrame(mopak.data)
df['dt_utc'] = pd.to_datetime(df.time, unit='s')
df.set_index('dt_utc', drop=True, inplace=True)
In [11]:
# Later, can drop time, dcl_date_time_string, and maybe dt_utc (not the index)
df.shape, df.columns
Out[11]:
In [12]:
df.dtypes
Out[12]:
In [13]:
df.tail(10)
Out[13]:
In [14]:
df[['acceleration_x', 'acceleration_y']].plot(figsize=(11,5));
In [15]:
df['angular_rate_z'].plot(figsize=(11,5));