Liisa Hirvonen edited Method.tex  almost 9 years ago

Commit id: baacb3dc142b15ea16bf3bd031f3a349af6b8109

deletions | additions      

       

The frames containing single photon events were processed with the ThunderSTORM \cite{Ovesny2014} superresolution imaging plug-in for ImageJ, whose settings are explored in this paper. Due to computational memory restraints, the USAF test pattern data was processed in six 5,000 image stacks. The data is first processed to detect the photon events from the noise background, before implementing a localization algorithm to detect the centre of the photon event. Subsequently, a sub-pixel localization algorithm is used on the data to fragment the pixels into 13$\times$13 grids to allow the centre of photon events to be located with greater resolution.  The camera default  parameters were set to a pixel used for the camera settings (pixel  size of 80.0nm, 80.0~nm,  3.6 photoelectrons per A/D count anda  base level of 414.0 A/D counts. counts).  For all the data, a wavelet (b-spline) image filter was applied with a b-spline order of 3 and a b-spline scale of 2.0. For the approximate localization algorithms, centroid of connected components (CoCC) was used with a peak intensity threshold (PIT) of 2*std(Wave.F1) and the watershed algorithm enabled, local maximum (LM) was applied using a PIT of 2.5*std(Wave.F1) with a 4-neighbourhood connectivity and non-maximum suppression (NMS) was used with a PIT of 2*std(Wave.F1) and 1 pixel dilation radius. For the sub-pixel localization algorithms, Gaussian  PSFgaussian  was set to a 2 pixel fitting radius, a 1.0 pixel standard deviation and the Maximum Likelihood (ML) fitting method;PSF  integrated gaussian Gaussian  (IG) PSF  was used with a 3 pixel fitting radius, a 1.6 pixel standard deviation and the least squares (LS) fitting method, radial symmetry was applied with a 2 pixel estimation radius and the centroid of local neighbourhood was set to a 1 pixel estimation radius. When the multiple-emitter fitting analysis (MFA) was applied in conjunction to the PSF gaussian or IG algorithms, it was used with a maximum of 1 molecule per fitting region, a model selection threshold (p-value) of 10-6 10^{-6}  and a limited intensity range of 1 to 800 photons.