I have a collection of Nx3 matrices in scipy/numpy and I\'d like to make a 3 dimensional scatter of it, where the X and Y axes are determined by the values of first and second c
AFAIK, legends for 3d scatter artists are not directly supported. See here: http://matplotlib.sourceforge.net/users/legend_guide.html#plotting-guide-legend
However, you can use a hack/work-around with a 'proxy artist', using something like:
p = Rectangle((0, 0), 1, 1, fc="r")
axis.legend([p], ["Red Rectangle"])
So the proxy artist won't be added to the axis, but you can use it to create a legend.