Matplotlib中的平行坐标图

使用传统的绘图类型，可以相对直接地查看二维和三维数据。即使具有四维数据，我们也经常可以找到一种显示数据的方法。但是，尺寸大于4的尺寸越来越难以显示。幸运的是，平行坐标图提供了一种查看较大尺寸结果的机制。

一些绘图包提供了平行坐标图，例如Matlab，R，VTK类型1和VTK类型2，但是我看不到如何使用Matplotlib创建一个。

1. Matplotlib中是否有内置的平行坐标图？我当然在画廊里没看到一个。
2. 如果没有内置类型，是否可以使用Matplotlib的标准功能来构建平行坐标图？

编辑：

基于下面的Zhenya提供的答案，我开发了以下支持任意数量轴的概括。按照我在上面原始问题中发布的示例的绘图样式，每个轴都有自己的比例尺。我通过标准化每个轴点处的数据并使轴的范围为0到1来实现此目的。然后，我向每个刻度线应用标签，以在该截距处给出正确的值。

该函数通过接受可迭代的数据集来工作。每个数据集被认为是一组点，其中每个点位于不同的轴上。中的示例以__main__两组30行的形式为每个轴获取随机数。在导致线聚集的范围内，线是随机的；我想验证的行为。

此解决方案不如内置解决方案好，因为您的鼠标行为不正常，并且我伪造标签上的数据范围，但是直到Matplotlib添加内置解决方案之前，它都是可以接受的。

#!/usr/bin/python
import matplotlib.pyplot as plt
import matplotlib.ticker as ticker

def parallel_coordinates(data_sets, style=None):

    dims = len(data_sets[0])
    x    = range(dims)
    fig, axes = plt.subplots(1, dims-1, sharey=False)

    if style is None:
        style = ['r-']*len(data_sets)

    # Calculate the limits on the data
    min_max_range = list()
    for m in zip(*data_sets):
        mn = min(m)
        mx = max(m)
        if mn == mx:
            mn -= 0.5
            mx = mn + 1.
        r  = float(mx - mn)
        min_max_range.append((mn, mx, r))

    # Normalize the data sets
    norm_data_sets = list()
    for ds in data_sets:
        nds = [(value - min_max_range[dimension][0]) / 
                min_max_range[dimension][2] 
                for dimension,value in enumerate(ds)]
        norm_data_sets.append(nds)
    data_sets = norm_data_sets

    # Plot the datasets on all the subplots
    for i, ax in enumerate(axes):
        for dsi, d in enumerate(data_sets):
            ax.plot(x, d, style[dsi])
        ax.set_xlim([x[i], x[i+1]])

    # Set the x axis ticks 
    for dimension, (axx,xx) in enumerate(zip(axes, x[:-1])):
        axx.xaxis.set_major_locator(ticker.FixedLocator([xx]))
        ticks = len(axx.get_yticklabels())
        labels = list()
        step = min_max_range[dimension][2] / (ticks - 1)
        mn   = min_max_range[dimension][0]
        for i in xrange(ticks):
            v = mn + i*step
            labels.append('%4.2f' % v)
        axx.set_yticklabels(labels)


    # Move the final axis' ticks to the right-hand side
    axx = plt.twinx(axes[-1])
    dimension += 1
    axx.xaxis.set_major_locator(ticker.FixedLocator([x[-2], x[-1]]))
    ticks = len(axx.get_yticklabels())
    step = min_max_range[dimension][2] / (ticks - 1)
    mn   = min_max_range[dimension][0]
    labels = ['%4.2f' % (mn + i*step) for i in xrange(ticks)]
    axx.set_yticklabels(labels)

    # Stack the subplots 
    plt.subplots_adjust(wspace=0)

    return plt


if __name__ == '__main__':
    import random
    base  = [0,   0,  5,   5,  0]
    scale = [1.5, 2., 1.0, 2., 2.]
    data = [[base[x] + random.uniform(0., 1.)*scale[x]
            for x in xrange(5)] for y in xrange(30)]
    colors = ['r'] * 30

    base  = [3,   6,  0,   1,  3]
    scale = [1.5, 2., 2.5, 2., 2.]
    data.extend([[base[x] + random.uniform(0., 1.)*scale[x]
                 for x in xrange(5)] for y in xrange(30)])
    colors.extend(['b'] * 30)

    parallel_coordinates(data, style=colors).show()

编辑2：

这是上面的代码在绘制Fisher的Iris数据时产生的示例。它不如Wikipedia的参考图像那么好，但是如果您只有Matplotlib并且需要多维图，则可以通过。