Two examples of astrometry

Michael Richmond
Jan 16, 2009

Today, our topic is astrometry, the science of measuring positions of celestial objects. You will have two different little projects to do today, to help illustrate two different facets of this field.

  1. In order to determine the parallax of an asteroid, you will need to make simple relative measurements of the position of an asteroid relative to a few nearby stars. Because the field of view is small, you can use simple plane geometry and get an answer which is good enough.
  2. In order to calculate the (RA, Dec) positions of a large number of stars in a large field, you'll need to deal with the conversion from a spherical coordinate system to a plane system and back again. Fortunately, there are plenty of tools you can use for this job.


Get connected

  1. Log into computers in the CIS lab.
  2. Start the "Terminal" program
  3. Type the following:
                ssh -Y -l student spiff.rit.edu
    
    When asked for a password, supply the value I will tell you.
  4. You should now be connected to a computer in my office, which has some images of a cluster of stars.
  5. Type the following commands:
                propinit
    
                cd grant         (or your own first name here)
    
                ls
    

You should now be in the directory with the images you'll need. The two images you'll use for the parallax experiment are

           rit.fit              an image taken from the RIT Observatory
           usna.fit             an image taken from the US Naval Academy

You can display one of the image by using the tv command, like so:

           tv rit.fit zoom=2

A new window should pop up like this:

The zoom=2 argument zooms into the image, making it twice as big as it really is. You can use zoom=0.5 to make an image appear smaller, or leave out the zoom keyword entirely to see the image at its normal size.

If you move your cursor into the image, you should see the coordinates (bottom of the window) change as you move the cursor. If the coordinates DON'T change, then click on the window's bar at the top. If they STILL don't change, left-click once while pointing anywhere in the image.

There are three keys you can use to measure stellar properties.

To destroy an image window, you can right-click while the cursor is inside it, or use the Mac red-color-button at the top left of the window.


Measure the parallax of an asteroid and compute its distance

One way to measure the distance to objects is via parallax, which involves measuring the apparent shift in position of an object when it is viewed from two different locations. You are probably most familiar with heliocentric parallax, in which measurements of a star are made 6 months apart and the orbit of the Earth provides the baseline.

In this case, however, two images were taken simultaneously. I took one from the RIT Observatory, and a colleague that US Naval Academy in Annapolis, Maryland, took the other. You can use the locations of our two observatories to determine the baseline for the measurements.


                             latitude        longitude
    ---------------------------------------------------
    RIT Observatory          43:04:33         77:39:53

    US Naval Academy         38:58:41         76:29:31

In order to determine the angle by which the asteroid appeared to shift, you need to measure the position of the asteroid relative to some other star in each image. You'll also need to determine the plate scale of each image -- in other words, the relationship between pixels and arcseconds. You can use this chart as a guide:

In this chart, stars B and D are 704 arcseconds apart.

  1. how many pixels apart are stars B and D in the RIT image?
  2. how many pixels are in one arcsecond in the RIT image?
  3. how many pixels apart are stars B and D in the USNA image?
  4. how many pixels are in one arcsecond in the USNA image?

The next step is to measure the location of the asteroid relative to some reference point in each image. One approach is to set up a coordinate system using the line connecting stars B and D. If you measure the vectors B-D and B-asteroid, you should be able to compute the x-component and y-component of the vector B-asteroid. Do that in both images, and you can then find the shift in the position of the asteroid.

  1. How far is the asteroid from star B in the RIT image? Express your answer in arcseconds.
  2. Break this distance in an X-component and Y-component, so you can write the distance between the asteroid and star B as a vector.
  3. How far is the asteroid from star B in the USNA image? Express your answer in arcseconds.
  4. Break this distance in an X-component and Y-component, so you can write the distance between the asteroid and star B as a vector.
  5. What is the shift in position of the asteroid between the two images? Write the magnitude of the shift in arcseconds.

Once you have the parallax angle θ, you can compute the distance to the asteroid.

What is your value for the distance to this asteroid?

Note that there are some complications that I've ignored in this exercise, which mean that your value is only a lower limit to the distance. Can you explain what these complications might be?


Computing the positions of many stars in a large field

If we take a picture showing a large area of the sky, then we run into a big problem: the positions of stars will be distorted, because our telescope projects the image of a curved, three-dimensional sky onto a flat, two-dimensional surface.

This distortion is the same reason that maps of large portions of the Earth's surface can look strange. This map make Greenland look as large as Africa, but it actually has less than one-tenth the area!


This Mercator projection of the Earth courtesy of Stefan Kuhn and Wikimedia Commons

So, one ingredient in a successful astrometrical recipe is a good method for accounting for the distortions involved in projecting the celestial sphere onto a plane.

Another ingredient is a good set of reference stars with known positions. In the optical regime, the best astrometric catalogs at the present time are

When you find yourself with an image of some portion of the sky, and you want to compute the positions of stars in the image, you can follow these steps:

  1. measure the positions of many stars in the image itself, in pixel coordinates (x, y)
  2. starting with one of the astrometric catalog, extract a list of stars with known positions in that area, with positions in (RA, Dec)
  3. project the (RA, Dec) positions of the catalog stars onto a plane, into what are called the "standard coordinates" or projected coordinates (xi, eta)
  4. match the catalog stars to stars detected in the image
  5. using the matched pairs, compute a coordinate transformation which transforms pixel coords (x, y) into projected coordinates (xi, eta)
  6. apply the transformation to ALL the detected stars, taking them all into the (xi, eta) system
  7. de-project the plane coordinate back into the spherical coordinate system (RA, Dec)

Today, I'd like simply to illustrate the procedure, so I'll do most of the work for you. We'll start with an image of this pretty galaxy: m51.fit.

I've measured the positions and instrumental magnitudes of the brightest stars in this image. You can find them in the file m51.pht. The columns are

  1. an index number
  2. row coordinate = y-coord (pixels)
  3. col coordinate = x-coord (pixels)
  4. local sky value
  5. uncertainty in local sky value
  6. instrumental magnitude
  7. uncertainty in instrumental magnitude
  8. flag indicating how "clean" the measurement is (0 = good)

I've also gathered the (RA, Dec) positions of many stars in this region from the USNO B1.0 catalog. You can find these positions in the file m51cat.dat. In this file, the important columns are


  0 - an index number 
  9 - the Right Ascension (decimal degrees)
 10 - the Declination (decimal degrees)
 11 - the R-band magnitude 

So, given all this preliminary work, your job is pretty easy.


For more information