Paolo ha scritto:Marco, sai cosa consiglia Ernst?

Paolo,
Ernst consiglia di utilizzare HRV-MM. Il tool che hai trovato sullo stesso sito credo che utilizzi il metodo della cross-correlation non il metodo del mirroring utilizzato da HRV-MM.
Buon fine settimana.
Marco Leonardi

Ciao Paolo e Marco.
Grazie Marco per le ulteriori dritte.
Ho fatto una prova cercando di togliere un po' le telluriche, che come dice Paolo potrebbero dare fastidio quando sono così intense come nel mio caso.
Poi ho attivato l'opzione "Advanced Parameters of spectral region" cercando di mettere un range che "ad occhio" ho reputato ragionevole.
A questo punto, seguendo la tua indicazione di cercare di vedere le righe orizzontali a sx e dx ho scelto i valori della "Correlated spectral region".

Anche in questo caso i valori cambiano al variare dei vari parametri, ma con un po' di prove ho visto che il risultato tende a convergere verso il valore di HRV = -26.6 km/s (fra l'altro identico a quello fornito dal Simbad) usando valori della "Correlated spectral region" che andavano da 11 a 15 Angstrom. Con valori più alti o più bassi, l risultato diverge.

Non so se questa sia una strada corretta, ma da questa prima prova direi che sembra buona.
Devo però fare altre riprove e non è detto che, data la relativa arbiitarietà nella scelta dei vari parametri, riesca sempre a trovare un risultato corretto e convergente.
Opinioni in merito?

Grazie
Nico

Facendo seguito al mio precedente messaggio, ho fatto una nuova prova su un vecchio spettro di Vega e dopo innumerevoli variazioni dei vari parametri direi che il risultato converge abbastanza nettamente, almeno per quel che vedo io, verso una HRV = -13.3 km/s.
Il problema però è un altro.
Il Simbad fornisce per Vega una Velocità radiale = -20.6 +/- 0.2 km/s. Fonte "Pulkovo Compilation of Radial Velocities for 35 495 Hipparcos stars in a common system" del 2006.

Wikipedia invece https://en.wikipedia.org/wiki/Vega ,mostra a Velocità radiale = -13.9 +/- 0.9 km/s.
Chi ha ragione?
Può essere che la VR di Vega sia leggermente variabile?
Se non ricordo male qualche tempo fa con Paolo ne avevamo parlato.... o forse ricordo male io?

Edit:
Ops... non ricordavo di aver aperto proprio io su questo forum l'anno scorso un topic dove si accenna alla "moltitudine" di valori in bibliografia per la VR di Vega.
viewtopic.php?f=22&t=3301&hilit=vega


Nico

Benissimo Nico! Usato in quel modo (regione relativamente stretta attorno al picco) mi pare di capire che sia equivalente a trovare la "Rv core" nell'analisi piu' estesa (di tutta la riga). Ok Marco, ho scaricato solo HRV-MM, se l'altro tool usa la cross-correlation su n-righe posso farne a meno perché abbiamo una analoga funzione integrata in ISIS.

Paolo

Ciao Nico,

sono molto contento che tu sia riuscito nell'intento e complimenti per la determinazione !

Per la comprensione di tutti riassumo le operazioni eseguite dal tool HRV-MM utilizzando due delle tue figure che si riferiscono al tuo test su Altair:
- è stato applicato il metodo delle bisettrici (sul quale si basa il metodo mirroring) sull'area del picco che hai selezionato con la linea rossa (utilizzo queste figure del mio post precedente per spiegare il modello delle bisettrici):

- è stato calcolata la variazione della RV lungo la linea delle bisettrici che tagliano l'area selezionata alle diverse altezze;
- è stata calcolata la BVS [differenza dei valori medi di velocità al top (Vtop) e al bottom (Vbottom) ];

- la RV e la BVS sono state correlate ed è stata calcolata la risultante RV.

Ci sono altre informazioni utili che risultano dall'utilizzo di questo tool che però non mi sono chiare. In particolare non mi è ancora comprensibile come interpretare il grafico presentato nel foglio (function of CC). Posso intuire che sull'asse delle ascisse c'è la RV determinata ma non riesco a capire cosa esprime l'asse delle ordinate. E' inoltre probabile che l'andamento del grafico ci aiuti a capire la correlazione RV e BVS, e di quanto siamo in prossimità del valore di RV di convergenza, ma non ne sono sicurissimo. Contatterò Ernst Pollmann per avere qualche dritta che condividerò immediatamente con il nostro gruppo non appena mi sarà possibile.

A presto

Marco Leonardi

Ciao Nico, Paolo,
come anticipato ho contatto Ernst Pollmann con il quale siamo rimasti in accordo di aprire un simile post anche in ARAS. Aprirò il post domani in modo che Ernst possa agganciarsi e di conseguenza raccogliere anche i pareri di altri esperti sull'argomento "calcolo velocità radiale" e considerazioni sui diversi metodi utilizzati.

Nell'email di Ernst viene chiarito come certi parametri vengano selezionati in funzione della dispersione del ns spettro e del potere risolutivo. Si capisce che la curva (function of CC) idealmente debba essere una gaussiana, questo particolare guida tutto. Incollo qui sotto la risposta che mi ha inviato e che mi sembra chiara. Comunque ditemi la vostra opinione/dubbi.

Un abbraccio !
Marco Leonardi

Dear Marco,
I am very pleased of your attempts, to calculate the RV by use of the PDM method.
I will try it too with your last pi Aqr spectrum. In case of that spectrum, I would prefer not to use selected sections of the profile, but to use the entire profile (see picture "spectrum").
You can recognize the corresponding parameter for that process as "advanced parameters of spectral region": high level = 2.3, low level = 0.
Corresponding your spectral parameters, like dispersion and spectral resolution power, you have to select certain values for "correlated spectral region" (I did choose 30 Angstr.) and for "velocity range of correlation" (I did choose 110 pixel). Laboratory wavelength = 6562.852 Angstr. With this pre-settings you can start the measurement "correlated line" (see picture "correlated line") and function of CC (see picture "function of CC"). Sometimes you have to "play" with this settings, depending on the line profile, in order to find out the best result of the CC.
The CC determines the size of the shift of the spectrum and/or the comparison (mirrored) spectrum (both leads to the same results, if one considers the direction of the shift), which is necessary, in order to obtain a maximum agreement of the correlated spectra (the measure of the agreement is defined by the coefficient of correlation).
An indicator for a "best result" is the CC curve itself. It has to be (as good as it possible) a symetrical distributed Gaussian function (curve).
Unfortunately there is no reference of a certain correlation coefficient in that program, so the only possibility you have is to judge only the CC curve.
I hope that this helps to understand the evaluation process and the CC a bit better.


Best wishes, Ernst
----------------------------------------------
International Working Group ASPA
Active Spectroscopy in Astronomy
http://www.astrospectroscopy.de
http://www.astronomie.de/astronomische- ... troskopie/

Ciao Marco, grazie per gli approfondimenti. Ho appena dato un'occhiata al forum ARAS e devo dire che la faccenda e' complessa anche piu' di quanto immaginavo. Sulle righe simmetriche non ci sono problemi mentre per la misura di righe asimmetriche bisogna necessariamente comprendere la natura delle sorgenti luminose (binarieta', condensazioni/asimmetrie del disco, ecc) altrimenti non si sa quello che si sta facendo. Comunque e' un bene avere a disposizione un tool con cui sperimentare, specie sotto la guida di persone piu' esperte.

A presto!
Paolo