Next Article in Journal
Intensity-Averaged Double Three-Step Phase-Shifting Algorithm with Color-Encoded Fringe Projection
Next Article in Special Issue
Blood Cell Analysis: From Traditional Methods to Super-Resolution Microscopy
Previous Article in Journal
Microwave Photonic Signal Generation in an Optically Injected Discrete Mode Semiconductor Laser
Previous Article in Special Issue
A Review of Raman-Based Technologies for Bacterial Identification and Antimicrobial Susceptibility Testing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

sCMOS Noise-Corrected Superresolution Reconstruction Algorithm for Structured Illumination Microscopy

1
State Key Laboratory of Membrane Biology, Beijing Key Laboratory of Cardiometabolic Molecular Medicine, Institute of Molecular Medicine, School of Future Technology, Peking University, Beijing 100871, China
2
Biomedical Engineering Department, Peking University, Beijing 100871, China
3
Chongqing Key Laboratory of Image Cognition, College of Computer Science and Technology, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
4
PKU-IDG/McGovern Institute for Brain Research, Beijing 100871, China
5
Beijing Academy of Artificial Intelligence, Beijing 100871, China
6
Shenzhen Bay Laboratory, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(3), 172; https://doi.org/10.3390/photonics9030172
Submission received: 23 February 2022 / Revised: 6 March 2022 / Accepted: 7 March 2022 / Published: 10 March 2022

Abstract

:
Structured illumination microscopy (SIM) is widely applied due to its high temporal and spatial resolution imaging ability. sCMOS cameras are often used in SIM due to their superior sensitivity, resolution, field of view, and frame rates. However, the unique single-pixel-dependent readout noise of sCMOS cameras may lead to SIM reconstruction artefacts and affect the accuracy of subsequent statistical analysis. We first established a nonuniform sCMOS noise model to address this issue, which incorporates the single-pixel-dependent offset, gain, and variance based on the SIM imaging process. The simulation indicates that the sCMOS pixel-dependent readout noise causes artefacts in the reconstructed SIM superresolution (SR) image. Thus, we propose a novel sCMOS noise-corrected SIM reconstruction algorithm derived from the imaging model, which can effectively suppress the sCMOS noise-related reconstruction artefacts and improve the signal-to-noise ratio (SNR).

1. Introduction

Superresolution (SR) microscopy enables biological researchers to see nanoscale images of intracellular structures. Various SR fluorescence microscopy techniques, such as stimulated emission depletion (STED) [1,2,3], photoactivated localization microscopy (PALM) [4,5,6], stochastic optical reconstruction microscopy (STORM) [7,8,9], and structured illumination microscopy (SIM) [10,11,12,13,14,15], have come to the fore during the past 20 years. These SR technologies can break the optical diffraction limit and achieve a spatial resolution of approximately 20~100 nm compared to the ~200 nm of conventional microscopes. SR SIM has a relatively high temporal resolution and long-term imaging ability and is a promising SR imaging technique for live-cell imaging [16,17,18].
Although an electron-multiplying charge-coupled device (EMCCD) can detect signals with sufficient SNR and sensitivity for SR SIM, its readout speed is slow and it has small sensor areas [19]. Therefore, scientific-grade complementary metal-oxide-semiconductor (sCMOS) cameras are a better choice for offering sufficient quantum efficiency and much faster readout speed, significantly increasing the data acquisition rate and improving the temporal resolution [20,21,22]. With the sCMOS camera, Hessian-SIM was developed to increase the temporal resolution and image rapidly moving vesicles or loops in the endoplasmic reticulum with a temporal resolution of 188 Hz [23].
Generally, images captured from sCOMS contain two kinds of typical noises: shot noise and readout noise. For the former, which comes from the arrival of photons at the camera sensor and is a stochastic process, we usually constitute shot noise following Poisson counting statistics. Shot noise equals the square root of the input signal and is always present due to its physical essence. For the latter, the readout noise comes from nonuniform pixel characteristics due to the process variation in the on-pixel amplifier circuitry and column-based readout structures [20]. This results in unique sCMOS pixel-dependent noise, which negatively impacts the image reconstruction and creates single-molecule localization errors in PALM/STORM and artefacts in SIM reconstruction. To correct the pixel-independent readout noise in the raw images, Huang et al. proposed a simplified sCMOS noise model incorporating single-pixel-dependent offset, gain, and variance. They developed a noise-corrected SMSN reconstruction algorithm [24] and a noise-correction algorithm for an sCMOS camera for a broad spectrum of microscopes [25]. By combining camera physics and sparsely layered filtering, Biagio et al. proposed a content-adaptive algorithm to automatically correct sCMOS-relatedd noise (ACsN) for fluorescence microscopy, which reduces the most relevant noise sources in the sCMOS camera while preserving the fine details of the signal [26]. Xue et al. developed a Hessian-based SMLM (Hessian-SMLM) method to correct the single-pixel variance, gain, and offset and effectively eliminated the pixel-dependent readout noise, especially under conditions of low signal-to-noise ratios [27]. Lin et al. presented a nonuniform noise model of sCMOS cameras that incorporates pixel-specific read-noise, offset, and sensitivity variation. Thus, they developed a new weighted least squared (WLS) fitting method designed to remove the heterogeneity of sCMOS pixels [20]. Li et al. developed the photon transfer curve (PTC) method to fully assess the performance of low-light cameras, including the sCMOS camera [28].
Despite all of these correction algorithms developed for the sCMOS cameras, correcting artefacts due to heterogeneous pixel noise in SIM reconstruction remains unexplored. However, the different offset and variance in each pixel causes artefacts in reconstructed SIM images (especially in hot pixels). To minimize these artefacts caused by readout noise, we propose an sCMOS noise model based on SIM imaging and present a novel noise-corrected algorithm for SIM reconstruction. The reconstruction is conducted in the spatial domain based on the stable bi-conjugate gradients descent algorithm (Bi-CGSTAB) [29] and split Bergman algorithm (SB) [30]. In all simulation experiments, the algorithm can effectively suppress the sCMOS pixel-dependent readout noise-related reconstruction artefacts, improve the SNR, and slightly increase the contrast of the reconstructed SR image using various regularized constraints.

