[Previous] [Next] [Up] [Top] [Edit]
1 slalink - programs available with the slalink extensions to tawk.
DESCRIPTION
The inclusion of several of the Star Link astrometry library functions in <#tawk#> allows the easy creation of scripts and programs that make complex astrometric transformations of the columns of starbase data tables.
PROGRAMS
*<#precess#> Precess the ra and dec columns from one equinox to another. *<#fk45#> Convert the ra and dec columns of an starbase data table from FK4 to FK5. *<#fk54#> Convert the ra and dec columns of an starbase data table from FK5 to FK4. *<#radecbox#> Search an starbase data table catalog for targets in an ra dec specified box. *<#radeccirc#> Search an starbase data table catalog for targets in an ra dec specified circle.
SEE ALSO
*<#slalink.4#> functions available with the slalink extensions to tawk. *<#slalib.3#> The Star Link Astrometry library. *<#starbase.1#> starbase data tables introduction (this page). *<#tawk#> awk language extensions implimented at SAO. *<#mawk#> an implimentation of the awk programming language by
Mike Brennan (brennanboeing.com).
4 slalink - functions available with the slalink extensions to tawk.
SYNOPSYS
preces(sys, ep0, ep1, ra, dec) convert from equinox0 to equinox1. fk45(r1950, d1950, bepoch) convert from FK4 to FK5. fk54(r2000, d2000, bepoch) convert from FK5 to FK4. r2d(radians) convert radians to degrees. d2r(degrees) convert degrees to radians. r2h(radians) convert radians to hours. h2r(hours) convert hours to radians.
DESCRIPTION
preces(sys, ep0, ep1, ra, dec) convert from equinox0 to equinox1.
Precession - either FK4 (Bessel-Newcomb, pre-IAU1976) or FK5 (Fricke, post-IAU1976) as required.
*sys precession to be applied: "FK4" or "FK5" *ep0,ep1 starting and ending epoch *ra,dec RA,Dec, mean equator & equinox of epoch ep0
RETURNS
*[ ra,dec ] RA,Dec, mean equator & equinox of epoch ep1
Notes
*1) The epochs are Besselian if sys='FK4' and Julian if 'FK5'.
For example, to precess coordinates in the old system from equinox 1900.0 to 1950.0 the call would be. #!
[ ra, dec ] = preces ("FK4", 1900.0, 1950.0, ra, dec) #
*2) This routine will not correctly convert between the old and
the new systems - for example conversion from B1950 to J2000. For these purposes see fk425, fk524, fk45z and fk54.
*3) If an invalid sys is supplied, values of -99.0,-99.0 will
be returned for both ra and dc.
*P.T.Wallace Starlink 22 December 1993
fk45(r1950, d1950, bepoch) convert from FK4 to FK5.
Convert B1950.0 FK4 star data to J2000.0 FK5 assuming zero proper motion in an inertial frame (double precision)
This routine converts stars from the old, Bessel-Newcomb, FK4 system to the new, IAU 1976, FK5, Fricke system, in such a way that the FK5 proper motion is zero. Because such a star has, in general, a non-zero proper motion in the FK4 system, the routine requires the epoch at which the position in the FK4 system was determined.
*r1950,d1950 B1950.0 FK4 RA,Dec at epoch (rad) *bepoch Besselian epoch (e.g. 1979.3)
RETURNS
*[ r2000, d2000 ] J2000.0 FK5 RA,Dec (rad)
Notes
*1) The epoch BEPOCH is strictly speaking Besselian, but
if a Julian epoch is supplied the result will be affected only to a negligible extent.
*2) Conversion from Besselian epoch 1950.0 to Julian epoch
2000.0 only is provided for. Conversions involving other epochs will require use of the appropriate precession, proper motion, and E-terms routines before and/or after FK425 is called.
*3) In the FK4 catalogue the proper motions of stars within
10 degrees of the poles do not embody the differential E-term effect and should, strictly speaking, be handled in a different manner from stars outside these regions. However, given the general lack of homogeneity of the star data available for routine astrometry, the difficulties of handling positions that may have been determined from astrometric fields spanning the polar and non-polar regions, the likelihood that the differential E-terms effect was not taken into account when allowing for proper motion in past astrometry, and the undesirability of a discontinuity in the algorithm, the decision has been made in this routine to include the effect of differential E-terms on the proper motions for all stars, whether polar or not. At epoch 2000, and measuring on the sky rather than in terms of dRA, the errors resulting from this simplification are less than 1 milliarcsecond in position and 1 milliarcsecond per century in proper motion.
Example:
#!
[ r2000, d2000 ] = fk45(r1950, d1950, 1950) #
fk54(r2000, d2000, bepoch) convert from FK5 to FK4.
Convert a J2000.0 FK5 star position to B1950.0 FK4 assuming zero proper motion and parallax.
This routine converts star positions from the new, IAU 1976, FK5, Fricke system to the old, Bessel-Newcomb, FK4 system.
*r2000,d2000 J2000.0 FK5 RA,Dec (rad) *bepoch Besselian epoch (e.g. 1950)
RETURNS
*[ r1950, d1950 ] B1950.0 FK4 RA,Dec (rad) at epoch BEPOCH
Notes
*1) The proper motion in RA is dRA/dt rather than cos(Dec)*dRA/dt.
*2) Conversion from Julian epoch 2000.0 to Besselian epoch 1950.0
only is provided for. Conversions involving other epochs will require use of the appropriate precession routines before and after this routine is called.
*3) Unlike in the slaFK524 routine, the FK5 proper motions, the
parallax and the radial velocity are presumed zero.
*4) It is the intention that FK5 should be a close approximation
to an inertial frame, so that distant objects have zero proper motion; such objects have (in general) non-zero proper motion in FK4, and this routine returns those fictitious proper motions.
*5) The position returned by this routine is in the B1950
reference frame but at Besselian epoch bepoch. For comparison with catalogues the bepoch argument will frequently be 1950.0.
*P.T.Wallace Starlink 30 October 1993
Example:
#!
[ r1950, d1950 ] = fk54(r2000, d2000, 1950) #
SEE ALSO
*<#slalink.1#> programs available with the slalink extensions to tawk. *<#slalib.3#> The Star Link Astrometry library. *<#tabtable.1#> starbase data tables introduction (this page). *<#tablefun.3#> starbase data tables introduction (this page). *<#tawk#> awk language extensions implimented at SAO. *<#mawk#> an implimentation of the awk programming language by
Mike Brennan (brennanboeing.com).