Calendario 2017

E mentre tutto sembra scomporsi: febbrili le convulsioni della ragione, che giace sbronza in un fosso; mute le urla che lacerano il silenzio; feriti gli occhi dai frammenti del tempo spaccato; straziato lo spazio, costipato in un punto; ed esangue il fuoco perpetuo dell’uguaglianza,… mentre succede tutto questo, zitta zitta, senza dar nell’occhio, la Terra la sua, di rivoluzione, ha già conclusa. Senza sostare di fronte ad alcuna delle biforcazioni del labirinto. Beffarda e inesorabile, scrivendo le sue balte nelle ombre dei monti.

Scarica il Calendario 2017.

Permalink link a questo articolo: http://www.giovanniceribella.eu/fuere/?p=1514

Transito di Mercurio (video)

Sebbene siano passati tre mesi e un paio di giorni, pubblico solo ora una piccola testimonianza del transito di Mercurio sul disco del Sole.

Il video ritrae, in maniera accelerata, il passaggio di Mercurio vicino ad un piccolo gruppo di macchie solari, attorno alle ore 16.10′ del 9 maggio scorso. La ripresa è allineata sul piccolo disco del pianeta, come se una telecamera lo seguisse; perciò si vede la superficie del Sole scorrervi sotto.

I fotogrammi di questo video, come le mie osservazioni, sono stati fatti da un piccolo comune in Baviera, Flintsbach am Inn. Il giorno del transito il tempo in Italia era così brutto che sono dovuto fuggire al di là delle alpi per poter vedere il fenomeno! Il cielo era completamente coperto fino al Brennero, mentre sopra Innsbruck si apriva qualche spiraglio. Alla fine, mio padre ed io ci siamo fermati poco oltre il confine tedesco, appena mezz’ora prima dell’inizio del transito. Siamo riusciti ad osservare tutta la sequenza d’ingresso e il lento passaggio fino al momento di massimo avvicinamento al centro del disco solare, quando anche lì sono sopraggiunte le nuvole (ore 16.50′ circa).

Come ho scritto nell’articolo precedente, c’è stato bisogno di scrivere un programma apposta per riallineare i fermi-immagine che compongono il breve video. Questo è anche il motivo per cui sin’ora non sono riuscito a trarre qualche altra bella fotografia dalle decine di pose che ho scattato quel giorno: i programmi che uso di solito non funzionano perché le immagini solari sono povere di dettagli. Pian piano, quando troverò il tempo per passare un paio di giornate al computer, dovrei poter aggiungere anche qualche altra ripresa.

Il transito di Mercurio è più frequente di quello di Venere. Mentre quest’ultimo avviene a coppie ogni 115 anni, il pianeta più piccolo del Sistema Solare passerà ancora davanti al Sole l’11 novembre 2019 e in altre dieci occasioni nel corso di questo secolo. Quindi non disperate! Se non l’avete visto, o se le nuvole ve l’hanno nascosto, ci sono ancora delle altre possibilità…

A presto!

Permalink link a questo articolo: http://www.giovanniceribella.eu/fuere/?p=1492

Allineamento immagini solari

Il 9 maggio scorso il tempo sull’Italia settentrionale era pessimo: era in transito una perturbazione di provenienza atlantica e sulle prealpi la copertura nuvolosa era pressoché totale. Per riuscire a vedere il transito di Mercurio, un fenomeno relativamente raro che consiste nel passaggio del piccolo pianeta davanti al Sole, ho fatto una breve gita in Germania, valicando le Alpi e sorpassando le nuvole in volata. Appena finirò di elaborare le immagini pubblicherò un resoconto dell’escursione.

Tornato a casa, si è posto il problema di come utilizzare il materiale raccolto. Mi sono concentrato su un breve video, della durata di 8 minuti circa, che riprendeva Mercurio mentre passava vicino a un gruppo di macchie solari. Il filmato era stato girato alla velocità di 50 fotogrammi al secondo, ma ciascuno di essi era disturbato dalla turbolenza atmosferica, che impediva di scorgere i dettagli dell’inquadratura. In casi come questo di solito si mediano più fotogrammi per ottenerne uno solo, in cui gli effetti di distorsione sono ridotti. Per fare ciò bisogna che tutti i fotogrammi siano prima riallineati nel modo più preciso possibile.

L’allineamento delle fotografie astronomiche che riprendono campi stellari è relativamente semplice: bisogna far corrispondere le immagini delle stelle in un’immagine alle immagini delle stesse stelle di un’altra. Con delle foto della superficie solare ad alto ingrandimento, la scarsità e la variabilità dei dettagli rendono l’operazione più difficile. Cercando in rete ho trovato un algoritmo generico abbastanza semplice che può essere d’aiuto. Funziona così: di ciascuna immagine si calcola la trasformata di Fourier bidimensionale discreta \mathcal{I}_i(\vec{k}\,) e si costruisce lo spettro di potenza incrociato:

{\rm P}_{i,j}(\vec{k}\,)\,=\,\mathcal{I}_i(\vec{k}\,)\cdot\mathcal{I}_j^{\mbox{*}}(\vec{k}\,)

Qui \vec{k}\,=\,(k_x,\,k_y) è il momento associato a \vec{r}\,=\,(x,\,y) che descrive le immagini bidimensionali.
Così come lo spettro di potenza di un segnale è la trasformata della funzione di autocorrelazione, si può mostrare che lo spettro di potenza incrociato è la trasformata della funzione di correlazione incrociata (o cross-correlazione):

