Oggi è 28 mar 2024, 12:07
   
Text Size
Login

Primo spettro di delCyg

Le nuove tecniche per studiare e comprendere la natura di tanti fenomeni celesti.

Messaggioda umberto » 13 gen 2019, 18:54

Tony, per questa operazione ho chiesto aiuto a tutti i fotografi del forum che usano maximDL. Se non avrò risposta chiederò al mio amico Joan oltre Pirenei, vedrai che riceveremo tutti i dettagli per l'operazione. Umberto
Ultima modifica di umberto il 13 gen 2019, 21:56, modificato 1 volta in totale.
umberto
Quasar Guru
 
Messaggi: 1590
Iscritto il: 11 gen 2011, 23:07

Messaggioda Paolo » 13 gen 2019, 19:33

Con il solo scopo di andare avanti nella tecnica possiamo anche usare lo spettro di libreria che dici. Gli spettri della libreria Miles hanno una risoluzione molto piu' alta rispetto a quanto si puo' ottenere con lo SA100, quindi bisogna lavorare un pochino nella fase di divisione. I profili della libreria Pickles hanno una risoluzione leggermente piu' bassa e arrivano a comprendere la parte di spettro pura (fino a 7800A circa). I Miles arrivano sui 7400A circa per cui, i tuoi profili verranno automaticamente tagliati a questa lunghezza d'onda.

Comunque anche Vspec permette la divisione dei pixel

Hai gia' perso tutti i conteggi sopra i 32767 ADU nel momento in cui apri l'immagine in VSpec! Bisogna operare prima, questo software non mantiene i valori di intensita' superiori a 32767 ADU (li livella a questo valore).

Appena riesco provo la correzione per la risposta strumentale del nostro profilo con VSpec. Se vuoi provare ti anticipo che e' meglio lavorare sul profilo croppato (3600-7800 circa), laddove cioe' esiste il segnale della stella, e non su quello completo che contiene l'ordine zero (se provi capirai il perche').

Paolo
Avatar utente
Paolo
Quasar Guru
 
Messaggi: 7747
Iscritto il: 16 gen 2006, 22:49
Località: L'Aquila

Messaggioda Paolo » 14 gen 2019, 13:02

Proviamo!

Abbiamo caricato in VSpec il nostro spettro calibrato solo in lunghezza d'onda (sopra) e lo spettro che hai proposto (sotto).

VS01.jpg

Selezioni il tuo spettro, poi dal menu Operation -> Divide profile by profile, fai doppio click sul nome del file di libreria (s0262.dat) e un click su Intensity. Poi OK.

VS02.jpg

Il risultato della divisione e' quella "robaccia" verde che termina con un picco

VS03.jpg

Occorre ritagliare la zona utile selezionandola (lascia fuori il picco, vedi selezione nell'immagine sopra) usando Edit -> crop

VS04.jpg

Il mio consiglio a questo punto e' quello di salvare la divisione grezza appena ottenuta. Siccome c'e' da lavorare un po' su questo profilo, in caso di errore basta ricaricarlo in VSpec (altrimenti devi rifare di nuovo tutti i passaggi).

Dal menu Edit -> Replace -> intensity, premi OK.

VS05.jpg

File -> Save as (chiamare risposta strumentale)

VS06.jpg


Ti prego di notare che la curva di risposta cosi' ottenuta e' grezza. I picchi che vedi in corrispondenza delle righe di Balmer non fanno parte della risposta del sensore (e di altri elementi della "catena osservativa") ma sono veri e propri artefatti dovuti principalmente alla differenza di risoluzione dei due spettri (il tuo e quello della libreria Miles). Pertanto vanno eliminati dalla curva di risposta che dobbiamo ricavare.

Oltre questi vistosi picchi puoi notare che il continuo non e' proprio "liscio" ma tutto "seghettato", in gergo diciamo che e' rumoroso. E' preferibile togliere anche questo rumore. In definitiva vogliamo una curva di risposta liscia, priva di picchi e discontinuita'. Se non la puliamo in questo modo, i picchi e il rumore saranno trasportati nel profilo finale, un bel guaio! Queste fasi si chiamano "estrazione del continuo" e "smussatura" (della curva di risposta).


continua...
Avatar utente
Paolo
Quasar Guru
 
Messaggi: 7747
Iscritto il: 16 gen 2006, 22:49
Località: L'Aquila

Messaggioda tonyven » 14 gen 2019, 14:54

Benissimo. Ho letto qualcosa sullo smhooting e vedere fatto in maniera diversa non può che farmi piacere. Non so come ringraziarti per la pazienza che hai con me.
tonyven
Quasar Dipendente
 
Messaggi: 621
Iscritto il: 10 gen 2019, 17:54

Messaggioda Paolo » 14 gen 2019, 15:05

Figurati, per me e' un piacere! L'interazione possibile su una piattaforma come questo (o altri) forum credo sia insostituibile. Perche' non sfruttarla?

Proseguendo...

In VSpec hai piu' di un modo di procedere per estrarre il continuo. Sono validi tutti, qui uso la rimozione di parti del profilo. Queste vengono poi ricostruite per interpolazione.

Dal menu Radiometry -> Compute Continuum

Uno alla volta, selezioni i picchi (non tenerti molto stretto, vedi il mio esempio sotto) e premi l'icona con la X. Il tratto di profilo nella selezione scompare.

VS07.jpg

Quando hai finito con tutti i picchi, premi l'icona con il "play". Dovresti vedere una linea arancione che attraversa il continuo

VS08.jpg

A questo punto devi applicare la smussatura. VSpec apre il piccolo tool con uno slide ed alcuni parametri da inserire. Prova i parametri che vedi qui sotto e cambiali se necessario fino a vedere una curva smussata (senza rumore) che segue il continuo in modo simile:

VS09.jpg

Premi OK e salva la preziosa curva di risposta. Ti ricordo che devi prima usare Edit -> Replace -> intensity, poi File -> Save as.

E' arrivato il momento di applicare la curva di risposta.

Carica in VSpec il tuo profilo calibrato solo in lunghezza d'onda e la curva di risposta. Seleziona il tuo profilo e dal menu Operation -> Divide profile by profile seleziona la risposta (io l'ho chiamata risposta_ok) e Intensity.

