An example of how to plot some data on the UW sermeq machine

sermeq-plot-ronne-rifts.py

import pickle
import pandas as pd
import numpy as np
import os

from matplotlib import cm
import matplotlib.pyplot as plt

shelf_name = 'ronne'
datapath = '/data/fast1/arc/'

atl06_file_name = os.path.join(datapath,'atl06/', shelf_name + '.pkl')

with open(atl06_file_name, 'rb') as handle:
    atl06_data = pickle.load(handle)
atl06_data = pd.DataFrame(atl06_data)



subset = pd.DataFrame(atl06_data[atl06_data.rgt==1009])
color=iter(cm.viridis(np.linspace(0,1,len(subset))))

fig,ax = plt.subplots(figsize=(10,10))
fig.patch.set_facecolor('white')
for i,row in subset.iterrows():
    
    c=next(color)
    plt.scatter(row.lat[row.quality==0],row.h[row.quality==0],
               color=c)
#     plt.plot(row.lat,row.h,'o')
    plt.ylim([20,60])
    plt.xlim([-75.5,-75.375])

The output of this code looks like this: