The objective of intrinsic images estimation is to decompose an input image into its intrinsic shading and reflectance components. This is a well-known under-constrained problem that has long been an open challenge. This paper proposes a novel approach for automatic intrinsic images decomposition that uses a new reflectance sparsity prior. On the basis of the observation that the reflectance of natural objects is commonly piecewise constant, we formalize this constraint on the entire reflectance image using the \(L_0\) sparse loss function that enforces the variation in reflectance images to be of high-frequency and sparse. This new sparsity constraint significantly improves the quality of Retinex intrinsic images estimation. It also functions effectively by combining a class of global sparsity priors on reflectance. Experimental results on MIT benchmark dataset as well as various real-world images and synthetic images demonstrate the effectiveness and versatility of our approach.