We reconstructed the super-resolution image with a single-emitter fitting method (ThunderSTORM), a multi-emitter fitting method (3D-DAOSTORM), and a compressed sensing-based approach (FALCON), all compared with WindSTORM. The single-emitter maximum likelihood Gaussian function fitting algorithm and wavelet filtering were used in ThunderSTORM, and the multi-emitter maximum likelihood Gaussian function fitting algorithm was used in 3D-DAOSTORM. The kernel width of the PSF used in WindSTORM was set to 1.5 pixels for the simulated dataset and our experimental PALM imaging dataset, 1.9 pixels for our experimental STORM imaging dataset, and 1.4 pixels for the open-access experimental dataset (http://bigwww.epfl.ch/smlm/datasets/index.html?p=real-hd). For iterative approaches (FALCON and ThunderSTORM), the above kernel widths of the PSF were used as the initial kernel. The drift correction was performed using cross-correlation method provided in ThunderSTORM (30 (link)). In this study, the GPU versions of WindSTORM and FALCON were executed on GPU (GTX1080, Nvidia), and the CPU version of WindSTORM, 3D-DAOSTORM, FALCON, and ThunderSTORM was executed on a Quad-core CPU (Core i7-4790, Intel).