ACsN combina calibración de cámara, estimación de ruido y filtrado disperso para corregir las fuentes de ruido más relevantes generadas por una cámara sCMOS (Fig. 1a y Notas complementarias 1 y 2.1). En particular, ACsN corrige primero el ruido de patrón fijo utilizando un mapa del desplazamiento y la ganancia de los píxeles de sCMOS. La presencia del ruido de patrón fijo en las cámaras sCMOS genera en diferentes píxeles (p) un número diferente de fotoelectrones del mismo número de fotones de impacto (Sp). Este efecto es proporcional al nivel de iluminación y se puede modelar como un factor multiplicativo yp aplicado al parámetro de la variable distribuida por Poisson Sp. Al mismo tiempo, durante la conversión analógica a digital (AD), el voltaje producido por cada píxel se lee como la diferencia con respecto a un nivel de referencia, que representa la ausencia de luz. En la práctica, a esta tensión de referencia se le asigna un valor positivo que es responsable de un sesgo (ßp) en los valores de intensidad medidos. Por lo tanto, la adquisición de una cámara sCMOS se puede modelar mediante la ecuación:20
donde Zp es el valor del píxel p, τ el tiempo de exposición, y N (0, σR) el ruido de lectura distribuido gaussiano de µR medio = 0 y desviación estándar σR. Teniendo en cuenta la practicidad de la microscopía de fluorescencia, en este modelo hemos omitido la contribución de la corriente oscura, que puede descartarse para tiempos de exposición inferiores a 1 s, y el ruido de cuantización debido a la conversión de ANUNCIO, que es insignificante en comparación con el ruido de lectura 3,21 (Nota complementaria 2.2).
Dado que el ruido de patrón fijo depende solo de los circuitos de la cámara, ßp y yp se pueden estimar a través de una calibración única (consulte Métodos). Sin embargo, una evaluación cuidadosa tanto del ruido de lectura distribuido gaussiano, N(0, σR), como de la fluctuación debida al ruido de disparo de fotones distribuido por Poisson, Pois{Sp(τ)}, es necesaria para obtener una estimación precisa de la señal subyacente Sp. Para realizar esta evaluación, diseñamos un modelo de ruido que permite una estimación conjunta de la varianza de ruido mediante el análisis de la respuesta de frecuencia del sistema de microscopía. Esto se basa en el hecho de que la distribución de Poisson del ruido de disparo de fotones se puede aproximar de manera factible mediante una distribución gaussiana cuando el flujo de fotones es >3 fotones por pixel22. En particular, el error introducido al aproximar la varianza de Poisson, \(\sigma _P^2\), con una varianza gaussiana, \(\sigma _G^2\), se convierte en <1% cuando el flujo de fotones es más de 5 fotones por píxel (Nota complementaria 2.3). En particular, las condiciones mencionadas anteriormente en el flujo de fotones generalmente se satisfacen para muchas aplicaciones en microscopía de fluorecencia23, 24. Por lo tanto, consideramos el ruido relacionado con la cámara como el resultado de la suma de dos variables aleatorias distribuidas por Gauss independientes, cuya varianza es \(\sigma _N^2 = \sigma _R^2 + \sigma _G^2\). Dicha distribución consiste en una densidad espectral de potencia constante, mientras que las señales procedentes de la muestra están contenidas en la función de transferencia óptica (OTF) 25. Por lo tanto, aprovechamos el conocimiento del sistema óptico para evaluar la fluctuación de píxeles fuera del OTF, que se debe solo al ruido, y luego usamos el valor obtenido para derivar σN en la imagen original (Nota complementaria 2.3).
A continuación, el algoritmo utiliza estas estadísticas de ruido para una evaluación no local de la auto-similitud de la muestra y para realizar un filtrado colaborativo disperso en la secuencia de entrada. A diferencia de implementaciones anteriores de filtrado colaborativo, adoptamos un enfoque por capas que sondea secuencialmente la auto-similitud de la imagen en el espacio y el tiempo para mejorar la corrección de ruido sin sacrificar la precisión y el tiempo de ejecución. En resumen, el filtro descompone la imagen en parches y los clasifica en grupos tridimensionales (3D) de acuerdo con su similaridad26. Luego, emplea una transformación 3D para procesar todos los grupos a la vez. La eliminación de ruido se realiza mediante umbral duro y se ve reforzada por el hecho de que, debido a la similitud entre los parches, la transformación 3D da como resultado una representación aún más escasa de los parches originales, mientras que el espectro de potencia de ruido permanece constante27. Posteriormente, los parches desprovistos se devuelven a sus ubicaciones originales para formar una imagen intermedia. En este punto, el filtro colaborativo se ejecuta por segunda vez, pero reemplazando el umbral mínimo por un filtro Wiener. El filtro se realiza utilizando tanto las imágenes ruidosas como las intermedias y genera la imagen desensibilada final (Nota complementaria 2.4). Cabe señalar que la variación espacial del ruido a lo largo de la imagen puede afectar al rendimiento del filtro Wiener. Sin embargo, esto se ve considerablemente mitigado por el uso de procesamiento basado en parches, que, en comparación con la imagen completa, mejora la uniformidad de intensidad dentro de los grupos de parches individuales, exhibiendo una gran estabilidad contra el ruido de variantes espaciales9.
Finalmente, se realiza otro filtro colaborativo buscando parches similares también en los marcos vecinos. De esta manera, el ruido persistente se puede reducir aún más aprovechando la auto-similitud de la muestra en el tiempo, preservando al mismo tiempo la resolución temporal18 (Nota complementaria 2.5).
Caracterización de ACsN
A continuación, caracterizamos el desempeño de ACsN utilizando datos numéricos y experimentales. En particular, el filtrado colaborativo de ACsN depende de la estimación de σN, así como de la elección de los parámetros en el algorítmo28, que se eligieron para optimizar tanto la corrección de ruido como el tiempo de ejecución (Nota complementaria 3.1). Observamos que nuestra estrategia puede atenuar significativamente el efecto perjudicial del ruido de la cámara, evitando la pérdida de resolución de la imagen, especialmente en presencia de ruido con variantes espaciales elevadas (Nota complementaria 3.2). Además, el ruido de la cámara puede inducir fluctuaciones temporales de los valores de píxeles que no están relacionados con la muestra, lo que afecta al análisis cuantitativo de los datos de lapso de tiempo. La eliminación de ruido ACsN reduce este efecto en aproximadamente un orden de magnitud, con fluctuaciones residuales comparables a las de una cámara ideal (Fig. 1b a g y Nota complementaria 3.3). Además, debe tenerse en cuenta que con recuentos bajos de fotones, los detalles de la muestra comienzan a ser comparables con las fluctuaciones de ruido y se vuelven más difíciles de recuperar. Por lo tanto, el rendimiento de la restauración de imágenes está intrínsecamente relacionado con el flujo de fotones de la imagen de entrada. Sin embargo, utilizando simulaciones y datos experimentales, verificamos una corrección de ruido ACsN robusta a niveles de poca luz de hasta 5-10 fotones por píxel (Nota complementaria 3.4).
Además, validamos el rendimiento de ACsN bajo varias tasas de muestreo normalmente adoptadas para microscopía de fluorescencia. En la práctica, una frecuencia de muestreo cercana al criterio de Nyquist representa un buen equilibrio entre la relación señal / ruido (SNR) y la preservación de los detalles. Aquí, examinando numérica y experimentalmente a través de una amplia gama de frecuencias de muestreo, demostramos la viabilidad de ACsN para SNR bajo con sobremuestreo y sin pérdida notable de señales con submuestreo (Nota complementaria 3.5).
A diferencia de las imágenes naturales, las imágenes fluorescentes de muestras biológicas son altamente especificadas, exhibiendo objetivos moleculares o estructuras en células etiquetadas con precisión. Por lo tanto, cada imagen fluorescente generalmente presenta objetos específicos recurrentes en todo el campo de visión, lo que proporciona suficiente auto-similitud no local para hacer que el algoritmo sea notablemente eficiente para la microscopía de fluorescencia. Con datos numéricos y experimentales, caracterizamos la dependencia del desempeño de ACsN en el uso de la auto-similitud de una imagen de entrada (Nota Complementaria 3.6). Además, como se muestra a continuación, evaluamos cuantitativamente una variedad de muestras no biológicas y biológicas para verificar la viabilidad del método, abarcando diversas dimensiones, morfología, aleatoriedad y densidad, como blancos de calibre, partículas fluorescentes, moléculas individuales, microtúbulos, filamentos de actina, mitocondrias, filopodios, lamelipodios y animales pequeños.
Microscopía de campo amplio
La microscopía de campo amplio, especialmente la microscopía de fluorescencia de reflexión interna total (TIRF, por sus siglas en inglés), es una de las técnicas más utilizadas en la creación de imágenes de cellas29. TIRF utiliza el fenómeno de la reflexión interna total de la luz en la interfaz vidrio/agua para crear una onda evanescente que se propaga solo por unos pocos cientos de nanómetros a través del cubreobjetos. Esto permite la excitación selectiva de las etiquetas fluorescentes en la parte inferior de la muestra (Fig. 1a). Sin embargo, en el caso de emisores fluorescentes débiles, poca intensidad de luz o un tiempo de exposición corto, el ruido relacionado con sCMOS se vuelve severo y deteriora la calidad de la imagen (Fig. 1b). La eliminación de ruido de ACsN puede reducir eficazmente dicha contribución y recuperar las señales no distorsionadas del ruido, lo que permite una adquisición más rápida sin comprometer la señal subyacente (Fig. 1c, d).
Se demostró la eliminación de ACsN de microscopía de campo amplio en configuraciones de epi-fluorescencia y TIRF utilizando varias muestras subcelulares fijas, vivas y multicolores, incluidos microtúbulos (Fig. 1 y Suplemento Fig. 1), mitocondrias (Fig. 2 y Películas suplementarias 1 y 2), y actina F (Fig. 2). El uso de ACsN puede mantener la misma calidad de imagen con un tiempo de exposición más corto (es decir, una mejor resolución temporal) y un nivel de excitación más bajo (es decir, menos daños fotográficos). El rendimiento está, por lo tanto, limitado principalmente por la fotofísica de los emisores fluorescentes. Usando métricas cuantitativas, mostramos que el método puede recuperar imágenes de campo amplio con un presupuesto de fotones dos órdenes de magnitud más bajo sin pérdida de calidad de imagen (Tabla complementaria 1).
Leave a Reply