def evaluate_features(output_file='features.pdf'):
    import algorithms
    blue = (0.0, 0.0, 1.0, 1.0)
    g = face(1, True, 18, 12)
    h = face(2, True, 18, 12)
    f = face(3, True, 18, 12)
    k = face(0, True, 18, 12)

    figure = pyplot.figure()

    subplot = pyplot.subplot(1, 5, 1)

    subplot.imshow(g.Z)

    for xx in range(1, 8):
        v = int((xx) * (g.height / 8))
        subplot.axhline(v, color=blue)

    subplot.set_xlabel('Different regions (N=8)')
    subplot.xaxis.set_visible(False)
    subplot.yaxis.set_visible(False)

    subplot = pyplot.subplot(1, 5, 2)
    subplot.imshow(g.features[1].reshape(g.features[2]))
    subplot.set_xlabel('Person 1a')
    subplot.xaxis.set_visible(False)
    subplot.yaxis.set_visible(False)

    subplot = pyplot.subplot(1, 5, 3)
    subplot.imshow(h.features[1].reshape(h.features[2]))
    subplot.set_xlabel('Person 1b')
    subplot.xaxis.set_visible(False)
    subplot.yaxis.set_visible(False)

    subplot = pyplot.subplot(1, 5, 4)
    subplot.imshow(f.features[1].reshape(f.features[2]))
    subplot.set_xlabel('Person 1c')
    subplot.xaxis.set_visible(False)
    subplot.yaxis.set_visible(False)

    subplot = pyplot.subplot(1, 5, 5)
    subplot.imshow(k.features[1].reshape(k.features[2]))
    subplot.set_xlabel('Person 2a') 
    subplot.xaxis.set_visible(False)
    subplot.yaxis.set_visible(False)   

    figure.savefig(output_file, format='pdf', dpi=600, orientation='landscape', bbox_inches="tight")
dataset/ETHPy150Open basilfx/Course-Face3D/face3d/utils.py/evaluate_features/original