def clipLog(im, lexp=1000, cmap=[-1.4914, 0.6273], scale=[-0.1, 10]):
    """
    Return normalized array like DS9 log
    """
    import numpy as np
    contrast, bias = cmap
    clip = (np.clip(im, scale[0], scale[1]) - scale[0]) / (scale[1] - scale[0])
    clip_log = np.clip((np.log10(lexp * clip + 1) / np.log10(lexp) - bias) *
        contrast + 0.5, 0, 1)
    return clip_log
