The task of separating an image into distinct components that represent different features plays an important role in many applications. Traditionally, such separation techniques are applied once the image in question has been reconstructed from measured data. We propose an efficient iterative algorithm, where reconstruction is performed jointly with the task of separation. A key assumption is that the image components have different sparse representations. The algorithm is based on a scheme that minimizes a functional composed of a data discrepancy term and the l(1)-norm of the coefficients of the different components with respect to their corresponding dictionaries. The performance is demonstrated for joint 2D deconvolution and separation into curve- and point-like components, and tests are performed on synthetic data as well as experimental stimulated emission depletion and confocal microscopy data. Experiments show that such a joint approach outperforms a sequential approach, where one first deconvolves data and then applies image separation.