-
Notifications
You must be signed in to change notification settings - Fork 16
/
from_ooi_json-pwrsys-Copy1.py
108 lines (67 loc) · 2.04 KB
/
from_ooi_json-pwrsys-Copy1.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
# coding: utf-8
# # Convert OOI Parsed pwrsys JSON to NetCDF file
# using CF-1.6, Discrete Sampling Geometry (DSG) conventions, **`featureType=timeSeries`**
# In[1]:
get_ipython().magic('matplotlib inline')
import json
import pandas as pd
import numpy as np
import glob
from pyaxiom.netcdf.sensors import TimeSeries
# In[2]:
def json2df(infile):
with open(infile) as jf:
df = pd.DataFrame(json.load(jf))
df['time'] = pd.to_datetime(df.time, unit='s')
return df
# In[3]:
path = '/sand/usgs/users/rsignell/data/ooi/endurance/cg_proc/ce02shsm/D00004/buoy/pwrsys/*.pwrsys.json'
odir = '/usgs/data2/notebook/data/nc'
ofile = 'ce02shsm_pwrsys_D00004.nc'
# In[4]:
df = pd.concat([json2df(file) for file in glob.glob(path)])
# In[5]:
df['depth'] = 0.0
df.index = df['time']
# In[6]:
df['solar_panel4_voltage'].plot();
# ### Define the NetCDF global attributes
# In[7]:
global_attributes = {
'institution':'Oregon State University',
'title':'OOI CE02SHSM Pwrsys Data',
'summary':'OOI Pwrsys data from Coastal Endurance Oregon Shelf Surface Mooring',
'creator_name':'Chris Wingard',
'creator_email':'cwingard@coas.oregonstate.edu',
'creator_url':'http://ceoas.oregonstate.edu/ooi'
}
# ### Create initial file
# In[9]:
ts = TimeSeries(
output_directory=odir,
latitude=44.64,
longitude=-124.31,
station_name='ce02shsm',
global_attributes=global_attributes,
times=df.time.values.astype(np.int64) // 10**9,
verticals=df.depth.values,
output_filename=ofile,
vertical_positive='down'
)
# ### Add data variables
# In[10]:
df.columns.tolist()
# In[11]:
for c in df.columns:
if c in ts._nc.variables:
print("Skipping '{}' (already in file)".format(c))
continue
if c in ['time', 'lat', 'lon', 'depth', 'cpm_date_time_string']:
print("Skipping axis '{}' (already in file)".format(c))
continue
print("Adding {}".format(c))
try:
ts.add_variable(c, df[c].values)
except:
print('skipping, hit object')
# In[ ]: