def hist_stretch(im, percentiles=(1, 99)): p2, p98 = np.percentile(im, percentiles) im = im *100000 #im = np.array(im, np.int64) return exposure.rescale_intensity(im, in_range=percentiles)