Aperture Photometry Overview - Revision history

Rebull at 02:23, 31 July 2020

2020-07-31T02:23:50Z

Rebull: Rebull moved page I'm ready to go on to a more advanced discussion of photometry to Aperture Photometry Overview

2020-07-31T02:23:33Z

Rebull moved page I'm ready to go on to a more advanced discussion of photometry to Aperture Photometry Overview

Rebull at 18:31, 17 June 2010

2010-06-17T18:31:06Z

Rebull: /* Applying aperture corrections */

2010-06-17T18:29:35Z

Applying aperture corrections

Rebull at 23:36, 10 May 2010

2010-05-10T23:36:43Z

Rebull: /* Aperture Photometry Overview */

2010-05-10T23:35:23Z

Aperture Photometry Overview

Rebull: /* Aperture Photometry Overview */

2010-05-10T23:33:30Z

Aperture Photometry Overview

Rebull at 19:38, 10 May 2010

2010-05-10T19:38:10Z

Rebull at 18:56, 10 May 2010

2010-05-10T18:56:54Z

Rebull at 22:24, 6 May 2010

2010-05-06T22:24:44Z

← Older revision		Revision as of 02:23, 31 July 2020
Line 73:		Line 73:

	==Go on to [[Aperture photometry using APT]].==		==Go on to [[Aperture photometry using APT]].==
−
−
−	~~=PSF photometry=~~
−	To come. Basic steps are essentially the same (detect objects, center up, determine what is background/source, sum up light for the source, check your results). More complicated than aperture photometry because you have many more free parameters than for aperture photometry.

@@ Line 35: / Line 35: @@
 #Check your numbers!
-==Detecting of objects==
+==Detecting objects==
 If you're doing photometry by hand, it's easy to recognize at least the brighter sources just using your own eye.  But, if you're doing photometry automatically (like with MOPEX), you need to teach the computer what to look for. I think it's probably easy to imagine the difference between finding bright single point sources on a flat black background, versus trying to find a wide variety of point source brightnesses, some crowded, inamongst complicated nebulosity.  If that's not clear, imagine just this - you probably can imagine having the computer look for peaks that are at least several times brighter than a dark background, but then what happens when you raise the background level by a constant amount?  The sources are less easily distinguished above the background.  MOPEX has many different parameters (check the online help files) that you can control to affect its ability to detect sources.

@@ Line 53: / Line 53: @@
 ==Applying aperture corrections==
-If you are calibrating your own data, then as long as you use the same aperture and annulus for everything (targets and calibrators), you don't need to apply aperture corrections.  However, for space-based data such as Spitzer (or any data that someone else calibrated for you), you need to be aware of what aperture and annulus were used to calibrate the instrument so that you can scale your results accordingly.  Both [http://ssc.spitzer.caltech.edu/irac/calib/ IRAC] and [http://ssc.spitzer.caltech.edu/mips/apercorr/ MIPS] have tabulated the aperture corrections they derive for the apertures and annuli they suggest you use.  Several suggested aperture/annulus pairs are tabulated on both pages, along with the corrections.  I recommend choosing one set of those predefined parameters, ideally such that the correction you need to use is not particularly large.  Note that you should ''not'' use the same parameters they use for their calibration objects, which are usually much brighter than most of the objects in your target frame.  You should choose one of the other tabulated options.  Note that if you choose to use other parameters entirely (not one of those listed), you will need to interpolate among the tabulated values to find the aperture correction value you should use. In all cases, you need to be aware of the pixel size of your mosaic and how it corresponds to arcseconds.  This information is in the FITS header of any mosaic.  To actually apply the aperture corrections, multiply the flux you have calculated by the aperture correction value.   Neither APT nor MOPEX will do this automatically!
+If you are calibrating your own data, then as long as you use the same aperture and annulus for everything (targets and calibrators), you don't need to apply aperture corrections.  However, for space-based data such as Spitzer (or any data that someone else calibrated for you), you need to be aware of what aperture and annulus were used to calibrate the instrument so that you can scale your results accordingly.  Both [http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/32/ IRAC] and [http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/50/ MIPS] have tabulated the aperture corrections they derive for the apertures and annuli they suggest you use (if those links are broken, go to the SSC website and google for IRAC or MIPS aperture corrections).  Several suggested aperture/annulus pairs are tabulated on both pages, along with the corrections.  Note that one lists NATIVE PIXELS and the other lists ARCSEC, and you should double-check exactly what you're using and what pixel size you have in order to make sure you are using the right thing.  I recommend choosing one set of those predefined parameters, ideally such that the correction you need to use is not particularly large.  Note that you should ''not'' use the same parameters they use for their calibration objects, which are usually much brighter than most of the objects in your target frame.  You should choose one of the other tabulated options.  Note that if you choose to use other parameters entirely (not one of those listed), you will need to interpolate among the tabulated values to find the aperture correction value you should use. In all cases, you need to be aware of the pixel size of your mosaic and how it corresponds to arcseconds.  This information is in the FITS header of any mosaic.  To actually apply the aperture corrections, multiply the flux you have calculated by the aperture correction value.   Neither APT nor MOPEX will do this automatically!
 ==Checking your photometry!==

