Import standard modules:


In [1]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from IPython.display import HTML 
HTML('../style/course.css') #apply general CSS


Out[1]:

In [2]:
from IPython.display import HTML
HTML('../style/code_toggle.html')


Out[2]:
The raw code for this notebook is by default hidden for easier reading. To toggle on/off the raw code, click here.

4.1 The Baseline and Its Representation in Space


Outline

4.1.1 The physical baseline

  • 4.1.1.A Local definition of the baseline vector

  • 4.1.1.B The (East, North, Up) reference frame: a global definition of the baseline vector

  • 4.1.1.C Expression in Horizontal Coordinates ($\mathcal{A}$, $\mathcal{E}$)

  • 4.1.1.D Expression in the Equatorial (X,Y,Z) Frame

4.1.2 The ($u$,$v$,$w$) Space


The baseline, which is the separation vector between two antenna element in an interferometric array, is a fundamental concept in interferometry. An array consists of several baselines (formed by every pair of antennas in the array), and thus the baseline configuration is determined by the configuration of antennas in the array. The concept of projected baseline will be introduced in this chapter, and expressions will be derived for the components of the baseline vector in various coordinate systems.

We will start by defining the notations along with the vector definition we will use in the entire course. This progression is designed to go from the local and instantaneous space of the observer (local $x$,$y$, $z$ on Earth) to the sky coordinates in the equatorial frame ($H$, $\delta$) by a series of coordinate transforms:

  1. From local ($x$,$y$,$z$) to the (East, North, Up)
  2. From local E, N, U to local azimuth/elevation ($\mathcal{A}$, $\mathcal{E}$)
  3. From local ($\mathcal{A}$, $\mathcal{E}$) to equatorial (($H$, $\delta$), or equivalently, ($X$, $Y$, $Z$))
  4. From the equatorial ($X$, $Y$, $Z$) to the ($u$, $v$, $w$)-space

4.1.1 The physical baseline

The physical baseline is a construct of the geometry which separates two antenna elements in a 3-D space, while the projected baseline is a mapping of the 3-D physical baseline to a 2-D plane determined by the direction of observation. While the physical baseline stays fixed in a terrestrial reference frame (ITRF is a common frame definition) the projected baseline changes as the Earth rotates relative to the point in the sky which is being observed.

4.1.1.A Local definition of the baseline vector

A position in 3-D can be described by a vector and defining an origin $O$ and orthogonal basis $\mathcal{B}$:

  • A local Cartesian reference frame attached to the local ground defined as $\mathcal{R}$ ($O$, $\mathcal{B}$($\mathbf{\hat{e}_x}, \mathbf{\hat{e}_y}, \mathbf{\hat{e}_z}$)), as illustrated in Fig. 4.1.1 ⤵.

  • Two antennas defined by their position vector in basis $\mathcal{B}$:

$$\textbf{a}_1=\vec{OA}_1=x_1 \mathbf{\hat{e}_x} + y_1 \mathbf{\hat{e}_y} + z_1 \mathbf{\hat{e}_z}$$$$\textbf{a}_2=\vec{OA}_2=x_2\mathbf{\hat{e}_x} + y_2 \mathbf{\hat{e}_y} + z_2 \mathbf{\hat{e}_z}$$

and defined by their components ($x_1$,$y_1$,$z_1$) and ($x_2$,$y_2$,$z_2$) in basis $\mathcal{B}$.

We assume that the distance between the antennas is small enough to consider the ground to be a ground plane on which the basis vectors can be translated without any rotation. The local physical baseline in Cartesian coordinates is defined as the difference vector

$$\mathbf{b} = \mathbf{a}_2 - \mathbf{a}_1 = \vec{A_1A_2} \qquad.$$

It is important to note that the physical baseline is a quantity which depends on the relative positions of two antenna elements and is independent of the origin of the reference frame.

Figure 4.1.1: Vectors $a_1$ and $a_2$ in a local 3-dimensional Cartesian coordinate system.

