Fluorescence microscopy method comprising the following steps: A. illuminating a scattering sample including one or more fluorescent emitters with a coherent excitation light beam having a wavefront with a first initial configuration C 1 for exciting said one or more fluorescent emitters comprised in the scattering sample; B. acquiring a first initial image F 1 corresponding to a background image F 1 = F B of the scattering sample, through a pixel imaging unit, obtaining a speckle grain image; C. selecting M target pixels P N, with N = 1,...,M, of coordinates [x N, y N] belonging to an area of the first initial image F 1 having fluorescence signals of at least one fluorescent emitter, each one of the M target pixels P N having a first initial intensity I PN 1; D. optimising the first initial configuration of the wavefront for each one of the M target pixels P N, for decreasing the speckle grain size of the first initial image F 1 and obtaining a final image F fin (I PN max) of the scattering sample with a local maximum I PN max of the fluorescence signal corresponding to a maximum intensity of the target pixels P N; E. subtracting the first initial image F B of background, obtained at step B, from the final image F fin (I PN max), obtained at step D, for each one of the M target pixels P N, obtaining an image F speckle having only fluorescent signals coming from a speckle grain located at pixel [x N, y N] with intensity I PN speckle, for each one of the M target pixels P N; F. fitting the image F speckle, obtained at step E, with a Gaussian function with free center coordinates, for each one of the M target pixels P N, thereby obtaining coordinates ( X N, Y N) and intensity I N in a xy plane; G. generating an image F N containing a Gaussian distribution centered at coordinates ( X N, Y N) and with intensity I N, obtained at step F, and with a waist equal to an average size S of speckle grain of the scattering sample, for each one of the M target pixels P N; H. generating a final image of the sample by summing the M images F N obtained at step G, for each one of the M target pixels P N. A. illuminating a scattering sample including one or more fluorescent emitters with a coherent excitation light beam having a wavefront with a first initial configuration C 1 for exciting said one or more fluorescent emitters comprised in the scattering sample; B. acquiring a first initial image F 1 corresponding to a background image F 1 = F B of the scattering sample, through a pixel imaging unit, obtaining a speckle grain image; C. selecting M target pixels P N, with N = 1,...,M, of coordinates [x N, y N] belonging to an area of the first initial image F 1 having fluorescence signals of at least one fluorescent emitter, each one of the M target pixels P N having a first initial intensity I PN 1; D. optimising the first initial configuration of the wavefront for each one of the M target pixels P N, for decreasing the speckle grain size of the first initial image F 1 and obtaining a final image F fin (I PN max) of the scattering sample with a local maximum I PN max of the fluorescence signal corresponding to a maximum intensity of the target pixels P N; E. subtracting the first initial image F B of background, obtained at step B, from the final image F fin (I PN max), obtained at step D, for each one of the M target pixels P N, obtaining an image F speckle having only fluorescent signals coming from a speckle grain located at pixel [x N, y N] with intensity I PN speckle, for each one of the M target pixels P N; F. fitting the image F speckle, obtained at step E, with a Gaussian function with free center coordinates, for each one of the M target pixels P N, thereby obtaining coordinates ( X N, Y N) and intensity I N in a xy plane; G. generating an image F N containing a Gaussian distribution centered at coordinates ( X N, Y N) and with intensity I N, obtained at step F, and with a waist equal to an average size S of speckle grain of the scattering sample, for each one of the M target pixels P N; H. generating a final image of the sample by summing the M images F N obtained at step G, for each one of the M target pixels P N.