Mapping Toolbox    

Astronomical Data

The Mapping Toolbox provides an interface to a set of astronomical data called FK5, the Fifth Fundamental Catalog of Stars. This catalog consists of positions, proper motions, and other characteristics for over 4500 stars. The data is provided by the Astronomical Data Center (ADC), located at the NASA Goddard Space Flight Center, and can be retrieved over the Internet at the ADC's home page at <http://adc.gsfc.nasa.gov/> or by anonymous FTP at <ftp://adc.gsfc.nasa.gov/>. Documentation and other information can be found at these sites as well.

The data comes in two files, `FK5.dat' and `fk5_ext.dat'. You can read both files using the readfk5 interface function. There are more than four thousand stars in the catalog, so reading and processing the data takes a while.

The star catalog is processed into MATLAB structures by the external interface function. The fields of the structures contain the different types of information for each star, while the structure elements represent the individual stars.

To view the fields of the structure, type the following:

If you just want to display the locations of the stars, you can find position and visual magnitude data in the stars workspace as part of the Mapping Toolbox atlas data. The remainder of this section illustrates the use of Mapping Toolbox functions to manipulate and transform star data for display.

The star positions are given in terms of right ascension and declination. You can display a star map by converting the positions to latitude and longitude coordinates. Right ascension is given in hours, minutes, and seconds, while declination is in degrees, minutes, and seconds.

You probably want the map to look like you are inside the celestial sphere looking out, so switch the sign of the longitudes:

It is customary to display the stars with sizes proportional to the visual magnitudes. Extract the magnitudes from the data structures:

Faint stars have large visual magnitude numbers, while bright ones have smaller values. The magnitudes of particular bright stars were used to define the visual magnitude scale. It was later determined that certain stars were even brighter than the reference stars, resulting in negative visual magnitudes. You will need to convert the values to relative intensities for plotting.

Stars in equatorial coordinates are typically displayed in the Equidistant Azimuthal projection. Here is the southern sky. You can see the southern cross at about 60° S, 175° E.

The sky is often depicted in galactic coordinates, with the center of our galaxy as the origin of a polar coordinate system. You can use the Mapping Toolbox to transform the star positions from equatorial coordinates to galactic coordinates. First define the location of the galactic pole and the galactic center.

Next create a new origin vector based on the location of the new pole in the old equatorial coordinate system. Use it to rotate the locations of the pole and the center of the galactic coordinates. Note that the origin provided by the newpole function results in a system centered on the location of the new pole in the old coordinate system. Rotate the star positions and bring the galactic center to zero longitude.

The stars in galactic coordinates are commonly displayed in a pseudo-cylindrical projection like the Hammer.


  U. S. Gridded Elevation Data Importing Other Data