Difference: 20160214MSciPhotometryAnalysisObs3part2 (1 vs. 2)

Revision 218 Feb 2016 - ElenaCukanovaite

Line: 1 to 1
 
META TOPICPARENT name="StewartBoogertPhotometry2015"

Testing the fitting

  • First we wanted to test that the function we are fitting is actually normalised.
Line: 72 to 72
 
$\nu_{\rm{tot}}$ 1.71e+06
$\chi^2$ 3222
Degrees of freedom 3591
Added:
>
>

Matching stars between images

  • Used Pythagoras to estimate the closest star image to image. Only used stars whose $\chi^2_{\rm reduced}$ were between 0.5 < $\chi^2_{\rm reduced}$ < 1.5
  \begin{equation} \rm{distance}_{ij} = (x_i - x_j)^2 + (y_i - y_j)^2  \end{equation}
  • Went through images and found the image with highest number of stars.
  • This was the base image that all other images were compared to.
  • It worked fine for images in focus and blue filter since they were taken close enough in time for individual stars to not have moved much.
focus_first.png
  • In red filter however, there was double matching, when more than one star was matched to the same base star as shown below:
red_first.pngred_first_closeup.png

red_first_closeup_two.png

  • Improved the algortihm by:
    • first calculate the minimum distances, some of them will be wronged as stars will be wrongly matched.
    • Use those matches to work out the median x and y shifts.
    • Then check between matches if the x and y shifts match the median x and y shifts to 1 pixel, if they do its a match.
  • Below is the same red image and matches:
red_second.png

red_second_closeup.png

red_second_closeup_two.png

HR diagram for Cr39

  • Using the matching and ML fits we got the following HR diagrams:
hr_b_v.png

hr_v_r.png

  • The stars in V-R HR that seem to have the same V-R are indexed 26 and 27. The star with V close to -8 is indexed 2, it appeared 3 times. Star with V close to -7.5 appeared 3 times.
  • The HR diagrams are terrible.
  • From SIMBAD they should look like this:
B_V_simbad.png

V_R_simbad.png

  -- ElenaCukanovaite - 14 Feb 2016
Line: 124 to 163
 
META FILEATTACHMENT attachment="latexe95444d662b6089e30433993bcc95746.png" attr="h" comment="" date="1455479439" name="latexe95444d662b6089e30433993bcc95746.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexb80262d0b0ab437f7d128be6dc25f7c0.png" attr="h" comment="" date="1455479440" name="latexb80262d0b0ab437f7d128be6dc25f7c0.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex92101250b7e949896dc7365a6d63e4e4.png" attr="h" comment="" date="1455479569" name="latex92101250b7e949896dc7365a6d63e4e4.png" user="zxap014" version="1"
Added:
>
>
META FILEATTACHMENT attachment="red_first.png" attr="" comment="" date="1455825044" name="red_first.png" path="red_first.png" size="525683" user="zxap014" version="1"
META FILEATTACHMENT attachment="red_first_closeup.png" attr="h" comment="" date="1455825265" name="red_first_closeup.png" path="red_first_closeup.png" size="351034" user="zxap014" version="2"
META FILEATTACHMENT attachment="red_first_closeup_two.png" attr="h" comment="" date="1455825198" name="red_first_closeup_two.png" path="red_first_closeup_two.png" size="445456" user="zxap014" version="2"
META FILEATTACHMENT attachment="focus_first.png" attr="h" comment="" date="1455825059" name="focus_first.png" path="focus_first.png" size="530273" user="zxap014" version="1"
META FILEATTACHMENT attachment="red_second.png" attr="" comment="" date="1455825970" name="red_second.png" path="red_second.png" size="98693" user="zxap014" version="2"
META FILEATTACHMENT attachment="red_second_closeup.png" attr="" comment="" date="1455825975" name="red_second_closeup.png" path="red_second_closeup.png" size="40623" user="zxap014" version="2"
META FILEATTACHMENT attachment="red_second_closeup_two.png" attr="h" comment="" date="1455825980" name="red_second_closeup_two.png" path="red_second_closeup_two.png" size="44952" user="zxap014" version="2"
META FILEATTACHMENT attachment="hr_b_v.png" attr="h" comment="" date="1455826422" name="hr_b_v.png" path="hr_b_v.png" size="31012" user="zxap014" version="1"
META FILEATTACHMENT attachment="hr_v_r.png" attr="h" comment="" date="1455826427" name="hr_v_r.png" path="hr_v_r.png" size="31135" user="zxap014" version="1"
META FILEATTACHMENT attachment="B_V_simbad.png" attr="" comment="" date="1455827326" name="B_V_simbad.png" path="B_V_simbad.png" size="50510" user="zxap014" version="1"
META FILEATTACHMENT attachment="V_R_simbad.png" attr="" comment="" date="1455827330" name="V_R_simbad.png" path="V_R_simbad.png" size="34826" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex9bf7a2667af9a16245d4de60644f5686.png" attr="h" comment="" date="1455827418" name="latex9bf7a2667af9a16245d4de60644f5686.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex555fe965e3be97eb5df382bf33d5eb4e.png" attr="h" comment="" date="1455827418" name="latex555fe965e3be97eb5df382bf33d5eb4e.png" user="zxap014" version="1"

