ACsN combina calibrazione della telecamera, stima del rumore e filtraggio sparso per correggere le sorgenti di rumore più rilevanti generate da una telecamera sCMOS (Fig. 1a e Note complementari 1 e 2.1). In particolare, ACsN corregge prima il rumore del pattern fisso utilizzando una mappa dell’offset e del guadagno dei pixel sCMOS. La presenza del rumore del modello fisso nelle telecamere sCMOS genera in pixel diversi (p) un numero diverso di fotoelettroni dallo stesso numero di fotoni che interferiscono (Sp). Questo effetto è proporzionale al livello di illuminazione e può essere modellato come un fattore moltiplicativo yp applicato al parametro della variabile Sp distribuita da Poisson. Allo stesso tempo, durante la conversione analogico-digitale (AD), la tensione prodotta da ciascun pixel viene letta come la differenza da un livello di riferimento, che rappresenta l’assenza di luce. In pratica, a questa tensione di riferimento viene assegnato un valore positivo che è responsabile di un bias (ßp) nei valori di intensità misurati. Pertanto, l’acquisizione di un sCMOS fotocamera può essere modellato mediante l’equazione:20
dove Zp è il valore del pixel p, t il tempo di esposizione, e N (0, σR) la Gaussiana distribuito rumore di lettura dei media µR = 0 e deviazione standard σR. Considerando la praticità della microscopia a fluorescenza, in questo modello abbiamo omesso il contributo della corrente oscura, che può essere ignorato per tempi di esposizione inferiori a 1 s, e il rumore di quantizzazione dovuto alla conversione AD, che è trascurabile rispetto al rumore di lettura3,21 (Nota integrativa 2.2).
Poiché il rumore del modello fisso dipende solo dalla circuiteria della fotocamera, ßp e yp possono essere stimati attraverso una calibrazione una tantum (vedi Metodi). Tuttavia, un’attenta valutazione sia del rumore di lettura distribuito da Gaussian, N(0, σR), sia della fluttuazione dovuta al rumore di fotone distribuito da Poisson, Poi{Sp (τ)}, è necessaria per ottenere una stima accurata del segnale sottostante Sp. Per eseguire questa valutazione, abbiamo ideato un modello di rumore che consente una stima congiunta della varianza del rumore analizzando la risposta in frequenza del sistema di microscopia. Questo si basa sul fatto che la distribuzione di Poisson del rumore del fotone può essere approssimata da una distribuzione gaussiana quando il flusso di fotoni è >3 fotoni per pixel22. In particolare, l’errore introdotto approssimando la varianza di Poisson, \(\sigma _P^2\), con una varianza gaussiana, \(\sigma _G^2\), diventa <1% quando il flusso di fotoni è superiore a 5 fotoni per pixel (Nota supplementare 2.3). In particolare, le condizioni sopra menzionate sul flusso di fotoni sono generalmente soddisfatte per molte applicazioni nella microscopia a fluorescenza23, 24. Pertanto, consideriamo il rumore correlato alla fotocamera come il risultato della somma di due variabili casuali distribuite gaussiane indipendenti, la cui varianza è \(\sigma _N^2 = \sigma _R^2 + \sigma _G^2\). Tale distribuzione consiste in una densità spettrale di potenza costante, mentre i segnali provenienti dal campione sono contenuti all’interno della funzione di trasferimento ottico (OTF)25. Pertanto, sfruttiamo la conoscenza del sistema ottico per valutare la fluttuazione dei pixel al di fuori dell’OTF, che è dovuta solo al rumore, e quindi usiamo il valore ottenuto per derivare σN nell’immagine originale (Nota integrativa 2.3).
Successivamente, l’algoritmo utilizza queste statistiche sul rumore per una valutazione non locale dell’auto-somiglianza del campione e per eseguire filtri sparsi collaborativi sulla sequenza di input. A differenza delle precedenti implementazioni del filtro collaborativo, abbiamo adottato un approccio a strati che sonda sequenzialmente l’auto-somiglianza dell’immagine nello spazio e nel tempo al fine di migliorare la correzione del rumore senza sacrificare la precisione e il runtime. In breve, il filtro scompone l’immagine in patch e le ordina in gruppi tridimensionali (3D) in base alla loro similità26. Quindi, impiega una trasformazione 3D per elaborare ogni gruppo tutto in una volta. Il denoising viene eseguito mediante hard-thresholding e migliorato dal fatto che, a causa della somiglianza tra le patch, la trasformazione 3D risulta in una rappresentazione ancora più sparsa delle patch originali, mentre lo spettro di potenza del rumore rimane costante27. Successivamente, le patch denoised vengono restituite alle loro posizioni originali per formare un’immagine intermedia. A questo punto, il filtro collaborativo viene eseguito una seconda volta, ma sostituendo l’hard-thresholding con un filtro Wiener. Il filtro viene eseguito utilizzando sia le immagini rumorose che quelle intermedie e genera l’immagine denoised finale (Nota complementare 2.4). Va notato che la variazione spaziale del rumore attraverso l’immagine può influenzare le prestazioni del filtro Wiener. Tuttavia, ciò è notevolmente mitigato dall’uso dell’elaborazione basata su patch, che, rispetto all’intera immagine, migliora l’uniformità dell’intensità all’interno dei singoli gruppi di patch, mostrando una grande stabilità contro il rumore spaziale9.
Infine, viene eseguito un altro filtro collaborativo alla ricerca di patch simili anche nei frame vicini. In questo modo, il rumore persistente può essere ulteriormente ridotto sfruttando l’auto-somiglianza del campione nel tempo preservando la risoluzione temporale18 (Nota complementare 2.5).
Caratterizzazione di ACsN
Successivamente, abbiamo caratterizzato le prestazioni di ACsN utilizzando sia dati numerici che sperimentali. In particolare, il filtraggio collaborativo ACsN dipende dalla stima di σN, nonché dalla scelta dei parametri nell’algoritmo28, che sono stati scelti per ottimizzare sia la correzione del rumore che il runtime (Nota integrativa 3.1). Abbiamo osservato che la nostra strategia può attenuare significativamente l’effetto dannoso del rumore della fotocamera, evitando la perdita di risoluzione dell’immagine, specialmente in presenza di rumore altamente spazialmente variante (Nota integrativa 3.2). Inoltre, il rumore della fotocamera può indurre fluttuazioni temporali dei valori dei pixel che non sono correlati al campione, influenzando così l’analisi quantitativa dei dati time-lapse. Il denoising ACsN riduce questo effetto di circa un ordine di grandezza, con fluttuazioni residue paragonabili a quelle di una fotocamera ideale (Fig. 1b-g e Nota complementare 3.3). Inoltre, va notato che a basso numero di fotoni, i dettagli del campione iniziano ad essere comparabili con le fluttuazioni del rumore e diventano più difficili da recuperare. Pertanto, le prestazioni del ripristino dell’immagine sono intrinsecamente correlate al flusso fotonico dell’immagine in ingresso. Tuttavia, utilizzando sia simulazioni che dati sperimentali, abbiamo verificato una robusta correzione del rumore ACsN a livelli di scarsa illuminazione fino a 5-10 fotoni per pixel (Nota supplementare 3.4).
Inoltre, abbiamo convalidato le prestazioni di ACsN sotto varie frequenze di campionamento normalmente adottate per la microscopia a fluorescenza. In pratica, una frequenza di campionamento vicina al criterio di Nyquist rappresenta un buon compromesso tra rapporto segnale-rumore (SNR) e conservazione dei dettagli. Qui, esaminando numericamente e sperimentalmente attraverso una vasta gamma di frequenze di campionamento, abbiamo dimostrato la fattibilità di ACsN per SNR basso con sovracampionamento e nessuna perdita evidente di segnali con sottocampionamento (Nota supplementare 3.5).
A differenza delle immagini naturali, le immagini fluorescenti di campioni biologici sono altamente specificate, esibendo bersagli molecolari o strutture etichettate con precisione nelle cellule. Di conseguenza, ogni immagine fluorescente caratterizza solitamente gli oggetti specifici che ricorrono attraverso il campo di vista, che fornisce l’auto-somiglianza non locale sufficiente per rendere l’algoritmo notevolmente efficiente per la microscopia di fluorescenza. Con dati numerici e sperimentali, abbiamo caratterizzato la dipendenza delle prestazioni ACsN dall’uso dell’auto-somiglianza di un’immagine di input (Nota supplementare 3.6). Inoltre, come mostrato di seguito, abbiamo valutato quantitativamente una varietà di campioni non biologici e biologici per verificare la fattibilità del metodo, che copre varie dimensionalità, morfologia, casualità e densità, come bersagli di calibro, particelle fluorescenti, singole molecole, microtubuli, filamenti di actina, mitocondri, filopodi, lamellipodi e piccoli animali.
Microscopia a ampio campo
La microscopia a ampio campo, in particolare la microscopia a fluorescenza a riflessione interna totale (TIRF), è una delle tecniche più utilizzate nell’immagine cellulare29. TIRF utilizza il fenomeno della riflessione interna totale della luce all’interfaccia vetro / acqua per creare un’onda evanescente che si propaga solo per poche centinaia di nanometri attraverso il vetrino coprioggetti. Ciò consente l’eccitazione selettiva delle etichette fluorescenti nella parte inferiore del campione (Fig. 1 bis). Tuttavia, in caso di emettitori fluorescenti deboli, bassa intensità luminosa o breve tempo di esposizione, il rumore correlato a sCMOS diventa grave e peggiora la qualità dell’immagine (Fig. 1 ter). ACSN denoising può ridurre efficacemente tale contributo e recuperare i segnali non distorti dal rumore, consentendo un’acquisizione più rapida senza compromettere il segnale sottostante (Fig. 1 quater, d).
Abbiamo dimostrato ACSN denoising di microscopia a ampio campo in entrambe le configurazioni di epi-fluorescenza e TIRF utilizzando vari campioni sub-cellulari fissi, vivi e multicolori, compresi i microtubuli (Fig. 1 e Fig. 1), mitocondri (Fig. 2 e Film supplementari 1 e 2), e F-actina (Fig. 2). L’uso di ACsN può mantenere la stessa qualità dell’immagine con un tempo di esposizione più breve (cioè, una migliore risoluzione temporale) e un livello di eccitazione più basso (cioè, meno danni alla foto). Le prestazioni sono, quindi, limitate principalmente dalla foto-fisica degli emettitori fluorescenti. Utilizzando metriche quantitative, abbiamo dimostrato che il metodo può recuperare immagini a campo largo con un budget di fotoni di due ordini di grandezza inferiore senza perdita di qualità dell’immagine (Tabella supplementare 1).
Leave a Reply