-- AaronAndrews - 01 Dec 2016

Observation on 25/11/2016:

Simulation: Simulation_test11.py.txt

  • Can generate background mean (or standard deviation) using taylor series:
\begin{equation} f\left ( x;a_{ij} \right ) = \sum_{i=0}^{N_{i}} \sum_{j=0}^{N_{j}} \frac{a_{ij}}{i!j!} \left ( x-x_{0} \right )^{i} \left ( y-y_{0} \right )^{j} \end{equation}
  • In order to get parameters, set up code to do a fit over an array of values using the expansion, up to a set number of terms.
  • First applied fit directly to image of M103 from observation 1, up to 25 terms:
  • LEFT: original image, M103. RIGHT: taylor series fit (25 terms, origin=(0,0))
  • M103_VIS_60S.pngM103_Simtest_BGmean_3.png
  • Then applied median analysis to frame first, then attempted fit (25 terms):
  • LEFT: Median filter for M103. RIGHT: taylor series fit (25 terms, origin=(0,0))
  • M103_Median_65box.pngM103_Simtest_BGmean_2.png
  • Attempted to fit up to 100 terms - failed.
Truncated Histogram:
  • Aim to make a mode filter over an image by defining a function that returned the mode of a truncated histogram
  • Stars populate the tail of the histogram with low frequencies --> set threshold frequency for mode filter as 1000/((1530*1020)/65^2) = 2
  • The mode is the pixel value with the highest frequency
  • Fit a gaussian to the truncated histogram, the mean of the gaussian corresponds to the mode of the truncated data
  • If there are not enough bins (less than 25), fit may not be executed, therefore we estimate the mode by finding the bin with the maximum frequency, return single value of mode
  • Run mode filter using output from truncation algorithm to obtain an array of modal pixel values
  • M103_Mode_full.png
Attachments Attachments
Topic attachments
I Attachment History Action SizeSorted ascending Date Who Comment
Texttxt M103_Simtest_BG_2.txt r1 manage 0.7 K 02 Dec 2016 - 00:07 AaronAndrews median background fit - parameters
Texttxt M103_Simtest_BG_3.txt r1 manage 0.7 K 02 Dec 2016 - 00:07 AaronAndrews direct background fit - parameters
PNGpng M103_Mode_full.png r1 manage 30.2 K 02 Dec 2016 - 13:10 AaronAndrews M103 sim. - mode for whole frame
Texttxt Simulation_test11.py.txt r1 manage 32.0 K 01 Dec 2016 - 23:05 AaronAndrews Simulation code (11)
PNGpng M103_Median_65box.png r1 manage 92.9 K 02 Dec 2016 - 00:03 AaronAndrews M103 median
PNGpng M103_Simtest_BGmean_2.png r1 manage 100.6 K 01 Dec 2016 - 23:59 AaronAndrews Background mean using median
PNGpng M103_Simtest_BGmean_3.png r1 manage 132.2 K 01 Dec 2016 - 23:59 AaronAndrews Background mean using raw image
PNGpng DoubleCluster_200s_blue_5stack.png r1 manage 1074.2 K 01 Dec 2016 - 09:53 AaronAndrews Double Cluster (blue) stacked
PNGpng M103_VIS_60S.png r1 manage 1533.2 K 02 Dec 2016 - 00:06 AaronAndrews M103 - real image
Unknown file formatfits M103_Simtest_starsBG_3.fits r1 manage 12195.0 K 02 Dec 2016 - 00:03 AaronAndrews M103 simulation plus background directly from real image
Compressed Zip archivezip Obs2_images_stacked.zip r1 manage 20147.7 K 01 Dec 2016 - 09:44 AaronAndrews Observation2 - stacked images
Edit | Attach | Watch | Print version | History: r6 < r5 < r4 < r3 < r2 | Backlinks | Raw View | Raw edit | More topic actions

Physics Webpages • RHUL Webpages • Campus Connect • Royal Holloway, University of London, Egham, Surrey TW20 0EX; Tel/Fax +44 (0)1784 434455/437520

Topic revision: r6 - 02 Dec 2016 - AaronAndrews

 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2023 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