domenica 23 dicembre 2012

DS9 Photometry

Utilizzando il software DS9 Saoimage, (gratuito), è possibile fare fotometria.

ho preparato un breve Tutorial nel quale mostro le procedure che ho adottato.
 Nel video riporto la varibilita della stella CY AQR, cefeide a breve periodo.

La procedura  è stata provata anche sul transito di un pianeta extrasolare  HATP36b ripreso da Libbiano Peccioli (AAAV) il risultato e sintetizzato nel seguente grafico, dove sono stati sottratti i valori ADU tra la stella di riferimento e la stella oggetto di transito.


Transito pianeta extrasolare
 di seguito il video esplicativo (spero)

 Video tutorial

 UPDATE: agosto 2013

nelle nuove versioni di ds9  il comando analisy è:
cat {/tmp/ds9ans1981141754.fits} |  funcnts -rp stdin physical;circle(732.62345,....

questo va cosi modificato:

funds9 funcnts "ds9" "K13P10J-001UF_alig.fit" " physical;circle(732.62345,415.68855,9.9999866);circle(576.75391,411.54652,9.9999866);circle(464.48249,424.81649,10);circle(692.3805,528.36309,9.9999866);circle(464.49663,640.50542,9.9999866);circle(565.08198,351.49805,9.9999866);circle(493.6906,371.49806,9.9999866);circle(756.35901,331.68787,9.9999866);circle(782.87574,486.80406,9.9999866);circle(814.58108,569.56309,9.9999866);circle(841.57469,349.55606,9.9999866) physical;circle(627.9635,455.57182,9.9999866) "
occhio lle vigolette " " 


UPDATE


2013 migliorato lo script adesso il graico viene coretto... ( vi era un errore nell'ordinare i dati in ordine di time  )






script bash  per estrapoalre i dati:
 
l # carica le immani con ds9
rm statdat.dat
rm ds9mag.txt
clear
echo " aprire ds9 selezionare in una immagine circle per la fotometria"
echo " fare una regione di backgrouns "
echo " lanciare il programma analisy funtolls analisy region"
echo " copiare lo script e lanciarlo per tutte le regioni "
echo  " funds9 funcnts ds9 SM0321_B-*.fit  physical;circle(528.51876,518.49495,20)"
echo " inserire asteriscoper selezionare le immagini "
echo " "
ls *.txt *.dat
echo " "
echo " DS9 fotometria "
echo -n " nome del file photometrico ds9 "
read nome
ak=` cat $nome |  sed -n '/net_counts/,/followin/P' | sed '/^$/d' | sed '/#/d'  | grep -e '.[0-9]  ' | sort -gk 1 |  awk '{ print $1,$2,$3,$4,$5, -2.5*log($2) }' `
# caloclo -magnitude = zeropoint - 2.5 log10 (flux)    
echo "$ak "
echo " Creato file ds9mag.txt "
echo "$ak " > ds9mag.txt
ste=`tail -1  ds9mag.txt | awk '{print $1}'`
echo " ----------------------- "
echo " statistiche immagini " > statdat.dat
echo " ------------------------" >> statdat.dat
echo " S - Media - Max - Min  "
echo " ------------------------"
for k in  `seq 1 $ste `
do      
# max min media ADU
cat ds9mag.txt | awk '{if ($1=='$k') print $0 }' |   awk '{if(min==""){min=max=$2}; if($2>max) {max=$2}; if($2> statdat.dat
# max min media magnitudinea
cat ds9mag.txt | awk '{if ($1=='$k') print $0 }' |   awk '{if(min==""){min=max=$6}; if($6>max) {max=$6}; if($6> statdat.dat
done
echo "------------- STATISTICHE -----------------------"
echo " max min medi ADU e magnitudine "
cat statdat.dat
echo " -----------------------------------------------"
cat ds9mag1.txt
echo "-------------------------------------------------"
# estrapola dai fits la data di osservazione
gethead -tpu -n 2 DATE-OBS *.fit | cut -f2 > timehead.txt
echo  " Creato file time.txe con gethead  OBS-DATE : "
# convete la data nel formato JD
for DATE in  `cat timehead.txt `
do
        YEAR=`echo $DATE | awk ' { print substr($0,1,4) } '`
        MONTH=`echo $DATE | awk ' { print substr($0,6,2) } '`
        DAY=`echo $DATE | awk ' { print substr($0,9,2) } '`
        ORA=`echo $DATE | awk ' { print substr($0,12,2) } '`
        MIN=`echo $DATE | awk ' { print substr($0,15,2) } '`
        SEC=`echo $DATE | awk ' { print substr($0,18,2) } '`

        A=$(( $DAY - 32075 + 1461 * ( $YEAR + 4800 - ( 14 - $MONTH ) / 12 ) / 4 + 367 * ( $MONTH - 2 + ( 14 - $MONTH ) / 12 * 12 ) / 12 - 3 * ( ( $YEAR + 4900 - ( 14 - $MONTH ) / 12 ) / 100 ) / 4  )) 

        K=` echo $DATE | awk -v OFMT=%7.5f '{ print ( '$A' + (3600 * '$ORA'  + 60 * '$MIN' + '$SEC' ) / 86400 - 0.5 ) }'`

        echo -e "$K \t $DATE" >> ds9timestar.txt
done
 echo " file ds9timestar.txt creato "
  tail -5 ds9timestar.txt
  cat ds9timestar.txt | wc
  echo " adesso creo il file con tempo e dati fotometrici "
#   ste=`tail -1  ds9mag.txt | awk '{print $1}'`
   echo " stelle esaminate in fototmetria : $ste "
   for n in `seq 1 $ste`
   do
            cat ds9timestar.txt  >> aaa.txt
    done
    echo " ----------------------------------------------- "
    echo " file       mag_time.txt     pronto per essere prottato "
    echo " ------------------------------------------------ "
    paste aaa.txt ds9mag.txt > mag_time.txt
    rm aaa.txt
    rm timehead.txt
echo "  -    JD     - Date head Fits        - S  -  ADU  - err -  Com -   err -     Mag. "
echo "------------------------------------------------------"
    tail -5 mag_time.txt | awk '{print $0}'
     cat mag_time.txt | awk '{ print  $1,$3,$4}' | grep " 1 " > zyzy.txt
for g in `seq 2 $ste`
do
cat mag_time.txt |  awk '{ print  $3,$4}' | grep "$g " |  cut -d ' ' -f2 >  $g-zy.txt
done
paste zyzy.txt *-zy.txt >> starcol.txt
rm *zy.txt
 echo "----------------------------------------------------------------"
 echo " creato file starcol.txt stelle in colonna  (nr $ste stelle )    "
 echo "----------------------------------------------------------------"
head -5 starcol.txt

# fine

C/2011 UF305 (LINEAR)

La cometa C/2011 UF305 è stata scoperta dall'osservatorio 704 Lincoln Laboratory ETS, New Mexico, coma da circolare M.P.E.C. 2011-Y51del 28/12/2011.

    C/2011 UF305 (LINEAR)
Epoch 2012 July 12.0 TT = JDT 2456120.5
T 2012 July 22.15968 TT                                 MPC
q   2.1382045            (2000.0)            P               Q
z  -0.0003117      Peri.  121.99192     -0.29610050     -0.35829972
 +/-0.0000006      Node   297.43591     +0.07008823     +0.91632627
e   1.0006665      Incl.   93.96064     +0.95258183     -0.17879452
From 1200 observations 2011 July 25-2012 Dec. 10, mean residual 0".4.
 
La Notte del 18/12/2012 è stata osservata da San Marcello, quando aveva una luminosità di ~ 14.8, D=1.966 r=2.738 UA.

Ecco l'immagine composta dallo stack di 17x120"

Clicca per ingrandire


sabato 22 dicembre 2012

(19367) Pink Floyd


Nel corso del follow-up di un asteroide nella NEOCP, nel campo di ripresa c'era un asteroide particolare.... almeno per il suo nome:



(19367) Pink Floyd, asteroide scoperto il 12(03/1997 dall'osservatorio di Caussols (MPC 010), è stato dedicato al grande gruppo musicale Pink Floyd.



Citazione:

(19367) Pink Floyd = 1997 XW3
     An influential English progressive rock group, Pink Floyd wrote
several astronomically themed songs, such as Interstellar Overdrive
and Astronomy Domine.  The group's 1973 album
Dark Side of the Moon is one of the classic rock recordings
 
Elementi orbitali 
 
Epoch 2012 Sept. 30.0 TT = JDT 2456200.5                MPC
M  10.38992              (2000.0)            P               Q
n   0.25767238     Peri.  305.50161     +0.79580772     -0.60213071             T = 2456160.17778 JDT
a   2.4458180      Node    91.60724     +0.57431811     +0.71686563             q =     2.0426651
e   0.1648336      Incl.    3.68552     +0.19196038     +0.35148581
P   3.83           H   14.4           G   0.15           U   0 

La notte del 22/12/2012 riesco a immortalare l'asteroide, che soddisfazione !!

Clicca per ingrandire

La canzone la dedico a Giovanni Sostero.


C/2012 Y1 (LINEAR)

Il Minor Plante Center con circolare M.P.E.C. 2012-Y18 del 21/12/2012, comunica la scoperta di una nuova cometa C/2012 Y1 (LINEAR) effettuata dall'osservatorio 704 Lincoln Laboratory ETS, New Mexico. la notte del 18/12/2012.

Quando l'oggetto era ancora nella NEOCP , è stato osservato da San Marcello confermando la natura cometaria (ma siamo arrivati in ritardo nell'invio delle misure)

ecco la nostra immagine  di conferma mag. ~ 18.5

clicca per ingrandire
Orbital elements:
    C/2012 Y1 (LINEAR)
T 2012 Dec. 8.59572 TT                                  MPC
q   2.1202862            (2000.0)            P               Q
                   Peri.  242.56443     +0.13525802     -0.98109896
                   Node   201.08074     +0.98918745     +0.14170371
e   1.0            Incl.   22.62749     +0.05668731     -0.13177590
From 30 observations 2012 Dec. 18-21.

C/2012 X1 (LINEAR)

Scoperta una nuova cometa dall'osservatorio 704 Lincoln Laboratory ETS, New Mexico, la notte del 08/12/2012.

Qunado l'oggetto era ancora nella NEOCP veniva osservato da San Marcello.

Il Minor Planet Center con circolare M.P.E.C. 2012-X70 del 12/12/2012, comunica la scoperta della nuova cometa.

Ecco l'immagine di conferma effettuata da 104, quando la cometa aveva una luminosità di ~ 19.0 mag.
Somma di 13 immagini da 120 secondi ciascuna.

clicca per ingrandire

Orbital elements:
    C/2012 X1 (LINEAR)
T 2014 Mar. 11.64562 TT                                 MPC
q   1.7441262            (2000.0)            P               Q
                   Peri.  128.68621     -0.23496663     +0.71451180
                   Node   114.25932     -0.95060712     -0.31038737
e   1.0            Incl.   46.28737     +0.20282207     -0.62700285
From 34 observations 2012 Dec. 8-12.

mercoledì 19 dicembre 2012

C/2012 K5 (LINEAR)

Ripresa dall'osservatorio astronomico di San Marcello la bella cometa C/2012 K5 (LINEAR), la notte del 18/12/2012.

La cometa scoperta da 704 Lincoln Laboratory ETS, New Mexico.  è stata osservata da LIBBIANO B33 quando ancora era nello NEOCP il 27/05/2012.

C/2012 K5 (LINEAR)
Epoch 2012 Sept. 30.0 TT = JDT 2456200.5
T 2012 Nov. 28.6881 TT                                  MPC
q   1.141790             (2000.0)            P               Q
z  +0.001314       Peri.  139.2908      -0.1510840      -0.0652734              T = 2456260.18812 JDT
 +/-0.000001       Node   279.0385      +0.4230551      +0.8975520              q =     1.1417902
e   0.998499       Incl.   92.8475      +0.8934193      -0.4360502     Earth MOID = 0.29415 AU
1/a(orig) = +0.001941 AU**-1, 1/a(fut) = +0.001801 AU**-1.


Adesso la cometa mostra una belliisima coda.

Al momento dell'osservazione la sua luminosità era di ~ 11.0 mag, con una distanza D=0.459  r=1.18.  e con un motion di ~ 5.40"/min.

Di seguito un animazione composta da 26 frame ognuno con esposizione di 45 secondi.


Animazione Clicca per ingrandire
 Somma di 20 immagini da 45 secondi ciascuna


Comet C2012 K5 clicca per ingrandire

la stessa immagini in falsi colori

clicca per ingrandire


C/2012 S1 ISON

Ripresa stanotte la cometa C/2012 S1 ISON dall'osservatorio di San Marcello (MPC 104)

La cometa ha un aspetto stellare,e non si è notato alcun accenno di coda....
Se le previsioni, rispettano le aspettative questa cometa potrebbe essere visibile di giorno il prossimo novembre 2013.

Al momento ha una luminosità "telescopia" di 16. mag. D=4.515 r=5.391.
(HA da poco passato l'orbita di Giove !!!)

Al momento della ripresa la cometa era rincorsa da alcuni asteroidi

Animazione clicca per ingrandire
Particolare della cometa ISON,  somma di 27X120" ,si nota un leggero allungamento verso ovest di ~21.5".




Particolare della cometa

Vedremo cosa ci riserverà in futuro.

lunedì 17 dicembre 2012

(99942) Apophis

L'asteroide (99942) Apophis = 2004 MN4 , è stato scoperto dal "nostro" Fabrizio Bernardi, Tucker, R. A., Tholen, D, la notte del 19/06/2004 dall'osservatorio astronomico di Kitt Peak. Conosciuto per essere stato il primo asteroide ad avere una certa probabilità di impatto con la Terra, con una probabilità stimata pari a 1 su 6250, Successive misure astronomiche hanno ridotto questa possibilità  a meno di 1 su 250.000.
Nel 2029 l'asteroide darà spettacolo tanto da poter essere visto ad occhio nudo, in quanto la sua luminosità raggiungerà la 3.3 mag. quando si avvicinerà a soli 0.0002563 UA ( a circa 38.000 km) dalla Terra.
Le dimensioni di Apophis sono stimate in ~ 350 metri.

Il giorno 17/12/2012, utilizzando il telescopio remoto H06 iTelescope Observatory, Mayhil USA, 0.25-m con risoluzione di 1.6"/pixell, sono riuscito a riprendere per la "MIA" prima volta il famoso asteroide.

Orbita di Apophis 

Ultime osservazioni


Dalle effemeridi del Minor Planet Center, noto che l'ultima osservazione di Apophis è stata fatta dall'osservatorio F51 PANSTARRS in data 08/05/2012.

L'oggetto decisamente basso si trovava  a 30° sull'orizzonte, di magnitudine stimata intorno alla 17.3 mag. con motion di 2.64"/min e PA 266.4.

Vale comunque la pena di provare a rintracciare per questa opposizione l'asteroide.

Effettuo una sequenza di 10 immagini da 10 secondi ciascuna.

Sommo le immagini acquisite sul moto di Apophis, e le blinko...

Sull'imagine un "puntino" si muove secondo le attese.

Davvero una grande soddisfazione, vedere il famoso asteroide Apophis


Animazione clicca per ingrandire

Sommo tutte le immagini acquisite per evidenziare l'asteroide (99942) Apophis.

Clicca per ingrandire


 Dal sito NEODyS si ricavano i prossimi passaggi ravvicinati dell'asteroide e le probabilità di impatto con la Terra.


passaggi ravvicinati



Probabilità di impatto




domenica 9 dicembre 2012

(4179) Toutatis

L'asteroide (4179) Toutatis, il giorno 11 dicembre 2012 passera "vicino alla Terra", ad una distanza di circa 0.046 UA (6 881 600 km) 17,92 LD. (Fonte MPC)


L'asteroide NEA ha dimensioni di circa  5.4 Km oggetto di studio da parte dei radiotelescopio di Goldstone


Immagine radar


Andrea Mantero dall'osservatorio C77 Bernezzo Observatory (Cuneo) , Nello Ruocco dall'osservatorio C82 Nastro Verde (Sorrento), ed io utilizzando il telescopio remoto Itelescope I89 RAS Observatory, Nerpio la sera del 9 dicembre hanno ripreso in contemporanea l'asteroide.

Utilizzando le immagini riprese in contemporanea dall'osservatorio C77 Cuneo e I86 Madrid (con una differenza di soli 2 secondi) è stato possibile apprezzare lo spostamento sulle immagini dovute alla parallasse.

I due osservatori distano 1014.15 km, permettendo di misurare un angolo di parallasse di circa 28.97"

Di seguito le immagini utilizzate per la determinazione dell'angolo di parallasse





SI evidenzia la parallasse


Misurazione dell'angolo di parallasse dell'asteroide (4179) Toutatis


I89 C77  differente posizione dell'asteroide (misura provvisoria)

Composizione delle immagini dove si evidenzia la differente posizione prospettica dell'asteroide (4179) Toutatis



Misura della parallasse migliorata.






 Successive analisi permetteranno di calcolare la distanza dell'asteroide (4179) Toutatis dalla Terra, dalle misure di posizione ottenute dai due osservatori.

Conclusioni preliminari


Al momento utilizzando le procedure illustrate sul numero della rivista (UAI) ASTRONOMIA n. 4 Luglio - Agosto 2012 l'articolo sulla parallasse dell'asteroide (433) Eros. e,

considerando la distanza dei due osservatori 1068 km è l'angolo di parallasse 29.5, risulta che (4179) Toutatis al momento dell'osservazione di trovava 7 457 381 km dalla Terra (0.0498 UA) compatibile con le previsione del NEODyS 7 405 095 (0.0495 UA) con una differenza di -52 286 km corrispondenti ad una differenza di -0,00035 UA, stima entro la deviazione standard

Animazione ottenuta con le immagini prese in remoto I89 Itelescope 0.25-mt, ripresa in Binnig 2x2, esposizione 15 secondi ciascuna, risoluzione 1,5" pixel, 1536x1024 pix.


Animazione clicca per ingrandire

Animazione clica per ingrandire






venerdì 7 dicembre 2012

Addio Giovanni

Oggi, 7 dicembre 2012,  mi è stato comunicata un triste notizia.
La scomparsa di Giovanni Sostero,  uno tra i migliori astrofili al mondo, il suo contributo è stato fondamentale per molti astrofili.

Per me era un esempio, un precursore da seguire.

Addio Giovanni

wikipedia Giovanni_Sostero

Giovanni Sostero INAF

http://remanzacco.blogspot.it/







venerdì 30 novembre 2012

M67 prove fotometriche

Mi sono divertito a fare alcune prove per automatizzare (o quasi) la ricerca dell'indice di colore in campi stellari acquisiti con filtri fotometrici RBVIU.

M67 VBR pseudo color
Premetto che non ho ancora calcolato il coefficiente di risposta dei filtri, per cui i dati qui presentati sono provvisori.

Mi sono collegato in remoto con ITELESCOPE utilizzando il telescopio T5 0.25-m, acquisendo immagini dell' ammasso aperto M67, utilizzando il filtri V R e B.

alcune stelle individuate con sextractor
Per individuare le stelle all'interno delle immagini utilizzo sextractor il quale mi crea un file per contenete la magnitudine strumentale, le coordinate AR e Dec (le immagini sono state precedentemente astrometrizzate)

 Nell'immagine a fianco un esempio delle stelle trovate da sextractor, filtrate con mag  17.

Una volta estrapolato dalla immagini i dati, ho tre file  contenenti la magnitudine strumentale per ogni filtro.




#   1 MAG_ISO         Isophotal magnitude                             [mag]
#   2 MAG_AUTO        Kron-like elliptical aperture magnitude         [mag]
#   3 MAG_BEST        Best of MAG_AUTO and MAG_ISOCOR                 [mag]
#   4 X_IMAGE         Object position along x                         [pixel]
#   5 Y_IMAGE         Object position along y                         [pixel]
#   6 ALPHA_J2000     Right ascension of barycenter (J2000)           [deg]
#   7 DELTA_J2000     Declination of barycenter (J2000)               [deg]
#   8 FWHM_IMAGE      FWHM assuming a gaussian core                   [pixel]
#   9 FLUX_MAX        Peak flux above background                      [count]
#  10 FLUX_BEST       Best of FLUX_AUTO and FLUX_ISOCOR               [count]
 10.9872  11.0272  10.9809    712.247     92.118 132.8639646 +12.1275237    27.32     65442.18      2556548
 11.7438  11.7824  11.7824   1095.791     74.601 132.6827019 +12.1279113     4.85     65451.47      1221912
 18.7755  17.9220  18.2863   1543.149     28.789 132.4707467 +12.1399173    12.84     78.95376     3058.259
 14.8945  14.9164  14.9164   2101.892     24.988 132.2071484 +12.1300568     4.95     6021.165      68144.2
 16.5219  16.5359  16.5359   2089.122     14.418 132.2129437 +12.1351981     4.21     1606.215     15334.09

 Vi è una correlazione tra il c.d. indice di colore  e la luminosità delle stelle che permette di individuare, la classe spettrale delle stelle come dal diagramma Hertzsprung–Russel (HR).

B-V e V-R temperatura
Dall'indice di colore è possibile indicativamente risalire alla classe spettrale, in rete si trovano moltissime informazione, io utilizzo questi dati per correlare i dati.
Nel grafico seguente sono plottare i valore dell'indice colore B-V e V-R, e la relativa temperatura superficiale delle stelle e quindi alla Classe Stellare.





 In pratica, se dalla nostra analisi troviamo un indice di colore B-V = 0.30 o V-R=-0.13, la stella ha una temperatura di 26500 corrispondente alla classe spettrale tipo B0V.

Utilizzando le immagini di M67 ( ammasso standard per la configurazione fotometrica) effettuate con i filtri B e V , e riportando per ogni stella la differenza tra la magnitudine strumentale del filtro B e quella del filtro V ( B-V), ottengo il seguente grafico, dove il colore del punto stellare corrisponde all'indice colore trovato e di conseguenza alla temperatura.


M67 Virtuale clicca per ingrandire

I dati plottati con indice di colore e magnitudine strumentale si ottiene un grafico che assomiglia al diagramma HR.

M67 HR clicca per ingrandire
Costruiamo adesso il grafico delle frequenze 
M67 Frequenza
Cosi possiamo notare che utilizzando i filtri B-V abbiamo la maggior parte di stelle hanno un indice colore tra 0.88 - 0.94 ( G8IV  5100 T - G8III   4870 T).

Con lo stesso procedimento mettiamo in grafico i dati relativi alla differenza tra il filtri V e R.


M67 V-R
 Il risultato ottenuto è in linea con quando in precedenza fatto da altri.
segnalo un interessante articolo sulla rivista online AstronomyOnline.net, che tratta in modo molto piu approfondito l'argomento. 

Concludendo, con il metodo semi-automatico di estrapolazione e analisi delle immagini astronomiche, nel video meglio descritto, sono riuscito ad ottenere un discreto risultato, nonostante non sia stata fatta alcuna calibrazione, che sicuramente avrebbe permesso di ottenere un risultato migliore.


Il Video mostra le procedure che ho utilizzato per ottenere i dati esposti.


 UPDATE

 Ho calcolato i coefficienti di correzione per il telescopio T5 NM (USA) Itelescope

Utilizzando le stelle standard dell'ammasso M67, (64 stelle), ho calcolato per ogni filtro il valore di slope utilizzando in ascisse i valori della mag strumentale ed in ordinate la magnitudine standar delle stelle in M67, come nel grafico a fianco.

Per i filtri ho trovato i seguenti valori :
 
  • Filtro V = 0.996(v)-2.813
  • Filtro B = 0.894(b)-2.940
  • Filtro R = 0.999(r)-2.675 
Dove v,b,r è la magnitudine strumentale estrapolata dalle imaggini con sextractor.


Dalla magnitudine strumentale, ho ricavato la magnitudine ridotta.
Confrontando l' Indice di colore delle 64 stelle standar in M67, e l'indice di colore ricavato utilizzando la magnitudine ridotta  delle stesse stelle individuate da sextractor ricavo il seguente grafico, che mostra la differenza dell'indice di colore dei filtri rispetto al sistema standard.

B-b, V-v R-r


Nel grafico sotto, la differenza tra l'indice di colore calcolato con la magnitudine strumentale e la magnitudine corretta.

Differenza tra mag strumentale e corretta B-V


Ed infine ecco il grafico HR dell'ammasso M67 nei filtri B-V

M67  B-V


 ed in ultimo il grafico V-R corretto




Adesso applichiamo la stessa procedura sull' ammasso NGC 1245.
In questo caso utilizzo solo le immagini ottenute con i filtri B e V, in quanto l'immagine in filtro R risulta mossa.

L'ammasso NGC 1245 è molto compatto pertanto ho deciso di estrapolare soltanto le stelle nelle vicinanze dell'ammasso, selezionandole con il software DS9 come nell'immagine.
Selezione delle stelle da esaminare in NCG 1245
Una volta selezionate le stelle da analizzare, in tutto 150, ho provveduto a calibrarle ed ecco il risultato.

NGC 1245 HR filter B-V
Adesso vediamo come le stelle dell'ammasso sono distribuite nello spazio secondo l'indice di colore B-V.
NGC 1245 B-V
Ed infine analizziamo la frequenza dei dati

Frequenza filtro B-V
NGC 1245 B-V > 1.9
Per finire un animazione filtrata per codice colore B-V con steep di 0.5

Animazione
Per maggiori informazioni sull'ammasso NGC 1245 consultare l'articolo NGC 1245 - an intermediate age open cluster.

Adesso vediamo l'ammasso M3
M3 BRV color
Anche in questo caso il diagramma è molto simile ad altri lavori svolti su questo ammasso globulare, come nell'artico di Sahana Kumar, ovviamente considerando la strumentazione utilizzata per effettuare questo lavoro.

M3 B.V


Stelle esaminate in M3 con indice colore B-V

Confronto tra l'ammasso M3 e NGC 1245, si nota che hanno una distribuzione similare, NGC 1245 risulta avere un indice colore più "freddo" rispetto a M3.


Sotto particolare del precedente diagramma, la riga di colore BLU rappresenta il diagramma HR. Si noto sia in M3 che in NGC 1245 la linea del turn-off.

M3 e NGC 1245 B-V.




Confrontiamo graficamente l'indice delle frequenze dei due ammassi globulari
M3 NGC 1245 frequneze

 
backman