DataWrangling_RunwayClassification.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>.container { width:75% !important; }</style>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from IPython.core.display import display, HTML\n",
    "display(HTML(\"<style>.container { width:75% !important; }</style>\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Description:\n",
    "In this notebook, we acquire one year of flight traffic data for EGLL airport. The dataset is located in [Skywise Foundry](https://core.skywise.com/workspace/dataset/preview/ri.foundry.main.dataset.33e61269-87a3-45ae-aaf2-af6dfdc18761/master). The output of this data wrangling is the removal of features not needed for the machine learning, feature engineering of new variables (eg: `hourofday`),removal of duplicates. Some basic Exploratory Data Analytics is also done."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "#Set some numpy print options for displaying numpy arrays to fit maximum width of cell\n",
    "np.set_printoptions(precision=8, edgeitems=30, linewidth=1000,formatter=dict(float=lambda x: \"%.8g\" % x))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Landing_WinAndGust_Cross_Class</th>\n",
       "      <th>Landing_WindAndGust_Front_Class</th>\n",
       "      <th>Landing_WindAndGust_Cross_kt</th>\n",
       "      <th>Landing_WindAndGust_Front_kt</th>\n",
       "      <th>Landing_Wind_Cross_kt</th>\n",
       "      <th>Landing_Wind_Front_kt</th>\n",
       "      <th>Wind_Speed_kt</th>\n",
       "      <th>Wind_And_Gust_kt</th>\n",
       "      <th>Landing_RWY_Mag_Bearing</th>\n",
       "      <th>Landing_Visibilty_Class_RVR</th>\n",
       "      <th>...</th>\n",
       "      <th>sky_level_4_coverage</th>\n",
       "      <th>sky_level_1_alt_ft</th>\n",
       "      <th>sky_level_2_alt_ft</th>\n",
       "      <th>sky_level_3_alt_ft</th>\n",
       "      <th>sky_level_4_alt_ft</th>\n",
       "      <th>present_weather_codes</th>\n",
       "      <th>present_weather_codes_meaning</th>\n",
       "      <th>AreaCode</th>\n",
       "      <th>RunwayTrueTrack</th>\n",
       "      <th>hourofday</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>(a) 5kt or less</td>\n",
       "      <td>(k) 0/+5kt</td>\n",
       "      <td>-0.693590</td>\n",
       "      <td>1.875882</td>\n",
       "      <td>-0.693590</td>\n",
       "      <td>1.875882</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>90.3</td>\n",
       "      <td>VMC</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>EUR</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>(a) 5kt or less</td>\n",
       "      <td>(l) +5/+10kt</td>\n",
       "      <td>3.032626</td>\n",
       "      <td>8.473676</td>\n",
       "      <td>3.032626</td>\n",
       "      <td>8.473676</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>90.3</td>\n",
       "      <td>VMC</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>EUR</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>(a) 5kt or less</td>\n",
       "      <td>(l) +5/+10kt</td>\n",
       "      <td>4.952684</td>\n",
       "      <td>8.687400</td>\n",
       "      <td>4.952684</td>\n",
       "      <td>8.687400</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.3</td>\n",
       "      <td>VMC</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>EUR</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>(b) 5/10kt</td>\n",
       "      <td>(l) +5/+10kt</td>\n",
       "      <td>-9.243315</td>\n",
       "      <td>7.652524</td>\n",
       "      <td>-9.243315</td>\n",
       "      <td>7.652524</td>\n",
       "      <td>12.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>270.4</td>\n",
       "      <td>VMC</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>800.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>EUR</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>17</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>(a) 5kt or less</td>\n",
       "      <td>(n) +15/+20kt</td>\n",
       "      <td>-3.067543</td>\n",
       "      <td>16.720950</td>\n",
       "      <td>-3.067543</td>\n",
       "      <td>16.720950</td>\n",
       "      <td>17.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>270.4</td>\n",
       "      <td>VMC</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>800.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>EUR</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>18</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 59 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "  Landing_WinAndGust_Cross_Class Landing_WindAndGust_Front_Class  \\\n",
       "0                (a) 5kt or less                      (k) 0/+5kt   \n",
       "1                (a) 5kt or less                    (l) +5/+10kt   \n",
       "2                (a) 5kt or less                    (l) +5/+10kt   \n",
       "3                     (b) 5/10kt                    (l) +5/+10kt   \n",
       "4                (a) 5kt or less                   (n) +15/+20kt   \n",
       "\n",
       "   Landing_WindAndGust_Cross_kt  Landing_WindAndGust_Front_kt  \\\n",
       "0                     -0.693590                      1.875882   \n",
       "1                      3.032626                      8.473676   \n",
       "2                      4.952684                      8.687400   \n",
       "3                     -9.243315                      7.652524   \n",
       "4                     -3.067543                     16.720950   \n",
       "\n",
       "   Landing_Wind_Cross_kt  Landing_Wind_Front_kt  Wind_Speed_kt  \\\n",
       "0              -0.693590               1.875882            2.0   \n",
       "1               3.032626               8.473676            9.0   \n",
       "2               4.952684               8.687400           10.0   \n",
       "3              -9.243315               7.652524           12.0   \n",
       "4              -3.067543              16.720950           17.0   \n",
       "\n",
       "   Wind_And_Gust_kt  Landing_RWY_Mag_Bearing Landing_Visibilty_Class_RVR  \\\n",
       "0               2.0                     90.3                         VMC   \n",
       "1               9.0                     90.3                         VMC   \n",
       "2              10.0                     90.3                         VMC   \n",
       "3              12.0                    270.4                         VMC   \n",
       "4              17.0                    270.4                         VMC   \n",
       "\n",
       "     ...    sky_level_4_coverage  sky_level_1_alt_ft  sky_level_2_alt_ft  \\\n",
       "0    ...                     NaN                 NaN                 NaN   \n",
       "1    ...                     NaN                 NaN                 NaN   \n",
       "2    ...                     NaN                 NaN                 NaN   \n",
       "3    ...                     NaN               800.0                 NaN   \n",
       "4    ...                     NaN               800.0                 NaN   \n",
       "\n",
       "   sky_level_3_alt_ft  sky_level_4_alt_ft present_weather_codes  \\\n",
       "0                 NaN                 NaN                   NaN   \n",
       "1                 NaN                 NaN                   NaN   \n",
       "2                 NaN                 NaN                   NaN   \n",
       "3                 NaN                 NaN                   NaN   \n",
       "4                 NaN                 NaN                   NaN   \n",
       "\n",
       "  present_weather_codes_meaning AreaCode RunwayTrueTrack hourofday  \n",
       "0                           NaN      EUR       89.713523         6  \n",
       "1                           NaN      EUR       89.713523        13  \n",
       "2                           NaN      EUR       89.713523        12  \n",
       "3                           NaN      EUR      269.713523        17  \n",
       "4                           NaN      EUR      269.713523        18  \n",
       "\n",
       "[5 rows x 59 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = pd.read_csv('data/2016_EGLL_Testdata_raw.csv',parse_dates = ['Landing_Time'], infer_datetime_format = True)\n",
    "df['hourofday'] = df.Landing_Time.dt.hour\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 221069 entries, 0 to 221068\n",
      "Data columns (total 58 columns):\n",
      "Landing_WinAndGust_Cross_Class     220349 non-null object\n",
      "Landing_WindAndGust_Front_Class    220349 non-null object\n",
      "Landing_WindAndGust_Cross_kt       220349 non-null float64\n",
      "Landing_WindAndGust_Front_kt       220349 non-null float64\n",
      "Landing_Wind_Cross_kt              220349 non-null float64\n",
      "Landing_Wind_Front_kt              220349 non-null float64\n",
      "Wind_Speed_kt                      221069 non-null float64\n",
      "Wind_And_Gust_kt                   221069 non-null float64\n",
      "Landing_RWY_Mag_Bearing            221069 non-null float64\n",
      "Landing_Visibilty_Class_RVR        221069 non-null object\n",
      "Landing_Visibilty_Class            221069 non-null object\n",
      "Landing_Coverage_Alt_ft            28078 non-null float64\n",
      "Landing_RVR_or_Vis_Quality         221069 non-null int64\n",
      "Landing_RVR_or_Vis_m               221069 non-null float64\n",
      "RVR_Mean                           1844 non-null float64\n",
      "METAR                              221069 non-null object\n",
      "METAR_Time                         221069 non-null object\n",
      "Landing_Time                       221069 non-null datetime64[ns]\n",
      "Landing_RWY                        221069 non-null object\n",
      "Landing_Airport                    221069 non-null object\n",
      "Landing_RWY_Quality                221069 non-null int64\n",
      "Dist_To_Threshold_m                221069 non-null float64\n",
      "Mean_DeltaTrack                    221069 non-null float64\n",
      "Nb_Snapshot_InRWY                  221069 non-null int64\n",
      "Mean_ProjDistance                  221069 non-null float64\n",
      "Nb_Pt_in_Last_min                  221069 non-null int64\n",
      "flight_id                          221069 non-null int64\n",
      "aircraft_registration              220928 non-null object\n",
      "equipment                          221067 non-null object\n",
      "msn                                147839 non-null float64\n",
      "manufacturer                       220760 non-null object\n",
      "aircraft_family                    220760 non-null object\n",
      "aircraft_type                      220750 non-null object\n",
      "airline_icao                       220974 non-null object\n",
      "Schedule_Landing_Airport           220726 non-null object\n",
      "ADBS_Landing_Airport               220720 non-null object\n",
      "temperature_c                      221069 non-null float64\n",
      "dew_point_c                        221069 non-null float64\n",
      "relative_humidity                  221069 non-null float64\n",
      "wind_direction_degrees             203342 non-null float64\n",
      "wind_speed_kn                      221069 non-null float64\n",
      "hour_precipitation                 221069 non-null float64\n",
      "pressure_altimiter_in              221069 non-null float64\n",
      "visibility_mi                      221069 non-null float64\n",
      "wind_gust_kn                       4801 non-null float64\n",
      "sky_level_1_coverage               220529 non-null object\n",
      "sky_level_2_coverage               43994 non-null object\n",
      "sky_level_3_coverage               9270 non-null object\n",
      "sky_level_4_coverage               773 non-null object\n",
      "sky_level_1_alt_ft                 141786 non-null float64\n",
      "sky_level_2_alt_ft                 37287 non-null float64\n",
      "sky_level_3_alt_ft                 4523 non-null float64\n",
      "sky_level_4_alt_ft                 37 non-null float64\n",
      "present_weather_codes              33033 non-null object\n",
      "present_weather_codes_meaning      32993 non-null object\n",
      "AreaCode                           221069 non-null object\n",
      "RunwayTrueTrack                    221069 non-null float64\n",
      "hourofday                          221069 non-null int64\n",
      "dtypes: datetime64[ns](1), float64(28), int64(6), object(23)\n",
      "memory usage: 97.8+ MB\n"
     ]
    }
   ],
   "source": [
    "df.dropna(axis=1, how='all',inplace=True)\n",
    "df.info()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "cols_of_interest = ['Landing_Time','hourofday','aircraft_registration','wind_direction_degrees','Wind_Speed_kt','Wind_And_Gust_kt',\n",
    "                    'Landing_RVR_or_Vis_m','temperature_c','dew_point_c','relative_humidity','pressure_altimiter_in',\n",
    "                    'RunwayTrueTrack','Landing_Visibilty_Class_RVR','present_weather_codes','Landing_RWY']\n",
    "df = df[cols_of_interest]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Landing_Time</th>\n",
       "      <th>hourofday</th>\n",
       "      <th>aircraft_registration</th>\n",
       "      <th>wind_direction_degrees</th>\n",
       "      <th>Wind_Speed_kt</th>\n",
       "      <th>Wind_And_Gust_kt</th>\n",
       "      <th>Landing_RVR_or_Vis_m</th>\n",
       "      <th>temperature_c</th>\n",
       "      <th>dew_point_c</th>\n",
       "      <th>relative_humidity</th>\n",
       "      <th>pressure_altimiter_in</th>\n",
       "      <th>RunwayTrueTrack</th>\n",
       "      <th>Landing_Visibilty_Class_RVR</th>\n",
       "      <th>present_weather_codes</th>\n",
       "      <th>Landing_RWY</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2016-01-01 06:15:58</td>\n",
       "      <td>6</td>\n",
       "      <td>GVIIC</td>\n",
       "      <td>70.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>8000.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>93.08</td>\n",
       "      <td>30.21</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2016-01-01 13:19:52</td>\n",
       "      <td>13</td>\n",
       "      <td>TCJJS</td>\n",
       "      <td>110.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>81.20</td>\n",
       "      <td>30.03</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2016-01-01 12:01:11</td>\n",
       "      <td>12</td>\n",
       "      <td>DAIZU</td>\n",
       "      <td>120.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>81.20</td>\n",
       "      <td>30.09</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2016-01-02 17:22:49</td>\n",
       "      <td>17</td>\n",
       "      <td>GEUPX</td>\n",
       "      <td>220.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>7000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>29.38</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2016-01-02 18:20:40</td>\n",
       "      <td>18</td>\n",
       "      <td>ECLXQ</td>\n",
       "      <td>260.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>87.47</td>\n",
       "      <td>29.41</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         Landing_Time  hourofday aircraft_registration  \\\n",
       "0 2016-01-01 06:15:58          6                 GVIIC   \n",
       "1 2016-01-01 13:19:52         13                 TCJJS   \n",
       "2 2016-01-01 12:01:11         12                 DAIZU   \n",
       "3 2016-01-02 17:22:49         17                 GEUPX   \n",
       "4 2016-01-02 18:20:40         18                 ECLXQ   \n",
       "\n",
       "   wind_direction_degrees  Wind_Speed_kt  Wind_And_Gust_kt  \\\n",
       "0                    70.0            2.0               2.0   \n",
       "1                   110.0            9.0               9.0   \n",
       "2                   120.0           10.0              10.0   \n",
       "3                   220.0           12.0              12.0   \n",
       "4                   260.0           17.0              17.0   \n",
       "\n",
       "   Landing_RVR_or_Vis_m  temperature_c  dew_point_c  relative_humidity  \\\n",
       "0                8000.0            2.0          1.0              93.08   \n",
       "1               10000.0            7.0          4.0              81.20   \n",
       "2               10000.0            7.0          4.0              81.20   \n",
       "3                7000.0           11.0         10.0              93.55   \n",
       "4               10000.0           11.0          9.0              87.47   \n",
       "\n",
       "   pressure_altimiter_in  RunwayTrueTrack Landing_Visibilty_Class_RVR  \\\n",
       "0                  30.21        89.713523                         VMC   \n",
       "1                  30.03        89.713523                         VMC   \n",
       "2                  30.09        89.713523                         VMC   \n",
       "3                  29.38       269.713523                         VMC   \n",
       "4                  29.41       269.713523                         VMC   \n",
       "\n",
       "  present_weather_codes Landing_RWY  \n",
       "0                   NaN       RW09R  \n",
       "1                   NaN       RW09L  \n",
       "2                   NaN       RW09L  \n",
       "3                   NaN       RW27L  \n",
       "4                   NaN       RW27L  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No. of Unique Weather Codes: 49 \n",
      "\n",
      "Unique Weather Codes: [nan 'VCFG' 'BR' '-RA' '-RADZ' 'TSRA' 'HZ' 'VCSH' '-SHRA' '-DZ' '-DZ BR' 'SHRA' 'MIFG' 'FG' 'RA' 'BR BCFG' 'PRFG VCFG' 'RADZ' 'VCTS -SHRA' 'VCTS' '-TSRA' '+SHRA' 'MIFG BR' 'FZFG' '-SN' '//////' '-SHRA VCTS' 'BR VCFG' 'BCFG' 'PRFG' '+RA' 'BR PRFG' 'SHRA VCTS' '-RADZ BR' '-SHRASN' 'RA HZ' '-SHRA TS' '-RA VCTS' 'TS' '-RADZ FG' 'DZ' '-RADZ HZ' 'VCTS VCSH' 'BR MIFG' '-DZ FG' '+RA VCTS' 'SHRA HZ' 'SG' 'BR VCTS' 'RADZ FG'] \n",
      "\n"
     ]
    }
   ],
   "source": [
    "print('No. of Unique Weather Codes:',df.present_weather_codes.nunique(),'\\n')\n",
    "print('Unique Weather Codes:',df.present_weather_codes.unique(),'\\n')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No. of Unique Visibility Classes: 5\n",
      "Unique Visibility Classes: ['VMC' 'Marginal VMC' 'CAT1' 'CAT2' 'CAT3']\n"
     ]
    }
   ],
   "source": [
    "print('No. of Unique Visibility Classes:',df.Landing_Visibilty_Class_RVR.nunique())\n",
    "print('Unique Visibility Classes:',df.Landing_Visibilty_Class_RVR.unique())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No. of Unique Aircraft Registrations landing in EGLL from 01-01-2016 to 31-12-2016: 3511\n"
     ]
    }
   ],
   "source": [
    "print('No. of Unique Aircraft Registrations landing in EGLL from 01-01-2016 to 31-12-2016:',df.aircraft_registration.nunique())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Landing_Time</th>\n",
       "      <th>hourofday</th>\n",
       "      <th>aircraft_registration</th>\n",
       "      <th>wind_direction_degrees</th>\n",
       "      <th>Wind_Speed_kt</th>\n",
       "      <th>Wind_And_Gust_kt</th>\n",
       "      <th>Landing_RVR_or_Vis_m</th>\n",
       "      <th>temperature_c</th>\n",
       "      <th>dew_point_c</th>\n",
       "      <th>relative_humidity</th>\n",
       "      <th>pressure_altimiter_in</th>\n",
       "      <th>RunwayTrueTrack</th>\n",
       "      <th>Landing_Visibilty_Class_RVR</th>\n",
       "      <th>present_weather_codes</th>\n",
       "      <th>Landing_RWY</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1627</th>\n",
       "      <td>2016-05-31 14:41:42</td>\n",
       "      <td>14</td>\n",
       "      <td>IBIXL</td>\n",
       "      <td>360.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>87.56</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2724</th>\n",
       "      <td>2016-04-30 16:02:36</td>\n",
       "      <td>16</td>\n",
       "      <td>GEUXL</td>\n",
       "      <td>360.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>46.59</td>\n",
       "      <td>30.21</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>VCSH</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3215</th>\n",
       "      <td>2016-09-30 18:46:00</td>\n",
       "      <td>18</td>\n",
       "      <td>GEUUT</td>\n",
       "      <td>210.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>81.99</td>\n",
       "      <td>29.85</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4353</th>\n",
       "      <td>2016-09-30 06:48:10</td>\n",
       "      <td>6</td>\n",
       "      <td>GEUPC</td>\n",
       "      <td>230.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>87.47</td>\n",
       "      <td>29.91</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7032</th>\n",
       "      <td>2016-01-31 17:17:31</td>\n",
       "      <td>17</td>\n",
       "      <td>GZBKE</td>\n",
       "      <td>250.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>81.99</td>\n",
       "      <td>29.77</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7768</th>\n",
       "      <td>2016-09-30 17:54:11</td>\n",
       "      <td>17</td>\n",
       "      <td>DAIZV</td>\n",
       "      <td>240.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>14.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>76.82</td>\n",
       "      <td>29.85</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14258</th>\n",
       "      <td>2016-06-30 13:14:47</td>\n",
       "      <td>13</td>\n",
       "      <td>GEUUX</td>\n",
       "      <td>230.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>63.82</td>\n",
       "      <td>29.83</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15237</th>\n",
       "      <td>2016-04-30 21:37:20</td>\n",
       "      <td>21</td>\n",
       "      <td>EIDSE</td>\n",
       "      <td>270.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>65.83</td>\n",
       "      <td>30.33</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17600</th>\n",
       "      <td>2016-05-31 09:06:02</td>\n",
       "      <td>9</td>\n",
       "      <td>GEUYU</td>\n",
       "      <td>350.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>29.97</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19672</th>\n",
       "      <td>2016-03-31 19:58:23</td>\n",
       "      <td>19</td>\n",
       "      <td>PHBGG</td>\n",
       "      <td>120.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>70.48</td>\n",
       "      <td>30.06</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20049</th>\n",
       "      <td>2016-07-31 12:31:28</td>\n",
       "      <td>12</td>\n",
       "      <td>GEUUG</td>\n",
       "      <td>280.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>21.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>43.13</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20341</th>\n",
       "      <td>2016-10-31 21:41:35</td>\n",
       "      <td>21</td>\n",
       "      <td>EIDEP</td>\n",
       "      <td>20.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>3300.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>30.27</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>BR</td>\n",
       "      <td>RW09R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22686</th>\n",
       "      <td>2016-11-30 19:14:57</td>\n",
       "      <td>19</td>\n",
       "      <td>GEUPZ</td>\n",
       "      <td>280.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>9000.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>80.51</td>\n",
       "      <td>30.53</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23457</th>\n",
       "      <td>2016-07-31 16:39:47</td>\n",
       "      <td>16</td>\n",
       "      <td>HL7628</td>\n",
       "      <td>280.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>43.42</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23814</th>\n",
       "      <td>2016-11-30 17:36:56</td>\n",
       "      <td>17</td>\n",
       "      <td>GSTBD</td>\n",
       "      <td>280.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>8000.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>69.88</td>\n",
       "      <td>30.53</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25074</th>\n",
       "      <td>2016-01-31 11:13:29</td>\n",
       "      <td>11</td>\n",
       "      <td>TFFII</td>\n",
       "      <td>210.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>3800.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>93.45</td>\n",
       "      <td>29.68</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>HZ</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25079</th>\n",
       "      <td>2016-01-31 20:34:13</td>\n",
       "      <td>20</td>\n",
       "      <td>PHBXI</td>\n",
       "      <td>250.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>81.99</td>\n",
       "      <td>29.83</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27001</th>\n",
       "      <td>2016-09-30 07:59:45</td>\n",
       "      <td>7</td>\n",
       "      <td>GTTOE</td>\n",
       "      <td>230.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>81.99</td>\n",
       "      <td>29.91</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27417</th>\n",
       "      <td>2016-02-29 12:35:00</td>\n",
       "      <td>12</td>\n",
       "      <td>VHOQB</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>53.00</td>\n",
       "      <td>30.33</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30846</th>\n",
       "      <td>2016-02-29 07:16:09</td>\n",
       "      <td>7</td>\n",
       "      <td>N177DN</td>\n",
       "      <td>350.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>92.97</td>\n",
       "      <td>30.33</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33087</th>\n",
       "      <td>2016-02-29 11:59:43</td>\n",
       "      <td>11</td>\n",
       "      <td>GTTOB</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>53.00</td>\n",
       "      <td>30.33</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34805</th>\n",
       "      <td>2016-08-31 07:29:32</td>\n",
       "      <td>7</td>\n",
       "      <td>N721AN</td>\n",
       "      <td>240.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>82.37</td>\n",
       "      <td>30.12</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35550</th>\n",
       "      <td>2016-04-30 08:06:56</td>\n",
       "      <td>8</td>\n",
       "      <td>GEUYF</td>\n",
       "      <td>330.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>70.68</td>\n",
       "      <td>30.15</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37635</th>\n",
       "      <td>2016-03-31 15:13:40</td>\n",
       "      <td>15</td>\n",
       "      <td>N966AM</td>\n",
       "      <td>30.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>40.54</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37637</th>\n",
       "      <td>2016-03-31 08:14:58</td>\n",
       "      <td>8</td>\n",
       "      <td>GEUYT</td>\n",
       "      <td>20.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>70.68</td>\n",
       "      <td>29.97</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39006</th>\n",
       "      <td>2016-05-31 11:20:23</td>\n",
       "      <td>11</td>\n",
       "      <td>GEUXK</td>\n",
       "      <td>350.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>4800.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>RA</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45492</th>\n",
       "      <td>2016-03-31 11:36:25</td>\n",
       "      <td>11</td>\n",
       "      <td>CSTTC</td>\n",
       "      <td>360.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>53.81</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54973</th>\n",
       "      <td>2016-09-30 09:07:17</td>\n",
       "      <td>9</td>\n",
       "      <td>CGSCA</td>\n",
       "      <td>240.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>15.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>72.01</td>\n",
       "      <td>29.91</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55731</th>\n",
       "      <td>2016-05-31 10:53:16</td>\n",
       "      <td>10</td>\n",
       "      <td>ODMEE</td>\n",
       "      <td>350.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>5000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55735</th>\n",
       "      <td>2016-05-31 17:54:04</td>\n",
       "      <td>17</td>\n",
       "      <td>EIDSA</td>\n",
       "      <td>360.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>87.56</td>\n",
       "      <td>30.06</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>180885</th>\n",
       "      <td>2016-10-31 14:26:36</td>\n",
       "      <td>14</td>\n",
       "      <td>BKQA</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>59.74</td>\n",
       "      <td>30.30</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>185093</th>\n",
       "      <td>2016-09-30 12:26:38</td>\n",
       "      <td>12</td>\n",
       "      <td>N223UA</td>\n",
       "      <td>230.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>55.37</td>\n",
       "      <td>29.88</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>185568</th>\n",
       "      <td>2016-02-29 12:30:01</td>\n",
       "      <td>12</td>\n",
       "      <td>4XECF</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>53.00</td>\n",
       "      <td>30.33</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>186788</th>\n",
       "      <td>2016-03-31 12:29:12</td>\n",
       "      <td>12</td>\n",
       "      <td>GEUUN</td>\n",
       "      <td>360.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>50.08</td>\n",
       "      <td>29.97</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>187145</th>\n",
       "      <td>2016-07-31 07:51:02</td>\n",
       "      <td>7</td>\n",
       "      <td>GEUYI</td>\n",
       "      <td>350.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>63.38</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>194749</th>\n",
       "      <td>2016-07-31 15:30:25</td>\n",
       "      <td>15</td>\n",
       "      <td>GVGEM</td>\n",
       "      <td>280.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>21.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>46.16</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>195032</th>\n",
       "      <td>2016-10-31 14:00:26</td>\n",
       "      <td>14</td>\n",
       "      <td>GEUYJ</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>9000.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>63.60</td>\n",
       "      <td>30.30</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>196337</th>\n",
       "      <td>2016-01-31 15:39:05</td>\n",
       "      <td>15</td>\n",
       "      <td>EIDES</td>\n",
       "      <td>250.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>14.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>76.82</td>\n",
       "      <td>29.71</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>196427</th>\n",
       "      <td>2016-02-29 19:16:46</td>\n",
       "      <td>19</td>\n",
       "      <td>GEUPE</td>\n",
       "      <td>180.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>56.48</td>\n",
       "      <td>30.27</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>197673</th>\n",
       "      <td>2016-04-30 09:16:23</td>\n",
       "      <td>9</td>\n",
       "      <td>GCIVU</td>\n",
       "      <td>310.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>57.52</td>\n",
       "      <td>30.15</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>198570</th>\n",
       "      <td>2016-02-29 09:17:07</td>\n",
       "      <td>9</td>\n",
       "      <td>GEUPX</td>\n",
       "      <td>60.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>70.08</td>\n",
       "      <td>30.36</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>199348</th>\n",
       "      <td>2016-10-31 13:38:16</td>\n",
       "      <td>13</td>\n",
       "      <td>GCIVD</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>7000.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>63.60</td>\n",
       "      <td>30.30</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>200046</th>\n",
       "      <td>2016-06-30 09:30:42</td>\n",
       "      <td>9</td>\n",
       "      <td>GBNLY</td>\n",
       "      <td>230.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>72.56</td>\n",
       "      <td>29.80</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>203128</th>\n",
       "      <td>2016-04-30 20:42:59</td>\n",
       "      <td>20</td>\n",
       "      <td>A7ACA</td>\n",
       "      <td>300.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>65.83</td>\n",
       "      <td>30.30</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>204302</th>\n",
       "      <td>2016-05-31 14:40:47</td>\n",
       "      <td>14</td>\n",
       "      <td>GXLEF</td>\n",
       "      <td>360.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>87.56</td>\n",
       "      <td>30.03</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>204391</th>\n",
       "      <td>2016-06-30 11:31:37</td>\n",
       "      <td>11</td>\n",
       "      <td>GYMMA</td>\n",
       "      <td>230.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>59.97</td>\n",
       "      <td>29.83</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>205027</th>\n",
       "      <td>2016-01-31 16:20:04</td>\n",
       "      <td>16</td>\n",
       "      <td>S2AFP</td>\n",
       "      <td>260.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>81.99</td>\n",
       "      <td>29.74</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>205267</th>\n",
       "      <td>2016-04-30 10:24:50</td>\n",
       "      <td>10</td>\n",
       "      <td>GEUPE</td>\n",
       "      <td>300.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>40.25</td>\n",
       "      <td>30.15</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>205445</th>\n",
       "      <td>2016-06-30 19:18:00</td>\n",
       "      <td>19</td>\n",
       "      <td>GEUOF</td>\n",
       "      <td>210.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>72.38</td>\n",
       "      <td>29.80</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>207560</th>\n",
       "      <td>2016-05-31 09:50:33</td>\n",
       "      <td>9</td>\n",
       "      <td>GEUOE</td>\n",
       "      <td>340.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>5000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>208566</th>\n",
       "      <td>2016-05-31 09:59:26</td>\n",
       "      <td>9</td>\n",
       "      <td>DAIUD</td>\n",
       "      <td>340.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>5000.0</td>\n",
       "      <td>11.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>93.55</td>\n",
       "      <td>30.00</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>208671</th>\n",
       "      <td>2016-06-30 12:14:15</td>\n",
       "      <td>12</td>\n",
       "      <td>N273AY</td>\n",
       "      <td>230.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>63.82</td>\n",
       "      <td>29.83</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>209517</th>\n",
       "      <td>2016-03-31 07:19:10</td>\n",
       "      <td>7</td>\n",
       "      <td>GCIVY</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>75.50</td>\n",
       "      <td>29.97</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>210471</th>\n",
       "      <td>2016-02-29 10:55:34</td>\n",
       "      <td>10</td>\n",
       "      <td>GXLEH</td>\n",
       "      <td>50.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>7.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>56.74</td>\n",
       "      <td>30.33</td>\n",
       "      <td>89.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW09L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>211327</th>\n",
       "      <td>2016-11-30 13:10:31</td>\n",
       "      <td>13</td>\n",
       "      <td>GMIDY</td>\n",
       "      <td>260.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>60.79</td>\n",
       "      <td>30.53</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>216402</th>\n",
       "      <td>2016-08-31 11:41:58</td>\n",
       "      <td>11</td>\n",
       "      <td>EIDVJ</td>\n",
       "      <td>260.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>23.0</td>\n",
       "      <td>12.0</td>\n",
       "      <td>49.90</td>\n",
       "      <td>30.12</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>218169</th>\n",
       "      <td>2016-04-30 19:48:44</td>\n",
       "      <td>19</td>\n",
       "      <td>B6075</td>\n",
       "      <td>NaN</td>\n",
       "      <td>3.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>49.81</td>\n",
       "      <td>30.27</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>219617</th>\n",
       "      <td>2016-08-31 08:34:43</td>\n",
       "      <td>8</td>\n",
       "      <td>GVIIJ</td>\n",
       "      <td>270.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>13.0</td>\n",
       "      <td>72.56</td>\n",
       "      <td>30.12</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27L</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>219711</th>\n",
       "      <td>2016-09-30 17:50:35</td>\n",
       "      <td>17</td>\n",
       "      <td>GEUPZ</td>\n",
       "      <td>240.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>14.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>76.82</td>\n",
       "      <td>29.85</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>-RA</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>219874</th>\n",
       "      <td>2016-11-30 11:12:06</td>\n",
       "      <td>11</td>\n",
       "      <td>GEUXF</td>\n",
       "      <td>240.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>10000.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>80.51</td>\n",
       "      <td>30.56</td>\n",
       "      <td>269.713523</td>\n",
       "      <td>VMC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>RW27R</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>129 rows × 15 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              Landing_Time  hourofday aircraft_registration  \\\n",
       "1627   2016-05-31 14:41:42         14                 IBIXL   \n",
       "2724   2016-04-30 16:02:36         16                 GEUXL   \n",
       "3215   2016-09-30 18:46:00         18                 GEUUT   \n",
       "4353   2016-09-30 06:48:10          6                 GEUPC   \n",
       "7032   2016-01-31 17:17:31         17                 GZBKE   \n",
       "7768   2016-09-30 17:54:11         17                 DAIZV   \n",
       "14258  2016-06-30 13:14:47         13                 GEUUX   \n",
       "15237  2016-04-30 21:37:20         21                 EIDSE   \n",
       "17600  2016-05-31 09:06:02          9                 GEUYU   \n",
       "19672  2016-03-31 19:58:23         19                 PHBGG   \n",
       "20049  2016-07-31 12:31:28         12                 GEUUG   \n",
       "20341  2016-10-31 21:41:35         21                 EIDEP   \n",
       "22686  2016-11-30 19:14:57         19                 GEUPZ   \n",
       "23457  2016-07-31 16:39:47         16                HL7628   \n",
       "23814  2016-11-30 17:36:56         17                 GSTBD   \n",
       "25074  2016-01-31 11:13:29         11                 TFFII   \n",
       "25079  2016-01-31 20:34:13         20                 PHBXI   \n",
       "27001  2016-09-30 07:59:45          7                 GTTOE   \n",
       "27417  2016-02-29 12:35:00         12                 VHOQB   \n",
       "30846  2016-02-29 07:16:09          7                N177DN   \n",
       "33087  2016-02-29 11:59:43         11                 GTTOB   \n",
       "34805  2016-08-31 07:29:32          7                N721AN   \n",
       "35550  2016-04-30 08:06:56          8                 GEUYF   \n",
       "37635  2016-03-31 15:13:40         15                N966AM   \n",
       "37637  2016-03-31 08:14:58          8                 GEUYT   \n",
       "39006  2016-05-31 11:20:23         11                 GEUXK   \n",
       "45492  2016-03-31 11:36:25         11                 CSTTC   \n",
       "54973  2016-09-30 09:07:17          9                 CGSCA   \n",
       "55731  2016-05-31 10:53:16         10                 ODMEE   \n",
       "55735  2016-05-31 17:54:04         17                 EIDSA   \n",
       "...                    ...        ...                   ...   \n",
       "180885 2016-10-31 14:26:36         14                  BKQA   \n",
       "185093 2016-09-30 12:26:38         12                N223UA   \n",
       "185568 2016-02-29 12:30:01         12                 4XECF   \n",
       "186788 2016-03-31 12:29:12         12                 GEUUN   \n",
       "187145 2016-07-31 07:51:02          7                 GEUYI   \n",
       "194749 2016-07-31 15:30:25         15                 GVGEM   \n",
       "195032 2016-10-31 14:00:26         14                 GEUYJ   \n",
       "196337 2016-01-31 15:39:05         15                 EIDES   \n",
       "196427 2016-02-29 19:16:46         19                 GEUPE   \n",
       "197673 2016-04-30 09:16:23          9                 GCIVU   \n",
       "198570 2016-02-29 09:17:07          9                 GEUPX   \n",
       "199348 2016-10-31 13:38:16         13                 GCIVD   \n",
       "200046 2016-06-30 09:30:42          9                 GBNLY   \n",
       "203128 2016-04-30 20:42:59         20                 A7ACA   \n",
       "204302 2016-05-31 14:40:47         14                 GXLEF   \n",
       "204391 2016-06-30 11:31:37         11                 GYMMA   \n",
       "205027 2016-01-31 16:20:04         16                 S2AFP   \n",
       "205267 2016-04-30 10:24:50         10                 GEUPE   \n",
       "205445 2016-06-30 19:18:00         19                 GEUOF   \n",
       "207560 2016-05-31 09:50:33          9                 GEUOE   \n",
       "208566 2016-05-31 09:59:26          9                 DAIUD   \n",
       "208671 2016-06-30 12:14:15         12                N273AY   \n",
       "209517 2016-03-31 07:19:10          7                 GCIVY   \n",
       "210471 2016-02-29 10:55:34         10                 GXLEH   \n",
       "211327 2016-11-30 13:10:31         13                 GMIDY   \n",
       "216402 2016-08-31 11:41:58         11                 EIDVJ   \n",
       "218169 2016-04-30 19:48:44         19                 B6075   \n",
       "219617 2016-08-31 08:34:43          8                 GVIIJ   \n",
       "219711 2016-09-30 17:50:35         17                 GEUPZ   \n",
       "219874 2016-11-30 11:12:06         11                 GEUXF   \n",
       "\n",
       "        wind_direction_degrees  Wind_Speed_kt  Wind_And_Gust_kt  \\\n",
       "1627                     360.0            9.0               9.0   \n",
       "2724                     360.0            8.0               8.0   \n",
       "3215                     210.0            6.0               6.0   \n",
       "4353                     230.0            7.0               7.0   \n",
       "7032                     250.0           16.0              16.0   \n",
       "7768                     240.0            8.0               8.0   \n",
       "14258                    230.0           13.0              13.0   \n",
       "15237                    270.0            6.0               6.0   \n",
       "17600                    350.0           10.0              10.0   \n",
       "19672                    120.0            7.0               7.0   \n",
       "20049                    280.0            7.0               7.0   \n",
       "20341                     20.0            2.0               2.0   \n",
       "22686                    280.0            4.0               4.0   \n",
       "23457                    280.0            7.0               7.0   \n",
       "23814                    280.0            3.0               3.0   \n",
       "25074                    210.0           12.0              12.0   \n",
       "25079                    250.0           17.0              17.0   \n",
       "27001                    230.0            7.0               7.0   \n",
       "27417                      NaN            2.0               2.0   \n",
       "30846                    350.0            4.0               4.0   \n",
       "33087                      NaN            3.0               3.0   \n",
       "34805                    240.0            5.0               5.0   \n",
       "35550                    330.0            4.0               4.0   \n",
       "37635                     30.0            6.0               6.0   \n",
       "37637                     20.0            5.0               5.0   \n",
       "39006                    350.0            8.0               8.0   \n",
       "45492                    360.0            5.0               5.0   \n",
       "54973                    240.0            8.0               8.0   \n",
       "55731                    350.0            9.0               9.0   \n",
       "55735                    360.0           10.0              10.0   \n",
       "...                        ...            ...               ...   \n",
       "180885                     NaN            1.0               1.0   \n",
       "185093                   230.0           13.0              13.0   \n",
       "185568                     NaN            2.0               2.0   \n",
       "186788                   360.0            6.0               6.0   \n",
       "187145                   350.0            4.0               4.0   \n",
       "194749                   280.0            7.0               7.0   \n",
       "195032                     NaN            2.0               2.0   \n",
       "196337                   250.0           16.0              16.0   \n",
       "196427                   180.0            7.0               7.0   \n",
       "197673                   310.0            6.0               6.0   \n",
       "198570                    60.0            5.0               5.0   \n",
       "199348                     NaN            1.0               1.0   \n",
       "200046                   230.0            6.0               6.0   \n",
       "203128                   300.0            4.0               4.0   \n",
       "204302                   360.0            9.0               9.0   \n",
       "204391                   230.0           12.0              12.0   \n",
       "205027                   260.0           17.0              17.0   \n",
       "205267                   300.0            7.0               7.0   \n",
       "205445                   210.0           13.0              13.0   \n",
       "207560                   340.0            9.0               9.0   \n",
       "208566                   340.0            9.0               9.0   \n",
       "208671                   230.0           12.0              12.0   \n",
       "209517                     NaN            3.0               3.0   \n",
       "210471                    50.0            6.0               6.0   \n",
       "211327                   260.0            4.0               4.0   \n",
       "216402                   260.0            8.0               8.0   \n",
       "218169                     NaN            3.0               3.0   \n",
       "219617                   270.0            6.0               6.0   \n",
       "219711                   240.0            8.0               8.0   \n",
       "219874                   240.0            4.0               4.0   \n",
       "\n",
       "        Landing_RVR_or_Vis_m  temperature_c  dew_point_c  relative_humidity  \\\n",
       "1627                 10000.0           12.0         10.0              87.56   \n",
       "2724                 10000.0           11.0          0.0              46.59   \n",
       "3215                 10000.0           13.0         10.0              81.99   \n",
       "4353                 10000.0           11.0          9.0              87.47   \n",
       "7032                 10000.0           13.0         10.0              81.99   \n",
       "7768                 10000.0           14.0         10.0              76.82   \n",
       "14258                10000.0           19.0         12.0              63.82   \n",
       "15237                10000.0            8.0          2.0              65.83   \n",
       "17600                10000.0           11.0         10.0              93.55   \n",
       "19672                10000.0            7.0          2.0              70.48   \n",
       "20049                10000.0           21.0          8.0              43.13   \n",
       "20341                 3300.0           11.0         10.0              93.55   \n",
       "22686                 9000.0            2.0         -1.0              80.51   \n",
       "23457                10000.0           22.0          9.0              43.42   \n",
       "23814                 8000.0            4.0         -1.0              69.88   \n",
       "25074                 3800.0            9.0          8.0              93.45   \n",
       "25079                10000.0           13.0         10.0              81.99   \n",
       "27001                10000.0           13.0         10.0              81.99   \n",
       "27417                10000.0            8.0         -1.0              53.00   \n",
       "30846                10000.0            0.0         -1.0              92.97   \n",
       "33087                10000.0            8.0         -1.0              53.00   \n",
       "34805                10000.0           16.0         13.0              82.37   \n",
       "35550                10000.0            8.0          3.0              70.68   \n",
       "37635                10000.0           12.0         -1.0              40.54   \n",
       "37637                10000.0            8.0          3.0              70.68   \n",
       "39006                 4800.0           11.0         10.0              93.55   \n",
       "45492                10000.0           11.0          2.0              53.81   \n",
       "54973                10000.0           15.0         10.0              72.01   \n",
       "55731                 5000.0           11.0         10.0              93.55   \n",
       "55735                10000.0           12.0         10.0              87.56   \n",
       "...                      ...            ...          ...                ...   \n",
       "180885               10000.0           19.0         11.0              59.74   \n",
       "185093               10000.0           17.0          8.0              55.37   \n",
       "185568               10000.0            8.0         -1.0              53.00   \n",
       "186788               10000.0           11.0          1.0              50.08   \n",
       "187145               10000.0           17.0         10.0              63.38   \n",
       "194749               10000.0           21.0          9.0              46.16   \n",
       "195032                9000.0           18.0         11.0              63.60   \n",
       "196337               10000.0           14.0         10.0              76.82   \n",
       "196427               10000.0            6.0         -2.0              56.48   \n",
       "197673               10000.0           10.0          2.0              57.52   \n",
       "198570               10000.0            5.0          0.0              70.08   \n",
       "199348                7000.0           18.0         11.0              63.60   \n",
       "200046               10000.0           18.0         13.0              72.56   \n",
       "203128               10000.0            8.0          2.0              65.83   \n",
       "204302               10000.0           12.0         10.0              87.56   \n",
       "204391               10000.0           20.0         12.0              59.97   \n",
       "205027               10000.0           13.0         10.0              81.99   \n",
       "205267               10000.0           11.0         -2.0              40.25   \n",
       "205445               10000.0           17.0         12.0              72.38   \n",
       "207560                5000.0           11.0         10.0              93.55   \n",
       "208566                5000.0           11.0         10.0              93.55   \n",
       "208671               10000.0           19.0         12.0              63.82   \n",
       "209517               10000.0            6.0          2.0              75.50   \n",
       "210471               10000.0            7.0         -1.0              56.74   \n",
       "211327               10000.0            6.0         -1.0              60.79   \n",
       "216402               10000.0           23.0         12.0              49.90   \n",
       "218169               10000.0           10.0          0.0              49.81   \n",
       "219617               10000.0           18.0         13.0              72.56   \n",
       "219711               10000.0           14.0         10.0              76.82   \n",
       "219874               10000.0            2.0         -1.0              80.51   \n",
       "\n",
       "        pressure_altimiter_in  RunwayTrueTrack Landing_Visibilty_Class_RVR  \\\n",
       "1627                    30.03       269.713523                         VMC   \n",
       "2724                    30.21       269.713523                         VMC   \n",
       "3215                    29.85       269.713523                         VMC   \n",
       "4353                    29.91       269.713523                         VMC   \n",
       "7032                    29.77       269.713523                         VMC   \n",
       "7768                    29.85       269.713523                         VMC   \n",
       "14258                   29.83       269.713523                         VMC   \n",
       "15237                   30.33       269.713523                         VMC   \n",
       "17600                   29.97       269.713523                         VMC   \n",
       "19672                   30.06       269.713523                         VMC   \n",
       "20049                   30.03       269.713523                         VMC   \n",
       "20341                   30.27        89.713523                         VMC   \n",
       "22686                   30.53       269.713523                         VMC   \n",
       "23457                   30.03       269.713523                         VMC   \n",
       "23814                   30.53       269.713523                         VMC   \n",
       "25074                   29.68       269.713523                         VMC   \n",
       "25079                   29.83       269.713523                         VMC   \n",
       "27001                   29.91       269.713523                         VMC   \n",
       "27417                   30.33        89.713523                         VMC   \n",
       "30846                   30.33        89.713523                         VMC   \n",
       "33087                   30.33        89.713523                         VMC   \n",
       "34805                   30.12       269.713523                         VMC   \n",
       "35550                   30.15       269.713523                         VMC   \n",
       "37635                   30.00       269.713523                         VMC   \n",
       "37637                   29.97       269.713523                         VMC   \n",
       "39006                   30.00       269.713523                         VMC   \n",
       "45492                   30.00       269.713523                         VMC   \n",
       "54973                   29.91       269.713523                         VMC   \n",
       "55731                   30.00       269.713523                         VMC   \n",
       "55735                   30.06       269.713523                         VMC   \n",
       "...                       ...              ...                         ...   \n",
       "180885                  30.30        89.713523                         VMC   \n",
       "185093                  29.88       269.713523                         VMC   \n",
       "185568                  30.33        89.713523                         VMC   \n",
       "186788                  29.97       269.713523                         VMC   \n",
       "187145                  30.03       269.713523                         VMC   \n",
       "194749                  30.03       269.713523                         VMC   \n",
       "195032                  30.30        89.713523                         VMC   \n",
       "196337                  29.71       269.713523                         VMC   \n",
       "196427                  30.27       269.713523                         VMC   \n",
       "197673                  30.15       269.713523                         VMC   \n",
       "198570                  30.36        89.713523                         VMC   \n",
       "199348                  30.30        89.713523                         VMC   \n",
       "200046                  29.80       269.713523                         VMC   \n",
       "203128                  30.30       269.713523                         VMC   \n",
       "204302                  30.03       269.713523                         VMC   \n",
       "204391                  29.83       269.713523                         VMC   \n",
       "205027                  29.74       269.713523                         VMC   \n",
       "205267                  30.15       269.713523                         VMC   \n",
       "205445                  29.80       269.713523                         VMC   \n",
       "207560                  30.00       269.713523                         VMC   \n",
       "208566                  30.00       269.713523                         VMC   \n",
       "208671                  29.83       269.713523                         VMC   \n",
       "209517                  29.97       269.713523                         VMC   \n",
       "210471                  30.33        89.713523                         VMC   \n",
       "211327                  30.53       269.713523                         VMC   \n",
       "216402                  30.12       269.713523                         VMC   \n",
       "218169                  30.27       269.713523                         VMC   \n",
       "219617                  30.12       269.713523                         VMC   \n",
       "219711                  29.85       269.713523                         VMC   \n",
       "219874                  30.56       269.713523                         VMC   \n",
       "\n",
       "       present_weather_codes Landing_RWY  \n",
       "1627                     NaN       RW27L  \n",
       "2724                    VCSH       RW27R  \n",
       "3215                     NaN       RW27R  \n",
       "4353                     NaN       RW27L  \n",
       "7032                     NaN       RW27L  \n",
       "7768                     -RA       RW27R  \n",
       "14258                    NaN       RW27R  \n",
       "15237                    NaN       RW27R  \n",
       "17600                    -RA       RW27R  \n",
       "19672                    NaN       RW27R  \n",
       "20049                    NaN       RW27R  \n",
       "20341                     BR       RW09R  \n",
       "22686                    NaN       RW27L  \n",
       "23457                    NaN       RW27L  \n",
       "23814                    NaN       RW27L  \n",
       "25074                     HZ       RW27R  \n",
       "25079                    NaN       RW27L  \n",
       "27001                    NaN       RW27L  \n",
       "27417                    NaN       RW09L  \n",
       "30846                    NaN       RW09L  \n",
       "33087                    NaN       RW09L  \n",
       "34805                    NaN       RW27L  \n",
       "35550                    NaN       RW27L  \n",
       "37635                    NaN       RW27R  \n",
       "37637                    NaN       RW27L  \n",
       "39006                     RA       RW27L  \n",
       "45492                    NaN       RW27L  \n",
       "54973                    NaN       RW27L  \n",
       "55731                    -RA       RW27R  \n",
       "55735                    NaN       RW27L  \n",
       "...                      ...         ...  \n",
       "180885                   NaN       RW09L  \n",
       "185093                   NaN       RW27L  \n",
       "185568                   NaN       RW09L  \n",
       "186788                   NaN       RW27L  \n",
       "187145                   NaN       RW27R  \n",
       "194749                   NaN       RW27L  \n",
       "195032                   NaN       RW09L  \n",
       "196337                   NaN       RW27L  \n",
       "196427                   NaN       RW27R  \n",
       "197673                   NaN       RW27L  \n",
       "198570                   NaN       RW09L  \n",
       "199348                   NaN       RW09L  \n",
       "200046                   NaN       RW27R  \n",
       "203128                   NaN       RW27R  \n",
       "204302                   NaN       RW27R  \n",
       "204391                   NaN       RW27R  \n",
       "205027                   NaN       RW27L  \n",
       "205267                   NaN       RW27L  \n",
       "205445                   NaN       RW27L  \n",
       "207560                   -RA       RW27R  \n",
       "208566                   -RA       RW27R  \n",
       "208671                   NaN       RW27R  \n",
       "209517                   NaN       RW27L  \n",
       "210471                   NaN       RW09L  \n",
       "211327                   NaN       RW27R  \n",
       "216402                   NaN       RW27L  \n",
       "218169                   NaN       RW27R  \n",
       "219617                   NaN       RW27L  \n",
       "219711                   -RA       RW27R  \n",
       "219874                   NaN       RW27R  \n",
       "\n",
       "[129 rows x 15 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Check for Duplicates in the Dataframe\n",
    "df[df.duplicated(keep='first',subset=['Landing_Time','aircraft_registration'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Grouping and sorting by ICAO values and postime\n",
    "df.drop_duplicates(keep='first',subset=['Landing_Time','aircraft_registration'], inplace=True)\n",
    "df.sort_values(by=['Landing_Time'], inplace = True)\n",
    "df.to_csv('data/2016_EGLL_Testdata_cleaned.csv', index = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x246136b2748>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAApUAAAJCCAYAAACRajwwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xm8XWV97/HPNxEFBAmDWkQUiiDKFCRSGaSoiNqqaEXBoQJV0dbZCxWHq6hFsXrrrUPVyEUcAVGoUdCgjMogBAhhEAEBK0QrMqRMIiS/+8daRzeHc072yTonOzt83q/Xep2113rW8/z2Tl7J7/ye9aydqkKSJEnqYsagA5AkSdLwM6mUJElSZyaVkiRJ6sykUpIkSZ2ZVEqSJKkzk0pJkiR1ZlIpSZI0hJIcneR3SS4f53ySfDrJtUkWJXlaz7kDklzTbgdMRTwmlZIkScPpGOD5E5x/AbBlux0MfB4gyQbAB4G/AnYGPphk/a7BmFRKkiQNoao6G7h1gib7AF+txvnArCQbA88DflRVt1bVbcCPmDg57YtJpSRJ0uppE+DXPa9vbI+Nd7yTh3XtQKufa3Z/3tB8d+fP//1Tgw6hb5tuNGvQIUzKMr/CddoM00d7zW9+N+gQ+nbxdb9efqNVyO5bbzHoEPq21sPXGHQIk7L39ltlZY43Xf9vbnXOqW+kmbYeMbeq5k6ii7E+h5rgeCcmlZIkSaugNoGcTBI52o3Apj2vHw8sbo/vOer4mR3GAZz+liRJ6iYzpmfrbh7w2nYV+DOAJVX1G2A+sHeS9dsFOnu3xzqxUilJktRFVupse8+wOZam4rhRkhtpVnSvAVBVXwBOAf4GuBa4GzioPXdrko8AF7ZdfbiqJlrw0xeTSkmSpCFUVa9czvkC3jzOuaOBo6cyHpNKSZKkDjJjMJXKVY33VEqSJKkzK5WSJEldTM2imqHnpyBJkqTOrFRKkiR1MaDV36sak0pJkqQuXKgDOP0tSZKkKWClUpIkqYM4/Q1YqZQkSdIUsFIpSZLUxQxrdGBSKUmS1I3T34DT35IkSZoCViolSZK6sFIJWKmctCSbJbl8JY/5iSRXJPnEGOfuXJmxSJIkjcVK5SogycOq6v4JmrwReHRV3buyYpIkSf2JC3UAK5UramaSL7XVw1OTrJVkdpLzkyxKclKS9QGSnJlkTru/UZIb2v0Dk5yQ5HvAqWl8IsnlSS5Lsl/bbh7wSOBnSfZLsnmS85JcmOQjIwElWSfJaUkubq/fpz3+kSRv72l3RJK3rawPSpKk1d6MGdOzDZnhi3jVsCXwuaraBrgdeBnwVeDdVbU9cBnwwT762QU4oKqeDfwdMBvYAdgL+ESSjavqxcA9VTW7qo4H/h34fFU9HfhtT19/AF5aVU8DngX8nzRPY/1/wAEASWYA+wPfGB1IkoOTLEiy4Ljf3jjZz0OSJD3EmVSumOuramG7fxGwBTCrqs5qj30F2KOPfn5UVbe2+7sDx1bV0qr6b+As4OljXLMbcGy7/7We4wE+mmQR8GNgE+CxVXUDcEuSHYG9gUuq6pbRnVbV3KqaU1Vz9v+Lx/cRuiRJApqFOtOxDRnvqVwxvfc2LgVmTdD2fv6cvK856txdPfuT+dtTYxx7NfBoYKequq+dZh8Z7yjgQOAvgKMnMY4kSVJfrFROjSXAbUme2b7+e5pKI8ANwE7t/r4T9HE2sF+SmUkeTVPpvGCMdufQTGFDk0iOWA/4XZtQPgt4Ys+5k4Dn01Q+5/f1jiRJUl+STMs2bKxUTp0DgC8kWRu4DjioPf5J4FtJ/h44fYLrT6K5x/JSmkrkP1fVb8do93bgm+3im+/0HP8G8L0kC4CFwFUjJ6rqj0nOAG6vqqUr9O4kSZImYFI5Se09itv2vP5kz+lnjNH+KmD7nkPvb48fAxzT066AQ9ttdB/r9OxfT5N8jjiyPf77Ucf/pF2g8wzg5eO8LUmStKJmDF9VcTqYVK7mkjwV+D5wUlVdM+h4JEla7cS7CcGkcrVXVVcCfznoOCRJ0urNpFKSJKkLp78BV39LkiRpCliplCRJ6mAYH/8zHUwqJUmSunChDuD0tyRJkqaAlUpJkqQuXKgDWKmUJEnSFLBSKUmS1EFmWKMDk0pJkqRuXP0NOP0tSZKkKWClUpIkqQsrlYCVSkmSJE0BK5WSJElduFAHsFIpSZKkKWClUg/y83//1KBD6NtT3v7OQYfQt1u/ctSgQ1htPWzmcP1+fPtd9ww6hL49au21Bh1C37Z74iaDDmFSbr97eP4e7HLm6YMOYXK232qlDud3fzdMKiVJkrrwG3UAp78lSZI0BaxUSpIkdRFrdGClUpIkSVPASqUkSVIXLtQBTColSZI6iQt1AKe/JUmSNAWsVEqSJHXh9DdgpVKSJElTwEqlJElSF373N2BSKUmS1ElMKgGnvyVJkjQFrFRKkiR14UIdwEqlJEmSpoCVSkmSpC6sVAImlZIkSd24UAdw+luSJElTwEqlJElSB3H6G5iCSmWSU5LMmkT7zZJcPon2ZyaZsyJjTdDn7CR/0/P6xUkO69rvBOMdk2Tf6epfkiRp0DpXKqvqb5bfamqMNVaaXw9SVcsm0dVsYA5wStvvPGDelAQ5jZI8rKruH3QckiSph5VKoI9KZZJ/TvK2dv9TSU5v95+T5OtJbkiyUVuB/HmSLyW5IsmpSdZq2+6U5NIk5wFvXs54ayU5LsmiJMcDa/WcGz3WfwAXA5sm2TvJeUkuTnJCknXaa56e5Nx2/AuSrAd8GNgvycIk+yU5MMln2/ZPTHJaO/5pSZ7QHj8myafbvq6bqPKYxmeTXJnkZOAxPed2SnJWkouSzE+ycU+ci9r38ImRam4b2wlJvgec2h47NMmFbfsP9fT9mvY9LkzyxSQz2+2YJJcnuSzJO5f3Zy5JkjRZ/Ux/nw08s92fA6yTZA1gd+Ano9puCXyuqrYBbgde1h7/MvC2qtqlj/H+Ebi7qrYHjgB2Gqfdk4GvVtWOwF3A+4G9quppwALgXUkeDhwPvL2qdgD2att+ADi+qmZX1fGj+v1s2+/2wDeAT/ec27h93y8EjpzgPby0jW874A3ArgDt5/YZYN+q2gk4un2P0HxGb2o/o6Wj+tsFOKCqnp1kb5rPeWeaiutOSfZI8hRgP2C3qprd9vHqts0mVbVtVW3XjvMgSQ5OsiDJgvknfmuCtyZJkh5gRqZnGzL9TH9fRJO4rAvcS1MZnEOTaL4NeE9P2+uramHPdZu1lcFZVXVWe/xrwAsmGG8P2kSuqhYlWTROu19V1fnt/jOApwLntDfLPhw4jyax+01VXdj29z+w3BtqdwH+rifWf+0595/tNPuVSR67nPdwbFUtBRaPVHfbeLYFftTGMBP4TXuf6LpVdW7b7ps0ieuIH1XVre3+3u12Sft6HZokc3uaBPzCtu+1gN8B3wP+MslngJNpq52jVdVcYC7AvIuurAnemyRJ6hUfpgN9JJVVdV+SG4CDgHOBRcCzgC2An49qfm/P/lKaxCbAZJOUftrf1bMfmsTrlb0Nkmy/AmNPFEvv+1verxBjjRvgitEV2yTrL6ev0e/1Y1X1xVF9vBX4SlW9h1GS7AA8j+bWg1cA/7Cc8SRJkial39T6bOCQ9udPgDcBC6tquQlbVd0OLEmye3vo1X2M9WqAJNvSVOCW53xgtyRPaq9bO8lWwFXA45I8vT2+bpKHAXcA647T17nA/j2x/rSP8cd6D/u39zNuTJOEA/wCeHSSXdp41kiyTVXdBtyR5Bltu/0f3OWfzAf+oeee0U2SPAY4Ddi33SfJBu39oRsBM6rqO8D/Bp62Au9HkiSNIzMyLduw6Xf190+A9wHnVdVdSf7Ag++nnMhBwNFJ7qZJiibyeeDL7bT3QuCC5XVeVTcnORA4Nskj2sPvr6qrk+wHfKZdNHQPzX2VZwCHJVkIfGxUd29rYz0UuLmNfbJOAp4NXAZcDZzVxvnHdoHPp9vbAh4G/F/gCuB1wJeS3AWcCSwZ572e2t4/eV47zX0n8JqqujLJ+4FTk8wA7qOpTN5D83mO/ALxoEqmJElSV+mj2KiVIMk6VXVnu38YsHFVvX0QsQzTPZVPefvwLGa/9StHDTqE1dbDZg7X/Uy33XnPoEPo2733D89TzG65467lN1qFDNP/vy+44LxBhzApjznkrSu1zLf43R+clj/Mx338Q0NVrvQbdVYdf5vkPTR/Jr8CDhxsOJIkqS8DfE5lkucD/06z+Peoqjpy1PlP8efb8NYGHlNVs9pzS2lmVQH+q6pe3CWWgSWVSZ4HfHzU4eur6qWDiGdFJNmOZoV4r3ur6q8m21f7aKPRjzeSJEkaU5KZwOeA5wI30jwBZl5VXTnSpqre2dP+rcCOPV3c0z6GcEoMLKmsqvks//7KVVpVXUbzHEhJkvQQtZxHFU6nnYFrq+q6No7jgH2AK8dp/0rgg9MVzHDdiCRJkvQQ0fvFJO128KgmmwC/7nl9Y3tsrL6eCGwOnN5zeM223/OTvKRrvN5TKUmS1MWM6anR9X4xyTjGKpGOt2hof+Db7RezjHhCVS1O8pfA6Ukuq6pfrmC4ViolSZI6SaZnW74bgU17Xj8eWDxO2/2BY3sPVNXi9ud1NI8z3PHBl/XPpFKSJGk4XQhsmWTzJA+nSRznjW6U5MnA+jRfYT1ybP2RZ3u3X5SyG+Pfi9kXp78lSZK6GNBCnaq6P8lbaBY+zwSOrqorknwYWFBVIwnmK4HjRn0T4lOALyZZRlNkPLJ31fiKMKmUJEkaUlV1CnDKqGMfGPX68DGuOxfYbipjMamUJEnqINO0UGfY+ClIkiSpMyuVkiRJXQzwaxpXJSaVkiRJXcwwqQSnvyVJkjQFrFRKkiR14fQ3YFKpMWy60axBh9C3W79y1KBD6NsGB7x+0CFMzoyZg46gb+ud/O1BhzAp66615qBD6NvNS+4cdAh9m7X2Wpx95bWDDqNvVy/+3aBD6NtzX3/QoEPQEDCplCStFoYpodTqxUcKNUwqJUmSuohJJbhQR5IkSVPASqUkSVIXPlIIsFIpSZKkKWClUpIkqYP4SCHApFKSJKkbF+oATn9LkiRpCliplCRJ6sKFOoCVSkmSJE0BK5WSJElduFAHsFIpSZKkKWClUpIkqYN4TyVgUilJktSNjxQCnP6WJEnSFLBSKUmS1IULdQArlZIkSZoCViolSZK6cKEOYFIpSZLUSWY48Qur0fR3kk8leUfP6/lJjup5/X+SvDfJtyfZ74FJPjvB+ScnOTPJwiQ/TzJ3xd5B3/HsmeT7E5w/PMkhYxw/MMnjpjM2SZL00LXaJJXAucCuAElmABsB2/Sc3xU4rar2neJxPw18qqpmV9VTgM9Mcf9T5UDApFKSpKmWGdOzDZnhi3h859AmlTTJ5OXAHUnWT/II4CnAbUkuhz9V7k5M8sMk1yT515GOkhyU5OokZwG7LWfcjYEbR15U1WU9/X+37f8XST7Y0/9rklzQVje/mGRme3zvJOcluTjJCUnWaY8/P8lVSX4K/F2/H0iSNyT5QZK/B+YA32jHXGuMtgcnWZBkwXe++fV+h5AkSQJWo3sqq2pxkvuTPIEmuTwP2ATYBVgCLAL+OOqy2cCOwL3AL5J8Brgf+BCwU3vdGcAlEwz9KeD0JOcCpwJfrqrb23M7A9sCdwMXJjkZuAvYD9itqu5L8h/Aq5OcArwf2Kuq7krybuBdbbL7JeDZwLXA8f18HkneAuwNvKSq7k3yOuCQqlowVvuqmgvMBbjkV4urnzEkSRIu1GmtNklla6RauSvwbzRJ5a40yeG5Y7Q/raqWACS5EngizbT5mVV1c3v8eGCr8Qasqi8nmQ88H9gHeGOSHdrTP6qqW9p+TgR2p0lad6JJMgHWAn4HPAN4KnBOe/zhNInx1sD1VXVN28/XgYOX8zn8PU319CVVdd9y2kqSpA7icyqB1S+pHLmvcjua6e9fA/8L+B/g6DHa39uzv5Q/fx6TqtRV1eK2/6Pb6fVtx+mngABfqar39J5I8iKaJPSVo47Pnmw8NO99NvB44PpJXitJkjRpq9M9ldBUKl8I3FpVS6vqVmAWzRT4eX328TNgzyQbJlkDePlEjdv7Hddo9/8C2BC4qT393CQbtPcwvqSN7zRg3ySPaa/ZIMkTgfOB3ZI8qT2+dpKtgKuAzZNs0fb5gKRzHJcAbwTm9az4vgNYt49rJUnSZCTTsw2Z1S2pvIxm+vr8UceWVNXv++mgqn4DHE6ThP4YuHg5l+wNXJ7kUmA+cGhV/bY991Pga8BC4DtVtaCqrqS5d/LUJIuAHwEbt9PtBwLHtsfPB7auqj/QTHef3C7U+VWf7+OnwCHtdRsBxwBfGG+hjiRJUhepck3GdEhyIDCnqt4y6Fgma5gW6vzx/qWDDqFvGxzw+kGHMDkzZg46gr6td/KkHj87cLfffc+gQ+jbzUvuHHQIfTv7ymsHHcKkXL34d4MOoW9HvOpFgw5hUh43a52VWua75aivTsv/mxu+/rVDVa5c3SqVkiRJGoDVbaHOtEnyPh58f+UJVXXEWO2r6hiaKefpiucg4O2jDp9TVW+erjElSdIYhvD+x+lgUtmnNnkcM4EchKr6MvDlQcchSdJDnY8Uajj9LUmSpM6sVEqSJHUxwxodWKmUJEnSFLBSKUmS1IX3VAImlZIkSd04/Q04/S1JkqQpYKVSkiSpg8xw+husVEqSJGkKWKmUJEnqwoU6gEmlJElSN3HiF5z+liRJ0hSwUilJktSBC3UaJpV6kGVVgw5h9TRj5qAjmJxlSwcdQd+WLls26BAmZdkQxVtD9O/BDP9jlwbKpFKSJKkLF+oA3lMpSZKkKWClUpIkqQtXfwMmlZIkSd14Py/g9LckSZKmgJVKSZKkDuJCHcBKpSRJkqaAlUpJkqQuvKcSMKmUJEnqZoYTv+D0tyRJkqaAlUpJkqQufE4lYKVSkiRJU8BKpSRJUgc+UqhhUilJktSFq78Bp78lSZI0BUwqJUmSukimZ+tr6Dw/yS+SXJvksDHOH5jk5iQL2+31PecOSHJNux3Q9WNw+luSJGkIJZkJfA54LnAjcGGSeVV15aimx1fVW0ZduwHwQWAOUMBF7bW3rWg8ViolSZK6yIzp2ZZvZ+Daqrquqv4IHAfs02fUzwN+VFW3tonkj4Dnr9D7b60SSWWSTyV5R8/r+UmO6nn9f5K8N8m3J9nvgUk+20e7S5McO7mo/3TtmUnmTHB+nSSfT/LLJJckuSjJG1ZwrJckeeqKxJPkvSsypiRJGowkBydZ0LMdPKrJJsCve17f2B4b7WVJFiX5dpJNJ3lt31aJpBI4F9gVIMkMYCNgm57zuwKnVdW+Uz1wkqfQfA57JHnkVPcPHAXcBmxZVTvS/BawwQr29RJgwqRyAiaVkiRNg8zItGxVNbeq5vRsc0cPPUY4Ner194DNqmp74MfAVyZx7aSsKknlObRJJU0yeTlwR5L1kzwCeApwW5LL4U8VyBOT/LC9ufRfRzpKclCSq5OcBezWx9ivAr4GnAq8uKefM5N8PMkFbX/PbI+vleS4NuM/HlhrvI6TbEFTmn5/VS0DqKqbq+rj7fk9k3y/p/1nkxzY7h+Z5Mp2nE8m2bWN7xPtjbZbTPSmksxI8pUk/5LkSGCt9rpv9PGZSJKkfg1uoc6NwKY9rx8PLO5tUFW3VNW97csvATv1e+1krRJJZVUtBu5P8gSa5PI84GfALjQ3kC4C/jjqstnAfsB2wH5JNk2yMfAhmmTyufRX1dsPOB44FnjlqHMPq6qdgXfQ3MwK8I/A3W3GfwR//sMZyzbApSMJZb/am2dfCmzTjvMvVXUuMA84tKpmV9UvJ+jiYcA3gKur6v1VdRhwT3vdq8cZ808l9hO/+fXJhCtJkgbjQmDLJJsneTiwP02u8CdtbjTixcDP2/35wN5tAW99YO/22ApblVZ/j1QrdwX+jWZef1dgCc30+GinVdUSgCRXAk+kmTY/s6pubo8fD2w13oBJng7cXFW/SnIjcHSS9XtWPp3Y/rwI2Kzd3wP4NEBVLUqyqN83mOR9wMuBx1TV4yZo+j/AH4CjkpwMfH+CtmP5IvCtqjqi3wvakvpcgItuuKlT+VuSpIeUGYOp0VXV/UneQpMMzgSOrqorknwYWFBV84C3JXkxcD9wK3Bge+2tST5Ck5gCfLiqbu0SzypRqWyN3Fe5Hc309/k0lcpdaRLO0e7t2V/KnxPkySRErwS2TnID8EvgUcDLxhijt//JjHElsEN7nyhVdURVzW7HgeYPuPfPYM223f000+bfobmP8od9jjfiXOBZSdac5HWSJGmIVNUpVbVVVW0xUkyqqg+0CSVV9Z6q2qaqdqiqZ1XVVT3XHl1VT2q3L3eNZVVKKs8BXgjcWlVL22x5Fk1ieV6fffwM2DPJhknWoKkKjqlN9F4ObF9Vm1XVZjTL8EdPgY92NvDqto9tge3Ha1hV1wILgH9pnyVFm+iN3CjxK+CpSR6RZD3gOW2bdYD1quoUmqn32W37O4B1lxMfwP8DTgFOSDKSDN/XfiaSJGkKJZmWbdisSknlZTTT1+ePOrakqn7fTwdV9RvgcJok9MfAxRM03wO4qapu6jl2Nk2St/E41wB8Hlinnfb+Z+CC5YT1emBD4NokF7VxvbuN99fAt2juGf0GcEl7zbrA99sxzgLe2R4/Dji0fTTRhAt1qurfaN7/19oEei6wyIU6kiRNsRkzpmcbMqny9jk90DDdU3n/0kmtgRqoDQ5646BDmJxlSwcdQd/WOfk7gw5hUpbcfc+gQ+jb7//nrkGH0LdzfnHdoEOYlKtu/O9Bh9C3I171okGHMCmPm7XOSi3z3XnmT6fl/8119tx9qMqVq9JCHUmSpOEzhFPV0+EhkVT2rLrudcJkVkf3McbPgEeMOvz3VXXZVI0xarzP8eDncP77VNxoK0mSNFkPiaSyTR6nLIEcZ4y/ms7+xxjvzStzPEmSNI4ZVirhIZJUSpIkTZf2yYEPeX4KkiRJ6sxKpSRJUhcu1AGsVEqSJGkKWKmUJEnqwoU6gJVKSZIkTQErlZIkSV24+hswqZQkSeokTn8DTn9LkiRpCliplCRJ6sJHCgFWKiVJkjQFrFRKkiR1YaUSMKnUkHvYzOEptq938rcHHcKkLF22bNAh9O3Ov33ZoEOYlA1+cNKgQ+jb1YtvHnQIfTvrimsHHcKkPGe7rQYdQt8W37Zk0CFMyuNmrbNSx8uM4fm/aDr5KUiSJKkzK5WSJEldWKkErFRKkiRpCliplCRJ6sKFOoBJpSRJUjd+ow7g9LckSZKmgJVKSZKkDhJrdGClUpIkSVPASqUkSVIXLtQBrFRKkiRpCliplCRJ6sLV34BJpSRJUjdOfwNOf0uSJGkKWKmUJEnqwEcKNfwUJEmS1JmVSkmSpC5cqAOYVEqSJHUzw4lfcPpbkiRJU8BKpSRJUgfxkULAJCuVSe6c6gCSHJ7kkHb/w0n2muL+lyZZmOTyJN9LMqs9fn2SJ49q+3+T/HOSPZMsSXJJkquSfHIqY+oZ75FJbkmy3qjj/5nkFUlenOSw6RhbkiRpKq1S099V9YGq+vEUd3tPVc2uqm2BW4E3t8ePA/YfaZTmeQD7Ase3h35SVTsCOwIvTLJblyCSPKgqXFV3AacCL+lptx6wO/D9qppXVUd2GVeSJE2zGTOmZxsynSNO8qIkP2urej9O8tj2+OFJjk5yZpLrkryt55r3JflFkh8DT+45fkySfdv9G5J8KMnFSS5LsnV7/NFJftQe/2KSXyXZqM9wzwM2afePpSepBPYAbqiqX/VeUFX3AAt7rhvrM9igrS4uSnJ+ku17PoO5SU4FvjrO5aPjeCnww6q6O8mBST7b9vXyttp6aZKzJ4jlwDaW77XV2LckeVf753N+kg3Gue7gJAuSLDjxm18fr3tJkjRaMj3bkJmKNPinwDPaqt5xwD/3nNsaeB6wM/DBJGsk2YkmidoR+Dvg6RP0/fuqehrweeCQ9tgHgdPb4ycBT+gnyCQzgecA8wCqahGwLMkObZP9aRK80detD2wJjJvIAR8CLqmq7YH38sAEcidgn6p61TjX/hDYKcmGE8UBfAB4XlXtALx4glgAtgVeRfO5HwHc3f75nAe8dqwLqmpuVc2pqjl/96rXLKd7SZKkB5qKpPLxwPwklwGHAtv0nDu5qu6tqt8DvwMeCzwTOKmq7q6q/6FN8sZxYvvzImCzdn93muSVqvohcNty4lsryULgFmAD4Ec9544F9m+npvcBTug598wki4Df0kxF/3aCMXYHvtbGdDqwYc99kvPaaueYquqPNJ/Bvm3FdTbNlPho5wDHJHkDMHOCWADOqKo7qupmYAnwvfb4Zfz5c5QkSVPBSiUwNUnlZ4DPVtV2wBuBNXvO3duzv5Q/rzavPvseub732sl+yvdU1WzgicDD+fM9ldAkla8A9gIWVdXves79pK08bgf8Y5LZE4wxVkwj7/GuPmIcmQLfF/huVd33oM6q3gS8H9gUWNhT2RxL7+e+rOf1MlzxL0mSpsFUJJXrATe1+wf00f5s4KVJ1kqyLvCiSY73U5pEkCR7A+v3c1FVLQHeBhySZI322C9pKphHMvaUM1V1NfAx4N0TdH828Oo2pj1ppu3/p5+4WmfQTLG/ebw4kmxRVT+rqg8Av6dJLiVJ0oBlRqZlGzaTTSrXTnJjz/Yu4HDghCQ/oUl2JlRVF9OssF4IfAf4ySRj+BCwd5KLgRcAvwHu6OfCqroEuJQHLow5lubez5MmuPQLwB5JNh/n/OHAnHa6/Ej6S65741pG81lsyPj3bn6iXbB0edvm0smMIUmSpklmTM82ZFLV70z0qiHJI4ClVXV/kl2Az7fT25oiF91w03D9pRgSm27YV1F9lbF02bJBh9C3O//2ZYMOYVLW+8FEv8OuWs6/+oZBh9C3L5z600GHMCnP2W6rQYfQt7/eZstBhzApczbfZKWW+e676TfT8v/mGptsPFTlymG8v+4JwLfa50r+EXjDgOORJEkPZUO4qGY6DF1SWVXX0DyO6E/aRSunjdH8OVV1y1SNneQg4O2jDp9TVW8eq/1KDIInAAAgAElEQVSoa58HfHzU4eur6qUrGMuU9idJktTF0CWVY2kTx2mfAq+qLwNfXsFr5wPzpzCWKe1PkiStoCFcVDMdhu8uUEmSJK1yVotKpSRJ0qBkCFdqTweTSkmSpC6c/gac/pYkSdIUsFIpSZLUwT1rPmJa+l13WnqdPlYqJUmS1JlJpSRJkjozqZQkSVJnJpWSJEnqzKRSkiRJnZlUSpIkqTMfKaQHqRp0BP27/a57Bh1C39Zda81BhzApy5YtG3QIfdvgBycNOoRJWfKClw46hL7tNOdpgw6hb18CXj1rk0GH0bcdNnv8oEPo2+133cOaa5gyaGJWKiVJq4VhSiiHjQml+mFSKUmSpM5MKiVJktSZSaUkSdKQSvL8JL9Icm2Sw8Y4/64kVyZZlOS0JE/sObc0ycJ2m9c1Fm+SkCRJGkJJZgKfA54L3AhcmGReVV3Z0+wSYE5V3Z3kH4F/BfZrz91TVbOnKh4rlZIkScNpZ+Daqrquqv4IHAfs09ugqs6oqrvbl+cD0/bYAZNKSZKk4bQJ8Oue1ze2x8bzOuAHPa/XTLIgyflJXtI1GKe/JUmSVkFJDgYO7jk0t6rm9jYZ47Ixnzad5DXAHOCvew4/oaoWJ/lL4PQkl1XVL1c0XpNKSZKkVVCbQM6doMmNwKY9rx8PLB7dKMlewPuAv66qe3v6X9z+vC7JmcCOgEmlJEnSINw3c41BDX0hsGWSzYGbgP2BV/U2SLIj8EXg+VX1u57j6wN3V9W9STYCdqNZxLPCTColSZKGUFXdn+QtwHxgJnB0VV2R5MPAgqqaB3wCWAc4IQnAf1XVi4GnAF9Msoxmjc2Ro1aNT5pJpSRJUgc15l2MK2vsOgU4ZdSxD/Ts7zXOdecC201lLCaVkiRJHSwbZFa5CvGRQpIkSerMSqUkSVIHZaUSsFIpSZKkKWClUpIkqQMrlQ2TSkmSpA5cqNNYLae/k8xK8k+DjmN5krwjydqDjkOSJKmr1TKpBGYBA08q05joM34HMKmkMonVZUmSViFV07MNm9U1qTwS2CLJwiSfSHJokguTLEryIYAkmyW5KslRSS5P8o0keyU5J8k1SXZu2x2e5GtJTm+Pv2FkkAn6/XmS/wAuBjZN8vkkC5Jc0dPubcDjgDOSnNEeu7On732THNPuH5Pk39p2H0/yyCRHt2NfkmSf8T6IJDOTfDLJZW2cb53ST1qSJInV957Kw4Btq2p2kr2BfYGdgQDzkuwB/BfwJODlwME035/5KmB34MXAe4GXtP1tDzwDeCRwSZKTgW2BLcfp98nAQVX1TwBJ3ldVtyaZCZyWZPuq+nSSdwHPqqrf9/GetgL2qqqlST4KnF5V/5BkFnBBkh9X1V1jXHcwsDmwY/t1Thv0+yFKkqTlc6FOY3WtVPbau90uoakcbk2TDAJcX1WXVdUy4ArgtGr+ZlwGbNbTx3er6p42+TuDJpGcqN9fVdX5Pde/IsnFbdttgKeuwPs4oaqW9rynw5IsBM4E1gSeMM51ewFfqKr7Aarq1rEaJTm4raYuOPHYr69AeJIk6aFsda1U9grwsar64gMOJpsB9/YcWtbzehkP/GxG/wpSy+n3rp7XmwOHAE+vqtvaKe01x4m1d5zRbXqrkAFeVlW/GKefB4Q0RvwPHrhqLjAXYMH1N/krlyRJfVq2/P9mHxJW10rlHcC67f584B+SrAOQZJMkj5lkf/skWTPJhsCeNFPl/fb7KJqEcEmSxwIvGCdOgP9O8pR2cc9LJ4hnPvDWJGnH3nGCtqcCbxpZ4OP0tyRJU6uqpmUbNqtlpbKqbmkX3FwO/AD4JnBem4PdCbwGWDpBF6NdAJxMM8X8kapaDCxO8pTl9VtVlya5hGZ6/TrgnJ7Tc4EfJPlNVT2L5l7Q7wO/Bi4H1hknno8A/xdY1CaWNwAvHKftUTT3Yy5Kch/wJeCz/b91SZKk5cswZsIrU5LDgTur6pODjmVlGabp79vvumfQIfTtCY9ef9AhTMqyZcsGHULfNljnkYMOYVKWvGCiiYhVy9pznjboEPr26lmbDDqESXn/vs8fdAh9W3ON4apB7f7kzbIyx7vhltun5f/NzTactVLfR1er6/S3JEmSVqLh+tVjAKrq8EHH0I8kzwM+Purw9VU1PCURSZKG0LJlQzPBN61MKlcTVTWfZgGPJElaibyTsOH0tyRJkjqzUilJktSBi54bViolSZLUmZVKSZKkDvxGnYZJpSRJUgdOfzec/pYkSVJnViolSZI6sFLZsFIpSZKkzqxUSpIkdeAX6jSsVEqSJKkzK5WSJEkdeE9lw6RSkiSpA5PKhtPfkiRJ6sxKpR7kmt/8btAh9O1Ra6816BD6dvOSOwcdwqQM02/eVy++edAhTMpOc5426BD6dveCiwcdQt+O+e67Bx3CpFx6/U2DDqFv6671iEGHsEpbNkT/Xk4nK5WSJEnqzEqlJElSB1YqGyaVkiRJHQzT7ULTyelvSZIkdWalUpIkqQOnvxtWKiVJktSZlUpJkqQOLFQ2TColSZI6cKFOw+lvSZIkdWalUpIkqQMX6jSsVEqSJKkzK5WSJEkdeE9lw0qlJEmSOrNSKUmS1IGFyoZJpSRJUgcu1Gk4/S1JkqTOrFRKkiR14EKdhpVKSZIkdWalUpIkqQPvqWxYqRxDksOTHDLgGB6X5Nt9tHvvyohHkiSNbVnVtGzDxqRyFVVVi6tq3z6amlRKkqSBM6lsJXlfkl8k+THw5PbYFkl+mOSiJD9JsnWSmUmuS2NWkmVJ9mjb/yTJk8bp//AkX0tyepJrkryhPZ4kn0hyeZLLkuzXHt8syeXt/oFJTmxjuSbJv7bHjwTWSrIwyTcmeG+vTbIoyaVJvjZOm4OTLEiy4LTvfmfFP0hJkh5iqmpatmHjPZVAkp2A/YEdaT6Ti4GLgLnAm6rqmiR/BfxHVT07ydXAU4HN23bPTPIz4PFVde0EQ20PPAN4JHBJkpOBXYDZwA7ARsCFSc4e49rZbXz3Ar9I8pmqOizJW6pq9gTvbRvgfcBuVfX7JBuM1a6q5rbvl2PPvWT4/iZLkqSBMqlsPBM4qaruBkgyD1gT2BU4IclIu0e0P38C7EGTVH4MeANwFnDhcsb5blXdA9yT5AxgZ2B34NiqWgr8d5KzgKcDi0Zde1pVLWnjuxJ4IvDrPt7bs4FvV9XvAarq1j6ukSRJfRrGquJ0MKn8s9F/I2YAt49TBfwJ8CbgccAHgEOBPYGxKowTjVFAxmo4hnt79pfS/59dxhhXkiRNkWX+Lwt4T+WIs4GXJlkrybrAi4C7geuTvBz+dO/jDm37n9FUMZdV1R+AhcAbaZLNieyTZM0kG9IkoRe2Y+/X3qv5aJoK6AWTiP2+JGtMcP404BXtmIw3/S1JktSFSSVQVRcDx9Mkh9/hz8nhq4HXJbkUuALYp21/L83U8/ltu58A6wKXLWeoC4CT2+s+UlWLgZNoprovBU4H/rmqfjuJ8OcCi8ZbqFNVVwBHAGe17+PfJtG3JElaDhfqNDKMQQ+jJIcDd1bVJwcdy/IM00KdR6291qBD6NusIYoVhuseoVvvvHvQIUzKTl8b8yEMq6S7F1w86BD6tsZ3jx90CJNy6fU3DTqEvv3F+o8adAiT8ldbbNrvrWVT4keXXTMt/2A+d7stV+r76Mp7KiVJkjoYpl/Cp5NJ5RRLchDw9lGHz6mqN0/zuBvS3D852nOq6pbpHFuSJMmkcopV1ZeBLw9g3FtonmUpSZJWomU+ZAUwqZQkSerE6e+Gq78lSZLUmZVKSZKkDnz4ecNKpSRJ0pBK8vwkv0hybZLDxjj/iCTHt+d/lmSznnPvaY//IsnzusZipVKSJKmDZQMqVSaZCXwOeC5wI3BhknlVdWVPs9cBt1XVk5LsD3yc5pv8ngrsD2xD87XTP06yVVUtXdF4rFRKkiR1MMBv1NkZuLaqrquqPwLH0X77X499gK+0+98GnpMk7fHjqureqroeuLbtb4WZVEqSJK2CkhycZEHPdvCoJpvQfG30iBvbY2O2qar7gSXAhn1eOylOf0uSJHUwXY8Uqqq5wNwJmoz1NY6jgxmvTT/XToqVSkmSpOF0I7Bpz+vHA4vHa5PkYcB6wK19XjspJpWSJEkdLKOmZevDhcCWSTZP8nCahTfzRrWZBxzQ7u8LnF5NaXUesH+7OnxzYEvggi6fg9PfkiRJHQzqG3Wq6v4kbwHmAzOBo6vqiiQfBhZU1Tzg/wFfS3ItTYVy//baK5J8C7gSuB94c5eV32BSqTFcfN2vl99oFbHdEzvdU7xSLbrhpkGHMCkzZox1u82q6awrrh10CJNy96zh+Xt7zHffPegQJuW+ffYbdAh9W/Kxjw46hL4tuecPbLjuIwcdhsZQVacAp4w69oGe/T8ALx/n2iOAI6YqFpNKSdJqYZgSymFjQjkxv/q74T2VkiRJ6sxKpSRJUgfLLFUCViolSZI0BaxUSpIkdTCo1d+rGpNKSZKkDkwqG05/S5IkqTMrlZIkSR24UKdhpVKSJEmdWamUJEnqwEplw6RSkiSpAxfqNJz+liRJUmdWKiVJkjpYZqESsFIpSZKkKWClUpIkqQPvqWyYVEqSJHVgUtlw+luSJEmdWamUJEnqwOdUNqxUSpIkqbPVIqlMcmaSOctp844ka/e8PiXJrCkaf7Mkl09FX2P0/eEke41xfM8k32/3X5zksHb/JUmeOh2xSJKkB6uanm3YDE1SmUaXeN8B/CmprKq/qarbu0c2varqA1X14+W0mVdVR7YvXwKYVEqSpJVqlU4q2wrgz5P8B3Ax8PdJzktycZITkqwzxjWfT7IgyRVJPtQeexvwOOCMJGe0x25IslGSjyf5p57rD0/yv9r9Q5NcmGTRSF8TmJnkS+24pyZZq+3jT1XUdrwb2v0Dk/xnku8luT7JW5K8K8klSc5PskHb7pgk+7b7z09yVZKfAn/XE/OBST6bZFfgxcAnkixMskWSi3vabZnkokn9IUiSpAlV1bRsw2aVTipbTwa+CjwXeB2wV1U9DVgAvGuM9u+rqjnA9sBfJ9m+qj4NLAaeVVXPGtX+OGC/ntevAE5IsjewJbAzMBvYKckeE8S5JfC5qtoGuB14WR/vbVvgVe0YRwB3V9WOwHnAa3sbJlkT+BLwIuCZwF+M7qyqzgXmAYdW1eyq+iWwJMnstslBwDFjBZLk4DYZX3Dp6fP7CF2SJEGzUGc6tmEzDEnlr6rqfOAZNNO65yRZCBwAPHGM9q9oq3OXANuwnKngqroEeEySxyXZAbitqv4L2LvdLqGpkm5NkziO5/qqWtjuXwRs1sd7O6Oq7qiqm4ElwPfa45eNcf3W7RjXVPPry9f76B/gKOCgJDNpkudvjtWoquZW1ZyqmrPDs5/XZ9eSJEmNYXik0F3tzwA/qqpXjtcwyebAIcDTq+q2JMcAa/YxxreBfWmqf8f1jPexqvpin3He27O/FFir3b+fPyfvo2PpvWZZz+tljP1nsyK/tnwH+CBwOnBRVd2yAn1IkqRxDONU9XQYhkrliPOB3ZI8CSDJ2km2GtXmUTRJ6JIkjwVe0HPuDmDdcfo+DtifJrH8dntsPvAPI/dtJtkkyWNWIO4bgJ3a/X1X4PoRVwGbJ9mifT1ecv2A91lVf6B5L58HvtxhfEmSpHENTVLZThEfCBybZBFNkrn1qDaX0kxXXwEcDZzTc3ou8IORhTqjrruCJhG7qap+0x47lWaq+Lwkl9Ekm+MlpRP5JPCPSc4FNlqB60di/ANwMHByu1DnV+M0PQ44tF3wM5KAfoOmynnqio4vSZLG5j2VjViyXf0lOQRYr6r+dz/tD/36vKH5S7HdEzcZdAh9u+mWVf4JVg8wY0YGHULfzrri2kGHMCl33/vHQYfQt2Pe+ppBh9C3+/bZb/mNViHnfeyjgw6hbxuu+8hBhzApfzt765X6D9jnTj1nWv7ffPPeuw3PP8QMxz2V6iDJScAWwLMHHYskSVp9mVROQpINgdPGOPWcVXUBTFW9dNAxSJK0OnPWt2FSOQlt4jh7uQ0lSZIeYkwqJUmSOrBQ2TCplCRJ6mAYV2pPh6F5pJAkSZJWXVYqJUmSOnChTsNKpSRJkjqzUilJktSBlcqGSaUkSVIHLtRpOP0tSZKkzqxUSpIkdWCdsmGlUpIkSZ1ZqZQkSerAeyobViolSZLUmZVKPcjuW28x6BD6dvvd9ww6hL5dvfh3gw5htfWc7bYadAiTssNmjx90CH279PqbBh1C35Z87KODDmFSdnnPewcdQt9OP/zwQYewSvORQg2TSkmSpA6WLTOpBKe/JUmSNAWsVEqSJHXg9HfDSqUkSZI6s1IpSZLUgY8UaphUSpIkdWBK2XD6W5IkSZ1ZqZQkSerAhToNK5WSJEnqzEqlJElSBy7UaZhUSpIkdeD0d8Ppb0mSJHVmpVKSJKkDp78bViolSZLUmZVKSZKkDixUNqxUSpIkqTMrlZIkSR24+rthUilJktSBC3UaTn9LkiSpM5PKSUiyZ5Lv9+zv2nPuTUleO8n+zm1/bpbkVVMY57lT1ZckSZrYsqpp2YbNwJLKJDNX4ljTMc2/J/CnpLKqvlBVX51MB1U1cv1mwKSSyok+v55+JUmSVoppSSrbyttVSb6SZFGSbydZO8kNST6Q5KfAy5NskeSHSS5K8pMkW7fXvzzJ5UkuTXJ2e2ybJBckWdj2uWU7zuU94x6S5PB2/8wkH01yFvD2JI9O8p0kF7bbbhPEv3OSc5Nc0v588uj3B7wJeGcbzzOTHJ7kkJ6xP5Xk7CQ/T/L0JCcmuSbJv/T0c2e7eyTwzLavdyaZmeQTbZyLkryxbb9nkjOSfBO4bIL47+xpf2b7+V+V5BtJMs41BydZkGTB/BNPGK9rSZI0SlVNyzZspnOhzpOB11XVOUmOBv6pPf6HqtodIMlpwJuq6pokfwX8B/Bs4APA86rqpiSz2uveBPx7VX0jycOBmcBjlxPDrKr663asbwKfqqqfJnkCMB94yjjXXQXsUVX3J9kL+CjwspGTVXVDki8Ad1bVJ9v+nzOqjz9W1R5J3g58F9gJuBX4ZZJPVdUtPW0PAw6pqhe2fR0MLKmqpyd5BHBOklPbtjsD21bV9ct57yN2BLYBFgPnALsBPx3dqKrmAnMBvrvgiuH7myxJ0oCsiglgkg2A42lmQ28AXlFVt41qMxv4PPAoYClwRFUd3547BvhrYEnb/MCqWjjRmNOZVP66qs5p978OvK3dHwl2HZrp4xN6imePaH+eAxyT5FvAie2x84D3JXk8cGKbiC4vhuN79vcCntpzzaOSrFtVd4xx3XrAV5JsCRSwxvIGGsO89udlwBVV9RuAJNcBmwK3jHchsDewfZJ9e+LZEvgjcMEkEkra9je2Yy+k+cv1oKRSkiStVg4DTquqI5Mc1r5+96g2dwOvbXOqxwEXJZlfVbe35w+tqm/3O+B0JpWj0/aR13e1P2cAt1fV7AddWPWmtnL5t8DCJLOr6ptJftYem5/k9cDVPHAKf81RXd3Vsz8D2KWq7ukj9o8AZ1TVS9up7jP7uGa0e9ufy3r2R14v73MP8Naqmv+Ag8mePPA9TSYOaH4L8TFSkiRNoWWrXqESYB+a9R8AX6HJZR6QVFbV1T37i5P8Dng0cDsrYDoX6jwhyS7t/isZVR2rqv8Brk/ycoA0dmj3t6iqn1XVB4DfA5sm+Uvguqr6NE0VcHvgv4HHJNmwnSZ+4QTxnAq8ZeRFW/Idz3rATe3+geO0uQNYd4I+JmN0X/OBf0yyBkCSrZI8corGkiRJQ6B3vUO7HTyJyx87Mkva/nzMcsbaGXg48Muew0e0azs+1eZZE5rOpPLnwAFJFgEb0MzZj/Zq4HVJLgWuoMmqAT6R5LJ2Ec7ZwKXAfsDl7RTu1sBXq+o+4MPAz4Dv09wLOZ63AXPaD+dKmns0x/OvwMeSnENz7+ZYvge8dGShzgR99WMRcH+ahUnvBI4CrgQubj+DL2KFUZKkVdJ0LdSpqrlVNadnm9s7bpIfp1nYPHrbZ7xYx5JkY+BrwEFVtaw9/B6afOvpNHnc6KnzB/czHTeXtlPG36+qbae8c027YVqoc/vd/dzNsGo48/JrBh3CamvbJ2w86BAmZYfNHj/oEPp21x/uXX6jVcSSe/4w6BAmZZf3vHfQIfTt9MMPH3QIk/LGvXZZ7qKLqfTOr/zntPy/+akDXrLC7yPJL4A9q+o3bdJ4ZlU9eYx2j6KZGv9YVY35+Jf29rs/LSgejw8/lyRJWv3MAw5o9w+geRLNA7RP0zmJZvb3hFHnNm5/BngJcPno60eblinVqroBWOWrlEkOAt4+6vA5Vf+/vTsPk7sq0///vlkDSNhFRRZBRFllX0VZRx0EEWRRRFFwB8EfuHydEUQUt2EGcQZUNCCgCCIiqBBA9lUCIUEQFxZFUESBRBYhcP/+OKfS1Z2q7grV3ed80s/ruurqrk+nO3c63VVPneU5/kiJPPND0grA5R0+tNOQdkUhhBBCGEOVnn7zJeAcSe8D/gi09rBsRmrneDCwD7A9sIKk9+TPa7UOOkvSSqTNw9MZftkgMMHX6dmeAkwpneOFyIXjcJuNQgghhDBB5TphaA9tbN8CHJzfP5PU9rHT5+84v3/nhC4qQwghhBD6VWPz8xJiTWUIIYQQQuhbjFSGEEIIIfSh0ubn4y6KyhBCCCGEPjw/t7XjxBbT3yGEEEIIoW8xUhlCCCGE0IfYp5PESGUIIYQQQuhbjFSGEEIIIfQhWgolUVSGEEIIIfSh0hN1xl1Mf4cQQgghhL7FSGUIIYQQQh9i+juJojLMY4nFFi0doWdbX/nL0hF6tsvBB5WOsMB68NHHS0eYL08/82zpCD1beonFS0fo2UuWm8wjs58oHaNnvzzmmNIR5suOTcq78yWlE0xIUVSGEEJYIDSpoGyaRhWUBcRIZRJFZQghhBBCH+KYxiQ26oQQQgghhL7FSGUIIYQQQh9i+juJkcoQQgghhNC3GKkMIYQQQujD88RIJcRIZQghhBBCGAUxUhlCCCGE0IdYU5lEURlCCCGE0Ifno6cQENPfIYQQQghhFMRIZQghhBBCH2L6O4mRyhBCCCGE0LcYqQwhhBBC6EMsqUyiqAwhhBBC6ENMfycx/R1CCCGEEPoWI5UhhBBCCH1wnKgDxEhlCCGEEEIYBTFSGUIIIYTQh+djTSVQ6UilpOckTZd0h6QLJS1bOhOApH/LuaZL+qeku/P73xuFr31K/lp3Snqq7e/Zs8+vO0nSI/3mCyGEEEJntsfk1jS1jlQ+Zfu1AJJOBz4CfKFsJLB9CXAJgKQrgSNt3zL0z0laxPac+fzaH8yf+0rgR61//2h87RBCCCGEsVblSOUQNwCrAEh6g6SLWh+Q9A1J78nv3yfpc5JulTRT0qvz9ZmSllXyd0kH5utnSNpZ0hqSrsmfd6ukbdo+vkfb33WWpN27hZR0sKSzc75f5K/9k7aPnyLpgPz+5pKukjRN0i8krTzcN0DSjZKOk3Q18CFJb5N0Ux7JvFjSivnPTc65Z0qaIWm3IV9nZUk3S9qll298CCGEEEb2vMfm1jRVF5WSFgZ2An7a46c8YnsT4GTgyHztOmBbYD3gHuB1+fpWwI3Aw8Au+fP2Bb6eP34qcFDOsQywDfDzEf7+rYF32e5atElaHDgR2Mv2psCZwOd7+LctZXt72ycBVwJb5dHMi4Aj8p/5PPCg7Q2AjfK/vfX3viz/2aNsX9rD3xdCCCGE0LNai8olJE0H/g4sD/RaBP04v50GrJHfvwbYPt9OBjaQtArwD9v/BBYFvi1pJnAusC6A7auAV0p6MbA/cF4P085TbT86wp95DanAvSz/Gz8FrNrDv+3stvdXBy7NmT+Wvx7AzsA3cn63ZZkEXAYclv9d85D0fkm3SLrl5z/6YQ9xQgghhACxprKl1qKytaZydWAx0ppKgDkMzjxpyOf9K799joH1oleTRidfRxrh+xuwN6nYhDTK91fSyN5m+e9rOQN4J2nEckoPuZ9oe79bVgEzbL823zaw/ab5/NonA1/JI5KHDvnanX4KnwFmAF1HUG1/y/Zmtjd789779hAnhBBCCGFArUUlALYfBw4DjpS0KHA/sK6kxfOU9E49fI0/ASsCa9u+B7iWNDXeKiqXAR6y/TzwLmDhtk8/DTg8f51fz2f8+4H1JC0maTlgx3z9TmAVSVsA5I+v1+2LdLEM8GdJAt7ddn0q8NH8dZX/XoDWv20LSYfP598VQgghhGHESGVSdVEJYPs24HZgv1wgnkMadTsLuK3HL3MT8Nv8/jWkjT/X5vv/B7xb0o3Aq2gbEbT9V+AuehulHJr7XuAnwEzge8Ct+fq/SCOlJ0i6Pf8btpzPL/9Z0vrIq4AHh1xfVdIdwHTSOtBWnmeBtwO7S3rf/P57QgghhNDZ8/aY3JpGTayEx4ukJUlF4SZ51HRCmDrjt435oXjt1EtKR+jZnIMPKh1hgfXgo8369Xz6mWdLR+jZoossPPIfqsQjs58Y+Q9V5IFHRlqCX48djzmmdIT5sva1l2g8/77dv3LqmDxv/vQTB4/rv6NftfapLE7SzsB3gRMmUkEZQgghhPnTxFHFsRBFZRe2LwNWK50jhBBCCKEJoqgMIYQQQuhDLCVMoqgMIYQQQuhD1JRJ9bu/QwghhBBC/WKkMoQQQgihD7FRJ4mRyhBCCCGE0LcYqQwhhBBC6ENs1EmiqAwhhBBC6ENMfycx/R1CCCGEEPoWI5UhhBBCCH2I6e8kRipDCCGEEELfYqQyhBBCCKEPMVCZxEhlCCGEEELoW4xUhhBCCCH0IXZ/J1FUhhBCCCH0ITbqJIpvRBgPkt5v+1ulc/SqSXmblBWalbdJWaFZeZuUFZqVt0lZoXl5Q3expjKMl/eXDjCfmpS3SVmhWXmblBWalbdJWfYxHM4AACAASURBVKFZeZuUFZqXN3QRRWUIIYQQQuhbFJUhhBBCCKFvUVSG8dK09TJNytukrNCsvE3KCs3K26Ss0Ky8TcoKzcsbuoiNOiGEEEIIoW8xUhlCCCGEEPoWRWUIIYQQQuhbFJUhhBBCCKFvUVSG0ECSFu/lWgghhDBeoqgMY0bSRyUtVzpHryR9TNJkJd+RdKukXUvn6uKGHq9VQdIZvVyrhaRlJR0m6QRJX2/dSudqJ+kjkpZtu7+cpA+XzLSgkXRL/j436XFsG0nvkHRg61Y6UyeSdupw7aASWcLoiaIyjKWXAL+SdI6kN0pS6UAjeK/tWcCuwErAQcCXykYaTNJLJG0KLCFpY0mb5NsbgCULxxvOeu13JC0MbFooSy9+DqwBzASmtd1qcojtx1p3bD8KHFIwz4gknd6hEP5uyUwj2A94Gelx7GxJ/1bz41h+ofY1YDtg83zbrGio7r4oabvWHUkfA95ZME8YBdFSKIyp/AC8K6lA2ww4B/iO7T8UDdaBpBm2N5R0InCl7fMl3WZ749LZWiS9G3gP6Xv5K6D1BDcLON32jwtF60jSp4H/BywBPNm6DDwDfNv2p0plG46kW21vUjrHcCTNADZyfhDPhfoM2+sN/5nldPp9qu13rBNJCwG7AScDzwPfBU60/Y+iwYaQdBewrhvwxC5pZeBC4FDSc8QWwNttP100WOjLIqUDhAWbbUv6C/AXYA6wHPAjSZfa/kTZdPOYJmkq8Arg05KWJj2BVMP26cDpkg6wfWb7xyQtXyhWV7aPB46XdLztT5fOMx/OkHQIcBHwr9bFyoqIS4BzJJ0CGPggcHHZSCNaSNJyeVS19TNb9fOQpA1JL4rfDJwHnEUaCfwl8NqC0Tq5gzRD9FDpICOx/VdJewCXAjOAt9p+rnCs0Keqf5lDs0k6DHg38AhwKnCU7Wfzq/7fAbUVle8jPUncY/tJSSuQnkxqtL+ks23PgTQtDvyMeqeUn22/k38GzrBd63TXM8BXgc+QCjby2zWLJZrXJ4H3Ax8ijf5OJf2e1ey/gOsl/Yj0/dwH+ELZSN1JmgY8BnwH+JTt1guMmyRtWy5ZVysCd0q6mcEvhnYvF2kwSX9j4HcK0izGS4GHJGH7xWWShdEQ099hzEg6ljTVfX+Hj73G9l0FYnWVp+rfCaxp+1hJqwEvsX1z4WjzyKNo/w7sBawK/BQ40vbUosG6kHQacLft4/Mu9XOBW20fUzRYF5L+AGxp+5HSWV4oSefZ3qt0jqEkrQvsSCqEL7d9Z9vH5o5i1kDSmrbvGXLtFbbvLZVpOJJe3+m67avGO0s3I3WpaCvcQwNFURnGnKQXA5Na923/sWCcriS11kvtaPs1ecfnVNubF47WkaSPAG8kbSj5gO3ryybqLhfsZ5E2vuwA/ML2f5dN1Z2knwL72X5yxD9cqSasVRyqtrWsnfJImma71hmBxpD0ZuA624/n+8sCW9v+RdlkoR8x/R3GjKS3ACeQdk8+DKwO3MWQncAV2dL2JpJug7SbVtJipUO1k/Tx9rukUcrpwFaStrJ9QplknUlqf0I+EfgmcB1wlaRNbN9aJtmIngOmS7qCwdOIh5WLNN+aOGJQxc5qSa8mPU4tI+ltbR+aTNsL5FpIutb2dpJmM/j/XaSl7ZMLRRvOF23PXZNq+zFJXwSiqGywKCrDWDoO2Aq4zPbGknYA9i+caTjP5h20rd20K1HZRh1g6SH3z+9yvRb/NeT+o8C6+bpJ06A1+km+hfFVSyG8Dmm397LAW9quz6bCtk22t8tvh30cqGx5QaeWhlGTNFxMf4cxI+kW25tJuh3Y2Pbzkm62vUXpbJ1IeiewL2mzy2nA3sB/2D63ZK4XQtJJtg8tnaNXkt6dd7Y3Qq3rFdvF9Hf/JG1tu9pDBeZXTd9fSd8D/gx8g/Ri4lBgVdsHFA0W+hKvCsJYekzSi4CrgbMkPUxqK1Ql22fl3Z47kaaN3lrbZqL5UOPO1OF8DGhMUUklu8AlLQGsZvvuDh/+5HjnGQW1TH9/wvZXgHdImmd2pWHLINpV8f3NPgJ8ntQaq9W94ENFE4W+RVEZxtIewNPAEaRd1csAxxZNNLIVgSdtT5G0Us07PRcwNT3Z9aL4FE9es/w1YDHgFZJeCxzbah9TYycASWfYftcw1+Y5uq+Q1ovJW4qmGH3Ff25bbM8GDi+dI4yuKCrDmLH9RNvd6kehJB1NOqlmHWAKsChwJs0b9Wuiap7sGuQY0ikkVwLYni5pjXJxejLscZ21NJe3fWF+W/3jVlPlxvdHkH4m2ruDvLlYqNC3KCrDqOuwA3GQSnciAuwJbAzcCmD7wXyqThM1beQv8s6/ObYfr/go6rnaj+uUNKt1mdRk/lvFgnUh6UKGfwyrppn4fKrph+VM0oEN+wCtgzL+XDRR6FsUlWHUtXYg5ubnfwHOID2YvZN6dykDPJOPlWzt/l6qdKA+nFg6QLtOywiGXLuuQKye5H6lq9qe0Xa5hvWKd0h6B7CwpLVJT8xV9ipt4HGdX8tv30Y69rB1JOr+wH0lAo0kn1I1w/b6w/yxWpYXQDpY4n8lfcD2JfmI3CtKhwr9id3fYcxIusn2liNdq4WkI4G1gV2A44H3At+3fVLRYG2aOoLStCbSkq4Edie98J4O/A24yvbHh/u88SRpSdIxkrvmS5cAx9l+ulyqkUlahdSzdu6ghu2ryyXqTtLVtrcf6VotJJ0FfLrWAybaSbrR9laSLgW+BDwI/Mx2FZvgwgsTI5VhLD2X2/ScTSqE9ic1la6S7a9J2gWYRVpX+VnblxaONVSjRlCa1kS6zTK2Z0k6GJhi+2hJM0b8rHGS1yJ+zvZRpMKyESR9CdgPuJOBxwKTOkTUaKX2oxolvQJYqXCm4bwU+HU++3vumvZKX2x+OZ+icxRwMukx4RNlI4V+RVEZxtI7SNOwJ5KeOK7L16qTn6Qvsb0zUFshOVfrDF9Jnx8yWnKhpBqfmBvVRLrNIpJeSlrvVV3RZvs5SVWO8o5gT2CdBp3vfARwpaTW+d9rAB8oF2dEnysdoBf58fZlth8jzQRsXThSGCVRVIZRJ+nLtj8JbG57j9J5epGfpJ+UtEzrLNrKNWIExfYFwAUNbCJ9LGk6+Vrbv5K0JvC7wpmGui2fUX4ug0elflwu0ojuIXVVaERRafvivF711fnSb2ouiG1fJWl1YG3bl+UlEguXzjVUfrx9O/C/pbOE0RVrKsOokzQT2AS4qZbTG3oh6RzSsZKXMvhJurpGx5LeSNo1O2gExfYlxUINQ9JXSMd2PgVcDGwEHG77zGE/MXQlaUqHy7b93nEP0yNJ55H+7y+n4jPVJe1o+5dDlmzMVWvhLukQ4P3A8rbXygXxKbZr2qADpNkWYEnS8qj2x9s7i4UKfYuRyjAWLgYeAZbK7UPa21i44pZCP8u36jVtBAXY1fYnJO0JPAC8nbTTs6qisnWSiqST6LAhqqbix/ZBpTO8AD/Nt9q9Hvglg5dstBiosqgknVKzBXATgO3fSXpx2Uhd7ZzfbtN2zUPuh4aJojKMurx54ChJFzRl+hua1eg4T2t9HFjd9iGS1pa0ju2LSmfrYtH89s3AD2z/o9L+io05SSWPVHYqfKsdqWzK75jto/PbphXu/7L9TOt3S9IiVHawgKQ9bF9gO9ZRLoCiqAxjxvYeklYGNs+XbrL9t5KZhpOn7Yc+AD9OKjCOs/338U/V1RRgGgML3B8gra2rtai8UNJvSNPfH5a0EukIz9osLmlSQ4qf9v/rSaRNMA8WyjIsSefY3qfL7xi2NywQa0R5d/KBpOUl7S2QqhmxHuIqSa0m87sAHwYuLJxpqKOBC0qHCGMj1lSGMZMXYn+NdIycgNcBR9n+Uclc3eR1f88B38+X9iPlfhzYznanqbAiJN1iezNJt9neOF+73fZGpbN1k5uIz8qL9JcEJtv+S+lc7SSdTzqW82LgB8BU29W2wWqXm19fZnvH0lmGkvRS2w/lTSTzsH3/eGfqhaTrgRuBmcDzreu1vujIPwPvI/UuFamjxbfLphqsU8/asOCIojKMGUm3A7vYfjjfX4n0pFdl4SPpOtvbdromaabtDUplGyo/2e0EXGd7E0lrkaaVtygcrSNJB3a6bvt7451lJJImk0b99iNtKrmA9L2tsWXTXJLWITWPfmXpLMMZMntxc+vxoUZNK4Akfcz2iSNdK0nSk8BvOn2ItOa+Md/vMK+Y/g5jaaEhTxh/BxYqFaYHL5K0pe2bACRtAbwof2xOuVgdHU0aTVs1n6KxLfCeoomGt3nb+5NIBfGtQHVFpe1ZwOnA6ZJWAPYGTpK0vO1Vy6YbIGk2g6eS/0Idx0d2JWkf4KsMzF6cJKna2QvgjLyj+iIG71b/R7lIw3o38x7R+p4O10q6l7RRLyyAoqgMY+liSZeQphEB9gV+XjDPSA4GviupVUjOBg5WOgP8+HKx5mX7Ukm3klogCfiY7UcKx+rK9qHt9yUtQzoTvlp5uv5tpJ/b5YHzyiYazPbSpTO8AJ8h9a8dNHsB1FpUPkMqgj/DQAFvoKqjBCXtTzpY4hW5d2nL0qQX8zV5xvYfSocIYyOKyjBmbB+V+7xtRyp8vmX7/MKxurL9K2CDXPAon/bQck6hWB0pbe98E7Cm7WMlrSZpC9s3l87WoydJ56xXRdLSwFtJx15uQmp/cxxwhStbKyTp8qH9Bztdq0zTZi8+Dryy5hds2fXAQ8CKwH+1XZ8NVHO8aHZj6QBh7ERRGcbEkGMPa+3pNkhe6/VF0vFhb5K0LrC17e8UjtbJ/5E2DuxIOv1lNmkkbfPhPqkUSRcyMNKzELAulRXq2b2kk3ROBi62/WzhPPOQNInUNHrFPJra6s00GXhZsWC9adrsxa9JL4Cqljc63U8Djju0/SGYO0p9HPBy2/+eH2+3sH1ayXyhP1FUhjHRwGMPAU4jteppnfX8W+CHQI1F5ZZ5g85tALYflbRY6VDD+Frb+3OA+20/UCrMMFazPWIRIek823uNR6AOPgAcTiogpzFQVM6i8mPvmjZ7QeoGMV3SFdR9AtC1trfrsM62tfmlxgMnTgPOYmAd8O9Ij7enFcoTRkEUlWEsPQ3MlFT9sYfZirbPkfRpANtzJNXaTubZPBpsmPuq//nhP6Uc21eVztCLXgrKrNiauryT90RJh9o+qVSOPlxPKtaeB35VOMtIfpJvVbO9XX7bpHW2L7b9fUlHAdh+tuLH29CjKCrDWGrMsYfZE3m3b6tQ24rUo7JGXwfOB1aW9AXSDuX/KBtpXpLeRzqH+Kv5/gOkaVoBn7B9csl8fSi+vtL2SZLWJy0lmNR2vbod9S2SDgY+SzoCsbX7+1jb3y2brLNa+1EOJy+JWJXBzdpvLZeoqyckLc/A4+3mpGU8ocGiT2UImaRNgJOA9YE7gJWAvW3XttAdAEmvJrXmAfil7buG+/MlSPoV8MbWaUStZu15XeBU29uXTfjC1NC/UNLRwBtIReXPSRu3rrW9d8lcw5F0N7BN28/DCsD1ttcpm6wzSffS+QSgqnZ/t0j6PKmF0D0MzFy40ob4m5FaHa0H3A6sQnq8nV40WOhLjFSGMdO0B2Tbt0p6PbAOaRTl7ho3arRZEmhNgS9ROEs3Cw053vJcANtPS6o1cy9qOLh8b1Jz9ttsH5Q3mp1aONNIHmDwaNRs4E+FsvRis7b3J5H6Ky5fKEsv9gHWsv1M6SAjsX2LpB2A15B+n+5sQu4wvCgqw1hqxANy3jjQyaskYbu63euSPkv6fp5HekCeIulc28eVTTaPZdrv2P4izD1OboUiiUaQ16qebvuAYf5YDU3Gn7L9vKQ5+RSgh6msf2IHfwZuknQB6cXQHsDNkj4OYPuEkuGGGvKCCOB/JF1LmsKv0R3AsqSfhapJeseQS+vmx9vvd/yE0AhRVIYx06AH5NaZ3i8GtiGt9wLYgXTyR3VFJamP4sa2nwaQ9CXSCTW1FZVTJR1ne+h6z2OBqSUCjSR3LlhJ0mLdRk5s15D9FknLAt8m7QL/J1B7n9I/5FvLBfltlRtM8pKYloVIL5SrzJodD9wm6Q4G71bfvVykrl7X9v4kUnu0aUAUlQ0WRWUYM015QLZ9EICki4B1bT+U77+Uelu03Ed6IH4631+cwU/WtTgKOFXS70nrpiBN2d5COsGoVvcB1+XTSdo7F1Qxkpab3x+fG/SfIuliYHKt639bbH9uuI9LOmno6UuFtTcSn0P6udinTJSenA58GZhJxd0gYKBfZUveYHRamTRhtERRGcZS0x6Q12gVlNlfgVeVCjOCfwG/zu2aDOwCXCvp61BP2ybbTwD7S1qTtCAf0tqpQQWwpPVs/3rcA3b3YL4tRJ0vhCzpJ8Cm+f59ZRONmm1LB2hne4fSGebTI7a/XjrECzSbeh9vQ4+iqAxjpoEPyFe2nfZhYD/girKRujo/31quLJSjJ7bvIe1I7eYM0rGIVWiNqOVjG237n4UjdXKjpM3z8aJhDEhaHNgLWIPBLXqOLZVpBNMkHU86XrR9+ru6lkKSzmfwKVvr0YCeoGF40VIojJl8hvbRQKttzFXAsTWfsCNpTwbyXl35aR9IWpTUAunPQ85UbpRWq6HSOVpy/8czGNhY9ghwYE2jqZLuJHUquI80Rd86PWXDkrn6UUOrpnZ5WcHjpLV+cxtz2/6vrp9UUD75Z6haWwq1n1HfOmXrvkJxwiiJojKMGUnnkXYjthoIvwvYyHa33dZVk3SD7aJn60o6BTjJ9q9z0X4D6clueeBI2z8Y9gtUqsJi4nrgM7avyPffAHzR9jZFg7WRtHqn6/kc6Eaq8MXFHbbXL51jQZM7LPzc9r+VzhJGV0x/h7G01pDzkT8nqcmNbSeN/EfG3OtsfzC/fxDwW9tvlfQS4BekqfvQv6VaBSWA7SslLVUy0FC275e0HbC27Sn5qM4Xlc41HElr5qUQ3Zw4bmF6c72kDWzPLB1kOJIOsH1mqzXTULVsMGvJHRaekTTZ9qzSecLoiaIyjKWnJG1n+1oASdsCTxXO1I8ahvXbW9zswkAz8b+kDcGNVVvT43sk/SdpChzgAODegnnmkU/U2Yw0BT4FWBQ4k8o2uwxxmqRVSGd+Xw1c016w2T6tVLB2kmaSft8XAQ6SdA9pjWKtSwxaL3g6bSqr4XGrk38Ct0uayuAOCx0L49AMUVSGsfQh4PQ8TQvwKPDugnkWBI9J2o3URHpb4H0AkhahwlN1hrSVmkdrA4HtrcYnUc/eC3yOgR6lV5NGhmuyJ7AxqT8pth/MG4uqZXt7SYsBm5OOmPyZpBfZru1QhN16+UOSlrP96FiHGYntb+Z3L7N9XfvH8ov5Gl2Wb2EBEkVlGEt3AV8B1iKd8vA48Fag6l56w6hhKPADwNeBlwCH2/5Lvr4T8LNiqbprbWiYRBpVu530fdwQuAnYrlCuYeVCoWtbpkr6KT6TWws5Z6pqer6TPF3/unxbFrgIuKZoqA7mY13q5VTUtQA4iXnzdLpWjKTTbL/H9ndKZwmjL4rKMJYuAB4jjaT8uXCW0fCu0gFs/xZ4Y4frlwCXtO5L+rTt48czWyettlKSzgbe35rqzLurjyyZrU81jP6cI+mbwLKSDiGNrn67cKaRXEVqfH88aaNGbcse5lcNLzSRtDXpNLCVhqyrnAwsXCZVV7UtHQijKIrKMJZebnueAqhW+QzwL5OOaxQD66cmk965o2C8+fV20hN3LV49ZO3cHZJeWzJQ09n+mqRdgFmkptGftX1p4VgjWYFUkG8PHCbpeeAG2/9ZNtYLVst6xcVIm7QWYfC6ylnA3kUSdbekpI3pUpDX2FMz9C6KyjCWGrFzss1XgLfYvqt0kFFQxQhKm7sknUraSGLSxpcF4ftc2kzSWlrn96tm+7G86WVV4OWk0bVFy6ZqPttXAVflqeXaW0qtQloW0+kxyqQzwENDRVEZRl0Dd062/HUBKSihnhGUloNIG7c+lu9fDZxcLk7fihftkg4GPgv8kpTnJEnH2v5u2WTdSfoDcDdwLXAKcFDDp8CL/xwASLqQ/DvfqQuE7d3HO9Mwfl9jM/YwOqL5eRh13Zoyt9T6SlrSiaQNMD9h8BFnP+76SZWqrYn0gkDS6q2fXUnvKd3+RtLdwDa2/57vrwBcb3udkrmGI2kh28+XzjE/OvUCtX1v/tjytv9RNiFIev1wH88jmVWIx6YFW4xUhlFXa9HYg8nAk8CubdfMQFuZJjm3dIB2ua3JMcDqDD5Dec1SmbrJmx5WIR3T+bCkDYFPkXYsrwrV9FN8AJjddn828KdCWXr1SkknAyvbXj9/b3e3fVzpYJ2M1Au0hoIS6ioae/BJgNwa7edNe5ERhhcjlSE0iKSTGGZq23bXNjglSfoNcATznqH892KhOpD0VVKPwunAK0ktbz4MfBH4pu2nC8YbRNL3gA1IXRYM7AHcDPwW6jtFBUDSVcBRpO/lxvlatUch5hPANgZubcs7o9YlPJLWJm3QW5e2E8AqffF2JrA1cB4wZQFaejShxUhlCJmkl5N6um1LepK+FviY7QeKBhvslvx2W9ITxw/z/beTCrZaPW77F6VD9ODfgY1tPy1pOeBBYEPbvyucq5M/5FvLBfltzQ3Ql7R985B1f3NKhelB03qBTgGOBv4b2IG0lrmKdZ9D2T5A0mRgf2BK/h5PAX5ge/bwnx1qFUVlCAOmAN8nFWiQdihPIR2HWAXbp0Na0wfsYPvZfP8UYGrBaCO5Io8C/pjB61Vrax/yVGs00vajku6utKDE9udKZ3gBHpG0FgObSvYGHiobaVhN6wW6hO3LJSkvQzpG0jWkQrM6tmdJOo/UweBw0ilRR0n6uu2TyqYLL0RMf4eQSZpu+7UjXatB3qSxdWtNVx5Vu7HWTRqSruhw2bXtApX0GGlnesv27fdr2kUraTPgM8y7TrXKqVkASWsC3yK1EnqUdJ76O2teh517ge5KGvG7pOZeoJKuI639/RGpK8CfgS/V+Lgg6S2kIn0t4Azg9LyGeUngLtvDbvgMdYqiMoRM0mXAacAP8qX9SS1PdioWqgtJB5E2vrSKtdcDx7RGMptA0sq2/1o6R7uG7aK9m7Q+cSYwd7NDjQXakFNeII1MLQQ8AXWu/wSQdARwbmVLYLqStDmp/+uywOdJmw+/avvGosE6yGuCT7V9dYeP7WT78gKxQp+iqAwhk7Qa8A3S4nED15PWVFb3JA0g6SXAlvnuTW3ngFdL0jLAXsA7gNfYXqVwpEEkHQn80Hbtu6iRdK3tKs9OHyrvooa0i3pz0vpPAW8h7bI/uFS24eTc+wD/AM4GflTbC6H5UcmZ9WEBFkVlCA0laRXmnfqc51V/aZKWAHYnFZKbkDaSvJVUTFTVTkTSf5OOtbuXNGJ9ru1HyqbqTNJOpNH0y2lIX1VJU4G9WhsxJC1N+h5XfZxrbn20L+kF0QO2dy4c6QWRdKvtTUrnAJC0FWlj5GtIx0wuDDzROhY3NFNs1AkTnqRP2P5Kt3Y9NbbpkfRl0pPcrxmY+jSD1wMWJ+ks0rrEqaRR4F+STtS4smSubmwfkadqtwf2A/5T0u2kAvP8ynalHgS8mtQ7sf1noNqiElgNaD9B5xlgjTJR5svDwF+AvwMvLpxlQfEN0u/YuaReoAeS2niFBouiMoSBM6hvGfZP1eWtwDq2/zXinyxrfdKGjLuA39h+rtWepVZO0zets5Q/CuwMfIl0rOCSJbMNsZHtDUqHmE9nADdLOp9UAO8JVLsOWNKHSC/eViJtfjnE9p1lUy04bP9e0sK2nyO1Fbq+dKbQnygqw4Rn+8L87pO2B51EI+ntHT6lBveQRqiqLiptbyTp1aSp78skPQwsLeklta8BlbQBaSRlX9II1f8rm2geN0pat0lFju0vSPoFaYcypI1wt5XMNILVgcNtTy8dZJTU1LPySUmLAdMlfYXUWqr2PqBhBLGmMoSs03qjmtYgtcu93TZi3vV01U3Vt8ttcPYn9QJ9wPY2hSMNkk8k2Y+U8TnS5owf2L6naLAOJN1FasdyL+lnQKSB1mpbCjWFpMm5h+LynT5ey/GM3UhayvYTHa4XP7O+RdLqwF9J6ymPAJYB/s/274sGC32JojJMeJLeBLyZtMvzh20fmgysa3uLIsGGIendna43paWQ0pEq27da9Ej6tO3jC8dC0j2k9ZNn255ZOs9w8pPyPGrtVtAkki6yvZuke0nT9O0jfK7x2EMASdsApwIvsr2apI2AD9j+cOFo85C0I6m37pOls4TRE0VlmPDyA+9rgWOBz7Z9aDZwhe1HiwSbQGodEW6RtDCwn+2zSmdpJ2k7YG3bUyStRCom7i2dK5Qh6SZS94Kf1n62eu5TuRVpack1+XZtPN42W6ypDBOe7duB2/PmgSfyovFWIbF40XBDSDrH9j6SZtJ5p3pTpz6rWOuVzyL+CLAK8FPgUuCjwJHAdKCaojL3UNyM1PtxCmmN7Zmkc+HDKJB0+dDDDzpdq4ntPw05W/25UlmGY/tAAEkvIxXC/wu8jKhLGi3+80IYMJW00/ef+f4S+VpN6/5uy6dm7Ak8WzrMKKplyuQM0m71G4CDSSfWLAbsUeFmjT2BjYFbAWw/mPs+hj5JmkTa6b9iPgK1VaVNJhU+tfpTngJ33gRzGAPdLaoi6QDShq0NgEdILYauKRoq9C2KyhAGTLLdKiix/c98Dm1NVgBOJPUnnEE69ec64IbaNw+MoIqRSmDNVpseSaeSnuxWq6w/Zcsztt1q0SQpds6Ong8Ah5MKyGkM/HzOIo2o1eqDpMeHVYAHSC+KP1I0UXf/A/yB1KrrCtv3lY0TRkMUlSEMeELSJrZvBZC0KfBU4UyD2D4SII9CbEYaRX0v8G1Jj9let2S+Ppw78h8ZF3NHf3NPzXsrLSgBzpH0TWBZSYeQfw4KZ1og2D4ROFHSobZPDAcV6QAAC9NJREFUKp2nF3m5zrtsv7N0ll7YXlHSeqSDBr6QOy/cbftdhaOFPkRRGcKAw4FzJT2Y77+U1KOwRkuQpuKWybcHgep2K3c7pail1QLJ9hfHLdTwNpI0K78vYIl8v9Wup6Yj5FoNuWeR1lV+lrR8I4wS2ydJWh9YF5jUdv175VJ1ll8E7QH8d+ksvcjrl1cj9QJdg/Q4VtWxrWH+xe7vENpIWpT0BC3SCTBVrVuU9C1gPdLO9JuAG0ltOarcMdnW+mhb0hNzq2XT24Fpto8oEqxPkpYr/T3v0ld1RoM3a1Unb4Z6A+ln9+fAm0g7lPcumasbSV8gFWc/BOb2qWzNvtRE0gzg2ny72vYDhSOFURBFZQhZXj/5cWB124fk6Zh1bF9UONpcki4GVgTuIK2nvAG4w5X/Iku6Ati1VaTn4n2q7R3KJnthSrZAykcHfhhYk7QmrWVp4DrbB5TItSDKXRY2Am7Lp0OtDJxq+y2Fo3WUf8+Gsu0dxz1MmJBi+juEAVNIi/K3zvcfIK31q6aotP3G3Dh8PdJ6yv8PWF/SP0ibdY4uGrC7l5GKntZmohdR9y7akZTcWPR94BfA8cCn2q7PbvhmrRo9Zft5SXPydO3DpGK+Sk16kSbpVaRWXWvQVotEAdxsUVSGMGAt2/tK2h/A9lMa0vCtBnlU8g5JjwGP59tuwBZArUXll0jtkFojKa8HjikXp2/FRoZtt/7P9y+VYQK5RdKypA1Q00jtxm4uG6k7SZ/tdN32seOdpQfnknZ+n0qlvTTD/IuiMoQBz0haglwwSFqLtnO1ayDpMNII5bakncrXkabAv0uFG3Va8okvvwC2zJc+ZfsvJTOFMJK24w1PyUtPJtueUTLTCNrP+55EerFZZZ9KYI7tk0uHCKMr1lSGkEnaBfgP0qL8qaTC7T22ryyZq52kE8i9KW0/VDrP/JC0CmmnZ/tU19XlEr1wkm5rHYMXFjyShl0vW+PGl04kLU46svHfSmcZStIxpOUE59P24j2WcDRbFJUhAHma++XAk6TzaEXaVf1I0WALCElfJrVn+jUDbUNse/dyqbqTtHyHy7PbNhotH09+C64uG15aGrPxJZ8GdLPttUtnGUpSpzPqbbvaNathZFFUhpBJmmZ709I5FkSS7gY2tF3VcoJuJN0HrEo6slHAssBDpJGVQ2xPK5cuhM7ybvXWk/rCpF6mx9r+RrlUYSKJNZUhDLhR0ua2f1U6yALoHmBRKlujOoyLgfNtXwIgaVfgjcA5wP8xsDY0LMDa2oytZvv9NbYZG2K3tvfnAH+1PadUmOFIOrDT9Roby4fexUhlCJmkO4FXAfeTFry3TlGJZtJ9knQeqd/f5QxeP3VYsVDDkHSL7c06XZM03fZrS2UL40fSD0m7vg+0vX7eyHdDrf//eXPhA7b/JekNwIbA92w/VjbZvPJpWy2TgJ2AW2ttLB96EyOVIQx4U+kAC7Cf5ltT/EPSJ4Gz8/19gUfz+cpxlNzE0Yg2Y23OAzaT9ErgO6Tfue8Dby6aqgPbh7bfl7QMcEahOGGURFEZJjxJk23PIh19GMaA7dNLZ5hP7yD1/PwJacT62nxtYWCfgrnC+Kq+zdgQz9ueI+ltwP/ks8tvKx2qR08C1W0oCvMnisoQ0iv53UjTXGbwaSmm4hM0aifpHNv7DNlAMFetSwvyrv9Du3z49+OZJZSRRyRPIa2vXVXSWeQ2YyVzjeDZPKp6INA6SnLRgnm6knQhA48JC5FauZ1bLlEYDbGmMoQwZiR9GriMtIv62aEft33/uIfqQRwhFyB1hAB2pSFtxiStC3yQtO7zB5JeAexr+0uFo81D0uvb7s4B7rf9QKk8YXREURkmvAWl0XGNJH2NdALQq4EZ5MbtpCe9avs8SrqdNEo1jbYj5KKV0MQi6X+B05rYESL3qFy18hOA5srrlfezfVbpLOGFi6IyTHhtjY4nAZsBt5NGJTYEbrK9XalsCwpJi5G+t9sAW+fbY7bXLRqsi+hZGqB5HSEkXQnsThpdnw78DbjK9sdL5monaTLwEWAV0kaiS/P9o4DptvcoGC/0KdZUhgnP9g4Aks4G3m97Zr6/PmkKNPRvCWAysEy+PUjFZ5UDF0r6MHGE3ETXtI4Qy9ieJelgYIrtoyXVNlJ5Bmk5zA3AwaRicjFgD9vTSwYL/YuRyhCyTv0HoydhfyR9C1iPtLP+JuBG0rq0R4sGG0EcIReaKG+I2xU4HfiM7V9JmlHTyKqkmbY3yO8vDDxCai4f3TcWADFSGcKAuySdCpxJ2pV4AHBX2UiNtxqwOPA74M/AA0B1jZiHsv2K0hlCeAGOBS4BrssF5Zqk372azN2wZ/s5SfdGQbngiJHKEDJJk4APAdvnS1cDJ9t+ulyq5sutWdYjrafcBlgf+Adps87RJbMNJWlH27/Mff7mYfvH450phAWJpOdI61MhrVFdgtSjsrVedXKpbKF/UVSG0CNJ59neq3SOppL0clKfv21IfUFXsL1s2VSDSfpcXoc2pcOHbfu94x4qhB7lVlgnAyvnYyU3BHa3fVzhaGGCiKIyhB5Jus32xqVzNImkw0hF5Lakaa/rSAv0rwNm2o4jD0MYJZKuIm18+WbrsUrSHbbXL5ssTBSxpjKE3sUrsPm3BvAj4AjbDxXOMiJJw7ZesX3CeGUJ4QVY0vbNQ44nn1MqTJh4oqgMIYyZmvrj9Wjp/HYdYHNSHz1IR95dXSRRCL17JJ9P3jqrfG+g+hdzYcER098h9CimvycOSVOBvVq7UiUtDZxr+41lk4XQXd7t/S3SkpNHgXuBd9Z6HGpY8MRIZQi9+2TpAGHcrAY803b/GdJUfghVkrQQsJntnSUtBSwUrXrCeIuiMkx4uWFw1yH7VuNg21PHLVQo7QzgZknnk3429gS+VzZSCN3Zfl7SR4FzbD8x4ieEMAZi+jtMeJJWz+9+JL89I799J/Ck7WPHP1UoTdKmQOvc96tt31YyTwgjkfSfwFPADxnoBRnHi4ZxE0VlCJmk62xvO9K1MDHkI+RWpm1Gx/YfyyUKYXj5eNF5ntTjeNEwXmL6O4QBS0nazva1AJK2AZYqnCkUIOlQ4Gjgr8Bz5NM+gGrOUA6hg3WBD5NG2A1cA5xSNFGYUGKkMoQsT3d+F1gmX3oMeK/tW8ulCiVI+j2wpe2/l84SQq8knQPMAs7Kl/YHlrW9T7lUYSKJojKEISRNJv1uPF46SyhD0hXALrajcXRoDEm3295opGshjJWY/g4hk7Q4sBepdcwirVMpYqPOhHQPcKWknwH/al2ME3VC5W6TtJXtGwEkbUk6EjWEcRFFZQgDLgAeB6bRVkiECemP+bZYvoXQBFsCB0pqbShbDbir1Tat1R4thLES098hZJLusL1+6RwhhPBCtLVH6yhO1gljLUYqQxhwvaQNbM8sHSSUJWkl4BPAesCk1nXbOxYLFcIIomgMpS1UOkAIFdkOmCbpbkkzJM2UNKN0qFDEWcBvgFcAnwPuA35VMlAIIdQupr9DyLpNHcWr/4lH0jTbm0qa0VqHJukq268vnS2EEGoV099hwpM02fYsYHbpLKEaz+a3D0n6d+BB4OUF84QQQvVipDJMeJIusr1b2xFnavuw44iziUfSbqTTSFYFTgImA8fYvrBosBBCqFgUlSFkks4Argausf2b0nlCXSQdbvt/SucIIYRaRVEZQiZpR9JmndcBawK3kQrME4sGC1WQ9Efbq5XOEUIItYqiMoQ2khYGNgd2AD4IPGX71WVThRpI+pPtVUvnCCGEWsVGnRAySZcDSwE3kNbTbW774bKpQkXiFXgIIQwjisoQBswANgXWJx3X+JikG2w/VTZWGC+SZtO5eBSwxDjHCSGERonp7xCGkPQi4CDgSOAlthcvHCmEEEKoXoxUhpBJ+ihpk86mwP3Ad0nT4CGEEEIYQRSVIQxYAjgBmGZ7TukwIYQQQpPE9HcIIYQQQujbQqUDhBBCCCGE5ouiMoQQQggh9C2KyhBCCCGE0LcoKkMIIYQQQt+iqAwhhBBCCH37/wE5OfK5P+OzrgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x576 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, ax = plt.subplots(figsize=(10, 8))\n",
    "corr = df.corr()\n",
    "sns.heatmap(corr, mask=np.zeros_like(corr, dtype=np.bool), cmap=sns.diverging_palette(220, 10, as_cmap=True),square=True, ax=ax)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Notes:\n",
    "`numerical_cols_of_interest = ['hourofday','aircraft_registration','wind_direction_degrees','Wind_speed_kt', 'Wind_And_Gust_kt','Landing_RVR_or_Vis_m','temperature_c','dew_point_c','relative_humidity','pressure_altimeter_in',\n",
    "'RunwayTrueTrack']`\n",
    "\n",
    "`categorical_cols_of_interest = ['Landing_Visibilty_Class_RVR','present_weather_codes']`\n",
    "\n",
    "\n",
    "* **Landing_RWY** is our main column to classify and predict.\n",
    "* Low Quality Data, i.e where Landing_RWY_Quality > 5 (is bad) has been removed in the raw dataset itself.\n",
    "* Landing time is used to sort the dataframe.\n",
    "* Many columns such as _Landing_Coverage_Alt_ft_, _Landing_RVR_or_Vis_Quality_ and manually calculated columns like *Landing_RWY_Quality, Dist_To_Threshold_m, Mean_DeltaTrack, Nb_Snapshot_InRWY, Mean_ProjDistance, Nb_Pt_in_Last_min*  have been dropped. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}