slaGeoc()function. When passed the geodetic latitude and height above sea-level, this returns the radial distance from the rotation axis, Rr, and the distance from the equatorial plane, denoted by
z. The distance from the center of the Earth is then given by
The altitude returned by GPS probably won't be accurate enough for tracking the moon to arcsecond precision. In particular it may (?) be referenced to a different reference spheroid than
slaGeoc(). If this ever becomes an issue, the altitude will have to be determined more accurately, perhaps by adjusting it interactively to produce stationary fringes while tracking the limb of the moon. For other, more distant sources, the accuracy of the GPS altitude should be more than sufficient.
The symbols shown in this diagram are used throughout this document. They are defined as:
Similarly the cosine formula of spherical trigonometery yields:
and the analogue to the cosine formula yields:
The above equations can be combined to find the horizon coordinates that correspond to given equatorial coordinates. The resulting equations are:
For the CBI we also need to know the rates at which A,E and P evolve per second of UT. For a sidereal source this is determined solely by the rotation of the earth, so the rates are related to the steady increase in the hour angle of the source. The rates due to precession are insignificant. For a nearby source an additional contribution comes from the apparent proper motion of the source. This appears as a steady change in and . For the computation of tracking rates the proper motion terms can be ignored for all standard astronomical sources except the moon. However it is possible that ephemerides could be used to specify unusual tracks, such as raster scans of the sky, so in the CBI, the proper motions of all sources that are implemented as ephemerides will be incorporated into the tracking rates. The proper motions of non-ephemeris objects will be ignored.
The tracking rates can be found by differentiating the above equations wrt hour angle, then multiplying the results by the sidereal rate.
where and represent the proper motion of the source (if any), and:
Note that d/dT is the ratio of the rotational period of the Earth to the length of the UT1 day. According to the explanatory supplement of the Astronomical Almanac, when the ratio of mean sidereal time to UT1, 1.002737909350795, is combined with the precession rate of right ascension, the rotation rate of the Earth is reasonably aproximated as 1.002737811906 rotations per UT1-day. In order to compute tracking rates in radians per UT1 second, the appropriate value for d/dT is thus.
Application of the sine rule to the above triangle reveals,
Solving this for the angle of parallax results in,
This correction will be applied to solar system objects only.
The CBI weather station will provide real-time measurements of the
local temperature and relative humidity. This will be used to estimate
conditions in the higher parts of the atmosphere and thus derive an
estimate of the elevation offset due to refraction. Rather than
devising an empirical model of the atmosphere, the CBI will, at least
initially, adopt the Starlink
slaRefro() function in
Patrick Wallace's slalib library. This was originally developed with
optical pointing in mind, but it was subsequently retrofitted for
radio wavelengths. It remains to be seen whether this will be adequate
for the CBI. If not, an empirical model will have to be
devised. Unfortunately the conditions in Pasadena are radically
different from those on Chajnantor, so initial testing at Caltech
won't necessarily resolve this issue.
Given information about local atmospheric conditions
slaRefro() returns two coefficients,
B. These relate the vacuum elevation to the observed
refraction via the following equation.
We need an equation that goes the other way, from the vacuum elevation to the observed elevation. Unfortunately the above equation is transcendental, so in Starlink User Note 67, Pat Wallace suggests the following approximation.
This can be expanded to
Because of this effect, the rotational velocity of the surface of the earth results in a small diurnal shift in the apparent positions of sources. The maximum aberration at the latitude of the CBI is about 0.3 arcseconds. Note that the much larger effect caused by the motion of the Earth around the Sun is assumed to have been taken care of during the precession calculations used to find the apparent Right Ascension and Declination of the source.
The tangential rotational velocity at the surface of the Earth is given by,
Where Rr was defined earlier to be the radial distance of the CBI from the center of the Earth, and 86400 is the number of seconds in a day.
To determine the adjusted azimuth and elevation, the following diagram shows a unit vector pointing at the source from the center of the celestial sphere, overlaid with a rectangular coordinate system that is aligned with the compass points and the zenith. A second unit vector illustrates the effect of diurnal aberration, and the accompanying equations show how to calculate its direction.
Nonetheless, there will be some flexure in the mount itself and in the azimuth bearing. This should be a function of elevation. Traditionally one assumes that flexure is proportional to the torque applied by gravity to a hypothetical telescope tube. This can be corrected using
where kf is an empirical constant to be fitted for. It remains to be seen whether this will be appropriate or necessary for the CBI.
From this diagram, given a topocentric source position (A,E,P) we need to find the new azimuth, , the new elevation, E' and the new paralactic angle, P'=P+P, that correspond to the tilted coordinate system of the telescope.
By applying the sine rule of spherical trigonometry to the lower left and rightmost triangles, and the 4-parts formula to the lower left triangle, one finds that:
from which one can solve for E'.
By applying the four-parts formula separately to the left and right triangles, one obtains two equations for tan(C), which in turn can be solved for the following equation for tan:
where the unknown, tan, is given by applying the 4-parts formula again to a different combination of parameters of the left triangle, resulting in:
Substituting this into the previous equation gives:
This can then be used to compute the new azimuth, ,
By applying the 4-parts formula again to two different sets of parameters of the rightmost triangle, plus the 4-parts formula to the lower left triangle, one can also obtain:
Using these equations to solve for P, gives:
which after having its denominator and numerator multiplied by results in:
Thus the apparent parallactic angle of the source is:
As mentioned above there are two main parts to the tilt.
To help clarify this diagram, note that to first order, tilting the front of the azimuth platform upwards corresponds to a positive tilt reading on the X-axis tilt meter, so the telescope ends up pointing at a higher elevation than expected.
In order to use the tilt meter readings to correct for the variable part of the tilt we first have to convert from the two orthogonal X,Y readings to their polar equivalents, T and T. The following diagram illustrates this.
Applying the four parts and cosine formulae to the highlighted triangle gives the following equations:
By applying, the spherical trigonometry sine, analogue to the cosine, and a variant of the cosine formula to the bottom left triangle, one obtains the following equations.
The CBI will need two separate collimation models.
The position angle of the deck will be denoted by . This increases as the deck rotates clockwise around the direction of the elevation vector and is here defined to point in the direction of increasing elevation when is zero. While tracking a source, the position angle of the deck will be rotated to keep the field of view from rotating. The deck angle will thus be given by,
= o + P.Where o is a user selected deck-angle offset and P is the corrected parallactic angle of the source. The collimation tilt will be defined by a position angle offset, c, and a magnitude c. The following diagram indicates the azimuth and elevation adjustments required to point at a source at the target azimuth and elevation.
The equations that relate the adjusted azimuth and elevation to the target elevation and azimuth are given by applying, respectively, the sine and cosine spherical trigonometry formulae to the green triangle.
Note that exploiting small-angle expectations to approximate cos2c to 1, and sin2c to 0 or c2, could result in significant errors, particularly near the zenith. For example, if +=90o, E=80o and c=0.1o, then the returned elevation will be wrong by almost 2 seconds of arc.
Also note that the apparent parallactic angle is not effected by collimation errors.
Note that terms involving have been replaced
by equations involving
cos of the
corrected A,E and P.
It remains to convert the positions and rates to encoder units. This involves multiplying by a fixed conversion factor and adding measured encoder offsets.
|kf||The constant of flexure (degree per cosine of elevation).|
|HA||The tilt of the azimuth axis in the direction of increasing hour angle.|
|A||The tilt of the azimuth axis in the direction of increasing latitude.|
|E||The tilt of the elevation axis perpendicular to the azimuth ring, measured clockwise around the direction of the azimuth vector. It makes sense to cache the values of sinE and cosE.|
|c||The magnitude of the collimation error. It makes sense to cache the values of sinc and cosc.|
|c||The deck angle at which the collimation error is directed radially outward.|
|Eo||The elevation encoder count at zero elevation.|
|Ao||The azimuth encoder count at zero azimuth.|
|Po||The deck encoder count that aligns the reference edge of the triangular platform to be horizontal. The reference edge should be chosen such that Po is close to zero.|
The following diagram illustrates the relationship between the X,Y offsets and the resulting azimuth and elevation corrections.
The corresponding equations are:
Note that since the offset depends on the azimuth and elevation of the pointing center, which changes as the source moves across the sky, this offset has to be recomputed every servo cycle.