{\rm R}_{i,j}(\vec{r}\,)\,=\,\mathcal{F}^{-1}\big[ {\rm P}_{i,j}(\vec{k}\,) \big]

Quando uno ha per le mani la funzione di correlazione incrociata “er ist fertig”, almeno in teoria. Infatti, in linea di principio, se le due immagini fossero l’una la traslazione rigida dell’altra la cross-correlazione avrebbe un picco di intensità alle coordinate corrispondenti al vettore di traslazione. In altre parole, quando si ha la cross-correlazione si sa di quanto spostare la seconda immagine per allinearla alla prima.

Tradurre in una tecnica pratica i risultati di qualche calcolo non è mai un’operazione immediata. Per allineare i fotogrammi ho scritto uno script in Python (potete scaricarlo al termine dell’articolo) che include questi passaggi:
 

  1. Le immagini sono ridotte in bianco e nero e la scala dei colori è invertita, in modo che il dettaglio più evidente (il disco di Mercurio) risulti molto luminoso su uno sfondo scuro. Le due immagini da allineare sono ricampionate alla metà della loro dimensione originaria (bin 2×2), per ridurre gli effetti della turbolenza.
  2. A ciascuna viene applicata una maschera gaussiana con 4px di sigma per distruggere eventuali difetti puntuali dell’immagine. Si trova il massimo d’intensità delle due immagini. Questo corrisponde pressapoco al centro del disco di Mercurio (X,Y).
  3. Attorno alle coordinate (2X,2Y) di ciascuna immagine originale (non ricampionata) si ritaglia un quadrato di 128 pixel per lato. Questa misura è comoda perché è una potenza di 2 ed è sufficiente, nel mio caso, a contenere l’intero Mercurio. Dei due quadrati si trova la funzione di correlazione incrociata, come descritto in precedenza.
  4. La funzione di cross-correlazione può essere interpretata essa stessa come un’immagine quadrata di dimensione 128×128 pixel. Se ne trovano le coordinate del massimo assoluto e tali valori sono usati per far partire un fit gaussiano del picco d’intensità, che fornisce la posizione del massimo con precisione sub-pixel (x,y).
  5. L’immagine originale è traslata della quantità non intera (-2X-x,-2Y-y) e ricampionata alla griglia dei pixel della prima immagine interpolandola con una spline cubica.

 

Schemata informativa di "Flecken"

La schermata informativa del programma “Flecken”. Con una scala di colori modificata, sono visibili l’immagine di referenza, quella da allineare, l’immagine di cross-correlazione col contorno del fit gaussiano.

 

È finita? Certamente no, a questo punto bisogna suddividere i fotogrammi in gruppi omogenei e mediarli. Chiaramente ciascun gruppo di fotogrammi non può contenerne più di un certo numero, perché col tempo che passa il pianeta si sposta! Nel mio caso i fotogrammi erano 24018 e li ho suddivisi in gruppi da 250 l’uno, ottenendo 96 immagini finali mediate. La media è stata fatta con una funzione che pure avevo scritto in precedenza e che fa parte di una minuscola libreria in Python dove ho raccolto un po’ di operazioni utili per manipolare le immagini, gcim.

Alla fine del circo, l’algoritmo di allineamento ha riallineato tutti i 24018 fotogrammi: in 5 casi (0,02% del totale) il fit gaussiano non è andato a buon fine e si è usata una traslazione con precisione al pixel intero; in 66 casi (0,27% del totale) il fit gaussiano è converso ad un valore abnorme, probabilmente a causa di un gradiente lineare sottostante al picco gaussiano. L’identificazione di questi 66 casi è stata possibile grazie al registro dove lo script salva tutti i dettagli di cosa fa. Sto ancora pensando a un modo per permettergli di riconoscere automaticamente questi casi.
 

Confronto fotogrammi

Confronto tra due fotogrammi non allineati. Il pianeta e le macchie solari sono in posizioni diverse.

Confronto fotogrammi allineati

Confronto tra i fotogrammi precedenti dopo l’allineamento. Il pianeta è ora fisso nella stessa posizione.

 

Alla fine ho usato un programma libero, ImPPg (che pure, volendo, implementa una procedura simile per l’allineamento delle immagini) per applicare dei filtri di deconvoluzione alle immagini risultanti. Ho quindi rimontato un piccolo video accelerato, che aggiungerò presto su questo sito.

Infine, ecco lo script di allineamento: ne ho prodotto due versioni, rispettivamente “Fleck” e “Flecken”. Il primo allinea due immagini, salvando una versione traslata della seconda per farla coincidere con la prima. Il secondo fa la stessa cosa ma in modo sequenziale, lavorando su una lista di immagini da allineare. Per eseguirli avete bisogno di un’installazione di Python con i moduli Numpy, Scipy, Matplotlib. Serve anche la mia piccola libreria gcim, che è già inclusa assieme agli script. Potete scaricare tutto da qui.

Francamente, dubito che il programma, così com’è, possa servire a qualcuno per allineare le sue immagini. Non so nemmeno come si comporti con dei fotogrammi diversi da quelli del video per il quale l’ho scritto. Ma di certo può essere un buon punto di partenza, se qualcuno ha voglia di lavorarci un po’ su.

A presto!

Permalink link a questo articolo: http://www.giovanniceribella.eu/fuere/?p=1389

Post precedenti «