Transform ECI to ECEF

Question

I need to transform ECI to ECEF coordinate. Here, I need the formula in math and algorithm. I have looking for the formula in the internet, but I can find one. I have x, y, z cartesian coord in ECI. And want to transform it to ECEF coord. Anybody knows the formula and algorithm?

Answer 1

The Explanatory Supplement to the Astronomical Almanac has all the equations you need. Take a look at chapters 3 and 4.

Keep in mind you need some clear definitions of what you mean by ECEF and ECI. Most people utilize WGS84 for ECEF, but that is not a requirement. Similarly ECI could be J2000 or ICRF

In general you will need 4 steps to convert from ECI to ECEF:

1. Calculate earth's precession
2. Calculate earth's nutation
3. Account for earth's rotation including UTC-UT1 offset
4. Account for polar motion
5. Convert to lan/lon if desired

The first two transformations are purely analytical, and can be obtained from books like the Explanatory Supplement. The third and fourth transformation rely on irregularly varying parameters that have to be updated for a specific date.

The help system for the Systems ToolKit (STK) has some excellent descriptions on various frames and transformation process, but not actual equations.

Answer 2

The ECI and ECEF frames have approximately the same origin and $z$ axis and only differ by an angular component on the $xy$ plane. This angular difference can be used to rotate a vector from one frame to the other.

The $x_{ECEF}$ axis rotates along the equatorial plane about $z_{ECI}=z_{ECEF}$. The angle that it makes with the $x_{ECI}$ axis is known as the Earth Rotation Angle (ERA) or Greenwich Sidereal Angle. To convert between ECI and ECEF frames, we must first gain the ERA denoted as $\gamma\in[0, 2\pi]$.

We begin at an epoch, say J2000, where the ERA is known to be $280.46$ degrees. Using the elapsed time since epoch, we can propagate the epoch ERA through time to a date of interest. Stated mathematically in the following block equation, the Earth rotation rate is given by $360.985...$ degrees per day, $\Delta T$ is the elapsed time in days since the epoch, and $280.46$ is the ERA at $\Delta T=0$.

\begin{equation}\label{eqERA} \gamma = 360.9856123035484\times\Delta T + 280.46 \end{equation}

To convert a vector between the ECI and ECEF frames requires rotating the vector about the $z$-axis by the Greenwich Sidereal Angle, $\gamma$, using the rotational matrix given in the following block equation. These conversions can be seen in equations the final two equations where $\mu_{C'}^C$ converts coordinates from frame $C$ to frame $C'$.

\begin{equation}\label{eqEz} E_3^\theta = \begin{pmatrix} \cos\theta & -\sin\theta & 0 \\ \sin\theta & \cos\theta & 0 \\ 0 & 0 & 1 \\ \end{pmatrix} \end{equation}

\begin{equation}\label{eqECEF2ECI} \vec{V}_{ECI} = \mu_{ECI}^{ECEF}\;\vec{V}_{ECEF} = E_3^\gamma\;\vec{V}_{ECEF} \end{equation} \begin{equation}\label{eqECI2ECEF} \vec{V}_{ECEF} = \mu_{ECEF}^{ECI}\;\vec{V}_{ECI} = E_3^{-\gamma}\;\vec{V}_{ECI} \end{equation}