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]:
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
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:
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.
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}$:
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.
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
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 ⤵).
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:
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}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:
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.
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.
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$,$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:
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.
where \begin{equation}
\end{equation}
and \begin{equation}
\end{equation}
Resulting in the transformation matrix: \begin{equation} \begin{pmatrix} u\ v\ w
\begin{pmatrix} X\\ Y\\ Z \end{pmatrix} \end{equation}
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 ➞