for k in range(class_count):
                norm = scipy.stats.norm(class_mean[k], class_std[k])
                pixel_class_prob[:, k] = class_prob[k] * norm.pdf(data)