Astronomy Lab
Determination of Stellar Magnitudes
from Digital
Images - Revised to use DSLR images
OBJECTIVES
- Learn and understand stellar
magnitude concept.
- Practice
calculating
stellar
magnitudes with
spreadsheet.
- Calculate
intrinsic
magnitudes of
stars based
on distance from the star.
Note: We
will do Procedure I on one day, then return later to Procedure II to
calculate
magnitudes.
INTRODUCTION.
Astronomers
describe the apparent brightness of stars in magnitudes
(see Comins p. 179f). Hipparchus developed the magnitude scale in
ancient
Greece. Arbitrarily, he said that a 6th magnitude star was 100 times
fainter
than a 1st magnitude star. A magnitude difference means
a
certain ratio in intensities. Let r be
the
ratio of intensities for stars differing by one magnitude. r x r
= r2
is the ratio of intensities for stars differing by two magnitudes, etc.
r5
= 100. Thus, 
The magnitude scale is logarithmic,
i.
e. the magnitude of a star is related to the exponent of the
intensity. Another strange characteristic of the magnitude scale is
that
it is backwards. A magnitude 6 star is 100 times fainter
than a magnitude 1 star. Thus magnitude is a logarithmic measure of
"faintness".
Suppose we
have a standard star of magnitude mo which yields an
intensity reading of Io from a photometer. Another
star
reads I on the same photometer. A simple formula for the
magnitude
m
is given by
m = m0 - 2.5*log(I/I0)
Although
this formula is rather cryptic involving logarithms, it presents no
problem
for a calculator or computer spreadsheet.
PROCEDURE
I. Measure brightness with digital imaging. The class has
previously
photographed a number of stars of varying brightness using a digital
SLR camera. Cassiopeia, Fall,
2008. You will load digital images of these stars on
machines in the Spidel Computer Laboratory or Physics Laboratory using
Matlab. Matlab is a sophisticated and powerful analytical package
that is programmable. In future years, a photometry package that
runs without the expensive Matlab program may be available.
Click here to examine a small portion of a
digital star image.
We will
measure selected stars in the constellation Cassiopeia -
photographed Oct. November 6, 2008 by Gracie McCarroll, Chandler Jones,
and Bayle Owens.
- Examine the image :Cassiopeia, Fall, 2008. Write
down the star name that you think is the brightest star by examining
the negative image. 2nd brightest star? Faintest of the
labeled stars? We will then measure the brightnesses.
- Start Matlab. This
is a big program, please be patient!
- From Matlab, just
type starbright. This is an application developed
by Donald Collins to measure the brightness of stellar images by "point
and click" techniques with the user's mouse. The camera used to
make the images is a digital SLR camera using raw images. Raw
images display a linear light response.
- starbright
asks for a file. If you are on WWC campus,
navigate to
G:\physics\Astronomy\2008Cas_Owens_McCarroll_Chandler and
select any of the three images that end in ".fit". starbright loads all three B, G, and
R images, and it doesn't matter which one the user selects. These
images were made by adding 20 digital images at 15 sec each.
Total exposure time of 10 x 15sec = 300 sec.
- Adjust the contrast to your
satisfaction. Use the contrast control on the image
window. Adjust the contrast
so that the faint stars are visible, but the image doesn't look too
noisy. You may want to try a Negative image. You can
re-adjust the contrast and pos/neg at any time throughout the session.
- NOTE: DO NOT ENTER NUMBERS FOR
THE DISTANCE, OR CALIB STAR . DO NOT PUSH THE "ENTER CALIB STAR
MAG" BUTTONS. THEY ARE FOR ANOTHER PURPOSE AND WILL INTERFERE
WITH TODAY'S MEASUREMENTS.
- Take the mouse and click on a
star in the
image. Note:
In order to view
the image better and see the stars better, you should zoom up with 1
click of the zoom button. Then you will have to pan the image to
get each star in the field of view. To calculate the
brightness of a star, the computer reads the sky background in an
annular region surrounding the star. The average background is
subtracted. Then the net value of all the pixels in a small
aperture that includes the star are added to give the brightness in
ADU's. In case the mouse position is not accurate, the computer
tries to
optimize
the brightness count by moving the aperture around in a small area of
the
image until it finds the largest brightness signal. The
brightness is displayed in another window as Analog Digital Units:
ADU. If the user misses a
star by such a large amount that the search algorithm fails,
nothing is displayed. The user should try
clicking more accurately.
- Record the starbrightness in a
table in your notebook along with the star name, distance, and Star
brightness (ADU). See this
link for an image of Cassiopeia (2008) and the star names and distances.
ADU
means "Analog Digital Units" and represents the sum of all the
pixels added up in the aperture (as was done with the pseudostar lab).
- When
finished
with all the stars, press the "Done" button in the star
image. You should also exit Matlab and log off the computer.
PART
II Use spreadsheet to calculate the magnitudes. After the
table
of the star distances and brightnesses has been completed, enter all
the data into a table in a spreadsheet. All computers can run the
spreadsheet (excel).
- In the spreadsheet you should enter your name in tht top left cell.
- Enter
column
headings for "starname", "distance (LY)", "brightness (adu)".
- Enter
the
starnames, distances, and brightness from starbright.
- Save
your
work!
- Enter
Reference
Star Delta Cas in
isolated row.. In a row several rows separated from
the main block of data, enter the
name, distance, and brightness for the reference star: Delta.
This is important because all the magnitudes will be calculated using
Delta as a standard. The magnitude for Delta is given as
2.7. The separate isolated row for the delta star is also needed
to sort the data and not mess-up the absolute addressing.
- Calculate
the
apparent magnitudes using Delta as a standard. This
will be done in a new column. In a new column next to brightness,
enter a label for "magnitude".
- Go to the first
cell of the new column to be calculated, type the equals '=' key
with
no return, then type the formula for calculating the magnitude by
entering numbers and pointing to cell values. This is best done
by means of example. The formula is:
m = mDelta - 2.5*log
(I/IDelta)
or m = 2.7 - 2.5*log
(I/IDelta)
For the "I" in the
above formula you point and click to the cell containing the ADU
brightness values. This will be demonstrated in class.
- Be very careful how you point to IDelta in
the above formula. Remember that you entered another row for IDelta.
You must use "absolute" addressing. When
you point (or highlight) the cell coordinates in the formula, you
should "absolutize" the coordinates. This is done using the
function key <F4> in Excel. You will notice that the cell
coordinates have a dollar sign ($) placed in front of the row and
column numbers. If the <F4> doesn't work you may
explicitly type the $ signs in front of both row number and column
number in the formula. The reason for this will become apparent
in the next step of replicating the formula.
- Once the
formula is entered correctly in the top cell, the cell may be easily
copied
down the rest of the column. The cell has a little "dragbox" in the
lower
right. Just click and drag the cell using the drag box, and the formula
gets copied for the whole column. This will require some practice.
- Check to make sure that the
formulas are correct.
- Format the calculated cells so
that the numbers are rounded to 2 decimal places.
where d
is the distance from the earth
in LY.
Enter the column heading and the formula for the first line of
data. Copy the formula down the column as before.
- Which star is the brightest
apparent magnitude?
Which star is the second brightest? Which star is the faintest?
How do these compare with your estimates at the beginning?
- Which star is the brightest
Intrinsic magnitude? Which star the second brightest?
Which star the faintest?
- Save your work!
- Calculate the Luminosity
relative to the sun. The luminosity of a star is the total
power radiated from the star. It is convenient to compare the
luminosity in solar units. If a star has a luminosity of 10
solar units, that star is 10 times brighter than the sun when the sun
and the star are viewed from the same distance. The intrinsic
magnitude of the sun is 4.83. For a star with an intrinsic
magnitude of M, the luminosity in solar units is given by thr formula:
Where 4.83 is the
intrinsic magnitude of the sun. In this formula, (4.83 - M)/2.5
is the exponent, or the power of 10. This is also taking the
antilog of a number. When typing this strange formula in the
spread sheet, you have to write it in the following format:
10^((4.83-M)/2.5). Use the pointing feature of the spreadsheet to
point to the intrinsic magnitude cell. You must be very careful
with the parentheses. Note that first you subtract numbers, then
divide the difference by 2.5, then rais 10 to the result! The
brightest star should be larger than 30,000 times brighter than the
sun. The "^" operator means "raise to a power".
- Which star has the largest
luminosity? How many times brighter than the sun is the brightest
luminosity?
- Which star has the smallest
luminosity? Which stars emit less total light then the Sun?
- Print the
spreadsheet. Make sure the spreadsheet is correct and your
name is at the top.
- Optional: Sort the spreadsheet
before printing in order of apparent magnitudes and again in order of
luminosity. You need to be careful here. Be sure to select
all the data, and only the data, before sorting. Make sure all
rows become sorted. You should certainly save the spreadsheet
before sorting so a mess-up can be restored by reloading the saved file.
Summarize
the data. Discuss the results, and any surprises
you observe.
Compare the brightest, second brightest, and faintest star (as viewed
from earth) with your "eyeball" predictions made earlier. Why is
the star with the greatest intrinsic magnitude not very bright in the
photograph? For the star(s) that is(are) fainter
than the sun where does(do) it(they) rank in the apparent
brightness?
Several stars in Cassiopeia are
very interesting. Gamma Cas is a variable star with very
interesting properties. See the link to the AAVSO
discussion. Beta Cas is another variable star. Bright
variable stars are very difficult for professional astronomers to
measure because they saturate the detectors. A possible research
project would be for a student to measure the brightness of some
of these stars over an extended period using the digital camera
technique. This would require special attention to detail,
correcting carefully for flat fielding, etc., but students should be
able to make a research contribution to science.
Cassiopeia

|
Star
|
Distance (LY)
|
alpha
|
229
|
beta
|
54.5
|
gamma
|
613
|
chi
|
204
|
delta
|
99
|
epsilon
|
442
|
eta
|
19.4
|
kappa
|
4129
|
lambda
|
355
|
mu
|
24.6
|
theta
|
137
|
upsilon1
|
229
|
upsilon2
|
206
|
zeta
|
597
|
|
Cassiopeia - Nov. 6, 2008. Sum of 20 exposures at 15
sec. using DSLR Camera "piggy-backed" onto Meade Clock-drive
telescope. Focal length: 55 mm. Photo made by Gracie
McCarroll, Chandler Jones, and Bayle Owens.