2. Methods

2.1. SIM Imaging Model

Compared to conventional wide-field fluorescence microscopy, SIM can double the spatial resolution by illuminating the observed sample with pattern illumination generated by the interference of two beams of light (Figure 1a). Without considering noise, the fluorescence emission distribution I n detected by the camera can be expressed as:
I n ( r k ) = [ s ( r l ) p n ( r l ) ] h ( r k r l ) ,
where r l and r k are the spatial coordinate vectors of the sample and camera, respectively, s ( r l ) is the fluorophore spatial distribution in the sample, p n ( r l ) is the intensity distribution according to the illumination pattern sequence n , h ( r k r l ) is the point spread function, and represents the convolutional operation. The model can be discretized and expressed as:
I k n = l s l p l n h k l
where l and   k are the discrete spatial coordinates as substitutes for the continuous spatial coordinate vectors r l and r k .

2.2. sCMOS Camera Noise Model

In SIM imaging, the acquired raw image contains noise from the camera. sCMOS camera noise primarily consists of shot noise and readout noise. While shot noise stems from the photon detection process, readout noise originates from the electronics built around the detector chip [24]. When photons are detected on the sensor chip, the camera’s analogue-to-digital unit (ADU) count output follows a probability distribution described here by the convolution of a Poisson distribution and a Gaussian distribution. The Poisson distribution represents the shot noise of photon detection, and the Gaussian distribution results from the readout noise (Figure 1b). Then, the pixel-dependent conditional probability density function (PDF) of the sCMOS camera can be expressed as:
P ( d k n | I k n ) = A q k n = 0 1 q k n e I k n ( I k n ) q k n 1 2 π σ k e ( d k n g k q k n o k ) 2 2 σ k 2
where d k n is the specific counts obtained by the sCMOS in the kth pixel when the sample is illuminated by the nth pattern (in units of ADU ). A is the normalization constant, q k n is the number of fluorescence photoelectrons as a random variable (in units of e ), and I k n is the number of expected fluorescence photoelectrons (in units of e ) and defined by Equation (2), σ k 2 is the readout noise variance in the kth pixel (in units of ADU 2 ), g k is the amplification gain of the kth pixel (in units of ADU / e ), and o k is the readout noise mean (offset) of the kth pixel (in units of ADU ).
For the convenience of later deduction, the distribution of random variables d k n can be equivalently expressed as:
d k n = P k n + G k + o k
where random variable P k n follows the Poisson distribution with mean value and variance value equal to g k I k n :
P k n ~ P ( g k I k n )
Random variable G k follows the Gaussian distribution with mean value equal to zero and variance value equal to σ k 2 :
G k ~ N ( 0 , σ k 2 )
Due to the noise, pixels on sCMOS cameras appear to flicker even when there are no expected incident photons. This noise level changes drastically across pixels from 1~2 ADU 2 to 1000~2000 ADU 2 [25]. This noise drastically reduces the image SNR and makes quantitative analysis a challenge.

2.3. sCMOS Camera Characterization