Revision 114 Feb 2016 - ElenaCukanovaite

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="StewartBoogertPhotometry2015"

Testing the fitting

  • First we wanted to test that the function we are fitting is actually normalised.
  • The skewed Gaussian function is:
  \begin{equation} f_G(x,y,\mu_x,\mu_y,\sigma_x,\sigma_y,\rho, \beta_1, \beta_2) = \frac{1}{2\pi\sigma_x\sigma_y\sqrt{(1-\rho^2)}} \exp \Bigg\{-\frac{1}{2(1-\rho^2)}\Bigg[\bigg(\frac{x-\mu_x}{\sigma_x}\bigg)^2 + \bigg(\frac{y-\mu_y}{\sigma_y}\bigg)^2 -2\rho\bigg(\frac{x-\mu_x}{\sigma_x}\bigg)\bigg(\frac{y-\mu_y}{\sigma_y}\bigg)\Bigg]\Bigg\} \times \Bigg[1 + \rm{erf} \frac{\beta_x x + \beta_y y}{\sqrt{2}}\Bigg]  \end{equation}
  • Below are some skewed Gaussians for different skews:
summed.png
  • A normalized function would have values between 0 and 1, but here the values are between 0 and 14. This is because the function is summed over all columns and rows.
  • If one looks at an individual column, it can be seen that the function is indeed normalised.
  • Middle column:
003_150.png
  • At the far columns one can see:
001_50.png

005_250.png

  • Also, the formula was extracted from this reference ( doi: 10.1093/biomet/83.4.715 ) which derives a normalized skewed bivariate Gaussian.
  • We wanted to look at a good star and see how the fit looks like in as much detail as possible, as the summation along rows and columns seems to distort the details of fitting.
  • The fit is done by minimizing the $\chi^2$ in the following way:
  \begin{equation} f(x,y,\vec{\theta}) = \frac{\alpha}{2\pi\sigma_x\sigma_y\sqrt{(1-\rho^2)}} \exp \Bigg\{-\frac{1}{2(1-\rho^2)}\Bigg[\bigg(\frac{x-\mu_x}{\sigma_x}\bigg)^2 + \bigg(\frac{y-\mu_y}{\sigma_y}\bigg)^2 -2\rho\bigg(\frac{x-\mu_x}{\sigma_x}\bigg)\bigg(\frac{y-\mu_y}{\sigma_y}\bigg)\Bigg]\Bigg\} \times \Bigg[1 + \rm{erf} \frac{\beta_x x + \beta_y y}{\sqrt{2}}\Bigg] + \frac{1-\alpha}{N*M}  \end{equation}  \begin{equation} \nu_{ij} = \nu_{\rm{tot}} \int \int_{ij} f(x,y,\vec{\theta}) dx dy  \end{equation}  \begin{equation} \chi^2_P = 2 \sum^{N,M}_{i, j=1} \left( n_{ij} \ln \frac{n_{ij}}{\nu_{ij}} + \nu_{ij} - n_{ij} \right)  \end{equation}  <br />
  • So we fit for 9 parameters.
  • We do not integrate, instead we evaluate the function at the pixel and multiply by pixel area.

  • We fitted this star:
