The goal of this tutorial is to setup the SBIG-6 CCD detector on the telescope. This tutorial is split into five parts: setup, connecting, calibrations, observations, and analysis. By the end of the tutorial, you will be familiar with the SBIG camera, how it operates, and how to use it. During these observations, we will characterize the pixel scale of the detector/telescope system by observing several pairs of binary stars and make observations of an open cluster to perform photometric analysis on it in the next tutorial.
Binary stars are all stars with known separation, which we will use to calculate the pixel scale for the system. The pixel scale is the conversion between on-sky separation of two objects and their separtion as they appear in an image. It is typically given in terms of arcsec/pixel.
The following command will display a list of binary stars in a separate browser window:
In [1]:
from astropy.table import Table
stars = Table.read('binarystars.csv')
stars.show_in_browser(jsviewer=True)
Here is a list of open clusters. We suggest selecting M6 during the winter months, M41 during the early summer, and M67 during the late summer months.
If you are only doing the binary star portion of this tutorial, the observations of the open cluster and the second observation of the binary stars can be skipped.
Prior to going to the observatory, there are a few things that will be required in order to successfully complete the tutorial. Please review the manual for the camera prior to using it.
Steps required to check while still in the lab:
If the filters are already installed, please make sure you have the list of filter associated with each position.
Plan your observing run prior to your trip to the telescope.
On the 14 in telescope, the detector has roughly a 18'x13' field of view. So your finding charts should be approxiamtely that size. In a 10 second exposure, the finding chart should have
If you want to perform the analysis while at the observatory, make sure the following software is installed on your computer:
When at the telescope, these are the steps that are required to connect the detector to the telescope. For operation of the telescope, please see the previous tutorial.
In order to process the images in the next tutorial, we will need to take several sets of calibration data. These data correct for systematic effects in the CCD observations. For details about these effects, see 'A User’s Guide to CCD Reductions with IRAF' by Phil Massey. Below we give a short description of each of the calibration images, how to take the images, and what they are correct.
If you are only doing the binary star portion of this tutorial, then this section can be skipped.
Description: Dark frames measure the amount of signal produced by the ambient heat of the CCD.
Process: Take 5 dark frames of 10s and 30s each.
Correction: Each dark frame is median combined into a master dark frame of a certain exposure time. Dark frames will be subtracted from our images. Frames should only be subtracted from frames with the same exposure time. If there is not a dark frame of a given exposure time, then a dark frame can be scaled by the ratio of the exposure times.
Description: Flatfields correct for pixel to pixel variations in the response of the CCD and can also correct for the response across the field of view of the CCD. Flatflieds can be taken either of the sky during twilight or on a flat, white surface that is illuminated by a source with a smooth continuum.
Process: For each filter used, take 5 flat fields. The exposure time should be adjusted such that the peak counts in the image should be ~30000. They should certainly be below the saturation rate of the image (typically 65000.)
Correction: Flat fields are median combined into a master flat field for each filter. Then, flat fields are divided into each of the images with the same filter.
The goal of the observations is to complete two tasks: 1) measure the pixel scale of the detector and 2) produce data to measure the HR diagram for an open cluster.
To measure the pixel scale, we will observe a series of binary stars and calculate the distance in pixels between the two stars and compare that to the known distance between the stars.
In order to measure the HR diagram (as part of the next tutorial), we will need to calibrate our observations and correct for the effect of atmospheric extinction. In order to do this, we will repeat the observations of the binary stars and use those objects as standard stars. By doing a second set of observations and multiple binary stars, we will be able to use the airmass) of different stars to calcute the correction for atmospheric extinction.
For the final part of this tutorial, we will measure the pixel scale of the system. This will require measure the distance in pixels between the two binary stars and comparing that to the known distance between the two stars.
The pixel scale is defined as: $$p = \frac{distance \ in \ arcsec \ between \ two \ stars}{distance \ in \ pixels \ between two \ stars}$$
Steps for mesuring the pixel scale:
A. First copy all of the images into the directory where you are running this notebook.
B. Start up a ds9 instance using imexam
In [21]:
import imexam
viewer=imexam.connect()
C. In ds9, load an image of the binary stars by going to 'File->open' and then select an image.
D. Start the image analysis tool using viewer.imexam()
In [22]:
viewer.imexam()
E. Now hover the circular marker over the image and select one of the two known binary stars. Once the marker is in position, press 'b' to produce an estimate of the centroid of the star. Repeat with the other star. Then press 'q' to quite the program.
F. Now save the x and y centroids calculated for the two stars
In [23]:
xc1=418.795523
yc1=372.080136
xc2=554.336772
yc2=299.007590
G. Calculate the distance between the two stars
In [28]:
d = ((xc1-xc2)**2+(yc1-yc2)**2)**0.5
d
Out[28]:
H. Compare to the known distance of the two stars to calculate the pixel scale.
I. Repeat steps C.-H. for each of the stars in your list. The distance only needs to be measured in one image for each binary star.
J. Create a numpy array with each of your measurements of the pixel scale and calculate the mean and standard deviation for the pixel scale.
In [ ]: