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

Nessun commento:

Posta un commento