To characterize the readout noise of the sCMOS camera, the offset, gain, and variance values need to be measured (Figure 2). The offset value describes a constant level of ADUs pre-engineered into the readout process to prevent negative ADUs caused by the readout noise [24]. After acquiring a series of images in an environment with zero expected incident photons, the offset value o k and the variance value σ k 2 of the kth sCMOS pixel can be calculated as:
o k = 1 M m = 1 M a k m
σ k 2 = 1 M m = 1 M ( a k m ) 2 ( 1 M m = 1 M a k m ) 2
where a k m is the ADU count at the mth frame for the kth pixel and M is the total number of dark frames acquired.
By illuminating the camera with quasi-uniform stationary intensity patterns and recording a series of image sequences at different average intensity levels, the gain value g k of the kth sCMOS pixel can be calculated as:
g k = a r g m i n g i { [ ( σ k i ) 2 σ k 2 ] g ( o k i o k ) }
where:
o k i = 1 M i m = 1 M i b k i m
The variance value ( σ k i ) 2 can be defined as:
( σ k i ) 2 = 1 M i m = 1 M i ( b k i m ) 2 ( 1 M i m = 1 M i b k i m ) 2
where b k i m is the ADU count at the ith illuminating intensity and the mth frame for the kth pixel, and M i is the total frame number acquired at the ith illuminating intensity.

2.4. sCMOS Noise-Corrected SIM Reconstruction Algorithm

Based on the deduction in Section 2.2., a new random variable Z k n can be defined and expressed as:
Z k n = d k n o k g k + σ k 2 g k 2
The mean value of the Gaussian distribution equals its variance value. When the mean of a Poisson distribution is large, it becomes similar to a Gaussian distribution. Then, the distribution of random variable Z k n approximately follows the Gaussian distribution and can be expressed as:
Z k n ~ N ( σ k 2 g k 2 + I k n , σ k 2 g k 2 + I k n )
For simplification, the I k n term in the variance value of Z k n is ignored, and then the approximate CPDF of Z k n can be expressed as:
P ( Z k n | I k n ) = g k 2 π σ k e ( g k 2 Z k n g k 2 I k n σ k 2 ) 2 2 g k 2 σ k 2
The approximate CPDF of d k n can be calculated by substituting Z k n into Equation (14) with Equation (12) and expressed as:
P ( d k n | I k n ) = g k 2 π σ k e ( d k n g k I k n o k ) 2 2 σ k 2
The sample fluorophore spatial distribution s is determined by maximizing the posterior probability P ( s | d ) defined as:
P ( s | d ) = n k P ( s | d k n )
where P ( s | d k n ) can be calculated with Bayes’ rule and expressed as:
P ( s | d k n ) = P ( d k n | s ) P ( s ) P ( d k n ) P ( d k n | s )
For operational convenience, the merit function E ( s ) is defined as the logarithm of P ( s | d ) and calculated by combining Equations (2) and (15)–(17). The simplified form of E ( s ) can be expressed as:
E ( s ) = n k [ g k 2 σ k 2 ( d k n o k g k l s l p l n h k l ) 2 ]
The regularization constraint can assist optimization convergence and prevent overfitting. The Tikhonov regularization term R T ( s ) , high-frequency suppression regularization term R H P ( s ) , and nonnegativity regularization term R N N ( s ) are appended to the merit function, and the whole constrained optimization problem can be expressed as:
m i n s n k [ g k 2 σ k 2 ( d k n o k g k l s l p l n h k l ) 2 ] + α R T ( s ) + β R H P ( s ) + R N N ( s )
where α and β are the weights of the corresponding regularization terms, and the definition of each regularization term can be expressed as:
R T ( s ) = l w s l 2
R H P ( s ) = k ( l s l h k l h p ) 2
R N N ( s ) = ( s 0 )
where w is the parameter of the Tikhonov regularization term and h h p is the parameter of the high-frequency suppression regularization term as a tunable high-pass filter. The matriculated representation of Equation (19) can be expressed as:
m i n s n A n s b n 2 2 + α s T C s + β s T D T D s + ( s 0 )
where s is the fluorophore spatial distribution in the sample as a column vector, A n is the illumination by the nth pattern and PSF convolution operation as a matrix, b n is the image obtained by the sCMOS as a column vector, C and D are the Tikhonov parameter matrix and high-pass filter matrix, respectively, and [ ] T represents the conjugate transpose operation.
The constrained optimization problem defined by Equation (23) can be solved by the combination of the stable bi-conjugate gradients descent algorithm (Bi-CGSTAB) and the split Bergman algorithm (SB). The Bi-CGSTAB algorithm can minimize convex differentiable functional such as the L2-functional in Equation (23) and avoid severe cancellation effects caused by the irregular convergence behavior of the conventional gradients’ descent algorithm [29]. The SB algorithm allows for the minimization of convex nondifferentiable functional efficiency and leads to a solution update for which L2-functional and the nonnegativity functions are decoupled and solved separately [30]. Using the Bergman method, the equivalent formulation of Equation (23) can be expressed as:
( s k + 1 , v k + 1 ) = a r g m i n s , v n A n s b n k 2 2 + δ s v + μ v k 2 2 + α s T C s + β s T D T D s + ( v 0 ) b n k + 1 = b n k + b n A n s k + 1 μ v k + 1 = μ v k + s k + 1 v k + 1
with b n 0 = b n , s 0 = 0 , v 0 = 0 . The variables s k + 1 and v k + 1 can be solved separately due to the decoupled relationship between them:
s k + 1 = a r g m i n s n A n s b n k 2 2 + δ s v k + μ v k 2 2 + α s T C s + β s T D T D s
v k + 1 = a r g m i n v δ s k + 1 v + μ v k 2 2 + ( v 0 )
The elements of v k + 1 in Equation (26) are solved independently using a shrinkage formula due to the decoupled relationship between them:
v k + 1 = m a x ( s k + 1 + μ v k , 0 )
To solve the variable s k + 1 in Equation (25) with the Bi-CGSTAB algorithm, the merit function can be expanded as:
a r g m i n s s T M s s T n + p
where M is a matrix, n is a column vector, and p is a scalar and can be calculated as:
M = n A n T A n + δ I + α C + β D T D
n = n A n T b n k + δ ( v k μ v k )
p = n { ( b n k ) T A n s + δ ( b n k ) T b n k [ ( μ v k ) T ( v k ) T ] ( s v k + μ v k ) }
where I is the unit matrix. According to convex optimization theory, the solution s ^ of Equation (28) satisfies the following equation:
M s ^ = n
which can be solved by the Bi-CGSTAB algorithm (Algorithm 1).
Algorithm 1. Pseudocode of the Bi-CGSTAB algorithm.
Bi-CGSTAB algorithm
Input: M , n
Initialization:
s 0 = 0 , r 0 = n M s 0 , r ^ 0 = r 0 , ρ 0 = ω 0 = α = 1 , v 0 = p 0 = 0
Iteration:
For k = 1 , 2 , 3  
ρ k = ( r ^ 0 , r k 1 ) , β = ( α ρ k ) / ( ρ k 1 ω k 1 )
p k = r k 1 + β ( p k 1 ω k 1 v k 1 )
v k = M p k
α = ρ k / ( r ^ 0 , v k )
b = r k 1 α v k
t = M b
ω k = ( t , b ) / ( t , t )
s k = s k 1 + α p k + ω k b
if s k s k 1 2 2 δ
Terminating.
Else
r k = b ω k t
End
Output: s = s k

