Site Map

Astrometry

Jump to:

Introduction

A detailed description of the astrometric calibration is given in Pier et al. (2003). Portions of that discussion are summarized here.

The r photometric CCDs serve as the astrometric reference CCDs for the SDSS. That is, the positions for SDSS objects are based on the r centroids and calibrations. The r CCDs are calibrated by matching up bright stars detected by SDSS with the UCAC astrometric reference catalogs.

Stars detected on the r CCDs are matched directly with stars in the United States Naval Observatory CCD Astrograph Catalog (UCAC2, Zacharias et al. 2000), which has a precision of 70 mas at its catalog limit of r= 16, and systematic errors of less than 30 mas. UCAC2 extends up to around a declination of 41 degrees north. Outside the UCAC2 area we use an "internal" UCAC data release known as "r14". Together UCAC2 and r14 cover the whole sky. There are approximately 2 - 3 magnitudes of overlap between UCAC and unsaturated stars on the r CCDs. The astrometric CCDs are not used.

The r CCDs are calibrated directly against the primary astrometric reference catalog. FRAMES uses the astrometric calibrations to match up detections of the same object observed in the other four filters. The accuracy of the relative astrometry between filters can thus significantly impact FRAMES, in particular the deblending of overlapping objects, photometry based on the same aperture in different filters, and detection of moving objects. To minimize the errors in the relative astrometry between filters, the u, g, i, and z CCDs are calibrated against the r CCDs. Each drift scan is processed separately. All six camera columns are processed in a single reduction. In brief, stars detected on the r CCDs , are matched to catalog stars. Transformations from r pixel coordinates to catalog mean place (CMP) celestial coordinates are derived using a running-means least-squares fit to a focal plane model, using all six r CCDs together to solve for both the telescope tracking and the r CCDs' focal plane offsets, rotations, and scales, combined with smoothing spline fits to the intermediate residuals. These transformations, comprising the calibrations for the r CCDs, are then applied to the stars detected on the r CCDs, converting them to CMP coordinates and creating a catalog of secondary astrometric standards. Stars detected on the u, g, i, and z CCDs are then matched to this secondary catalog, and a similar fitting procedure (each CCD is fitted separately) is used to derive transformations from the pixel coordinates for the other photometric CCDs to CMP celestial coordinates, comprising the calibrations for the u, g, i, and z CCDs.

Note: At the edges of pixels, the quantities objc_rowc and objc_colc take integer values, contrary to standard practice.

[Back to top]

Proper Motions

The SAS and CAS include proper motions for objects derived by combining SDSS astrometry with USNO-B positions, recalibrated against SDSS (Munn et al. 2004). These are given in the ProperMotions table in the CAS, and in the "PM" external catalogs directory in SAS.

The proper motions in DR9 and later correct an error in DR7 proper motions for stars at low Galactic latitude; see the QA discussion below for more details.

[Back to top]

Calculating Errors for Individual Objects

The calibrations are performed in great circle coordinates. The estimated errors in the calibrations are given on a per-frame basis. The calibration errors in great circle longitude and latitude are given by the attributes muErr and nuErr, respectively (in arcseconds). These are in the photoField files in the SAS. These should be added in quadrature with the centroiding errors for individual objects to give the estimated total error in the position of a given object. The centroiding errors in great circle longitude and latitude are given by the attributes objc_rowcErr and objc_colcErr, respectively (in pixels; these should be multiplied by the focal plane scale of 0.396 arcseconds/pixel to convert to arcseconds). These attributes are in the photoObj files in the SAS.

[Back to top]

Calibration Equations

Astrometric calibrations are generated as a separate set of equations for each frame, converting frame row (x), frame column (y), and star color to catalog mean place great circle longitude (μ) and latitude (ν), in degrees:

for color < (color)0:

x' = x + g0 + g1 y + g2 y2 + g3 y3 + px color
y' = y + h0 + h1 y + h2 y2 + h3 y3 + py color

for color > (color)0:

x' = x + g0 + g1 y + g2 y2 + g3 y3 + qx
y' = y + h0 + h1 y + h2 y2 + h3 y3 + qy

μ = a + b x' + c y' ν = d + e x' + f y'

Note that in these equations, for DR8 we did not account for the color term at all, which results in 10 to 20 mas systematic errors. However, for DR9 and subsequent releases, we accounted for this color term correctly.

The transformation from (x, y) to (x', y') corrects for optical distortions (which, in TDI mode, are a function of column only) and differential chromatic refraction (DCR). For u and g frames, DCR is modeled as a linear function of color (u-g for u frames, g-r for g frames) for blue stars [(color)0 = (u-g)0 = 3.0 for u frames, (color)0 = (g-r)0 = 1.5 for g frames], and a constant for redder stars. For r, i, and z frames, DCR is modeled as a linear function of color (r-i) for all stars [(color)0 = (r-i)0 >> 1]. (The DCR corrections are mis-stated in Pier et al. [2003], where [r-i]0 appears in the equations rather than the correct [color]0, and where the wrong value for [color]0 is given for u frames.) The corrected frame coordinates (x', y') are then transformed to catalog mean place great circle coordinates (μ, ν) using an affine transformation.

The calibration coefficients may be found in the photoField files in the DAS, where the attribute names are different than given in the transformation equations above; (color)0 is called riCut; g0, g1, g2, and g3 are called dRow0, dRow1, dRow2, and dRow3, respectively; h0, h1, h2, and h3 are called dCol0, dCol1, dCol2, and dCol3, respectively; px and py are called csRow and csCol, respectively; and qx and qy are called ccRow and ccCol, respectively.

[Back to top]

Transformation from Great Circle Coordinates to J2000 Celestial Coordinates

The calibration equations above yield catalog mean place in great circle coordinates. To convert these to J2000 celestial coordinates you need to know the right ascension and inclination of the ascending node of the scan great circle with respect to the J2000 celestial equator. These are given as the header keywords "NODE" and "INCL", respectively, in the "photoField" file. The celestial coordinates are then

tan(α2000 - α0) = [sin(μ - α0)cos ν cos i - sin ν sin i]/[cos(μ - α0)cos ν]
sin δ2000 = sin(μ - α0)cos ν sin i + sin ν cos i

where μ and ν are great circle longitude and latitude, α0 and i are the right ascension and inclination of the ascending node of the great circle with respect to the J2000 celestial equator, and α2000 and δ2000 are J2000 right ascension and declination.

Astrometry QA files

We have implemented a new astrometry quality assurance system in order to identify errors in the SDSS imaging astrometry better. The astromqa data model fully describes all of the files that are produced in this process. Here we outline the major results and plots that are produced, both the global plots and those appropriate for each run.

The astrometric QA is defined with respect to a set of reference catalogs:

The astrometry QA results are summarized in the following way:

Caveats on DR8 astrometric calibration

The DR8 astrometric calibrations were substantially degraded relative to the DR7 astrometric calibrations, particularly at declinations northward of about 41 deg. These errors have a much smaller, but non-zero, effect on the DR8 proper motions. We recommend using the DR9/DR10 astrometry and proper motions, which are available in the DR9 or DR10 releases, or in the astromDR9 and properMotionsDR9 tables of the DR9 CAS.

A detailed description of these errors was released as part of the DR8 Paper Erratum, and is also available on the DR8 version of this algorithms page.

Note in particular that the proper motions tabulated in the CAS were only mildly affected by these problems. The primary effects on the proper motions are to introduce an additional systematic error with color of order 0.5 mas/yr, and to introduce an additional source of error for stars with declinations above 41 deg, of order 1 mas/yr.

[Back to top]