4.1.1.B The (East, North, Up) Reference Frame: A Global Definition of the Baseline Vector

The problem with the previous formulation is that directions of $\mathbf{\hat{e}_x}, \mathbf{\hat{e}_y}, \mathbf{\hat{e}_z}$ are arbitrary. Instead the coordinate basis can be attached to an absolute reference in relation to the Earth. This can be achieved by defining a new Cartesian basis defined by the cardinal points and use a terrestrial reference frame. To do this we define a new origin and basis

  • We define the origin $O$ at some point on Earth. This origin can be arbitrary, but it is convenient to chose the 'centre' of the array layout or at the position of one of the array elements.
  • A global right-handed orthonormal basis $\mathcal{B}'$($\mathbf{\hat{e}_E}, \mathbf{\hat{e}_N}, \mathbf{\hat{e}_U}$) pointing respectively towards the geographical East and North, and towards the local vertical (pointing to the Zenith away from the Earth, e.g. Up).

The baseline vector is still expressed in a local reference frame but is coincidental with the geographical North and East (see Fig. 4.1.2 ⤵).

4.1.1.C Expression in Horizontal Coordinates ($\mathcal{A}$, $\mathcal{E}$)

As seen in $\S$ 3.4 ➞, an observer located somewhere on Earth can define a direction in the sky in terms of local azimuth $\mathcal{A}$ and elevation $\mathcal{E}$. Similarly, a baseline can also be expressed with these two coordinates. We now define:

  • Antenna 1 as the origin of the reference frame on the ground.
  • The North, East and Up axes (see $\S$ 4.1.1.B)

The azimuth is the angle contained in the plane of the local ground measured clock-wise from North to East. The elevation is the angle measured from the horizon to the local zenith.

Figure 4.1.2: Relation between the horizontal frame ($\mathcal{A}$, $\mathcal{E}$) and the ($E$, $N$, $U$) Cartesian frame.

The baseline vector is expressed in the $\mathcal{B}'$ basis:

\begin{equation} \mathbf{b}_{\text{ENU}} = \lvert \mathbf{b} \rvert \begin{bmatrix} \sin \mathcal{A} \cos \mathcal{E}\\ \cos \mathcal{A} \cos \mathcal{E}\\ \sin \mathcal{E} \end{bmatrix} \end{equation}
Note: If Antenna 2 and Antenna 1 share the same component along $\mathbf{\hat{e}_U}$ (e.g. Antenna 1 and Antenna 2 lying on the ground plane), then the elevation of the baseline is $\mathcal{E}=0^\circ$.
Warnings:
  • Geographers and the Astronomers use different conventions for the origin and variation of the azimuth. We will use the Astronomer convention as defined previously in [Chapter 3 ➞](../3_Positional_Astronomy/3_0_Introduction.ipynb).
  • The horizontal coordinate system looks like a 2D spherical coordinates system with different conventions in terms of angle origin and direction of variation.
  • The latitude $L_a$ was not expressed here but will be required later in the chapter.

4.1.1.D Expression in the Equatorial (X,Y,Z) Frame

The baseline is fully described in the $ENU$ system using the azimuth $\mathcal{A}$ and elevation $\mathcal{E}$. As shown in $\S$ 3.2 ➞, at the observer location at latitude $L_a$, the extension of the direction of the baseline defines a position on the sky, which can be associated with Equatorial coordinates ($H$, $\delta$), with $H$ being the hour angle and $\delta$ the declination.

To generalize the baseline further (and to ease the subsequent derivations), we need to define a set of reference frames which will associate the baseline definition to sky coordinates on the celestial sphere. For that purpose, we will define a new set of terrestrial reference frame $XYZ$ (defined by the basis ($\mathbf{\hat{e}_X}$,$\mathbf{\hat{e}_Y}$,$\mathbf{\hat{e}_Z}$)), attached to the Earth, as an intermediate reference frame. We can position these axes with respect to the Equatorial coordinates $(H; \delta)$. We defined the axes of the $XYZ$ frame as:

  • the $X$-axis points towards $(H=0^\textrm{h}, \delta = 0^{\circ})$ ($\gamma$ the point where the vernal equinox crosses the local meridian)
  • the $Y$-axis towards $(H=-6^\textrm{h}, \delta = 0^{\circ})$ (East)
  • the $Z$-axis towards the North Celestial Pole (NCP).