complete_fit.png crop.png
  • Column by column the fit looks like this:
column_15.png column_16.png
column_17.png column_18.png
column_19.png column_20.png
column_21.png column_22.png
column_23.png column_24.png
column_25.png column_26.png
column_27.png column_28.png
column_29.png column_30.png
column_31.png column_32.png
column_33.png column_34.png
column_35.png column_36.png
column_37.png column_38.png
column_39.png column_40.png
column_41.png column_42.png
column_43.png column_44.png
column_45.png  
  • The plot parameters are shown below:
$\rho$ 0.0509
$x_0$ 33.5
$y_0$ 33.8
$\sigma_x$ 5.021
$\sigma_y$ 5.373
$\beta_x$ -1.184
$\beta_y$ -1.384
$\alpha$ 0.0014
$\nu_{\rm{tot}}$ 1.71e+06
$\chi^2$ 3222
Degrees of freedom 3591

-- ElenaCukanovaite - 14 Feb 2016

META FILEATTACHMENT attachment="summed.png" attr="" comment="" date="1455477180" name="summed.png" path="summed.png" size="142981" user="zxap014" version="1"
META FILEATTACHMENT attachment="003_150.png" attr="" comment="" date="1455477197" name="003_150.png" path="003_150.png" size="105656" user="zxap014" version="1"
META FILEATTACHMENT attachment="005_250.png" attr="" comment="" date="1455477202" name="005_250.png" path="005_250.png" size="84196" user="zxap014" version="1"
META FILEATTACHMENT attachment="001_50.png" attr="" comment="" date="1455477207" name="001_50.png" path="001_50.png" size="83770" user="zxap014" version="1"
META FILEATTACHMENT attachment="crop.png" attr="" comment="" date="1455477684" name="crop.png" path="crop.png" size="59764" user="zxap014" version="1"
META FILEATTACHMENT attachment="complete_fit.png" attr="" comment="" date="1455477691" name="complete_fit.png" path="complete_fit.png" size="114384" user="zxap014" version="1"
META FILEATTACHMENT attachment="column_15.png" attr="h" comment="" date="1455478119" name="column_15.png" path="column_15.png" size="38923" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_16.png" attr="h" comment="" date="1455478126" name="column_16.png" path="column_16.png" size="41241" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_17.png" attr="h" comment="" date="1455478154" name="column_17.png" path="column_17.png" size="39862" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_18.png" attr="h" comment="" date="1455478158" name="column_18.png" path="column_18.png" size="43716" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_19.png" attr="h" comment="" date="1455478163" name="column_19.png" path="column_19.png" size="42753" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_20.png" attr="h" comment="" date="1455478167" name="column_20.png" path="column_20.png" size="43058" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_21.png" attr="h" comment="" date="1455478171" name="column_21.png" path="column_21.png" size="45588" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_22.png" attr="h" comment="" date="1455478175" name="column_22.png" path="column_22.png" size="45782" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_23.png" attr="h" comment="" date="1455478178" name="column_23.png" path="column_23.png" size="50482" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_24.png" attr="h" comment="" date="1455478183" name="column_24.png" path="column_24.png" size="44552" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_25.png" attr="h" comment="" date="1455478187" name="column_25.png" path="column_25.png" size="47111" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_26.png" attr="h" comment="" date="1455478192" name="column_26.png" path="column_26.png" size="48647" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_27.png" attr="h" comment="" date="1455478196" name="column_27.png" path="column_27.png" size="49285" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_28.png" attr="h" comment="" date="1455478200" name="column_28.png" path="column_28.png" size="52047" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_29.png" attr="h" comment="" date="1455478204" name="column_29.png" path="column_29.png" size="52552" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_30.png" attr="h" comment="" date="1455478208" name="column_30.png" path="column_30.png" size="53859" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_31.png" attr="h" comment="" date="1455478212" name="column_31.png" path="column_31.png" size="53637" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_32.png" attr="h" comment="" date="1455478215" name="column_32.png" path="column_32.png" size="53943" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_33.png" attr="h" comment="" date="1455478220" name="column_33.png" path="column_33.png" size="53177" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_34.png" attr="h" comment="" date="1455478223" name="column_34.png" path="column_34.png" size="51721" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_35.png" attr="h" comment="" date="1455478227" name="column_35.png" path="column_35.png" size="51468" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_36.png" attr="h" comment="" date="1455478230" name="column_36.png" path="column_36.png" size="56885" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_37.png" attr="h" comment="" date="1455478234" name="column_37.png" path="column_37.png" size="45663" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_38.png" attr="h" comment="" date="1455478238" name="column_38.png" path="column_38.png" size="51061" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_39.png" attr="h" comment="" date="1455478242" name="column_39.png" path="column_39.png" size="48392" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_40.png" attr="h" comment="" date="1455478246" name="column_40.png" path="column_40.png" size="44720" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_41.png" attr="h" comment="" date="1455478249" name="column_41.png" path="column_41.png" size="43718" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_42.png" attr="h" comment="" date="1455478254" name="column_42.png" path="column_42.png" size="42595" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_43.png" attr="h" comment="" date="1455478257" name="column_43.png" path="column_43.png" size="40670" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_44.png" attr="h" comment="" date="1455478262" name="column_44.png" path="column_44.png" size="41082" user="zxap014" version="2"
META FILEATTACHMENT attachment="column_45.png" attr="h" comment="" date="1455478266" name="column_45.png" path="column_45.png" size="42211" user="zxap014" version="2"
META FILEATTACHMENT attachment="latexe69822a46019b6a4e4cd8d38f0184b29.png" attr="h" comment="" date="1455478614" name="latexe69822a46019b6a4e4cd8d38f0184b29.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex2f796cc9e36e7a28f035aa71bffa8204.png" attr="h" comment="" date="1455479436" name="latex2f796cc9e36e7a28f035aa71bffa8204.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexc1ee5e4b7bbc4b3de0498005a0d5d5dd.png" attr="h" comment="" date="1455479436" name="latexc1ee5e4b7bbc4b3de0498005a0d5d5dd.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex2c1db3ab700e3511de9f696230f6dca6.png" attr="h" comment="" date="1455479437" name="latex2c1db3ab700e3511de9f696230f6dca6.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex609b1f9aaf260c1506b685af6cec7404.png" attr="h" comment="" date="1455479437" name="latex609b1f9aaf260c1506b685af6cec7404.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexe7b9e3aacf9c90156004f9381aba5cc4.png" attr="h" comment="" date="1455479437" name="latexe7b9e3aacf9c90156004f9381aba5cc4.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex2c74254b08260f2622cce7e27b09e566.png" attr="h" comment="" date="1455479438" name="latex2c74254b08260f2622cce7e27b09e566.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex44fd95b4604601049034012851eedcf0.png" attr="h" comment="" date="1455479438" name="latex44fd95b4604601049034012851eedcf0.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexa023df2e6bfceb5d5e63e9f68f847525.png" attr="h" comment="" date="1455479439" name="latexa023df2e6bfceb5d5e63e9f68f847525.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexe95444d662b6089e30433993bcc95746.png" attr="h" comment="" date="1455479439" name="latexe95444d662b6089e30433993bcc95746.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latexb80262d0b0ab437f7d128be6dc25f7c0.png" attr="h" comment="" date="1455479440" name="latexb80262d0b0ab437f7d128be6dc25f7c0.png" user="zxap014" version="1"
META FILEATTACHMENT attachment="latex92101250b7e949896dc7365a6d63e4e4.png" attr="h" comment="" date="1455479569" name="latex92101250b7e949896dc7365a6d63e4e4.png" user="zxap014" version="1"
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding RHUL Physics Department TWiki? Send feedback