3. Results

Compared with the conventional Wiener algorithm, and not considering the sCMOS noise, our proposed sCMOS noise-corrected SIM reconstruction algorithm can suppress the sCMOS noise-related spatial punctiform artefact in the SR reconstructed image, especially in the background area (Figure 3a–c). The punctiform artefact within the SR image reconstructed by the conventional algorithm is widely distributed (Figure 3a) and can be observed within the zoomed-in image (left in Figure 3c). After the proposed sCMOS noise-correction algorithm, the punctiform artefact is eliminated globally (Figure 3b) and invisible to human eyes (right in Figure 3c). The statistical box plot (Figure 3d) shows that the SNR of the SR image reconstructed by the conventional algorithm is ~16.8 and improved to ~17.7 by the proposed sCMOS noise-corrected SIM reconstruction algorithm, which is consistent with the phenomenon of punctiform artefact elimination above.
Furthermore, the intensity fluctuation standard deviation in each pixel over time, which reflects the temporal pixel fluctuation map, is shown in Figure 3e–g. The bright and punctiform structures within the temporal pixel fluctuation map corresponding to the conventional algorithm (Figure 3e and left in Figure 3g) indicate that the intensity fluctuation standard deviation of the corresponding pixels is relatively large. That is, the temporal fluctuation noise level of these pixels is relatively high. These structures reflecting high noise levels almost completely disappear (Figure 3f, and right in Figure 3g) with the proposed sCMOS noise-corrected SIM reconstruction algorithm. Hence, the temporal fluctuation noise is suppressed effectively. The statistical histogram of the temporal pixel fluctuation map (Figure 3h) confirms the conclusion above. The pixel intensity fluctuation standard deviation distribution corresponding to the conventional algorithm moves left globally after the sCMOS noise correction. The mode interval of the pixel intensity fluctuation standard deviation corresponding to the noise-corrected algorithm is ~4.75, which is 30% lower than the ~6.75 of the conventional algorithm.
The punctiform artefact and high pixel intensity fluctuation standard deviation may result in negative impacts on the data analysis. For example, the dynamic analysis of the organelle structures and functions in the living cell is likely to be misleading when the sCMOS noise-related pixel intensity fluctuation is severe in the SR reconstructed images. Such possible mistakes can be avoided by the proposed sCMOS noise-corrected SIM reconstruction algorithm.
To further verify the robustness of our method, we validated our approach with different SNRs. We found that the SNR of the SR image reconstructed by the proposed algorithm is robust under various shot noise and readout noise levels (Table 1). The average improvement percentage of the SNR is ~22%, and the maximum improvement percentage of the SNR is ~30% when the shot noise level and the readout noise level are equal to 100 and 30, respectively. The robustness of the proposed algorithm guarantees broad applicability in the reconstruction of the raw data at various noise levels. The simulated ground truth image is identical to Figure 3. The shot noise level is quantified by the parameter representing the total photon count of each emitter on the simulated samples. The read out noise level is quantified by the parameter multiplied to the variance map of the sCMOS camera. For example, if the read out noise level is set as L, the variance of the ith pixel in the sCMOS camera variance map is V. Then, the variance of Gauss noise is added to the ith pixel is set as LV when we generate the simulation SIM raw noisy data. For all of the reconstructions, only the Tikhonov regularized constraint is used and the regularized term weight parameter α is set as 2.
The proposed algorithm can also combine regularization to further suppress the noise. The SR images reconstructed by the noise-corrected SIM reconstruction algorithm with nonnegativity, high-frequency suppression, and Tikhonov regularized constraint are compared in Figure 4.
The ground truth image contains the chirped line structure (Figure 4a). The shot noise level and read out noise level of SIM raw data are 150 and 10, respectively. Within all of the reconstructed SR images (Figure 4b–d), the reconstruction with a nonnegativity regularized constraint exhibits the highest contrast (Figure 4b). The grey value range of the profile corresponding to the nonnegativity regularized constraint is from ~96 to ~174 (Figure 4f), which is ~30% and ~16% wider than the high-frequency suppression regularized constraint (Figure 4g) and Tikhonov regularized constraint (Figure 4h). However, the reconstruction with the nonnegativity regularized constraint still produces severe artefacts. Thus, the nonnegativity regularized constraint is usually not used alone as denoising regularization in the reconstruction of low SNR raw images.
In contrast, the reconstruction of the high-frequency suppression regularized constraint has the lowest artefact level compared with those of the other two methods (Figure 4c). The smoothness of profiles provides additional evidence of artefact suppression (Figure 4g). However, the two lines are indistinguishable in reconstructing the high-frequency suppression method (Figure 4c, white arrow), which means the resolution decreased for reconstructing the high-frequency suppression regularized constraint. In addition, the average full width at half maximum (FWHM) of the ten spikes on the right of the profile corresponding to the high-frequency suppression regularized constraint (Figure 4g) is ~3 pixels, which is ~33% larger than the ground truth profile (Figure 4e).
The resolution, artefact level, and contrast level of the SR image reconstructed with Tikhonov regularization (Figure 4d) are the most moderate within all three reconstructed SR images (Figure 4b–d). Additionally, the profile is also the closest to the ground truth profile (Figure 4e,h). Therefore, the Tikhonov regularized constraint is a more general choice than the nonnegativity regularized constraint and high-frequency suppression regularized constraint.

