def bow_histogram(data, codebook, pts=None, shape=None): code, dist = vq(data, codebook) code_hist = bow(data, code, codebook.shape[0]) return code_hist