Gradient domain high dynamic range compression