ACsN łączy kalibrację kamery, szacowanie szumów i rzadkie filtrowanie, aby skorygować najbardziej istotne źródła szumów generowane przez kamerę sCMOS (rys. 1a oraz uwagi uzupełniające 1 i 2.1). W szczególności, ACsN najpierw koryguje szum o stałym wzorze za pomocą mapy przesunięcia i wzmocnienia pikseli sCMOS. Obecność szumu o stałym wzorze w aparatach sCMOS generuje w różnych pikselach (p) inną liczbę fotoelektronów z tej samej liczby uderzających fotonów (Sp). Efekt ten jest proporcjonalny do poziomu oświetlenia i może być modelowany jako współczynnik mnożnikowy YP zastosowany do parametru zmiennej rozproszonej Poissona Sp. Jednocześnie, podczas konwersji analogowo-cyfrowej (AD), napięcie wytwarzane przez każdy piksel jest odczytywane jako różnica w stosunku do poziomu odniesienia, co oznacza brak światła. W praktyce temu napięciu odniesienia przypisuje się wartość dodatnią, która odpowiada za odchylenie (ßp) w zmierzonych wartościach natężenia. Zatem pozyskanie kamery sCMOS można modelować za pomocą równania:20
gdzie ZP jest wartością piksela p, τ czasu ekspozycji i N (0, σr) rozkładu Gaussa szumu odczytu średniego µr = 0 i odchylenia standardowego σr. Biorąc pod uwagę praktyczność mikroskopii fluorescencyjnej, w tym modelu pominęliśmy udział ciemnego prądu, który może być pominięty dla czasów ekspozycji poniżej 1 s, oraz szum kwantyzacji spowodowany konwersją AD, który jest znikomy w porównaniu do szumu odczytu3,21 (Uwaga uzupełniająca 2.2).
ponieważ szum o stałym wzorze zależy tylko od obwodu kamery, ßp i yp można oszacować za pomocą jednorazowej kalibracji (patrz metody). Jednak dokładna ocena zarówno szumu odczytu rozproszonego Gaussa, N(0, σR), jak i fluktuacji spowodowanej szumem fotonu rozproszonego Poissona, Pois{SP(τ)}, jest konieczna do uzyskania dokładnego oszacowania sygnału wyjściowego Sp. W celu przeprowadzenia tej oceny opracowaliśmy model szumu, który pozwala na wspólne oszacowanie wariancji szumu poprzez analizę odpowiedzi częstotliwościowej układu mikroskopowego. Jest to oparte na fakcie, że rozkład Poissona szumu fotonu może być prawdopodobnie przybliżony przez rozkład Gaussa, gdy strumień fotonu wynosi >3 fotony na piksel22. W szczególności, błąd wprowadzony przez przybliżenie wariancji Poissona, \(\sigma _p^2\), z wariancją Gaussa, \(\sigma _g^2\), staje się <1%, gdy strumień fotonów jest większy niż 5 fotonów na piksel (Uwaga uzupełniająca 2.3). W szczególności, wyżej wymienione warunki dotyczące strumienia fotonów są zwykle spełnione dla wielu zastosowań w mikroskopie fluorescencyjnym23, 24. Dlatego uważamy, że szum związany z kamerą jest wynikiem sumy dwóch niezależnych zmiennych losowych rozproszonych Gaussa, których wariancja wynosi \(\sigma _n^2 = \ sigma _r^2 + \ sigma _g^2\). Taki rozkład składa się ze stałej gęstości widmowej mocy, podczas gdy sygnały pochodzące z próbki są zawarte w funkcji transferu optycznego (OTF)25. Dlatego wykorzystujemy znajomość układu optycznego do oceny fluktuacji pikseli poza OTF, która jest spowodowana tylko szumem, a następnie wykorzystujemy uzyskaną wartość do uzyskania σN na oryginalnym obrazie (Uwaga uzupełniająca 2.3).
następnie algorytm wykorzystuje te statystyki szumu do nielokalnej oceny samopodobieństwa próbki i do przeprowadzenia wspólnego filtrowania na sekwencji wejściowej. W przeciwieństwie do poprzednich implementacji filtrowania opartego na współpracy, przyjęliśmy podejście warstwowe, które sekwencyjnie bada samopodobieństwo obrazu w przestrzeni i czasie w celu poprawy korekcji szumów bez poświęcania dokładności i czasu pracy. Krótko mówiąc, filtr rozkłada obraz w plastry i sortuje go w trójwymiarowe (3D) grupy zgodnie z ich podobieństwo26. Następnie wykorzystuje transformację 3D, aby przetworzyć każdą grupę naraz. Denoising jest wykonywany przez twarde progowanie i wzmocniony przez fakt, że ze względu na podobieństwo między łat, transformacja 3D powoduje jeszcze rzadszą reprezentację oryginalnych łat, podczas gdy widmo mocy szumu pozostaje stałe27. Następnie odbarwione łaty są zwracane do swoich pierwotnych lokalizacji, tworząc pośredni obraz. W tym momencie filtr współpracujący jest uruchamiany po raz drugi, ale zastępuje twardy próg filtrem Wienera. Filtr jest wykonywany przy użyciu zarówno głośnych, jak i pośrednich obrazów i generuje końcowy odbarwiony obraz (Uwaga uzupełniająca 2.4). Należy zauważyć, że przestrzenna zmienność szumu na obrazie może wpływać na działanie filtra Wienera. Jest to jednak znacznie złagodzone przez zastosowanie przetwarzania opartego na łatach, które w porównaniu z całym obrazem zwiększa jednolitość intensywności w poszczególnych grupach łatek, wykazując dużą stabilność w stosunku do przestrzennego wariantu szumu9.
wreszcie wykonywany jest kolejny filtr współpracujący szukający podobnych łat również w sąsiednich klatkach. W ten sposób można jeszcze bardziej ograniczyć utrzymujący się hałas, wykorzystując samopodobieństwo próbki w czasie, przy zachowaniu rozdzielczości czasowej18 (Uwaga uzupełniająca 2.5).
charakterystyka ACsN
następnie scharakteryzowaliśmy wydajność ACsN zarówno na podstawie danych numerycznych, jak i eksperymentalnych. W szczególności, wspólne filtrowanie ACsN zależy od estymacji σN, jak również od wyboru parametrów w algorytmie28, które zostały wybrane w celu optymalizacji zarówno korekcji szumu, jak i czasu pracy (Uwaga uzupełniająca 3.1). Zauważyliśmy, że nasza strategia może znacznie złagodzić szkodliwy wpływ szumu kamery, unikając utraty rozdzielczości obrazu, zwłaszcza w przypadku zakłóceń o dużym wariancie przestrzennym (Uwaga uzupełniająca 3.2). Co więcej, szum kamery może wywoływać czasowe wahania wartości pikseli, które nie są związane z próbką, wpływając w ten sposób na ilościową analizę danych poklatkowych. Odazotowanie ACsN zmniejsza ten efekt o około jeden rząd wielkości, z resztkowymi fluktuacjami porównywalnymi do idealnego aparatu fotograficznego (rys. 1b-g i Uwaga uzupełniająca 3.3). Ponadto należy zauważyć, że przy niskiej liczbie fotonów szczegóły próbki zaczynają być porównywalne z fluktuacjami szumów i stają się trudniejsze do odzyskania. W ten sposób wydajność przywracania obrazu jest nieodłącznie związana ze strumieniem fotonów obrazu wejściowego. Niemniej jednak, wykorzystując zarówno symulacje, jak i dane eksperymentalne, zweryfikowaliśmy solidną korekcję szumów ACsN przy słabym oświetleniu do 5-10 fotonów na piksel (Uwaga uzupełniająca 3.4).
ponadto zwalidowaliśmy wydajność ACsN przy różnych częstotliwościach próbkowania normalnie przyjętych do mikroskopii fluorescencyjnej. W praktyce częstotliwość próbkowania zbliżona do kryterium Nyquista stanowi dobrą kompromis między stosunkiem sygnału do szumu (SNR) a zachowaniem szczegółów. Tutaj, badając numerycznie i doświadczalnie w szerokim zakresie częstotliwości próbkowania, wykazaliśmy żywotność ACsN dla niskiego SNR z oversamplingiem i bez zauważalnej utraty sygnałów z niedostatecznym próbkowaniem (Uwaga uzupełniająca 3.5).
w przeciwieństwie do naturalnych obrazów, fluorescencyjne obrazy próbek biologicznych są wysoce określone, wykazując precyzyjnie oznaczone cele molekularne lub struktury w komórkach. Dlatego każdy obraz fluorescencyjny zwykle zawiera określone obiekty powtarzające się w polu widzenia, które zapewnia wystarczające nielokalne samopodobieństwo, aby algorytm był szczególnie skuteczny w mikroskopii fluorescencyjnej. Na podstawie danych numerycznych i eksperymentalnych scharakteryzowaliśmy zależność działania ACsN od zastosowania samopodobieństwa obrazu wejściowego (Uwaga uzupełniająca 3.6). Ponadto, jak pokazano poniżej, oceniliśmy ilościowo wiele niebiologicznych i biologicznych próbek w celu weryfikacji żywotności metody, obejmujących różne wymiarowość, morfologię, losowość i gęstość, takie jak cele kalibru, cząsteczki fluorescencyjne, pojedyncze cząsteczki, mikrotubule, włókna aktyny, mitochondria, filopodia, lamellipodia i małe zwierzęta.
mikroskopia szerokokątna
mikroskopia szerokokątna, w szczególności mikroskopia TIRF (Total internal reflection fluorescence), jest jedną z najczęściej stosowanych technik w obrazowaniu komórkowym29. TIRF wykorzystuje zjawisko całkowitego wewnętrznego odbicia światła na styku szkło-woda w celu stworzenia fali ulotnej, która rozchodzi się tylko przez kilkaset nanometrów w poprzek pokrywy. Umożliwia to selektywne wzbudzenie etykiet fluorescencyjnych na dole próbki(dodatkowe rys. 1a). Jednak w przypadku słabych emiterów fluorescencyjnych, słabego natężenia światła lub krótkiego czasu ekspozycji szum związany z sCMOS staje się poważny i pogarsza jakość obrazu (dodatkowe rys. 1B). Odszumianie ACsN może skutecznie zmniejszyć taki wkład i odzyskać niezakłócone sygnały z szumu, umożliwiając szybszą akwizycję bez narażania sygnału bazowego(dodatkowe rys. 1c, d).
zademonstrowaliśmy denoizację mikroskopii szerokokątnej ACsN zarówno w konfiguracjach epi-fluorescencji, jak i TIRF przy użyciu różnych stałych, żywych i wielokolorowych próbek subkomórkowych, w tym mikrotubul (Fig. 1 i dodatkowe rys. 1), mitochondria (rys. 2 i dodatkowe filmy 1 i 2) oraz F-aktyna (rys. 2). Zastosowanie ACsN może utrzymać tę samą jakość obrazu przy krótszym czasie ekspozycji (tj. lepszej rozdzielczości czasowej) i niższym poziomie wzbudzenia (tj. mniej uszkodzeń zdjęć). Wydajność jest więc ograniczona przede wszystkim przez fotofizykę emiterów fluorescencyjnych. Korzystając z metryk ilościowych, pokazaliśmy, że metoda może odzyskać obrazy o szerokim polu przy budżecie fotonu o dwa rzędy wielkości mniejszej bez utraty jakości obrazu (tabela uzupełniająca 1).
Leave a Reply