13.07.2015 Views

Using Rotations to Build Aerospace Coordinate Systems - Defence ...

Using Rotations to Build Aerospace Coordinate Systems - Defence ...

Using Rotations to Build Aerospace Coordinate Systems - Defence ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

DSTO–TN–0640ZNDEN 0E 0D 0 N 1E 1ωαYXD 1Figure 5: Constructing local north–east–down axes N, E, D at a given poin<strong>to</strong>f known latitude and longitude. Start with the three vec<strong>to</strong>rs N 0 , E 0 , D 0 .Construct the intermediate set N 1 , E 1 , D 1 by rotating the original setthrough the longitude ω about N 0 . Then rotate N 1 , E 1 , D 1 about −E 1 by thelatitude α, according <strong>to</strong> the right hand rule, <strong>to</strong> give the final set N, E, D. Orequivalently, rotate the intermediate set about E 1 through minus the latitude.The resulting vec<strong>to</strong>rs are the local north–east–down axesthe act of rotating each of the vec<strong>to</strong>rs that represent the axes of some frame.There are two basic tasks <strong>to</strong> consider in this type of aerospace calculation. First, givena place on Earth (usually in lat–long–height coordinates), we need <strong>to</strong> construct the localgeographic frame’s axes, such as NED. Second, an aircraft can be introduced in<strong>to</strong> thisNED frame, and its orientation found relative <strong>to</strong> that frame.4.1 Constructing NED Axes for the Local Geographic FrameGiven the lat–long–height of a place on or near Earth, the first thing we wish <strong>to</strong> do is findthe local NED axes. The technique used is the primary one of this report. We build aninitial set of NED axes in a place where it’s simple <strong>to</strong> determine them, and then we rotatethese around <strong>to</strong> the required place.The simplest place <strong>to</strong> construct an initial set of NED axes is on the junction of theEqua<strong>to</strong>r and the prime meridian, since this has a latitude and longitude of α = ω = 0 ◦ .Represent each axis by a unit vec<strong>to</strong>r in the ECEF frame, which all calculations are beingdone within:⎡ ⎤0⎡ ⎤0⎡ ⎤−1N 0 = ⎣0⎦ , E 0 = ⎣1⎦ , D 0 = ⎣ 0⎦ . (4.1)1 0 0The scenario is shown in Figure 5. Now rotate each of the N 0 , E 0 , D 0 vec<strong>to</strong>rs in two steps.12

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!