domenica 24 febbraio 2008

Come calcolare il gain e il readout noise

Ora che conosciamo il significato di gain e readout noise per una camera CCD, vediamo come possiamo determinarne il valore nel caso queste quantità non fossero riportate nella documentazione che accompagna il nostro dispositivo.

Esiste un metodo molto semplice e rapido da effettuare, anche se non molto preciso, che ci darà comunque una stima affidabile di queste due quantità.

Una volta installata la camera CCD sul telescopio, è sufficiente riprendere due FLAT FIELD, [F1] e [F2] e due BIAS FRAME [B1] e [B2] .(Nota: utilizzeremo sempre questa convenzione per indicare un "frame" della camera CCD, ovvero quella di rinchiudere l'identificativo tra parentesi quadre: si tratta infatti di una matrice di numeri e non di un singolo valore numerico; questo ci sarà utile per una più agevole lettura delle formule matematiche utilizzate di seguito).

L'unica importante avvertenza riguarda i flat field: questi non potranno essere dei sky-flat in quanto ben difficilmente i due flat field avranno lo stesso identico illuminamento poichè saranno stati ripresi in tempi diversi e quindi con una luminosità del cielo differente. E' necessario riprendere due flat field con la tecnica del dome-flat ovvero utilizzando una sorgente luminosa stabile e costante.

Una volta ripresi e salvati su disco i due flat field e bias frame, utilizzando Astroart possiamo aprire le finestre della statistica per ogni frame ed annotarne il valore medio: indicheremo questi valori ripettivamente con avg[F1], avg[F2], avg[B1] e avg[B2].

A questo punto sottraiamo un flat field con l'altro, cioè utilizzando il comando Matematica > Sottrai di Astroart eseguiamo l'operazione [F1]-[F2]. Otterremo un frame dall'aspetto completamente nero come se tutti i pixel avessero un valore nullo: in realtà sono le soglie di visualizzazione dell'immagine che debbono essere regolate e per farlo velocemente e in modo automatico è sufficiente fare click con il mouse sulla barra di stato dell'immagine. Dovreste ottenere una rappresentazione della "mappa del rumore" presente nel flat field come nell'immagine qui sotto.

Tuttavia non è tanto la rappresentazione grafica che ci interessa quanto il suo valore numerico: torniamo ad aprire la finestra delle statistiche ed annotiamoci il valore della deviazione standard: la indicheremo con ds([F1]-[F2]). Ripetendo gli stessi passaggi per i due bias frames otterremo il valore ds([B1]-[B2]).

Ora abbiamo tutte le informazioni necessarie per applicare la seguente formula per il calcolo del gain:

mentre applicheremo la seguente per il calcolo del readout noise:

Ricordiamo sempre che il gain è espresso in elettroni/ADU mentre il readout noise è espresso in elettroni.