DS9 campi stellari

Prove tecniche per evidenziare oggetti in campi stellari affollati.

Teoricamente potrebbe essere utile, lo vedremo alla prova dei fatti.

Elaborazione dell'orbita dell'asteroide 2016 NK22 effettuata con Solex10

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