4. Discussion

Although quite a few sCMOS noise-correction algorithms for various imaging technologies have been developed, the sCMOS noise-correction algorithm in SIM reconstruction has never been presented to date as far as we know. The experiment shows that the proposed sCMOS noise-corrected SIM reconstruction algorithm can suppress the sCMOS readout noise-related spatial punctiform artefact and temporal fluctuation noise caused in the SR reconstructed images and improve the SNR by 22%, which is validated by the simulation experiments.
We combined the proposed sCMOS noise-corrected reconstruction with three regular terms. The combination with the nonnegativity regularized constraint can be suitable for raw data with a relatively high SNR, while the combination with high-frequency suppression regularization is more appropriate for raw data with a low SNR. As the most general method, the sCMOS noise-corrected SIM reconstruction algorithm with Tikhonov regularization provides the most moderate reconstructed SR images in terms of the resolution, artefact level, and contrast level.

5. Conclusions

A nonuniform sCMOS noise model in the SIM imaging process is established and a novel sCMOS noise-corrected SIM reconstruction algorithm is proposed based on the imaging noise model. The sCMOS noise causes artefacts in the reconstructed SR image and the novel algorithm can effectively suppress the artefacts and improve the SNR.

Author Contributions

Conceptualization, X.H., J.F. and B.Z.; methodology, J.F. and B.Z.; software, B.Z. and J.F.; validation, L.C., J.F., X.H. and B.Z.; formal analysis, B.Z. and J.F.; investigation, B.Z.; resources, L.C. and J.F.; data curation, B.Z.; writing—original draft preparation, B.Z. and J.F.; writing—review and editing, L.C., J.F. and B.Z.; visualization, B.Z. and J.F.; supervision, L.C., J.F. and X.H.; project administration, L.C. and J.F.; funding acquisition, L.C. and J.F. All authors have read and agreed to the published version of the manuscript.

