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/