import numpy as np
import zlib
import pylab as plt

yt = np.random.rand(2000)
N = 1000
ids = np.stack([np.isclose(yt,np.linspace(yt.min(),yt.max(),N)[i], rtol=0.01) for i in range(N)]).argmax(0)
dec = np.linspace(yt.min(),yt.max(),N)[ids]
enc = zlib.compress(ids)
enc_org = zlib.compress(yt)
ratio = len(enc)/len(enc_org)
mse = np.mean((dec-yt)**2)
print("ratio :", ratio, "MSE :", mse)