Fig. 4.1.3 ⤵ shows the plane of the array elements (red) can be related to the the plane of the celestial sphere via a coordinate transformation.

Figure 4.1.3: The (X,Y,Z) Reference frame. $T$: local zenith. ($X$,$Y$) = Celestial plane (blue) $\perp$ to the NCP. ($N$,$S$,$E$,$W$) Antenna plane (red) $\perp$ to the local zenith (T).

Warning:
  • Although we use Equatorial coordinates to define the $XYZ$-coordinate system, it represents a terrestrial reference frame, as the $X$, $Y$ and $Z$ axes never change orientation for a local observer on Earth.
  • However, as the declination refers to the celestial equator, the latitude of the observer $L_a$ is required.

We can convert ($\mathcal{A}$, $\mathcal{E}$, $L_a$) into $(H, \delta)$ in this new frame:

\begin{equation} \begin{bmatrix} X\\ Y\\ Z \end{bmatrix} = \begin{bmatrix} \lvert \mathbf{b} \rvert \cos \delta \cos H\\ -\lvert \mathbf{b} \rvert \cos \delta \sin H\\ \lvert \mathbf{b} \rvert \sin \delta \end{bmatrix} = \lvert \mathbf{b} \rvert \begin{bmatrix} \cos L_a \sin \mathcal{E} - \sin L_a \cos \mathcal{E} \cos \mathcal{A}\nonumber\\ \cos E \sin \mathcal{A} \nonumber\\ \sin L_a \sin \mathcal{E} + \cos L_a \cos \mathcal{E} \cos \mathcal{A}\\ \end{bmatrix} \end{equation}

Equation A: Conversion from baseline vector $\mathcal{E},\mathcal{A}$ to $XYZ$ . $\mathbf{b}$: amplitude of baseline vector. $H$: Hour angle, $\delta$: Declination, $L_a$ latitude of the array.

Note: The complete derivation of Eq. A from the spherical triangle ZTS can be found in [Appendix A.1 ➞](../0_Introduction/2_Appendix.ipynb)

4.1.2 The ($u$,$v$,$w$) Space

Now that we have defined a baseline in a terrestrial XYZ frame we are ready for the final transformation to the celestial uvw reference frame. Let ($H_0$, $\delta_0$) be the point on the celestial sphere in the direction of of a source $\mathbf{s_0}$. We now define a new set of axes:

  • the $u$-axis lies in the celestial equatorial plane, and points toward the hour angle $H_0-6^\text{h}$.
  • the $v$-axis lies in the plane of the great circle with hour angle $H_0$, and points toward the declination $\frac{\pi}{2}-\delta_0$.
  • the $w$-axis points in the direction of $\mathbf{s_0}$.

The ($u$,$v$) axes form a 2-D plane perpendicular to $\mathbf{s_0}$.

