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 |
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
# 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
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
