def normalize(x): y = np.copy(x); minv = -1000; y[x<minv]=minv; maxv = 200; y[x>maxv]=maxv return ((y*1.0-minv)/(maxv-minv)*255).astype(np.uint8);#0-255, to save disk space