15

I am trying to vary the color of a line plotted from data in two arrays (eg. ax.plot(x,y)). The color should vary as the index into x and yincreases. I am essentially trying to capture the natural 'time' parameterization of the data in arrays x and y.

In a perfect world, I want something like:

fig = pyplot.figure()
ax  = fig.add_subplot(111)
x   = myXdata 
y   = myYdata

# length of x and y is 100
ax.plot(x,y,color=[i/100,0,0]) # where i is the index into x (and y)

to produce a line with color varying from black to dark red and on into bright red.

I have seen examples that work well for plotting a function explicitly parameterized by some 'time' array, but I can't get it to work with raw data...

Tom Stephens
  • 475
  • 4
  • 15

1 Answers1

22

The second example is the one you want... I've edited it to fit your example, but more importantly read my comments to understand what is going on:

import numpy as np
from matplotlib import pyplot as plt
from matplotlib.collections import LineCollection

x   = myXdata 
y   = myYdata
t = np.linspace(0,1,x.shape[0]) # your "time" variable

# set up a list of (x,y) points
points = np.array([x,y]).transpose().reshape(-1,1,2)
print points.shape  # Out: (len(x),1,2)

# set up a list of segments
segs = np.concatenate([points[:-1],points[1:]],axis=1)
print segs.shape  # Out: ( len(x)-1, 2, 2 )
                  # see what we've done here -- we've mapped our (x,y)
                  # points to an array of segment start/end coordinates.
                  # segs[i,0,:] == segs[i-1,1,:]

# make the collection of segments
lc = LineCollection(segs, cmap=plt.get_cmap('jet'))
lc.set_array(t) # color the segments by our parameter

# plot the collection
plt.gca().add_collection(lc) # add the collection to the plot
plt.xlim(x.min(), x.max()) # line collections don't auto-scale the plot
plt.ylim(y.min(), y.max())
Ethan Coon
  • 751
  • 5
  • 16
  • Thank you for pointing out what is happening with the reshape and concatenation. This is working well. – Tom Stephens Apr 21 '12 at 02:01
  • If you want a smoother transition between line segments you can do `segs = np.concatenate([points[:-2], points[1:-1], points[2:]], axis=1)` instead. – shockburner Oct 20 '17 at 17:58