Source reconstruction of SLACS lenses, assuming analytical mass and lens light models from Shajib et al. (2020b). Left to right: imaging data with lens-light subtracted, normalised residuals with shapelet model, shapelet source model, normalised residuals with starlet model, starlet source model (r_pix = 3). The maximal polynomial order n_max used for shapelet models are, from top to bottom: 6, 15, 12, 15. We refer to the text for more details, specifically concerning the case of SDSS J0959−0410. Dark grey areas in model panels correspond to non-positive pixel values (which can be negative with shapelets). All source reconstructions are shown at the supersampled resolution corresponding to the chosen r_pix. Colour maps for flux and residuals retain their meaning from Fig. 6 (https://github.com/aymgal/SLITronomy-papers/blob/master/paper_I/SLACS_source_reconstruction.ipynb).