features = np.zeros([len(imlist),162])
for i,f in enumerate(imlist):
    im = np.array(Image.open(f))
    h=mahotas.features.tas(im)
    features[i] = h.flatten()