@@ Line 35: / Line 35: @@
 #Check your numbers!
-==Detection of objects==
+==Detecting of objects==
 If you're doing photometry by hand, it's easy to recognize at least the brighter sources just using your own eye.  But, if you're doing photometry automatically (like with MOPEX), you need to teach the computer what to look for. I think it's probably easy to imagine the difference between finding bright single point sources on a flat black background, versus trying to find a wide variety of point source brightnesses, some crowded, inamongst complicated nebulosity.  If that's not clear, imagine just this - you probably can imagine having the computer look for peaks that are at least several times brighter than a dark background, but then what happens when you raise the background level by a constant amount?  The sources are less easily distinguished above the background.  MOPEX has many different parameters (check the online help files) that you can control to affect its ability to detect sources.
@@ Line 43: / Line 43: @@
 Centroiding is finding the center of the object.  If your aperture is not centered exactly right on the object, you will derive incorrect photometry. (You can empirically discover this in APT by manually changing the center (even the initial center guess) and watching how the photometry changes.)  Usually the computer can centroid to within a small fraction of a pixel.  This is easier for the computer (or a human) when the star is well-sampled, e.g., when many pixels define the point sources.  For images like those for IRAC-1, usually only 2 pixels define a detection, so the results are more sensitive to your ability to centroid.  Neither MOPEX nor APT have many parameters that you can control to affect this algorithm in their respective routines.  APT provides a way to interactively graph the radial profile of the PSF, which will give an indication of how well the centroiding is doing.
-==Determine background==
+==Determining background==
 You need to define an annulus (a ring) that is far enough from the star that you have avoided most of the light in the object, but close enough that the flux you measure is still representative of the true background under your source.  The annulus needs to be thick enough that you have enough pixels included that the noise from the background estimate doesn't adversely affect the total error in the measured flux.  You need to decide if you should take the mean, median, or mode of the pixels in the background, or some other estimator.  Because an annulus might include a nearby source (or cosmic ray), you need to decide how many (if any) pixels should be thrown out of the annulus background calculation, and the criteria for doing so.  APT provides some tools to interactively explore some of these options.  You need to decide how to scale the background you measure to be appropriate for the size of the annulus centered on your source (are the total pixel areas the same in the annulus and aperture?). If you are using ground-based data, you should use the same annulus (and aperture, see next item) for every one of your calibration and target objects.  If you are using Spitzer data, you need to be aware of what annulus (and aperture) were used by the instrument teams to calibrate the instrument so that you can scale your results accordingly.  (See information on aperture corrections below.)  MOPEX has options for all of this (and more); APT has options for most of this.  In all cases, you need to be aware of the pixel size of your mosaic and how it corresponds to arcseconds.  This information is in the FITS header of any mosaic.