VS10.jpg

Il risultato e' il profilo verde:

VS11.jpg

Come vedi, va ritagliato. Finalmente il profilo finale!

VS12.jpg

Finezza, se vuoi. Normalizzare la zona vicino alla riga h-alfa (e' una procedura molto comune fra gli spettroscopisti).

Selezioni la regione (vedi sotto), poi dal menu Operation -> normalize:

VS14.jpg

Nota la scala delle ordinate a sinistra del grafico (i valori erano prima molto alti). Adesso il profilo ha l'intensita' pari a 1 nella regione che abbiamo selezionato e raggiunge circa 5 nel blu.

Fammi sapere se ottieni risultati simili o incontri qualche problema. La parte piu' delicata e' l'estrazione del continuo e la smussatura. E' la fase purtroppo piu' soggettiva dell'intero processo, in grado di alterare in modo anche pesante il risultato finale.

Volevo farti notare infine che qui non abbiamo sottratto il fondo nell'immagine dello spettro. Senza questa (semplice) operazione, la curva di risposta non e' applicabile ad altri spettri (vale solo per la stella di riferimento). Ne parleremo.

La procedura e' da ripetere con un nuovo spettro di migliore qualita'. Aspettiamo le nuove osservazioni!

Paolo
Avatar utente
Paolo
Quasar Guru
 
Messaggi: 7747
Iscritto il: 16 gen 2006, 22:49
Località: L'Aquila

Messaggioda tonyven » 14 gen 2019, 23:14

ALLELUIA! ALLELUIA! A me sembra molto simile al tuo
sirio-014 div2calibratacorretta.png
sirio-014 div2calibratacorretta.png (4.75 KiB) Osservato 2398 volte


E questa è la risposta strumentale
risposta strumentaleok.png
risposta strumentaleok.png (3.6 KiB) Osservato 2398 volte


Come puoi vedere il profilo calibrato e corretto è ancora verde.Anche qui Edit->Replace ?
(So troppo contento =p~ =p~ )

Ora non vedo l'ora di acquisire un'altra stella da usare come riferimento. Ovviamente la risposta strumentale è valida solo per la sera in cui e' stata acquisita? Cioè ogni serata ha una sua risposta strumentale?
Proverò ad acquisire qualche stella come mi è stato suggerito da Umberto in questo post.
RAGAZZI CHE SODDISFAZIONE!!!! Mi si passi il francesismo ma qui è meglio di una trom...ta!
tonyven
Quasar Dipendente
 
Messaggi: 621
Iscritto il: 10 gen 2019, 17:54

Messaggioda umberto » 15 gen 2019, 0:53

:YMAPPLAUSE: Hip Hip Hurraaaaaaaaaa!!!!!!!.... "Che cavalcata!"
umberto
Quasar Guru
 
Messaggi: 1590
Iscritto il: 11 gen 2011, 23:07

Messaggioda Paolo » 15 gen 2019, 11:42

ALLELUIA! ALLELUIA! A me sembra molto simile al tuo

Certamente! Le uniche differenze sono dovute al maggiore smussamento che hai applicato alla curva di risposta.

Qui sotto vedi la differenza nelle nostre curve.

comp_risp.jpg

In pratica hai "tirato" maggiormente lo slide sulla sinistra con i suoi parametri. Considera che se tutte le righe sono state correttamente rimosse con il tool, puoi permetterti di smussare un po' meno (vedi la mia curva di risposta). Questo ti darebbe il vantaggio di attenuare le oscillazioni (ripple) introdotte dalla tua camera digitale.

Comunque va bene anche cosi' come hai fatto!

Oltre le stelle che ti ha consigliato Umberto, quando hai tempo prova a prendere lo spettro di qualche stella peculiare come gamma Cassiopeiae (mag 2.2) e P-Cygni (mag 4.8, bassa sull'orizzonte ovest) che hanno la riga h-alfa in emissione.

Se osserverai stelle "fredde" come Betelgeuse o Aldebaran (classi M, K), ti faremo vedere come si procede per calibrare in lunghezza d'onda (con lo Star Analyser non vedi le righe di Balmer in questi spettri).

Adesso comincia il divertimento!

Paolo
Avatar utente
Paolo
Quasar Guru
 
Messaggi: 7747
Iscritto il: 16 gen 2006, 22:49
Località: L'Aquila

PrecedenteProssimo

Torna a Spettroscopia, Radioastronomia

Chi c’è in linea

Visitano il forum: Nessuno e 6 ospiti

cron

Chi c’è in linea

In totale ci sono 6 utenti connessi :: 0 iscritti, 0 nascosti e 6 ospiti (basato sugli utenti attivi negli ultimi 5 minuti)
Record di utenti connessi: 595 registrato il 22 dic 2022, 1:59

Visitano il forum: Nessuno e 6 ospiti

Login Form