Funding

Please add: This research was funded by the National Natural Science Foundation of China grant nos. (81925022, 92054301, 91750203, 31821091, 62103071), the National Key Research and Development Programme of China grant nos. (SQ2016YFJC040028), the Beijing Natural Science Foundation grant nos. (Z200017, Z201100008420005, Z20J00059), the National Science and Technology Major Project Programme grant nos. (2016YFA0500400), the Natural Science Foundation of Chongqing grant nos. (cstc2021jcyj-msxmX0526), and the Science and Technology Research Programme of Chongqing Municipal Education Commission grant nos. (Grant No. KJQN202100630).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hell, S.; Wichmann, J. Breaking the Diffraction Resolution Limit by Stimulated-Emission: Stimulated-Emission-Depletion Fluorescence Microscopy. Opt. Lett. 1994, 19, 780–782. [Google Scholar] [CrossRef] [PubMed]
  2. Dyba, M.; Jakobs, S.; Hell, S. Immunofluorescence stimulated emission depletion microscopy. Nat. Biotechnol. 2003, 21, 1303–1304. [Google Scholar] [CrossRef] [PubMed]
  3. Gao, P.; Prunsche, B.; Zhou, L.; Nienhaus, K.; Nienhaus, G.U. Background suppression in fluorescence nanoscopy with stimulated emission double depletion. Nat. Photonics 2017, 11, 163–169. [Google Scholar] [CrossRef]
  4. Betzig, E.; Patterson, G.H.; Sougrat, R.; Lindwasser, O.W.; Olenych, S.; Bonifacino, J.S.; Davidson, M.W.; Lippincott-Schwartz, J.; Hess, H.F. Imaging Intracellular Fluorescent Proteins at Nanometer Resolution. Science 2006, 313, 1642–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hess, S.T.; Girirajan, T.P.K.; Mason, M.D. Ultra-high resolution imaging by fluorescence photoactivation localization microscopy. Biophys. J. 2006, 91, 4258–4272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Errico, C.; Pierre, J.; Pezet, S.; Desailly, Y.; Lenkei, Z.; Couture, O.; Tanter, M. Ultrafast ultrasound localization microscopy for deep super-resolution vascular imaging. Nature 2015, 527, 499–502. [Google Scholar] [CrossRef]
  7. Huang, B.; Wang, W.; Bates, M.; Zhuang, X. Three-dimensional super-resolution imaging by stochastic optical reconstruction microscopy. Science 2008, 319, 810–813. [Google Scholar] [CrossRef] [Green Version]
  8. Linde, S.; Löschberger, A.; Klein, T.; Heidbreder, M.; Wolter, S.; Heilemann, M.; Sauer, M. Direct stochastic optical reconstruction microscopy with standard fluorescent probes. Nat. Protoc. 2011, 6, 991–1009. [Google Scholar] [CrossRef]
  9. Rust, M.J.; Bates, M.; Zhuang, X. Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM). Nat. Methods 2006, 3, 793–795. [Google Scholar] [CrossRef] [Green Version]
  10. Gustafsson, M.G.L. Surpassing the lateral resolution limit by a factor of two using structured illumination microscopy. Short Commun. J. Microsc. 2000, 198, 82–87. [Google Scholar] [CrossRef] [Green Version]
  11. Cragg, G.E.; So, P. Lateral resolution enhancement with standing evanescent waves. Opt. Lett. 2000, 25, 46–48. [Google Scholar] [CrossRef] [PubMed]
  12. Frohn, J.T.; Knapp, H.F.; Stemmer, A. True optical resolution beyond the Rayleigh limit achieved by standing wave illumination. Proc. Natl. Acad. Sci. USA 2000, 97, 7232–7236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Günther, B.; Hehn, L.; Jud, C.; Hipp, A.; Dierolf, M.; Pfeiffer, F. Full-field structured-illumination super-resolution X-ray transmission microscopy. Nat. Commun. 2019, 10, 2494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. York, A.; Parekh, S.; Nogare, D.; Fischer, R.; Temprine, K.; Mione, M.; Chitnis, A.; Combs, C.; Shroff, H. Resolution doubling in live, multicellular organisms via multifocal structured illumination microscopy. Nat. Methods 2012, 9, 749–754. [Google Scholar] [CrossRef]
  15. Mudry, E.; Belkebir, K.; Girard, J.; Savatier, J.; Le Moal, E.; Nicoletti, C.; Allain, M.; Sentenac, A. Structured illumination microscopy using unknown speckle patterns. Nat. Photonics 2012, 6, 312–315. [Google Scholar] [CrossRef]
  16. Kner, P.; Chhun, B.B.; Griffis, E.R.; Winoto, L.; Gustafsson, M.G.L. Super-resolution video microscopy of live cells by structured illumination. Nat. Methods 2009, 6, 339–342. [Google Scholar] [CrossRef] [Green Version]
  17. Shao, L.; Kner, P.; Rego, E.H.; Gustafsson, M.G.L. Super-resolution 3D microscopy of live whole cells using structured illumination. Nat. Methods 2011, 12, 1044–1046. [Google Scholar] [CrossRef]
  18. Fiolka, R.; Shao, L.; Rego, E.H.; Davidson, M.W.; Gustafsson, M.G.L. Time-lapse two-color 3D imaging of live cells with doubled resolution using structured illumination. Proc. Natl. Acad. Sci. USA 2012, 109, 5311–5315. [Google Scholar] [CrossRef] [Green Version]
  19. Diekmann, R.; Till, K.; Müller, M.; Simonis, M.; Schüttpelz, M.; Huser, T. Characterization of an industry-grade CMOS camera well suited for single molecule localization microscopy-high performance super-resolution at low cost. Sci. Rep. 2017, 7, 14425. [Google Scholar] [CrossRef]
  20. Lin, R.; Clowsley, A.H.; Jayasinghe, I.D.; Baddeley, D.; Soeller, C. Algorithmic corrections for localization microscopy with sCMOS cameras-characterisation of a computationally efficient localization approach. Opt. Express 2017, 25, 11701–11716. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, Z.; Wang, Y.; Piestun, R.; Huang, Z. Characterizing and correcting camera noise in back-illuminated sCMOS cameras. Opt. Express 2021, 29, 6668–6690. [Google Scholar] [CrossRef] [PubMed]
  22. Saurabh, S.; Maji, S.; Bruchez, M.P. Evaluation of sCMOS cameras for detection and localization of single Cy5 molecules. Opt. Express 2012, 20, 7338–7349. [Google Scholar] [CrossRef] [PubMed]
  23. Huang, X.; Fan, J.; Li, L.; Liu, H.; Wu, R.; Wu, Y.; Wei, L.; Mao, H.; Lal, A.; Xi, P.; et al. Fast, long-term, super-resolution imaging with Hessian structured illumination microscopy. Nat. Biotechnol. 2018, 36, 451–459. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, F.; Hartwich, T.; Rivera-Molina, F.; Lin, Y.; Duim, W.; Long, J.; Uchil, P.; Myers, J.; Baird, M.; Mothes, W.; et al. Video-rate nanoscopy using sCMOS camera-specific single-molecule localization algorithms. Nat. Methods 2013, 10, 653–658. [Google Scholar] [CrossRef]
  25. Liu, S.; Mlodzianoski, M.; Hu, Z.; Ren, Y.; McElmurry, K.; Suter, D.; Huang, F. sCMOS noise-correction algorithm for microscopy images. Nat. Methods 2017, 14, 760–761. [Google Scholar] [CrossRef]
  26. Mandracchia, B.; Hua, X.; Guo, C.; Son, J.; Urner, T.; Jia, S. Fast and accurate sCMOS noise correction for fluorescence microscopy. Nat. Commun. 2020, 11, 94. [Google Scholar] [CrossRef] [Green Version]
  27. Xue, F.; He, W.; Xu, F.; Zhang, M.; Chen, L.; Xu, P. Hessian single-molecule localization microscopy using sCMOS camera. Biophys. Rep. 2018, 4, 215–221. [Google Scholar] [CrossRef]
  28. Li, L.; Li, M.; Zhang, Z.; Huang, Z. Assessing low-light cameras with photon transfer curve method. J. Innov. Opt. Health Sci. 2016, 9, 1630008. [Google Scholar] [CrossRef] [Green Version]
  29. van der Vorst, H.A. Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems. Siam J. Sci. Stat. Comput. 1992, 13, 631–644. [Google Scholar] [CrossRef]
  30. Abascal, P.; Chamorro-Servent, J.; Aguirre, J. Fluorescence diffuse optical tomography using the split Bregman method. Med. Phys. 2011, 38, 6275–6284. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Diagram of sCMOS camera noise in SIM. (a) Schematic diagram of the SIM setup. PBS, polarization beam splitter; AOTF, acousto-optic tunable filters; HWP, half wave-plate; DM, dichroic mirror; SLM, spatial light modulator; PR, polarization rotator; L1–L5, lenses. (b) SIM imaging model and sCMOS camera noise model.
