test data - Excel spreadsheet



   Lymm Observatory

     optimal photometry

     how it works

     test data 1

     test data 2

     real data

    Download spreadsheet

    Let's put the following star into the spreadsheet:

        sky = 40

        amplitude = 40

        fwhm = 15

    If you think of pixels 1 unit wide, then the fwhm is 15 pixels, and the central pixel has 40 photons from the sky and about 40 from the star. Alternatively, you could think of this as a sky background of 4000 photons per pixel, with a star fwhm of 1.5 pixels (although the image will then be processed in 0.1-pixel-width rings, which we can't do in practice, of course). The spreadsheet tells you the total flux from this star is 10208 photons, regardless of actual pixel size. Ignore the two comparison stars for now.

    The extraction PSF should have an FWHM equal to the actual star PSF. The camera gain is assumed to be 1 (1 count = photon) and the errors are computed based on photon statistics and pixel sampling.

    Notice the S/N achieved with the different methods. The most precise is obtained with well sampled optimal photometry, followed by 'poorly sampled' optimal photometry (reducing the image to 10 unit wide annuli), and finally aperture photometry. At best, optimal photometry gives about 10% better S/N than aperture photometry, although there is a suggestion optimal extraction suffers with a poorly sampled image.

    Some comments are in order about aperture photometry. First, in this case, the highest S/N is obtained with an aperture that samples about 80% of the total light from the star. There's a simple reason for this - the remainder of the light is thinly dispersed against the sky background, and inclusion of this part of the image actually reduces the overall S/N. Second, in practice, aperture photometry is difficult to do with a small FWHM. I'll explore this further later (when discussing results form synthetic images).

    Now try varying PSF size, whilst keeping the total flux constant, with the following parameters,

        amplitude = 20   fwhm = 20.213

        amplitude = 14.4   fwhm = 25

        amplitude = 10   fwhm = 30.

    Notice how S/N drops as the image is made larger and the starlight diluted against a greater area of noisy sky. There is also a suggestion optimal extraction performs closer to the theoretical ideal with a wider (better sampled) PSF.

    Now, let's consider the two comparison stars. Enter some values for the flux (but brighter than the target star).

    One important issue is what happens if the we can't model the star's PSF exactly - which is why I provide the option to use a different extraction PSF. Try changing the width of the extraction PSF slightly - and notice the following:

    (1) The extracted flux changes (rather similar to changing the aperture size in aperture photometry). (2) The effect is the same for target and comparison stars, so flux ratios (differential magnitudes) are still correct (again like aperture photometry).

    It is useful that the differential magnitudes do not depend critically upon the extraction PSF, since we are unlikely to be able to precisely model the actual PSF anyway!

    Finally, there is the question of which star to use to generate the weight mask. Again there is an analogy with aperture photometry - fainter stars yield highest S/N with smaller apertures than brigher stars, but we must use the same aperture for all stars. With optimal photometry, we use the faintest star of interest to generate the weight mask. You can demonstrate this with the spreadsheet. For example, let the target (and weighting star) have an amplitude of 40, and the comparison stars fluxes of 400 to 4000 units. The spreadsheet calculates the error in the flux ratio (differential magnitude) between the two comparison stars. Now increase the flux of the target star - this error remains more or less constant until the flux reaches that of the fainter comparison star.