Position angle
Position angle, usually abbreviated PA, is the convention for measuring angles on the sky in astronomy. The International Astronomical Union defines it as the angle measured relative to the north celestial pole (NCP), turning positive into the direction of the right ascension. In the standard (non-flipped) images this is a counterclockwise measure relative to the axis into the direction of positive declination.
In the case of observed visual binary stars, it is defined as the angular offset of the secondary star from the primary relative to the north celestial pole.
As the example illustrates, if one were observing a hypothetical binary star with a PA of 135 degrees, that means an imaginary line in the eyepiece drawn from the north celestial pole to the primary (P) would be offset from the secondary (S) such that the NCP-P-S angle would be 135 degrees.
When graphing visual binary orbits, the NCP line is traditionally drawn downward—that is, with north at bottom—and PA is measured counterclockwise.
Also the direction of the proper motion can for example be given by its position angle.
The definition of position angle is also applied to extended objects like galaxies, where it refers to the angle made by the major axis of the object with the NCP line.
Contents
1 Nautics
1.1 Global geocentric coordinate system
1.2 Brief Derivation
1.3 Long Derivation
2 Astronomy
2.1 Brief Derivation
2.2 Long Derivation
3 See also
4 Further reading
5 References
6 External links
Nautics
The concept of the position angle is inherited from nautical navigation on the oceans, where the optimum compass course is the course from a known position s to a target position t with minimum effort. Setting aside the influence of winds and ocean currents, the optimum course is the course of smallest distance between the two positions on the ocean surface. Computing the compass course is known as the inverse problem of geodesics.
This article considers only the abstraction of minimizing the distance between s and t traveling on the surface of a sphere with some radius R. In which direction angle p relative to North should the ship steer to reach the target position?
Global geocentric coordinate system
Detailed evaluation of the optimum direction is possible if the sea surface is approximated by a sphere surface. The standard computation places the ship at a geodetic latitude φs and geodetic longitude λs, where φ is considered positive if north of the equator, and where λ is considered positive if east of Greenwich. In the global coordinate system centered at the center of the sphere, the Cartesian components are
- s=R(cosφscosλscosφssinλssinφs){displaystyle {mathbf {s} }=Rleft({begin{array}{c}cos varphi _{s}cos lambda _{s}\cos varphi _{s}sin lambda _{s}\sin varphi _{s}end{array}}right)}
and the target position is
- t=R(cosφtcosλtcosφtsinλtsinφt).{displaystyle {mathbf {t} }=Rleft({begin{array}{c}cos varphi _{t}cos lambda _{t}\cos varphi _{t}sin lambda _{t}\sin varphi _{t}end{array}}right).}
The North Pole is at
- N=R(001).{displaystyle {mathbf {N} }=Rleft({begin{array}{c}0\0\1end{array}}right).}
The minimum distance d is the distance along a great circle that runs through s and t. It is calculated in a plane that contains the sphere center and the great circle,
- ds,t=Rθs,t{displaystyle d_{s,t}=Rtheta _{s,t}}
where θ is the angular distance of two points viewed from the center of the sphere, measured in radians. The cosine of the angle is calculated by the dot product of the two vectors
- s⋅t=R2cosθs,t=R2(sinφssinφt+cosφscosφtcos(λt−λs)){displaystyle mathbf {s} cdot mathbf {t} =R^{2}cos theta _{s,t}=R^{2}(sin varphi _{s}sin varphi _{t}+cos varphi _{s}cos varphi _{t}cos(lambda _{t}-lambda _{s}))}
If the ship steers straight to the North Pole, the travel distance is
- ds,N=Rθs,N=R(π/2−φs){displaystyle d_{s,N}=Rtheta _{s,N}=R(pi /2-varphi _{s})}
If a ship starts at t and swims straight to the North Pole, the travel distance is
- dt,N=Rθt,n=R(π/2−φt){displaystyle d_{t,N}=Rtheta _{t,n}=R(pi /2-varphi _{t})}
Brief Derivation
The cosine formula of spherical trigonometry [1] yields for the
angle p between the great circles through s that point to the North on one hand and to t on the other hand
- cosθt,N=cosθs,tcosθs,N+sinθs,tsinθs,Ncosp.{displaystyle cos theta _{t,N}=cos theta _{s,t}cos theta _{s,N}+sin theta _{s,t}sin theta _{s,N}cos p.}
- sinφt=cosθs,tsinφs+sinθs,tcosφscosp.{displaystyle sin varphi _{t}=cos theta _{s,t}sin varphi _{s}+sin theta _{s,t}cos varphi _{s}cos p.}
The sine formula yields
- sinpsinθt,N=sin(λt−λs)sinθs,t.{displaystyle {frac {sin p}{sin theta _{t,N}}}={frac {sin(lambda _{t}-lambda _{s})}{sin theta _{s,t}}}.}
Solving this for sin θs,t and insertion in the previous formula gives an expression for the tangent of the position angle,
- sinφt=cosθs,tsinφs+sin(λt−λs)sinpcosφtcosφscosp;{displaystyle sin varphi _{t}=cos theta _{s,t}sin varphi _{s}+{frac {sin(lambda _{t}-lambda _{s})}{sin p}}cos varphi _{t}cos varphi _{s}cos p;}
- tanp=sin(λt−λs)cosφtcosφssinφt−cosθs,tsinφs.{displaystyle tan p={frac {sin(lambda _{t}-lambda _{s})cos varphi _{t}cos varphi _{s}}{sin varphi _{t}-cos theta _{s,t}sin varphi _{s}}}.}
Long Derivation
Because the brief derivation gives an angle between 0 and π which does not reveal the sign (west or east of north ?), a more explicit derivation is desirable which yields separately the sine and the cosine of p such that use of the correct branch of the inverse tangent allows to produce an angle in the full range -π≤p≤π.
The computation starts from a construction of the great circle between s and t. It lies in the plane that contains the sphere center, s and t and is constructed rotating s by the angle θs,t around an axis ω. The axis is perpendicular to the plane of the great circle and computed by the normalized vector cross product of the two positions:
- ω=1R2sinθs,ts×t=1sinθs,t(cosφssinλssinφt−sinφscosφtsinλtsinφscosλtcosφt−cosφssinφtcosλscosφscosφtsin(λt−λs)).{displaystyle mathbf {omega } ={frac {1}{R^{2}sin theta _{s,t}}}mathbf {s} times mathbf {t} ={frac {1}{sin theta _{s,t}}}left({begin{array}{c}cos varphi _{s}sin lambda _{s}sin varphi _{t}-sin varphi _{s}cos varphi _{t}sin lambda _{t}\sin varphi _{s}cos lambda _{t}cos varphi _{t}-cos varphi _{s}sin varphi _{t}cos lambda _{s}\cos varphi _{s}cos varphi _{t}sin(lambda _{t}-lambda _{s})end{array}}right).}
A right-handed tilted coordinate system with the center at the center of the sphere is given by the
following three axes: the
axis s, the axis
- s⊥=ω×1Rs=1sinθs,t(cosφtcosλt(sin2φs+cos2φssin2λs)−cosλs(sinφscosφssinφt+cos2φssinλscosφtsinλt)cosφtsinλt(sin2φs+cos2φscos2λs)−sinλs(sinφscosφssinφt+cos2φscosλscosφtcosλt)cosφs[cosφssinφt−sinφscosφtcos(λt−λs)]){displaystyle mathbf {s} _{perp }=omega times {frac {1}{R}}mathbf {s} ={frac {1}{sin theta _{s,t}}}left({begin{array}{c}cos varphi _{t}cos lambda _{t}(sin ^{2}varphi _{s}+cos ^{2}varphi _{s}sin ^{2}lambda _{s})-cos lambda _{s}(sin varphi _{s}cos varphi _{s}sin varphi _{t}+cos ^{2}varphi _{s}sin lambda _{s}cos varphi _{t}sin lambda _{t})\cos varphi _{t}sin lambda _{t}(sin ^{2}varphi _{s}+cos ^{2}varphi _{s}cos ^{2}lambda _{s})-sin lambda _{s}(sin varphi _{s}cos varphi _{s}sin varphi _{t}+cos ^{2}varphi _{s}cos lambda _{s}cos varphi _{t}cos lambda _{t})\cos varphi _{s}[cos varphi _{s}sin varphi _{t}-sin varphi _{s}cos varphi _{t}cos(lambda _{t}-lambda _{s})]end{array}}right)}
and the axis ω.
A position along the great circle is
- s(θ)=cosθs+sinθs⊥,0≤θ≤2π.{displaystyle mathbf {s} (theta )=cos theta mathbf {s} +sin theta mathbf {s} _{perp },quad 0leq theta leq 2pi .}
The compass direction is given by inserting the two vectors s and s⊥ and computing the gradient of the vector with respect to θ at θ=0.
- ∂∂θs∣θ=0=s⊥.{displaystyle {frac {partial }{partial theta }}mathbf {s} _{mid theta =0}=mathbf {s} _{perp }.}
The angle p is given by splitting this direction along two orthogonal directions in the plane tangential to the sphere at the point s. The two directions are given by the partial derivatives of s with respect to φ and with respect to λ, normalized to unit length:
- uN=(−sinφscosλs−sinφssinλscosφs);{displaystyle mathbf {u} _{N}=left({begin{array}{c}-sin varphi _{s}cos lambda _{s}\-sin varphi _{s}sin lambda _{s}\cos varphi _{s}end{array}}right);}
- uE=(−sinλscosλs0);{displaystyle mathbf {u} _{E}=left({begin{array}{c}-sin lambda _{s}\cos lambda _{s}\0end{array}}right);}
- uN⋅s=uE⋅uN=0{displaystyle mathbf {u} _{N}cdot mathbf {s} =mathbf {u} _{E}cdot mathbf {u} _{N}=0}
uN points north and uE points east at the position s.
The position angle p projects s⊥
into these two directions,
s⊥=cospuN+sinpuE{displaystyle mathbf {s} _{perp }=cos p,mathbf {u} _{N}+sin p,mathbf {u} _{E}},
where the positive sign means the positive position angles are defined to be north over east. The values of the cosine and sine of p are computed by multiplying this equation on both sides with the two unit vectors,
- cosp=s⊥⋅uN=1sinθs,t[cosφssinφt−sinφscosφtcos(λt−λs)];{displaystyle cos p=mathbf {s} _{perp }cdot mathbf {u} _{N}={frac {1}{sin theta _{s,t}}}[cos varphi _{s}sin varphi _{t}-sin varphi _{s}cos varphi _{t}cos(lambda _{t}-lambda _{s})];}
- sinp=s⊥⋅uE=1sinθs,t[cosφtsin(λt−λs)].{displaystyle sin p=mathbf {s} _{perp }cdot mathbf {u} _{E}={frac {1}{sin theta _{s,t}}}[cos varphi _{t}sin(lambda _{t}-lambda _{s})].}
Instead of inserting the convoluted expression of s⊥, the evaluation may employ that the triple product is invariant under a circular shift
of the arguments:
- cosp=(ω×1Rs)⋅uN=ω⋅(1Rs×uN).{displaystyle cos p=(mathbf {omega } times {frac {1}{R}}mathbf {s} )cdot mathbf {u} _{N}=omega cdot ({frac {1}{R}}mathbf {s} times mathbf {u} _{N}).}
If atan2 is used to compute the value, one can reduce both expressions by division through cos φt
and multiplication by sin θs,t,
because these values are always positive and that operation does not change signs; then effectively
- tanp=sin(λt−λs)cosφstanφt−sinφscos(λt−λs).{displaystyle tan p={frac {sin(lambda _{t}-lambda _{s})}{cos varphi _{s}tan varphi _{t}-sin varphi _{s}cos(lambda _{t}-lambda _{s})}}.}
Astronomy
The parallactic angle of astronomy is the position angle of the zenith (considered the target position) as seen from the star (considered the starting position), where the North Celestial Pole (NCP) of the celestial sphere defines the reference direction.
The sign of the angle is basically kept, north over east in both cases,
but as astronomers look at stars from the inside of the celestial sphere,
the definition uses the convention that the q is the angle in an image that turns the direction to the NCP counterclockwise into the direction of the zenith.
Brief Derivation
The formulas are special cases of the formulas shown above. In the equatorial system of right ascension α and declination δ of a star
one considers the spherical triangle with an angle q at the position of the star,
with an angle h, the hour angle, at the NCP, and with an angle π-A at the
zenith, where A is the azimuth in the south-over-west convention. The side lengths of the spherical triangle are (i) π/2-φ as the distance from the zenith to the NCP, where φ is the observer's geographic latitude, (ii) π/2-δ as the distance between the star and the NCP, and (iii) z=π/2-a, the zenith distance, complement to the star elevation, as the distance between the zenith and the star.
The cosine formula (one of the three cosine formulas) is
- cos(π/2−φ)=cos(π/2−δ)cosz+sin(π/2−δ)sinzcosq.{displaystyle cos(pi /2-varphi )=cos(pi /2-delta )cos z+sin(pi /2-delta )sin zcos q.}
The sine formula is
- sin|q|sin(π/2−φ)=sin|h|sinz=sin(π−A)sin(π/2−δ).{displaystyle {frac {sin |q|}{sin(pi /2-varphi )}}={frac {sin |h|}{sin z}}={frac {sin(pi -A)}{sin(pi /2-delta )}}.}
The 4-parts formula is
- cosq=−coshcos(π−A)+sin|h|sin(π−A)cos(π/2−φ).{displaystyle cos q=-cos hcos(pi -A)+sin |h|sin(pi -A)cos(pi /2-varphi ).}
As in the brief derivation above, these formulas hide the sign of q and h, which
we indicated by placing the absolute values in arguments of the sines, which are sign-dependent.
Long Derivation
The explicit vector algebra to recover the correct sign is equivalent to the calculation of the
long derivation for the compass course. In the global equatorial system the star is at
- s=(cosδcosαcosδsinαsinδ).{displaystyle mathbf {s} =left({begin{array}{c}cos delta cos alpha \cos delta sin alpha \sin delta end{array}}right).}
In the same coordinate system the zenith is found by inserting a=π/2, cos a=0
into the transformation formulas
- z=(cosφcoslcosφsinlsinφ),{displaystyle mathbf {z} =left({begin{array}{c}cos varphi cos l\cos varphi sin l\sin varphi end{array}}right),}
where φ is the observer's geographic latitude, a the star's altitude,
z=π/2-a the zenith distance, and l the local sidereal time. The North Celestial Pole is at
- N=(001).{displaystyle mathbf {N} =left({begin{array}{c}0\0\1end{array}}right).}
The normalized cross product is the rotation axis that turns the star into the direction of the zenith:
- ωz=1sinzs×z=1sinz(cosδsinαsinφ−sinδcosφsinl−cosδcosαsinφ+sinδcosφcoslcosδcosφsin(α−l)).{displaystyle mathbf {omega } _{z}={frac {1}{sin z}}mathbf {s} times mathbf {z} ={frac {1}{sin z}}left({begin{array}{c}cos delta sin alpha sin varphi -sin delta cos varphi sin l\-cos delta cos alpha sin varphi +sin delta cos varphi cos l\cos delta cos varphi sin(alpha -l)end{array}}right).}
Finally ωz X s is the third axis of the tilted coordinate system and the direction into which the star is moved on the great circle towards the zenith.
The plane tangential to the celestial sphere at the star is spanned by the unit vectors to the north,
- uδ=(−sinδcosα−sinδsinαcosδ),{displaystyle mathbf {u} _{delta }=left({begin{array}{c}-sin delta cos alpha \-sin delta sin alpha \cos delta end{array}}right),}
and to the east
- uα=(−sinαcosα0).{displaystyle mathbf {u} _{alpha }=left({begin{array}{c}-sin alpha \cos alpha \0end{array}}right).}
These are orthogonal:
- uδ⋅uα=0;uδ2=uα2=1.{displaystyle mathbf {u} _{delta }cdot mathbf {u} _{alpha }=0;quad mathbf {u} _{delta }^{2}=mathbf {u} _{alpha }^{2}=1.}
The parallactic angle q is the angle of the initial section of the great circle
at s, east of north,
- ωz×s=cosquδ+sinquα.{displaystyle omega _{z}times mathbf {s} =cos q,mathbf {u} _{delta }+sin q,mathbf {u} _{alpha }.}
- cosq=(ωz×s)⋅uδ=1sinz(cosδsinφ−sinδcosφcosh),{displaystyle cos q=(omega _{z}times mathbf {s} )cdot mathbf {u} _{delta }={frac {1}{sin z}}(cos delta sin varphi -sin delta cos varphi cos h),}
- sinq=(ωz×s)⋅uα=1sinzsinhcosφ.{displaystyle sin q=(omega _{z}times mathbf {s} )cdot mathbf {u} _{alpha }={frac {1}{sin z}}sin hcos varphi .}
The values of sin z and of cos φ are positive, so using atan2 functions one may
divide both expressions through these without losing signs; eventually
- tanq=sinhcosφcosδsinφ−sinδcosφcosh=sinhcosδtanφ−sinδcosh{displaystyle tan q={frac {sin hcos varphi }{cos delta sin varphi -sin delta cos varphi cos h}}={frac {sin h}{cos delta tan varphi -sin delta cos h}}}
yields the angle in the full range -π ≤ q ≤ π. The advantage of this expression is that it
does not depend on the various offset conventions of A; the uncontroversial offset
of the hour angle h
takes care of this.
For a sidereal target, by definition a target where δ and α are not time-dependent,
the angle changes with a period of a sidereal day Ts.
Let dots denote time derivatives; then the hour angle changes as
- h˙=2πTs{displaystyle {dot {h}}={frac {2pi }{T_{s}}}}
and the time derivative of the tan q expression is
- q˙1cos2q=cosφ[coshcosδsinφ−sinδcosφ](cosδsinφ−sinδcosφcosh)2h˙;{displaystyle {dot {q}}{frac {1}{cos ^{2}q}}={frac {cos varphi [cos hcos delta sin varphi -sin delta cos varphi ]}{(cos delta sin varphi -sin delta cos varphi cos h)^{2}}}{dot {h}};}
- q˙=cosφ[coshcosδsinφ−sinδcosφ]sin2zh˙=cosφcosacosAsin2zh˙=cosφcosAsinzh˙.{displaystyle {dot {q}}={frac {cos varphi [cos hcos delta sin varphi -sin delta cos varphi ]}{sin ^{2}z}}{dot {h}}={frac {cos varphi cos acos A}{sin ^{2}z}}{dot {h}}={frac {cos varphi cos A}{sin z}}{dot {h}}.}
See also
- Angular distance
Further reading
D. Scott Birney; Guillermo Gonzalez; David Oesper (2007). Observational Astronomy. Cambridge University Press. p. 75. ISBN 0-521-85370-2..mw-parser-output cite.citation{font-style:inherit}.mw-parser-output q{quotes:"""""""'""'"}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-lock-limited a,.mw-parser-output .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}
Taff, Laurence G. (1981). Computational spherical astronomy. Wiley. Bibcode:1981csa..book.....T. ISBN 0471-873179.
Karttunen, Hannu; Kröger, P.; Oja, H.; Poutanen, Markku; Donner, Karl Johan (1987). Fundamental Astronomy. Springer. Bibcode:2003fuas.book.....K. ISBN 0-387-17264-5.
References
^ Abramowitz, Milton; Stegun, Irene Ann, eds. (1983) [June 1964]. "Chapter 4.3.149". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications. ISBN 978-0-486-61272-0. LCCN 64-60036. MR 0167642. LCCN 65-12253.
External links
- The Orbits of 150 Visual Binary Stars, by Dibon Smith (Accessed 2/26/06)