Figure 1. Diagram of sCMOS camera noise in SIM. (a) Schematic diagram of the SIM setup. PBS, polarization beam splitter; AOTF, acousto-optic tunable filters; HWP, half wave-plate; DM, dichroic mirror; SLM, spatial light modulator; PR, polarization rotator; L1–L5, lenses. (b) SIM imaging model and sCMOS camera noise model.
Photonics 09 00172 g001
Figure 2. Diagram of the simulation pipeline.
Figure 2. Diagram of the simulation pipeline.
Photonics 09 00172 g002
Figure 3. Performance improvement of the sCMOS noise-corrected SIM reconstruction algorithm. (a) The conventional algorithm reconstructs the SR image. (b) The noise-corrected algorithm reconstructs the SR image. (c) Zoomed-in images of selected subregions i and ii in (a,b). (d) The SNR of the SR image reconstructed by conventional and noise-corrected algorithms. (e) Temporal pixel fluctuation map (standard deviation in each pixel over time) over 80 SR images reconstructed by the conventional algorithm. (f) Temporal pixel fluctuation map over 80 SR images reconstructed by the noise-corrected algorithm. (g) Zoomed-in ROI of selected subregions i and ii in (e,f). (h) The distribution of pixel fluctuation standard deviation corresponds to conventional and noise-corrected algorithms.
Figure 3. Performance improvement of the sCMOS noise-corrected SIM reconstruction algorithm. (a) The conventional algorithm reconstructs the SR image. (b) The noise-corrected algorithm reconstructs the SR image. (c) Zoomed-in images of selected subregions i and ii in (a,b). (d) The SNR of the SR image reconstructed by conventional and noise-corrected algorithms. (e) Temporal pixel fluctuation map (standard deviation in each pixel over time) over 80 SR images reconstructed by the conventional algorithm. (f) Temporal pixel fluctuation map over 80 SR images reconstructed by the noise-corrected algorithm. (g) Zoomed-in ROI of selected subregions i and ii in (e,f). (h) The distribution of pixel fluctuation standard deviation corresponds to conventional and noise-corrected algorithms.
Photonics 09 00172 g003
Figure 4. Performance comparison of the noise-corrected SIM reconstruction algorithm with different regularized constraints. (a) Ground truth image. (b) The SR image reconstructed by the noise-corrected algorithm with a nonnegativity regularized constraint. (c) The SR image reconstructed by the noise-corrected algorithm with a high-frequency suppression regularized constraint. (d) The SR image reconstructed by the noise-corrected algorithm with the Tikhonov regularized constraint. (e) The horizontal profile averaged along the vertical direction of (a). (f) The horizontal profile averaged along the vertical direction of (b). (g) The horizontal profile averaged along the vertical direction of (c). (h) The horizontal profile averaged along the vertical direction of (d).
Figure 4. Performance comparison of the noise-corrected SIM reconstruction algorithm with different regularized constraints. (a) Ground truth image. (b) The SR image reconstructed by the noise-corrected algorithm with a nonnegativity regularized constraint. (c) The SR image reconstructed by the noise-corrected algorithm with a high-frequency suppression regularized constraint. (d) The SR image reconstructed by the noise-corrected algorithm with the Tikhonov regularized constraint. (e) The horizontal profile averaged along the vertical direction of (a). (f) The horizontal profile averaged along the vertical direction of (b). (g) The horizontal profile averaged along the vertical direction of (c). (h) The horizontal profile averaged along the vertical direction of (d).
Photonics 09 00172 g004
Table 1. SNR of the SR reconstructed image under various noise levels.
Table 1. SNR of the SR reconstructed image under various noise levels.
Shot Noise
Level
Readout Noise LevelSNR
ConventionalNoise-Corrected
1001011.16713.349
2010.33612.942
309.62712.521
1501012.28014.338
2011.43213.865
3010.85213.496
2001011.14913.363
208.73410.562
305.2946.460
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, B.; Huang, X.; Fan, J.; Chen, L. sCMOS Noise-Corrected Superresolution Reconstruction Algorithm for Structured Illumination Microscopy. Photonics 2022, 9, 172. https://doi.org/10.3390/photonics9030172

AMA Style

Zhou B, Huang X, Fan J, Chen L. sCMOS Noise-Corrected Superresolution Reconstruction Algorithm for Structured Illumination Microscopy. Photonics. 2022; 9(3):172. https://doi.org/10.3390/photonics9030172

Chicago/Turabian Style

Zhou, Bo, Xiaoshuai Huang, Junchao Fan, and Liangyi Chen. 2022. "sCMOS Noise-Corrected Superresolution Reconstruction Algorithm for Structured Illumination Microscopy" Photonics 9, no. 3: 172. https://doi.org/10.3390/photonics9030172

APA Style

Zhou, B., Huang, X., Fan, J., & Chen, L. (2022). sCMOS Noise-Corrected Superresolution Reconstruction Algorithm for Structured Illumination Microscopy. Photonics, 9(3), 172. https://doi.org/10.3390/photonics9030172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop