Abstract We consider estimating a piecewise-constant image, or gradient-sparse signal on general graph, from noisy linear measurements. propose and study an iterative algorithm to minimize penalized least-squares objective, with penalty given by the “ℓ0-norm” of signal’s discrete graph gradient. The method uses non-convex variant proximal gradient descent, applying alpha-expansion procedure app...