3_Positional_Astronomy/3_4_direction_cosine_coordinates.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "***\n",
    "\n",
    "* [Outline](../0_Introduction/0_introduction.ipynb)\n",
    "* [Glossary](../0_Introduction/1_glossary.ipynb)\n",
    "* [Positional Astronomy](3_0_introduction.ipynb)\n",
    "    * Previous: [Horizontal Coordinates](3_3_horizontal_coordinates.ipynb)\n",
    "    * Next: [Further Reading](3_x_further_reading_and_references.ipynb) \n",
    "\n",
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Import standard modules:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       "\n",
       "p.prerequisites {\n",
       "    border:5px;\n",
       "    border-style:none;\n",
       "    background-color:#F0F0F0;\n",
       "    box-shadow: 5px 5px 2px #888888;\n",
       "    border-color:#FF0000;\n",
       "    padding: 1em;\n",
       "    }\n",
       "\n",
       "p.conclusion {\n",
       "    border:5px;\n",
       "    border-style:none;\n",
       "    background-color:#c8ebca;\n",
       "    box-shadow: 5px 5px 2px #888888;\n",
       "    border-color:#c8ebca;\n",
       "    padding: 1em;\n",
       "    }\n",
       "\n",
       "div.warn {    \n",
       "    background-color: #fcf2f2;\n",
       "    border-color: #dFb5b4;\n",
       "    border-left: 5px solid #dfb5b4;\n",
       "    padding: 0.5em;\n",
       "    }\n",
       "\n",
       "div.advice {    \n",
       "    background-color: #ecf8ed;\n",
       "    border-color: #a2dda5;\n",
       "    border-left: 5px solid #a2dda5;\n",
       "    padding: 0.5em;\n",
       "    }\n",
       "\n",
       "</style>\n",
       "\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 71,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "from IPython.display import HTML \n",
    "HTML('../style/course.css') #apply general CSS"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Import section specific modules:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<script>\n",
       "    code_show=false; \n",
       "    function code_toggle() {\n",
       "        if (code_show){\n",
       "            $('div.input').hide();\n",
       "        } else {\n",
       "            $('div.input').show();\n",
       "        }\n",
       "        code_show = !code_show\n",
       "    } \n",
       "    $( document ).ready(code_toggle);\n",
       "</script>\n",
       "\n",
       "The raw code for this notebook is by default hidden for easier reading. To toggle on/off the raw code, click <a href=\"javascript:code_toggle()\">here</a>.\n",
       "\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import HTML\n",
    "HTML('../style/code_toggle.html')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "## Direction Cosine Coordinates"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "There is another useful astronomical coordinate system that we ought to introduce at this juncture, namely the *direction cosine coordinate system*. The direction cosine coordinate system is quite powerful and allows us to redefine the fundamental reference point on the celestial sphere, from which we measure all other celestial objects, to an arbitrary location (i.e. we can make local sky-maps around our own chosen reference point; the vernal equinox need not be our fundamental reference point). Usually this arbitrary location is chosen to be the celestial source that we are interested in observing. We generally refer to this arbitrary location as the *field centre* or *phase centre*.\n",
    "\n",
    "<div class=advice>\n",
    "<b>Note:</b> The direction cosine coordinate system is useful for another reason, when we use\n",
    "it to image interferometric data, then it becomes evident that there exists a Fourier relationship between the sky brightness function and the measurements that an interferometer makes (see <a href='../4_Visibility_Space/4_0_introduction.ipynb'>Chapter 4 &#10142;</a>).\n",
    "</div><br>\n",
    "We use three coordinates in the direction cosine coordinate system, namely $l$, $m$ and $n$. The coordinates $l$, $m$ and $n$ are dimensionless direction cosines, i.e."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "\\begin{eqnarray}\n",
    "l &=& \\cos(\\alpha) = \\frac{a_1}{|\\mathbf{a}|}\\\\\n",
    "m &=& \\cos(\\beta) = \\frac{a_2}{|\\mathbf{a}|}\\\\\n",
    "n &=& \\cos(\\gamma) = \\frac{a_3}{|\\mathbf{a}|}\n",
    "\\end{eqnarray}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<a id='pos:fig:cosines'></a> <!--\\label{pos:fig:cosines}--><img src='figures/cosine.svg' width=35%>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "*Figure 3.4.1*: Definition of direction cosines. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "The quantities $\\alpha$, $\\beta$, $\\gamma$, $a_1$, $a_2$, $a_3$ and $\\mathbf{a}$ are all defined in <a class='pos_fig_cos_dir'></a> <!--\\ref{pos:fig:cos}-->. Moreover, $|\\cdot|$ denotes the magnitude of its operand. The definitions above also imply that $l^2+m^2+n^2 = 1$. When $|\\mathbf{a}|=1$ then we may simply interpret $l$, $m$ and $n$ as   Cartesian coordinates, i.e. we may simply relabel the axes $x$, $y$ and $z$ (in <a class='pos_fig_cos_dir'></a><!--\\ref{pos:fig:cos}-->) to \n",
    "$l$, $m$ and $n$. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "So the question now arises, how do we use $l$, $m$ and $n$ to uniquely identify a location on the celestial sphere? The direction cosine coordinate system (and the relationship between it and the celestial coordinate sytem) is depicted in <a class='pos_fig_dirconversion_dir'></a><!--\\ref{pos:fig:dirconversion}-->. Note that the $n$-axis points toward the field center (which is denoted by $\\boldsymbol{s}_c$ in <a class='pos_fig_dirconversion_dir'></a><!--\\ref{pos:fig:dirconversion}-->. It should be clear from <a class='pos_fig_dirconversion_dir'></a><!--\\ref{pos:fig:dirconversion}--> that we can use $\\mathbf{s} = (l,m,n)$ to uniquely idnetify any location on the celestial sphere."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<a id='pos:fig:convert_lmn_ra_dec'></a> <!--\\label{pos:fig:convert_lmn_ra_dec}--><img src='figures/conversion2.svg' width=40%>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "*Figure 3.4.2*: The source-celestial pole-field center triangle; which enables us to derive the conversion equations between direction cosine and equatorial coordinates. The red plane represents the fundamental plane of the equatorial coordinate system, while the blue plane represents the fundamental plane of the direction cosine coordinate system. We are able to label the orthogonal fundamental axes of the direction cosine coordinate system $l$,$m$ and $n$, since the radius of the celestial sphere is equal to one."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "We use the following equations to convert between the equatorial and direction cosine coordinate systems: "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<p class=conclusion>\n",
    "  <font size=4><b>Converting between the equatorial and direction cosine coordinates (3.1)</b></font>\n",
    "  <br>\n",
    "  <br>\n",
    "\\begin{eqnarray}\n",
    "l &=&  \\sin \\theta \\sin \\psi = \\cos \\delta  \\sin \\Delta \\alpha \\nonumber\\\\\n",
    "m &=& \\sin \\theta \\cos \\psi = \\sin \\delta \\cos \\delta_0 - \\cos \\delta \\sin \\delta_0 \\cos\\Delta \\alpha \\nonumber\\\\\n",
    "\\delta &=& \\sin^{-1}(m\\cos \\delta_0 + \\sin \\delta_0\\sqrt{1-l^2-m^2})\\nonumber\\\\\n",
    "\\alpha &=& \\alpha_0 + \\tan^{-1}\\bigg(\\frac{l}{\\cos\\delta_0\\sqrt{1-l^2-m^2}-m\\sin\\delta_0}\\bigg)\\nonumber\n",
    "\\end{eqnarray}\n",
    "</p>\n",
    "**<a id='pos_eq_convertlmnradec'></a><!--\\label{pos:eq:convertlmnradec}-->**"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<div class=advice>\n",
    "<b>Note:</b> See <a href='../0_Introduction/2_Appendix.ipynb'>Appendix &#10142;</a> for the derivation of the above relations.\n",
    "</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "We can obtain the conversion relations above by applying the spherical trigonemetric identities in <a href='../2_Mathematical_Groundwork/2_13_spherical_trigonometry.ipynb'>$\\S$ 2.13 &#10142;</a> to the triangle depicted in <a class='pos_fig_dirconversion_dir'></a><!--\\ref{pos:fig:dirconversion}--> (the one formed by the source the field center and the NCP)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "There is another important interpretation of direction cosine coordinates we should\n",
    "be cognisant of. If we project the direction cosine position vector $\\mathbf{s}$ of a celestial body onto the $lm$-plane it's projected length will be equal to $\\sin \\theta$, where $\\theta$ is the angular distance between your field center $\\mathbf{s}_c$ and $\\mathbf{s}$ measured along the surface of the celestial sphere. If $\\theta$ is small we may use the small angle approximation, i.e. $\\sin \\theta \\approx \\theta$. The projected length of $\\mathbf{s}$ is also equal to $\\sqrt{l^2+m^2}$, implying that $l^2+m^2 \\approx \\theta^2$. We may therefore loosely interpret $\\sqrt{l^2+m^2}$ as the angular distance measured between the source at $\\mathbf{s}$ and the field-center $\\mathbf{s}_c$ measured along the surface of the celestial sphere, i.e. we may measure $l$ and $m$ in $^{\\circ}$. \n",
    "The explenation above is graphically illustrated in <a class='pos_fig_proj_dir'></a> <!--\\ref{pos:fig:proj}-->. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<a id='pos:fig:understand_lm'></a> <!--\\label{pos:fig:understand_lm}--><img src='figures/conversion2b.svg' width=40%>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "*Figure 3.4.3*: Why do we measure $l$ and $m$ in degrees? "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<p class=conclusion>\n",
    "  <font size=4><b>Three interpretations of direction cosine coordinates</b></font>\n",
    "  <br>\n",
    "  <br>\n",
    "&bull; **Direction cosines**: $l$,$m$ and $n$ are direction cosines<br><br>\n",
    "&bull; **Cartesian coordinates**: $l$,$m$ and $n$ are Cartesian coordinates if we work on the \n",
    "    unit sphere<br><br> \n",
    "&bull; <b>Angular distance</b>: $\\sqrt{l^2+m^2}$ denotes the angular distance $\\theta$, $(l,m,n)$ is from the field center (if $\\theta$ is sufficiently small).\n",
    "</p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Example"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Here we have a couple of sources given in RA ($\\alpha$) and DEC ($\\delta$):\n",
    "* Source 1: (5h 32m 0.4s,60$^{\\circ}$17' 57'') - 1Jy\n",
    "* Source 2: (5h 36m 12.8s,61$^{\\circ}$ 12' 6.9'') - 0.5Jy\n",
    "* Source 3: (5h 40m 45.5s,61$^{\\circ}$ 56' 34'') - 0.2Jy\n",
    "\n",
    "The field center is located at $(\\alpha_0,\\delta_0) = $ (5h 30m,60$^{\\circ}$). The first step is to convert right ascension and declination into radians with"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "\\begin{eqnarray}\n",
    "\\alpha_{\\textrm{rad}} &=& \\frac{\\pi}{12} \\bigg(h + \\frac{m}{60} + \\frac{s}{3600}\\bigg)\\\\\n",
    "\\delta_{\\textrm{rad}} &=& \\frac{\\pi}{180} \\bigg(d + \\frac{m_{\\textrm{arcmin}}}{60}+\\frac{s_{\\textrm{arcsec}}}{3600}\\bigg)\n",
    "\\end{eqnarray}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "In the above equations $h,~m,~s,~d,~m_{\\textrm{arcmin}}$ and $s_{\\textrm{arcsec}}$ respectively denote hours, minutes, seconds, degrees, arcminutes and arcseconds. If we apply the above to our three sources we obtain"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "RA (rad) of Sources and Field Center =  [ 1.43989663  1.44865237  1.46700741  1.48683872]\n",
      "DEC (rad) of Sources =  [ 1.04719755  1.05241899  1.06817495  1.08110542]\n"
     ]
    }
   ],
   "source": [
    "RA_rad = (np.pi/12) * np.array([5. + 30./60, 5 + 32./60 + 0.4/3600, 5 + 36./60 + 12.8/3600, 5 + 40./60 + 45.5/3600])\n",
    "DEC_rad = (np.pi/180)*np.array([60., 60. + 17.0/60 + 57./3600, 61. + 12./60 + 6.9/3600, 61 + 56./60 + 34./3600])\n",
    "Flux_sources_labels = np.array([\"\", \"1 Jy\", \"0.5 Jy\", \"0.2 Jy\"])\n",
    "Flux_sources = np.array([1., 0.5, 0.1]) #in Janskys\n",
    "print \"RA (rad) of Sources and Field Center = \", RA_rad\n",
    "print \"DEC (rad) of Sources = \", DEC_rad"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Recall that we can use <a class='pos_eq_convertlmnradec_dir'></a><!--\\label{pos:eq:convertlmnradec}--> to convert between equatorial and direction cosine coordinates, in terms of the current example this translates into the python code below. Note that before we can do the conversion we first need to calculate $\\Delta \\alpha$."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "l (degrees) =  [ 0.          0.24855826  0.74818685  1.26458942]\n",
      "m (degrees) =  [ 0.          0.30010768  1.21061225  1.96811494]\n"
     ]
    }
   ],
   "source": [
    "RA_delta_rad = RA_rad-RA_rad[0] #calculating delta alpha\n",
    "\n",
    "l = np.cos(DEC_rad) * np.sin(RA_delta_rad)\n",
    "m = (np.sin(DEC_rad) * np.cos(DEC_rad[0]) - np.cos(DEC_rad) * np.sin(DEC_rad[0]) * np.cos(RA_delta_rad))\n",
    "print \"l (degrees) = \", l*(180./np.pi)\n",
    "print \"m (degrees) = \", m*(180./np.pi)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Plotting the result."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/trienko/NASSP/fundamentals/lib/python2.7/site-packages/ipykernel/__main__.py:8: MatplotlibDeprecationWarning: pyplot.hold is deprecated.\n",
      "    Future behavior will be consistent with the long-time default:\n",
      "    plot commands add elements without first clearing the\n",
      "    Axes and/or Figure.\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEOCAYAAACTqoDjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGv1JREFUeJzt3X+UVPWd5vH3I6gEccjRaZiEHyGKzYCKKGA07qitKIqK\nGjRKOhJQT4cc43HOJrrjMiHGyOqOWdfJkIxhNDAZCNh7DMqg+ANt4rCOUVCC4q+OoyjEjebsjrGH\nqAE++0cVnQZud1XB7b63up7XOXW6v7e+dfvpVvrp+6PuVURgZma2pwOyDmBmZvnkgjAzs0QuCDMz\nS+SCMDOzRC4IMzNL5IIwM7NEuSsISX0kPS9pZdZZzMxqWe4KArgOeDnrEGZmtS5XBSFpKHAecHfW\nWczMal3frAPs4U7gBuDQziZIagKaAPr16zd++PDhPRRt3+3cuZMDDshVFydyzvRUQ0ZwzrRVS87X\nXnvttxFRV3JiROTiAZwP/LD4+enAylKvqa+vj2rQ0tKSdYSyOGd6qiFjhHOmrVpyAuuijN/Leaq6\nU4Cpkt4ElgFnSFqcbSQzs9qVm4KIiBsjYmhEjAAuB56IiC9nHMvMrGblpiDMzCxf8naQGoCIWAOs\nyTiGmVlN8xaEmZklckGYmVkiF4SZmSVyQZiZWSIXhJmZJXJBmJlZIheEmZklckGYmVkiF4SZmSVy\nQZiZWSIXhJmZJXJBmJlZIheEmZklckGYmVkiF4SZmSVyQZiZWSIXhJmZJXJBmJlZotwUhKR+kp6R\n9EtJmyR9J+tMZma1LE/3pP4IOCMi2iQdCKyVtCoins46mJlZLcpNQUREAG3F4YHFR2SXyMystuVm\nFxOApD6SNgDvAo9FxC+yzmRmVqtU+MM9XyR9ElgOXBsRL+7xXBPQBFBXVze+ubk5g4SVaWtrY8CA\nAVnHKMk501MNGcE501YtORsaGtZHxISSEyMilw9gLvDNrubU19dHNWhpack6QlmcMz3VkDHCOdNW\nLTmBdVHG7+Hc7GKSVFfcckDSJ4CzgFeyTWVmVrtyc5Aa+BTwj5L6UDg20hwRKzPOZGZWs3JTEBGx\nETg+6xxmZlaQm11MZmaWLy4IMzNL5IIwM7NELggzM0vkgjCrMQ8//DCjRo1i5MiR/PSnP02cc8cd\ndzBmzBjGjh3LmWeeyebNmxPnVcObwmzfuSDMasiOHTu45pprWLVqFS+99BKPP/44L7300l7zjj/+\neNatW8fGjRu55JJLuOGGGzJIa1lzQZjVkGeeeYaRI0dyxBFHcNBBB3HGGWfwwAMP7DWvoaGB/v37\nA3DSSSexZcuWLtc7Y8YM7r///vZxY2Nj4nqturggzGrI1q1bGTZsWPu4rq6OrVu3dvmae+65h3PP\nPbfLOVdddRWLFi0C4P333+epp57ivPPO2++8li0XhJl1avHixaxbt47rr7++y3mnnXYara2tvPfe\neyxdupRp06bRt29u3odr+8j/Bc1qyJAhQ3j77bfbx++99x5HHHFE4tzVq1czb948fv7zn3PwwQeX\nXPeMGTNYvHgxy5YtY+HChalltux4C8KshkycOJHW1lbeeOMNPv74Y5544gmmTp2617znn3+er371\nq6xYsYJBgwaVte6ZM2dy5513AjBmzJhUc1s2vAVhVkP69u3L/PnzmTx5Mjt27KChoYGjjz4agLlz\n5zJhwgSmTp3K9ddfT1tbG5deeikAw4cPZ8WKFbuta/v27bttWQwePJjRo0dz0UUX9dw3ZN3KBWFW\nY6ZMmcKUKVMAWLNmTfvym2++uf3z1atXl1zPpk2bOPLII9vH27Zto7W1lenTp6cX1jLlXUxmVrG7\n7rqL6dOnc8sttwCFQhk9ejTXXnstAwcOzDidpcVbEGZWsdmzZzN79uz28aRJkzp9t7VVL29BmJlZ\nIheEmZklckGYmVkiF4SZmSXKTUFIGiapRdJLkjZJui7rTGZmtSw3BQFsB74REWOAk4BrJPntmGZZ\nW7IERoyAAw4ofFyyJOtE1kNyc5prRLwDvFP8/ANJLwNDgL0vVm9mPWPJEmhqgm3bCuPNmwtjgMbG\n7HJZj1BEZJ1hL5JGAE8Cx0TE7/Z4rgloAqirqxvf3Nzc4/kq1dbWVhV33nLO9FRDRiid86TLL6ff\nb36z1/IPBw/m6WXLujPabnrLzzMvGhoa1kfEhJITIyJXD2AAsB74Qqm59fX1UQ1aWlqyjlAW50xP\nNWSMKCOnFAF7P6QeybdLr/l55gSwLsr4fZynYxBIOhC4D1gSET/LOo9ZzRs+vLLl1qvkpiAkCbgH\neDki7sg6j5kB8+ZB8daj7fr3Lyy3Xi83BQGcAlwBnCFpQ/ExJetQZll4+OGHGTVqFCNHjuS2225L\nnLNo0SLq6uoYN24c48aN4+67706ct1/7xBsbYcEC+MxnQCp8XLDAB6hrRJ7OYloLKOscZlnbsWMH\n11xzDY899hhDhw5l4sSJTJ06NfEmPJdddhnz58/v3kCNjS6EGpWnLQgzA5555hlGjhzJEUccwUEH\nHcTll1/OAw88sN/rnTFjBvfff3/7uLGxkbVr1+73eq33ckGY5czWrVsZNmxY+3jo0KFs3bo1ce59\n993H2LFjueSSS3a713SSq666ikWLFgHw/vvv89RTT3HyySenltt6HxeEWZW64IILePPNN9m4cSNn\nnXUWX/nKV7qcf9ppp9Ha2sp7773H0qVLmTZtGn369OmhtFaNXBBmOTNkyJDdtga2bNnCkCFD9pp3\n+OGHt98T+uqrr2b9+vUl1z1jxgwWL17MwoULufLKK9MLbb2SC8IsZyZOnEhraytvvPEGH3/8McuW\nLWPq1Kl7zXvnnXfaP1+xYgWjR48uue6ZM2dy5513AiQe9DbrKDdnMZlZQd++fZk/fz6TJ09mx44d\nXHnllRx99NEAzJ07lwkTJjB16lS+//3vs2LFCvr27cthhx3Wfnyho+3bt7dvZQAMHjyY0aNHc9FF\nF/XUt2NVzAVhlkNTpkxhypS93wZ08803t39+6623cuutt3a5nk2bNnHkkUe2j7dt20ZrayvTp09P\nL6z1Wt7FZNZL3XXXXUyfPp1bbrkFgNWrVzN69GiuvfZaBg4cmHE6qwbegjDrpWbPns3s2bPbx5Mm\nTWLz5s0ZJrJq4y0IMzNL5IIwM7NELggzM0vkgjDrTXz/aEuRD1Kb9Ra+f7SlzFsQZr3FnDl/LIdd\ntm0rLDfbBy4Is97irbcqW25WggvCrLfw/aMtZS4Is97C94+2lOWqICT9WNK7kl7MOotZ1fH9oy1l\nJc9iknRYGevZGRH/nkKeRcB84CcprMus9vj+0Zaick5z/XXxoS7m9AH2e0dnRDwpacT+rsfMzPaf\nIqLrCdLzEXH8/s4pO1ChIFZGxDGdPN8ENAHU1dWNb25uTuPLdqu2tjYGDBiQdYySnDM91ZARnDNt\n1ZKzoaFhfURMKDkxIrp8AP3SmFPuAxgBvFjO3Pr6+qgGLS0tWUcoi3OmpxoyRjhn2qolJ7Auyvgd\nW/IgdUR8CCDpUkmHFj//lqSfSTqh4xwzS3bllVcyaNAgjjkmccMYgJtuuonvfe97PZjKrGuVnMX0\nrYj4QNJ/As4E7gH+vntimfUuM2fO5OGHH846hllFKimIHcWP5wELIuJB4KA0w0haCvwrMErSFklX\npbl+s6yceuqpHHZYOScEwuuvv84JJ5zQPm5tbd1tbNZTKimIrZJ+BFwGPCTp4ApfX1JETI+IT0XE\ngRExNCLuSXP9ZtXgyCOPZODAgWzYsAGAhQsXMmvWrIxTWS2q5Bf8F4FHgMlReM/DYcD13ZLKrMZd\nffXVLFy4kB07dnDvvffypS99KetIVoPKLoiI2BYRP4uI1uL4nYh4tPuimdWuadOmsWrVKlauXMn4\n8eM5/PDDs45kNajsglDBlyXNLY6HSzqx+6KZ1a5+/foxefJkvva1r3n3kmWmkl1MPwROBqYXxx8A\nP0g9kVkvNH36dE4++WReffVVhg4dyj337H14bfv27Rx88MHt48bGRg444ADOPvvsnoxq1q6SO8p9\nLiJOkPQ8QET8P0mpnsVk1lstXbq05JxNmzbx+c9/vn28du1aZs2aRZ8+fbozmlmnKimIP0jqAwSA\npDpgZ7ekMqsxxx57LPX19e1bCxdffDGvv/46TzzxRMbJrJZVUhDfB5YDgyXNAy4B/rpbUpnVmBde\neGG38fLlyzNKYvZHZRdERCyRtJ7Cu6gFXBQRL3dbMjMzy1RFZzEBE4DDI2I+8B8+i8nMrPfyWUxm\nZpaokoL4XERcA3wIhbOYSPlaTGY1ZckSGDECDjig8HHJkqwTme3GZzGZZWHJEmhqgm3bCuPNmwtj\n8C1DLTcq2YLYdRbToOJZTGuB/9Ytqcx6uzlz/lgOu2zbVlhulhNlbUEUD1A/CfgsJrM0vPVWZcvN\nMlBWQURESHooIo4FXunmTGa93/Dhhd1KScvNcqKSXUzPSZrYbUnMasm8edC//+7L+vcvLDfLiYrO\nYgKelvS6pI2SXpC0sbuCmfVqjY2wYAF85jMgFT4uWOAD1JYrlZzFNLnbUpjVosZGF4LlWiVbENMS\nHpMkjUsrjKRzJL0q6VeS/iqt9Zr1hL/5G2hp2X1ZS0thuVk1qqQgJgCzgSHFx1eBc4B/kHTD/gYp\nvsfiB8C5wBhguqQx+7tes54ycSJ88Yt/LImWlsJ4oo/cWZWqZBfTUOCEiGgDkPRt4EHgVAqnv+7v\n30knAr+KiH8rrn8ZcCHw0n6u16xHNDRAc3OhFM49dwSrVhXGDQ1ZJzPbN5UUxCDgow7jPwCDI+L3\nkj7q5DWVGAK83WG8hcKB8d1IagKaAOrq6lizZk0KX7p7tbW1OWeK8pxTKpTDP/3TCK644k2kN8lp\nVCDfP8uOnDMblRTEEuAXkh4oji8AfirpEHrwr/yIWAAsABg1alScfvrpPfWl99maNWtwzvTkOWdL\nC6xaBVdc8SarVo1g1qwRud6CyPPPsiPnzEYl94P4rqRVwCnFRbMjYl3x8zROxdgKDOswHlpcZlYV\ndh1zaG4G6U1mzRrRPs5zSZh1ptL7Qfw5MDAi/hZ4N+X7QTwLHCXps8V7XV8OrEhx/Wbd6tlndy+D\nXccknn0221xm+6qSXUw/pHD11jOAmyncD+I+IJVzNCJiu6SvA48AfYAfR8SmNNZt1hNuSDiXr6HB\nWw9WvSopiM9FxAmSnofC/SCKf+mnJiIeAh5Kc51mZrZvKnkfhO8HYWZWQ3w/CDMzS1TJWUxLJPl+\nEGZmNaKSYxBExCv4fhBmZjWhZEFI+s9dPR8Rd6QXx8zM8qKcLYhDix9HUTilddd7Ey4AnumOUGZm\nlr2SBRER3wGQ9CSFi/V9UBzfROFifWZm1gtVchbTYODjDuOPi8vMzKwXquQg9U+AZyQtL44vAhal\nnsjMzHKhktNc5xUv1vcXxUWzIuL57ollZmZZq/Q01+eA57opi5mZ5UjJYxCSShZCOXPMzKy6lLMF\nMVrSxi6eFzAwpTxmZpYT5RTEn5cxZ8f+BjEzs3wp530Qm3siiJmZ5Usl74MwM7MaUlFBSDqwu4KY\nmVm+lH2aq6R/AM6XtB34NbAR2BgRf9dd4czMLDuVbEGcCgyNiGHAFyjcPOiQNEJIulTSJkk7JU1I\nY51mZrZ/Knmj3C+Aw4F3I2IrsJX07h/9IoXS+VFK6zMzs/1USUH8CPi5pHsolMXGiHg/jRC77kwn\nKY3VmZlZChQR5U2U3gAWAAEcB4wF+kXEkamFkdYA34yIdV3MaQKaAOrq6sY3Nzen9eW7TVtbGwMG\nDMg6RknOmZ5qyAjOmbZqydnQ0LA+Ikrvzo+Ish7AvyQsO7iC16+msCtpz8eFHeasASaUu876+vqo\nBi0tLVlHKItzpqcaMkY4Z9qqJSewLsr4HVvJLqYNkq6LiL/tUC4flfviiJhUwdcyM7OMVVIQg4FJ\nkv4LhSu6/hLYEBH/q1uSmZlZpso+zTUivhgRo4HPAnOB14AT0wgh6WJJW4CTgQclPZLGes3MbN9V\ndD8IaN+tlOp9ISJiOYX3VZiZWU74WkxmZpbIBWFmZolcEGZmlsgFYWZmiVwQZmaWyAVhZmaJXBBm\nZpbIBWFmZolcEGZmlsgFYWZmiVwQZmaWyAVhZmaJXBBmZpbIBWFmZolcEGZmlsgFYWZmiVwQZmaW\nyAVhZmaJclEQkm6X9IqkjZKWS/pk1pnMzGpdLgoCeAw4JiLGAq8BN2acx8ys5uWiICLi0YjYXhw+\nDQzNMo+ZmYEiIusMu5H0z8C9EbG4k+ebgCaAurq68c3NzT0Zb5+0tbUxYMCArGOU5JzpqYaM4Jxp\nq5acDQ0N6yNiQsmJEdEjD2A18GLC48IOc+YAyykWV6lHfX19VIOWlpasI5TFOdNTDRkjnDNt1ZIT\nWBdl/I7t2z39lFhEk7p6XtJM4HzgzOI3YGZmGeqxguiKpHOAG4DTImJb1nnMzCwnB6mB+cChwGOS\nNki6K+tAZma1LhdbEBExMusMZma2u7xsQZiZWc64IMzMLJELwszMErkgzMwskQvCzMwSuSDMzCyR\nC8LMzBK5IMzMLJELwszMErkgzMwskQvCzMwSuSDMzCyRC8LMzBK5IMzMLJELwszMErkgzMwskQvC\nzMwSuSDMzCxRLgpC0nclbSzej/pRSZ/OOpOZWa3LRUEAt0fE2IgYB6wE5mYdyMys1uWiICLidx2G\nhwCRVRYzMyvom3WAXSTNA2YA7wMNGccxM6t5iuiZP9YlrQb+LOGpORHxQId5NwL9IuLbnaynCWgC\nqKurG9/c3NwdcVPV1tbGgAEDso5RknOmpxoygnOmrVpyNjQ0rI+ICSUnRkSuHsBw4MVy5tbX10c1\naGlpyTpCWZwzPdWQMcI501YtOYF1Ucbv2Fwcg5B0VIfhhcArWWUxM7OCvByDuE3SKGAnsBmYnXEe\nM7Oal4uCiIhpWWcwM7Pd5WIXk5mZ5Y8LwszMErkgzMwskQvCzMwSuSDMzCyRC8LMzBK5IMzMLJEL\nwszMErkgzMwskQvCzMwSuSDMzCyRC8LMzBK5IMzMLJELwszMErkgzMwskQvCzMwSuSDMzCyRC8LM\nzBK5IMzMLFGuCkLSNySFpD/NOouZWa3LTUFIGgacDbyVdRYzM8tRQQD/E7gBiKyDmJkZ9M06AICk\nC4GtEfFLSaXmNgFNxeFHkl7s7nwp+FPgt1mHKINzpqcaMoJzpq1aco4qZ5IieuYPdkmrgT9LeGoO\n8F+BsyPifUlvAhMiouQPWdK6iJiQbtL0OWe6qiFnNWQE50xbb8vZY1sQETEpabmkY4HPAru2HoYC\nz0k6MSL+T0/lMzOz3WW+iykiXgAG7RpXsgVhZmbdJ08HqffFgqwDlMk501UNOashIzhn2npVzh47\nBmFmZtWl2rcgzMysm7ggzMwsUa8piLxfpkPSdyVtlLRB0qOSPp11pj1Jul3SK8WcyyV9MutMSSRd\nKmmTpJ2ScndKoaRzJL0q6VeS/irrPEkk/VjSu3l/H5GkYZJaJL1U/G9+XdaZkkjqJ+kZSb8s5vxO\n1pk6I6mPpOclrSw1t1cURJVcpuP2iBgbEeOAlcDcrAMleAw4JiLGAq8BN2acpzMvAl8Answ6yJ4k\n9QF+AJwLjAGmSxqTbapEi4Bzsg5Rhu3ANyJiDHAScE1Of54fAWdExHHAOOAcSSdlnKkz1wEvlzOx\nVxQEVXCZjoj4XYfhIeQwa0Q8GhHbi8OnKbwnJXci4uWIeDXrHJ04EfhVRPxbRHwMLAMuzDjTXiLi\nSeD/Zp2jlIh4JyKeK37+AYVfbEOyTbW3KGgrDg8sPnL3b1zSUOA84O5y5ld9QXS8TEfWWUqRNE/S\n20Aj+dyC6OhKYFXWIarQEODtDuMt5PAXWjWSNAI4HvhFtkmSFXfdbADeBR6LiDzmvJPCH9M7y5mc\n+RvlylHOZTp6NlGyrnJGxAMRMQeYI+lG4OvAt3s0IKUzFufMobBpv6Qns3VUTk6rHZIGAPcBf7nH\n1nhuRMQOYFzx2N1yScdERG6O8Ug6H3g3ItZLOr2c11RFQVTLZTo6y5lgCfAQGRREqYySZgLnA2dG\nhm+SqeBnmTdbgWEdxkOLy2wfSTqQQjksiYifZZ2nlIj4d0ktFI7x5KYggFOAqZKmAP2AP5G0OCK+\n3NkLqnoXU0S8EBGDImJERIygsDl/Qh6v4STpqA7DC4FXssrSGUnnUNj8nBoR27LOU6WeBY6S9FlJ\nBwGXAysyzlS1VPjL7x7g5Yi4I+s8nZFUt+usP0mfAM4iZ//GI+LGiBha/F15OfBEV+UAVV4QVeY2\nSS9K2khhl1geT9ebDxwKPFY8HfeurAMlkXSxpC3AycCDkh7JOtMuxYP8XwceoXBAtTkiNmWbam+S\nlgL/CoyStEXSVVln6sQpwBXAGcX/JzcU/wLOm08BLcV/389SOAZR8jTSvPOlNszMLJG3IMzMLJEL\nwszMErkgzMwskQvCzMwSuSDMzCyRC8LMzBK5IMzMLJELwmqOpDMlLe4wHiHp98ULrSXNv0nSN3su\nYWmSPlF809jHeb0HilU/F4TVouOAPcvg9eK9OrqNClL5NxcRvy/m/XUa6zNL4oKwWnQc0OXl4SXN\nkfSapLXAqD2e+3Lx7mEbJP2oeJMgJH2reCe5tZKWSvpmcevkVUk/oXDhtmFdvL6z5YdIerB4t7IX\nJV3WHT8Usz25IKwWJW1BtJM0nsLFzMYBU4CJHZ4bDVwGnFL8C34H0ChpIjCtuO5zgY63Qj0K+GFE\nHA307+T1iestvv4c4NcRcVxEHAM8vJ/fv1lZquJy32ZpKV46emBEvNfFtL8Alu+6oq2kjldjPRMY\nDzxbvMT8JyjcIOYw4IGI+BD4UNI/d3jN5oh4usTr/6ST5QAvAP9D0n8HVkbEv+zL925WKReE1ZrR\nlHk/3k4I+MeI2O1+3ZL+sovX/EcZr782aTlARLwm6QQKWzO3SHo8Im7e5+/ArEzexWS1puTxB+BJ\n4KLimUKHAhd0eO5x4BJJgwAkHSbpM8D/Bi6Q1K9497PzO1l3Z6/vbDmSPg1si4jFwO3ACZV/22aV\n8xaE1ZrjKFyvv1MR8ZykeykUybsd50fES5L+Gni0eEbSH4BrIuLp4q6ojcBvKOwWej9h3V29fq/l\nwGbgWOB2STuLy7+2fz8Cs/L4fhBW8ySNoLBv/5j9XM+AiGiT1J/CVkhTRDyXQsSuvuabwISI+G13\nfh2rTd7FZFY4Y2hgZ2+Uq8CC4jqeA+7rznLY9UY54EBgZ3d9Hatt3oIwM7NE3oIwM7NELggzM0vk\ngjAzs0QuCDMzS+SCMDOzRC4IMzNL5IIwM7NE/x+Q+lpiCc+kqQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f334fddbbd0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(111)\n",
    "plt.xlim([-4., 4.])\n",
    "plt.ylim([-4., 4.])\n",
    "plt.xlabel(\"$l$ [degrees]\")\n",
    "plt.ylabel(\"$m$ [degrees]\")\n",
    "plt.plot(l[0], m[0], \"bx\")\n",
    "plt.hold(\"on\")\n",
    "plt.plot(l[1:]*(180/np.pi), m[1:]*(180/np.pi), \"ro\")\n",
    "counter = 1\n",
    "for xy in zip(l[1:]*(180/np.pi)+0.25, m[1:]*(180/np.pi)+0.25):                                              \n",
    "    ax.annotate(Flux_sources_labels[counter], xy=xy, textcoords='offset points',horizontalalignment='right',\n",
    "                verticalalignment='bottom')  \n",
    "    counter = counter + 1\n",
    "\n",
    "plt.grid()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "*Figure 3.4.4*: "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "***\n",
    "\n",
    "Next: [Further Reading](3_x_further_reading_and_references.ipynb)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<div class=warn><b>Future Additions:</b></div>\n",
    "\n",
    "* figure: projection plot similar to white book figure 2-9"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}