Canopus Userguide
Canopus Userguide
Canopus Userguide
Copyright 1993-2012
Bdw Publishing
Table of Contents
CORE OPERATIONS
2. OPENING IMAGES
11
3. SCALING/INVERTING AN IMAGE
13
15
5. AUTOMATCH
17
6. MANUAL MATCHING
23
7. BLINKING IMAGES
27
29
35
ASTROMETRY IN CANOPUS
37
1. MEASURING AN ASTEROID - I
39
2. MEASURING AN ASTEROID II
45
49
51
PHOTOMETRY IN CANOPUS
55
1. A FIRST LIGHTCURVE
61
95
Table of Contents
USING PHOTORED
99
103
109
111
2. TRANSFORMS ERRORS
125
127
137
139
147
155
157
159
165
181
185
187
189
193
4. ALCDEF EXPORT
197
207
211
215
Table of Contents
8. ALCDEF IMPORT
SUPPLEMENTAL TUTORIALS
221
225
227
245
253
4. EXOPLANET LIGHTCURVE
265
271
277
281
283
289
Table of Contents
Introduction
The Canopus/PhotoRed Reference Manual are combined into a single Adobe PDF
document (Canopus_PhotoRedV10.PDF). By default, it is installed in \MPO\DOCS.
If you installed the MPO software to a different base directory, look for the DOCS
directory in that location.
Introduction
to learning even after years of practice and still then, using a little black magic to make
the final determination.
It is strongly recommended that you install the example images since they are reference throughout the Users Guide. Attempting to use your own images when first
leaning may mean some things are lost in translation and the learning curve will be
a vertical cliff.
The images for lightcurve determinations were dark framed but not flat-fielded. The
absence of flat-fields had no adverse affect on the determination of a lightcurve. In fact, a
very clean, smooth curve was obtained. This will not always be the case and the use of
flat-fields for all work is strongly encouraged.
After installation (assuming you selected the default destination of \MPO), the astrometry
images are located in
\MPO\EXAMPLES\ASTMTRY
The lightcurve photometry examples are in
\MPO\EXAMPLES\LTCURVES
There are subdirectories below these with images from separate nights.
The photometry reduction examples are in several directories under
\MPO\EXAMPLES
The Moving Object Search images are in
\MPO\EXAMPLES\MOVINGOBJECTSEARCH
The Variable Star Search images are in
\MPO\EXAMPLES\VARIABLESTARSEARCH
All the images are real world, i.e., they were measured in Canopus and/or PhotoRed
and the results, when appropriate, submitted. In the case of the astrometry images, the
results went to the Minor Planet Center. The photometry results were published in the
Minor Planet Bulletin. At the Palmer Divide Observatory, more than 600 lightcurves
have been measured and then published in the MPB using Canopus. Other observers using Canopus have contributed hundreds of curves and the number is growing with each
quarterly publication. In other words, it can be done. All you need is a little practice,
patience, and, of course, Clear Skies!
Introduction
Meaning
<Key> <button>
Click
Right-click
Shift+click
Press and hold down either shift key and click the left
mouse button. Release the shift key after the click.
<Ctrl+Key>
Press and hold down either Ctrl key and press then release
the indicated key or mouse button. Then release the Ctrl
key. For example, <Ctrl+3> means to hold down the Ctrl
key and press the 3 key on the main keyboard (not number pad)
If the key is a capital letter, the lower case version is OK,
i.e., dont hold down the Shift key unless its specifically
part of the keystroke command. For example, if told to
press <Ctrl+Alt+I>, dont hold down or even press the
Shift key as part of the keystroke combination.
Drag
Depress the left mouse button and, while keeping it depressed, move the mouse cursor to a specified location,
e.g., to a list in a box.
In some cases, you depress the right mouse button. This is
not common and will be specifically mentioned. If no
mouse button is specified, use the left button.
Some mice have features that allow dragging by pressing
a single or special button. Consult the mouses user manual or on-line help.
When asked to enter a string, e.g., Examples, enter the value in the appropriate field
without the quotes unless specifically told to do so.
Mac Users
Canopus and PhotoRed are Windows/Intel (Wintel) programs. They will run on a Mac
that is running a virtual machine that loads the Windows operating system, e.g., Parallels.
There is no provision made for Mac-specific needs or requirements and it is presumed
that you are familiar with running programs under Windows.
You must have files on a local drive or, if on a network drive, that drive must be
mapped to a local drive letter. You cannot, for example, load an image via a URL
(\\remotemache\remotedirectory\Afile.FIT). While the image will load, the AutoMatch routine will likely fail with a 120 error.
Inverted Images
The charts and images in the tutorials are often inverted (dark stars, light sky) even
though this may not be the normal appearance in the program. This is done to make see-
Introduction
ing fine details a little easier as well as to improve printing quality. Large areas of black
do not always print well, especially if there are fine details within the area, such as stars.
Instead of noting every instance where the chart and/or image has been inverted for this
Users Guide, you should presume that a screen shot (or part of it) has been inverted when
you see, for example, dark (black) stars on a light (white) sky. Do not be concerned that
the screen shot is not an exact match for what you see if the only difference is white on
black versus black on white.
Core Operations
Core Operations
The following lessons cover things that you will use at one time or another regardless of
the task you are running. These include loading and scaling images, setting the configuration, adjusting the apertures for astrometry and photometry (the areas used to measure the
target and sky background), and so on. You should be thoroughly familiar with these
concepts before trying to move on to a specific task.
As noted in the introduction, none of these tutorials contain blow-by-blow instructions
to the level of how to click on a star in an image. That is a basic Windows function that
you are presumed to know.
It will be helpful to follow the tutorials in this section in order. In some cases, a
given lesson presumes that youre familiar with information from a previous tutorial.
Core Operations
1.
Several of the settings in the configuration form are important for making astrometry and
photometry easier and accurate. In some cases, it means the difference between getting
and not getting the program to match the image to a set of reference stars, whats called
an AutoMatch.
See the Configuration section in the Canopus Reference Manual for details on the individual settings. It is important that you understand each value, how it is used, and in
some cases how it interacts with other values in ways that allow the programs to work
as intended or prevent them from doing so.
The configuration form allows you to save multiple configurations under different
names. The settings are stored in the Windows registry. Once created and saved,
you can quickly recall the settings for a given telescope/camera/preferences combination by selecting the previously saved profile.
1.
Open the Configuration form by clicking the configuration speed button on the
top tool bar, using the keyboard combination of <Shift+Ctrl+S>, or selecting
File | Configuration from the main menu.
2.
The Configuration form has several tabs, each with a specific purpose. Go to the
General page by clicking on its tab.
3.
Enter a profile name in the Profile field or select one from the drop down list.
4.
Enter your longitude, latitude, and elevation (in meters). 1 foot = 0.3048 meters.
5.
This is not the difference between your computer clock and Universal Time. It is the
difference between the time in the FITS headers and UT. Under most circumstances,
this value is 00:00:00, meaning the time obtained from the FITS header is in UT.
6.
Enter the focal length. Use the radio buttons immediately to the right of this field
to indicate if the value is inches or millimeters.
7.
Enter the eADU value. This should be a positive, non-zero number that is the
conversion factor from electrons to ADU as given by the camera manufacturer.
Enter 2.3 if you dont know the value.
8.
Enter the number of columns (Cols) and rows (Row) in the images. The program
usually reads these from the FITS header but it is a good idea to enter the correct
values regardless.
9.
Enter the pixel sizes in microns or arcseconds. Use the radio button group to specify the units. See the Canopus Reference Manual for additional information.
10.
Select the appropriate Header Exposure Time. This is usually Start, meaning
that the time taken from the FITS header is the exposure start time.
11.
If you have trouble getting the programs to AutoMatch, check Ignore OBJECT
KW and/or Ignore FL/Pix Size. If you do this, the values for focal length and
numbers for Col, Row, and pixel sizes must be accurate.
2.
Check the Reverse E/W and Reverse N/S boxes and enter an appropriate
value in the Rotate field so that a chart generated by the program matches the
approximate scale and orientation of your images. Usually, all thats needed is to
set the two check boxes as needed. However, if using an AO or other device that
affects the image in one axis but not the other, it may be necessary to set the chart
rotation as well.
3.
The Max Scale Diff setting may affect if you can get an AutoMatch. After an
initial AutoMatch, Canopus and PhotoRed check the original scale of the chart (in
degrees per side) versus the value after the match and find the ratio of the smaller
value to the largest (meaning the ratio is always 1.0 or greater). If the ratio exceeds this value, then the AutoMatch is forced to fail.
If the configuration is set correctly and/or Canopus can get the needed values
form the FITS header, then the scale ratio should be in the range of 1.0 to 1.5. Its
recommended that you dont go less than 2.0, which would account for a difference of having an image binned at 1x1 but the configuration set to pixel sizes that
are equal to 2x2 binning.
This Max Scale Diff setting is a safety valve in that a very bad mismatch between the derived and true plate scales can make the program appear to hang (or
actually do so). For example, the original scale might be 0.3 degrees per side but
the derived scale was 12 degrees per side. This could mean tens of thousands of
stars are being added to the chart and Canopus has to work its way through all of
them. That can take a long time.
1.
This setting is important when using TrueMags for analysis (these will be covered in
later tutorials). You want to keep the overall standard deviation as low as possible
to reduce the scatter in the derived data.
2.
Check the Ignore Saturation box if you want Canopus to measure and use stars
that are above the Saturation value. This is not recommended unless you have
good reason and are aware of the potential consequences. The program cannot
find an accurate centroid or magnitude for a star that is in the non-linear response
of the camera or full saturated.
About Saturation
Having the right saturation value is critical to some operations in Canopus, e.g., when
doing photometry using magnitudes based on the Magnitude/Intensity Relationship
(M/IR, see the Reference Manual and tutorial in the Photometry chapter). If you set the
value too high, so that stars that are in the non-linear response of your camera, or are
actually saturated, then the MI/R solution can be off significantly and your results badly
skewed. When in doubt, underestimate the value and never assume that it is the maximum allowed by your cameras ADU converter. Just because the ADU is 16-bit doesnt
mean that the camera is linear all the way up to 65535 ADU counts.
Note that the preceding assumes 16-bit images, usually the native output of cameras.
If youre using processed images, they may be 32-bit (integer or float) or 64-bit
(float) and so the saturation number will be on the order of 2 billion for 32-bit integer images. If youre working with 32 or 64-bit images, check Ignore Saturation.
The odds of any pixel actually reaching the maximum value for these images are extremely remote unless the images are the result of stacking several 16-bit images
and a given star was non-linear or saturated on one or more of those images.
Saving Changes
Click <OK> to save the current settings.
Click <Cancel> to close without saving changes and revert to the settings before the form
was opened.
Click <To Text> to save the settings a text file. This is handy if you need tech support.
Some of the tutorial sets, e.g., astrometry and lightcurve photometry, include short
tutorials that set certain fields to specific values. What was covered here are the
core values that MUST be set or at least considered in order for the programs to
work properly.
10
2.
Opening Images
Canopus can load FITS (Flexible Image Transport System), SBIG, BMP, and JPG files.
However, FITS and SBIG are the only two that allow you to AutoMatch an image and for
which you should use for science work with FITS the standard format.
It should be noted now that there are many FITS standards. MPO software does
its best to adhere to those whose definition can be found at
http://fits.gsfc.nasa.gov/fits_home.html.
The Canopus/PhotoRed manual has an appendix of recognized keywords, i.e.,
names for values that are found in FITS/SBIG headers. Many programs add keywords to the header to define other values or use standard values in a way not
expected. Canopus/PhotoRed may or may not recognize them or use them properly.
The only guarantee is that every effort will be made to have the MPO programs
work within the standards defined on the above link and the keywords listed in the
Reference Manual appendix.
Select Image | Open image from the main menu or press (<Ctrl+O>, the letter,
not the number). This displays a Windows file open dialog where you can select
the image to be loaded.
A list of the five most recently opened images is found on the Image menu. Use
this to reload one of those images without having to navigate through the Windows file open dialog.
2.
Select Image from the main menu and chose a previously opened image from the
MRU (most recently used) file list.
11
12
3.
Scaling/Inverting an Image
Canopus and PhotoRed provide some basic image manipulation to make it easier to see
faint targets or bring out details in other ways.
The scaling does NOT change the original data. It simply modifies the way the image data are displayed. You can save the modified display as a BMP or PNG but the
underlying FITS data does not change.
1.
2.
Right-click over the image and choose Image scaling from the pop-up menu to
display the Image Scaling form.
3.
Select a method of scaling and adjust the values to get the desired effect.
4.
2.
Right-click over the image and select Image scaling from the pop-up menu to
display the scaling form.
3.
4.
5.
13
All images opened after setting a default use those settings. Some images may not
scale well with the default values, showing as all black, all white, or with extremely
high contrast. You should use custom scalings carefully and remember that they
could make some images worse than better.
Right-click the over the image and select Invert from the pop-up menu.
or
2.
Each time you select Invert, the display is set to the opposite state.
If an image is open, it is set to dark stars, light sky. If it is already set to inverted, it
does not switch to normal mode.
Aperture Colors
Even when an image is inverted, the aperture color settings are retained. See the tutorial
on setting aperture sizes and colors (pg. 15).
14
4.
The measuring apertures are used to define the areas on the image used to find the brightness of a star, asteroid, or other object what will be called the target and the brightness of the nearby sky. There three apertures, centered on a common point and concentric
with one another. The screen shot below shows a sample set of apertures.
The inner-most circle is the target aperture. Pixel values from within this region define
the brightness of the target.
Just outside that is another circle that forms an annulus with the target aperture. The area
of this annulus is called the dead zone. Pixels within this region are ignored entirely.
Its a buffer zone between the region that defines the target and the region that defines the
sky background so that pixels within the annulus are not counted twice.
The outer-most circle defines the outer boundary of the sky annulus with the dead zone
aperture defining the inner boundary. The pixels within this area are used to compute the
sky background. The average value is then subtracted from the value of every pixel in the
target aperture so that what remains is the result of the only the target and not the target +
sky.
Click <Apertures> on the main form tool bar to display the apertures form.
2.
15
The height and width values can be different but both must be odd numbers, e.g., 3,
7, etc.
You can define aperture sets for three types of objects: the target (e.g., an asteroid
or variable star being measured for time-series work), the comparison stars (when
using the lightcurve wizard and transform wizards), and astrometry (when measuring field stars during an AutoMatch).
2.
The Width and Height settings for a given object type are for the target aperture, the inner-most aperture. The Dead and Sky settings are for the width
of the given zone. For example, assume the target aperture is 11x11 pixels. A
Dead setting of 2 means the dead zone is two pixels wide and so has a total
width of 13 and radius of 2 + 5, or 7 pixels. Keep in mind that because of the
forced odd number of pixels in the target aperture of 11 pixels there are five pixels on either side of the center column (or above/below the center row). A
Sky setting of 11 means the outer aperture has a radius (assuming circular apertures) of 11 + 2 + 5, or 18 pixels.
3.
You can set the target apertures independently of the comparisons apertures. This
is useful if working a fast moving object and so the target and/or the comparisons
are trailed.
4.
The apertures for a given object type do not have to be symmetrical, i.e., round or
square. They can have different heights and widths. All the apertures (target, dead
zone, sky) for a given object type are the same shape since the latter two are defined in terms of the target aperture.
5.
You can set the rotation angle of a non-symmetrical aperture set. For example, if
using elliptical apertures, where the height and width are not the same, you can
rotate the ellipse to match the trailed image of a fast moving object.
6.
You can use different colors for each object type, e.g., yellow for the target and
red for comparisons.
When using the Lightcurve and Photometry wizards in Canopus and PhotoRed
where you are measuring a series of images in which an anchor star is used to
find all other objects, the anchor star (one of the Comparisons group) always
uses bright green (lime) apertures. You should select a color for the Comparisons
other than lime and that used for the Target type so that you can easily spot the
anchor star.
7.
You can copy the Target type aperture set to the Comparison and Astrometry types by pressing <F9> or clicking <Copy Target>.
See the set of Astrometry and Photometry tutorials for examples of setting the apertures.
16
5.
AutoMatch
One of the most common and important operations when using Canopus and PhotoRed is
AutoMatch. This is where the program reads the information from the FITS/SBIG header
and/or configuration settings, generates a star chart centered on the presumed center of
the image, and matches the charts center, rotation, and scaling to the image. In more
specific terms, the chart is created so that for any given star in the image, the corresponding star in the chart is at the same distance and angle from the chart center. Mathematically, the star has the same polar coordinates in the image and chart, both of which have
the same origin point in the sky.
When this is done, so-called plate constants are found. These allow converting an X/Y
coordinate on the image to an RA and Declination, thus allowing, for example, Canopus
to determine the coordinates of an asteroid to be sent to the Minor Planet Center. Also
found during an AutoMatch is whats called the Magnitude-Intensity Relationship
(M/IR). This converts the measured total ADU count for an object in the image to a magnitude based on those taken from the one or more catalogs used to generate the chart.
AutoMatch Fundamentals
The Canopus Reference Manual contains details on the AutoMatch process such as how
the programs read the header to get the necessary information and when/if those values
are superseded by the configuration settings. On the assumption that the necessary values
are available, it is important that you have a basic idea of how the AutoMatch works. For
the sake of brevity, the following steps assume Canopus is being used. The core AutoMatch code is used in all MPO programs and so a general description applies to all.
1.
Canopus reads the FITS header and determines the approximate RA and Declination of image center as well as the scale (degrees/pixel) and orientation (e.g.,
north up, east left) for the chart.
2.
3.
The chart is generated and the star data sorted to find a fixed number of the
brightest stars (the charts constellation).
4.
Canopus forms a number of two-star pairs using the image constellation and compares each pair to every possible two-star combination among the chart constellation.
During the AutoMatch process, you will see a set of colored lines on the chart.
These represent the pairs that Canopus uses to try to find a match.
17
For each comparison, the chart data is translated, rotated, and scaled so that the
separation and angle of the test pair from the chart matches the pair from the image.
6.
Canopus then checks how many matches can be made between the new X/Y coordinates of the chart constellation to the image constellation.
7.
The chart pair that results in the most number of matches to the image is used redraw the entire chart and attempt to find an astrometric solution that has an acceptable error. That error is found by back-fitting the derived plate constants to
find the RA/Declination of the chart stars and comparing them to their catalog positions. The mean error and standard deviation of the comparisons are used to determine if the solution is acceptable.
8.
If Canopus can read the FITS header such that it knows that it was approximately
centered on a given asteroid, it will use the plate constants to find the asteroid on
the image and automatically measure it (AutoMeasure). This simulates the actions
of clicking on the asteroid in the image, which derives the RA/Declination from
the plate constants, right clicking over the image and selecting Set Asteroid Position from the pop-up menu.
AutoMatch Tutorial
The following steps will take you through an AutoMatch. These general steps apply to
any image. Where the details may change is the configuration and aperture settings.
1.
18
Open the Apertures form and duplicate the settings in the screen shot below.
Open the Configuration form, duplicate the settings shown below, and save the
changes. The name you give the profile is not important.
The presumption is that you installed the MPOSC3 catalog but not UCAC2 or
UCAC3 and that you do not have the 10/11 CD set of the USNO-A catalog. Even
if you did, you should generally use only one catalog at a time, assuming that
there are enough stars in the catalog for the field of view being covered.
When doing astrometry, the UCAC2 or 3 catalog is the one of choice. For photometry in Canopus, use MPOSC3 since it contains data from the CMC-14 and SDSS
catalogs. Where available, the Sloan r magnitudes are much more internally consistent over the BVRI magnitudes derived from a 2MASS to BVRI conversion (Warner, 2007; Minor Planet Bulletin 34, 113-119).
19
Note that the Sloan r magnitudes were chosen on the Photometry tab for this
exercise. In this case, it is OK because there are a sufficient number of stars in the
field with r magnitudes. This fact will be used to demonstrate an important feature
in Canopus for doing lightcurve photometry. In other cases, you may want to or
have to choose another photometric band.
3.
4.
5.
You will see a set of colored lines on the chart. After a few seconds, those go
away and are replaced by numbers. On the image, a set of measuring apertures are
displayed. One set, those about the asteroid, is in yellow.
6.
Look at the upper right of the Canopus form. Youll see the summary of the M/IR
solution in a tool bar giving the number of the stars in the solution, the M/IR offset, and the standard deviation of that offset. See the Reference Manual for additional information.
The Refs number may not match the number of stars used in the astrometric solution, which determines the plate constants that allow computing the RA and Declination of a star based on its X/Y coordinates on the image.
7.
20
Click <Reductions> on the tool bar at the top of the Canopus form (it has an
icon that resembles a spread sheet or table).
8.
Press <Ctrl+3>
This page shows the derived position for the asteroid at the bottom of the page.
The table shows the data for the reference stars used to find the plate constants.
The numbers in the far left column of the table correspond to the numbers on the
chart back on the Measurements page. The RARes/DCRes values are the differences between the catalog position and the one computed from the plate constants
and are given in arcseconds. The standard deviations of these residuals are shown
in the S.D. fields at the bottom of the form. In this case, the match is very good,
with the deviations ~0.15 arcseconds.
As youll see in a later tutorial, the M-C values are the differences between the
measured and calculated position. The absolute value of these numbers is not as
important as the trend of the numbers when comparing results from a number of
images taken on the same night.
AutoMatch in PhotoRed
In PhotoRed, you go through the same basic steps as above. The difference is that after an
AutoMatch, there is no Reductions page and you wont see the apertures on the image.
Instead, a simple message appears that says AutoMatch Successful.
21
22
6.
Manual Matching
Sometimes Canopus cannot perform an AutoMatch. This can for several reasons, usually
because it cannot interpret information from the FITS header correctly and so cannot
determine the approximate center of the image or approximate scale and/or orientation
for drawing the chart. When you encounter such a case, you can still get Canopus to
match the image and so get the plate constants and M/IR solutions needed for other
things.
The following steps refer to these two screen shots. The top one is the CCD image while
the bottom is a chart drawn by Canopus that approximately matches the image.
To do a manual match requires that you first have Canopus draw a chart that approximately matches the scale and orientation of the image. From there, you will do manually
what Canopus does automatically, indicate a pair of stars on the chart that correspond to
the same stars on the image. From there, Canopus can do the rest.
23
After you open an image, right click over it and select one of the three Zoom
items from the popup menu. Each one provides a different level of magnification,
200%, 300%, or 400%. You can also display the form by pressing <Ctrl+Shift+2>
(200%), <Ctrl+Shift+3> (300%), or <Ctrl+Shift+4> (400%).
2.
Move the box to a convenient location on the screen so that it doesnt get in the
way of doing things.
Note that the box also shows the current set of apertures. If you change the aperture settings (sizes), the box automatically updates. The apertures are always red.
If you open the configuration form, the box is closed and you have to open the box
again.
2.
Select Image | Generate manual chart from the main menu. This displays a message.
3.
Click <No>.
24
4.
If you do an AutoMatch, it generates a file that has data on the stars in the image
that was measured. You may want to generate the chart for that image again but not
do an AutoMatch. Also, that file might still be needed for some other functions. By
answering No, that file is not deleted. Of course, the presumption is that youll be
working with the same image from which the data file was created. If you change
images, then that file is no longer valid and you should delete it.
After you click < No>, the chart generation form appears.
If possible, Canopus fills in the fields for you. In this case, it was able to determine that the image was of 771 Libera and set the form accordingly.
5.
Click <OK> to generate the chart. Your screen should look similar to this:
6.
Using the reference images at the beginning of this tutorial, locate star A on the
image and click on it. You should see a set of measuring apertures centered on it.
7.
Right-click over the image and select Select Image Star | Add Star 1 from the
popup menu.
25
Find star A on the chart. Center the mouse cursor over the star and right-click.
9.
Select Set chart star | Add Star 1 from the popup menu.
10.
Repeat these steps for Star B, except that you select Add Star 2 from the two
popup menus.
When selecting the two stars to be used in this process they should be 1) as far apart
as possible but not too close to the edge and 2) not be on a vertical or horizontal
line, i.e., a line joining the two would not be nearly parallel to the vertical or horizontal edges of the image. Following these guidelines improves the accuracy of the
match.
11.
Right-click over the chart and select Match to Image | Auto measure from the
popup menu.
After a few moments, you should see something similar to what followed an AutoMatch
in that tutorial: a number of aperture sets about stars on the image and reference numbers
next to a number of stars on the chart.
What you see after this match may not match the screen shot of the AutoMatch exactly,
i.e., a different set of stars may have been used in the final solution. That could be due to
the AutoMatch finding a different star pair that had a larger separation (or angle) that
produced a solution with smaller residuals. The main point is that you do see the apertures and reference numbers and when you go to the Reductions page that you see something similar to what followed in the AutoMatch tutorial.
26
7.
Blinking Images
There may be times when you want to compare two or more images to see if there is an
unexpected moving object in your image or to check a reference image of a galaxy
against a more recent image to see if a supernova may have occurred. The Moving Object
Search (MOS) utility in Canopus can search for moving objects in a number of sets of
images but that is overkill for a simple check to make sure you know where an asteroid is
or to locate unexpected targets. This is where the Blinker in Canopus comes in handy.
The idea behind the blinker is simple: two or more images are aligned and then each one
is shown in quick succession. All non-moving objects that are in every image appear to
stand still while a moving object changes location from image to image. It was a similar
concept that was used by Clyde Tombaugh to discover Pluto, but he was using large photographic plates and a mechanical device.
2.
Click <Open> and load a file. For this tutorial, locate and open
\MPO\EXAMPLES\LTCURVES\A771\SEP18\A771A001.FIT
3.
The image windows are cascaded so that you can see all (or most) of one image and the
caption bars of the other images underneath.
You can load more than one image at a time. When the file open dialog appears
after clicking <Open>, select as many files as you want (be reasonable) using the
standard Windows techniques for selecting multiple files.
4.
On the first image, click on a brighter star near the middle of the image.
5.
Click the caption bar of the second image to bring it to the front and then click on
the same star on that image.
6.
If you want to view the images inverted (dark stars and light sky), which can help
finding fainter objects, click <Invert> on the Blinker page (its immediately to the
left of the Clear button).
7.
Click <Start>. This minimizes the two images and presents a new window where
the images are alternately displayed. The caption on the button changes to
Pause.
27
Unless you picked the asteroid (its happened), all the stars in the image should be stationary while one star the asteroid appears to move back and forth.
You can change the rate at which the two images are displayed by changing the value in
the Rate (ms) field by entering a value directly or using the spinner control.
If you click <Pause>, the blinker stops on the current image (and the caption reverts to
Start.) Use the radio buttons at the top of the form to select which image is displayed.
If one or fewer images are open, the Start/Pause button is disabled.
Use <Clear> to close all images. The button is disabled while the blinker is operating.
You can close a window and open one or more others to use for blinking. Make sure that
you click on the same reference star in each image before starting the blinker. If not, the
images are not aligned and you could get dizzy as the blinker displays those misaligned
images.
28
8.
Saving Plots
All plots can be saved using a menu option or button provided for that purpose. The default is to save as a PNG since these are much smaller than BMP files and have very little
compression loss. When a file dialog appears, select the type of file before saving. Regardless of the extension you give the file, MPO software forces the extension to the
selected file type.
Here is a screen shot of Canopus after generating the lightcurve of an asteroid. The box
has been added to indicate the approximate region on which we want to zoom.
Zooming In
1.
Move the mouse cursor a little above and to the left of the region that you want to
see zoomed.
2.
Drag the mouse cursor a little below and to the right of the region and then release
the mouse button.
29
Now you see the region filling the entire plot area.
You can zoom in as much as you want. If you want to zoom in some more, just repeat the
process of selecting a region but on the zoomed plot. It is not necessary to reset the plot
to zoon in more.
Zooming Out
1.
To zoom back to normal view, reverse the process, i.e., place the mouse cursor
anywhere on the plot area, drag up and to the left a few pixels, and release the
mouse. The plot returns to normal scaling.
Panning
Panning means to move the entire set of data in any direction. You can pan the plot
regardless of the zoom level.
1.
Move the cursor anywhere on the plot and drag using the right mouse button instead of the left.
This screen shot shows the original plot after panning to the right.
30
The most common use of panning it to move about the data set after zooming in.
Reset Panning
Use the same procedure to reset the plot as you did for zooming out, i.e., drag up and to
the left from anywhere on the plot and then release the left mouse button.
Floating Range
Check this radio button to have Canopus scale the vertical axis such that the curve is
about 70-80% of the plot area height, but with no less than ~0.12 mag total height.
Fixed Range
Check this radio button and enter the vertical range of the plot. If the amplitude of the
curve exceeds the range, Canopus automatically switches internally to use the floating
range option but the setting does not change on the tool bar.
Click <Replot>, immediately to the right of the period slider bar, to replot the data if you
make any changes to the range settings.
Move the slider bar left or right of center. As you do, youll see the new period
value and the amount of change in the labels below the slider.
2.
This is a handy tool for getting an estimate of the true period error. The one reported by
Canopus is the formal error, which should often be increased to provide a more reasonable error estimate. If you change the period by 0.01 h and the plot falls apart or the
31
Click on the sliders thumb (button) and press one of the number keys on the
main keyboard (not the number pad). This resets the slider to the middle and
changes the total range.
Key
1
2
3
4
5
6
7
8
9
0
Range (hours)
5.000
1.000
0.500
0.200
0.100
0.067
0.050
0.033
0.025
Reset to default (0.100)
Clicking <Find> resets the slider range and position to its defaults.
Plotting Options
Click <Options> (immediately below the <Replot> button) to see a menu with several
options. Most of the menu options are toggles, meaning if you select one, it goes from off
to on or vice versa.
Zoom
This is a special implementation of zooming the lightcurve.
Use the slider bars to change the vertical and or horizontal scaling of the plot. If the
Track together box is checked, the sliders move in unison. Uncheck the box to change
the scaling independently.
You can save the rescaled plot as long as you dont close the zoom form or click
<Reset>.
Show Legend
This menu item toggles the legend on the lightcurve plot. When publishing the lightcurve, you should include the legend. Use a graphics program to move the legend within
the data area. This makes for a much better presentation.
32
Flip Y-axis
Lightcurves should be presented so that the object is brightest when the data points are at
the top. If youve been given data where the asteroid is brightest at the bottom, use this
option to present the lightcurve correctly.
Show Errors
This toggles the error bars on the data points in the plots.
Save LC Plot
Click this to display a file dialog. Plots can be saved as PNG or BMP files.
The Reference Manual covers the details of the plotting options controlled on the Configuration | Photometry | Plotting Options tab.
One worth mentioning here is the ability to generate custom-sized plots. This is useful for
when publishing your lightcurves in a journal such as the Minor Planet Bulletin, Journal
of the AAVSO (JAAVSO), and so on. Another is to plot the data showing more than one
cycle (the data are repeated when going past the first cycle). The MPB Phased Plots
forces the phase range from -0.05 to +1.05 cycles. Some journals want a wider range. If
you uncheck the MPB Phased Plots box, you can set the range from 1.0 2.0 cycles.
These custom-phased plots always start at 0.0 phase.
33
34
9.
The MPO Server is a Windows program that runs automatically whenever you start an
MPO program (MPO Canopus, Connections, MPO Viewing Guide, or MPO LCInvert).
The server is one side of whats called a client-server database system. The MPO programs are the clients and, of course, the MPO Server is the server.
The idea is that the client applications never access data file directly. Instead, they send a
request to the server, which then accesses the files, and returns the requested data. When
editing or adding data, the server is the one that does the actual modifying, using the data
supplied from the client. Think if the server as a bank teller who takes and gives you
money from the banks vault (Im showing age by assuming that not all people use electronic banking).
You should almost never have to worry about the MPO server application. However,
there may come a time when you need to start and restart it. For example, if one of the
MPO client programs shuts down unexpectedly or in a way not intended, the server may
keep a hold on one or more files such that when the client application tries to access the
file again, the server says that its not available. The teller is waiting to get the managers
approval to release your funds.
You will get error messages on occasion that tell you specifically to stop and restart the
MPO server or that a certain action cannot be done. The latter is usually when trying to
make a backup of the main files as Canopus closes. Should you suspect that the teller is
being stingy, its easy to get the managers approval and conduct your transactions.
1.
2.
On the far right of the Windows task bar, usually at the bottom of your screen, is
the icon tray..
This holds the icons of programs that are running in the background, e.g., a virus
scanner. Among the icons you should see is one that looks like an Earth globe
with FF superimposed. A white box has been added to the screen shot above to
indicate the MPO Server icon.
3.
Note the green LED at the bottom right of the server form.
Do not change any options or settings on the MPO Server unless instructed in this
Guide or by technical support.
35
Click the black square at the upper left. This stops the MPO Server.
5.
Wait until the black square becomes disabled (turns gray) and the button to its left
displays a green, right-pointing arrow. The LED should now be red.
It can take ten or more seconds for the MPO Server to stop.
6.
Click the green arrow and wait until it turns gray, the square is black, and the
LED is back to green. This should take only a few seconds at the most.
7.
Click the X at the far upper right of the form. This does not close the application but simply returns it to the icon tray.
Warning
As noted above, do not change any options or settings in the MPO Server without being
told to do so by technical support. Otherwise, the MPO client programs may no longer
work.
36
Astrometry in Canopus
Astrometry in Canopus
The following sections cover the steps used in Canopus to measure the position of an
asteroid (or other target) using the AutoMatch feature. As seen in the core lesson on
AutoMatch, Canopus also attempts to locate an asteroid based on information in the FITS
header and automatically measure its position and magnitude as well. It does not and
cannot measure multiple targets in the same field. You must measure additional items
separately. That will be covered in a subsequent tutorial. Targets other than asteroids,
e.g., variable stars, also are not measured automatically after an AutoMatch. Some tutorials in the Photometry section of this guide cover measuring such objects.
As part of the AutoMatch, Canopus automatically measures the reference stars on the
image and eliminates those that have excessive residuals. You can override the set of
reference stars that Canopus selects and the position that it uses for the asteroid. With
practice and a field that is not overly crowded, you can go through the entire process of
measuring an image in less than 10 seconds.
The tutorials in this section assume that you have gone through the core lessons and understand the basics of setting the configuration, setting aperture sizes and colors, and how
to AutoMatch an image (or, if necessary, do a manual match).
The objectives of the tutorials are to:
1.
Teach the basics of doing astrometry, which is the measuring of an objects position and, as a secondary goal, its brightness.
2.
Understanding the results and, if necessary, working with the data set to get the
best possible answer.
3.
Reporting your results to the Minor Planet Center or other institution or person.
4.
Practicing the core functions of opening images, setting the configuration, and doing an AutoMatch of an image.
37
Astrometry in Canopus
38
1.
Measuring an Asteroid - I
The purpose of this tutorial is to measure the position of an asteroid. It is essentially the
same as the AutoMatch tutorial (pg. 17). The difference is that youll use different images.
Be sure to keep Canopus open after running this lesson so that you can move immediately to the subsequent lessons on reviewing and manipulation the results and then to
generating a report for the Minor Planet Center.
1.
Open the configuration form and duplicate the settings on the next page. If you
did the AutoMatch tutorial in the Core Operations section, you would just select
the profile that you created for that tutorial from the drop down list but confirm
that the settings match all pages.
2.
If necessary, open the apertures form and set the apertures as shown in the screen
shot.
3.
4.
39
The bottom screen shot is of the Save Options subtab on the Photometry tab of
the configuration form.
40
5.
If not already, open the Zoom box at 400% and move the mouse cursor over the
asteroid (yellow apertures in the screen shot) so that the mouse cursor is centered
on the inner aperture (the measuring aperture). Do not click on the asteroid.
6.
If you look closely, youll see that there is a very faint star just below the asteroid.
The screen shot of the zoom box shows a highly-processed section of the image.
The asteroid is the brighter (darker) of the two objects within the aperture while
the star is the smaller (fainter) object.
7.
Switch over to the Reductions page (press <Ctrl+3>) and look at the bottom part
of the form, where the M-C values are reported.
8.
Since the asteroid and star were both within the measuring aperture, the derived
position is not that of the asteroid but that of the combined objects. You would not
to submit these results as-is, unless you also indicated that the measurement was
affected by a field star.
41
Open the apertures form and change the target aperture settings to use a width of
9 pixels and height of 5 pixels. Click <OK> to save the changes.
This causes the apertures on the image to disappear. If you dont remember where
the asteroid was, do another AutoMatch.
11.
Using the zoom box to help locate the mouse cursor, click on the asteroid.
The screen shot above shows the placement of the elliptical aperture before clicking on
the image to get an updated position.
12.
After clicking on the asteroid in the image with the new apertures, go back to the
Reductions page (<Ctrl+3>) and look at the revised data for the asteroid.
13.
Note that the RA and Declination have changed significantly and that the residuals are much smaller.
14.
Before going on, restore the target aperture height and width to 11 pixels each.
42
If you want to generate a report for the MPC or someone else, you must save the results.
These include not just the position of the object but all the information on the Reductions
page. This allows you to load the file in the future and generate a new report or check the
results.
15.
16.
Assuming that your results from the AutoMatch are still there, click <Save> on
that page. This displays a Windows file dialog.
17.
By default, the program assigns a file name based on the number or name of the
object plus an auto-incremental number. For example, A1577_0001.AST where
the leading A indicates an asteroid. Its recommended that you change the name
to include the date of the observations, e.g., A1577_19990107_0001.AST would
show that the observations were made on 1999 January 7.
The program uses a default path of \MPO\CANOPUS\ASTFILES. You can
change this but its recommend that you keep all astrometry data this directory,
which was created when you installed Canopus.
18.
19.
Proceed immediately to the next tutorial, where youll measure a second image of
the same asteroid and well take a more detailed look at the results.
43
44
2.
Measuring an Asteroid II
The previous tutorial gave us a statistic of one, meaning that we have one position for
an asteroid with specific measured versus computed residuals. What happens if we measure a second image? Do we get the same approximate errors or something completely
different? What if we measure two more images taken on a different night and compare
the residuals to those from the first night? Most important of all, do we know that the
residuals are even remotely correct? The answer to that last question might surprise you.
This tutorial will also take a closer look at analyzing the results on the Reductions page,
beyond just those for the asteroid, and what they tell you about the overall quality of the
solution.
1.
Confirm that you are using the same configuration and initial aperture settings
that you used at the start of the previous tutorial.
2.
3.
4.
If you look carefully, youll see that the yellow set of apertures, those for the asteroid, are slightly off-center from the asteroid. Thats because a nearby star encroached into the aperture and skewed its placement.
5.
Note the large M-C value in Declination, no doubt due to including part of that
star above the asteroid.
6.
Were going to try cheating as we did in the previous tutorial. Reset the height
and width of the target aperture to 9 pixels. Do not change any other settings.
45
8.
The star is still just on the edge of the measuring aperture for the target but notice
how the M-C residuals have gone down considerably. Also note the significant
difference in the measured RA and Declination.
9.
You can see that there appears to be plenty of sky around the asteroid with the aperture set to 9x9. So, reset the apertures to 7x7 and do the AutoMatch again.
10.
The star is no longer affecting the asteroid and so this is as good as were going to
get. The residuals are much smaller and the position, especially in Declination has
changed yet again.
This is as far as you want to go. As noted earlier, you cant really see the entire
stellar profile and so to go smaller might mean cutting off a significant portion of
the wings of the asteroids image.
11.
Save the results under a unique name. Do not overwrite the results from the previous tutorial.
Do not close Canopus or do another AutoMatch. You want the current results available for the next tutorial.
Theres a saying that Its not cheating if you dont get caught. Cheating is cheating and to quote the famous phrase, Theres no such thing as a free lunch. When
you use the off-center technique, be aware that the centroiding algorithm has to
work overtime and that you may not really be cutting off the interfering star completely. Both of these can affect the final results to one degree or another. However,
sometimes, there is no alternative and so you must accept the tainted results or
have none at all.
46
RA M-C
-0.764
-0.552
Dec M-C
+0.205
+0.432
What is important here is not so much the actual values but the trend of the values. In this
case, both of the residuals in RA are negative and approximately 0.5 arcseconds. The
Declination values are both positive and are centered on about 0.3 arcseconds. If we had
a third image for the night, we would hope to get residuals that are similar, i.e., about -0.5
arcseconds for RA and +0.3 arcseconds for Declination.
The M-C values are affected by many things. For example, you may not have the most
recent orbital elements from the Minor Planet Center and so the calculated position is not
right. Your longitude and latitude settings may be off, which affects the correction from
geocentric positions to topocentric (from the surface of the Earth at your location). This
correction can be significant for asteroids very close to Earth.
Its also possible that the times in your FITS header are wrong and so youll get a larger
error along the line of the asteroids path. If this is the case, then when you submit your
report to the Minor Planet Center, no doubt youll get a return message asking you to
check your measurements, including the times.
Again, the important consideration is the trend of the residuals. If a third set showed the
RA residual to be +0.5, then that should cause you to recheck the results. If they are the
same and you can find no reason for the difference, report the data any way. It is often
worse not to report what appears to be strange data than to hold onto it. There are many
cases in history where the unusual proved to be real and normal to be wrong.
Proof of Concept
To explore these concepts, repeat the two tutorials using images of the same asteroid
from a second night. These are found in
\MPO\Examples\ASTMTRY\JAN0899
Try using apertures of 11x11 to start. If necessary, use the techniques described above to
reduce the interference from any field stars.
Save your results and compare them to those from this and the previous lessons. What do
you find? If the trends are not similar, check for faint field stars.
47
48
3.
This lesson covers aspects of the Reductions page not covered in detail in the previous
two tutorials, most important, how to improve the astrometric solution.
To restate what was mentioned in the introduction of the Users Guide, your results
may vary some from what you see in this lesson. Do not get hung up on the fact that
the numbers are not exactly the same. The point here is how to interpret and manipulate the results to get the best possible solution. Naturally, if your results are
significantly different, you should double check your steps to make sure that you
have followed along as expected.
1.
Assuming that you still have the results from the previous tutorial, or that rerun it,
go to the Reductions page (<Ctrl+3>)
2.
Note the number in the upper-left corner of the table. This is the number of stars
used in the astrometric solution. Canopus limits the number to 75. This is more
than enough to get a good solution. Keep an eye on this number throughout the
tutorial.
3.
If you scroll down the table, you may see one or more rows where the RARes
and/or DCRes cells have green, blue, or red backgrounds. These indicate that the
RA and Declination for the star that is found from the astrometric solution differs
a significant amount from the catalog RA and Declination.
Green difference of 0.25 < x < 0.5 arcseconds.
Blue difference of 0.5 < x < 0.75 arcseconds.
Red difference of x > 0.75 arcseconds.
4.
If you see a row with colored cells, click on the cell in the U (used) column so
that the X goes away. As you do, keep an eye on the Measured position, the
S.D. (standard deviation of backfit residuals), and M-C values. Usually, youll see
very little change in these values. Those of primary concern are the S.D. values.
The lower these values, the better the astrometric solution. As mentioned in the
previous tutorial, the trend of the M-C values and not their actual values should be
taken as a whole over several images.
49
Note that the number of stars in the upper-left corner has changed as have the
RARes and DCRes values. The latter two are now 99.99999, indicating that the
star is not part of the calculations. If you move off the row, youll see that these
cells have a gray background, providing another visual cue about the stars status.
6.
Locate other rows that have colored cells and remove those stars from the calculations and compare the S.D. values to what they were when you started.
These are the results found running this tutorial. Note that the M-C values are the
same as those at the end of the previous tutorial. However, the S.D. values have
gone down by about 0.04 arcseconds in RA and Declination. Thats not a significant difference, but it shows how removing a few bad objects affects the overall
solution. The number of stars in the final solution, that number in the upper-left
corner of the table, is 56, meaning 9 stars were removed from the initial solution.
Other Things
The Reference Manual covers the rest of this page in more detail. For the moment, note
the three drop down lists under Observation Codes. Usually, you should not change
these, most important of all, Note 2 which is set to CCD and so indicates the observations were made from CCD images.
If you cannot completely eliminate a field star from the observation, you can use the
Note 1 drop down list and select I involved with star.
The Mag Code defaults to Do not use, meaning the derived magnitude that is reported in the Mag field on the page is not included. However, if you use the MPOSC3
or UCAC catalogs as the source for magnitudes, it may be appropriate to set this field to
the photometric band defined by Configuration | Photometry | Default Filter.
If you didnt already, be sure to save the final results (see the previous tutorials in
this section).
50
4.
The ultimate goal of doing the previous lessons is to generate a report for the Minor
Planet Center and then submit it.
Do NOT submit the results from these tutorials to the Minor Planet Center.
1.
2.
The settings in the screen shot above (except the TYP field) are the defaults immediately after installation. Before you actually submit a report to the Minor Planet Center, be sure to enter correct values in all the fields on this tab and save them.
3.
Assuming that you saved the results from the previous tutorials, go to the Reductions page (<Ctrl+3>).
4.
Click <MPC>. This displays a Windows file dialog where you give the name for
the report (text) file. All such files are forced to use the extension TXT.
5.
6.
51
This form is filtered for AST files, which are created when you save the results on
the Reductions page. This is a typical Windows multi-select form where you can
select one or more files, contiguous or non-contiguous within the list.
Note the dropdown list box Mag Code. This allows you to force using magnitudes
in the MPC report with a specific band, not using magnitudes at all, or having each
file determine how its magnitude data are written. Generally, youll select U use
original on the presumption that you saved the original file with the magnitude and
band setting that you wanted to be used for that data.
7.
Select the files that will be used to generate the MPC report. In this case, select
the files you created as part of the tutorials in this section and then click <OK>.
This presents the MPC Report form.
If you have TYP on the configuration form (see Step 1) set to Unidentified, a
small form appears after you click <OK> in Step 7 asking for additional information. See the Reference Manual
8.
The form contains a memo (text) control that you can change if necessary. The
control comes up in overwrite mode so that you do not change the columns in
which each data field in the report lines starts.
Be very careful when editing this report. You must keep data in their correct columns.
9.
52
The file automatically saved automatically as the form appeared. If you make
changes, then click <Save>.
Click <Email> to invoke your default email program and start a new message.
The To address is automatically set to the To email setting on the Configuration | MPC tab.
11.
Attach the just created file and send the message to the MPC.
53
54
Photometry in Canopus
Photometry in Canopus
The primary goal of photometry in Canopus is to measure asteroid (or variable star) magnitudes in dozens if not hundreds of images and then to determine the period and amplitude of the lightcurve. The process can be broken in to two steps.
Measure the series of images
Find the lightcurve characteristics (period and amplitude)
At first, the number of steps and detail required might seem daunting but, again, accurate
photometry is neither simple nor simply done. Go through each step carefully, refer to the
Reference Manual if you have questions about why, and before long the hard part will
be waiting to get the next series of images and not the processing of measuring them.
Once again, its important that you follow the steps below in order. You cant measure
images before you tell Canopus how to determine magnitudes and you cant determine a
lightcurve until you measure images.
There are other goals within the photometry features, for example, finding the color index
of an object or finding the transforms for your system for all-sky photometry or higherprecision/accuracy work. These will be covered as well.
A Matter of Magnitudes
Canopus always records the instrumental magnitude of a comparison or target. This is
computed directly from the sum of the ADU counts for the object and then normalized to
a 1-second exposure. This means that even images with different exposure times can be
directly compared among one another, assuming the same system was used.
The instrumental magnitude is independent of any catalog used. It is strictly a function of the system that took the image.
Canopus can also record a true magnitude (TrueMag from here on). This is based on a
solution using the measured instrumental magnitude versus the catalog magnitude for up
to 75 images in the star. For each star, the difference between instrumental magnitude and
catalog magnitude is computed (m-M). The mean and standard deviation of these results
comprise what is called the Magnitude/Intensity Relationship (M/IR from here on). If you
think about it for a moment, in a perfect world, the instrumental magnitude should vary in
direct proportion to catalog magnitudes, e.g., a star that is one magnitude brighter in a
catalog should have an instrumental magnitude that is also one magnitude brighter. Of
course, this is not the case.
Errors in the catalog values, measurement of the instrumental magnitudes, and differences in the color of the stars used in the data set contribute to some degree of uncertainty. These are all wrapped up into the standard deviation of the mean of the m-M differences, or at least to a first order approximation. These are only the internal factors.
External conditions such as variable observing conditions mean that the solution for one
image may not be the same as for another image, e.g., if a slight haze moves in, the instrumental magnitudes are all fainter, and so the offset (m-M) is different.
The MPOSC3 catalog provided with MPO software includes a large subset of the Carlsberg Meridian Catalog (CMC-14) as well as from the Sloan Digital Sky Survey (SDSS).
More specifically, the catalogs strongly favor stars that about the same color as the Sun.
These magnitudes are mostly SDSS r with some g and i when available. Almost all
stars in MPOSC3 also have BVRI magnitudes derived from 2MASS J-K magnitudes (see
55
Photometry in Canopus
Warner, 2007; Minor Planet Bulletin 34, 113-119). These have an internal consistency of
~0.05 mag for V and 0.03 mag for R. Many tests have been run using the magnitudes
from the MPOSC3 catalog, mostly to link data sets from many nights. Generally, individual sessions fall into line to within 0.02 mag when using a calibration method involving the Comparison Star Selector and instrumental magnitudes.
Differential Photometry
This is where the difference between a target and comparison (or average of several comparisons) is found for each image. It is this differential value that is used in analysis.
Since, in most cases, any variations caused by external factors are the same across the
entire image, this approach nullifies those factors. What remain are the internal factors,
the differences in color between the object and comparison stars usually dominating
things.
You must take care not to use stars that are beyond the linear response of the camera, let alone that are saturated even if its only one pixel out of dozens.
Instrumental Photometry
For MPO software, this will mean that instrumental magnitudes are used for differential
photometry and that no offset is applied to put the differential values on a standard system. These will be called Instrumental in this Guide.
Catalog-based Photometry
This approach uses the M/IR to compute the magnitude of the target directly and ignores
the comparisons altogether. If the M/IR is computed image by image, then this reduces
but does not eliminate external factors. Under good conditions, this approach produces
good catalog-based magnitudes. However, one cannot be sure if good conditions prevailed throughout a session. These will be called TrueMags in this Guide.
56
Photometry in Canopus
That standard deviation incorporates the uncertainty in the measurement of the target and
each comparison, the uncertainty in the catalog value, the uncertainty in the correction for
color differences (if used), and so on. It is not a perfect predictor of the error budget, but
it is better than much simpler methods.
The DerivedMags method is the recommended approach for almost all measuring
for time-series photometry. It gives catalog-based values while retaining the (considerable) virtues of differential photometry using instrumental magnitudes.
All three of these methods will include corrections for changing geometry and/or
phase (if appropriate) when doing period analysis. These corrections are mandatory
in order to do proper analysis. The Transformed Absolute method does not include these corrections.
mo
Tr
V-R
Zo
When using the Comp Star Selector (covered in the tutorials), Canopus can transfer
the B-V or V-R color indices of the comparisons from the MPOSC3, making them
available for applying transforms on-the-fly. As for the target, that is taken from the
current PhotoRed transforms data. This is covered in more detail in the tutorials on
finding a lightcurve.
57
Photometry in Canopus
The Differentials
Lets go a step further and look at the internal workings of Canopus when computing the
instrumental magnitude difference between the target and a comparison star. A subscript
of t indicates the target and one of c indicates the comparison. The o subscripts in
the formula above are presumed. CI represents the chosen catalog color index, e.g., V-R,
for the target or comparison.
m = (mt - mc) + Tr(CIt CIc) + (Zt Zc)
Obviously, the Z values are identical on a given night and so that differential is 0 and
drops out. This leaves only the difference in the instrumental magnitudes plus the correction for the difference in the color indices of the target and comparison which, as noted
above, can be transferred into the current session when using the Comp Star Selector. The
net result is that Canopus computes the catalog color index differential and applies it to
the difference in instrumental magnitudes to derive a transformed differential magnitude.
Dont let this confuse you too much. The tutorials will demonstrate these principles.
Whats important to understand is that by taking the time to find the transforms for your
system (if you have the necessary filters), you can increase the accuracy of your results
and, within some limits, have a little more flexibility when choosing comparison stars by
not having to keep them as close to the same color (same color index) as the target.
Deprecated Methods
The Transformed and Transformed Absolute methods have been retained in the current version of Canopus but should be considered deprecated, meaning that they are no
longer recommended and used very infrequently. The main reason for retaining them is to
work with data that has been converted to standard magnitudes in PhotoRed and then
imported back into Canopus. In those cases, you do not want to make further use of the
comparison star set since that was done in the PhotoRed processing.
As mentioned above, you should use the Transformed method if working with the
TrueMags found during AutoMatch (see below).
58
Photometry in Canopus
The advantage to this approach is if the tracking isnt very good or you deliberately dithered around the target. In this case, Canopus may not be able to find the set of comparisons and target. The AutoMatch routine knows the RA and Declination of whats
called the anchor star (see the first tutorial in this chapter). After an AutoMatch, it can
find that star and, using that as the reference point, the remaining comparisons and target.
Remember: you do not have to use the resulting catalog magnitudes for the target
for analysis. Since the instrumental magnitudes are measured regardless, you can
use the Instrumental or DerivedMags methods. The main (really, only) reason for
the AutoMatch processing is to let you do other things on a different computer or
away from the processing computer while Canopus is measuring images. You should
not try to do things on the processing computer during auto-processing since the
program is constantly trying to get focus.
Make it a Habit
You should always use the Comp Star Selector. The most important reason is that it records the catalog magnitudes of the comparison stars and can transfer those into Canopus
so that they become available for using DerivedMags approach. If doing Instrumental
analysis, the average of the comparison star catalog magnitudes can be used to set the
nightly zero point for each set of data. Linking sessions via this approach has been quite
successful in work on very long period asteroids.
Finally, while you may hope to use TrueMags, those many external factors mentioned
earlier may come into play such that non-differential photometry produces very poor
results. With the DerivedMags approach, you still get catalog-based magnitudes that are
much more reliable and stable.
The benefits of the CSS are many while the time to set it up it is minimal. Use it!
59
Photometry in Canopus
Why have separate sessions? Its to account working the general case of a moving target
(I think of variable stars as a special case where the object is not moving which can be
approximately true for an asteroid near one of its stationary points). In the general case,
youre not assured of using the same comparison stars each night, sometimes over the
course of single night. Since the target magnitude on each image is based on the comparison star set one way or another, any shift in the comp star set average value affects the
target magnitudes throughout out the session. Even if you are using the same comparison
star set, many factors including just bad luck (aka Murphy) can mean a slight shift
in the data from session to session. This is called zero point shift.
By having separate sessions, you can set the zero point for each session separately and so
account for minor variations in the zero point over the duration of the observations.
When using the DerivedMags method with SDSS magnitudes, especially if applying
color corrections, it is rare to have to apply any correction to the zero points. If so,
its on the order of 0.01-0.02 mag. When using the 2MASS to BVRI values, the shifts
are usually on the order of < 0.05 mag, though 0.1 mag is not impossible. This demonstrates the importance of using native catalog magnitudes that have a good internal consistency.
The second important reason behind sessions is the most critical when working asteroids.
All other things being equal, an asteroid changes brightness as its distances from the Sun
and Earth change. Canopus can correct for these changes so that all magnitudes are based
on the notion that the average magnitude of the asteroid is always the same, no matter
how much it fades or brightens due to changing distance.
Asteroids also change brightness with changing phase angle, which is the apparent separation of the Earth and Sun as seen from the asteroid. For example, when the asteroid is at
opposition and the Earth and Sun are in a line with the Earth on the asteroid side of the
Sun, the phase angle is 0. There is a very distinct change in rate of change of the magnitude as the asteroid phase angle is near opposition (the opposition effect). This is most
pronounced with phase angles less than about 7. The relationship is complex and calculated using the H (absolute magnitude) and G (phase slope parameter) values that you
find the Minor Planet Centers MPCORB file (or Lowell Observatorys ASTORB,
among others). Canopus computes the correction due to changing phase and adjusts the
objects magnitude accordingly. You can even use a custom value for G should it be determined by measurements or to test the affect of different values on period analysis.
As you might have guessed by now, the phase correction varies constantly and so a new
value is needed for each session. The issue is even more complex than whats been explained before. See the Reference Manual for additional information.
A third reason, related more to processing observations than making them, is that separate sessions are needed for observations in different filters for proper processing in
PhotoRed. Such processing will be covered in later tutorials.
Just remember: one object on one night using one filter and one set of comparison stars
equals 1 session and youll be on your way to understanding sessions and why they are
necessary.
The concept of sessions still holds even if using DerivedMags or TrueMags. The
corrections for changing distance and phase must still be included in order to do period analysis. What you dont have to worry about is the different values for the
comparison star set average. Even so, you may still have to tweak the zero point a
little because of the inherent errors in the catalog magnitudes for the comparisons,
among other things.
60
1.
A First Lightcurve
This tutorial will take you through finding the lightcurve of an asteroid on a single night.
At times, youll have several forms open at once. Position them the best you can so that
you can see as much of each one as possible.
Remember that DerivedMags is the recommended method for analyzing data. This tutorial will generate the necessary data to use that method or instrumental magnitudes (both
still taking advantage of differential photometry). AutoMatch will not be used here so
that you get data as quickly as possible. When working the second night of data for this
asteroid, AutoMatch will be used, first to demonstrate its use and, second, because the
field of view jumped around considerably that night and youd be resetting the location of
comparison star apertures a number of times. In the end, it would still take less time to
work without AutoMatch but the frustration factor and being able to have the measurements take place unattended after they are started are both worth considering.
1.
There are several options when measuring images that affect what data are available for analysis. With such flexibility comes complexity. Take your time and try to
understand the why and much as the what that are covered in this and following tutorials.
Open the configuration form and match the settings on the screen shots (the profile name can be different if you want).
61
Note the bottom two shots are of the subtabs on the Photometry tab of the configuration form. These subtabs include options that control the plotting of the data when analyzing a lightcurve.
62
Set the apertures to use 11x11 for target, comparisons, and astrometry.
3.
Open the first image for the first night for 771 Libera from the Examples folder,
\MPO\Examples\LTCURVES\A771\SEP18\A77A001.FIT
In general, the image you open just before creating a session and running the lightcurve wizard must be the first image that you will use in the lightcurve wizard. The
image does not have to be the first image that will eventually measure. This is explained in more detail later in this tutorial.
There is another important consideration for the first image. It should be of higher
quality (good guiding, no clouds, etc). This is particularly important if you will be
using the star subtraction feature (StarBGone) since you need to be able to see and
select all potential subtraction stars.
4.
<Ctrl+Alt+S>
<Ctrl+Alt+C>
moves the splitter so that only the chart is seen on the Measurements page.
You should always AutoMatch the first image to be measured just before creating
a session. The plate constants are used to find the asteroid on the images opened in
the Lightcurve Wizard (see below) and allow using the Comp Star Selector. It also
allows data about the M/IR to be transferred to the sessions Notes field.
5.
6.
63
8.
Enter 771 Libera (remember, without the quotes) in the Object field. The drop
down list contains up to the last ten objects worked, but you can also type in the
name of a new object.
64
9.
It is very important that you use the exact same name for a given object for all sessions. Canopus groups sessions by name. The grouping is not case-sensitive (so
Fred and FRED will be in the same group) but it is sensitive to any other differences.
The Mid-Date should already be set to 1999/09/17. This was taken from the FITS
header. Unfortunately, the header was not written in UT and so the date is off by
one day.
Change the date to 1999/09/18.
10.
Generally, all date fields in Canopus use ANSI format, yyyymmdd (year, month,
date), always with four-digit years. Time fields always use 24-hour format, e.g.,
13:05 and not 1:05 PM.
The Mid-Time field indicates the approximate middle UT of the observations to
be made. Dont get caught up on having this value exactly the middle. Somewhere close is good enough.
This field defaults to 07:00 UT, which is reasonable for data taken in the continental U.S. and good for this tutorial.
11.
Enter information in the Telescope and Camera fields (or use the drop down
list if there are previous entries).
12.
The F.L. (focal length) can be given in millimeters or inches. This is an informational field and is not used elsewhere.
13.
The Temp is the CCD operating temperate (in C). This is an informational field
and not used elsewhere.
14.
The Exp is the exposure in seconds. If available, Canopus takes this from the
FITS header of the image that is open when the sessions form is opened.
15.
If the Filter setting is not C (clear filter), select C filter in the drop down list.
If you do change to C, a form appears to set the Mag Band value.
If you are going to apply transforms on-the-fly, be sure to select the filter that was
used to take the images and that you have good transforms for that filter. Canopus
reads this field to determine which transforms from the PhotoRed transforms data to
use.
Note that the Filter combo box is automatically set to the Default Filter from
the Configuration | Photometry tab.
16.
17.
Click <Calc M/D/P>. This displays a form where you pick the asteroid being
measured.
65
If youre not working an asteroid, you should still enter the RA and Dec values. See
below.
18.
Locate the asteroid in question by scrolling down the table or using the search
controls on the right side and bottom of the form.
19.
Once the asteroid is highlighted, click <OK>. This automatically enters data in
the fields within the same area as <Calc M/D/P>. The values in these fields are
used to correct for changing distances and phase angle, light-time travel, and to
compute the air mass of the field on-the-fly since the AIRMASS FITS keyword is
not assumed to be available. See the Reference Manual for additional information.
20.
If you used the Asteroid Lookup form or are working a variable star, check the
entries in the RA and Dec fields. This is the approximate RA and Declination of
the center of the field being worked. The data are used to compute air mass and,
for a variable star, the heliocentric JD. See the Reference Manual.
21.
22.
Click <OK> to close the sessions form and make the newly created session the
default session. The default session is the one to which any new data are added or
in which existing observations are edited.
23.
Open the Canopus Lightcurve Wizard (<Ctrl+Shift+W>) or Photometry | Lightcurve wizard from the main menu. This displays a dialog which is a Its 2 AM
and the operator may doing something that he really doesnt want to do warning.
24.
Click <Yes>. This displays the first page of the Lightcurve Photometry Wizard
(called the wizard or lightcurve wizard for the rest of the Canopus photometry tutorials).
66
Even thought it is already open, click the button next to the top entry field and
load \MPO\Examples\LTCURVES\A771\SEP18\A771A001.FIT.
26.
Check the Star Subtraction box. Make a habit of doing this if working a moving
target. In the end, the target may not go near any field stars and so star subtraction
is not needed. Its easy to change your mind after you know for certain that you
wont need to use star subtraction (This feature is called StarBGone in the Reference Manual; be sure to read that section).
27.
Click <Next> to see a page of informational text. Click <Next> again to get to the
page where you define the comparison stars on the first image.
28.
Click <Show Path>. This displays a form similar to the one used in the Sessions
form to pick an asteroid.
29.
67
30.
You should see something similar to above, i.e., a small circle around the asteroid
and a line off in one direction. At the end of the line is the length of the line, in
hours of asteroid motion.
Make a habit of using the Show Path feature. Its most important use is picking the
comparison star set on the first image so that you wont pick comparisons that are
too close to the asteroids path. Also, you cannot show the path on the second image
unless youve activated it on the first image. You must have done an AutoMatch on
the image being used for Image 1 for the Show Path feature to work.
31.
Click on the asteroid, the star at the end of the line, surrounded by a small circle, and away from the +12hrs label.
32.
33.
Click <Selector> the wizard. This displays the Comp Star Selector (CSS from
here on).
Reference the screen shot below for the next few steps.
34.
68
If you are youre going to use TrueMags and the Transformed method, then its
very important that you select the same filter as the Default Filter on the
Configuration | Photometry tab. The configuration setting determines which catalog magnitudes are used to establish the M/IR (magnitude/intensity relationship).
You want to pick comp stars based on the same photometric band.
If youre using Instrumental or DerivedMags but havent found or dont want to apply
transforms, the Default Filter and the filter you select here should be the same and
come closest to the band of your system. For example, if using the Clear filter, many
systems using the Kodak chips (KAF series) come close to matching the R band. So the
R filter is probably a good choice.
If at all possible, use the r magnitudes for unfiltered systems. These are much more
consistent and you will get tighter results (lower standard deviation). Unfortunately, the SDSS and CMC-14 catalogs do not cover the entire sky and so you may
be forced to use the 2MASS to BVRI magnitudes. Which ever you chose, make sure
to use the same setting for every session for the object. Otherwise you will have unpredictable zero point shifts that can and will affect period analysis.
The Clear filter has a significant color dependency and so it is critical that you use
comparisons that match the color of the target as closely as possible.
Put another way, the transform for a given filter is 0 for a perfect match to the standard
system, meaning that there is no difference in the color correction over a wide range of
colors. This is not true for the Clear filter. A typical example at my location is Tc ~ 0.3 *
(V-R). So, if the differential color index between target and comparison is only 0.1 mag,
the error by not applying a transform will be 0.03 mag (0.3 * 0.1). This can be and often
is a significant error when trying to get consistent data and linking data from one night to
another night.
Continuing On
35.
Click <Plot Comps>. This finds stars in the chart that was generated during the
AutoMatch done just before staring the wizard. If you dont do an AutoMatch,
youll get an error message about there not being enough stars (youll also get that
warning if you pick a magnitude band for which values are available for very few,
if any, stars in the field.
36.
You should see a plot similar to the one above. A perfect world solution would be
a slope of 1.000 and S.D. of 0.000 mag. That never happens.
69
If you click the Solar radio button on the panel towards the bottom, the plot is redrawn
to include only those stars with 0.5 B-V 0.9, or approximately the same color as most
asteroids, which are simply reflecting sunlight but affected some by surface minerals.
When you click the Solar radio button, the Show on Image check box is enabled. If
you then check the box, or its checked when you go to solar-colored only stars, circles
are placed around the solar-colored stars in the image as shown in the screen shot below.
If using a dark background, the circles are bright green. On an inverted image, as in the
screen shot, the circles are magenta.
70
If you uncheck the box, or select All, the circles are removed.
The speed button next to the check box forces the circles to be hidden. Internal flags can
get out of step and hide the circles even though the CompStarSelector thinks they are
visible. If you click the button, youre asked if the circles are being displayed. Give the
appropriate answer and the CSS then hides the circles and unchecks the Show on Image box. The button is enabled only when plotting only solar-colored stars.
For this set of images, finding a good set of comp stars was not hard. If the field had been
much more crowded, then using the solar-colored star only feature makes picking comparisons much easier when working targets like asteroids and stars of about the same
color as the Sun.
71
37.
72
38.
Note that as you click one of the CSS comp buttons, the X/Y location of the star is
recorded in the lightcurve wizard fields as well. Do NOT click <Star X> buttons on
the wizard. While the X/Y position will be recorded, the magnitude information that
the CSS would retrieve is not recorded.
After you have measured the fourth comparison star, click the Comps tab at the
top of the CSS notebook control. This displays the magnitude data for the chosen
comparisons.
Canopus Transfer
Its just a good idea to keep a permanent but handy record of all the contents of the CSS
comps memo. For example, if you are doing Instrumental photometry and it turns out that
one of the comparisons is variable, you must first eliminate that comparison from the
calculations and then find a new average magnitude for the comparisons to use in the
sessions DeltaComp field. If you have the CSS data available, this is very easy to do.
The CSS can transfer the average value to the DeltaComp field, the individual magnitudes of the comparisons to the respective CM and CMCI (catalog magnitude, color index) fields, and the contents of the CSS memo to the Notes field of the current session
without having to go back to the Canopus main form and opening the sessions form.
39.
Do not check the Avg to DeltaComp box. This is because were using the DerivedMags method. If we were using the Instrumental method, then you would
check this box.
40.
41.
42.
Check the V-R button. This transfers the V-R magnitudes to the sessions form.
If you did (or will do) transforms based on B-V, then you would check that button.
If you are applying transforms, you must select the color index that you used for
finding the transform. For example, if when finding the transform for the Clear filter
you used V-R catalog magnitudes, then you must select V-R here. Making a selection does not automatically mean that transforms will be applied; it means only
which values are used should you chose to apply transforms.
73
Click <Transfer>. This displays a confirmation message that the data were transferred or, if not, the nature of the problem.
In case youre curious, the screen shots above show the session form just after the
transfer. The C1-C5 names have not yet changed (they will when you close the
wizard) but the individual catalog magnitudes for the selected comps have been
entered. Since only four comparisons are being used, the C5 catalog magnitude is
0.0. The Use box is still checked, but will be cleared when you finish the wizard.
The small form to the right is the second tab of the catalog magnitudes and shows
the color index values that were transferred. The 0.45 for the ObjCI came from
the PhotoRed transforms data and is the approximate default V-R for solarcolored objects.
Check the Mag Band setting. Remember that when you created the session, that
field was set to the Default Filter on the Configuration | Photometry tab. If you
used a different filter in the Comp Star Selector, the Mag Band setting is changed
to match the CSS filter so that you know the photometric band of the comp star
magnitudes.
Do not edit the Sessions form manually at this time. Finish the wizard before making
any changes.
74
Click <Save to Text> on this page of the CSS. This saves the contents of the
memo field to a text file with the default directory being the same as where the
image was found.
This is somewhat redundant to transferring the text to the session Notes field, but
it does provide a second copy in a separate file that you could send to another observer.
Make a habit of saving the CSS data and saving the data in the separate text file.
You will need it more often than you think and the only way to get it back if you
dont save it one place or another is to run the wizard again.
Moving On
45.
After you close the CSS, the wizard should look something like this.
46.
47.
Use the speed button next to the entry field to bring up a Windows file dialog and
load \MPO\Examples\LTCURVES\A771\SEP18\A771A045.FIT.
75
When the second image appears, a set of numbers and T appear superimposed.
Numbers 2 and higher and T are located in reference to the number 1 star (the
anchor) based on the offsets in X and Y measured in the first image.
Many times, the labels will not be exactly next to the stars and target. This is particularly true if you repositioned the telescope during the observing run to the asteroid near center. This is not a problem and is easily fixed.
If you invert the images in the wizard, it may be a good idea to select Black numbers so that the labels are easier to see. If you have trouble seeing the numbers in
this exercise, you can use the Back button to get to the first page, check that box,
and then return to this page (Image #2) of the wizard.
49.
As you can see, the numbers are not next to their corresponding comparison stars.
You can reposition the entire set of apertures by moving the mouse cursor over
the anchor star (Comp 1) and then <Ctrl+Click>. Note how all the labels move to
their now positions.
This provides a visual check that you have correctly identified the anchor star on
the second image, i.e., all the labels except the T for a moving target should
76
If the labels are close to where they should be or you dont care to reposition the
labels, just <Click> on the anchor star so that its X/Y coordinates are recorded in
memory and available for the next step.
Sometimes a label may be on top of the anchor star and you wont be able to click
on the star (the cursor disappears as you move over the star). Click the speed button
to the right of the Clear button so that it appears raised (it should appear depressed). This hides the labels so that you can click on the anchor star.
The screen shot above shows the target apertures around the anchor star and labels next to their respective comparisons.
50.
Click <Star 1> button on the wizard (make sure you are on the Image #2 page).
Note how not just the anchor star position is entered but those of the other comparisons as well. These will be updated a little (to account for image rotation) before the final results are derived.
51.
Where is the target? Click <Show Path>. As with Image 1, a line appears that
shows the path of the asteroid but this time there is a small circle at a point along
the path where Canopus predicted the asteroid would be. See the screen shot below.
52.
53.
Click <Target> on the wizard. The X/Y values are entered on the wizard.
77
Click <Next>. This displays information about star subtraction because you
checked that box on the first page. If you did not check that box, then the following steps can be skipped. In fact, you must skip them because the relevant wizard
pages will not be displayed.
If after reviewing the paths from Show Path in the previous steps you determine
that the moving target does not come close to any field stars, you can either go
through these steps and specify a random star as the subtraction star (dont make it
a comp star!) or go back to the first page of the wizard and uncheck Star subtraction. Then click <Next> to get to the final page of the wizard to close it out.
StarBGone does not have infinite power. It cannot remove a 10th magnitude star as
the 16th magnitude asteroid goes over it. Statistics of random noise will show that at
some point the noise in a given stars measurement is equal to or greater than the
entire contribution of the target. Do not expect major miracles with StarBGone but
do expect to recover data in a number of cases you thought might not be possible.
55.
Click <Load Image 1> to load the same image that was used as Image 1 earlier in
this tutorial.
56.
A line appears going through the target that is derived from the positions you gave
for it on the first and second image. Any field star that touches this line or lies sufficiently close to it should be marked for subtraction.
It is much easier to see faint field stars when you invert the image (black stars, light
background).
78
Click on comparison star 2 on the image and then click <Set Ref Star 1> on the
wizard.
58.
Move the mouse cursor to X = 220, Y = 94. This is a star near a brighter star, the
latter of which is close to the asteroid path.
59.
Click on the star and then click <Set Ref Star 2>.
60.
If it isnt still open, display the zoom box. This helps see faint stars along the asteroid path and position the cursor.
61.
Locate stars near the path. For each one, click on the star and then click <Add> on
the wizard. To delete an entry, click on it in the list, and then click <Delete>.
79
You do not have to mark every star near the path for subtraction, only those that are
close enough that they might intrude into the measuring area of the target aperture.
The sky background algorithm in Canopus can handle fairly bright stars that fall
into the sky aperture. However, it cant hurt to mark somewhat brighter stars that
would intrude into the sky annulus in order to reduce the stress on the algorithm.
There is no point in trying to remove an excessively bright star. Simply skip those
images where the target and star are too close.
62.
After you have marked all the subtraction stars, click <Next> to go to the final
page of the wizard.
63.
Click <Save Comps>. This saves a copy of Image 1 with the number labels included as either a PNG or BMP file. This gives you a quick visual guide to the
comparisons that you selected. The default directory is the same as for Image 1.
64.
80
The Canopus Image List provides the interface to measuring the selected images. The
Target and Comp buttons determine which set of apertures are used when you click
on an image while processing is paused. In most cases, you will use Target since youll
be tweaking the position of the target apertures if they begin to wander away from the
target.
A quick note on the Beep On control. If checked, Canopus sounds a short beep after an
image is measured. This can be both comforting and annoying. You can check or uncheck the box at almost any time.
A Bit of Philosophy
Before describing the various measuring methods available, let it be said that fullyautomated software has never been the goal or philosophy of Canopus. It is a good idea
to review every image at least quickly before it is measured and, if not up to standards,
reject it. If one trusts the computer without question, this could result in a set of data with
numerous bad points that must be reconciled one way or another. By reviewing the
measuring process and maintaining at least some control, you can reduce the number of
outlier data points considerably and so proceed to lightcurve analysis much sooner and
with a greater confidence in the data set and subsequent results.
81
The AutoTO field changes the default timeout, in seconds, allowed for Canopus to extract the stars from an image and do an AutoMatch. This is needed for cases where an
image was very badly guided and so the AutoMatch generates a chart that is way out of
scale to the image and has thousands of stars. This works in conjunction with the
Configuration | Charting option for Max Scale Diff to avoid the program appearing
to hang during an AutoMatch.
Make sure that you dont set the timeout so short that your computer cant extract
the stars from a crowded field before the timeout occurs.
82
66.
Position the Image List and zoom box so that you can see the middle part (if not
all) of an image as it is loaded by the Image List.
67.
68.
Match the settings in the screen shot immediately above. Make sure that the
AutoMatch box is not checked.
If you have doubts that the images can be extracted and matched in 10 seconds,
increase the value in the AutoTO to something you think is reasonable.
69.
Click <Auto>. This displays a form that allows you to select the auto-processing
mode.
70.
Check the Simple (AutoMatch optional) radio button and then click <OK>.
71.
Note that the Auto button caption changes to Abort while an image is being
loaded and measured. If you want to top auto-processing, click the button to display
a confirmation message. Click <Yes> to abort processing.
The first image should load and be measured automatically. After the second image is loaded, you should see an error message. As you can see from the screen
shot immediately below, the apertures are no longer centered on the stars and target.
The green-colored (lime) aperture set represents the anchor star. The other apertures represent the target (default yellow) and comparisons 2-5 (default red).
83
72.
73.
Click on the anchor star in the image. The apertures should all move in unison
such that they are centered on the comparisons and target again.
If the anchor and other comparison star apertures are well-positioned but the target
aperture is not, <Ctrl+Click> on the target on the image. This repositions only the
target apertures.
If you accidentally click on the wrong star for the anchor star, double-click the images name in the Image List files list. This reloads the image so that you can click
on the correct star.
Do not click on the star again without reloading the image. This does not reposition
the other apertures.
74.
The highlight in the files list should be on the second image. If so and the apertures are all where they need to be, click <Auto> again. You will see the same
query form as before. Check Simple (AutoMatch optional).
75.
Since the highlight is not on the first image in the list, you see a new confirmation
message.
76.
Click <Yes> to start with the second (highlighted) image. If you click <No>, the
processing starts with the first image in the list.
Continue as above until all images are measured. You will have to reposition the apertures a few times but, for the most part, Canopus will measure merrily along.
Do not check the AutoMatch box during this tutorial, e.g., when paused to reposition
apertures.
77.
84
Once the last image is processed, click <Close> on the Image List.
Aborting Auto-processing
If you need to stop auto-processing, for example, the target aperture is drifting away from
the target and needs to be reset, click <Abort>. When you click <Abort>, a message appears confirming if you want to abort. Click <Yes>. If you have AutoMatch turned on,
the processing may not stop until after the match is performed.
79.
Enter these settings in the analysis fields that are found at the top left of the
Analysis page.
Reduced Magnitudes
Reduced magnitudes, for the purposes of Canopus in this case, are those that have been
corrected to unity distances at a fixed phase angle, i.e., the magnitude that the object
would have if it were at 1 AU from the Sun and 1 AU from the Earth at the given phase
85
Even when using Reduced magnitudes, the correction for phase angle is included
based on the value for the phase slope parameter (G) thats entered in the sessions
form.
80.
81.
Make sure the All button is checked. Otherwise, you see only those sessions
where the Filter was set in the sessions from to the filter selected here.
82.
Check the session in which youve been working (the current default session).
If this is the first time youve measured an object with the name you used in the
Object field of the sessions form, there will be only one item. If youve done other
measurements, more than one item appears. Select (check) only those sessions
you want to analyze.
This is a typical Windows multi-select list where you use click, Shift+Click,
and/or Ctrl+Click to select one or more items.
This shows the importance of giving the same name in the Object field on the sessions form for observations involving the same object. Canopus uses the Object field
to search for all sessions involving a given target. Only those that match the name in
the Object field of the default session are found.
83.
Once you select the session on which were working, click <OK>.
84.
86
Two other forms are also displayed. These are not shown or discussed here. Well
look at those in detail in a following tutorial. Close those forms for now.
This is a plot of the raw data, meaning that no attempt has been made to fit the
data to a specific period. Since this is a raw plot, the period is not displayed.
Phase angle, the angle between the Sun and Earth as seen from the asteroid, is designated
with the Greek letter (alpha). Since theres no assurance that the character will be
properly displayed, it is spelled out on the plot. This value tells you the reference phase
angle for all observations.
Remember that all data must be corrected to remove effects not only of changing distance
but changing phase angle. The corrections are calculated using the value for G (phase
slope parameter) entered in the sessions data and a complex formula that predicts magnitude based on that value. The phase angle for the first session is the zero point for all
other corrections all data are corrected to this phase angle. Correcting to 0 requires
accurate knowledge of G, which is not often available.
If the phase angle for the first session (usually earliest by date) has not be computed, then
? replaces the value in parentheses.
Since Reduced magnitudes were selected, the title includes Reduced. Had that
box not been checked before clicking <Find>, the title would have been just Magnitude plus the phase angle information.
The industry standard for Reduced magnitudes is the V band. If you do not use V
magnitudes, be sure to note which band was used when you publish your results
(you DO publish your results, dont you?).
87
Click <Notepad> (notepad icon) at the far upper right of the Canopus main form.
This displays the sessions form in edit mode on the current default session.
Before moving on, note that the C1-C4 names were changed from A, B, etc.
They are now a combination of the RA and Declination of their respective positions
without the field delimiters. This allows you find those comps on any image (since
their X/Y values may change) or in an on-line search. C5 has been set to blank and
the Use box is no longer checked.
86.
On the Session Data tab, click <Calc DC>. This may display a message (it
shouldnt in this case if youve followed the instructions carefully).
87.
Click <No>. This causes Canopus to recompute some internal values regarding
the comparison star set so that you can plot their data.
In most cases, click <No>. You may have spent some time finding the correct
value for the DeltaComp value. You dont want to override it without first being
asked.
The DeltaComp values are used with differential photometry (when using DerivedMags or Instrumental methods) to correct for the difference in the average magnitude of the comparison star set from one session to another.
88.
89.
If necessary, click on the Data | Raw and Comp | 1 radio buttons on the right
side of the form. You should see something similar to the screen shot below.
88
90.
Click on the other Comp radio buttons to see the raw data for the other comparisons.
Note that the data are tight, there is little scatter, but that they are not flat. Instead
they rise and then level off. This is normal. You are seeing the effect of extinction
that causes a star to get brighter as the field rises higher above the horizon. Near
the meridian, the magnitude changes very little. If the session had gone longer, the
curve would have started going down.
As with M-C values in astrometry (see those tutorials), you are looking for a similar trend among all the comparisons. If they all show about the same scatter and
trend, then all is good.
91.
You can confirm the quality of the comparison star data further by clicking the
Data | Avg radio button. The plot will change slightly (in this case).
92.
Now you see (as you would hope) the data on a nearly flat, horizontal line.
Here the program is plotting the differences between the instrumental magnitude
of the select comparison and the average instrumental magnitude of the remaining
comparisons. This is akin to the old C-K values that showed the stability of the
measurements.
At the top of the plot, you see the S.D: 0.007 (in this example). This is the standard deviation, in magnitudes, of the average of all the data points. You want
something on the order of 0.01 mag. Here, its even better. Sometimes, however,
89
Fixing Things
If plotting the Avg and the line is not flat but, for example shows a sinusoidal shape or
obvious eclipse, you have used a variable star for a comparison! If you select other stars,
the trend will reverse itself for those stars as the variable becomes part of the average of
the remaining comparisons.
Instrumental Method
If you are working the Instrumental method, you must change the DeltaComp value to
the average of the catalog magnitudes of the comp stars still being used. If you transferred the magnitudes to the session when using the Comp Star Selector, this is very easy
to do.
As you check or uncheck the Use box on the Session Data, the average of
the catalog magnitudes is displayed immediately below the C5 CM entry
field. Enter that value in the DeltaComp field.
Go to the Comparison Plots tab and repeat the exercise to be sure that no other
star is a variable.
If you did not transfer the values but did save the Comp Star Selector data to the NOTES
field:
Write down the appropriate catalog values for the comps (all of them so that
you do this once).
DerivedMags Method
All that is necessary is to uncheck the Use box of the offending comparison. You
should leave the DeltaComp value at 0.000, or whatever the value is after you did some
tweaking to align sessions.
90
The data points are not deleted from the data table, just excluded from being used.
You can restore data points (or exclude them) by going to the Observations tab
and either checking or unchecking the Use box for a given observation.
93.
Once you are satisfied with the comp star data, return to the Session Data tab,
click <Save> and then click <OK>.
94.
Click <Find> on the Lightcurve Analysis page of Canopus and repeat the steps
from above. The plot reflects the changes.
If you just click on a point, information about it appears on the lower status bar.
This includes the JD (light-time corrected or Heliocentric JD) and the session number from which the data point was taken. The latter is handy if you are plotting data
from more than one session.
95.
This feature can be used whether using Instrumental or DerivedMags but only if you
transferred from the Comp Star Selector or manually entered the catalog magnitudes for the comparisons in the CM field so on the sessions form.
Click on the Catalog Check tab to display the DerivedMags data in a different
way.
91
Raw
This screen shot shows the Raw view of the Catalog Check. Here, the Y-values are the
derived magnitude for the asteroid when using one of the comparisons. In a perfect
world, the plots would like exactly on top of one another since the derived magnitude for
the asteroid should be the same regardless of which comparison is used. If the plots are
not well-aligned, that indicates, mostly likely, that the catalog magnitude for the corresponding comparison is bad.
In this case, it appears that Comp1s catalog magnitude is a little too low since its curve is
below the others, but only slightly.
Note that the CI Corr button is set to None. If you were applying transforms,
you would want to check Apply to see how applying the transforms affects the
plots. This is just for the sake of comparisons here. You must still choose to apply
transforms when finding a lightcurve to apply the color corrections during period
analysis.
Average
The average plot shows the data for each comparison as the difference between the derived magnitude using the selected comparison star catalog magnitude minus the average
of the derived magnitudes for the object using the remaining comparisons. Here again, in
that perfect world the plots would sit exactly on top of one another.
92
With Comp 2 used in the calculations, the DeltaComp for this session had to be changed
considerably from 0.00 to make it mate with sessions from the same night. When Comp 2
was removed, the session lined up perfectly with the others when DeltaComp was set to
0.00.
On the left is the Raw plot with no color correction included. On the right, the color index
correction has been applied (V-R in this case). You can see that the five curves come
much closer to matching and so the benefit of using transforms, assuming that good color
index values and transforms are available.
The color index values in the sessions form for the comparison and object must be
based on the same one used to find the transforms.
93
In Conclusion
This concludes this tutorial. Its been long but, I hope, productive. The data stored here is
used in a following tutorial, so be sure not to delete this session.
Check the Reference Manual for more details about specific fields or methods demonstrated here. Once you understand the process and get comfortable with it, youll find that
you can measure a nights set of images in a very short order. If the target is cooperative,
period analysis goes just as quickly, maybe even faster.
94
2.
This tutorial measures the second night of images of 771 Libera. The process is nearly
identical to what was described in the previous tutorial in this chapter. The main difference is that AutoMatch is used in order to demonstrate its use and also because the field
of view jumped around quite a bit the second night and you would have to reset the apertures a number of times.
Even when having to reset the measuring apertures when not using AutoMatch, it
would take much less time to go through this tutorial. AutoMatch is handy but definitely not for the impatient or those with limited time.
1.
If necessary, open the configuration and select the profile created in the second tutorial, the one used for looking at DerivedMags values. If you forget, it is not important as long as the profile is one of the two created in the first two tutorials.
You can switch just before doing the analysis. The main point is to make sure that
Canopus had the right information so that it can AutoMatch images.
2.
Open \MPO\Examples\LTCURVES\A771\SEP22\A771001.FIT.
3.
4.
5.
Run the lightcurve wizard, making sure to use the Comp Star Selector.
6.
Using star subtraction is optional (it wont be needed). However, you should
make a habit of checking that box on the first page of the wizard so that you have
the option of using it should a moving target come close to field stars. See Lesson
1 in this chapter and the Reference Manual for additional information.
Use \MPO\Examples\LTCURVES\A771\SEP22\A7710001.FIT for the first image in the wizard and \MPO\Examples\LTCURVES\A771\SEP22\A7710045.FIT
for the second image.
95
9.
Get in the habit transferring data to the sessions form and saving the CSS data to a
text file. Its easier to delete a text file no longer needed that to go through the steps
to recreate it.
10.
Save the Image 1 comps reference image when you get to the last page of the
wizard.
11.
12.
13.
Canopus should get through all the images without stopping to have you reset apertures.
Yes, thats all there is to implementing AutoMatch, unless it cant find the comps. If you
had selected Full Auto (AutoMatch Forced), then even if there is a problem, the process will not stop but a memo form appears after all images have been processed that displays which images failed and why.
14.
Go to the Lightcurve Analysis page and plot the raw data (check the Raw box)
for this session only.
15.
96
17.
Set the period search values to the settings shown in the screen shot immediately
above.
Experience teaches you a likely period search range for an initial search. For one,
asteroid lightcurves with an amplitude > 0.3 mag are almost always bimodal,
meaning that they have two maximums and minimums per rotation (think of a spinning potato). The raw lightcurve covered almost a full bimodal cycle over about five
hours. So, we made sure that a period between 5-7 hours was covered.
Click <Find> again, this time selecting this session and the one created in the second tutorial, the one that also was the result of using AutoMatch.
This shot shows the results of the period search when running this tutorial. The match is
close and the good news is that the two sessions agree very closely vertically. You may
see something a little different.
Remember that the error in the MPOSC3 catalog is ~ 0.03 mag for R. In addition,
were using Clear filter and no transforms. It wont be unusual to see slight displacements, but they should be 0.05 or less if you make sure to use comparisons of the same
color (same color index) as the target.
97
98
3.
Many books and papers can be (and have been) written on period analysis. Youll find a
chapter on it in my book, A Practical Guide to Lightcurve Photometry and Analysis
(Springer, 2006; available on Amazon and Barnes&Noble and may other sites).
Canopus uses the FALC (Fourier Analysis for Lightcurves) algorithm developed by Alan
Harris. The algorithm in the industry standard for asteroid period analysis. Petr Pravec
modified it to do both additive and non-additive period analysis. The former works for
finding the orbital period of a binary system and the rotation period of the satellite. Tumbling asteroids have complex curves that can be analyzed only by the non-additive routines such as Pravecs. Canopus includes the additive dual period analysis. There is a
tutorial on finding the orbital period of a binary asteroid in this Guide.
For now, we will concentrate on a simple case where there is only one period, the curve
is relatively symmetric, the data are of high quality, the amplitude is fairly large, and the
period is reasonably short. If you work with an asteroid with such characteristics, consider yourself lucky. Of the more than 600 that Ive done to date, less than 100 fit that
entire set of characteristics.
This tutorial presumes that you have run the previous three in this chapter and that
you have TrueMags data for two nights on 771 Libera.
1.
Open the sessions form and select one of the two sessions that contain the DerivedMags data generated in Lessons 1 and 2.
2.
Go to the Lightcurve Analysis page (<Ctrl+4>) and enter the search parameters
shown in the screen shot below.
See the Reference Manual for a more detailed explanation of the relation between
the search type, Min, Size, and Steps fields.
99
The Analysis
3.
Click <Find>. The result from running this tutorial is shown below. You may get
something else. The next step explains why.
4.
Note the small form that also appeared. This is the Comp Adjust form and will
be covered in a moment. Another form also appeared.
This is the Period Spectrum form. It has two tabs. The first is a plot of the results of the period search. The period (hours or days) is the X-axis and the RMS
error (in units of 0.01 mag) is the Y-axis.
A lower RMS error indicates a more likely period solution, with some caveats.
For example, note how many periods have low RMS values that are extremely
similar. Any one of those might be the true period and the program could latch
onto wrong one very easily.
The second tab displays the data in a table and shows three error values:
100
PE (Probable Error):
3-sig Error:
This is the error that would cause the last data point in
the set (by date order) to be shifted by 2% (or, an error of
about 7 rotation). It uses a formula similar to the one
used to estimate the period error when doing spin and
shape axis modeling.
What is Right?
Asteroid lightcurve analysis has been called a bit of science and a bit of black magic.
There are many possibilities to consider. We can cover only the basics here.
First, asteroid lightcurves tend to be bimodal, meaning that they have two maximums and
two minimums per rotation. If you picture a spinning potato (Spuds in Space!), you can
see why this would be the case assuming that you were seeing the potato with its spin
axis nearly perpendicular to the line of sight. However, if youre looking down on one of
the potato poles, then you wont see a bimodal curve but more likely one that is
monomodal and has a very low amplitude.
If the asteroid is at a very high phase angle, all bets are off. You can have very complex
curves with three or more maximums and minimums.
As a general rule of thumb, if the phase angle is low, the asteroid is near opposition, and
the amplitude is < 0.10 mag, then you probably cannot say with certainty whether or not
a monomodal or bimodal solution, the latter having a period about double that of the
monomodal curve, is correct. If the amplitude is 0.10 0.20 mag, then the tendency is
towards accepting a bimodal solution. As the amplitude increases, it is more likely that a
bimodal solution is correct, especially with amplitudes of more than 0.3 mag.
If you look at the period spectrum plot, you see a number of solutions between 2.5 3.25
h and another set centered on about 5.7 h. That first group represents monomodal solutions, which have a period about one-half that of the bimodal solution. This fact can be
put to good use, as will be shown later.
Remember, these are all rules of thumb and you could end up all thumbs if you accept them as absolute. Fans of a certain series of movies will recognize the line,
We consider them to be more of a set of guidelines.
5.
In the first solution, a period of about 5.7 h was found with a bimodal lightcurve.
Since the amplitude of the curve is more than 0.3 mag, we will presume that a bimodal solution is correct and what needs to be done is refine the search.
6.
Reset the search parameters as shown here. Note that the period method has been
switched to Period (reg).
101
When the error in the period (shown at the bottom of the plot) is equal to the step
size, it is not a bad idea to refine the search a little more by setting the minimum
period and a smaller step size that centers the new search on the period that was
just found.
9.
Set Size to 0.0001 h and the Min to 5.880 h. Leave Steps at 300.
10.
Check the Fourier box. This is available only if you do not have the Raw box
checked.
11.
12.
This is about the best you can do with the available data. If you want to experiment a little, try changing to higher and lower values for Orders. Youll see that
Orders = 8 probably gives a slightly better fit to the Fourier curve and higher precision.
Do not close Canopus just yet. Move immediately to the next tutorial to explore
using Instrumental magnitudes and the Comp Adjust form.
102
4.
This tutorial covers using the Instrumental method and the Comp Adjust form. Even
when using DerivedMags or TrueMags, the data from various sessions may not agree
perfectly. There are several reasons. For example, there is the inherent error in the catalog
magnitudes or the M/IR solution is not exactly the same for every image. Another may be
color differences between the bulk of the stars used to find the M/IR and the target.
If you recall, Canopus corrects data from a given session to make it directly comparable
to the first session. These corrections are based on using a specific value for the phase
slope parameter, G. In most cases, G = 0.15 is used. This is not always correct and under
extreme circumstances, one can get a better fit of the data by changing the assumed value
of G and recalculating the corrections. The sessions form stores the value for G to be
used to compute the corrections and so you can change that value and see if it makes for a
better fit.
All sessions involving the same target must use the same value of G. If you change
one, you must change them all to agree with the first session for that target.
The Comp Adjust form is used to change the DeltaComp values in sessions on-the-fly
and visually so that you can immediately see the effect of changing the value for one
session on the overall solution.
1.
Open the configuration form and select the profile that you created in Lesson 1 of
this chapter and change from the Photometry Magnitudes | Method from Derived to Instrumental.
2.
Go to the Lightcurve Analysis page (<Ctrl+4>) and set the search parameters as
shown below. Make sure that the Fourier box is not checked.
Were going to cheat by using parameters that are close to the final solution since
we have a good idea of that final solution. In real life, you would start with the
auto search as we did in the previous tutorial.
3.
Click <Find>.
4.
Select the sessions you created in the first tutorials of this chapter.
103
5.
Scroll down until you find the CSS data and note the average R magnitude
for the comparisons. In this case, it is 12.850.
Enter that value in the DeltaComp field.
Click <Save>.
Highlight the second session and find/enter its average R magnitude (13.453
when running this tutorial)
6.
104
Click <OK> to close the sessions form after you have saved the changes for the
second session.
An Alternate Method
You may have noticed by now that the average of the comparison magnitudes is immediately available on the Sessions form; its shown just below the five comp star values.
If you manually change any of the values or change the status of the Use box for a
comparison, the average is updated immediately.
Instead of opening the notes form, you can simply enter the average value into the DeltaComp fields.
If you should enter values for a comparison star directly, e.g., after finding more
accurate values or finding them from another source, make sure you enter 0.000 for
any star that is not used at all. Then, even if you check the Use box for the star, it
does not affect the average. However, it still affects the calculations for the object
during period analysis if Use is checked.
Change the search parameters so that the search is about centered on that value.
8.
105
The spinner control next to the DeltaComp field can be used to change the value in that
field, or you can enter a new value directly. The spinner is easier since you can use the
mouse all the time.
The legend just above the list of sessions tells how to change the amount the DeltaComp
changes for each mouse click. Clicking on the spinner changes DeltaComp by 0.01 mag.
<Shift+Click> causes DeltaComp to change by 0.1 mag while <Ctrl+Click> causes DeltaComp to change by 0.001 mag.
A given session always moves in the direction of the arrow, regardless if using TrueMags
or instrumental magnitudes for analysis. So, if you click <Up>, all the data in the selected
session move up.
You must click <Replot> to have the change take affect.
Restoring Values
You can restore original values as long as you do not click <Find> again. Right-click
over the sessions list and choose from the popup menu to restore the highlighted or all
DeltaComps to their original values. Canopus automatically replots the data with the
restored values you do not have to click <Replot>.
106
10.
The second session appears to be a little too high against the first session. So,
highlight the second session in the Comp Adjust list and click the down arrow 3
or 4 times (do not hold down the Shift or Ctrl key).
Click <Replot>.
Now the two sessions are almost exactly on top of one another and the period solution has a much better precision.
You are now a steely-eyed lightcurve analyst with honors!
In Conclusion
An adjustment of < 0.05 mag (in other terms, < 0.05 mag agreement between nightly zero
points) is very good. If you get < 0.03 mag shifts (not uncommon) then the results are
excellent. Most people doing all-sky photometry under average conditions cant do much
better, at least consistently. Of course, there are those who achieve much better results
but usually in above average locations and only with very careful work. Our goal here is
not to derive magnitudes for high-precision work but to get good enough agreement
among sessions for proper period analysis.
107
108
Using PhotoRed
Using PhotoRed
PhotoRed is a utility program that helps you perform a variety of special photometric
calculations:
Determination of your systems transformation coefficients, which convert your
instrumental magnitudes to the standard Johnson-Cousins BVRI or Sloan (SDSS)
gri systems.
Determination of first and second order extinction terms.
Determination of color indices.
Determination of the nightly zero point.
Conversion of instrumental magnitudes to exoatmospheric standard magnitudes.
Automated observations of variable stars for submission to the AAVSO.
A detailed discussion of photometry and the derivation of the routines to derive the transformation values is beyond the scope of this Guide and even the Reference Manual.
Readers should refer to the several excellent references listed in the Credits section of the
Canopus/PhotoRed Reference Manual.
A Diminishing Need
PhotoRed was written in large part to provide the necessary steps to convert instrumental
magnitudes taken during time-series work to standard magnitudes. By reducing to standard magnitudes, results can be directly related to other such data and make it easier for
some studies since all data would be on a common, well-defined system.
With the TrueMags method in Canopus, all of this almost becomes moot when using
DerivedMags or TrueMags, especially if using the SDSS magnitudes in MPOSC3. However, there are instances where a more detailed and rigorous process is necessary.
Since the TrueMags method does not do differential photometry, it is still subject
to variations due to changing conditions (more so than, for example, changing extinction). When differential photometry using instrumental magnitudes produces
better results, then if the observer requires the data can still be put on a standard system.
Until an all-sky, well-calibrated catalog is available, the DerivedMags and TrueMags values can still have some significant errors. That catalog is coming but its
not here yet.
For DerivedMags and AAVSO batch processing, you can apply color corrections
(based on the color index of the target and a given comparison star). These improve the quality of the data by making them conform better to the standard system and help remove scatter (errors) caused by color differences.
The need for PhotoRed may not be as great as before, but there is still enough need such
that it is worth your time to work through some quick tutorials if for no other reason than
to understand the basics of calibrated photometry.
109
Using PhotoRed
110
1.
Theres no point in repeating steps when not necessary. This is best demonstrated by the
fact that the Transforms and Modified Hardie First Order Extinction (MHFOE) methods
can be run using the data obtained from a single combined set of images taken of two
reference fields.
As covered in the opening sections of the PhotoRed Reference Manual, the MHFOE
method relies on having transforms. This is backwards to most peoples thinking but it
really does work. The transforms value for a given filter is the slope of a linear regression
using the instrumental magnitudes in that given filter versus the standard (catalog) color
indices of the stars that were measured. It turns out that the transform values (slopes) are
independent of first order extinction. If you change the FOE values when finding transforms, the only thing affected is the zero point of the linear regression.
Furthermore, the only way the MHFOE method can work is to use the transforms against
standard color index values for the stars being used in the analysis. This converts the
instrumental magnitudes close to the standard system, on which Hardies original method
relied. The MHFOE method provides poor results at best if this is not true.
The above are very important concepts in PhotoRed. Understanding them makes understanding and using PhotoRed much easier. See the First Order Extinctions Is This Trip
Necessary? and the following section in the PhotoRed Reference Manual for a detailed
explanation.
The Modified Hardie First Order Extinction (MHFOE) method is preferred over the
comparison star method for several reasons.
It takes only a few minutes, thus leaving more time for observing the program field.
It avoids problems with changing photometric conditions through the night.
It uses standard (or near standard) stars, so the catalog values are well known.
Getting Started
In this tutorial, youll use some images taken of a field at air mass 1.06 and another at
2.0.
1.
Open the configuration form in Canopus and match the settings on the screen
shots below. PhotoRed does not allow you to create or change the profile setting,
so you must make sure you have the settings entered in Canopus first. The reason
for this restriction is because more than one instance of Canopus can be started at
a time. The currently selected profile name in Canopus is used to keep track of
settings among the multiple instances.
Note that youre creating a new profile name in this tutorial. You wont necessarily
create a new profile name each time you run the program. Youre creating this one
so that the settings match the system used to take the images youre about to measure. If you have several setups, creating a separate profile makes using Canopus
and PhotoRed much easier. Instead of changing values on one or more page of the
configuration form when you measure images taken with a given system, you can
simply select a profile of preset values.
111
Turn off (uncheck) AutoSave raw data on the Save Options subtab on the Photometry page.
2.
Once you have matched the settings, click <OK> to save the new configuration.
3.
Measuring Field 1
4.
Set the apertures to 11x11/2/11 for target, comps, and astrometry. See Setting the
Apertures tutorial in the Core Operations chapter.
5.
Select Image | Open from the main menu and open the file
112
Select Image | Auto match/measure from the main menu (<Ctrl+A>). The program matches the chart to the image.
7.
8.
Select Yes. This displays the first page of the Transforms Wizard.
9.
10.
Click <Next>. This displays another page of explanatory text. After youve read
the text, click <Next>. This displays the USER catalog embedded in the Photometry Wizard. The goal from here is to measure the several USER stars in the image.
Use the screen shot below as a guide for the new few steps.
113
Youre about to have the PhotoRed form, the Wizard, and one other form on the
screen. Arrange things as you need. Use the zoom box. It helps seeing faint companions near potential candidate stars.
11.
<Ctrl+Click> on the star in the image that has apertures around it.
12.
Note that in the Find field of the Wizard, the table has been automatically positioned to the star and that its name is in the Find field.
13.
Click <Add> on the Photometry Wizard. The data for the star appears in the list
on the Wizard.
14.
Add more stars to the list by repeating the three previous steps, using a different
star each time, until you have at least five, preferably more, stars. If you click on a
star that is not in the UserStar catalog, the program beeps. If you select a star that
is already in the list and then click <Add>, the program displays a message and
does not add the star again.
When selecting stars, avoid those that have close companions and/or are very faint.
This helps avoid bad data points later on.
The screen shot above shows the wizard after a number of stars have been added.
114
Click <Next> to display the final page of the wizard and then click <Finish>. This
displays a file selection list.
16.
Select all the images for AK GEM and then click <Open>. This displays a confirmation message.
17.
Click <Yes>.
18.
19.
Since the images were taken so that all of those of one filter were taken back-toback, the selected images are grouped by filter regardless of your choice. However,
this will not always be the case and which option you chose can have a profound affect on the final results. See the Reference Manual for additional information.
Double-click on the first image in the list. You should see something like the
screen shot above for the image.
You may have selected different stars, so the number and placement of apertures
may be different. The important point is that you see the anchor star surrounded
by green (lime) apertures.
Set the controls as shown in the screen shot above.
115
By setting Group to 1 and selecting Increment | None, all the images will be
put into the same group, i.e., Group 1. See the Reference Manual for additional information about groups.
The Read Filter feature (activated by checking the Read Filter box) requires that
the FILTER key word be used in the FITS header and that the values assigned to
that key word use the standard names for BVRIC filters, i.e., B(lue), V(isual), R(ed),
I(nfrared), and C(lear). For SDSS filters, the entries need to be g, r, or i. Your
imaging software must follow these conventions when writing data to the FITS
header. MPO Connections does include the FILTER key word and allows the names
to be mapped to meet the naming requirements.
20.
The apertures should be centered over the stars. If they are not, click on the anchor star. This should reposition all the apertures so that they are centered over
their respective stars.
21.
Click <Accept>. This measures the image and automatically loads the next image
in the list, which is the second V image.
22.
For each successive image, confirm that the apertures are centered after the image
is loaded. If not, reposition them. Confirm that the Filter radio button setting
matches the filter for the image (the name of the sample files contains the filter).
Then click <Accept>.
23.
24.
Click <Close>.
25.
The file save box appears as soon as you close the wizard. This helps assure that
you do not forget to save your data. Save the file under the name
PHOTORED_LESSON1.OBS
Of course, you can select about any name you want but future lessons will refer to this
specific filename.
Measuring Field 2
26.
Select Image | Open from the PhotoRed main menu and open the file
\MPO\EXAMPLES\PHOTORED_FOE\ LW_CAS_REF_C_00001.FIT
27.
Set the apertures to 9x9/2/11. These images have much better seeing and the
available stars have closer companions that can be excluded by slightly smaller
apertures.
28.
29.
Start the Photometry Wizard as you did for the first image. This displays the confirmation message.
30.
Click <No>. You want to add the observations youre about to measure to those
from the first field. If you click <Yes>, the observations from the first field are
deleted.
116
31.
<Ctrl+Click> on the star on the image that has apertures around it. This is your
anchor star for the second set of images.
32.
As you did with the first image, select an additional set of stars with no less than
five stars and preferably 10 or more.
33.
When done defining the stars to be used, click <Next> to the last page of the Wizard and then click <Finish>. This presents the file dialog.
34.
35.
Use the screen shot below for the next few steps.
117
36.
Double-click on the first file in the list. This loads the image and displays the
measuring apertures. Again, you may have selected different stars. Make sure that
the anchor star has the green (lime) apertures.
37.
38.
As you did for the first set of images, as each image is loaded confirm that the
measuring apertures are centered and reposition them if necessary. Also confirm
that the Filter radio button matches the filter for the image (the filter is included
in the file name).
39.
Click <Accept> to record the data and move to the next image.
40.
41.
Click <Close>.
42.
The file save box comes up as soon as you close the wizard. Save the file under
the same name you previously saved the data, i.e.,
PHOTORED_LESSON1.OBS
PhotoRed warns that you are about to overwrite the file. Click <Yes> to overwrite the
file, which now contains the data for both images.
You can take a break here if you want. However, youre just a few seconds away from
seeing the results of your efforts so far. You may have a little work to do, e.g., removing
some bad stars from the solution set, but that wont take long either.
Transforms First
Contrary to some conventions and for the reasons discussed at the start of this lesson,
well find the transforms first. The assumption is that you have not previously measured
images for the night and so must use assumed first order extinction values when finding
the transforms.
43.
118
From the PhotoRed main menu, select Photometry | Transforms values (or
<Ctrl+T>). This opens the transforms form.
44.
Go to the Extinction page and enter the values shown above for the First Order
V, R, and C fields. These are just estimates. Dont worry about the other values
on the tab.
When entering assumed values, the V and C extinctions should be similar while the
R should be a little less than V, so that V-R is a small positive number on the order
of 0.05mag.
45.
46.
Remember that your data may not be exactly the same since you may have picked
different reference stars.
47.
Click on the Grp column header so that the data is sorted by group number. If
Group 2 is at top, click the column header again.
48.
Make sure you are scrolled all the way to the top of the data list and then click on
the first observation.
49.
Without clicking on the data list, scroll down until you see the first observation in
Group 2.
119
<Shift+Click> on the last observation in Group 1. This should highlight all observations in Group 1, which should be the lower field and have an air mass value
near 2.0.
51.
Right-click over the data area and select Set Use to False from the popup menu.
You should see the entries in the Use column of the table change to False for
the selected observations.
52.
Select the Transforms method under Reduction Method and check the boxes
for Visual, Red, and Clear.
53.
Click <Reduce>. This displays a form asking which color index to use.
54.
Normally, you would think that because you took images in V and R that you would
automatically select V-R for the Standard Color Index, and youd be right. However, strange as it seems, PhotoRed allows finding a standard color index different
from the filters used, say B-V versus v-r. This mixing of color indices is not generally recommended because of secondary extinctions and other factors. However,
tests have shown that good to excellent results judged by comparisons to reduced
magnitudes of stars in reference fields can be obtained. See the PhotoRed Reference Manual for additional information.
Very Important! Even though you can select an off-color index, there must still be
catalogs magnitudes for the filters for which you are trying to find transforms. For example, the program uses r-R (instrumental minus catalog) magnitude when finding the R
transform. It cannot use B, or V, or any other magnitude to find the transform for R.
55.
56.
120
When transforming Clear observations, you can have them referenced to any one of
the supported filters. Usually, this is V and so PhotoRed would find the transform
that converts a Clear instrumental magnitude to standard V. Regardless of your
choice, and as noted above, there must be catalog values for which ever filter you
chose.
57.
58.
Using the data measured for this tutorial, the standard deviation for the linear regression line is 0.018 mag.
59.
Click on the R and C tabs at upper right to display the results for those filters.
121
If youre satisfied with the results, click <Save> on the Transforms form and then
click <Close>.
62.
Make sure that you save the observations data again if you made any changes,
e.g., marked the Use flag some stars to False.
k
0.200
0.150
0.200
Tx(V-R)
-0.132
-0.112
+0.344
ZP
21.903
21.833
22.718
At the end of the lesson, youll see how the Tx value (transform for a given filter) does
not change (or very little) but the ZPx (zero point value) does when actual first order extinction values are used instead of the assumed values.
64.
If you changed the sort order of the data for some reason, click on the Grp column header until the observations for Group 1 are at the top of the data list.
65.
Select all the observations for Group 1 and use the popup menu to change the
Use column back to True.
66.
Select First Order (Hardie) in the Reduction method and confirm that only the
Visual, Red, and Clear boxes are checked.
67.
Click <Reduce>.
68.
Select V-R.
This method uses only catalog magnitudes and so it is not necessary to select
which instrumental color index to use.
122
You must have observations in the two filters of the selected color index and there
must be catalog values for those filters for this method to work.
69.
Click <OK>. You are asked if you want to set the nightly zero points.
70.
71.
As before, youre asked which standard filter magnitudes to use for the Clear filter. Click <V> for Visual.
Remember that finding the transforms earlier also found nightly zero points based
on the assumed FOE values. Since we are finding the true FOE values, which in
this case are not the same as those assumed values, we must accept the new
nightly zero points. Otherwise, all the reduced standard magnitudes will be in error.
See Appendix A of this Guide where this lesson is repeated and you do not set the
nightly zero points.
Again, your results may be different. Note the tight grouping of the data at each end. The
group on the right is the data from the high air mass (lower) field while the data at upper
left is from the low air mass (higher) field. If the data are spread widely at either end, you
may have some bad data points and need to remove them from the calculations as described above.
72.
The standard deviations in the derived values were on the order of 0.02 mag,
which is very good for my location.
The kV-R was +0.063. This is good because it means that the V first order extinction is greater than for R and that the two are similar. Remember that this value
should be on the order of +0.05.
73.
74.
If satisfied with the results, click <Save> on the Transforms form and then close
that form.
123
124
2.
Transforms Errors
This lesson assumes you did the previous tutorial and saved the observations data from
that tutorial. Here, well presume you saved the data under the name
\MPO\UDATA\PHOTORED_LESSON1.OBS
The point of this tutorial is to allow you to review the quality of the transforms and zero
points you found earlier by using them against the original data. In a perfect world, using
those transforms and zero points against the original data should produce standard magnitudes that exactly match the catalog values.
You dont have to run this routine every time you find transforms but its a good idea to
do spot checks from time to time, just to be sure your results are valid. If nothing else,
this routine helps spot a bad star that might be adversely affecting the transforms values.
You wont be measuring any images in this lesson, so you dont need to worry about
configuration settings.
1.
If you closed PhotoRed after the previous tutorial, you need to open PhotoRed
and load the data. If youve moved directly here from the previous lesson, you can
skip to step 5.
2.
To load a PhotoRed OBS file, select Photometry | Load observations from the
PhotoRed main menu. This opens a Windows file dialog. Locate and open the appropriate file.
3.
Remember that the data are for two fields, one high and one low, and that you
found the transforms using only the higher field.
4.
Click on the Grp column header to sort the data by group number.
5.
Select all the data from Group 1 and use the popup menu to set the Use flag to
False.
6.
Set the Reduction Method and select the filters as shown in the screen shot above.
125
8.
Select V-R since this is color index that you used to find the transforms.
9.
Click <OK>. This displays the form to select to which filter to convert Clear magnitudes.
10.
Select V since that is what you used when finding the transforms.
In this and other reduction methods when you are offered a choice of color index or
Clear to X filter, you must select the same values that you did when finding the
transforms. It will not do to find transforms in terms of B-V and later select V-R.
11.
Remember that your results may be different. The point here is to check the mean
value it should be 0.000 or very close to it and the standard deviation of the
errors. In this case, it is 0.021 mag for the V filter (R had 0.014 mag and Clear
had 0.019 mag standard deviations).
12.
Check Again
Its not really valid to run the errors test on the same data used to find the transforms. Try
running the test on the data for the lower field by setting the Use flag to True for the
Group 1 stars only. Here are the filters, mean errors, and standard deviations found when
running this tutorial
V: 0.006 0.021
R: 0.009 0.015
Given the large air mass of the lower field, these results are very good.
126
3.
One of two methods of getting data into PhotoRed is outlined here. In either case, the
general approach requires that you follow a target field over a large range of air mass
values in order to determine first order extinction (FOE from here on; Im tired of typing
first order extinction). Changing conditions for all but the best sites often makes getting
a usable value difficult.
The Modified Hardie First Order Extinction (MHFOE) method, covered in the previous
lesson, is preferred for finding first order extinction values. The comparison star method
requires that conditions be not only stable but photometric quality for the duration of
the run. This occurs very rarely for most sites, especially backyard observatories.
Open the configuration form in Canopus and match the settings in the screen
shots below.
127
2.
Switch back to PhotoRed and, if you want, use Canopus | Hide Canopus on the
PhotoRed main menu to hide the Canopus form.
Use the blinker in Canopus with the following two images to determine which
star is the asteroid.
\MPO\EXAMPLES\PHOTORED_ASTEROID\A36_C_0001.FIT
\MPO\EXAMPLES\PHOTORED_ASTEROID\A36_C_0124.FIT
128
5.
6.
Click <Yes>. This displays the first page of the Differential Photometry Wizard.
This wizard works identically to the one you used in the Photometry chapter to set up
measuring images to find an asteroids lightcurve in Canopus.
7.
Load \MPO\EXAMPLES\PHOTORED_ASTEROID\A36_C_0001.FIT
8.
Click <Next> until you get to the page to select the comparisons for Image 1.
Use the reference image above for the next few steps.
129
Set the two comparison stars and target on the wizard by clicking on each one and
then the appropriate button on the page. You should see something like this:
10.
Click <Next> to go to the next page and load the second image,
\MPO\EXAMPLES\PHOTORED_ASTEROID\A36_C_0124.FIT
11.
Click <Next> until you get to the page to set the comparisons and target on the
second image.
12.
13.
Click <Next> until you get to the end of the wizard and then click <Finish> to
display a Windows file dialog.
14.
Select all the Clear (A36_C_*.FIT) files except A36_C_0124 and A36_C_0127.
The asteroid is too close to a field star for those images to be used.
15.
130
Depending on how the file names are created by your imaging program, the images you select for measuring in PhotoRed may not always be in chronological
order or grouped by filter (if the program includes the filter in the file name). This
option allows you to try to sort the files in the list so that they make measuring as
easy as possible by not having to switch filters on the Images List about to come
up. Experience will tell you how to answer this question for a given situation.
16.
17.
Check the Increment | List Index radio button. This automatically assigns the
Group number to the list index + 1 (the list index starts with 0. Since 0 is not a
valid group number, PhotoRed automatically adds 1 to the list index).
18.
Do not check Read Filter. The images do not have the FILTER key word in the
header and so PhotoRed cannot automatically make the filter selection.
19.
AutoMatch and auto-processing are not available on the PhotoRed Image List.
20.
Double-click on the first image in the list. This loads the image and automatically
measures the comps and target, as shown in the screen shot above. Note that the
Group box now shows 1.
21.
If the apertures are centered on the appropriate objects, click <Accept> (or press
<Enter> if the Image List form has focus). This records the data for the observation and automatically loads the next image.
131
Proceed carefully in the next step. If you make a mistake and click <Accept> with
the apertures in the wrong place, you can still double-click on the mismeasured image in the list, position the apertures, and then click <Accept>. However, unlike
Canopus, PhotoRed always appends data and so it will not overwrite the bad data
with good. When you go to reduce the data, you will need to delete the offending observation or set Use to False.
22.
Repeat until all the images are measured. If the apertures are not properly centered, click on the anchor star (A in the reference image) to reposition them.
23.
After you have measured the last image, click <Close>. This causes a file save
form to appear. Save the file so you dont lose your data. For this tutorial, use the
file name
\MPO\UDATA\PHOTORED_LESSON3.OBS
The data from this lesson are not required for future lessons but you should get in
the habit of saving your observations data. The one time you dont is the one time
you really should have.
24.
Select Photometry | Open image list from the PhotoRed main menu. This presents a message.
25.
Click <OK>. Were using the same comps and target but a different set of images.
Another message appears.
26.
Click <Yes>. The PhotoRed Reference Manual gives more information on the
available choices and when youd use each one. The Image Selection list is displayed as before.
27.
28.
Sort the files by Date (click <Yes> on the message that appears). This displays the
PhotoRed Image List.
132
29.
Select List Index and set the R radio button. Do not check Read Filter.
30.
As you did before, measure all the images, making sure to reset the apertures if
necessary.
31.
32.
33.
Once youve measured all the V images, save the observations again into
\MPO\UDATA\PHOTORED_LESSON3.OBS.
34.
133
Click the Grp column header to sort the observations by group number until
Group 1 is at the top of the list.
36.
Select First Order (Comps) and check the Visual, Red, and Clear boxes, then
click <Reduce>. This displays a small form that allows you to choose which
comparison star to use for the calculations. The default is 1.
You can enter a number from 1 to 5 but, of course, there must be a corresponding
comparison star. For example, if you have only two comparisons, as in this case,
only 1 or 2 is a valid entry.
37.
Do not click <Save> on the Transforms form. These are very bad results. Youll
see why in a moment.
38.
Use the tabs at the upper right of the plot to switch to those for the R and C filters.
39.
40.
Click <No>.
134
135
136
4.
This lesson teaches you how to import data from Canopus for additional processing in
PhotoRed, usually to convert Canopus instrumental magnitudes to standard magnitudes
on the Johnson-Cousins or Sloan system.
Given the DerivedMags method available in Canopus, especially when the AAVSOs
APASS catalog is fully on-line, PhotoReds future lies mostly in finding the color index of unknown objects and batch processing of variable star observations for the
AAVSO.
For this tutorial, you will use the data from the first lesson of the Photometry chapter
where you measured images for 771 Libera to find its lightcurve. Since instrumental
magnitudes are always measured and saved in Canopus, the data from any one of those
lessons can be used.
1.
Select Utilities | Import Canopus session | from PHSESS from the PhotoRed
main menu. If there are observations on the data page, a confirmation message is
displayed.
2.
This is a multiselect form, meaning that you can select more than one session.
This would make sense only if the sessions were from the same night. They might
be with different filters or the result of having to create several sessions because
of a fast moving object.
In this case, were going to select one and only one session.
3.
Locate the session Sept. 18, 1999, session for 771 Libera and click on it. This
highlights the session (the session number in the Sess column may not be the
same as shown above).
137
Click <OK>. The cursor changes to an hourglass (or whatever your Wait cursor is) while the data are imported. This may be so fast that you dont see the cursor change.
5.
6.
Select First Order (Comps) and check the Clear filter box.
7.
Click <Reduce> and select Comp 1 in the form that appears. This displays a plot
similar to the one below.
Note that the comp star does what it should, get brighter with lower air mass. The
range of air masses is also much larger than before, going from nearly 1.0 to just
over 2.0. This helps assure a good liner regression solution.
Something that should be quickly apparent is that if you follow a field all night, going from eastern horizon to western, is that the preponderance of data will be at
lower air masses. The air mass changes quickly when near the horizon but more
slowly when the field is near the meridian. If you want a good range of air masses
for this method, make sure that you get a number of images when the field is only
30-40 high.
A later lesson will show you how to export the data from PhotoRed to files that can be
read by Canopus.
138
5.
Second order extinction is a term that is dependent on both air mass and color. Youre
probably familiar with the concept of how objects near the horizon are redder than when
seen overhead. This change in color (due to scattering) affects the correction to the raw
instrumental magnitudes to make them exoatmospheric, i.e., as seen from outside the
Earths atmosphere.
By definition, second order extinction (SOE) is zero for the V and R filters. Most tests
have indicated that this is true. However, you should confirm it for your system and, if
for some reason it is not, try to figure out why and if necessary apply it as needed.
The B and Clear filters do have second order terms that cannot be ignored. Keeping all
comparisons and the target (and check) similar in color can reduce the correction to insignificant amounts but if youre finding the standard magnitudes of comparisons to generate your own sequence, then you should consider second order terms.
The most effective way to determine SOE is to follow a pair of stars over a wide range of
air mass. These stars should be significantly different in color, one being red and the
other blue, in order to give the largest possible difference and so make for a more accurate determination of SOE.
Reference Image
Use the screen shot below for reference in the tutorial.
139
Finding SOE
1.
140
Set the configuration to match the shots below. The only items of actual interest
on the General tab are the Header time information (to determine mid-exposure),
the focal length, and pixel sizes. The name of the profile is not important save to
make it easy for you to know why it was created.
2.
3.
Open SDSS_RB_SA92_G_00017.FIT
4.
5.
6.
Select Extinction/Transforms for the purpose and User catalog in the drop
down list on the first page.
141
7.
Add stars A and B from the reference image to the list of targets.
8.
Finish the wizard. When the file selection dialog appears, select all SDSS*.FIT
files.
9.
142
Set up the Image List as shown above, specifically, make sure to select Index
DIV X and set Mod to 3.
One image in each filter was taken per set. Sets were taken throughout the night in
order to follow the field over a wide range of air mass. You want the images in
the each set to have the same group number. By sorting by date, the first three images in the list are each filter in the first set, the second three in the list are the
three images one in each filter for the second set, and so on.
The DIV operator combined with Mod = 3, means that the first three images will
be in group 1, the second in group 2, etc. because DIV returns the integer part of
X/Y and so, 1 div 3 = 0, 5 div 3 = 1 and likewise down the list.
11.
Double-click on the first image, confirm the apertures are correctly placed, and
click <Accept> (or press <ENTER> if the Image List has focus). This automatically loads the next image. Continue until all images are measured.
12.
Save the observations data set before running the reductions routine.
Proceed slowly. You may have to reset the apertures a number of times during the
measuring.
No Do-overs
Unlike the Canopus Image List, with the PhotoRed Image List you cannot re-measure an
image and have the new data replace the old. If you make a mistake, reload the image
that was mismeasured (double-click on it in the list), reset the apertures, and then click
<Accept>. You will have a bad data point in the set that you can usually spot by the
SNR of the stars being dramatically lower compared to the correctly measured image
or by an outlier data point on the plots.
14.
Check the g, r, and i filter boxes. Do not check the Clear filter box at this time.
15.
Click <Reduce>.
16.
143
For each filter, PhotoRed plots the raw data for the first star in the list. This lets
you check that the overall data appear to be good and that there are no outlier
points throwing off the solution.
17.
Review the plots for the first star for the three filters. They should all be reasonably good. The standard deviation for all three was 0.007 mag when this tutorial
was run.
18.
This tab displays the intermediate and final results. For each star, you see the
slope of the raw data (instrumental magnitude versus air mass), the standard color
index, and the standard deviation of the least squares solution for the raw data.
19.
144
This tab shows the plots for the three filters with the X-axis being the standard
color index of the star and the Y-axis the slope of the raw data least squares
analysis.
What you hope to see is a value of 0.0 for the second order slope. This would
mean that there is no color dependency in extinction for the given filter.
For the g filter, SOE = +0.03. This is a little strange since it implies that red stars
fade faster. That is obviously not the case. When you get a strange result, and
even in general, you should get data from additional nights and confirm what
youre getting and take an average. This may have been an off night that needs
verification. However, that is not the point to this tutorial.
When there are only two data points (stars), as in this case, PhotoRed adds an artificial point that is the average of the two points. This is required in order to get a
least squares solution with standard deviation. The slope (SOE term) does not
change as a result.
145
146
6.
The underlying approach to PhotoRed for finding standard magnitudes for a target relies
on knowing the standard color index of the comparisons and target beforehand. To make
the final step to standard magnitudes for the target, the standard magnitudes of the comparisons must be known. The next two lessons are dedicated to finding those values.
In case your curiosity is up, the final reductions use a differential photometry formula
(based on a paper by Richard Miles of the BAA), that accounts for the difference in the
standard color index of the target and a given comparison star. If more than one comparison is involved, PhotoRed finds the derived standard magnitude differential value between the target and a given comparison and applies the derived standard magnitude for
that comparison. The average of the derived standard magnitudes of the target is then
reported as the final value.
Since Practice makes perfect (better yet Perfect practice makes perfect), youre
going to repeat many of the steps in the first part of a previous lesson except on a different set of images and with a slightly different but critical difference. Youll also learn a
trick that save you lots of time for the next lesson, which is finding the standard magnitude of the comparisons.
Check that you are using the profile created in PhotoRed Lesson 3 (First Order
Extinction Comp Star I). That lesson called the profile PhotoRed FOE. If
you didnt run that lesson, go to it now and create a profile with the settings
shown in that lesson.
2.
3.
147
4.
Set the values on the Extinction and Transforms tabs as shown in the two screen
shots immediately below and then click <Save>.
The values were determined off-line using the same set of images and the modified Hardie method.
5.
6.
7.
Click <Yes> on the confirmation message to clear all observations. The wizard
then appears.
8.
148
9.
For each of the numbered stars (the comparisons), click on the star on the image
and then click the same numbered button. For example, click on star 1 and then
click <Star 1> on the wizard. Do this for all five comparisons.
10.
Set the position for the target by clicking on it on the image and then click
<Target> on the wizard. You should see something similar to this.
11.
The X/Y coordinates may be slightly different in the third or fourth decimal place.
Click <Next> on the wizard and load the second image,
\MPO\EXAMPLES\PHOTORED_M67\M067_C_0003.FIT.
12.
Click <Next> until you get to the page to set the positions on the second image
13.
Set the positions for the comps and target on the second image. Your results
should look similar, but maybe not exactly, as these:
149
14.
Click <Next> once. This takes you to the final page where the button is now
called Finish.
15.
Click <Finish>. This closes the wizard and displays a Windows file dialog.
16.
17.
Click <Yes> to the question about sorting the files by date or name.
Depending on how the file names are created by your imaging program, the images
you select for measuring in PhotoRed may not always be in chronological order or
grouped by filter (if the program includes the filter in the file name). This option allows you to try to sort the files in the list so that they make measuring as easy as
possible by not having to switch filters on the Images List. Experience will tell you
how to answer this question for a given situation.
18.
19.
Check the Increment | Index mod X radio button. The combination of these two
settings allows you to measure all the images without having to change the Group
number manually. As you'll see, each of the three images for a given filter is assigned to a different group.
20.
Check the C radio button (assuming the first image is as shown in the screen
shot below).
21.
Do not check the Read Filter box. These images do not have the FILTER key
word in the header.
22.
Double-click on the first file name in the list. This loads the image into the program. The measuring apertures should be centered on the comparisons and target.
150
23.
Proceed carefully from here on. The following steps assume that the files are in
the order seen in the screen shot.
Measure the three Clear images by clicking <Accept>. Reposition the apertures if
necessary.
Note how the value in the group number automatically increases as each Clear
image is loaded
24.
When the first R (red) image is loaded, do not immediately click <Accept>, even
if the apertures are aligned.
25.
Change the filter to R (not r) and note that the Group field has gone back to 1.
This is as it should be.
When running the tutorial, it was necessary to click on Comp 1 (the anchor star) for
each R image since the apertures were slightly off center. This may or may not be
necessary for you.
26.
27.
Pause again when the first V (visual) image is loaded and set the Filter to V.
151
28.
29.
Click <Close>. This displays a file save form. Save the file as
PHOTORED_COLORINDEX.OBS
It is not necessary to enter the OBS extension. PhotoRed forces the file to have that
extension when the file is saved. This is required since, when trying to load previously saved files, the file dialog looks only for OBS files.
30.
A Stitch in Time
Here comes a time-saving trick. The Color Index (Comps/Target) method requires three
groups of data in order to calculate a standard deviation of the derived values. On the
other hand, the Comps Standard Mags method, which uses the exact same data, must
have the data arranged so that there is one and only one group but with at least three observations in each filter per group to compute the standard deviations. You dont want or
need to remeasure the images, where you would force the group number to 1.
152
Click on the Grp column header on the data list until the data is sorted by group
number with all Group 1 data at the top of the list.
32.
33.
Right-click over the data list and select Set group number from the popup
menu. This displays an entry form.
34.
Enter 1 and click <OK>. All the data should have the group number set to 1.
35.
Select Photometry | Save Observations from the PhotoRed main menu and save
the data under the name PHOTORED_CMPSTDMAG.OBS.
Select Photometry | Load Observations and load the original data that was
saved under PHOTORED_COLORINDEX.OBS
37.
38.
Clear is not a standard filter, so you do not check that box. If you do, PhotoRed
ignores it.
39.
40.
Select V-R since this is the color index you used to find the transforms and
FOE values.
41.
This is important! You must select the same color index that you used in finding the
transforms and first order extinction (if you used the modified Hardie method). You
also must have observations in the two filters involved in the color index and those
two filters must be checked in the Filters section.
Click <OK>. This displays a confirmation message.
153
42.
The results in listed in the Transforms form as well as at the top of the plot. The Sigma
values are the standard deviations of the average V-R for each comparison and the target.
When running this tutorial, all the sigmas were 0.005 mag or less.
43.
Comparing Results.
How good are the results? Heres a comparison of the derived versus Henden catalog
values for the five stars and target.
Object
Comp1
Comp2
Comp3
Comp4
Comp5
Target
154
Henden
V-R
0.249
0.355
0.141
0.532
0.568
0.049
Derived
V-R
0.247
0.366
0.136
0.538
0.561
0.050
Diff.
H-D
+0.002
-0.011
+0.005
-0.006
+0.007
-0.001
7.
The final step before finding target standard magnitudes is to find the standard magnitudes of the comparisons. This is done by combining all the values found so far, i.e.,
transforms nightly zero points, extinction, and standard color indices. Once the standard
magnitudes of the comparisons are known, then the color corrected differential magnitude between each comparison and the target are found and applied against the comparisons standard magnitude to give a reduced standard for the target.
To repeat one more time, you should try to use a reference field close to the target
field, maybe a couple of degrees or less. Also, if you use reference and comparison
stars that are similar in color to the target, then the extinction and color index values are of little or no consequence in the final reductions. The reason youve
learned how to find those values is because you cant always have a reference field
nearby or have the luxury of finding reference and/or comparison stars that sufficiently match the color of the target.
1.
From the PhotoRed main menu, select Photometry | Load observations. Load
the PHOTORED_CMPSTDMAG.OBS that you created in the previous lesson.
2.
3.
Select Comps Standard Mags and check the Visual, Red, and Clear filter boxes.
4.
5.
Select V-R.
6.
Remember that you must use the same color index throughout the process, from
transforms to target standard magnitudes.
Click <OK>. This displays the Plots page and Transforms form.
155
7.
8.
Internally, the standard deviations are less than 0.01m for each comparison. Lets compare results against catalog values again. H stands for the Henden magnitude while P
stands for the PhotoRed derived magnitude.
Obj
Comp1
Comp2
Comp3
Comp4
Comp5
VH
VP
11.546 11.544
11.284 11.290
10.929 10.912
12.123 12.131
11.242 11.250
Err: -0.001
+0.011
RH
RP
CH(v)
CP(v)
Vp-Cp
11.297 11.297 11.546 11.547 -0.003
10.924 10.925 11.284 11.295 -0.005
10.788 10.778 10.929 10.907 +0.005
11.590 11.593 12.123 12.128 +0.003
10.692 10.689 11.242 11.242 +0.008
Err: -0.000
Err: +0.001
+0.002
0.002
0.013
0.005
These results will do for many cases. The standard deviations of for the errors are a little
large; youd prefer them to be more on the order of 0.01 mag and not 0.02 mag. The good
news is that the Clear-to-V reduction shows the magnitudes to be nearly identical to the
directly reduced V magnitudes, which are acceptably close to the Henden values.
156
8.
In the final step, youre going to find the standard magnitude of the target from the last
lesson. The approach of using the data from the Color Index measurements, i.e., a limited
number of images, is not uncommon. As previously mentioned, you might visit a target
field only once or twice in a night to keep tabs on a long period variable or some other
project. In this case, you wouldnt need a protracted time-series and theres no need to
clutter the Canopus sessions data table with a number of sessions with very few observations in each.
Its assumed for this tutorial that you did the lessons in the Photometry chapter
where you measured the images for 771 Libera on two nights.
The results that you get here will be meaningless. This is because there were no images in
V and R of Libera to find the true color indices of target and comparisons or the standard
magnitudes of the comparisons. Once again, the goal is not to reproduce a set of results
exactly but to learn the processes in PhotoRed well enough to use them on your own
images.
1.
Select Photometry | Load Observations from the PhotoRed main menu and load
PHOTORED_COLORINDEX.OBS created in the Color Indices tutorial.
By using the data where the observations were split into three groups, you can get
an average and valid standard deviation for the target. If you were to use the data in
with all the data in Group 1, you would get an average but the reduction routine
could not find a valid standard deviation.
2.
3.
Set the Reduction Method to Target Standard Mags and check the Visual, Red,
and Clear filters.
4.
5.
Click <No>. The errors are going to be above average (or what you hope is average) and so the error bars will only fill the plot needlessly.
6.
As before, select V-R since that is the color index that has been used up to this
point. This displays the V plots page, which should show three data points.
157
7.
Click on R and C tabs at upper left to see the plot for those filters.
8.
Click the Std. Mags tab to display a text file with the data broken out by filter.
9.
You can click <Save> at lower right to open a file save form and save the text
file, which can then be used in a spreadsheet or sent to other observers who might
not have Canopus.
The standard V and R magnitudes from the Henden catalog for the target star are 11.267
and 11.218. PhotoRed found 11.272 (+0.005 mag difference) and 11.219 (+0.001 mag
difference) respectively. The Clear-to-V reduction found 11.265, or 0.002 from the
Henden value and 0.007 mag from the PhotoRed V value.
158
9.
This tutorial is going to use a very important trick within PhotoRed in order to get reasonable results. This trick is good to know if you already know the standard magnitudes
and color indices for the comparison stars, which is often the case when working a variable star field that has comps that have been well-calibrated. Well use an average V-R
value for asteroids since there no images available to determine its true value.
The images of 771 Libera were taken with an entirely different system from those used to
find the transforms, and so the color corrections will not be exact. Dont let all this stop
you from getting comfortable with the process. When you get your own data, then you
can concentrate on the absolute certainty of the results.
Select Photometry | Backup saved transforms from the PhotoRed main menu.
This displays a Windows file dialog. Save the current set of Transforms under any
name you like, e.g., TUTORIAL_TRANSFORMS_1.PRR. The PRR extension is
forced when you save the file name.
You can restore these transforms to the Transforms form (make them the current
Transforms) by selecting Photometry | Restore to saved transforms and loading
the file that you just saved.
What youre going to do is enter known color index and standard magnitude values for
the comparisons. Again, youre working a field where such data are available, you can
skip all the steps used to find the color indices and standard magnitudes for the comparisons and just enter as them as youre about to do. The only step you would have to take,
and only if the data are not available or a reasonable guess is not good enough, is to
find the color index of the target.
2.
159
The values youre entering were found in Canopus by doing an AutoMatch on the
first image of the set and then clicking on the comp stars on the chart. This displays
the Object Info form which includes, among other things, the BVRI magnitudes derived from the 2MASS J-K magnitudes. The inherent errors in those values migrate
through to this exercise and so increase the overall error of each data point.
3.
On the Color Index tab, enter the values shown above in the V-R column.
4.
On the Std Mags tabs, enter the values shown for the V and C columns.
5.
Click <Save> so that the changes become the active set of Transforms and then
click <Close>.
6.
Select Utilities | Import Canopus data | from PHSESS from the PhotoRed main
menu. This displays a confirmation message.
7.
160
Locate the session Sept. 18, 1999, session for 771 Libera and click on it. This
highlights the session and changes the symbol in the far left, fixed column to a
combined arrow and dot. This indicates that not only is the highlighted row the
current row but also has been selected. If the highlighted row has no symbol in
the column, it has not been selected.
9.
Click <OK>. The cursor changes to an hourglass (or whatever your Wait cursor is) while the data are imported.
10.
Go to the Reductions page (<Ctrl+2>). The observations should be available. Remember: it may not be exactly the same as what you have but should be close.
11.
Select Target Standard Mags for the reduction method and check only the Clear
filter box.
12.
13.
14.
You can click <Yes> because you cheated by manually entering comparison star
data and an assumed V-R for the asteroid in the Transforms form.
Select V-R (why?) and click <OK>. This displays the plot.
161
Now you have a real lightcurve with standard magnitudes. These are V magnitudes since the Clear transforms were to the V band. The average magnitude from
this plot is about 13.75. The predicted magnitude based on data from the Minor
Planet Center was 13.81, or 0.06 mag fainter.
Dont necessarily write this off due to errors in the reduction process. That predicted magnitude is based on two values: H (absolute magnitude) and G (phase
slope parameter). A slight error in either one or both could make up a large part of
the different here. Checking the accuracy of the H-G parameters is one of the reasons for reducing data to the standard system.
In this case, however, the inherent errors in the 2MASS conversions and the fact
that the transforms were found using data from a different telescope/camera combination had more to do with the error, which is still not excessive.
15.
Click the Std. Mags tab to review and, if you want, save the text file for experimenting with another program.
Return to Sender
16.
Had this not been a test, you would normally want to export the data back to
Canopus export files (also called saved session files) where you could do further period analysis or more detailed plotting or send them to another Canopus
user.
17.
18.
162
If the data you import from Canopus has only V and/or C measurements that were
converted to V in PhotoRed, then you want to force the session information in the
exported files to indicate that the V filter was used. In this case, click <Yes>.
If you have observations in other filters or the Clear transform converts the Clear
instrumental magnitudes to a band other than V, you should click <No>, otherwise the sessions file carries the wrong filter information.
19.
Click <Yes>. The cursor changes to the wait condition while the data are written to the hard drive. When the export is complete, a Windows select directory
dialog appears.
The default is \MPO\UDATA. You can save the files here or any other directory,
including one that you create on-the-fly. It is recommended that you store all export files in \MPO\UDATA or a subdirectory of it.
20.
Once you have selected an existing directory or entered the path to a new one,
click <OK>. This displays a confirmation message.
21.
Click <OK>.
Whats in a Name?
Note the naming convention used when creating the files.
First is the object name as listed in the Object field of the Canopus Sessions form (up to a
maximum of 30 characters). This is followed by _X_ where X is the filter for the given
session, and then by an assigned number so that you could export additional files, e.g.,
from a second night. For example, the exported files from this tutorial were named
163
The _OBS.FF2 and _SESS.FF2 that end the file name are mandatory. If you
change the name of the files, you must keep these ending characters and maintain
the FF2 extension. Otherwise, Canopus cannot work with the files.
Furthermore, what precedes these two mandatory suffixes must also be identical.
For example, if the OBS file in the files above had V_2 instead of V_1,
Canopus would not be able to import the observations. In short, change file names
with great care and only if you know what youre doing.
Whenever youre working with imported Canopus data, many of the menu items under
Photometry on the PhotoRed main men are disabled. While the data appear to be normal in relation to what youve seen in most of the lessons, there are some underlying
differences important enough to have PhotoRed try to prevent some actions. Before you
can work with PhotoRed generated data, you must clear the Canopus data.
Cleaning Up
22.
From the PhotoRed main menu, select Utilities | Clear Lists. This displays a
confirmation message.
23.
Click <Yes>. This clears the data lists and you can work with PhotoRed generated
data again.
Of course, it does little good to create these files if you dont know how to use them or
pass them around to other Canopus users. Data Importing and Exporting are a vital part
of the Canopus/PhotoRed system and so a separate set of tutorials is devoted to those
tasks.
164
10.
AAVSO Batch Processing allows you to measure a number of images for a non-moving
target automatically and find the standard magnitude of a target based on differential
photometry using pre-known comparison stars. This is the approach taken by the
AAVSO in its CCD reporting format. In previous times, there would be one and only one
comparison, a check star, and the target. That has changed in that more than one comparison can be used and the check stars role has changed when more than one comparison is
used in ensemble photometry.
Before, the role of the check star was to subtract its magnitude from that of the comparison star. The constancy and scatter of that differential was a measure of the quality of the
target observations and to assure that the comparison was not variable.
The new role of the check under ensemble photometry is as a control target. Just as
with the target, the magnitude of the check star is found for each comparison. The average of the derived magnitudes of the check and standard deviation of the average are used
for the check star value in the report. If all is well, the reduced magnitude will match the
adopted magnitude of the check. If the reduced magnitude does not match the adopted
value, then the AAVSO or future researcher can apply the difference to the target magnitude to get the true value. Of course, the derived value for the check should be consistent, even if its not exactly equal to its adopted value. Also, should the adopted magnitude of the check change, then the difference can be factored into the target values.
Process Overview
The process of measuring images for AAVSO work (or any non-moving target) is as
follows.
1.
Take images of the target field in one or more standard filters or Clear.
2.
If youre going to transform the values, i.e., apply color corrections, youll need to
find the transforms for your system. See the PhotoRed tutorials in this Guide and
the PhotoRed Reference Manual for more information.
Since AAVSO Batch Processing relies strictly on differential photometry, the nightly
zero points and first order extinction terms are not critical. What is important is that
the difference in the first order extinction be about right. For example, if relying on
V-R color index, then kvr should be about +0.05.
The only remaining terms of importance are the transforms and second-order corrections. The latter are usually negligible for differential work and in small fields.
This leaves the transforms, which allow PhotoRed to correct for color differences
between the target and comparisons.
3.
Create a Batch Reference File for the field. This specifies the location and magnitudes of the target, check, and one or more comparison stars.
4.
Create a Batch Definition File. This is a list of objects to be measured with each
item in the list using a Batch Reference File and measuring a specific set of files.
5.
6.
165
Since the path name is long, and more so after adding one of two subdirectories
and a file name, the path to files will be shortened to, for example,
\..\20071023\
1.
2.
3.
166
4.
Select Utilities | Generate batch reference file from the main menu. This displays the Batch Reference File Generator (called the Generator from here on).
5.
6.
7.
8.
Click <Get> on the Generator. This puts the stars information into the entry
fields.
167
The magnitudes are from the 2MASS J-K conversion (see Warner, 2007, Minor
Planet Bul. 34, 113-119). However, the reduced magnitudes using the PhotoRed reduction methods found different values for the stars to be used. Those values will be
used instead.
Changing the catalog values will not be uncommon. If the AAVSO chart for the field has
adopted magnitudes for comparisons, you must use those instead. The main advantage of
using this approach is that the Generator automatically enters the correct RA and Declination for you, thus avoiding a critical error.
9.
Change the V magnitude to 13.352 and the CI (Color Index) magnitude to 0.481.
10.
There is no AUID (AAVSO Unique Identifier) for any of the stars in the field
(except the variable, of course), so leave the Name field alone. A comment entry
in the batch processing file lets the AAVSO know these are 2MASS identifiers.
Do NOT change the group number. It must be -1 for this exercise.
11.
Check the Comp radio button and then click <Add>. The information is placed
into the list.
When you click <Add>, all the entry fields are set to default values that the Generator recognizes as no star. To see the values in the fields for an item, click on that
item in the list.
12.
<Ctrl+Click> on star B. Change the data to the values below and add the star to
the list.
V = 12.319
CI = 0.381
Type = CHECK
13.
<Ctrl+Click> on star T (the target) and repeat the steps for adding the star to the
list, using these values in the Generator before you add the star.
V = 99.999
CI = 0.486
168
Naturally, if you are working a target where one or more of the other magnitudes are
known in other bands for the comps and/or check, and you are going to find reduced
magnitudes in those bands, you want to enter the proper values in the appropriate fields.
Many variables change color as they vary. You can chose to enter a compromise value,
which introduces a maximum error about one-half the total range of color index values,
or create different Batch Reference Files for maximum, mid-range, and minimum. This
would imply, of course, that you know where on its cycle the star was when you measure
a set of images.
14.
Enter W20111220 (without quotes) into the Chart (Sequence) ID field. This
identifies the source of the comp magnitudes. If you were importing data from the
AAVSO Variable Star Plotter (VSP), this field should be automatically field with
the ID assigned by the VSP. In this case, the ID is simply a last name initial and
the date the file was generated. You can use any scheme you want, but keep the
entry to 10 or fewer characters.
15.
16.
Click <Save> and save the file in the default directory of \MPO\UDATA. For this
lesson, use the name
ET_PER_TUTORIAL.TXT
The Batch Reference Files are simple text files that can be edited in a text editor. See the
PhotoRed manual for the formatting of the.
17.
This completes this part of the process. The next step is to create a Batch Definition File
that uses this information.
169
Never create a definition that measures images using different filter and/or taken on
different dates or of different targets. The strict rule is for a definition is: one target, one filter, one date per definition. The good news is that you can use the same
Batch Reference file for all definitions involving the same target, regardless of filter
or date, just as long as the comps, check, and target information do not change.
This tutorial is going to use transforms, i.e., apply color corrections to account for different color indices for the target, comp, and check stars. This is not a required step and
should be done only if you know youve established good transforms for your system.
The derived magnitudes are based on differential photometry, so the absolute first
order extinction (FOE) and nightly zero point values are not critical. However, it is
important to have FOE values such that the difference between kv and kr is accurate. Otherwise, the correction based on color index will not be right.
18.
If you want to save whatever transforms values you have, select Photometry |
Backup saved transforms from the main menu before proceeding. You can reload the original transforms when needed.
19.
20.
Match the screen shots below for V and R FOE and the V and R transforms (including hidden).
170
21.
22.
23.
Go to the Reductions page (<Ctrl+2>) and select AAVSO Batch for the reduction method.
24.
25.
Click <Add> on the form. This puts it into edit mode and places the cursor in the
Target Data | Name.
26.
Click the folder button next to the Ref File field. This displays a file open dialog. Locate and select the Batch Reference File that you created above. The form
should now look similar to the screen shot below. Of course, the path to the reference file may be different.
171
Note that the AAVSO and Chart fields have been filled using information
from the Batch Reference File. Be sure to confirm that this information is correct.
27.
The Name field is not the AAVSO name but a unique identifier for the batch definition. You could use anything that makes sense to you.
172
Click <Add> under the Images field (a list box). This displays a message.
29.
Click <Yes>. This displays the Windows directory selector. Locate and select
\..\20071023. Remember that \..\ is short-hand for the path to the example images.
30.
Saving Time
Selecting a folder instead of a specific set of images can save you a considerable amount
of time by not having to edit the batch definitions every time you measure images of the
same target. For example, say you regularly observe ET PER. You may save the images
for a given night, maybe even filter, in separate directories. However, if you move the
images to be measured even temporarily to a work directory, e.g., ETPER_WORK, you
can specify this temporary directory in the batch definition. Measure the images and then
move them to their permanent location. Unless something else changes about the definition, you need to create it only once.
Repeat the above steps above to add a new definition. Keep all information the
same except:
Name: ET PER 2
Images: \..\20071026 (click <Yes> to Select an entire folder?).
173
32.
Click <Save>. Save the definition in the default directory, \MPO\UDATA\, using
the name
ET_PER_TUTORIAL.BFD
Measuring Images
33.
34.
35.
36.
Select V-R (remember it was V-R magnitudes that were entered in the Batch
Reference File) and then click <OK>. This starts batch process.
37.
The batch processing takes some time since each image must be opened, the stars
extracted, and then PhotoRed does an AutoMatch. This is why is important to
have the correct configuration settings. PhotoRed extracts the stars in the first
pass for each definition and then does the AutoMatch. PhotoRed uses the plate
constants and the positions in the Batch Reference File to find the target, check,
and comparison and then measure their magnitudes.
174
Do not disturb the program during batch processing. You can try to work on other
things as long as they do not involve MPO files, e.g., you can do email. However, as
each image is loaded and AutoMatch invoked, PhotoRed grabs focus from any
other program. The best thing to do is sit back, keep an eye on things, and have a
cup of coffee.
Remember that youre processing two sets of images, so dont think things are done
after the first set of extractions and measurements.
Monitor the bottom status line of the Batch form. It tells you which image is being
processed and the image count versus total image count.
38.
After all the images have been processed, the AAVSO Report form should appear.
39.
If the first tab is blank. Click <Make List>. This generates a summary of all observations, showing the name of each star, the date, and number of observations
on that date. It does not subdivide by filter.
40.
Click <Save List> to save the text file to your hard drive.
The setting of Auto-generate is remembered each time you run the program. It is
off by default since generating the list can take some time if you have thousands of
observations. This is something that you can and should avoid. See the Data Import/Export chapter for lessons on managing the AAVSO tables. Later in this tutorial, youll export data for use in Canopus. The Archiving functions (located on the
Batch Edit/Delete/Archive tab) are covered in tutorials in the Data Import/Export
or Supplemental chapters.
Click on the AAVSO Reporting tab. You should see your data (and more if this
is not the first object youve measured).
42.
ET PER
2007/10/23
2007/10/27
Checked
175
XXX
| - pipe
The search for the Name field is case-sensitive, so entering ET Per would not find
any records in this case. If you do not enter a name, then all records within the
start-end range are included.
All data you want to include in the report must fall on or between the two dates at
00:00 UT.
43.
The default directory is the last one used when saving a file. If there is no previous
default, \MPO\UDATA is the default.
44.
Save the file. This is the file you would email or upload to the AAVSO if this
were live data.
As a check, the data were imported into a spread sheet where the mean magnitude and
standard deviation of the Check star were found: V = 12.312 0.012. Remember that the
Check star had a catalog value of 12.319, or a difference of 0.007 m. This provides another check on the quality of the data and accuracy of the results. Should the check star
actually be a different magnitude, the difference between its true and adopted values can
be applied to the target magnitudes to improve the accuracy of those results.
176
46.
In the Export to Canopus SESS/OBS section, set the entry fields as follows.
Name:
Start:
End:
V:
ET PER
2007/10/23
2007/10/23
Checked
00:00
18:00
The start and end date/times are in UT. The start and end dates can be different but the
total time between the combined start date/time and end/date time must be less than 24
hours.
47.
48.
Keep in mind the rules for a Canopus session: one target, same comparisons, same
filter, same date. The last is a little flexible since many locations will have data from
a single run that cross dates, e.g., those in the Eastern Hemisphere who start observing before 00:00 UT and keep observing until after 00:00 UT the next day.
Click <Generate> to export the data to a pair of export files (SESS/OBS).
PhotoRed generates a unique base name for the pair, so previous exports, even if of
the exact same data set, are not overwritten. Your file names may be different.
Repeat the steps above in this section, with the following exception
Start:
End:
49.
2007/10/26
2007/10/26
00:00
18:00
177
Open the Canopus configuration form and go to the Photometry tab and match
the settings on the sub-tabs as show below. Save the revised settings.
51.
Select Photometry | Load saved session(s) from the Canopus main menu and
load the first of the two export file sets you created in the previous section.
52.
Select Photometry | Import into saved session | from file from the Canopus
main menu and import the second set of export files from the previous section of
this tutorial.
53.
Youll see a warning message that you are working with a temporary copy. If you
want to save the combined set, select Utilities | Backup current sessions/observations and save the combined set. They will be saved with the same
name as the first files you loaded but you can save them to a different directory.
54.
55.
Why set the orders so high? As youll see, this is an RR Lyr star that has a very fast rise
from minimum to maximum. The higher number of orders allows the Fourier curve to
follow that sharp rise and rapid roll-over towards minimum more closely as well as a
small bump towards minimum. If this were a typical contact binary star, you could use
lower orders. An Algol-type variable would also probably require higher orders since
they, too, have relatively sharp rises and falls during the primary eclipse.
56.
178
Click <Find> and select both sessions. After a few seconds, you should see a plot
similar to this.
57.
You can try to refine the period search by setting the Min to something a little less
that 0.394 and decreasing the step size. The full data set for this star included
more than 700 data points and led to finding a very precise period.
In Conclusion
As youve seen, PhotoRed can process large numbers of images on any number of targets
from a given night. You can generate reports to submit to the AAVSO (ET PER was one
of those submitted), and do period analysis in Canopus.
179
180
11.
The Quick Binzel method is designed to use a small number of images in two standard
filters, V and R, to find the V magnitude of an asteroid or other object as well as the V-R
of that target. Its primary use is to obtain V magnitudes over a range of phase angles to
determine the H and G parameters of an asteroid.
This tutorial does not include detailed step-by-step instructions. See previous tutorials on
using the Photometry Wizard for Transforms for a general description of how to match an
image and select stars from a catalog using the matched image.
2.
3.
Select Quick Binzel on the first tab of the wizard and MPOSC3 from the drop
down list for the catalog.
4.
Select as many reference stars from the catalog as possible that are on the image.
Try to avoid excessively faint stars (low SNR) or those that are non-linear or saturated. Also try to avoid those with close companions such that a one pixel error in
the placement of the aperture would change the reading of the star.
5.
Measure the targets position on the initial image. In Quick Binzel mode, the
wizard includes a separate page where you set the target name and location.
6.
Measure the images. Put all data into Group 1, making sure that the correct filter
is selected on the image list before you click <Accept>.
7.
8.
Run the Quick Binzel reduction method. If you select the Clear filter, youll see a
warning message.
181
9.
Uncheck the Find offset only box so that PhotoRed can find the average difference between instrumental and catalog magnitudes.
11.
12.
13.
Click <OK>. Since you checked Enter CI, an entry form appears.
14.
Enter the color index of the objects. The default is 0.45, the average V-R for asteroids.
15.
This option allows you to reset the nightly zero point used for the transforms and
all-sky photometry. It uses the current value for first order extinction to calculate
the result.
16.
182
The plot above shows the individual c-V values for each star in each image. In this case,
there were 11 stars and 3 images. You can click on a given data point to see which star
produced the data point. The above data is acceptable, with a standard deviation of 0.034
mag. If using the LONEOS or User Star catalogs, you can usually get this down to 0.02
mag or less.
The Reduced Target tab shows the reduced V magnitude of the target.
Since there were three images, PhotoRed found the V magnitude of the target three times
by applying its corrected instrumental magnitude (for color and extinction) to the Average Field Offset value shown above. The three values are plotted as are the average of the
three values and the standard deviation.
If you had measured V and R images, you could have obtained the V-R value for the
target. You dont see this page because you chose to enter the value manually. Had those
images been available, the V-R would have been computed by making three pairs of v-r
instrumental magnitudes (corrected for extinction) and applying the hidden transform to
each v-r value. The average of these three and the standard deviation are shown in the
header of the plot along with the average JD of the set of images and the reduced V magnitude of the target.
Note the error for the reduced magnitude. The error given immediately after the magnitude is the standard deviation of the three measurements. The value in the parentheses is
183
184
It is extremely important that you not change the names of the files or, if you do, that
the base file name is exactly the same for both files and that you keep the
_SESS.FF2 and _OBS.FF2 endings..
Canopus relies on the SESS file name to find the corresponding OBS file. You cannot
specify them separately when loading data into Canopus. For example, in the examples
above, if the SESS file was named 771_LIBERA_1_SESS.FF2 and the OBS file was
771_LIBERA_2_OBS.FF2, Canopus would not be able to find the OBS file.
The file names are not case-sensitive. Otherwise, they must follow these rules.
Text Files
As youll see in the tutorials, you can export to and import from text files. Certain information must be available for at least the OBS file (the observations data). You will have
to create the SESS file or use an existing one. This is covered in the tutorials.
185
Version 10 cannot convert data from v8.x and before on-the-fly. Those data files
must be converted to v10 using the CONVERT10 program found in
\MPO\PROGRAMS. See the CONVERT10 documentation and Installation Guide for
additional information.
186
1.
It is very easy to create a Canopus export set from data in your PHSESS/PHOBS files.
You would create an export set for a number of reasons.
For archiving purposes so that you have the data for a given object in its own
files.
To share with other Canopus users who want to incorporate your data into their
data sets.
To create a Master export set that will be used to import data from other users
in a collaboration. This keeps the data for the given object isolated from your data
set and allows you to resend the merged set back to the other Canopus users.
The easiest way to exchange data among Canopus users is with export sets. Sending
text files requires additional work and can lead to errors.
1.
2.
Select the first session for the object for which you want to export data. In the example above, it is 771 Libera.
3.
4.
For this tutorial, we want on only Canopus export files (an export set), so check
the Canopus export files radio button.
5.
187
The JD Correction setting has no meaning for an export set. The raw dates are
exported so that the receiver can have Canopus apply or not apply corrections as he
wants. By never applying an corrections in an export set, there is never any confusion about whether or not the Julian Dates have been adjusted for light-time correction or to heliocentric JD.
Do not check the Delete after export box! If you do, you will get a warning
message that allows you back out before you permanently remove the data from
the PHSESS/PHOBSS files.
The Delete after export option should be used VERY carefully. The Delete
button on the sessions form accomplishes the same goal, but one session at a time.
7.
Click <OK>. This displays the (what should be) familiar Select Sessions form.
8.
Select the All radio button so that all sessions, regardless of filter setting in the
sessions, are included in the output.
If you select another filter, then only those sessions with that filter are included in
the list and will be in the export set.
9.
10.
After you select the directory, you should see a confirmation message that the
files were saved.
11.
Use Windows explorer to check the directory for the two files in the export set. If
you were to repeat the process above selecting sessions with the same Object
name, you would get a separate set of export files with the auto-increment number
one higher.
If sending to another Canopus user, be sure to send both files of an export set.
188
2.
Not every person with whom you work will be using Canopus. You can still provide
them with as much or as little of the information within the SESS and OBS files as you
want to send or they need to receive.
When exporting to a text file, the raw data are output, meaning that no corrections for phase angle or distance are included.
This tutorial will be very brief, since it follows the same steps as the previous lesson up to
the point where you get to the Output Options form.
Open the configuration form and set the type of magnitudes to be used when calculating the average of the comparisons and the differential magnitude of the target.
The Method is whats used for the calculations. The Plot Method setting as
no effect.
If you select Derived make sure that you have magnitudes entered in the CM
field for each comparison that is used. Otherwise, you will get exaggerated errors.
Do not select the deprecated Transformed or Transformed Absolute method.
2.
For this tutorial, select Instrumental, which is also the choice you should make
if catalog magnitudes are not entered in the CM fields for the comparisons on
the sessions form.
Open the sessions form and select the first session for the object for which you
want to export data.
4.
5.
189
7.
Select the sessions to included and click <OK>. This presents a file dialog so that
you can change the default name and path if you want.
8.
Save the file. If successful, you will see a confirmation message giving the full
file name.
2009/11/11 07:00:00
C
0.00
14.47
1.609
0.704
011546.19
011524.53
011546.39
011551.63
011541.68
+464250.6
+464234.2
+464430.3
+464420.4
+464106.1
14.230
13.716
14.334
14.187
14.302
Y
Y
Y
Y
Y
Note the header includes a line regarding JD correction. In this case, the JD values
were corrected for light-time, the time light took to go from the asteroid to the Earth
(in days). The JD values for this session were all corrected by subtracting 0.004063
days from the raw JD. Light-time correction is always subtracted from the raw JD.
This is followed by a detailed listing of the observation data, with the instrumental magnitude for the target and each comparison. At the far right of each line are three columns
that give data based on the Photometry Method setting.
190
CA
O-CAvg
Err
-------------------0.000 +14.396 0.068
The errors can be considerably larger with the DerivedMags, mostly because of inconsistencies within the catalog and color differences between each comp star and the target.
The Instrumental errors are based solely on the 1.0857/SNR conversion for each comparison and target.
At the end of the output file is a condensed listing that gives the JD, differential (or
DerivedMags) value, and error.
Condensed Observation Data for Session: 3511
--------------------------------------------U
J.D.
Obj-CAvg
Err
---------------------------------------Y 2455146.85369
+0.242
0.003
Y 2455146.85679
+0.228
0.003
191
192
3.
The previous tutorial covered exporting the session and observations data for an object to
a text file using a fixed format. While this usually works in most cases, you can also output the data in a custom format should a collaborator want the data in specific format.
The custom output generator provides a considerable amount of flexibility in this regard
but, of course, cannot possibly account for every contingency.
If you have not run the previous tutorial on exporting to a text file, at least read the introductory material.
1.
You do not have to select the Photometry Magnitudes method before running a
custom report. The Custom Output form allows you to make that selection.
2.
Open the sessions form (<Shift+Ctrl+S>) and highlight a session of the object for
which you are exporting data. Click <To File>. This displays the Output Options
form.
3.
4.
Check the Used Observations Only box. This setting can be overridden on the
Custom Output form.
5.
Set Every X Obs to 1 (meaning all used observations in this case). If this field
is, e.g., 2, then only every other observation that would otherwise be included
(based on the Used Observations Only box) is sent to the output file.
6.
Select the JD Correction | Light-time option for this tutorial. Your needs may
require a different selection.
7.
Click <OK> to display the Select Sessions form. Select one more sessions and
then click <OK> on that form to display the Custom Output form.
193
On the left is a list of stored data fields. These are values that are stored in the OBS file.
The lists and combo boxes to the right of the form are calculated fields that can be included in the report. The middle list is the fields, stored and calculated that have been
selected, the Select Fields list.
See the Reference Manual for details about the individual fields.
Using the instructions above, move the JD and Obj IMag fields from the
Fixed Fields list to the Selected Fields list.
9.
From the Calculated Fields list, select CAvg but do not move it just yet.
194
When using CAvg, you must also select which comp stars to use in building the average. This is done by selecting one or more items in the CAvg Stars list. If you
select a star that is not used, it is ignored.
10.
11.
12.
Move Obj-CAvg and Obj-CAvg Err to the Selected Fields list using a single move (select both items before moving).
13.
14.
15.
16.
Many people like to have a K-C magnitude as a check on the quality of the night. If
the value, which is the difference between two selected comp stars, is stable, then the
observations can be considered trustworthy.
Select Delta Mags | Instrumental. This means only the instrumental magnitudes
are used in computing the CAvg, Obj-CAvg, Obj-CAvg Err, and C1-K2 values.
If you select Derived, the CAvg value is 0.0. Obj-CAvg value is found by first
computing the instrumental difference between the target and the comparisons
catalog magnitude to get a derived magnitude for the target. This is repeated for
each comparison. The Obj-CAvg value is the mean of those derived magnitudes and
the error is the standard deviation of the mean.
If Transformed is selected, the C-K value uses the M/IR derived value for the comparisons.
17.
Set the remaining options to match those in the screen shot above.
18.
Click <Save> to save these settings when the Windows file dialog appears. The
file is a simple text file using the Windows INI file structure.
195
Click <OK> to run the report. When done, a Windows file dialog appears. Save
the file. The default name is the entry in the Object field of the sessions with
Custom included to indicate that it is not a standard output format text file.
Here is a short section of the file generated when running this tutorial on data for 1750
Eckert.
Observation Data Fields: JD,Obj IMag,CAvg,Obj-CAvg,Obj-CAvg Err,C1-K2
Observation Data for Session: 3511
---------------------------------2455146.853691
-7.823
-8.065
2455146.856792
-7.866
-8.094
2455146.859926
-7.806
-8.038
2455146.863077
-7.715
-7.961
+0.242
+0.228
+0.232
+0.246
+0.004
+0.004
+0.004
+0.005
+0.501
+0.502
+0.508
+0.494
The first line gives the names of the columns. If more than one session had been selected,
each session would be in its own block.
The non-data lines can easily be removed before the receiver imports the data. They are
put there to document the data, which should always be done to one degree or another.
A Custom Header
Sometimes a research requests a particular bit of information at the top of the file, something different from the standard Canopus header. If that is the case, you can create a text
file that has the header information thats needed and then specify Header | Custom on
the Custom Output form. You would also enter the full name of the file before generating
the report. You can use the speed button next to the Header Options radio button group
to locate the custom header file and enter its full name into the field.
196
4.
ALCDEF Export
At the 2010 October Division of Planetary Sciences (DPS) meeting, Stephens et al. proposed the adoption of the Asteroid Lightcurve Data Exchange Format (ALCDEF) as a
universal standard for reporting asteroid lightcurve data. Along the same lines as astrometry data submitted to the Minor Planet Center and the FITS standard for astronomical
images, the intent is to have a common data format so that asteroid lightcurve data can be
easily exchanged among observers and researchers.
With the standard in place, a person could write a single parsing program that would read
any ALCDEF file from any program and convert the data as the individual wants instead
of having to go through a number of steps to prepare the incoming data, as will be shown
in the importing tutorials to follow. Furthermore, the standard allows including any number of lightcurves for any number of asteroids in a single file instead of forcing, at the
least, one asteroid per file.
If you are in a collaboration where everyone is using MPO Canopus, it is recommended that you use Export Sets for both exporting and importing data since they
contain additional information not included among the mandatory data fields in the
ALCDEF.
This tutorial will be very similar to the previous one (Exporting to a Text File - II) except
that certain settings will be forced in order to conform to the ALCDEF standard.
It is presumed that you did the tutorials on measuring the lightcurve and finding the period for 771 Libera. If not, the principles are the same regardless of which sessions you
have available.
Open the sessions form and select one of the sessions for 771 Libera by clicking
on it in the sessions list.
197
Click <To File> button. This displays the Output Options form.
For maximum compatibility, always use uncorrected Julian Dates. Let the person
using the data determine if applying a correction is required. The same will apply to
magnitudes in the steps below.
3.
Set the form as show in the screen shot and then click <OK>. This displays the
Select Sessions form.
4.
Select all sessions for the asteroid and then click <OK>. This displays the Custom
Output form.
If you had multiple sessions for the asteroid that were not part of a single group,
e.g., data from several different apparitions, you could still select all sessions. You
do not necessarily have to generate a separate file for each apparition. This will be
covered a little more in steps below.
198
If you havent generated a custom report before, the form appears as above.
5.
Click the <Header | ALCDEF> radio button. Note the changes in the form. Do not
change any of the settings at this time.
The ALCDEF standard uses two mandatory and two optional items for the actual
lightcurve data. The required data are the JD and magnitude. The optional data are
the magnitude error and/or air mass. The JD, magnitude, and error are automatically selected by default. Do not add any other fields to this list. As shown below, whether or not air mass is included is handled a special way.
All comp stars have been selected. Canopus will include only those that are
marked Used in the sessions form for a given session.
199
Just as important, checking the Force DM TRUE box forces the DIFFERMAGS
keyword in the ALCDEF file to be TRUE and the STANDARD keyword to be NONE.
Setting these two keyword values correctly is vital to the proper interpretation and
use of the data.
When you select Derived, the normal way to compute the target magnitude is to
find the target value using
Instrumentaltarget Instrumentalcomp + Catalogcomp
for each comp star. The average of the derived values is the value for the target
and the standard deviation of the mean is the error.
This is fine if you did the original work in Canopus and used the Comp Star Selector to set the catalog values for each comparison. If, for example, youre working with imported data from another observer that included only JD/mag pairs,
there are no instrumental or catalog magnitudes to be used. Or, if working with
old Canopus data where you did not set the M/IR, then there will be M/IR derived
magnitudes for the target and comps but no instrumental magnitudes. In these and
similar cases, you must check Force DM TRUE. Canopus will then find differences as if all the data were strictly instrumental magnitudes and use those to find
the targets value.
If you select Derived and do not check Force DM TRUE, then DIFFERMAGS=FALSE and STANDARD=<NONE/INTERNAL/TRANSFORMED>
(see the description below for the settings on the ALCDEF tab).
If you select Transformed for the Delta Mags method, the usual procedure is
to use only the value in the ObjMag field. Unless these values truly are transformed to a standard system, you should check Force DM TRUE. If not
checked, the DIFFERMAGS and STANDARD keywords are set as with Derived.
Back to Business
The delimiter is forced to Pipe ( | ) since this is one of only two delimiters allowed in the ALCDEF standard. The comma and space are not allowed/
Remove LTC is not checked by default. If checked, Canopus will compute the
light-time correction for the mid-date/time defined in the session data. It will then
200
The main reason this option is provided is for when using data that were imported
from Canopus CONDENSED (or other) files where the JD were light-time corrected
using a constant value. If the JD were corrected point-by-point, you should not try to
work with the data at this point but find a way to remove the corrections by some
means and reimport the data into Canopus.
Point-by-point corrections (PxP corrections) are not allowed and so the checkbox is disabled. Allowing PxP corrections would require that a person wanting
uncorrected data would have to go through unreasonable efforts to remove the
corrections since each data point would have a unique correction value. For objects that were very near the Earth during a session, the correction value can
change dramatically over a short time.
Unity corrections by applying 5*log(rR) are allowed. If selected, a single value
is used so that a person wanting to remove the correction can do so. Before unity
corrections can be applied, you must load the orbital elements for the object. That
is covered in steps below.
6.
It is strongly recommended that you do not apply a unity correction. Let the enduser apply those, if desired.
Click on the ALCDEF tab.
Note how the longitude and latitude from the current profile have been automatically entered. You can change these if necessary.
7.
Do not check the Revised Data box at this time. You would do so if you were
generating a report to be resubmitted to a central database that modifies previously submitted data.
8.
Click the button to the right of the Object Name field. This displays a form that
loads the asteroids orbital elements.
201
9.
10.
Type 771 (without the quotes) into the Search for field. After a brief pause,
the list should reset to 771 Libera.
11.
Click <OK>. This loads the orbital elements into memory and sets the name in the
Object Name field.
12.
The ALCDEF standard requires a contact name and info (usually email). Enter
those in the appropriate field.
13.
Set the Output Band to the magnitude band used in the Comp Star Selector (or,
if you entered comp star magnitudes, the band for the magnitudes).
Note that this not the same as the filter setting in the sessions form. You may have
used a Clear filter but used R (or V) magnitudes in the Comp Star Selector.
14.
202
The magnitudes are on an arbitrary system. If using Instrumental for Delta Mags or you have check the Force
DM TRUE box, you must use this setting.
The magnitudes are on an internal system that approximates catalog values of the Output Band but have not
been formally transformed to that system. For example, you
would use this setting if you used the Comp Star Selector in
Canopus and the R magnitudes from the MPOSC3 catalog
(based on 2MASS to BVRI conversions).
The publication and BinCode fields are not required. However, if the data have
been used to generate a lightcurve that was published (formally, not just on a personal web page), enter the brief citation in this field.
16.
ALCDEF also requires the observers. These are not necessarily the same as the
contact person, e.g., if you are submitting data on behalf a collaboration. Use the
format indicated on the form.
Proper attribution of those involved in obtaining the data goes beyond common
courtesy and just being fair. Give credit when its due; youd want nothing less for
yourself.
17.
The comments field entry must apply to all sessions or, at the least, you must be
specific as to which session (date, observer) the comments apply.
18.
Click on the Core tab and then click <OK>. The program generates the report
and then presents a file save dialog. The file is given a default name of
ALCDEF_<SessionObjectName>.TXT
in the \MPO\UDATA directory. You can change the directory and/or name but
the file extension is forced to TXT.
19.
Click <Save> once you have set the file path and name. This saves the file and a
confirmation message appears.
203
If all the data in the PHSESS/PHOBS files are strictly your work product, then
proceed as above and list yourself as the only observer.
2.
If you worked with collaborators and their sessions are part of your
PHSESS/PHOBS files, then select all sessions and make sure to enter each observers name in the OBSERVERS entry field before you generate the file.
3.
The alternative to #2 is to report only your data and make sure the others do the
same. However, this could mean some lightcurves are not submitted and it breaks
the paper train of the combined set being a collaboration.
204
If you move or copy your data from PHSESS/PHOBS to export sets, you will
need to load those and generate the reports accordingly along the same lines as
above, i.e., fill in the OBSERVERS field as needed.
Where things may be left in a bit of confusion is which observer contributed a specific
lightcurve within a set. If using Canopus 10.2.0.5 or later, the NOTES field of a new
session includes a section that carries the core information from the Configuration | MPC
tab. When generating an ALCDEF report, Canopus looks for this block and, if present,
includes the contents of the MPC CONTACT1 line as a COMMENT in the ALCDEF
file. For example:
COMMENT=Canopus MPC Contact1: CONTACT1=B. D. Warner, 17995 Bakers Farm
Rd., Colorado Springs, CO, USA [[email protected]]
If your data were created before 10.2.0.5, this information is not included. You can either
manually edit (carefully) the ALCDEF report or, more likely the best course of action,
make no special effort unless it is to try to use the COMMENT entry field on the Custom
Output forms | ALCDEF tab. However, that could be prone to error. In this case, it may
be best to provide a little less information rather than spend extra time and take the
chance of providing wrong information.
205
206
5.
Importing export sets is relatively easy and straightforward. Where the trick likes is planning a routine when working in collaboration with other Canopus users. For example, if
you are the primary, then all data comes to you. As mentioned earlier, you really dont
want to merge other peoples data into your PHSESS/PHOBS files. Here is just one possible outline of how you would proceed.
Create an export set that includes all your sessions and observations for the target.
Rename the files to indicate they are the master files. For example, say your export files had the original names
771_LIBERA_1_SESS.FF2
771_LIBERA_1_OBS.FF2
Rename them to
771_LIBERA_MASTER_SESS.FF2
771_LIBERA_MASTER_OBS.FF2
Its also a good idea to create a new subdirectory under \MPO\UDATA to hold all
files that youll be getting, e.g., \MPO\UDATA\771_DATA.
The total length for each of the two export set files, excluding path, cannot exceed
32 characters. This is a limit of the data base engine. If you try to load files with file
names that are too long, the engine, via Canopus, reports an error and the data are
not loaded.
If you make additional observations, continue putting them into the
PHSESS/PHOBS files so that you have a complete record of your work.
Load the master set and import the new data from your PHSESS/PHOBS files and
collaborators, be it their export data sets or text files.
The following tutorials show you how to do handle all these basic steps.
Avoid Duplication
If you or your collaborators will be sending new data from time-to-time, it is very helpful
if there is no duplication of data. For example, say you get an export set from your collaborator that has data from Nov. 1 and you have incorporated that data into your master
set. He then observes on Nov. 3. Make sure that when he creates the export set to send
you that he selects only the Nov 3 session, not both. Otherwise, Canopus loads both sessions as new data, meaning that youll have two sessions from Nov 1 from your collaborator.
The same applies if youre creating export sets to send to someone else. Select only new
data each time you create a new export set. Remember that Canopus automatically creates new export sets with unique base names by using an auto-increment number.
You want to be careful if collecting sets from others. They may send you files that have
the same file name as some on your computer. I put all incoming files into a temporary
directory, immediately rename the two files in the set, often including the initials of the
observer in the new name, and then move the renamed files to the working directory.
207
Temporary Data
In many cases when working with export sets, you are not working on the original files
but a temporary copy. If you do not backup the files (covered in the tutorials), then you
will lose all changes when you close Canopus.
Assuming that you are not already working with an export set, select
Photometry | Load export set from the Canopus main menu. This displays a
Windows file dialog that is filtered to show only those files that end with
_SESS.FF2.
If you did the tutorials on creating an asteroid lightcurve and then exporting the
export set for those sessions, there should be at least one _SESS.FF2 file in the
dialog somewhere on your hard drive.
2.
If there is at least one session with data in the file, Canopus automatically sets the first
session as the default session. If there are no sessions, an error message appears and no
sessions are available.
3.
All the sessions are (or should be) for the same object. If there sessions with
slightly different names, this suggests that you received a file that included data
from several observers that had been created by steps that are outlined in the tuto-
208
You cannot be editing an existing record or adding a new one to use these features.
Press <Ctrl+F>. This displays an entry form where you can enter all or part of the
name that appears in the Object. The search is not case-sensitive, so Fred and
FRED will be found. The search looks for the search string anywhere within a
name, so 1503 will find (21503) Beksha or 1503 Kuopio
If no match is found, the table is positioned on the last session.
You cannot edit the table entries directly. You must use the entry fields below the
table.
Select Photometry | Clear session from the Canopus main menu. This clears the
export set and makes PHSESS/PHOBS the default photometry files. It also resets
209
The default session is the one to which new observations are appended or existing
data are edited. Its Object field serves as the filter when selecting sessions for doing period analysis or exporting data.
210
6.
As outlined earlier in this chapter, if you are the primary in a collaboration of several
observers, it is best to work with a master export set and not import data from several
observers into your primary PHSESS/PHOBS files. This tutorial presumes that you have
already loaded that master export set and now want to add additional data.
The first case is where you are in collaboration with yourself. Thats not as strange as it
sounds. There are some Canopus users who operate out of more than one observatory and
work with a copy of Canopus at each location. In this case, they have one computer
where they tend to do most of the work, the master, and create export sets on the computers to copy to the master where they can create a merged data set. This is covered in
the first part of this tutorial.
The second, and more likely case, is where you are brining in data from other Canopus
observers who have sent you one or more export sets. This is covered in the second part
of this tutorial
Load an export set, e.g., the master being used in a collaboration where the data
from several observers are being merged. For this lesson, an export set from the
771 Libera lightcurve tutorials (see the Photometry chapter) was loaded.
2.
Select Photometry | Import into export set | from PHSESS from the Canopus
main menu. This displays the Select Sessions form.
The form lists only those sessions from the PHSESS file that have the same name
as that in the Object field of the currently active session from the master export
set.
3.
Select one or more sessions and click <OK>. You will see some information messages as sessions and observations are imported and renumbered. Canopus assures
that there are no sessions with duplicated session numbers since those uniquely
identify a given session.
211
5.
6.
Open the sessions form (<Ctrl+S>) and confirm that the correct data were imported and that things look OK.
7.
If so, close the sessions form and immediately select Utilities | Backup current
sessions/observations from the Canopus main menu. This displays a select directory dialog.
If you want, you can save the updated export set to the original directory. If you
select that directory, a warning message appears that the files already exists. Click
<Yes> to overwrite the existing files.
8.
After the files are saved, a confirmation message is shown that gives the path to
which they were saved.
You can now work with the master set just as would the PHSESS/PHOBS files.
Select Utilities | Backup current sessions/observations files before you clear the
session or exit Canopus. Otherwise your changes will be lost.
2.
Select Photometry | Import into export set | from file from the Canopus main
menu. This displays a Windows file dialog that filters to show only files ending
with _SESS.FF2.
You can select one and only one export set per merge.
3.
212
Select the SESS file to merge and follow the steps outlined above, including
backup the merged data set.
If you have additional export sets to merge, repeat the steps above.
Be sure to backup the export set after you have imported all additional data.
There is no rule against merging sets that cover different objects. You would work
with the merged, multi-object set just as you do the PHSESS/PHOBS files.
213
214
7.
Many times you will need to import lightcurve data that did not come from Canopus.
Since that data can be in many forms, this requires a specialized import process. However, its not hard and can be easily learned.
There must be at least a date (JD or calendar date/time) and magnitude for each observation.
The import file must be simple ASCII text and have one and only one observation
per line.
If there are header lines (non-data), those must be removed before trying to import
the data.
The magnitudes can be instrumental or on a standard system.
Data must be imported into an existing, empty session.
If there are data for multiple nights, you must create a separate session for each
night. If necessary, parse the data file into several files, each one having the data
from one and only one night.
1.
Two sample files are available for this tutorial. They are found in
\MPO\EXAMPLES\DATAIMPORT. Only the file names (without this path) will be
given in the tutorial. These were for the asteroid 4265 Kani and kindly supplied by
Richard Miles of the British Astronomical Association
Open the sessions form (<Shift+Ctrl+S>) and click <New> to create a new session.
The only critical information thats needed is the Object, Mid-Time, Band,
and distance/phase data.
2.
Match the settings for these critical items as shown above. To get the distance/phase data, you need to click <Calc M/D/P> and locate 4265 in the look up
list.
3.
Click <Save> to save the session data but do not close the sessions form.
215
Because of internal structures, you cannot import the data while creating a new
session. So, you must save the session and then go immediately into edit mode.
4.
5.
Click <Import>. This displays the Observations import form (the Import form).
The Import form allows you to import data from two sources: MPC-compliant reports such as you create when you generate an MPC report from astrometry files
(see the Astrometry chapter) or Custom, meaning you specify the format of the
incoming data.
6.
The left side of the form is where you define the incoming data.
7.
Select the Delimiters | Tab button since the incoming data use tabs between the
values in a given line.
It is important that one and only one tab be used between fields. See the Reference
Manual for additional details on delimiters in the data file.
8.
216
Check JD in the fields list since that is one of the fields in the incoming data.
Locate Obj IMag and <Ctrl+Click> on it. If you click, you will lose the JD
selection since this is a typical Windows multi-select list box.
10.
Use the up arrow below the list to move Obj IMag just below JD.
11.
Located Obj Error in the list, <Ctrl+Click> on it, and move it to just below Obj
IMag.
You dont have to move the items in the list. This is just so that you see all the
checked items at the top and dont forget to select a field that you want.
12.
Set the Base JD to 2450000.0. The data in the file are Julian Date 2450000.0,
so you need to add that constant to the data so that the correct JD are stored in the
file.
13.
Do not check either of the two check boxes. The data are such that the object is
brightest at minimum value and the errors are in magnitudes.
If using DerivedMags you would still select the Obj IMag field. This is because
that method reads the Instrumental magnitude field in the observations file to do
the calculations. Only the Transformed or Transformed Absolute photometry
methods would use the Obj Mag field.
If we had selected Obj Mag, the data would go into the true magnitude field in the
observations file. Then, if you try to do analysis in the Instrumental or Derived mode,
there would be no data in the instrumental magnitude field and you get a blank plot along
with a complaint about no data.
14.
Click <Load>. This opens a Windows file dialog. Locate and open
A4265_20081025.TXT
217
If you dont see any data in the memo form, look at the status line at the very bottom
of the form. It will likely say something about an incorrect number of fields. This is
usually caused by setting the wrong delimiter or indicating more fields than there
actually are.
If you get an error about a bad date (probably a negative one), you forgot or improperly
set the Base JD value.
15.
You should see something like the screen shot above. Note the O Mag and O
IM columns. The O Mag column stores the true magnitude value. Here it
99.99 since we did not import data into that field. The O IM column holds the
data even though, in this case, the data look to be DerivedMags values. Again, the
DerivedMags photometry method reads the instrumental magnitude values from
the observations file to generate a final catalog magnitude for the target.
This goes back to section above about selecting a magnitude for the input data.
Canopus does not combine data from these two columns during analysis; its one
or the other. When merging data sets that include instrumental differential values
and true magnitudes, this means working with the DeltaComp values in the sessions form to get the sessions to merge properly.
16.
Click <OK>.
17.
You will see a message asking you to confirm that you want to import the data.
Click <Yes>.
18.
After a second or two, you should see a confirmation message saying that 104 of
104 records were imported. Close that message (click <OK>). This returns you to
the sessions form.
19.
218
Here is a plot of the data after doing a period search and before tweaking the nightly
zero points.
If you have a large number of data sets to import, this process can get tedious. However,
it is the only way to get it done within Canopus. A special utility program, the Lightcurve
Analysis Engine, shipped with MPO Canopus, allows importing large number of data sets
into the PHSESS/PHOBS files or, preferred, new export sets.
A tutorial on the LAE is forthcoming. In the meantime, consult the on-line help file for
the program.
219
220
8.
ALCDEF Import
At the 2010 October Division of Planetary Sciences (DPS) meeting, Stephens et al. proposed the adoption of the Asteroid Lightcurve Data Exchange Format (ALCDEF) as a
universal standard for reporting asteroid lightcurve data. Along the same lines as astrometry data submitted to the Minor Planet Center and the FITS standard for astronomical
images, the intent is to have a common data format so that asteroid lightcurve data can be
easily exchanged among observers and researchers.
With the standard in place, a person could write a single parsing program that would read
any ALCDEF file from any program and convert the data as the individual wants instead
of having to go through a number of steps to prepare the incoming data, as will be shown
in the importing tutorials to follow. Furthermore, the standard allows including any number of lightcurves for any number of asteroids in a single file instead of forcing, at the
least, one asteroid per file.
If you are in a collaboration where everyone is using MPO Canopus, it is recommended that you use Export Sets for both exporting and importing data since they
contain additional information not included among the mandatory data fields in the
ALCDEF.
Just to be Safe
A detailed description of the ALCDEF standard will not be given here. Instead visit
http://www.MinorPlanetObserver.com/ALCDEF/index.htm
In brief, however, each session and its data is called a lightcurve block that consists of
two subblocks: metadata, which is essentially the session data, and data, which is the
actual lightcurve observations. A single file can have any number of lightcurve blocks
involving any number of objects and sessions.
Canopus allows importing information from one or more ALCDEF files at a time. Each
file can include data for only one asteroid or for multiple asteroids. Since things can go
wrong during an import, Canopus does not load the file directly into the
PHSESS/PHOBS files. Furthermore, the import feature is disabled if you are working
with an export set. Instead, Canopus creates a new, empty export set and then imports the
data into those two files.
As it does so, Canopus keeps track of any errors, including the specific file and lightcurve
block within the file. At the end of the import process, Canopus displays a Success
message or a memo form giving the details of the errors.
Even if there are no errors, you should then load the newly created export set and confirm
that the information is as you expect and want. From there, you can manage the individual sessions as described in the previous tutorials and the manuals.
Getting Started
The following presumes that you did the ALCDEF Export tutorial on page 197and so
have an ALCDEF-compatible file handy.
1.
Select Photometry | Import ALCDED to empty export set form the main menu.
This displays a file open dialog that allows you select more than one text file.
221
2.
In this case, select one and only one file, the one created in the ALCDEF export
tutorial and then click <Open>.
3.
The program imports the data into the empty export set. If there are no errors, a
file save dialog appears with the default directory set to \MPO\UDATA and the
file name to
ALCDEF_1_SESS.FF2
It is very important that you do not change the suffix and extension, i.e.,
_SESS.FF2. If you recall, all export sets have two files, one ending in _SESS.FF2
and the other _OBS.FF2. If you change the base name (whats before the suffix), the
program automatically changes the base name for the _OBS.FF2file.
4.
After youve set the file path and base name, click <Save>. If the file already exists, youre asked if you want to replace it.
5.
6.
Click <OK>
7.
Select Photometry | Load export set from the main menu and check that the data
were properly imported. This includes checking the settings for each session (filter, date, etc.
222
223
224
Supplemental Tutorials
Supplemental Tutorials
The tutorials that follow provide insights into using specialized features of MPO Canopus
and PhotoRed. Among these are the dual-period search in Canopus when working with
binary asteroids and the AAVSO-driven methods in PhotoRed.
The supplemental tutorials are designed for advanced users in the sense that they do not
walk through procedures literally step-by-step, down to individual keystrokes, mouse
clicks, and menu selection. Its assumed the user is familiar with the Windows operating
system and doesnt need detailed instructions when told to open a given image, to select a
range of items in a list box, or save a set of data. Furthermore, its assumed the user has a
solid familiarity with the MPO programs, enough so that they are able to find menu items
on their own and have a good grasp of the fundamentals of how to open and measure
images using the various wizards. The tutorials in the main section of this Users Guide
provide somewhat more detailed steps to acquire this familiarity and should be run not
only on the sample images and data but on the users data as well. Throughout, however,
it is assumed that the user is comfortable with the Windows operating system and its
common controls, e.g., file dialogs, multi-select lists, and so on.
By not going into the greater detail of the main tutorials (which takes considerable time
to document), its hoped that this Users Guide will expand over time to include many
more supplemental tutorials that allow you a more complete understanding of MPO
Canopus and PhotoRed. This will let you get the most benefit and enjoyment from what I
consider to be an extremely powerful toolbox for research on asteroids, variable stars,
and double stars, among others.
The tutorials are offered in no particular order of importance. You get to decide which
ones matter most to you.
Brian D. Warner
Palmer Divide Observatory
2011 December
225
Supplemental Tutorials
226
1.
MPO Canopus incorporates Alan Harris FALC program, a Fourier analysis algorithm
that has become the industry standard for asteroid lightcurve period analysis when only
a single period is involved. With the discovery of binary asteroids, the need for finding
two or more periods within an asteroid lightcurve was required. That capability has now
been added, but still using the FALC routines as the underlying engine.
For more complete details about binary period analysis, refer to the paper by Pravec et al.
(2006) in Icarus. A brief outline follows here. It is not meant to be even slightly comprehensive or thorough in all details.
In a typical asynchronous binary system, a small satellite rotates about the primary asteroid, its rotation period being tidally locked to the orbital period, like the Earths moon,
but this is not always the case. The primary spins at a faster rate. For NEAs, the average
primary rotation rate is between 2 and 5 hours with the orbital periods of the satellite
going from 12 to 40 hours. In some cases, careful analysis of the lightcurve shows that
the satellite contributes a small amount of variation on top of that of the primary as it
rotates when the satellite is not tidally locked to the period of its orbit or if it is locked but
is sufficiently elongated so that the out of eclipse time shows a upward hump. Those
who work certain types of eclipsing variable stars are familiar with this behavior.
If the viewing aspect is just right, occultations and eclipse events can be seen in the lightcurve as unexpected drops in the overall brightness. If there are events, then youll be
finding the primary rotation and, most likely, orbital periods since those will dominate
the data. You can then try to find the rotation period of the secondary, but keep in mind
that it may not be rotating independently of its orbital period and/or the amplitude of the
lightcurve is too small in comparison to the average scatter of the data.
If there are no events, then youll find the primary rotation period and can then try to find
a secondary rotation period due to the rotation of the satellite. Again, it may be too small
in comparison to the scatter of the data or be tidally locked.
Complicating matters even more is if the phase angle is large enough, conditions may be
right such that the shadow of the satellite on the primary is visible independently, in
whole or in part, of the satellite itself. Observers of Jupiters four largest satellites are
familiar with the difference between a transit of the moon and the transit of its shadow.
In MPO Canopus, the search for the various periods and amplitudes within the total lightcurve is an iterative process that proceeds as follows:
Using the standard procedure, an initial primary period is found and that curve
is subtracted from the data, leaving only the variations due to the satellite and
events. Unless the data are of high quality, the low amplitude lightcurve of the
satellite is lost in the noise of the data. However, the event data should dominate
whats left after subtraction and so allow finding of the orbital period.
The secondary period is subtracted from the data, allowing a refinement search for
the period of the primary. Once that is found, the Fourier coefficients are again
saved and used to subtract the primary curve from the data.
The secondary period is again analyzed, with the final results being saved.
Once the primary and orbital lightcurve Fourier coefficients have been found, a
search can be made for the satellites curve. The depth of the events can be used
to derive an estimate of the size ratio of the two objects, given as D2/D1 and in the
range of 0 < x < 1.
227
Again, as you step through this lesson, your results may not exactly match those
presented. They should be very similar if you follow the instructions carefully. Keep
your eye on the big picture, which is to learn the concepts behind and the steps
used in a dual period search.
Getting Started
1.
On the Configuration | Photometry tab, check Instrumental for the photometry magnitudes method and Range for the plot method.
2.
As noted above, make sure the UTOffset is 00:00:00. Save the changes if necessary.
If the UTOffset for the profile is not normally 00:00:00, be sure to reset this value to
its original state after youre done with this tutorial.
3.
Load the SESS/OBS files for the asteroid using Photometry | Load export set on
the Canopus main menu
4.
On the Photometry page (where you do period searches) set the initial search parameters to the following:
5.
The reason this period range was select is that visual inspection of the raw plots
showed that after accounting for obvious events, the period was very likely
around 2-4 hours. A search going further out does find a very strong solution at
6.05 hours. However, the curve is dominated by the presumed events and forms
228
Use a high number of orders when doing a dual period search. This helps by allowing the Fourier solution to cling tighter to the actual data.
You can check the Fourier box on the period search page when finding the primary curve and there is sufficient data to fill in the curve. You should not check the
box when doing the initial searches for the orbital period since missing data and
high orders often result in very wild Fourier curves that greatly expand the vertical range of the plot. Once you start narrowing in on things, you can elect to include
the Fourier curve for one or both curves as appropriate.
7.
Click the down arrow button next to the Find button. This displays a submenu:
8.
Select Dual period search. The Find period is bold because that is the menu
item selected if you just click <Find>. You should see a confirmation message.
9.
Should you have the dual period search form open and click <Find> or select the
Find period menu button, youll see a warning message that asks if you want to
close the Dual Search Parameters form and do a single period search.
Select the Clear Info option.
See the Reference Manual for an explanation of the options.
10.
11.
12.
After a few moments, Canopus computes the period and presents the plot, which
should look something like the one below.
229
Before going on, take a look at the period spectrum plot on the Period Search Results form.
There are three or four possible solutions, obviously none of them very outstanding since they show as only slight dips below the rest of the data.
Which one, if any, is correct? In this case, go with the period suggested by
Canopus but keep in mind for the future that the data due to events influences the
period solution, sometimes dramatically. In the lightcurve plot, you can clearly
see the data points that might be due to eclipse/occultation events.
From here on, the Period Search Results form will be called just Results form.
14.
230
In addition to the Results form you see when finding a period, the Dual Period Parameters form now appears.
15.
This form remains open unless you deliberately close it. Dont do that right now.
Note that it is empty. That will also soon change.
16.
17.
Click <Save>. The default directory is \MPO\UDATA but you can save the file
anywhere. For this tutorial, use the name
A5477_P1.TXT
This convention tells you the name of the asteroid (or number) and that the data
are for Period 1.
18.
After you click <Save>, another form appears since you have the Dual Period Parameters form open (DPP for short).
231
19.
This form allows you to load the Fourier values directly into the DPP, assigning
them to either Period 1 or Period 2.
20.
Check Import to Period 1, click <OK>, and then look at the DPP.
Note that the Period 1 memo field now has the Fourier values.
As with almost all tutorials in the Users Guide, dont worry if your results are not
exactly the same. As long as they are in the ballpark, keep moving ahead. The
goal is to learn the process, not duplicate results. Unlike a chemistry professor I
had, you dont get a reduced grade for not getting exactly the same results he did
even when the grad student confirmed my observations.
21.
Check the Subtract | Period 1 radio button and then click <Data | Save>. Save
the DPP data as
A5477_ALL.TXT
22.
Click <OK> on the Results form to get it out of the way. Move the DPP to some
convenient location.
You can also double-click on the DPP caption bar to roll it up reduce it to just
the caption bar and maybe a small bit of the form under it. Double-click the caption
bar again to restore the form to its normal size.
Do not click <Close> on the DPP at this time.
232
24.
Confirm that the Subtract | Period 1 button is checked. This subtracts the theoretical curve defined by the Period 1 Fourier coefficients from the data before the
period search is started.
25.
Click the down arrow next to <Find> and select Dual period search from the
drop down menu. This displays the Dual Period Info options form.
26.
Select Use current info. This uses the information in the DPP form and subtracts Period 1 from the data set.
27.
Select all sessions in the Select Sessions form and click <OK>. You should see
results similar to these.
The events are better defined, though there is still some oscillation outside the
events. This is a possible clue that you have not locked onto the true primary period.
28.
Since were in the neighborhood, click <OK> on the Results form to clear it.
233
Set the period search parameters to the following and re-run the dual period
search.
30.
Youll see a plot similar to the one above, but the period may be slightly different.
When running this exercise, the new period was 24.39 0.02 hr.
Dont try to go with higher precision for now; we need to be sure of the primary
period.
31.
Go to the Data tab on the Results form and save the Fourier values using the name
A5477_P2.TXT
This displays the Dual Period Import form (Import form from here on).
32.
234
Select Import to Period 2 and click <OK>. Note how the values are now in the
Period 2 memo field on the DPP.
33.
Check the Subtract | Period 2 box and uncheck the Subtract | Period 1 box.
34.
35.
37.
Set the period search parameters as follows and run a dual period search, selecting
Use current info and all sessions.
Why these parameters? First they include the initial primary period of about 3.05
hours. However, they also include a period that was almost as good that was at
around 2.99 hours. Remember that the initial period was influenced by the then unsubtracted event data. When this is the case, the true period is often a little less
that what is initially found, especially if the event data are strong.
The period spectrum shows that the 2.99 hour solution is now favored but by the
slightest amount.
235
When you encounter this situation, i.e., two periods very close to one another, its
best to include both periods in future searches to see which one eventually prevails.
Again, because the secondary data period may not be perfect at this time, it doesnt
hurt to presume the shorter period due to the influence of the event data on the solution
(makes it longer than it might be). Here is the plot of the data at this point.
38.
Once again, save the Fourier values from the Data tab on the Results form and
import them into the DPP, making sure to select Period 1. Save the data under
the previous name, i.e., A5477_P1.TXT
39.
Check Subtract | Period 1, uncheck Subtract | Period 2, and save the DPP data
under the previous name, i.e., A5477_ALL.TXT.
40.
Set the search parameters as follows and do another dual period search.
41.
236
42.
Save the Fourier values in the Results form, import them into Period 2 on the
DPP, check Period 2, uncheck Period 1, and Save the DPP data.
You should be starting to get into the rhythm of switching between which period to
search and saving the results after each search.
43.
This will be the last pass at finding the primary. The secondary period does not
need to be examined again, unless additional data are available.
44.
Set the period search parameters as follows and do the dual period search.
45.
Make sure that you save the Fourier data from the Results form, import it into the
DPP, and then save the DPP data set.
237
There is no longer any ambiguity! The 3.05 hour solution is now a relatively minor blip.
In the sense of full disclosure, the search can be very sensitive to very minor changes in
either period. If too precise of a solution for the second period had been tried, it may have
been entirely possible that the longer period near 3 hours prevailed.
To repeat, when you have two solutions that are very close to one another and one
or the other appears to be the truth. you should keep working on refinements until
you can eliminate one or the other with as much certainty as possible.
46.
Again, there is no real need to work on the second period again. However, you
can see if you can get a slightly better solution on the primary period.
47.
Leave DPP settings as they were, i.e., subtracting Period 2, and do a dual period
search using these search parameters:
Note that the plot range (at the far right of the parameters tool bar, has been reset to
0.40. This expands the plots vertical scale so that you can see the curve in more detail. In general, you want to keep the value at 0.8 so that 1 magnitude in amplitude is
about the same scale as 1 cycle of the lightcurve. This avoids giving a false impression of amplitude changes in the curve, which can be interpreted incorrectly as
having significant physical meaning.
238
The precision is as low as you can go. Note that there are still some data points from
sessions that might have had events. Manual elimination of these points, at least temporarily, might improve things but doing so would take more effort that resulting gain and
would be based at least on some guesswork.
Start the search fresh using the initial search parameters, selecting Clear info,
and selecting all sessions.
49.
When the Dual Search Parameters form appears, check the Show Plot box towards the lower left.
239
If necessary, move the form so that you can see all of it, including the plot.
51.
The plot shows all the data as black dots, except for those associated with the session selected in the drop down list at upper right. As a quick aside, change the selection in the list and see how the plot changes.
52.
In the screen shot above, the event data is very apparent in the form of data well
below the general curve and that some of the data points from session 1 are
probably due to an event. We want to eliminate these from the solution.
53.
Click on the left-most green data point that is in the event area. Then click
<Set Start>. The status bar at the bottom shows the date/time of that data point.
54.
Click the right-most green data point that is in the event area. This is not always
obvious but make the best estimate you can.
55.
Click <Set End>. If the end point you selected is later in time than the start point,
the end date and time are shown in the status bar and the two data points and all
those that were in between are changed to red. The screen shot below shows the
results of excluding points from session 1.
You may get a warning message that the end date is before the start date. This is a
phased plot, so data points are not necessarily in chronological order. Keep working to find the two end points that exclude the most number of event data points.
56.
Repeat the process for all sessions. Note that not all sessions have event data. Do
not exclude data from those.
57.
Once you have reviewed all sessions, click <Save> to save the entire set.
Use a different file name for the ALL set, e.g., A5477_ALL_WEXCLUDE.TXT so
that you wont lose the results of whats been done until now.
58.
Uncheck the Show Plot box to shrink the form and then check only the Excluded box. This ignores the excluded data points when doing the period search
(this is not the same as subtracting them).
59.
Run the dual period search again without changing the search parameters, making
sure to select Use current info and all sessions. Heres the resulting plot.
240
60.
If you did a good job picking the event data points, there will be none or only a
few below the general plot of the rest of the data. Note that the period found is
2.99 hr. By eliminating the event data, the period search was able to find the approximate correct period on the first try. Heres the period spectrum.
These are the same plots. The one on the right is zoomed a little to show the differences better.
61.
The alternate 3.05 hr period value is significantly higher though not completely
lost in the noise.
62.
To find the orbital period, you want to uncheck the Excluded box so that the
event data are available for analysis and then check the Period 1 box to subtract
that information from the data set.
241
In this case, the period was close to the orbital period. However, theres nothing significant in this solution since the data are fairly noisy. In the analysis by Petr Pravec of a
more complete data set, P2 was not found, indicating either the satellite was tidally
locked to the primary or that the amplitude of the curve was too small to be reliably detected.
242
= [1 10(0.4 * 0.16)]1/2
= [1 0.8630]1/2
= 0.37 0.02
A Second Try
As mentioned at the start of this tutorial, there is a second set of files available in
\MPO\EXAMPLES\DUALPERIODSEARCH
These are for 6265 1985 TW2, another binary asteroid discovered by David Higgins,
Hunters Hill Observatory, Australia, who provided his original data.
This set is cleaner than that for 5477 and should provide an excellent second chance at
finding the primary and orbital periods. The steps wont be covered here, nor will the
final results so that you dont have any preconceptions. Besides, you can find the values
on the Internet. It wont hurt to give you a clue: the primary period is about 1/6 the orbital
period and the orbital period is less than a day.
A Final Thought
All of the discussions have presumed an asynchronous system where the satellite is
tidally locked to the orbit or is maybe spinning at a rate independent of the orbital period.
Also, the primary in the system has a period in the range of 2-5 hours and the orbital
period is 15-40 hours. These are not absolute limits. A discovery at the Palmer Divide
Observatory in late 2009 of revealed what is probably a binary asteroid of unusual type.
Here the primary has a period of almost 600 hours and the secondary has a period of
about 2.5 hours. The amplitude of the primary is ~0.40 mag while that of the secondary is
~0.04 mag. The orbital period is unknown since there were no mutual events.
In short, it never hurts to look for two additive periods in lightcurve data if there are signs
of such. In the end, the second period may be due to nothing more than the analysis latching onto noise. However, it just may be evidence of something real and lead to a new
discovery.
243
244
2.
The Moving Object Search allows you to search for moving objects in your images. For
those hunting new asteroids, its an ideal tool and beats using the blinker and having to
stare at the images, especially if you have dozens of image to review.
The basic principle is to take a set of three images, align and stack them, and the remove
any object that is the same place on all three images. Anything that remains is a potential
candidate. However, what remains may not be a real object but clumps of noise, hot pixels, and so on. The MOS looks at the remaining objects for such things as being on all
three images, the position change between images 1 and 2 and then images 2 and 3 is
proportional to the time changes between images, and that the three points lie close to a
line.
You can set various parameters to determine how faint the search goes, the maximum
motion allowed, and how close to a line the three points must lie. When done, you are
presented with the results in visual and text form. This allows you to confirm that what
the MOS found is real. The MOS also generates an MPC-ready report that you can edit as
needed before sending it to the Minor Planet Center.
No search engine is going to be perfect. Its going catch some things it shouldnt and
miss some things it shouldnt. Protecting against the former is why there is the visual
inspection tool. Unfortunately, there is no tool that finds excluded objects. In this case,
you make the search parameters less restrictive knowing that you will get more false hits
than before. At some point, you have to decide on the cutoff point between catching as
many new targets as possible and spending lots of time rejecting false hits.
The number of images is fixed at three. In the future it may be possible to use additional
images (you dont want to use less than three).
File Location
The files for this tutorial are found in
\MPO\EXAMPLES\MOVINGOBJECTS
which will be referred to as the MOS directory from here on.
The Moving Object Search form will be called the MOS from here on.
The images used in this tutorial were kindly supplied by Larry Robinson of Sunflower Observatory, Olathe, KS. They have been processed with darks and flats.
Open the configuration form and create a profile with the settings shown below.
You can access the MPC tab on the configuration form directly by clicking <MPC>
on the MOS form. The other pages are also available. You should confirm the MPC
settings before running any search.
245
MPC Settings
The Typ on the MPC page of the configuration affects the MPC report that is generated. What you select is not always straightforward. Generally, the Minor Planet Center
likes to make the identifications and so you indicate all targets as unidentified. In this
case, all the fields are of known main-belt objects, and so the MBA option in the drop
down list is not inappropriate.
The NET field is set to MPOSC3. This tells the MPC which catalog was used for finding
positions and magnitudes. Since the MPOSC3 astrometry is from the 2MASS catalog,
you could also enter 2MASS, which is more recognized by the MPC.
246
Regardless of which type you select, the MOS doesnt know which object youre
measuring, especially if its a new one. So, it uses a base name with sequential numbers to identify all objects in the report. The base name changes based on the option
you select. If necessary and desired, you can edit the MPC report so that it identifies
the specific object.
Consult the Minor Planet Center web sites astrometry guide for additional insights.
The settings on the other pages and tabs are not critical to this tutorial.
2.
Open the MOS by selecting Utilities | Moving Object Search from the Canopus
main menu, or press <Shift+Ctrl+M>.
3.
Copy the settings seen in the screen shot below, except the Save Dir field.
The MOS uses AutoMatch to extract the star data and stores the extracted data for
each image in its own file that can be read for the extraction process. The files are
stored in the Save Dir location.
The Mag Code drop down list includes the most common supported filters. If
you want to include magnitude estimates in the report, set this value to match the
instrumental band of the catalog magnitudes (the Default Band setting on the
Canopus Configuration | Photometry tab). If using a Clear filter, R is probably the best choice for most Kodak (KAF) blue-enhanced chips.
You can also select Do not use, but this means the magnitudes are not included at
all. Its easier to delete them in the editor if you change your mind than it is to go
through the search again. So, its recommended that you select an appropriate filter
when doing a search.
The FWHM setting is in pixels. If youre not sure of the value, move the MOS
aside and open one of the images. Click on a non-saturated star and review the
statistics in the bottom status bar of the Canopus form.
247
The SN:85/2.16 in this example indicates the star has a signal-to-noise ratio
(SNR or SN) of 85 and the full width-half max of the star profile is 2.16 pixels.
Click on several stars and use a mental average to enter into the FWHM field.
Its OK to go a little larger. Some of the images in the examples have slightly
greater FWHM and so 3.0 was used as a one size fits all estimate.
4.
Manually enter a valid directory in the Save Dir field or use the speed button
next to the field to display a directory selector.
5.
6.
Click <Add>. This displays a Windows file dialog that allows you to select more
than one file at a time.
7.
Locate the MOS directory and select the three files with a base name of A56. The
MOS Groups form should look something like this.
Click <+> or <> next to the group name to expand or collapse the tree that
shows which files are in the group.
8.
248
Repeat adding groups for all the images in the MOS directory. When done, you
should have something like the screen shot below.
Make sure that each group has exactly three images and that the three images are of
the same target (they have the same base name).
Dont worry if the files in the group arent in increasing date order. The MOS sorts the
files before doing the search.
9.
10.
The screen shot above shows the results of running this tutorial.
249
Click on the name of a group in the top list. This changes the image to the first
one in that group.
12.
Click <+> next to a group in the top list to show the targets in that group. The
numbering within the group starts at 1 but the number to the right is an incrementing number that corresponds to the objects in the MPC report.
13.
If there is more than one target in the group, clicking on each target loads the first
image of the group and the image is updated to show the position of the selected
target.
As you change targets in the top list, the bottom list is updated to show the three
images in the group. For each image, the date and time are displayed along with
the X/Y coordinate of the target on that image.
14.
Click on the name of a file in the bottom list to show its image and have the green
aperture re-center on the target.
16.
The information at the top of the report is taken from the MPC tab in the configuration form.
The memo can be edited and the form comes up in overwrite mode. This means that
if you start typing, you will replace existing characters and not force the characters
to the right of cursor to push right. This helps assure that you do not change the
structure of the data lines, which MUST have the data in specific columns.
17.
250
Note that the file has already been saved using the name given at the top of the
page. This is the file that you would attach if you click <Email> to send the report.
See the Minor Planet Centers web site for its Astrometry guide and the details of
this report format.
Search Recommendations
Dont set the search floor too low (closer to 0). Otherwise almost every bit of noise might
show up as a potential target.
Keep Min sigma at 1.0 or greater. This is the degree to which a pixel rises above the
general background. If you set it below 1.0, you quickly start adding noise as potential
targets.
Keep the tolerance at 2-3 pixels. Anything more from a straight line or away from the
predicted position on the image (part of the check for a real target), and you might get
false hits, especially in crowded fields.
The Max motion is for the full range of time among the images. For example, if the
first and last image are one hour apart, the object must move less than Max. motion
pixels over that time. If the images are 512 pixels on a side, a setting of 200 means the
object moves no more than about 2/5 the field over that hour.
251
252
3.
The Variable Star Search (VSS) is another example of putting images taken for one purpose to use for another. The images that are used in this tutorial were taken to generate
the lightcurve of an asteroid. The VSS was used to see if any of the stars in the field happened to be variable. If so, then it would be possible to work them further and generate
lightcurves and periods for them and even publish the results if worthwhile. As it turns
out, there was at least one previously unreported variable and a couple of other potential
candidates.
The concept of the VSS is similar to the Moving Object Search (MOS) in the previous
tutorial: start with a series of images, match them, and in this case look only for stars
that dont move and then auto-generate a lightcurve for every star on the images that
meets the initial search criteria. Those that show variability of a minimum level are reported for visual inspection.
You can give the profile any name you want. Note that its the same profile that was used
in the PhotoRed tutorial for finding first order extinction using the Modified Hardie
method. If you ran that tutorial, these settings should be the same, or very close to it.
The only setting on note on the Photometry page is for Circular apertures. This exercise does not worry about the types of magnitudes since they wont be stored in the observations file. If youre interested, the instrumental magnitudes of the stars as extracted
from the image are used for analysis.
This utility has been used to find at least a dozen new variable stars, some of which lead
to a published paper. It was inspired by talks with Robert Koff, Antelope Hills Observatory, Bennett, CO, as we drove from Colorado to California for the annual Society for
Astronomical Sciences meeting (http://www.SocAstroSci.org).
1.
253
2.
Open the VSS (<Ctrl+Shift+V>) or select Utilities | Variable star search from
the Canopus main menu.
3.
Match the search parameters that you see above in that section on the VSS, except
the Save Dir. Manually enter or use the speed button to the right of the field to
select a valid directory for the temporary files.
See the Reference Manual for a discussion on each of the parameters.
Unlike the MOS, these files can be used again during the same VSS session to do a
follow up search with different search parameters or comparisons.
4.
5.
Select all the files in the \MPO\EXAMPLES\VARIABLESTARSEARCH directory. This starts the search where the data are extracted from the images and moving objects are removed.
254
6.
When the data extraction is done, the Select Comparison Star form appears. Use
the screen shot above as a guide.
7.
Click on each of the numbered stars. As you do, its X/Y coordinates appear in the
list at right. If you click on the wrong star, highlight the entry in the list and click
<Delete>.
Picking Comps
Unless you know the field for some reason, you really dont have an idea of which stars
are variable and which are not. It is very possible when picking these three stars that one
or more is variable. The VSS provides a way to check, which is covered below.
For a general rule, try to pick isolated stars that are reasonably bright but not near saturation or non-linearity.
8.
9.
255
Oops! One of the stars appears to be variable! All three lines should be flat. Instead, Comp 2 goes a different direction from the other two. It is the most likely
culprit.
11.
Click <New Comps>. This displays the same form that you used before to select
comparison stars.
12.
Select the new set of comps (only #2 changes) and click <OK>. This time the
search is much shorter because all the original star data are still intact. The only
thing new is to compute the differential magnitudes between the comps and stars.
256
Click on the Results (Mean/SD) tab. This displays the search results using the
Max/Min S.D. parameters in the search parameters, meaning that the standard deviation of the average magnitude of each candidate was between the two limits (in
magnitudes).
At the upper left is the first image in which the variable appeared (it may not always be the first image in the set if the field of view changes significantly during
the run). The comparison stars are surrounded by red circles (cyan on an inverted
image) and the target surrounded by a green circle (magenta on inverted image).
At the upper right is a plot of the variable star. Below is a plot (SD Plot) of the
257
You can zoom and pan the SD Plot, which is often necessary to see individual data
points in crowded plots.
Note in the screen shot above that one of the stars in the SD Plot is a red triangle. This is
the currently selected star.
15.
Click on any point in the SD Plot; that point turns into a red triangle and the previous one returns to a black diamond. The image is reloaded and the newly selected stars location is indicated.
The lightcurve plot changes to show the selected stars data and the stats box at
lower right is updated to reflect the new data.
16.
If not already, check the Force Scale box above the lightcurve plot after you set
the value to 1.000 in the entry field next to it. This forces the lightcurve to have a
1-magnitude vertical scale. You do this when reviewing plots in order to avoid
seeing what are really small amplitude changes as being significant.
As you review the candidates, you will see general trends that repeat in the lightcurves, e.g., they all have a bump about half way through or have about the same
amount of scatter. Those that follow the trend are most likely not variables. Those
that show unusual behavior deserve closer attention.
17.
In the plot shows in the screen shot above, the lightcurve is relatively flat. This is
probably not a variable.
18.
258
Now you have something that appears to be legitimate! As it turns out, this was a
previously unreported variable star. Additional data were obtained and a model
was generated in David Bradstreets Binary Maker 3.0.
Note that the image was scrolled so that you could see the variable at the far right
of the image.
19.
Click on the Results (Targets) page. You should see the same lightcurve.
20.
21.
You can save this a text file for later reference. After youve reviewed the search
results, you can use the information in this file to relocate the stars and do proper
lightcurve photometry in Canopus. This is often necessary if you have a target
that is borderline case. The more accurate tools in Canopus will determine for certain if the star is variable and, if so, allow you to do a period analysis on the data.
259
As you saw, when you select a target on the Mean/SD page, it is also selected on the
Targets page. This reverse is true as well: when you select a target in the upper list on the
left of the page, it is selected on the Mean/SD page.
The screen shot above shows the Targets page with the Image tab selected. This allows
you to view an image just as before.
22.
Click on an item in the upper list on the Data tab. As you pick different objects,
their information appears in the bottom list and the lightcurve changes.
As before, look for general trends among lightcurves and then take a close look at
those that dont follow the trend.
23.
24.
If necessary, refer back to the beginning of the tutorial where the first set of comp
stars was selected, the set where the second comp was variable.
25.
Click on that star on the image. The selected item in the list changes and the information in the bottom list is updated.
26.
Click on the Data tab. You should see the lightcurve for that star.
260
While not as dramatic as the previous variable, the lightcurve on this star is
smooth. This is very likely a variable star.
Note that the information includes the RA and Declination of the star. You can
use the International Variable Star Index (VSX) on the AAVSO web site
(http://www.aavso.org) to look for known stars based on coordinates.
A search on the VSX found that the star is AX CAS, an Algol-type eclipsing binary. The magnitude range of the star is given at 0.5 mag. The dip here was only
0.15 mag so I happened to be lucky and caught the secondary eclipse.
261
Click on the Mag header until the values closest to zero (the faintest stars) are at
the top of the list. The instrumental magnitudes are usually lower (fainter) than
the floor you set. This is because all stars are extracted, regardless of the floor setting, which controls which stars are included in a search.
This means that you can use the same data for multiple searches, changing the
floor, minimum amplitude, and/or standard deviation ranges. This is different
from just selecting a new set of comparison stars because one is variable. Once
you have a good set of comparisons, use them again for different search parameters.
29.
30.
Click <Search>. Since you have one or more files in that list select, a message appears.
31.
Click <Yes> to use the extracted data. If you click <No> then the VSS will extract the data from the images again.
The extracted star data are not saved, so you cannot rerun searches with different
parameters. However you can use a different set of comparison stars provided the
image that was used to select the comparisons is still available in the same location
as when you ran the initial search.
32.
Click <Save> on the Select/Extract page. This presents a Windows file dialog.
Enter a file name and save the data. The extension is forced to VDF. The files
can be fairly large (by some standards). The file for this tutorial was about 2 MB.
33.
Click <Load> to reload the file. There are no files in the Extracted Files list, but
you are able to review the previous results and select a new set of comparisons, if
that original image is available.
You will not able to see the images unless they are available in the same location as
when the original search was run.
If you click <Search> after loading a VDF, the VSS will try to do a full search by asking
for a set of images. The original data files are no longer available, so the VSS must recreate the data.
262
This is when you know that youve picked some very bad comparisons.
263
264
4.
Exoplanet Lightcurve
Exoplanets are big news these days. More than 400 have been found as of this writing.
This tutorial is nearly identical in concept to the one where you found the lightcurve of an
asteroid and so will not go into great detail on the individual steps, except to direct you to
the necessary files and configuration settings.
The images for the exoplanet, WASP-1, were donated by Cindy and Jerry Foote.
Cindy is a noted backyard astronomer who has contributed to the discovery of
several exoplanets. The original images were cropped using the Crop function the
Batch Image Processor of Canopus.
Use the image below to identify the two comparisons and variable (exoplanet).
265
1.
Set the configuration in Canopus to match the settings above. You can give the
profile any name that you want.
2.
3.
From here on, only a file name will be given when needed and the WASP directory will refer to the above path.
4.
266
Create a new session in Canopus (<Shift+Ctrl+S>). Match the settings in the entry fields as shown below.
When working a variable star, heliocentric JD values are preferred, thus referencing all observations to when the distant stars light reaches the sun. If not used, timing errors of greater than 16 minutes can occur. Canopus always stores the raw JD
based on the FITS header. It is only when doing period analysis or when you export
data asking for such that either light-time or HJD corrections are applied.
You can tell Canopus to use HJD by selecting that option on the Photometry page of the
configuration.
Click <Save> and then <OK> after you have made the entries into the sessions
form. Since the two distance fields are 0, a warning message is displayed.
267
7.
Click <Yes> to continue. This message is more for those working asteroids who
forgot to click <Calc M/D/P> while creating the new session.
8.
9.
Finish the wizard and select all the WASP1_R*.FIT images in the WASP1 directory.
10.
Turn off AutoMatch in the Canopus Image List (uncheck the box) and select the
Auto | Simple method.
It took about 90 seconds to measure all the images on the tutorial machine, an older
one with only 512 MB of RAM running XP SP3.
11.
12.
268
Note the switch to Floating Range. The amplitude of the lightcurve is so small
that the default setting of Fixed Range = 0.8 mag shows a flat line.
You can clearly see the dip in the lightcurve caused by the exoplanet.
13.
14.
Reset the zoom to show the plot normally and rerun the search, but change the
Bin setting from 1 to 3. After the plot appears, zoom in again.
269
The Bin feature finds the average of Bin-successive data points with each point
no more than Max Diff minutes between the points in the bin.
This is not a running average binning. No given data point is used more than once.
15.
Before you move on, be sure to reset the Bin field to 1. Its easy to overlook
this and then be wondering whats wrong when you start plotting your next target.
The pair was less than 30 above the horizon on a not so good night and so the data are
noisier than I might have gotten otherwise given the brightness of the two.
270
5.
In the PhotoRed chapter, the AAVSO Batch Processing tutorial takes you through the
steps to get data into the Batch Reference File Generator (the Generator) using an image
of the variables field and catalog values from the MPOSC3 or other source.
If youre working a variable with an established sequence of comparisons, one that you
can retrieve from the AAVSOs Variable Star Plotter (VSP) then you can capture that
data and import it into the Generator in a few simple steps. This makes for much quick
work and assures that youre working with accepted and adopted values for comparison
stars.
1.
271
For this tutorial, well retrieve data for R AND. Not all stars will have a photometric sequence. The AAVSO is working on that as fast as they can.
2.
3.
4.
Select the Do you want a chart or list of field photometry | Photometry Table
radio button.
5.
Enter 20 for the Field of View (FOV). This is the approximate FOV, in arcminutes, for the telescopes at my location. You would normally enter a value appropriate to the system that you are using.
272
Select the How would you like the output | Printable radio button. This is important. The HTML option generates different underlying HTML code that the
page parser does not understand.
7.
Since youre not generating a chart for this tutorial, the other settings dont matter. However, if you want to generate a chart afterwards, at the least set the two
orientation buttons so that the chart matches your images.
8.
Its a good idea to generate and save the chart after getting the photometry data.
This way you can compare the Comp and Check assignments made by
PhotoRed and change them if necessary.
Click <Plot Chart>. This generates a photometry sequence table (only a portion is
shown in the screen shot).
Note that the table does not list the variables data other than name (if the chart was
generated by name) and RA and Declination. There are no magnitudes for the variable, which means that its color index will be set to 0.000.
9.
Use your browsers options to save this page as a simple HTM/HTML file. Do
not save it as a simple text file nor as an MTH (Web page archive), nor as Web
page, complete (IE option)..
10.
Close or minimize the browser and open the Generator in PhotoRed (Utilities |
Generate batch reference file from PhotoRed main menu).
11.
Set the Color Index to B-V. If your transforms were found using a different
color index, select the one you used. If youre not going to apply transforms, this
setting doesnt matter but you should get in the habit of selecting something.
12.
273
This gives you the option of including stars that do not have magnitudes in both
of the bands in the color index you selected.
13.
Click <No>, indicating that you want only stars that have both V and R magnitudes. A Windows file dialog is displayed.
14.
15.
This is a reminder that the target name may not have been found, e.g., you generated the chart by position, and that the color index for the target is was not set.
16.
17.
The Group should be -1. The Star Source setting is not important since youre
not working with an image and PhotoRed generated chart.
274
Chart ID
Note that, if possible, the chart (sequence) ID is parsed and entered into the form. This is
important information. Confirm that it was correctly entered. When you save the file, the
chart ID is saved as well.
Batch Reference Files created before version 10.4.0.11 do not have the chart ID. The
program will still read legacy files. You can update the file to include the chart ID by
loading it into the Generator, entering the chart ID, and then saving the file.
Always confirm the imported data. Parsing HTML code can be tricky, and theres always
a chance that the format of the output file can change. If it appears the parsing is no
longer working, please notify tech support, providing the information you used on the
VSP page to generate the photometry data.
275
276
6.
Canopus includes enhanced features for measuring double stars, saving the results, and
reporting those results. The summary reports can even reduce the saved raw instrumental
magnitudes in different filters into standard magnitudes, greatly increasing the scientific
value of your observations.
Heres an outline of what you do.
Use the Double Stars List to select all the images you want to measure. This List
is easier than having to use the Image | Open and locating the next file after
each measurement.
For each image, specify the reference (primary) star, specify the secondary star,
and then add the information to the DoubleStarsManagement table.
After all images have been measured, use the Double Stars Measurements form to
review the data and generate a report of the raw information for all data or that
matching a name and/or date range, or a summary report that can combine observations of a given star/filter/date into a single observation with errors and reduced
standard magnitudes.
Canopus must be able to AutoMatch the images (or you will have to match them
manually). This is so it has the plate constants to convert X/Y positions to RA and
Declination. See the Reference Manual for how Canopus looks for certain FITS keywords in order to generate a chart centered on the image center.
No Samples
There are no example images for double star measurements. These tutorials will be read
along only but the process is fairly simple and straightforward and should be easily applied to your images.
Double star measuring requires using AutoMatch, so open the configuration form
if necessary and create or select a profile that fits the images you will be measuring, e.g., the charting orientation is approximately correct and the core values on
the General tab match the focal length and pixels for the system.
2.
Set the measuring apertures as needed. You can change them while measuring, if
necessary.
277
Select Utilities | Double stars | Set double star list from the Canopus main menu
(or press <Shift+Ctrl+L>). This displays a Windows file dialog.
4.
Select one or more images to be measured. When done, the Double Stars List (the
List) is displayed.
5.
Double-click on the first file in the list (or highlight it and click <Load>). This
loads the image into Canopus. Each time an image is loaded, the highlight automatically moves to the next item in the list so that once you've measured the just
selected image, you can quickly load the next image in the list.
If successful, youll see the typical apertures around a number of stars on the image and numbers next to stars on the chart. If not successful, you will have to
match the chart manually. See the tutorials in the Core Operations chapter.
To clear the astrometry apertures so that you can see things better, right-click over
the image and select Toggle AA boxes from the popup menu.
6.
Once the image is matched, <Alt+Click> on the primary star on the image. This
prevents the Object Info box from appearing but still measures the stars position.
7.
Right-click over the image and select Set Dist/PA reference from the popup
menu.
8.
9.
Right-click over the image and select Save Dist/PA Data from the popup menu
to display the Double Stars Input form.
278
Enter the name for the double star in the Name field and select the appropriate
filter from the Filter drop down list.
11.
Check the Load DS List Image box if you want to have Canopus automatically
load the next image in the double star list.
In most circumstances, you want this box checked, which is the default, when using the Double Stars List. The exception would be when you have a multiple star
on a set of images, e.g., an AB, AC, and BC pair. In this case, uncheck the box so
that when you add a measurement, the current image remains loaded with its astrometric solution still valid.
For example, if you set the A star as the reference and measured the B as the secondary, you can then measure the C component as the secondary without having
to reset the A as the reference. Just be sure to change the entry in the Name field
to match which pair is being measured. For the B-C pair, you would have to reset
the B star as the reference and the C as the secondary.
12.
Repeat the steps above until the last image has been measured.
13.
The data were previously entered using the Double Stars list. Again, there are no sample
images, so you wont be able to perform these steps directly until you get your own images and have you own data.
For this tutorial, well generate a full report for A 1613AC.
14.
Check Define Filter | Name. This allows filtering the report by double star
name.
15.
Enter A 1613AC (without the quotes) in the Name field and then check the
Filter Active box.
279
Set the other controls as shown. The WAB in the Obs field is mine. You
would use another identifier for your observations.
17.
Click <Report>. This generates the text file and presents a Windows file dialog so
that you can change the name and or path where the file is stored.
18.
Save the file. This displays the report on the Report page of the DSM form.
19.
In this case, a summary report is being generated that applies PhotoRed transforms to the instrumental magnitudes to find the standard magnitudes for the
measured stars. This report also finds the mean and standard deviation of the position angle (PA, degrees) and separation (Dist, arcseconds).
20.
Click <Report>. Save the file and then review the results on the Report page.
This data set found PA = 139.2 0.19 and Dist = 60.88 0.094
There are many combinations to generate reports. This tutorial was designed to give you
an introduction into this useful tool that can used for simple but very important research.
280
7.
After finding a reasonably good solution for a lightcurve, you can use the data generated
when finding the period to generate a lightcurve for a future date. This is very helpful if
youre trying to fill in a specific part of the curve or planning observations around a
maximum or minimum. The Lightcurve Ephemeris tool (the LCE) in Canopus provides
this capability
This tutorial presumes that you did the lightcurve lessons in the Photometry chapter and
so have data for 771 Libera.
1.
Before you can run the LCE, you must do a period search and save the Fourier
data.
2.
Click <Save> on the Period Search Results form and save the file. I recommend a
name that includes Fourier, e.g., Libera_Fourier.txt.
3.
Do not change the contents of the file in any way. Otherwise, Canopus may not be
able to read it correctly.
Press <Shift+Ctrl+E> in Canopus (Utilities | Lightcurve ephemeris from the
main menu). This displays the LCE.
281
Click <Load Data>. Locate and load the Fourier data file that you saved.
5.
If the file is successfully loaded, the entry fields should change to reflect the data
in the file.
6.
Enter 771 Libera in the Name field, which is not automatically changed when
you load the data.
The Days field is set to generate approximately two complete cycles as based on
the period.
No matter how well the period has been determined, it is not reasonable (for an asteroid
at least) to generate a predicted lightcurve many years after the fact. For one, the period
you find for an asteroid is synodic, meaning that it will change slightly depending on
viewing aspects. Another is that even the slightest error can quickly amount to a significant difference over time. This tool is intended to be used during a current campaign to
plan observations.
However, you can turn things around and see how close the prediction comes to reality.
This may allow you to correct the original period by increasing the precision some. This
is what variable star people do all the time. Again, however, asteroids do not have constant synodic periods (and some even have changing sidereal periods) and so you must
keep things in perspective.
7.
Check the in UT box. This makes the horizontal axis in time instead of lightcurve phase.
8.
Click <Generate>. This displays the predicted lightcurve starting at the current
date and extending for about 0.5 days.
9.
Click <Save> to save the plot. This forces the LCE to full screen and then presents a Windows file dialog. After you save the file, the LCE returns to its original state.
282
8.
After youve done an AutoMatch to find an asteroid (see the AutoMatch tutorial on
page 17), you may wonder if there are any other asteroids in the image. You could use
the Moving Object Search (MOS) utility (see page 245) but there is a faster way to do a
preliminary check. This check can also be used to cross-check the MOS results. For example, you may have the overall sensitivity of the MOS too low and are overlooking
some asteroids.
1.
2.
Once you have the AutoMatch, right-click over the image to display the image
popup menu and select Other Asteroids | Locate
3.
This automatically starts looking at every asteroid in the MPCORB_EX file (compressed version of the MPCORB table). If the asteroid can reach at least magni-
283
Actually, two positions are found in that final computation. The first is for the midexposure of the image and the second is for 6 hours (0.25 day) later. If one or both
of the positions are on the image, the asteroid is plotted.
Heres an example of an image that had several asteroids enter or leave the field
of view between the time the image was taken and six hours later.
During the search, the bottommost status bar shows the progress.
On a Win7 64-bit machine, the search of almost 500,000 asteroids took less than
20 seconds.
284
4.
The screen shot above shows the final results. For each asteroid, the initial position is indicated by a square with the number (if not 0), name, and estimated V
magnitude above the square. If the square is too close to the top, the text is below
the square.
The +6 hour position is indicated by a circle. A line between the two shows the
path. In some cases, both markers may not be shown, the asteroid either leaving
or entering the image field of view over the 6 hours. In either case, the ID text appears above or below the one end marker that is shown.
285
If you want to keep this information for an astrometry report, click <Save> (see the tutorial on astrometry). Now you can measure and save the data for the other asteroid in the
image.
1.
Return to the Measurements page and click on the other asteroid that is within the
rectangle (this star was confirmed to be moving with the Blinker).
Note that the information in the status bars at the bottom update to show the objects
SNR and measured coordinates and position.
2.
Right-click over the image and select Set Asteroid Centroid from the popup
menu.
3.
Go to the Reductions page. Note the M-C values. They are huge. Thats because the computed position is not for the newly selected asteroid. You need to
update the information, including the name and MPC ID.
286
4.
In the Fixed Data section of the Reductions page, click the button next to
Number field, which currently has 771. This displays the asteroid picker.
5.
Make sure the Sort by combo box is set to Number and then type 5083
(without the quotes) in the Search for field. Let the table reposition to that asteroid (theres about a 0.25 second pause after you stop typing) and then click
<OK>.
6.
Click <Recompute>.
7.
The information on the Reductions page is not automatically updated after you pick
a new asteroid. You must click <Recompute> to complete the calculations and update the Object section data fields.
Review the information in the Object section. Note that the name and MPC
fields have been updated as have the magnitude and computed SNR. Most important is to note that the M-C values are much better.
287
As you did earlier, click <Save> and save the new data. By default, Canopus
should use a base name of A5083 for the file. Confirm this is the case so that
you know which asteroid is represented in the new data file.
If you had an image with several asteroids, measure the position and save the data
for each one following the steps above.
Copy to Clipboard
If you select Copy to Clipboard from the image popup menu while the markers are
displayed, the markers are included.
Lost List
Since the marker positions are very specific to the image and chart following an AutoMatch (or manual match), the asteroids list is cleared every time you generate a chart on
the Measurements page (either due to an AutoMatch attempt or generating the chart
manually from the menu) or if you load an image. It doesnt matter if the image is the
same image or the chart is for the same image on which you did the match. To keep
bookkeeping to a minimum, Canopus takes the blunt approach and clears the list.
288
9.
You can change the G or Mag Band data for multiple sessions involving the same
target in a set of simple steps. The general outline for doing so is the same for each: select
which parameter to edit, set the new value, select the sessions to edit, and make the edits
happen.
The G Field
The exact value for the Phase Slope Parameter (G) in the H-G magnitude system can
have an effect when finding the period of an asteroids lightcurve. When several sessions
in Canopus are combined, the data are normalized the phase angle of the asteroid at the
time of the first (usually earliest) session. To do this requires using the G parameter. So,
if the wrong value is used, the correction can be seriously affected.
The default value of G is 0.15. This is sufficient for most asteroids but some asteroids
have considerably different values of G. For example, high albedo objects (those of type
E, V, and R) have G ~ 0.43. Very dark objects (low albedo) have low values of G, ~0.12
for C, D, P, etc.
If the actual value of G becomes available after youve created a number of sessions for a
given object, its a good idea to enter the revised value and recheck the period analysis. It
can also help in period analysis if using the default value of G initially to try other values
and see if that improves the solution, e.g., better RMS fit. If you a have number of sessions, it can be very tedious, and prone to error, to edit each session manually. This is
where the batch editing feature for changing the value of G comes in handy.
Remember that the G parameter is not used for variable stars. Use this feature only
if the sessions involve an asteroid. Every session for a given object must use the
same value of G.
289
2.
3.
Click the button next to the G entry field in the session data section of the session form. This displays the asteroid selection form.
While its possible that Canopus might be able to determine the asteroid from the
entry in the Object field, there are too many possibilities and not enough certainties of how that field is managed to assure the correct interpretation. This is
why the Asteroid Lookup form is displayed.
4.
290
Locate the same asteroid as that in the session and click <OK>. This displays the
Sessions Selection form.
Select all sessions for the asteroid and click <OK>. This displays a data entry
form.
Remember! All sessions for a given asteroid must use the same value of G for analysis. The asteroids physical properties cannot and do not change from one night to
the next!
The initial value that appears in the entry field is the value of G in the session that
you highlighted in step 2.
6.
0.12 0.08
0.20 0.07
S, Q
0.24 0.11
E, V, R
0.43 0.08
The larger scatter in S, Q is because the various sub-classes of the S type each
have their own values of G. See the Icarus paper for a more detailed breakout.
7.
From here, you can redo the period analysis to see what effect, if any, changing the value
of G had on the data. Since the period solution can change significantly with different
values of G, be sure to do a search on an appropriate range of periods. Dont presume that
the new solution will be within a few percent of the old solution.
Open the sessions form and locate one of the sessions for the object for which you
want to set the Mag Band field.
Do not go into edit mode since batch editing is not allowed if you are editing
a specific session.
291
Click the button next to the Mag Field drop down list. This displays the Select
Sessions form.
3.
Check the boxes for those sessions where you want to set the Mag Band field
and then click <OK>. This displays a form to select the photometric band.
4.
Select the desired Catalog Mags (in this case, R) and then click <OK>.
5.
Check that the Mag Band field was updated to the new value in all the selected
sessions.
292
Notes
293
Notes
294