@@ Line 51: / Line 51: @@
 You need to define an aperture that is close enough to the star that it includes most of the flux, but not so far away that it adds a lot of sky, because you want to count just the flux from the object, not increase the noise by counting up lots of sky.  For bright objects, the PSF can appear quite large, while for faint sources, the PSF can appear quite small.  However, the PSF is a function of the telescope+CCD+filters, etc., and in reality is ALWAYS the same shape.  You're just seeing different fractions of it for bright and faint sources.  To make things easier (and minimize errors), you should pick the same aperture for all of your data, bright and faint sources.  Thus, you need to pick an aperture (and annulus, see previous item) that works reasonably well for both bright and faint sources.   The graph that describes how the flux changes as your aperture changes is called the "curve of growth", and APT provides an interactive way to explore these settings.   If you are using ground-based data, you should use the same aperture (and annulus, see previous item) for your calibration and target objects.  If you are using Spitzer data, you need to be aware of what aperture (and annulus) were used by the instrument teams to calibrate the instrument so that you can scale your results accordingly.  (See information on aperture corrections below.)  MOPEX has options for all of this (and more); APT has options for most of this.  In all cases, you need to be aware of the pixel size of your mosaic and how it corresponds to arcseconds.  This information is in the FITS header of any mosaic.
-==Apply aperture corrections==
+==Applying aperture corrections==
 If you are calibrating your own data, then as long as you use the same aperture and annulus for everything (targets and calibrators), you don't need to apply aperture corrections.  However, for space-based data such as Spitzer (or any data that someone else calibrated for you), you need to be aware of what aperture and annulus were used to calibrate the instrument so that you can scale your results accordingly.  Both [http://ssc.spitzer.caltech.edu/irac/calib/ IRAC] and [http://ssc.spitzer.caltech.edu/mips/apercorr/ MIPS] have tabulated the aperture corrections they derive for the apertures and annuli they suggest you use.  Several suggested aperture/annulus pairs are tabulated on both pages, along with the corrections.  I recommend choosing one set of those predefined parameters, ideally such that the correction you need to use is not particularly large.  Note that you should ''not'' use the same parameters they use for their calibration objects, which are usually much brighter than most of the objects in your target frame.  You should choose one of the other tabulated options.  Note that if you choose to use other parameters entirely (not one of those listed), you will need to interpolate among the tabulated values to find the aperture correction value you should use. In all cases, you need to be aware of the pixel size of your mosaic and how it corresponds to arcseconds.  This information is in the FITS header of any mosaic.  To actually apply the aperture corrections, multiply the flux you have calculated by the aperture correction value.   Neither APT nor MOPEX will do this automatically!
-==Check your photometry!==
+==Checking your photometry!==
 After you do all of this, the next steps are to check your work. You should check it in at least two ways. Overlay your detections on images, make sure that you are detecting real sources and not noise.  Merge your photometry such that measurements for the same object across many bands are tabulated, and make some color-color plots and color-magnitude plots and see if (a) the plots look like you expect, and (b) if the scatter in various parts of the diagram looks like you expect.  '''Just because the computer  comes up with a number does not mean it is right!!'''  For example, if you plot up [3.6] vs. [3.6]-[4.5] for a galactic observation, the bulk of the objects you detect should be unadorned photospheres, so they should have zero color, with larger scatter at the fainter end.  If the locus of most of your points is not in fact centered on zero, and/or if it's not straight (until saturation at the bright end), then you've done something wrong and you should go back and fiddle with some of your parameters until it is straight and centered on zero, or you understand why it's not.

@@ Line 3: / Line 3: @@
 ''Much of the first several paragraphs on aperture photometry from ground-based telescopes was originally ruthlessly copied without permission from the document entitled "Photometry Using IRAF" by Lisa A. Wells, from 1994, found on the [http://iraf.noao.edu/docs/photom.html IRAF photometry documentation page].''
-There are many techniques involved in doing aperture photometry  and these methods vary from one astronomer to another. Some observers use large apertures for their measurements to account for seeing, tracking, and focus variations (which vary for ground-based observations), while others use small apertures and apply aperture corrections. The sky algorithm used may vary according to the chip characteristics and the data  There are a number of ways to do the standard calibration so be sure to observe standards in a way that is compatible with the calibration package you wish to use.
+There are many techniques involved in doing aperture photometry, and these methods vary from one astronomer to another. Some observers use large apertures for their measurements to account for seeing, tracking, and focus variations (which vary with time for ground-based observations), while others use small apertures and apply aperture corrections. The sky algorithm used may vary according to the chip characteristics and the data. There are a number of ways to do the standard calibration, so if you are planning ground-based observations, be sure to observe standards in a way that is compatible with the calibration package you wish to use.
-Space-based data does not have to worry about variations in seeing and focus (generally), and the sky background is a function primarily of the direction in which you are looking.  Moreover, space-based data usually arrive on your desktop already calibrated.  So, space-based data are much easier to process than ground-based data.  However, you still need to use care in selecting your photometry reduction parameters, because it is very easy to shoot yourself in the foot.
+Space-based data does not have to worry about variations in seeing and focus (generally), and the sky background is a function primarily of wavelength and the direction in which you are looking.  Moreover, space-based data usually arrive on your desktop already calibrated.  So, space-based data are much easier to process than ground-based data.  However, you still need to use care in selecting your photometry reduction parameters, because it is very easy to shoot yourself in the foot.
 Some references on the theory and techniques of aperture photometry are

← Older revision	Revision as of 02:23, 31 July 2020
(No difference)

@@ Line 34: / Line 34: @@
 #Summing up the light in the object (define the size of the aperture to use, subtracting off the background).
 #Apply aperture corrections, if necessary.
 ==Detection of objects==
@@ Line 77: / Line 78: @@
 =PSF photometry=
-To come.
+To come.  Basic steps are essentially the same (detect objects, center up, determine what is background/source, sum up light for the source, check your results).  More complicated than aperture photometry because you have many more free parameters than for aperture photometry.

← Older revision		Revision as of 18:56, 10 May 2010
Line 67:		Line 67:
	[[image:checkphotom.jpg]]		[[image:checkphotom.jpg]]

		+
		+	==Bonus question==
		+	Knowing what you do from the above discussion, why is the computer having trouble finding the fainter sources around these bright sources?
		+
		+	[[image:bright+faintsources.png]]

	==Go on to [[Aperture photometry using APT]].==		==Go on to [[Aperture photometry using APT]].==

@@ Line 1: / Line 1: @@
-=Aperture Photometry=
+=Aperture Photometry Overview=
 ''Much of the first several paragraphs on aperture photometry from ground-based telescopes was originally ruthlessly copied without permission from the document entitled "Photometry Using IRAF" by Lisa A. Wells, from 1994, found on the [http://iraf.noao.edu/docs/photom.html IRAF photometry documentation page].''
@@ Line 26: / Line 26: @@
 If you are working with ground-based data, you need to worry about extinction corrections to the data, so extinction stars need to be observed and reduced along with the object data.  The extinction stars should be observed at airmasses corresponding to the range in airmass of the program objects.  Color and zero point corrections are often applied to the instrumental magnitudes as well to put them on the standard system defined by a set of observed standard stars -- these same standard stars can also be used as the extinction stars.  These stars should be chosen prior to observing so that their colors bracket those of the program objects -- a good rule of thumb is to have at least a 0.5 magnitude range in the colors of the standards to determine reasonable calibrations.
-=Steps=
+=Aperture Photometry Steps=
 '''Once you have calibrated data (ground- or space-based), the basic series of steps for doing aperture photometry are as follows''', with many different options and parameters for each step: