Fast live cell imaging at nanometer scale using annihilating filter-based low-rank Hankel matrix approach

Localization microscopy such as STORM/PALM can achieve a nanometer scale spatial resolution by iteratively localizing fluorescence molecules. It was shown that imaging of densely activated molecules can accelerate temporal resolution which was considered as major limitation of localization microscopy. However, this higher density imaging needs to incorporate advanced localization algorithms to deal with overlapping point spread functions (PSFs). In order to address this technical challenges, previously we developed a localization algorithm called FALCON1, 2 using a quasi-continuous localization model with sparsity prior on image space. It was demonstrated in both 2D/3D live cell imaging. However, it has several disadvantages to be further improved. Here, we proposed a new localization algorithm using annihilating filter-based low rank Hankel structured matrix approach (ALOHA). According to ALOHA principle, sparsity in image domain implies the existence of rank-deficient Hankel structured matrix in Fourier space. Thanks to this fundamental duality, our new algorithm can perform data-adaptive PSF estimation and deconvolution of Fourier spectrum, followed by truly grid-free localization using spectral estimation technique. Furthermore, all these optimizations are conducted on Fourier space only. We validated the performance of the new method with numerical experiments and live cell imaging experiment. The results confirmed that it has the higher localization performances in both experiments in terms of accuracy and detection rate.

[1]  Keith A. Lidke,et al.  Simultaneous multiple-emitter fitting for single molecule super-resolution imaging , 2011, Biomedical optics express.

[2]  Suliana Manley,et al.  3D high-density localization microscopy using hybrid astigmatic/ biplane imaging and sparse image reconstruction. , 2014, Biomedical optics express.

[3]  J. Lippincott-Schwartz,et al.  Imaging Intracellular Fluorescent Proteins at Nanometer Resolution , 2006, Science.

[4]  Jong Chul Ye,et al.  A General Framework for Compressed Sensing and Parallel MRI Using Annihilating Filter Based Low-Rank Hankel Matrix , 2015, IEEE Transactions on Computational Imaging.

[5]  Jong Chul Ye,et al.  Annihilating Filter-Based Low-Rank Hankel Matrix Approach for Image Inpainting , 2015, IEEE Transactions on Image Processing.

[6]  Yingbo Hua.  Estimating two-dimensional frequencies by matrix enhancement and matrix pencil , 1992, IEEE Trans. Signal Process..

[7]  José M. Bioucas-Dias,et al.  Fast Image Recovery Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

[8]  Michael J Rust,et al.  Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM) , 2006, Nature Methods.

[9]  S. Holden,et al.  DAOSTORM: an algorithm for high- density super-resolution microscopy , 2011, Nature Methods.

[10]  Suliana Manley,et al.  Quantitative evaluation of software packages for single-molecule localization microscopy , 2015, Nature Methods.

[11]  Nicolas Olivier,et al.  FALCON: fast and unbiased reconstruction of high-density super-resolution microscopy data , 2014, Scientific Reports.

[12]  Marc Moonen,et al.  An efficient subspace algorithm for 2-D harmonic retrieval , 1994, Proceedings of ICASSP '94. IEEE International Conference on Acoustics, Speech and Signal Processing.

[13]  B. Krauskopf,et al.  Proc of SPIE , 2003 .

[14]  Lei Zhu,et al.  Faster STORM using compressed sensing , 2012, Nature Methods.

[15]  Michael D. Mason,et al.  Ultra-high resolution imaging by fluorescence photoactivation localization microscopy. , 2006, Biophysical journal.

[16]  Nathan Halko,et al.  Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions , 2009, SIAM Rev..