The transformation from the ($X$,$Y$,$Z$) frame to the ($u$, $v$, $w$) frame requires two succesive rotations applied to the equatorial coordinates of the baseline:

  1. $\mathcal{R}_\text{Z}^{\frac{\pi}{2}-H_0}$: the rotation of the $X$,$Y$ axes around the $Z=Z'$ axis of angle $\frac{\pi}{2}-H_0$ making the new axes $X'$, $Y'$ (Green axes in Fig. 4.1.4 ⤵)
  2. $\mathcal{R}_\text{X'}^{\frac{\pi}{2}-\delta_0}$: the rotation of the $Y'$,$Z'$ axes around the $X'=u$ axis of angle $\frac{\pi}{2}-\delta_0$ making the new axes $u$, $v$ (red axes in Fig. 4.1.4 ⤵).

Figure 4.1.4: Relation of the (X,Y,Z) frame to the ($u$,$v$,$w$) frame. $Z$: local zenith. ($X$,$Y$) = Celestial plane $\perp$ to the NCP.

The conversion from $XYZ$ to $uvw$ is expressed as following: \begin{equation} \begin{pmatrix} u\ v\ w

\end{pmatrix}

\mathcal{R}_\text{X'}^{\frac{\pi}{2}-\delta0} \mathcal{R}\text{Z}^{\frac{\pi}{2}-H_0} \begin{pmatrix} X\\ Y\\ Z \end{pmatrix} \end{equation}

where \begin{equation}

\mathcal{R}_\text{X'}^{\frac{\pi}{2}-\delta_0}

\begin{pmatrix} 1& 0& 0\\ 0 & \cos (\frac{\pi}{2}- \delta_0) & \sin (\frac{\pi}{2}- \delta_0) \\ 0 &-\sin (\frac{\pi}{2}- \delta_0) & \cos (\frac{\pi}{2}- \delta_0) \\ \end{pmatrix}

\begin{pmatrix} 1& 0& 0\\ 0 & \sin \delta_0 & \cos \delta_0 \\ 0 &-\cos \delta_0 & \sin \delta_0 \\ \end{pmatrix}

\end{equation}

and \begin{equation}

\mathcal{R}_\text{Z}^{\frac{\pi}{2}-H_0}

\begin{pmatrix} \cos (\frac{\pi}{2}- H_0) & \sin (\frac{\pi}{2}- H_0) & 0 \\ -\sin (\frac{\pi}{2}- H_0) & \cos (\frac{\pi}{2}- H_0) & 0 \\ 0 & 0 & 1 \end{pmatrix}

\begin{pmatrix} \sin H_0 & \cos H_0 & 0 \\ -\cos H_0 & \sin H_0 & 0 \\ 0 & 0 & 1 \end{pmatrix}

\end{equation}

Warning: The sequence of operations applied to vectors should be read from right to left as it is the rule in linear algebra.

Resulting in the transformation matrix: \begin{equation} \begin{pmatrix} u\ v\ w

\end{pmatrix}

\begin{pmatrix} \sin H_0 & \cos H_0 & 0\\ -\sin \delta_0 \cos H_0 & \sin\delta_0\sin H_0 & \cos\delta_0\\ \cos \delta_0 \cos H_0 & -\cos\delta_0\sin H_0 & \sin\delta_0\\ \end{pmatrix}

\begin{pmatrix} X\\ Y\\ Z \end{pmatrix} \end{equation}

Warning: We usually express the components ($u$, $v$, $w$) in units of wavelength (of the observation) by dividing the baseline vector by $\lambda$:
\begin{equation} \mathbf{b_\lambda} = \frac{\mathbf{b}}{\lambda} = \begin{pmatrix} u\\ v\\ w \end{pmatrix} = \frac{1}{\lambda} \begin{pmatrix} \sin H_0 & \cos H_0 & 0\\ -\sin \delta_0 \cos H_0 & \sin\delta_0\sin H_0 & \cos\delta_0\\ \cos \delta_0 \cos H_0 & -\cos\delta_0\sin H_0 & \sin\delta_0\\ \end{pmatrix} \begin{pmatrix} X\\ Y\\ Z \end{pmatrix} \end{equation} In the following sections and chapters, we will drop the $\lambda$ underscript, as we will always assume that the $u$, $v$, $w$ components are in wavelength unit.

Conclusion

This section was dedicated to the expression of the baseline vector is astronomy-friendly reference frames. Later on, we will use the results of this section to quickly link the physical length of the baseline between two receivers, to a quantity measured on the celestial sphere. In $\S$ 4.2 ➞, we will focus on a simple one dimension 2-element interferometer. We will address the correlation of two signals collected from a remote point source and construct the corresponding visibility function which will be addressed with more detailed in $\S$ 4.3 ➞