%pylab inline
import numpy as np
import matplotlib.pylab as plt
from nipy.modalities.fmri import design_matrix as dmtx
highpass = dmtx._cosine_drift(2, np.arange(0,180,1))
highpass.shape
plt.imshow(highpass, interpolation="nearest")
plt.plot(highpass[:,80])
plt.plot(highpass[:,0:8])