Week8/PandasExamples.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Dealing with other data types\n",
    "Python can deal with basically any type of data you throw at it. The community have provided many packages that make things easy, today we will look at the \"pyshp\" (for dealing with shapefiles) and \"pandas\" (great for tables nad time series) packages.\n",
    "\n",
    "Data for this exercised was downloaded from http://www.bom.gov.au/water/groundwater/explorer/map.shtml"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Shapefiles\n",
    "Shapefiles are a very common file format for GIS data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Load the required modules\n",
    "import shapefile\n",
    "\n",
    "#NOTE: Weirdly and confusingly, this package is called \"pyshp\" but you call it via the name \"shapefile\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on module shapefile:\n",
      "\n",
      "NAME\n",
      "    shapefile\n",
      "\n",
      "FILE\n",
      "    c:\\users\\nbutter\\anaconda3\\envs\\pygeol\\lib\\site-packages\\shapefile.py\n",
      "\n",
      "DESCRIPTION\n",
      "    shapefile.py\n",
      "    Provides read and write support for ESRI Shapefiles.\n",
      "    author: jlawhead<at>geospatialpython.com\n",
      "    date: 2017/08/24\n",
      "    version: 1.2.12\n",
      "    Compatible with Python versions 2.7-3.x\n",
      "\n",
      "CLASSES\n",
      "    exceptions.Exception(exceptions.BaseException)\n",
      "        ShapefileException\n",
      "    Reader\n",
      "    Writer\n",
      "        Editor\n",
      "    \n",
      "    class Editor(Writer)\n",
      "     |  Methods defined here:\n",
      "     |  \n",
      "     |  __init__(self, shapefile=None, shapeType=1, autoBalance=1)\n",
      "     |  \n",
      "     |  balance(self)\n",
      "     |      Adds a corresponding empty attribute or null geometry record depending\n",
      "     |      on which type of record was created to make sure all three files\n",
      "     |      are in synch.\n",
      "     |  \n",
      "     |  delete(self, shape=None, part=None, point=None)\n",
      "     |      Deletes the specified part of any shape by specifying a shape\n",
      "     |      number, part number, or point number.\n",
      "     |  \n",
      "     |  point(self, x=None, y=None, z=None, m=None, shape=None, part=None, point=None, addr=None)\n",
      "     |      Creates/updates a point shape. The arguments allows\n",
      "     |      you to update a specific point by shape, part, point of any\n",
      "     |      shape type.\n",
      "     |  \n",
      "     |  select(self, expr)\n",
      "     |      Select one or more shapes (to be implemented)\n",
      "     |  \n",
      "     |  validate(self)\n",
      "     |      An optional method to try and validate the shapefile\n",
      "     |      as much as possible before writing it (not implemented).\n",
      "     |  \n",
      "     |  ----------------------------------------------------------------------\n",
      "     |  Methods inherited from Writer:\n",
      "     |  \n",
      "     |  bbox(self)\n",
      "     |      Returns the current bounding box for the shapefile which is\n",
      "     |      the lower-left and upper-right corners. It does not contain the\n",
      "     |      elevation or measure extremes.\n",
      "     |  \n",
      "     |  field(self, name, fieldType='C', size='50', decimal=0)\n",
      "     |      Adds a dbf field descriptor to the shapefile.\n",
      "     |  \n",
      "     |  line(self, parts=[], shapeType=3)\n",
      "     |      Creates a line shape. This method is just a convienience method\n",
      "     |      which wraps 'poly()'.\n",
      "     |  \n",
      "     |  mbox(self)\n",
      "     |      Returns the current m extremes for the shapefile.\n",
      "     |  \n",
      "     |  null(self)\n",
      "     |      Creates a null shape.\n",
      "     |  \n",
      "     |  poly(self, parts=[], shapeType=5, partTypes=[])\n",
      "     |      Creates a shape that has multiple collections of points (parts)\n",
      "     |      including lines, polygons, and even multipoint shapes. If no shape type\n",
      "     |      is specified it defaults to 'polygon'. If no part types are specified\n",
      "     |      (which they normally won't be) then all parts default to the shape type.\n",
      "     |  \n",
      "     |  record(self, *recordList, **recordDict)\n",
      "     |      Creates a dbf attribute record. You can submit either a sequence of\n",
      "     |      field values or keyword arguments of field names and values. Before\n",
      "     |      adding records you must add fields for the record values using the\n",
      "     |      fields() method. If the record values exceed the number of fields the\n",
      "     |      extra ones won't be added. In the case of using keyword arguments to specify\n",
      "     |      field/value pairs only fields matching the already registered fields\n",
      "     |      will be added.\n",
      "     |  \n",
      "     |  save(self, target=None, shp=None, shx=None, dbf=None)\n",
      "     |      Save the shapefile data to three files or\n",
      "     |      three file-like objects. SHP and DBF files can also\n",
      "     |      be written exclusively using saveShp, saveShx, and saveDbf respectively.\n",
      "     |      If target is specified but not shp,shx, or dbf then the target path and\n",
      "     |      file name are used.  If no options or specified, a unique base file name\n",
      "     |      is generated to save the files and the base file name is returned as a \n",
      "     |      string.\n",
      "     |  \n",
      "     |  saveDbf(self, target)\n",
      "     |      Save a dbf file.\n",
      "     |  \n",
      "     |  saveShp(self, target)\n",
      "     |      Save an shp file.\n",
      "     |  \n",
      "     |  saveShx(self, target)\n",
      "     |      Save an shx file.\n",
      "     |  \n",
      "     |  shape(self, i)\n",
      "     |  \n",
      "     |  shapes(self)\n",
      "     |      Return the current list of shapes.\n",
      "     |  \n",
      "     |  zbox(self)\n",
      "     |      Returns the current z extremes for the shapefile.\n",
      "    \n",
      "    class Reader\n",
      "     |  Reads the three files of a shapefile as a unit or\n",
      "     |  separately.  If one of the three files (.shp, .shx,\n",
      "     |  .dbf) is missing no exception is thrown until you try\n",
      "     |  to call a method that depends on that particular file.\n",
      "     |  The .shx index file is used if available for efficiency\n",
      "     |  but is not required to read the geometry from the .shp\n",
      "     |  file. The \"shapefile\" argument in the constructor is the\n",
      "     |  name of the file you want to open.\n",
      "     |  \n",
      "     |  You can instantiate a Reader without specifying a shapefile\n",
      "     |  and then specify one later with the load() method.\n",
      "     |  \n",
      "     |  Only the shapefile headers are read upon loading. Content\n",
      "     |  within each file is only accessed when required and as\n",
      "     |  efficiently as possible. Shapefiles are usually not large\n",
      "     |  but they can be.\n",
      "     |  \n",
      "     |  Methods defined here:\n",
      "     |  \n",
      "     |  __init__(self, *args, **kwargs)\n",
      "     |  \n",
      "     |  iterRecords(self)\n",
      "     |      Serves up records in a dbf file as an iterator.\n",
      "     |      Useful for large shapefiles or dbf files.\n",
      "     |  \n",
      "     |  iterShapeRecords(self)\n",
      "     |      Returns a generator of combination geometry/attribute records for\n",
      "     |      all records in a shapefile.\n",
      "     |  \n",
      "     |  iterShapes(self)\n",
      "     |      Serves up shapes in a shapefile as an iterator. Useful\n",
      "     |      for handling large shapefiles.\n",
      "     |  \n",
      "     |  load(self, shapefile=None)\n",
      "     |      Opens a shapefile from a filename or file-like\n",
      "     |      object. Normally this method would be called by the\n",
      "     |      constructor with the file name as an argument.\n",
      "     |  \n",
      "     |  record(self, i=0)\n",
      "     |      Returns a specific dbf record based on the supplied index.\n",
      "     |  \n",
      "     |  records(self)\n",
      "     |      Returns all records in a dbf file.\n",
      "     |  \n",
      "     |  shape(self, i=0)\n",
      "     |      Returns a shape object for a shape in the the geometry\n",
      "     |      record file.\n",
      "     |  \n",
      "     |  shapeRecord(self, i=0)\n",
      "     |      Returns a combination geometry and attribute record for the\n",
      "     |      supplied record index.\n",
      "     |  \n",
      "     |  shapeRecords(self)\n",
      "     |      Returns a list of combination geometry/attribute records for\n",
      "     |      all records in a shapefile.\n",
      "     |  \n",
      "     |  shapes(self)\n",
      "     |      Returns all shapes in a shapefile.\n",
      "    \n",
      "    class ShapefileException(exceptions.Exception)\n",
      "     |  An exception to handle shapefile specific problems.\n",
      "     |  \n",
      "     |  Method resolution order:\n",
      "     |      ShapefileException\n",
      "     |      exceptions.Exception\n",
      "     |      exceptions.BaseException\n",
      "     |      __builtin__.object\n",
      "     |  \n",
      "     |  Data descriptors defined here:\n",
      "     |  \n",
      "     |  __weakref__\n",
      "     |      list of weak references to the object (if defined)\n",
      "     |  \n",
      "     |  ----------------------------------------------------------------------\n",
      "     |  Methods inherited from exceptions.Exception:\n",
      "     |  \n",
      "     |  __init__(...)\n",
      "     |      x.__init__(...) initializes x; see help(type(x)) for signature\n",
      "     |  \n",
      "     |  ----------------------------------------------------------------------\n",
      "     |  Data and other attributes inherited from exceptions.Exception:\n",
      "     |  \n",
      "     |  __new__ = <built-in method __new__ of type object>\n",
      "     |      T.__new__(S, ...) -> a new object with type S, a subtype of T\n",
      "     |  \n",
      "     |  ----------------------------------------------------------------------\n",
      "     |  Methods inherited from exceptions.BaseException:\n",
      "     |  \n",
      "     |  __delattr__(...)\n",
      "     |      x.__delattr__('name') <==> del x.name\n",
      "     |  \n",
      "     |  __getattribute__(...)\n",
      "     |      x.__getattribute__('name') <==> x.name\n",
      "     |  \n",
      "     |  __getitem__(...)\n",
      "     |      x.__getitem__(y) <==> x[y]\n",
      "     |  \n",
      "     |  __getslice__(...)\n",
      "     |      x.__getslice__(i, j) <==> x[i:j]\n",
      "     |      \n",
      "     |      Use of negative indices is not supported.\n",
      "     |  \n",
      "     |  __reduce__(...)\n",
      "     |  \n",
      "     |  __repr__(...)\n",
      "     |      x.__repr__() <==> repr(x)\n",
      "     |  \n",
      "     |  __setattr__(...)\n",
      "     |      x.__setattr__('name', value) <==> x.name = value\n",
      "     |  \n",
      "     |  __setstate__(...)\n",
      "     |  \n",
      "     |  __str__(...)\n",
      "     |      x.__str__() <==> str(x)\n",
      "     |  \n",
      "     |  __unicode__(...)\n",
      "     |  \n",
      "     |  ----------------------------------------------------------------------\n",
      "     |  Data descriptors inherited from exceptions.BaseException:\n",
      "     |  \n",
      "     |  __dict__\n",
      "     |  \n",
      "     |  args\n",
      "     |  \n",
      "     |  message\n",
      "    \n",
      "    class Writer\n",
      "     |  Provides write support for ESRI Shapefiles.\n",
      "     |  \n",
      "     |  Methods defined here:\n",
      "     |  \n",
      "     |  __init__(self, shapeType=None)\n",
      "     |  \n",
      "     |  bbox(self)\n",
      "     |      Returns the current bounding box for the shapefile which is\n",
      "     |      the lower-left and upper-right corners. It does not contain the\n",
      "     |      elevation or measure extremes.\n",
      "     |  \n",
      "     |  field(self, name, fieldType='C', size='50', decimal=0)\n",
      "     |      Adds a dbf field descriptor to the shapefile.\n",
      "     |  \n",
      "     |  line(self, parts=[], shapeType=3)\n",
      "     |      Creates a line shape. This method is just a convienience method\n",
      "     |      which wraps 'poly()'.\n",
      "     |  \n",
      "     |  mbox(self)\n",
      "     |      Returns the current m extremes for the shapefile.\n",
      "     |  \n",
      "     |  null(self)\n",
      "     |      Creates a null shape.\n",
      "     |  \n",
      "     |  point(self, x, y, z=0, m=0, shapeType=1)\n",
      "     |      Creates a point shape.\n",
      "     |  \n",
      "     |  poly(self, parts=[], shapeType=5, partTypes=[])\n",
      "     |      Creates a shape that has multiple collections of points (parts)\n",
      "     |      including lines, polygons, and even multipoint shapes. If no shape type\n",
      "     |      is specified it defaults to 'polygon'. If no part types are specified\n",
      "     |      (which they normally won't be) then all parts default to the shape type.\n",
      "     |  \n",
      "     |  record(self, *recordList, **recordDict)\n",
      "     |      Creates a dbf attribute record. You can submit either a sequence of\n",
      "     |      field values or keyword arguments of field names and values. Before\n",
      "     |      adding records you must add fields for the record values using the\n",
      "     |      fields() method. If the record values exceed the number of fields the\n",
      "     |      extra ones won't be added. In the case of using keyword arguments to specify\n",
      "     |      field/value pairs only fields matching the already registered fields\n",
      "     |      will be added.\n",
      "     |  \n",
      "     |  save(self, target=None, shp=None, shx=None, dbf=None)\n",
      "     |      Save the shapefile data to three files or\n",
      "     |      three file-like objects. SHP and DBF files can also\n",
      "     |      be written exclusively using saveShp, saveShx, and saveDbf respectively.\n",
      "     |      If target is specified but not shp,shx, or dbf then the target path and\n",
      "     |      file name are used.  If no options or specified, a unique base file name\n",
      "     |      is generated to save the files and the base file name is returned as a \n",
      "     |      string.\n",
      "     |  \n",
      "     |  saveDbf(self, target)\n",
      "     |      Save a dbf file.\n",
      "     |  \n",
      "     |  saveShp(self, target)\n",
      "     |      Save an shp file.\n",
      "     |  \n",
      "     |  saveShx(self, target)\n",
      "     |      Save an shx file.\n",
      "     |  \n",
      "     |  shape(self, i)\n",
      "     |  \n",
      "     |  shapes(self)\n",
      "     |      Return the current list of shapes.\n",
      "     |  \n",
      "     |  zbox(self)\n",
      "     |      Returns the current z extremes for the shapefile.\n",
      "\n",
      "FUNCTIONS\n",
      "    b(v)\n",
      "    \n",
      "    calcsize(...)\n",
      "        Return size of C struct described by format string fmt.\n",
      "    \n",
      "    is_string(v)\n",
      "    \n",
      "    pack(...)\n",
      "        Return string containing values v1, v2, ... packed according to fmt.\n",
      "    \n",
      "    signed_area(coords)\n",
      "        Return the signed area enclosed by a ring using the linear time\n",
      "        algorithm. A value >= 0 indicates a counter-clockwise oriented ring.\n",
      "    \n",
      "    test(**kwargs)\n",
      "        # Begin Testing\n",
      "    \n",
      "    u(v)\n",
      "    \n",
      "    unpack(...)\n",
      "        Unpack the string containing packed C structure data, according to fmt.\n",
      "        Requires len(string) == calcsize(fmt).\n",
      "\n",
      "DATA\n",
      "    MISSING = [None, '']\n",
      "    MULTIPATCH = 31\n",
      "    MULTIPOINT = 8\n",
      "    MULTIPOINTM = 28\n",
      "    MULTIPOINTZ = 18\n",
      "    NULL = 0\n",
      "    POINT = 1\n",
      "    POINTM = 21\n",
      "    POINTZ = 11\n",
      "    POLYGON = 5\n",
      "    POLYGONM = 25\n",
      "    POLYGONZ = 15\n",
      "    POLYLINE = 3\n",
      "    POLYLINEM = 23\n",
      "    POLYLINEZ = 13\n",
      "    PYTHON3 = False\n",
      "    __version__ = '1.2.12'\n",
      "\n",
      "VERSION\n",
      "    1.2.12\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(shapefile)\n",
    "#Or check out the help pages https://github.com/GeospatialPython/pyshp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Set the filename\n",
    "boreshape='../data/shp_torrens_river/NGIS_BoreLine.shp'\n",
    "\n",
    "#read in the file\n",
    "shapeRead = shapefile.Reader(boreshape)\n",
    "\n",
    "#And save out some of the shape file attributes\n",
    "recs    = shapeRead.records()\n",
    "shapes  = shapeRead.shapes()\n",
    "fields  = shapeRead.fields\n",
    "Nshp    = len(shapes)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "7635\n"
     ]
    }
   ],
   "source": [
    "print(Nshp) #print the Number of items in the shapefile"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('DeletionFlag', 'C', 1, 0),\n",
       " ['HydroID', 'N', 10, 0],\n",
       " ['HydroCode', 'C', 30, 0],\n",
       " ['BoreID', 'N', 10, 0],\n",
       " ['TopElev', 'F', 19, 11],\n",
       " ['BottomElev', 'F', 19, 11],\n",
       " ['HGUID', 'N', 10, 0],\n",
       " ['HGUNumber', 'N', 10, 0],\n",
       " ['NafHGUNumb', 'N', 10, 0],\n",
       " ['SHAPE_Leng', 'F', 19, 11]]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fields[:]#print the fields"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[32001999, '652800645', 30027773, 6.74, -74.26, 31000043, 1042, 104005, 0.0]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "recs[0] #print the first record, then this is a list that can be subscripted further"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[(591975.5150000006, -3816141.8817), (591975.5150000006, -3816141.8817)]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "shapes[0].points #print the point values of the first shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Shapefiles are not a native python format, but the community have developed tools for exploring them. The package we have used \"pyshp\" imported with the name \"shapefile\" (for some non-consistent weird reason), is one example of working with shapefiles. Alternatives exist."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## More table manipulation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Import the module\n",
    "import pandas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "#read in the data\n",
    "log_data=pandas.read_csv(\"../data/shp_torrens_river/NGIS_LithologyLog.csv\",\\\n",
    "                         header=0,sep=',',skipinitialspace=True,quotechar ='\"',\\\n",
    "                         usecols=list(range(0,13)),\\\n",
    "                         skiprows=[453,456,458,460,689,697,720,723,726,839,880,884,885,890,898,934])\n",
    "\n",
    "#This data was weird because it has quotation marks to signify inches inside comments within the file, \n",
    "#making automatic reading of it tricky"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "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>OBJECTID</th>\n",
       "      <th>BoreID</th>\n",
       "      <th>HydroCode</th>\n",
       "      <th>RefElev</th>\n",
       "      <th>RefElevDesc</th>\n",
       "      <th>FromDepth</th>\n",
       "      <th>ToDepth</th>\n",
       "      <th>TopElev</th>\n",
       "      <th>BottomElev</th>\n",
       "      <th>MajorLithCode</th>\n",
       "      <th>MinorLithCode</th>\n",
       "      <th>Description</th>\n",
       "      <th>Source</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1769789</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>18.0</td>\n",
       "      <td>19.5</td>\n",
       "      <td>39.25</td>\n",
       "      <td>37.75</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1769790</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>19.5</td>\n",
       "      <td>22.0</td>\n",
       "      <td>37.75</td>\n",
       "      <td>35.25</td>\n",
       "      <td>ROCK</td>\n",
       "      <td>None</td>\n",
       "      <td>Rocks and sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1769791</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>22.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>35.25</td>\n",
       "      <td>33.25</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1770725</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>0.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>No sample</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1770726</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>6.0</td>\n",
       "      <td>15.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>-11.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - orange-red grey, mottled; stiff-sticky....</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>1770727</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>15.0</td>\n",
       "      <td>18.0</td>\n",
       "      <td>-11.0</td>\n",
       "      <td>-14.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay - red brown, grey, green, white mot...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>1770728</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>18.0</td>\n",
       "      <td>21.0</td>\n",
       "      <td>-14.0</td>\n",
       "      <td>-17.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - red green mottled; stiff. fine sand, qtz.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>1770729</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>21.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>-17.0</td>\n",
       "      <td>-20.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand - orange green mottled; stiff. fin...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>1770730</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>24.0</td>\n",
       "      <td>28.5</td>\n",
       "      <td>-20.0</td>\n",
       "      <td>-24.5</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - red-green grey, mottled; stiff.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>1770731</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>28.5</td>\n",
       "      <td>30.0</td>\n",
       "      <td>-24.5</td>\n",
       "      <td>-26.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand - red green, mottled; sticky. fine...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>1770732</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>30.0</td>\n",
       "      <td>34.5</td>\n",
       "      <td>-26.0</td>\n",
       "      <td>-30.5</td>\n",
       "      <td>SAND</td>\n",
       "      <td>GRVL</td>\n",
       "      <td>Sand/gravel - fine sand - fine gravel, clear, ...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>1770733</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>34.5</td>\n",
       "      <td>36.0</td>\n",
       "      <td>-30.5</td>\n",
       "      <td>-32.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - grey red mottled; soft - plastic.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>1770734</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>36.0</td>\n",
       "      <td>44.0</td>\n",
       "      <td>-32.0</td>\n",
       "      <td>-40.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Sand - fine-medium, clear, opaque, fe stained ...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>1770735</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>44.0</td>\n",
       "      <td>48.0</td>\n",
       "      <td>-40.0</td>\n",
       "      <td>-44.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - grey green; stiff - plastic; minor fine...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>1770736</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>48.0</td>\n",
       "      <td>51.0</td>\n",
       "      <td>-44.0</td>\n",
       "      <td>-47.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>No sample.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>1770737</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>51.0</td>\n",
       "      <td>54.0</td>\n",
       "      <td>-47.0</td>\n",
       "      <td>-50.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - red green; plastic, sticky; minor fine-...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>1770738</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>57.0</td>\n",
       "      <td>64.0</td>\n",
       "      <td>-53.0</td>\n",
       "      <td>-60.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - green grey; stiff. minor silt/fine sand...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>1770739</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>64.0</td>\n",
       "      <td>66.0</td>\n",
       "      <td>-60.0</td>\n",
       "      <td>-62.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay - cream; fine sand.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>1770740</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>66.0</td>\n",
       "      <td>75.0</td>\n",
       "      <td>-62.0</td>\n",
       "      <td>-71.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand - cream-pale green grey, fine sand...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>1770741</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>75.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>-71.0</td>\n",
       "      <td>-86.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>GRVL</td>\n",
       "      <td>Sand/gravel - coarse (1-3mm), predominantly op...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>1770742</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>90.0</td>\n",
       "      <td>98.5</td>\n",
       "      <td>-86.0</td>\n",
       "      <td>-94.5</td>\n",
       "      <td>SILT</td>\n",
       "      <td>SAND</td>\n",
       "      <td>Silt/sand - grey; fine sand (0.05-0.1mm), qtz,...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>1770743</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>98.5</td>\n",
       "      <td>114.0</td>\n",
       "      <td>-94.5</td>\n",
       "      <td>-110.0</td>\n",
       "      <td>LMST</td>\n",
       "      <td>None</td>\n",
       "      <td>Limestone - shell and coral fragments (drill c...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>1770744</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>114.0</td>\n",
       "      <td>117.0</td>\n",
       "      <td>-110.0</td>\n",
       "      <td>-113.0</td>\n",
       "      <td>SDST</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandstone - yellow grey; interbedded with sand...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>1770745</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>117.0</td>\n",
       "      <td>125.0</td>\n",
       "      <td>-113.0</td>\n",
       "      <td>-121.0</td>\n",
       "      <td>SDST</td>\n",
       "      <td>LMST</td>\n",
       "      <td>Sandstone/limestone - yellow grey; carbonate +...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>1770746</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>125.0</td>\n",
       "      <td>147.3</td>\n",
       "      <td>-121.0</td>\n",
       "      <td>-143.3</td>\n",
       "      <td>SDST</td>\n",
       "      <td>LMST</td>\n",
       "      <td>Sandstone/limestone - grey white and grey; fin...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>1770747</td>\n",
       "      <td>30141910</td>\n",
       "      <td>662816624</td>\n",
       "      <td>4.0</td>\n",
       "      <td>NGS</td>\n",
       "      <td>147.3</td>\n",
       "      <td>148.3</td>\n",
       "      <td>-143.3</td>\n",
       "      <td>-144.3</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay - steel grey-blue; sticky. (munno para cl...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>1771871</td>\n",
       "      <td>30141911</td>\n",
       "      <td>662816625</td>\n",
       "      <td>4.0</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>4.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay - yellow brown; minor sand component.</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>1771873</td>\n",
       "      <td>30141911</td>\n",
       "      <td>662816625</td>\n",
       "      <td>4.0</td>\n",
       "      <td>UNK</td>\n",
       "      <td>6.0</td>\n",
       "      <td>21.0</td>\n",
       "      <td>-2.0</td>\n",
       "      <td>-17.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay (gravel) - grey and brown, mottled;...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>1771874</td>\n",
       "      <td>30141911</td>\n",
       "      <td>662816625</td>\n",
       "      <td>4.0</td>\n",
       "      <td>UNK</td>\n",
       "      <td>21.0</td>\n",
       "      <td>36.0</td>\n",
       "      <td>-17.0</td>\n",
       "      <td>-32.0</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay (gravel) - grey and brown, mottled;...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>1771875</td>\n",
       "      <td>30141911</td>\n",
       "      <td>662816625</td>\n",
       "      <td>4.0</td>\n",
       "      <td>UNK</td>\n",
       "      <td>36.0</td>\n",
       "      <td>42.0</td>\n",
       "      <td>-32.0</td>\n",
       "      <td>-38.0</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand - orange brown; fine (0.1mm) - coa...</td>\n",
       "      <td>SAGeodata</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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70127</th>\n",
       "      <td>4120275</td>\n",
       "      <td>30300938</td>\n",
       "      <td>662829132</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.8</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SILT</td>\n",
       "      <td>None</td>\n",
       "      <td>Gravelly silt</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70128</th>\n",
       "      <td>4120276</td>\n",
       "      <td>30300938</td>\n",
       "      <td>662829132</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.8</td>\n",
       "      <td>7.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70129</th>\n",
       "      <td>4120277</td>\n",
       "      <td>30300938</td>\n",
       "      <td>662829132</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>7.5</td>\n",
       "      <td>8.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70130</th>\n",
       "      <td>4120278</td>\n",
       "      <td>30300938</td>\n",
       "      <td>662829132</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>8.5</td>\n",
       "      <td>10.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Gravelly/silty clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70131</th>\n",
       "      <td>4120279</td>\n",
       "      <td>30300938</td>\n",
       "      <td>662829132</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>10.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70132</th>\n",
       "      <td>4120280</td>\n",
       "      <td>30301335</td>\n",
       "      <td>662829148</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>16.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70133</th>\n",
       "      <td>4120281</td>\n",
       "      <td>30301339</td>\n",
       "      <td>662829152</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11.2</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70134</th>\n",
       "      <td>4120282</td>\n",
       "      <td>30301339</td>\n",
       "      <td>662829152</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>11.2</td>\n",
       "      <td>12.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>GRVL</td>\n",
       "      <td>None</td>\n",
       "      <td>Gravels</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70135</th>\n",
       "      <td>4120283</td>\n",
       "      <td>30301339</td>\n",
       "      <td>662829152</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>12.0</td>\n",
       "      <td>15.2</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70136</th>\n",
       "      <td>4120284</td>\n",
       "      <td>30301339</td>\n",
       "      <td>662829152</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>15.2</td>\n",
       "      <td>17.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay, clayey sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70137</th>\n",
       "      <td>4120302</td>\n",
       "      <td>30302253</td>\n",
       "      <td>662829174</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.3</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand dark brown</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70138</th>\n",
       "      <td>4120303</td>\n",
       "      <td>30302253</td>\n",
       "      <td>662829174</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.3</td>\n",
       "      <td>1.1</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Sand clay orange brown</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70139</th>\n",
       "      <td>4120304</td>\n",
       "      <td>30302253</td>\n",
       "      <td>662829174</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>1.1</td>\n",
       "      <td>4.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand orange grey</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70140</th>\n",
       "      <td>4120305</td>\n",
       "      <td>30302252</td>\n",
       "      <td>662829173</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.3</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand dark brown</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70141</th>\n",
       "      <td>4120306</td>\n",
       "      <td>30302252</td>\n",
       "      <td>662829173</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.3</td>\n",
       "      <td>1.1</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Sand clay orange brown</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70142</th>\n",
       "      <td>4120307</td>\n",
       "      <td>30302252</td>\n",
       "      <td>662829173</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>1.1</td>\n",
       "      <td>4.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand orange grey</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70143</th>\n",
       "      <td>4120312</td>\n",
       "      <td>30303742</td>\n",
       "      <td>662829205</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>FILL</td>\n",
       "      <td>None</td>\n",
       "      <td>Pug hole no returns</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70144</th>\n",
       "      <td>4120313</td>\n",
       "      <td>30303742</td>\n",
       "      <td>662829205</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>8.0</td>\n",
       "      <td>17.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>UKN</td>\n",
       "      <td>None</td>\n",
       "      <td>Cuttings lost into pug hole</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70145</th>\n",
       "      <td>4120314</td>\n",
       "      <td>30303990</td>\n",
       "      <td>652802870</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>FILL</td>\n",
       "      <td>None</td>\n",
       "      <td>Fill material, gravel, sand, bricks, rubble</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70146</th>\n",
       "      <td>4120315</td>\n",
       "      <td>30303990</td>\n",
       "      <td>652802870</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>2.0</td>\n",
       "      <td>2.2</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay grey brown, medium plasticity</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70147</th>\n",
       "      <td>4120316</td>\n",
       "      <td>30303990</td>\n",
       "      <td>652802870</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>2.2</td>\n",
       "      <td>5.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Sand, grey to black, medium to coarse grained</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70148</th>\n",
       "      <td>4120317</td>\n",
       "      <td>30303994</td>\n",
       "      <td>652802874</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>FILL</td>\n",
       "      <td>None</td>\n",
       "      <td>Fill - cobble sized lumps of slag with gravel ...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70149</th>\n",
       "      <td>4120318</td>\n",
       "      <td>30303994</td>\n",
       "      <td>652802874</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>2.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Silty sand - grey to black medium to coarse gr...</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70150</th>\n",
       "      <td>4120343</td>\n",
       "      <td>30304039</td>\n",
       "      <td>662829228</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>FILL</td>\n",
       "      <td>None</td>\n",
       "      <td>Fill</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70151</th>\n",
       "      <td>4120344</td>\n",
       "      <td>30304039</td>\n",
       "      <td>662829228</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>1.0</td>\n",
       "      <td>9.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70152</th>\n",
       "      <td>4120345</td>\n",
       "      <td>30304039</td>\n",
       "      <td>662829228</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>9.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Sandy clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70153</th>\n",
       "      <td>4120346</td>\n",
       "      <td>30304039</td>\n",
       "      <td>662829228</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>10.0</td>\n",
       "      <td>12.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70154</th>\n",
       "      <td>4120347</td>\n",
       "      <td>30304050</td>\n",
       "      <td>652802882</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.3</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>FILL</td>\n",
       "      <td>None</td>\n",
       "      <td>Fill</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70155</th>\n",
       "      <td>4120348</td>\n",
       "      <td>30304050</td>\n",
       "      <td>652802882</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.3</td>\n",
       "      <td>0.8</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Clayey sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70156</th>\n",
       "      <td>4120349</td>\n",
       "      <td>30304050</td>\n",
       "      <td>652802882</td>\n",
       "      <td>None</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0.8</td>\n",
       "      <td>3.5</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>SAND</td>\n",
       "      <td>None</td>\n",
       "      <td>Sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>70157 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       OBJECTID    BoreID  HydroCode RefElev RefElevDesc  FromDepth  ToDepth  \\\n",
       "0       1769789  30062892  662815923   57.25         NGS       18.0     19.5   \n",
       "1       1769790  30062892  662815923   57.25         NGS       19.5     22.0   \n",
       "2       1769791  30062892  662815923   57.25         NGS       22.0     24.0   \n",
       "3       1770725  30141910  662816624     4.0         NGS        0.0      6.0   \n",
       "4       1770726  30141910  662816624     4.0         NGS        6.0     15.0   \n",
       "5       1770727  30141910  662816624     4.0         NGS       15.0     18.0   \n",
       "6       1770728  30141910  662816624     4.0         NGS       18.0     21.0   \n",
       "7       1770729  30141910  662816624     4.0         NGS       21.0     24.0   \n",
       "8       1770730  30141910  662816624     4.0         NGS       24.0     28.5   \n",
       "9       1770731  30141910  662816624     4.0         NGS       28.5     30.0   \n",
       "10      1770732  30141910  662816624     4.0         NGS       30.0     34.5   \n",
       "11      1770733  30141910  662816624     4.0         NGS       34.5     36.0   \n",
       "12      1770734  30141910  662816624     4.0         NGS       36.0     44.0   \n",
       "13      1770735  30141910  662816624     4.0         NGS       44.0     48.0   \n",
       "14      1770736  30141910  662816624     4.0         NGS       48.0     51.0   \n",
       "15      1770737  30141910  662816624     4.0         NGS       51.0     54.0   \n",
       "16      1770738  30141910  662816624     4.0         NGS       57.0     64.0   \n",
       "17      1770739  30141910  662816624     4.0         NGS       64.0     66.0   \n",
       "18      1770740  30141910  662816624     4.0         NGS       66.0     75.0   \n",
       "19      1770741  30141910  662816624     4.0         NGS       75.0     90.0   \n",
       "20      1770742  30141910  662816624     4.0         NGS       90.0     98.5   \n",
       "21      1770743  30141910  662816624     4.0         NGS       98.5    114.0   \n",
       "22      1770744  30141910  662816624     4.0         NGS      114.0    117.0   \n",
       "23      1770745  30141910  662816624     4.0         NGS      117.0    125.0   \n",
       "24      1770746  30141910  662816624     4.0         NGS      125.0    147.3   \n",
       "25      1770747  30141910  662816624     4.0         NGS      147.3    148.3   \n",
       "26      1771871  30141911  662816625     4.0         UNK        0.0      6.0   \n",
       "27      1771873  30141911  662816625     4.0         UNK        6.0     21.0   \n",
       "28      1771874  30141911  662816625     4.0         UNK       21.0     36.0   \n",
       "29      1771875  30141911  662816625     4.0         UNK       36.0     42.0   \n",
       "...         ...       ...        ...     ...         ...        ...      ...   \n",
       "70127   4120275  30300938  662829132    None         UNK        0.0      0.8   \n",
       "70128   4120276  30300938  662829132    None         UNK        0.8      7.5   \n",
       "70129   4120277  30300938  662829132    None         UNK        7.5      8.5   \n",
       "70130   4120278  30300938  662829132    None         UNK        8.5     10.0   \n",
       "70131   4120279  30300938  662829132    None         UNK       10.0     17.0   \n",
       "70132   4120280  30301335  662829148    None         UNK        0.0     16.0   \n",
       "70133   4120281  30301339  662829152    None         UNK        0.0     11.2   \n",
       "70134   4120282  30301339  662829152    None         UNK       11.2     12.0   \n",
       "70135   4120283  30301339  662829152    None         UNK       12.0     15.2   \n",
       "70136   4120284  30301339  662829152    None         UNK       15.2     17.0   \n",
       "70137   4120302  30302253  662829174    None         UNK        0.0      0.3   \n",
       "70138   4120303  30302253  662829174    None         UNK        0.3      1.1   \n",
       "70139   4120304  30302253  662829174    None         UNK        1.1      4.5   \n",
       "70140   4120305  30302252  662829173    None         UNK        0.0      0.3   \n",
       "70141   4120306  30302252  662829173    None         UNK        0.3      1.1   \n",
       "70142   4120307  30302252  662829173    None         UNK        1.1      4.5   \n",
       "70143   4120312  30303742  662829205    None         UNK        0.0      8.0   \n",
       "70144   4120313  30303742  662829205    None         UNK        8.0     17.0   \n",
       "70145   4120314  30303990  652802870    None         UNK        0.0      2.0   \n",
       "70146   4120315  30303990  652802870    None         UNK        2.0      2.2   \n",
       "70147   4120316  30303990  652802870    None         UNK        2.2      5.0   \n",
       "70148   4120317  30303994  652802874    None         UNK        0.0      2.0   \n",
       "70149   4120318  30303994  652802874    None         UNK        2.0      5.0   \n",
       "70150   4120343  30304039  662829228    None         UNK        0.0      1.0   \n",
       "70151   4120344  30304039  662829228    None         UNK        1.0      9.0   \n",
       "70152   4120345  30304039  662829228    None         UNK        9.0     10.0   \n",
       "70153   4120346  30304039  662829228    None         UNK       10.0     12.5   \n",
       "70154   4120347  30304050  652802882    None         UNK        0.0      0.3   \n",
       "70155   4120348  30304050  652802882    None         UNK        0.3      0.8   \n",
       "70156   4120349  30304050  652802882    None         UNK        0.8      3.5   \n",
       "\n",
       "      TopElev BottomElev MajorLithCode MinorLithCode  \\\n",
       "0       39.25      37.75          CLYU          None   \n",
       "1       37.75      35.25          ROCK          None   \n",
       "2       35.25      33.25          CLYU          None   \n",
       "3         4.0       -2.0          None          None   \n",
       "4        -2.0      -11.0          CLYU          None   \n",
       "5       -11.0      -14.0          CLYU          None   \n",
       "6       -14.0      -17.0          CLYU          None   \n",
       "7       -17.0      -20.0          SAND          None   \n",
       "8       -20.0      -24.5          CLYU          None   \n",
       "9       -24.5      -26.0          SAND          None   \n",
       "10      -26.0      -30.5          SAND          GRVL   \n",
       "11      -30.5      -32.0          CLYU          None   \n",
       "12      -32.0      -40.0          SAND          None   \n",
       "13      -40.0      -44.0          CLYU          None   \n",
       "14      -44.0      -47.0          None          None   \n",
       "15      -47.0      -50.0          CLYU          None   \n",
       "16      -53.0      -60.0          CLYU          None   \n",
       "17      -60.0      -62.0          CLYU          None   \n",
       "18      -62.0      -71.0          SAND          None   \n",
       "19      -71.0      -86.0          SAND          GRVL   \n",
       "20      -86.0      -94.5          SILT          SAND   \n",
       "21      -94.5     -110.0          LMST          None   \n",
       "22     -110.0     -113.0          SDST          None   \n",
       "23     -113.0     -121.0          SDST          LMST   \n",
       "24     -121.0     -143.3          SDST          LMST   \n",
       "25     -143.3     -144.3          CLYU          None   \n",
       "26        4.0       -2.0          CLYU          None   \n",
       "27       -2.0      -17.0          CLYU          None   \n",
       "28      -17.0      -32.0          CLYU          None   \n",
       "29      -32.0      -38.0          SAND          None   \n",
       "...       ...        ...           ...           ...   \n",
       "70127    None       None          SILT          None   \n",
       "70128    None       None          CLYU          None   \n",
       "70129    None       None          CLYU          None   \n",
       "70130    None       None          CLYU          None   \n",
       "70131    None       None          CLYU          None   \n",
       "70132    None       None          CLYU          None   \n",
       "70133    None       None          CLYU          None   \n",
       "70134    None       None          GRVL          None   \n",
       "70135    None       None          CLYU          None   \n",
       "70136    None       None          CLYU          None   \n",
       "70137    None       None          SAND          None   \n",
       "70138    None       None          SAND          None   \n",
       "70139    None       None          SAND          None   \n",
       "70140    None       None          SAND          None   \n",
       "70141    None       None          SAND          None   \n",
       "70142    None       None          SAND          None   \n",
       "70143    None       None          FILL          None   \n",
       "70144    None       None           UKN          None   \n",
       "70145    None       None          FILL          None   \n",
       "70146    None       None          CLYU          None   \n",
       "70147    None       None          SAND          None   \n",
       "70148    None       None          FILL          None   \n",
       "70149    None       None          SAND          None   \n",
       "70150    None       None          FILL          None   \n",
       "70151    None       None          CLYU          None   \n",
       "70152    None       None          CLYU          None   \n",
       "70153    None       None          SAND          None   \n",
       "70154    None       None          FILL          None   \n",
       "70155    None       None          SAND          None   \n",
       "70156    None       None          SAND          None   \n",
       "\n",
       "                                             Description     Source  \n",
       "0                                                   Clay  SAGeodata  \n",
       "1                                         Rocks and sand  SAGeodata  \n",
       "2                                                   Clay  SAGeodata  \n",
       "3                                              No sample  SAGeodata  \n",
       "4      Clay - orange-red grey, mottled; stiff-sticky....  SAGeodata  \n",
       "5      Sandy clay - red brown, grey, green, white mot...  SAGeodata  \n",
       "6       Clay - red green mottled; stiff. fine sand, qtz.  SAGeodata  \n",
       "7      Clayey sand - orange green mottled; stiff. fin...  SAGeodata  \n",
       "8                 Clay - red-green grey, mottled; stiff.  SAGeodata  \n",
       "9      Clayey sand - red green, mottled; sticky. fine...  SAGeodata  \n",
       "10     Sand/gravel - fine sand - fine gravel, clear, ...  SAGeodata  \n",
       "11              Clay - grey red mottled; soft - plastic.  SAGeodata  \n",
       "12     Sand - fine-medium, clear, opaque, fe stained ...  SAGeodata  \n",
       "13     Clay - grey green; stiff - plastic; minor fine...  SAGeodata  \n",
       "14                                            No sample.  SAGeodata  \n",
       "15     Clay - red green; plastic, sticky; minor fine-...  SAGeodata  \n",
       "16     Clay - green grey; stiff. minor silt/fine sand...  SAGeodata  \n",
       "17                        Sandy clay - cream; fine sand.  SAGeodata  \n",
       "18     Clayey sand - cream-pale green grey, fine sand...  SAGeodata  \n",
       "19     Sand/gravel - coarse (1-3mm), predominantly op...  SAGeodata  \n",
       "20     Silt/sand - grey; fine sand (0.05-0.1mm), qtz,...  SAGeodata  \n",
       "21     Limestone - shell and coral fragments (drill c...  SAGeodata  \n",
       "22     Sandstone - yellow grey; interbedded with sand...  SAGeodata  \n",
       "23     Sandstone/limestone - yellow grey; carbonate +...  SAGeodata  \n",
       "24     Sandstone/limestone - grey white and grey; fin...  SAGeodata  \n",
       "25     Clay - steel grey-blue; sticky. (munno para cl...  SAGeodata  \n",
       "26      Silty clay - yellow brown; minor sand component.  SAGeodata  \n",
       "27     Silty clay (gravel) - grey and brown, mottled;...  SAGeodata  \n",
       "28     Silty clay (gravel) - grey and brown, mottled;...  SAGeodata  \n",
       "29     Clayey sand - orange brown; fine (0.1mm) - coa...  SAGeodata  \n",
       "...                                                  ...        ...  \n",
       "70127                                      Gravelly silt  SAGeodata  \n",
       "70128                                         Silty clay  SAGeodata  \n",
       "70129                                         Sandy clay  SAGeodata  \n",
       "70130                                Gravelly/silty clay  SAGeodata  \n",
       "70131                                         Silty clay  SAGeodata  \n",
       "70132                                         Silty clay  SAGeodata  \n",
       "70133                                         Silty clay  SAGeodata  \n",
       "70134                                            Gravels  SAGeodata  \n",
       "70135                                         Sandy clay  SAGeodata  \n",
       "70136                            Sandy clay, clayey sand  SAGeodata  \n",
       "70137                             Clayey sand dark brown  SAGeodata  \n",
       "70138                             Sand clay orange brown  SAGeodata  \n",
       "70139                            Clayey sand orange grey  SAGeodata  \n",
       "70140                             Clayey sand dark brown  SAGeodata  \n",
       "70141                             Sand clay orange brown  SAGeodata  \n",
       "70142                            Clayey sand orange grey  SAGeodata  \n",
       "70143                                Pug hole no returns  SAGeodata  \n",
       "70144                        Cuttings lost into pug hole  SAGeodata  \n",
       "70145        Fill material, gravel, sand, bricks, rubble  SAGeodata  \n",
       "70146                 Clay grey brown, medium plasticity  SAGeodata  \n",
       "70147      Sand, grey to black, medium to coarse grained  SAGeodata  \n",
       "70148  Fill - cobble sized lumps of slag with gravel ...  SAGeodata  \n",
       "70149  Silty sand - grey to black medium to coarse gr...  SAGeodata  \n",
       "70150                                               Fill  SAGeodata  \n",
       "70151                                               Clay  SAGeodata  \n",
       "70152                                         Sandy clay  SAGeodata  \n",
       "70153                                          Clay sand  SAGeodata  \n",
       "70154                                               Fill  SAGeodata  \n",
       "70155                                        Clayey sand  SAGeodata  \n",
       "70156                                               Sand  SAGeodata  \n",
       "\n",
       "[70157 rows x 13 columns]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data           # print the first 30 and last 30 rows"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# add a new column as a function of existing columns\n",
    "log_data['Thickness'] = log_data.ToDepth - log_data.FromDepth"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.frame.DataFrame"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(log_data)     # see what Python type the DataFrame is"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "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>OBJECTID</th>\n",
       "      <th>BoreID</th>\n",
       "      <th>HydroCode</th>\n",
       "      <th>RefElev</th>\n",
       "      <th>RefElevDesc</th>\n",
       "      <th>FromDepth</th>\n",
       "      <th>ToDepth</th>\n",
       "      <th>TopElev</th>\n",
       "      <th>BottomElev</th>\n",
       "      <th>MajorLithCode</th>\n",
       "      <th>MinorLithCode</th>\n",
       "      <th>Description</th>\n",
       "      <th>Source</th>\n",
       "      <th>Thickness</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1769789</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>18.0</td>\n",
       "      <td>19.5</td>\n",
       "      <td>39.25</td>\n",
       "      <td>37.75</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "      <td>1.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1769790</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>19.5</td>\n",
       "      <td>22.0</td>\n",
       "      <td>37.75</td>\n",
       "      <td>35.25</td>\n",
       "      <td>ROCK</td>\n",
       "      <td>None</td>\n",
       "      <td>Rocks and sand</td>\n",
       "      <td>SAGeodata</td>\n",
       "      <td>2.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1769791</td>\n",
       "      <td>30062892</td>\n",
       "      <td>662815923</td>\n",
       "      <td>57.25</td>\n",
       "      <td>NGS</td>\n",
       "      <td>22.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>35.25</td>\n",
       "      <td>33.25</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "      <td>2.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   OBJECTID    BoreID  HydroCode RefElev RefElevDesc  FromDepth  ToDepth  \\\n",
       "0   1769789  30062892  662815923   57.25         NGS       18.0     19.5   \n",
       "1   1769790  30062892  662815923   57.25         NGS       19.5     22.0   \n",
       "2   1769791  30062892  662815923   57.25         NGS       22.0     24.0   \n",
       "\n",
       "  TopElev BottomElev MajorLithCode MinorLithCode     Description     Source  \\\n",
       "0   39.25      37.75          CLYU          None            Clay  SAGeodata   \n",
       "1   37.75      35.25          ROCK          None  Rocks and sand  SAGeodata   \n",
       "2   35.25      33.25          CLYU          None            Clay  SAGeodata   \n",
       "\n",
       "   Thickness  \n",
       "0        1.5  \n",
       "1        2.5  \n",
       "2        2.0  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.head(3)    # print the first 3 rows"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "RangeIndex(start=0, stop=70157, step=1)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.index     # “the index” (aka “the labels”). \n",
    "#Pandas is great for using timeseries data, where the index can be the timestamps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index([u'OBJECTID', u'BoreID', u'HydroCode', u'RefElev', u'RefElevDesc',\n",
       "       u'FromDepth', u'ToDepth', u'TopElev', u'BottomElev', u'MajorLithCode',\n",
       "       u'MinorLithCode', u'Description', u'Source', u'Thickness'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.columns   # column names (which is “an index”)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "OBJECTID           int64\n",
       "BoreID             int64\n",
       "HydroCode          int64\n",
       "RefElev           object\n",
       "RefElevDesc       object\n",
       "FromDepth        float64\n",
       "ToDepth          float64\n",
       "TopElev           object\n",
       "BottomElev        object\n",
       "MajorLithCode     object\n",
       "MinorLithCode     object\n",
       "Description       object\n",
       "Source            object\n",
       "Thickness        float64\n",
       "dtype: object"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.dtypes    # data types of each column"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(70157, 14)"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.shape     # number of rows and columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[1769789L, 30062892L, 662815923L, ..., 'Clay', 'SAGeodata', 1.5],\n",
       "       [1769790L, 30062892L, 662815923L, ..., 'Rocks and sand',\n",
       "        'SAGeodata', 2.5],\n",
       "       [1769791L, 30062892L, 662815923L, ..., 'Clay', 'SAGeodata', 2.0],\n",
       "       ...,\n",
       "       [4120347L, 30304050L, 652802882L, ..., 'Fill', 'SAGeodata', 0.3],\n",
       "       [4120348L, 30304050L, 652802882L, ..., 'Clayey sand', 'SAGeodata',\n",
       "        0.5],\n",
       "       [4120349L, 30304050L, 652802882L, ..., 'Sand', 'SAGeodata', 2.7]],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "log_data.values    # underlying numpy array — df are stored as numpy arrays for effeciencies."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "#log_data['MajorLithCode']         # select one column\n",
    "##Equivalent to \n",
    "#log_data.MajorLithCode \n",
    "##and\n",
    "#log_data.iloc[:,9]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "pandas.core.series.Series"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "type(log_data['MajorLithCode'])   # determine datatype of column (e.g., Series)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "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>OBJECTID</th>\n",
       "      <th>BoreID</th>\n",
       "      <th>HydroCode</th>\n",
       "      <th>RefElev</th>\n",
       "      <th>RefElevDesc</th>\n",
       "      <th>FromDepth</th>\n",
       "      <th>ToDepth</th>\n",
       "      <th>TopElev</th>\n",
       "      <th>BottomElev</th>\n",
       "      <th>MajorLithCode</th>\n",
       "      <th>MinorLithCode</th>\n",
       "      <th>Description</th>\n",
       "      <th>Source</th>\n",
       "      <th>Thickness</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>7.015700e+04</td>\n",
       "      <td>7.015700e+04</td>\n",
       "      <td>7.015700e+04</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157.000000</td>\n",
       "      <td>70157.000000</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157</td>\n",
       "      <td>70157.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>unique</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>5068</td>\n",
       "      <td>4</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>27777</td>\n",
       "      <td>27878</td>\n",
       "      <td>81</td>\n",
       "      <td>42</td>\n",
       "      <td>33598</td>\n",
       "      <td>39</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>top</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>None</td>\n",
       "      <td>NGS</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>CLYU</td>\n",
       "      <td>None</td>\n",
       "      <td>Clay</td>\n",
       "      <td>SAGeodata</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>freq</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>18514</td>\n",
       "      <td>44946</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>18514</td>\n",
       "      <td>18514</td>\n",
       "      <td>25857</td>\n",
       "      <td>62797</td>\n",
       "      <td>4603</td>\n",
       "      <td>70119</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>2.505842e+06</td>\n",
       "      <td>3.018201e+07</td>\n",
       "      <td>6.624491e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>24.942443</td>\n",
       "      <td>30.626594</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>5.684151</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>9.276598e+05</td>\n",
       "      <td>8.068098e+04</td>\n",
       "      <td>2.130462e+06</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>45.435866</td>\n",
       "      <td>48.609957</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>9.943264</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>1.769789e+06</td>\n",
       "      <td>3.002715e+07</td>\n",
       "      <td>6.528000e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.010000</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>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>1.932799e+06</td>\n",
       "      <td>3.014558e+07</td>\n",
       "      <td>6.628129e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.800000</td>\n",
       "      <td>4.000000</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>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>1.999036e+06</td>\n",
       "      <td>3.018487e+07</td>\n",
       "      <td>6.628196e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>11.580000</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>2.800000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>3.967159e+06</td>\n",
       "      <td>3.025487e+07</td>\n",
       "      <td>6.628248e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>25.900000</td>\n",
       "      <td>34.750000</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>7.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>4.120349e+06</td>\n",
       "      <td>3.030405e+07</td>\n",
       "      <td>6.728042e+08</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>610.300000</td>\n",
       "      <td>620.160000</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>300.500000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            OBJECTID        BoreID     HydroCode RefElev RefElevDesc  \\\n",
       "count   7.015700e+04  7.015700e+04  7.015700e+04   70157       70157   \n",
       "unique           NaN           NaN           NaN    5068           4   \n",
       "top              NaN           NaN           NaN    None         NGS   \n",
       "freq             NaN           NaN           NaN   18514       44946   \n",
       "mean    2.505842e+06  3.018201e+07  6.624491e+08     NaN         NaN   \n",
       "std     9.276598e+05  8.068098e+04  2.130462e+06     NaN         NaN   \n",
       "min     1.769789e+06  3.002715e+07  6.528000e+08     NaN         NaN   \n",
       "25%     1.932799e+06  3.014558e+07  6.628129e+08     NaN         NaN   \n",
       "50%     1.999036e+06  3.018487e+07  6.628196e+08     NaN         NaN   \n",
       "75%     3.967159e+06  3.025487e+07  6.628248e+08     NaN         NaN   \n",
       "max     4.120349e+06  3.030405e+07  6.728042e+08     NaN         NaN   \n",
       "\n",
       "           FromDepth       ToDepth TopElev BottomElev MajorLithCode  \\\n",
       "count   70157.000000  70157.000000   70157      70157         70157   \n",
       "unique           NaN           NaN   27777      27878            81   \n",
       "top              NaN           NaN    None       None          CLYU   \n",
       "freq             NaN           NaN   18514      18514         25857   \n",
       "mean       24.942443     30.626594     NaN        NaN           NaN   \n",
       "std        45.435866     48.609957     NaN        NaN           NaN   \n",
       "min         0.000000      0.010000     NaN        NaN           NaN   \n",
       "25%         0.800000      4.000000     NaN        NaN           NaN   \n",
       "50%         7.000000     11.580000     NaN        NaN           NaN   \n",
       "75%        25.900000     34.750000     NaN        NaN           NaN   \n",
       "max       610.300000    620.160000     NaN        NaN           NaN   \n",
       "\n",
       "       MinorLithCode Description     Source     Thickness  \n",
       "count          70157       70157      70157  70157.000000  \n",
       "unique            42       33598         39           NaN  \n",
       "top             None        Clay  SAGeodata           NaN  \n",
       "freq           62797        4603      70119           NaN  \n",
       "mean             NaN         NaN        NaN      5.684151  \n",
       "std              NaN         NaN        NaN      9.943264  \n",
       "min              NaN         NaN        NaN      0.000000  \n",
       "25%              NaN         NaN        NaN      1.000000  \n",
       "50%              NaN         NaN        NaN      2.800000  \n",
       "75%              NaN         NaN        NaN      7.000000  \n",
       "max              NaN         NaN        NaN    300.500000  "
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#describe the data frame\n",
    "log_data.describe(include='all')     "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count    70157.000000\n",
       "mean        24.942443\n",
       "std         45.435866\n",
       "min          0.000000\n",
       "25%          0.800000\n",
       "50%          7.000000\n",
       "75%         25.900000\n",
       "max        610.300000\n",
       "Name: FromDepth, dtype: float64"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# summarise a panda Series\n",
    "log_data.FromDepth.describe()   # describe a single column"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "24.9424428068475"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#calculate mean of 5th column (\"FromDepth\")\n",
    "log_data.iloc[:,5].mean()      "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "24.9424428068475"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#alternate method to calculate mean of FromDepth column (the 5th one)\n",
    "log_data[\"FromDepth\"].mean()    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Count how many Lith Codes there are\n",
    "lithCounts=log_data.MajorLithCode.value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "CLYU    25857\n",
       "SAND    12772\n",
       "SLAT     4069\n",
       "FILL     4020\n",
       "SDST     3207\n",
       "TPSL     2767\n",
       "GRVL     2445\n",
       "QTZT     1891\n",
       "SHLE     1601\n",
       "SCHT     1576\n",
       "ROCK     1572\n",
       "LMST     1571\n",
       "SILT     1436\n",
       "None      874\n",
       "SLST      530\n",
       "MDST      528\n",
       "STON      474\n",
       "SOIL      441\n",
       "QUAR      251\n",
       "BITM      205\n",
       "SPST      190\n",
       "SHEL      181\n",
       "DLOM      169\n",
       "LOAM      143\n",
       "CLYS      119\n",
       "CORL      117\n",
       "MARL      115\n",
       "GNSS       99\n",
       "RUBL       90\n",
       "UKN        82\n",
       "        ...  \n",
       "BASE       13\n",
       "MANM       13\n",
       "PEAT       12\n",
       "TRAV       10\n",
       "COAL        8\n",
       "WKST        7\n",
       "BSLT        6\n",
       "CAAR        6\n",
       "PYRT        5\n",
       "SLAG        5\n",
       "GPSM        5\n",
       "CAVE        3\n",
       "WOOD        2\n",
       "PEGM        2\n",
       "CLYW        2\n",
       "CLYR        2\n",
       "VOLC        2\n",
       "OPCA        2\n",
       "CLLV        2\n",
       "HFLS        1\n",
       "ARKS        1\n",
       "ALDG        1\n",
       "DIOR        1\n",
       "LMSD        1\n",
       "DUST        1\n",
       "FOSS        1\n",
       "REGO        1\n",
       "SCOR        1\n",
       "CALU        1\n",
       "SANN        1\n",
       "Name: MajorLithCode, Length: 81, dtype: int64"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Print the lithcodes, use .index or .values \n",
    "lithCounts"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x9f13990>"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#plot a bar chart of the lith codes\n",
    "lithCounts.plot.bar(rot=90,figsize=(15,5))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0xa410b70>"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2oAAAFJCAYAAADqjN66AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XvcZXPd//H3xww55TAZh4ghh5IcapLoriTdahSlQoV+KXW7lUPRoJOSJilC6VZkOjBEoqZEjrdDMg7jnOM4FBlFdKuEz++Pz3fPta49+/Bde63ruta4Xs/HYx6zT+u7v9fea6/1/XwPn2XuLgAAAABAcyw21hUAAAAAAAxHoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANM3E032yllVbyKVOmjOZbAgAAAEBjXHvttY+6++R+rxvVQG3KlCmaM2fOaL4lAAAAADSGmd2X8zqmPgIAAABAwxCoAQAAAEDDZE19NLN5kp6U9KykZ9x9qplNknS6pCmS5kl6n7s/NjLVBAAAAIDxo8yI2tbuvqm7T033p0u60N3Xk3Rhug8AAAAAqKjK1McdJM1Mt2dK2rF6dQAAAAAAuYGaSzrfzK41s73SY6u4+0OSlP5feSQqCAAAAADjTW56/q3c/U9mtrKkC8zs9tw3SIHdXpK05pprDlBFAAAAABhfskbU3P1P6f9HJJ0taXNJfzaz1SQp/f9Il21PdPep7j518uS+13UDAAAAgHGvb6BmZsuY2QtbtyW9VdLNks6VtEd62R6SzhmpSgIAAADAeJIz9XEVSWebWev1p7r7eWZ2jaQzzGxPSfdLeu/IVRMAAAAAxo++gZq73yNpkw6P/0XSNiNRKQAAAAAYz3KTiYyIKdNn933NvBnTRqEmAAAAANAcVa6jBgAAAAAYAQRqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwE8e6AlVNmT6772vmzZg2CjUBAAAAgHowogYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA0zcawr0ARTps/u+5p5M6aNQk0AAAAAgBE1AAAAAGgcAjUAAAAAaBgCNQAAAABoGAI1AAAAAGgYkonUhIQkAAAAAOrCiBoAAAAANEx2oGZmE8zsejP7Zbq/tpldbWZ3mtnpZrbEyFUTAAAAAMaPMiNq+0q6rXD/a5KOdvf1JD0mac86KwYAAAAA41VWoGZma0iaJun76b5JerOkM9NLZkracSQqCAAAAADjTe6I2jGSDpL0XLr/IkmPu/sz6f6DklbvtKGZ7WVmc8xszvz58ytVFgAAAADGg76BmpltL+kRd7+2+HCHl3qn7d39RHef6u5TJ0+ePGA1AQAAAGD8yEnPv5Wkd5rZ2yUtKWk5xQjbCmY2MY2qrSHpTyNXTQAAAAAYP/qOqLn7we6+hrtPkbSLpIvc/QOSLpb0nvSyPSSdM2K1BAAAAIBxpMp11D4j6QAzu0uxZu2keqoEAAAAAONbztTHBdz9EkmXpNv3SNq8/ioBAAAAwPhWZUQNAAAAADACCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYSaOdQUwZMr02X1fM2/GtFGoCQAAAICxxIgaAAAAADQMgRoAAAAANAyBGgAAAAA0DIEaAAAAADQMgRoAAAAANAyBGgAAAAA0DIEaAAAAADQMgRoAAAAANAyBGgAAAAA0DIEaAAAAADTMxLGuAOo3Zfrsns/PmzFtlGoCAAAAYBCMqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMNMHOsKoJmmTJ/d8/l5M6aNUk0AAACA8YcRNQAAAABomL6BmpktaWa/N7O5ZnaLmR2WHl/bzK42szvN7HQzW2LkqwsAAAAAz385I2r/kvRmd99E0qaStjOzLSR9TdLR7r6epMck7Tly1QQAAACA8aNvoObh7+nu4umfS3qzpDPT4zMl7TgiNQQAAACAcSZrjZqZTTCzGyQ9IukCSXdLetzdn0kveVDS6l223cvM5pjZnPnz59dRZwAAAAB4XssK1Nz9WXffVNIakjaX9PJOL+uy7YnuPtXdp06ePHnwmgIAAADAOFEq66O7Py7pEklbSFrBzFrp/deQ9Kd6qwYAAAAA41NO1sfJZrZCur2UpLdIuk3SxZLek162h6RzRqqSAAAAADCe5FzwejVJM81sgiKwO8Pdf2lmt0qaZWaHS7pe0kkjWE8AAAAAGDf6BmrufqOkzTo8fo9ivRrQ0ZTps3s+P2/GtFGqCQAAALBoKbVGDQAAAAAw8gjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYQjUAAAAAKBhCNQAAAAAoGEI1AAAAACgYSaOdQWAXqZMn93z+Xkzpo1KGQAAAMBoYkQNAAAAABqGQA0AAAAAGoZADQAAAAAahkANAAAAABqGZCJABhKSAAAAYDQxogYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADTNxrCsAjBdTps/u+5p5M6aNQk0AAADQdIyoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAw0wc6woAyDdl+uy+r5k3Y9oo1AQAAAAjiRE1AAAAAGiYvoGamb3EzC42s9vM7BYz2zc9PsnMLjCzO9P/K458dQEAAADg+S9nRO0ZSZ9y95dL2kLSf5vZhpKmS7rQ3deTdGG6DwAAAACoqG+g5u4Puft16faTkm6TtLqkHSTNTC+bKWnHkaokAAAAAIwnpdaomdkUSZtJulrSKu7+kBTBnKSVu2yzl5nNMbM58+fPr1ZbAAAAABgHsgM1M1tW0lmS9nP3J3K3c/cT3X2qu0+dPHnyIHUEAAAAgHElK1Azs8UVQdpP3P1n6eE/m9lq6fnVJD0yMlUEAAAAgPElJ+ujSTpJ0m3u/s3CU+dK2iPd3kPSOfVXDwAAAADGn5wLXm8laTdJN5nZDemxQyTNkHSGme0p6X5J7x2ZKgIAAADA+NI3UHP3yyVZl6e3qbc6AEbalOmz+75m3oxpo1ATAAAAdFMq6yMAAAAAYOQRqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDEKgBAAAAQMNMHOsKAFj0TJk+u+9r5s2YNgo1AQAAeH5iRA0AAAAAGoZADQAAAAAahkANAAAAABqGQA0AAAAAGoZADQAAAAAahkANAAAAABqGQA0AAAAAGoZADQAAAAAahkANAAAAABpm4lhXAMD4NGX67L6vmTdj2ijUBAAAoHkYUQMAAACAhiFQAwAAAICGIVADAAAAgIYhUAMAAACAhiFQAwAAAICGIesjgEVWHZkj68o+2a8cMlgCAIAyGFEDAAAAgIYhUAMAAACAhiFQAwAAAICGIVADAAAAgIYhmQgANEQdCUlIagIAwPMDI2oAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAkEwEADENCEgAAxh4jagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMARqAAAAANAwBGoAAAAA0DAEagAAAADQMH0DNTM72cweMbObC49NMrMLzOzO9P+KI1tNAAAAABg/ckbUTpG0Xdtj0yVd6O7rSbow3QcAAAAA1KBvoObul0n6a9vDO0iamW7PlLRjzfUCAAAAgHFr4oDbreLuD0mSuz9kZit3e6GZ7SVpL0lac801B3w7AMCiZMr02T2fnzdj2ijVBACARdOIJxNx9xPdfaq7T508efJIvx0AAAAALPIGDdT+bGarSVL6/5H6qgQAAAAA49uggdq5kvZIt/eQdE491QEAAAAA5KTnP03SVZI2MLMHzWxPSTMkbWtmd0raNt0HAAAAANSgbzIRd9+1y1Pb1FwXAAAAAIBGIZkIAAAAAKAcAjUAAAAAaBgCNQAAAABoGAI1AAAAAGiYvslEAAAYC1Omz+77mnkzpo1CTQAAGH2MqAEAAABAwxCoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDkPURAPC8ReZIAMCiihE1AAAAAGgYAjUAAAAAaBgCNQAAAABoGAI1AAAAAGgYkokAANBDHQlJSGoCACiLETUAAAAAaBgCNQAAAABoGAI1AAAAAGgYAjUAAAAAaBiSiQAAsAhoSlITEqMAwOhgRA0AAAAAGoZADQAAAAAahkANAAAAABqGQA0AAAAAGoZADQAAAAAahkANAAAAABqGQA0AAAAAGoZADQAAAAAahkANAAAAABqGQA0AAAAAGmbiWFcAAACMP1Omz+75/LwZ00apJgDQTIyoAQAAAEDDEKgBAAAAQMMQqAEAAABAwxCoAQAAAEDDkEwEAAAskupISEIZ5crIKaeOMgAwogYAAAAAjUOgBgAAAAANQ6AGAAAAAA1DoAYAAAAADUMyEQAAADRKU5KaNKWMnHIoo5ll5JbTCSNqAAAAANAwlQI1M9vOzP5gZneZ2fS6KgUAAAAA49nAgZqZTZD0bUlvk7ShpF3NbMO6KgYAAAAA41WVEbXNJd3l7ve4+9OSZknaoZ5qAQAAAMD4VSVQW13SA4X7D6bHAAAAAAAVmLsPtqHZeyX9p7t/JN3fTdLm7v6JttftJWmvdHcDSX/oUexKkh4dqEL1ltGkulAGZSwKZTSpLpRBGeOljCbVhTIoY1Eoo0l1oYzxXcZa7j65b0nuPtA/Sa+T9JvC/YMlHTxoeamMOVW2r6uMJtWFMihjUSijSXWhDMoYL2U0qS6UQRmLQhlNqgtlUEbOvypTH6+RtJ6ZrW1mS0jaRdK5FcoDAAAAAKjCBa/d/Rkz20fSbyRNkHSyu99SW80AAAAAYJwaOFCTJHf/laRf1VQXSTqxIWXUVQ5lUMZ4KaOuciiDMihj9MuhDMoYL2XUVQ5lUMZIlyGpQjIRAAAAAMDIqLJGDQAAAAAwAgjUAAAAAKBhCNQAALUxs/eY2ZJjXQ8AABZ14z5QM7MjRrj8Fczs0LGui5lVShxT8r3WHIX3WCLzdVuMdF2awMzuH2CbyWbW/2KL5ctdvO4yu7zPfmb24tF4rz712Ges61CnGn6/H5B0v5n90MzeZmYT6qhXUe7+Xscx1cy+XLWMPuUvU+K1549kXUaLmb2scPsFbc+Ni2N23czsNWb2tg6Pv9PMXl3Xe2S+rvKx2cyWLp5LzGwDM9vfzN5dpdzxromfn5ktW+Y4mFnmVnWWN5bGLJmImf1AUrc3d3ffM6OMG7s9lcrYOKOM69z9Vf1el1HOSyR9TtKLJf1c0qmSvixpN0mnufu+I10XM/uFpH3c/b62x98i6Rh33yizjK47hbu/M6OMuj7Tz7n7Qo0kM1te0jnu/qbRqIuZ3avhn4kV7ru7vzSjjAPaHnLFVesvd/d7q9Qvlf+Au78k43Um6QuS9lH8HYtJekbSce7+pQrvb5K2lvR+Se9w91X6vH4tSY+7+9/S/a0l7SjpPknHu/vTGe95nKR3S/qDpNMknenujw1Y/7Uk/Z+7P5oaiq+XdLe7n52xbV37+0HufmS6/V53/2nhuSPc/ZCMMnpey3K0fr9mtpykdymur7mJpHMUx8HLqpRbKD93f6/jb6nr+11d0mqSbnT3p81sZUn7SfqQu2c1as3senffrEId1pA0xd0vT/cPkLRsevpUd78ro4xa97H2z7fM5506IN8mqRX43SbpPHd/Jmf7VEYdf8+1ki6XdKWkK9z9wdz3zyh7a0kHuftCQVjb6y5R7Evz2h5fV9KJ7v7mAd9/Q8XveFdJf3P3qRnbVD42m9llkvZ09zvT3/B7ST+RtKGk37v7wRllfFTSJakMk3SypJ0kzVN8VtdllLGnpEnu/vV0/4+SXqg4fx7k7ieU+bs6lL9tKmfbzNdX2udrOiYuLenf7v7vdH8DSW+XdJ+7/6xEOXtLmi5pGcXn+aSkr7n7dzK3nyDpfZJWV3wGN5vZ9pIOkbRUzrGyQ5u31Ta72N1/nFmPnp2b7l66I71o1EZZOvhlh8fWVJy4cntgn1N8qKdK+oWkfwxQjwlmtqJiJ1mIu/81s5wfSrpU0lmStpP0O0m3SNrY3R8epbrMknSxmZ0k6UhJkyUdo/hc98isw1Hpf5P0PUkfydyuqGP9B/AfZvYVd18wImlmqyqu3XdWTe+Ro/3EtJji4PBpSddnlvHCDo9NkXSomX3R3WcNXj1JPYLrNvtJ2krSa1oBopmtI+kEM9vf3Y8u86Zm9lpFcPYuSZMk/bekAzM2PSNt8zcz21TSTyV9VdGo/44y9jt3/4SZ7asIEHeR9GUzu0bRMDjH3f8v82/4nKQPSXIzmyXpLZIukTTNzN7o7vvllFODXRS/W0k6WPGZtGynOPn08zpJDyg+g6s12G+x8u/X3Z+QNFPSTDN7kaT3SDrOzCblBFg5b5H5ujqO7xPMrNUw61TGE/0KMLP9JB0q6S5JLzCzb0n6puK8UWa0Y/lePeIZjaSvKxq7LR9TpJFeWtJhitHQfurex9q3zyovjdhcLOkhxXHYJG0v6RtmtrW7/ymzLnX8PXtK2lLSOyR9NTWmr1AEble6+7X9CjCzN0o6QUOdvV9V/IaWkvSVjDq8qD1IkyR3vyv9BrOljqtd079nJK0laWqn8jup6di8orvfmW7voejo+YTFbJprFcfIfvaVdEq6vaukjSWtLWkzSd+S9B8ZZXxccfxtecTdV7eY3n2+4jvry8zeLOm7Gvp+j1D8/k1532+d+3xV5yn2+VYQfZXiuLK9mb0mM4j+rOI38yZ3vyc9to6kb6XzxOEZ9ThJ0ksUQfyxZnaf4vc83d1/nvm3HNXhsUmSPmhmG7n79IwyZivOScVjhyva4CsrP6bpqBHp+dOXc4ikN0g6WtJJOT3qaduXKX6A75B0qyJoO79E78K/JP1RnQ/O7u7rZJYz1903Kdz/s6Q13f1fOdvXVZc02nSkpG0kLa44AHzPB/iiB+29NbNHFEFjR+7+ycxylpR0pqQ73P0AM1tP0q8lfd3d/yezjMclde3Fz+ktLZS1mGKE9EBJN0g6wt1vzd2+S5mTJP02p4erw6jcgqckHerukzLKuF7Stu7+aNvjkxW/m6zv28y+oghW71ecfM+WNMfd187c/sbWiLeZHSXpOXc/KH3GN+SMhncoc6Kk/5R0uKT13T1rKoWZ3SppU0Vj9X5Jq7r7U6m8G/qNRJvZM5Ke6vSU4ne7XGY9Fvze2n97ub/F1MO4rYYaJbMVDZxbcuqQyqjl95vKWlERpO0qaT1JZ+UGvjXt73UcU/8l6c9tZbROyu7ufaeKpn3s9e7+19T7epekN7j77/pt21bOXxSjk93+ng/32b599Kq4z/2vu/dtuNa0j1UeUTOzUxS/z2PaHv+kpFe7e1bnZB1/T4cyV5H0XkXH2Nru3rehlo7NByoavW9TNEIPc/dvZr7nXe6+btnnOrz2SknLK44Bs9Jo1L25x/YuZZY+NredI65QnPd/nu4Pa2/1KOMGd9803T5V0tXu/q10P3c/u9bdX124f4i7H5FuX+PuudNBr5e0v4a+3x9K+lyrPpllnKKK+7yZPaU4/iz0lPJnot3k7q9Mt7+sGHH871YQ3XquTxl/kLSJu/+z7fGlJM119/UzyrhZMRjyXGovPipp3RKDI73KnqD4WzYdYNspkj6j6PQ91t2Pq1KXsRxRk5m9XNHLuJmip+/juQFWi7vfrpjK9QUz21mx838tlZfj1kGCkU7aem4flrS0pXm3mT23ddRlQ0mbK3oYpkpaRfE9/3uAsgaN4v+h6PGqxN3/aWbvkjQrjXS8TtJ+njEdrWC+pG9UqYfFPPkPKw6yl0vawd3vrlJmS2q45fbgdhqVa8k92C/eHqSlesy3cmvL9lJMazlB0i/Td1Vmfyn+zW9W6h1NB9wSxaTC4liyS/r3d8UxIdc/U8fQ02Z2t7s/leryjJnldBjdVNMxpH36Rbfnuhfg/qyit/M8i7U/u0q6xMy+VOJkUen3m0afdkzv/SpJ5yoaaBeX7DCqY3+v45haRxn/bJ0D3P1+M7ujbJCW3NcvGOujPcnLNoXbWaMuNe1ja5jZsYrjQOu20v3VM8vYwt0/1KF+x6ZGYJY6/p50DN9EMUKwlaT1FcH9KYqGeW5dfptunmlmRyo6rnP9NnWgfbb4OzOzwyRdVKKc+ZLWULQdJku6U4O3Baocm29MnXh/lLSuYvRKZrZCibd/zsxWk/SYYl8vjlwtlVnG8sU7hSBtMWX+ZoY29UvS7Z+b2fwyQVpSxz5/r2Jgo4ri/vBmpfa2x5Tu57ILaQvS0mP/KFHG0+7+XKusdFytHKSl8p4t2xZJgwmHSnqtot35SU/TQ6sYs0DNzH6qCCSOUjSAn5W0XOuDyQxsWvP+d1FMo3oslVWmIV+X5RWNm+I325r/7JKyRuaqMLPvKxpGe7v7VSlIPEzSXDPbz937LkRPIzwtC00byvxe/uLuM0tWv1NdWj3qv5d0kKT/lbR26/HMnsYn3f3SilW5VzH94xjFiMsmZragN89LzMluZzEdInfu/m8GbNwV9Qo8skaxk1UlvVXRqDnGzC6WtJSZTczsbLnIzM5QTOFYUakhkU6quaPpa2to7cRESacr1sfdUeLvkKQVLKaTmeIY1JpaZmo7SY+wTczsifS+S6XbrXpkZ1FMjc1pis9liqRjJZXZR6v+fu9VTE8+QbFuYKATlbsf1u05y0xq0CDFYESSVi7eLzFKWXVa6pNmtn7rN9I6nlvMTPl7biE17GPF6dFz2p5rv99Nr6UOnUa4u6rh73lCQx1XX/ShKXtlLG9m7TM83lFoE/VcSyfpU5K+L+kuM7shPbaJ4vPMXsLg7jtYzMzZSdJhFlPbVjCzzd399zll1HRs/qhi6uIUSW9tdaApOqNzO+I/r/j7J0g6tzVKajHN9J7MMs43s8Pd/bNtj39JKXjMtIINn7ZsxfuZ7Yg69vmnvS2PwQDqCKIfNLNt3P3C4oOpTfRQZhkvs6FcFSbppel+mdHBTjMzVpS0u2L5Ul9mtpEiQHuFYkbbnqkDqBZjmUxkngoJGdQ2rcTzpqRcquh1PUMxRW5YEJETVJjZh9z9lLbHVlQkOhjVD6dqXcxsf8Uw67Ntj79S0nc8b1pLK3FGlalCv3P3ypm7zKxnz1uvhlyhjJ+5e6UsR2m6Qa/EN317uM3spg5lTJL0J0m7p5HhfmVcrxS0ekrCUZaZPSup0/oAk7Sku5fO2JimHGyvOCm/XtKF7v7+PtuYpJ0VyRXOcPc/psc3k7Syu/8m433v09D0nNy1gp3K+UGv5939//XZfsFUmLFmZjMlbaSYIjzL3W8eoIxKv18zW7rQqCo+/hJJu3halD9AuYMkNah8fDezPd39pC7PHeXun84oo+eUpNzA2GLNxELfqcU0nV3c/ScdNiu+bjtFEPIVDXUkvlqx9GBfd/91Rh3q2MfeozQSX3bbQhn3KNYJL/SUpCM9I8lTKqeOv2c3xYyPzST9U3GcvkrSVe7+58wyftTjaXf33TPLWUfRYJSkW9z9HjN7sQ+4fski6c3Oit/dSzwviU8tx+Ye5Z/u7jtnvG5xxXn3hV5IZpI6sddx95syylhGEQC/RtLc9PCCANjdszo4+pxnctsRlfd5Mzve3StlKraYnriv4vx9srvPTY9vKeml7t5rX26V8QrFNO7LFYMcrviMt1LMWuobJFmspewqJyDt0OZtJRO5RNLhnrcG+VnFOtfZikGn9npkLxfoWP5YBWp16BDsLXhK+UHF5xUNxdtTr9p5ih/hM5LeX5iKkFOf9mw8typGQXLXy1Wqi5kt1hoG7vDcij5gRryybIQz4KT3WMYzk0V02T4705KZLd8tMLJYOHtNRhntBxRXjFxk/w0WUy0+KWlvSV/OORiONotpb+8edEQmt8GZUc4V7j4q6XktMpz1ypSauyazWyateZ453ddiykhrn+p0TMxaL9eh3JcqGmq7eEb22MJ2KynW6uyqmNJ2dk5gU9i+UlKDOo/vXcq/3zPWqPXYfknFSMNP+744Xr+cImHP6ooppRcoMrh+WrF2ZYeMMjZSzFBY0JhXNPKyApQ69jEzO1vRIDtPscb1/LI90FU7WArl1PqbMbNlJW2hmAa5m6TFMhvQL3H3B8q8V4k6Ze+naZ98obvPb3t8FcVapNsq1qXysTn37zGzXysa/U+3Pb6JIqnJlBLvWQyAb/Walj+UUcc+b2bvUGSevS/d/7xi9PQ+RWdN6SzUKSDeSNIf3f2REtstqUhG9grF7+0WST8p24GTRnBfofj93uYpOUlVZvZad78643UfUu82QLUZZu4+Jv/p2dJGAAAgAElEQVQUQ9HLFu5voUgm8gbFQaJq+atnvu4WDQWseyky6kyQ9HJFCtjc93uxYurDJYp55ccoskD+QdKLR6Muih7S13Z4/COS7smsw1qSli/c31qxHmR/SUtklnGTpBvT/zcV7j8k6dmy36NiiuwS6f7KimxJf8rcfmtJdyim9fxYMW1ijqIH592ZZcxRZKBqf3xbSQ9klrG0Yn1Y6/4G6TN91wD79oaS/qZIZftE6//MbSe1/Vuxtc+VrMMBiuH99sc/oVhH2G/75RTr0o5XTKG0tO19ihNoqfp0KD/reym8foKklQr3l0i/wdsytt2j178SdbhM0nrp9rqKGQLHSbpQ0oyqn8kAn+FqimQIv1eMEnxB0isztnuhYtrIeYrpRd+Q9OAA73+l4pj4ucLncm/JMmo5vte1nxX2tVYygT8rUpfnbnuOYt3TxxQzSS5QnGc2rfh3LCnpvaO8fy2XfiO/VswsOEGRYGXU6lDz37Okov0yXZGF+hHFue+7mdvfrQi4J4xA3bL3U0UW0IXOjYqMoCeMZl16lHF/5usOT8fPpQuPvUkx+rFtifebqFjXdWD6t72kiQPUeyNFJs85kq5Jt/seU2veF25sfR7p77hDMar+EcXAQk4Z35X0inR7ecWgxE2KqZC7juLfslw6Dt6jmKp8drr9U0nLjdZ+NuJ/55i9caxNO6hw/950cLtAcR2FUfmAJV1fuH2WpI8V7l9X4v1OUYcGqmIEZOZo1EUx7WyuIq3+JMU0jKvSTrtGZh2uVgosFZnwHlXMfZ8p6fsDfhdTFCfhOyV9osR2+ykWNl+lCEL3kPQXRSC8Wu5nmg7ML1AkOHhC0WtUpv4fVWR5nFx47P1pn904s4xejfCvlqjLnulz3EeDBVj3pgPZvYV/j0r6reL6Srnl3KwOgXv6nG/M2H5EGpyF8rMPsIrpdH9TNBYvVQT3DyoO+q+qWI+1Srz2psLtL0v6drq9RPG5kf6X9veLFCfwwxWZ8O4tsf0/0uf4HxoKkrI6ijrsI/crgvktBylHNRzftXDnRuvfi1QiAFU04r+raCSepZRwquTfU9xHJijWtw7UsakBA0bV2OlUKONF6VgwV+WCitobvorOwTXTv6zGeHrvv6bj6JcV6dyXL/m+KyjOCde39ve6/pU8Ht7a47lbRqsuivX2nf69WtJDJd7vUMX69mUVI0f3K0blc7ev3AmfytlBce7+cDqmbpJu36kY9csp4x0qnFMUa/DmKkbX184sY27h9smSPlO4n3tMvKVwez9JP0+3V1XhmNunjD0lHVi4/6CGOp7/K7OMUyR9UTFq3XrM0ufywxr21dyO+B+kz7LTv5Mq16NqARU+gOuLB8HWl5s+5MtH8QP+neJgPzkdaNcuPPeHEu93e4/nssqpoy6Kk+/hihP4A4pFuGU+txsLt49STImR4tphfRvgbWWtl35Ityl6axYvuf2tiqkWUpw0n1ZkPSpTxnVt9+8ecH/aTdFj1BpluF3lApvKjXDFCMOpivTxA/82upT9bkXih9J/T5nnunweAzU4U507/dtJ0vwS5dysSOkrRUPgXyrZ6FSsT3mPYn2dFCfiU3OPQ2mb4m/vCkk7Fu7PzSyjOMr6ZOH+U5KeySzjaUVDZGrhsewASdFovzp9rodIemmZ7dvKWl7RmLlA0anwmKTNS2zf65ja9ZjdVkanzo3Wv9yZCg+m3+9urf1cJUcH0zbtx7PszsTCNpUCRtU88qsY1W91Djwg6ZjM7So3fFM5B0v6fOH+/Ypj/e2SDi7xmdYyEqZYp/O4onPwOkVbKaej9jjF+sP2f8cpc8ZFKqfrLIJez7W9rvKxWTH63fVfyc/0AEVgdYfSsb7EtqeoYid8ev1cdWgzKDqyc4/vdYyG3agIWhdTzGApHue7BultZRQ7wGYrLh6+0HN9yrhGce2/YdspRqYvyyzjzkGeK/Gd5XYq7NTh3/7p8y09m6T931im51/Mh6/d+owUE8LTHO+qPPN1+yoSkUyWdLQPXQT47RpabJ2jjmw8ddSltR7kBMU1HHY2szmef+HuXmnTcy9EWlcGnDpSWteRaUnu/iMz+6fipHm/pK3c/S8l6lHcHwdNZ/sFd7+gxHtmc/efWVyAMpuZreJtC+XTOoYcCzIBeqTBvdfdnyzz/uqdYviXJcp52t3vSnW5LtUlO3OsmX1dcdK8QdJnzOyXinWERygajrkqZ9Jy92Ep7dOawb0VoxW5f9OLFceRb6bv8wzFNRlz63C0pKPTmo5dFRd3fbGZfUaxRi0785vH2tCTJZ2c6rKzIstoVlIDRadKt2NqVpIDr3D9qIKzFCP6O0t61szO0WApz1uZQSUNyw6atZ7KzB5UHL9OUPRmP5n29zJZEvtdiLjvBWKtnks4fEkxfW1e4bG5ZnaRYjT2nMxy3qvhFz7+i7tvltbLXqq48HQ/P5B0gpkdPeC5TpJkkY3wOEVw8G1J2anO1TtbZm4mTUl6xDpkeLTItDq/yzbtKh+b3X3rzPfqysx+oaFEEZMV1w/7pg1l0sy5jmotl4FQdFTP61DOPMu/NI4XfqvvVozWXCvpWjPbO7OMYxTnqicUgfccaUEir9xsi4+b2faKc9VWitGxVp6G3MseLNbWhvqptCDNfm4ZVbPgFveRTmXnXrLkrEJ5xetCz1BcD7GSsQzUljCzF7YaZ55Sx1ukhc1KRd1jEb8pphH05bFQ8GUdHv9ViZ1FitS6nbILmmIe7YjXxcx+qwgY3+Lu95rZoYopcteY2dfc/cSMavRKm567wHOuhjLgbC5p82KM5/kZcOpIaX2php80ivddGSmYbShjoymm/bxI0sUpcHXPuzhzHels35kWAndU4nNdSOocWazEJl+XNNvMPqXh2eOOVIzE9lOpwSn1XjhtZjtl1KFlZRt+ceVli/e9/2UgpknaLJ1gVlRModzYy6fo7pWOOuczXSDtV/sp1oqdKuk1uR0LHtfZO0HR8FxDMTX0ETO7TRFoHZJZzj2K7IJfscg8+37FeqSsbHxtf8+ykv7u7sdKOrZDcp5ubnP3bsfUrEZneq/HU9AoM9taEWTMU4yM972chLvva2b7KabV7qr4/SxnZu+T9CvPzBznGRdO7qOOgLGOTqd7Vf0SDnU0fFvbFJM6fSs99myJNsCrFIHmHDP7b3e/ssz7S5KZ/VhxGZ/d3f2Gfq9v54WEBen34j5Ywq0DJZ1hke24dT3FqYpjyS6Zdal8bDazI1rHGjPbdsBOyqO63C6jrstA/NvM1vS2ZGrp+JJ7/WBL3+1TiuvCfafwXFa72d1PNrPfKNb7zy089bCkrAQ8io6/Y5VmGPnQtcu2UbT7ctRxfborUjKULxc7eMzsc4rZFDl67RfZ+4zVcF3ormXnd17VKzWE3qL4Y+5Pj62lOHBf6O59L1JsNaU87lF+mUxJP+hTl9wfwMB1MbN3dRoJMLNVJX3D3T+Q8T7d0qa/QdIPPC97VS3fy0h/v7n6NQo9LwVsHelsi5/HYWq7cGjO59EWkLSsKOmdko539+/1K6NQ1tsUveetLIA3K6Y+9U3xPdJK/na/0ONpd/cv9dn+Wnd/deH+De6+aWZVc+qXndbeIsvipxS/4ZMlHecDXsqhQ9kbSNq53+fRY/uJir/jxyW22Vuxj7VmWTypWMP8ne5bDdt+jmLU5bG2x9+q6I3OSTV+tWIq7J/MbFPFOqSvKqbb/dvds69RVShzccX6sF0UQflKmdv9SnGdzHll37NQhmkoYHy7oiNxT2UGjCmgeFjRIfEZxXTSp1LnwKXuvknPAqKMypdwMLO5ioyZnRq+v8jsQJOZ3aFIjPDvtsdfIOlmd18vp5y0zWsU03TnKUbDWh1Pr8rY9r8UiUcWapSZ2U7FXvs+ZRwsaZn00N9V4vdSKGdlRXbR1rH9FsX5ITujX4+yczM2Xtf63Iq362JmW7n7FRmvq+syEDsqOjKP0PB09NMV68R+nlHGhxWjNU9IesTdt0uPbybpKHffptf26bXtn6NLetRHKONoj3p8R9Jfve36dGZ2uCKx18czylhOMWL1KsUooSsCpesVl054vEL9yhyLiteFPkNtKfo9f0Zb5/LHKlCTJDP7uGKnKx5UZrj7CWNWqQIzeyDnRJ5RTtZBtmpdzGzZbidaM3upl0wnmxol75f0PkUP6M/c/bgyZXQoc62cwKYuXQKTBTJGS0Zc7gmjbZvr3X2zAd6rPShxRYKWyzzjmjIjxeLC9a3Rgj9V7Ymq8bfb9/ILZva4Yt1OyxvS/VYjLWd6TXuZA6W1N7P/U0xP+oEiqBkmZ383s4Pc/ch0+71eSB1f7OXusX23FPKfUqzF6JtCPpXzWUWK833S6FxrWsm3JF3t7odnlPHRVJdtPaUbN7P3K0b6dnD3G3ttn15/Y6vRbzEq/py7H5R6fm/IDQi6lL2mYn1HVvCbRuAOVyTNOHLAkahieYsrEl/sqsyAsdDptKqi8650p1NbeYPu65UbvqmcI9Lfsk8reLS4dtbxkh5294Mzy2lNW7xIbdMWy557O5Sd01Fb+fdSKGszxcj3LV4xHX+HsrOOzXUEahbTV9+n2K/Oc/ebLabsHSJpqZxzaJ2d8BaXBfiUhqejP6r1G8osY3Wl0TBPl2OymPG0eHunRZftL+7w8CTFevldPXM0t9BR20qLf6uiU+BXmdv3uj7dR73EUgiLy8dsqPSZDvp7q3AsmqeFLxW24LpsnnGpsJ7lj2WgtqASMZRrZb6YtN3rFRct/GG6f6Zih5PiQnUXVaxXpWvk1FlO5oH6bsXi5zMKjy0p6bOKnvC+PYNmtr6GLir7F0mnS/q0u+dONWqV8zrFjn6Zuz9iZhsrftT/UaYBnUaR9lVkFpMiMcmxre88Y/tiYPIxSf9TfN7zLpp9rzpcX2eoiKxRxsonjLbyau9hLPn+la4dZmYHK04sX0r371dkXlxcsUA7Z11Ir/IH/s1ZyQsrpwZauwUHa3e/NPN9XyjpXYrOkfUVa8p2dvc1StT9i+r+vXhOQNCrgZSz31lMp3tMka11G8WI7RKKbKvZU7os1n9s4m3X1EmBwlx3Xz+znN0U1w17q2Kk8eOStssdlTKzm9z9len2dYpj7G/S/QVBXK5BGwOF7ZdRZDXbTtKPNDwgGLjjyczOcPf3Zb62UkO+jn09lVNHw3eCInD/iGLxvxTJq06S9NmcTiMbmra4d5l9vEQdczpq6/q9fF7SBxXB72sVWYmzZ1pklJ87ovagpG8qvtf90+0FMjudTpH0EsUlRl6r+H5fJ2l6biDfp/zKnfAl36/nsdfdy+RVaC97qqRvuvsbMl77UUV76iANrX+cqliT9X3PW2bTKmvg69OZ2X8qkjOd2fb4BxQjjn2ny9Z1LBpJY7ZGrdNIhw1fx5RzwjlMce2llg0kfUgxQneI0vqqPvVorT9a6ClJuYkR+r5N1ouq1+Wtko5PP6L/Uuz8RykW8+cGArcrUtm+w1OCBTPbP3NbpdfXklzBzHZXrLM5QLEOyhRD3F83M+UEa8VAzMx2zAnMOmhvqC+mCLo+rcyEBIqTfuuEcayZ1XrCyGVm5/Z6vsToT5nF6Z1UXsBf52/Xql1YeQXF5S++ncr6vWLhuislScr0iGL/+Kwi862b2btKbC93/2K35yymZeWwLrc73e9knUJg833F5R/WLNsRJ8XC8g6P/cPy10LVkQio17rdvuvT0ms7NQbWGbAx8G/FBZpfoLhmXZmEE71skfOitob8kWY2SEO+8r4uSSkg290qrMnySP4x3cwOU6wflqS73L3X2qR2V0jazTv0fNfUmM/qUa/j96LozNjUYzrrixTXQyz1/dZ0bP6eYv9uv13GVMV64edSp/WjiqyPD/fZLtfRinWffVn3pBWSss+9cxSdEa31tcXjsSvWjA7E3edYfiK//SW93odP6bsojbJdrrgWX09m9kF3/7G732Nmq3lhVpGZ7ePux2fU4zB1TlxzoeIYm7OusZZjUbs0yreLYpRyo36v72Usk4n0+tHlDvMt5+63Fu7f6ZEBR2aW2yP/bsWBo31+7lqKOfh1yP17tq/0JtET8TYzO1ARcD0s6T/d/ZYSxeyk2LkuNrPzJM1SZqBZUFdyhb0Va0PmFR67yGIx8izF9X/KGGj4uNWos5jqtJtisfUNkqa17X+9VD5hmNmTGkpq0kq+oXTfPSMBhyI4fEDSaYoU6gNlTfIa1gd69QX8nX4vJmkNRUdNFjO7UrGweZak97j7nRaZ8OZlFnGQhi+yX0LxfS+jmIL4004bdXBIKucESaea2emZ23XVPjqohTsdOvEutzvd76SOjJ6S9KCZbePuFxYfNLNtlJmdzOpJBLSfhtbtvt6Hphuuqlg8nqOWxoCZbacYWThXcY2/MskM6lK5Ia+a9nUbWsO4TLpfek2Wpam+KaB5mZec6itJ3nu5RlZjvobgptvv5c3Kz+YnRbblp6Q496XzXlmV21UDdqq2e9rT9MDUHrmjxiBNKnf+HDSZSdGnFG20fyjOV2d7ZiKifiyy6ua2kcw7rLtK+0vuWx4gqbVe+ThFJ3zLhxVTj/tZ2tOU9rZ6PJxmHuSo7bybOu92VnTIbazocN510PJaxixQ6/UjtMiOlWNYxjx3L2ZdzO21OVrSId62bsrMJqfneqWZLb6+cg9Sex1SuSspRhv6/oAsFusfqFgYvrdiofixZra3u2elkfVIRnJ22sl3VPScrGJmJygOCudnFPOPVs+euz9mZn8YIEiTIhCf16GO8yzWwYwKi3UcH1Z8Fpcr1raUnQNd+YThbanXB7SqpG0VB4/3KzI0nVYymK9jZG5ZM1u81eh191NSuS9QfpbUBb8XW3g9ZZke7PmK4G4VxUjYnSoX1C/hwxdiX56C+7+UOFnUlta+4uigNJSRs1OHQE5msU2tYkbP5JOSzjGzyzV8DdJWimto5ThZ8Zt9TIUAsox07J1VfCwdl2/IOS4ndTUGDlV0JuR2EA3TY+qUKf8SDJUb8j329YMUF8/tu6/b0JqsN3nbmiwzm+T5a7J2Uax1kyIRR7FjZTuV6PTpVtXM11UNbj4h6dyKvxdJemnb8b11v8ya28rtKhue8Xkhnpfl+GVm1lqHaoq/5cZ0+znPSHzTR/Z5wgtT4NPnoE5BRp8yWr+btRW/mwvT7JwjPH9tWadlC5MUv6V9M6vyhJlt4m1TjC2mIud2ylWduSFJS5rZRG+bnpzabFkdvjUdiz6atl1DkUzkI5LOqamzoRlr1NpZ/hzmXyiyJM1ue3x7xZXNp2WUcXO3YUkrrE3IKKeOzIBbKOb4/lVxUeQfSVpJMdVud3c/r8/2NymmjR3qQ6mkt1f05Pwsp2ewS7mTFFPVdnb3vkPrVlNyBWvLppf7XNvrigH0uorrqCyQ06NuMVf+GcX1RxZarOsZ12Izs6cK722K9R13aegzyanHkor1NesqLlp5cvsBqowUFLVShX/JSySKsUht3nVkzvusy7LuC/i/Lekhz1jAbzWtp0xlLa/oqdxV8fmuoBiN/n3PDWPbu9x93S7P3e2ZWcG6bN9Ka/++nHLaRgdnFUYH67gWWBYbMNFNh3LWVewj62v4GqQ7Jf0xp7PEIvnHlorLntyouOj0FZKu6tQj3KWMSsfltrJajYFdJK2nWGuW1RhI2y8taS/FPnqTInNl9jHAIpFAa4RRalt76xnXrmo7vptiCnOp43v6blfx4dOdNlYcY9/oGZchsPrWZC3YX9v33Tr25RLtmV8qgpsb2x6fqriGZs/gxsyOVxyPN9Dw38tP2j+jPuW01twupdhHn5N0t1Ka+n7H9lRG5XaV1ZPluNO5YMGsC3d/e0YZvTrh13f3F/Qro1DWFxQBtSmOH88oMvOWzqRrZq9QHEd2k3SQF3IT9NmuPZt2K6HYNZ6Z1dMiP8RPFDNGip0Ce0j6oLtfnlFGpbXQ6XUzFJ0b+3iaoZPaEccpLqzed+lBTceipxVrsj/lQ9elu8crJhFZUH5DA7XcrEDrKkYErtTwazltKWn7zEi4VyOr63M5yoyGpdfPUfTeLa+Y4/s2d/+dmb1MMerR84RhZm9w98s6PL6kpNmekbq1DlZfcoVicDPsKcX6jr6jFWa2nnr0Unpah9enjFbWp/aEIlI0TPquu6spkD9dMSrwv4rU3ve5e24PWLGcFyimp+6quGbXuYqg748lypigoZG5jVVyZM7qWcD/nOKz2NOH1lNWPjja0IWVd5HU98LKZvYTSZd42zodM/uYord/4KkPAxxDzlGsRz1X0qnufmXZz6Rqh0DuSTajnEoN17ZtWtNRt1RM/32d4tpoG2ZsW+m4nMqo3BhI21Q6BpjZ5pIecPeH0v09FB0U8yR9MSd4reP43uO7fY3iu+27DMBipsYGXZ673TtcP6/La+toMFZuzFcNbsxsX8UxazVFp9VpPkBikzQa8RXFLJL7NRTYnKL4zvqOTNfdrqopYF5o1oVnrIOqI9hL5eyvmOW0l7vfmx5bR0PXEjw6o4x1FN/xDoo2zSxJvywTiLeVN9DIXtp2FUVG3WKnwLc9c6ZQoX1X7LiWyrXvJiqy4HZqR3wuc1+t41hUTBC1imJU7UM5cUyOpgZquT1Q7T2uUuwsdym/x/U0SRd1aGTtqUhXvHNmnSv3ulrh+ktmdpu7v7zwXN8DlcX1Pr6ryNzzTHpsFUnfkLSBu+cmE6jEzHZQj+QKXlgD0KecOoKbyo09i4s6t/dCz1dMcbu33/Z1seHZ5yZK+n3ZBrGZzVRcG+fXihGXm2uoV+mRuXQQfFDS44qA4E2K6TC3K7/B+C7FSWtLxTqZWYqMU7WNHlnG5SQsrjn0c0n/0vAOoxdI2tHd/5z5XrWM3FiF0cG0fdVgoJWtrSPPzExYR6984fXLK4KzrdL/K0i6yTNSa1c9LqfXVW4MpNdXOgZYZK18i7v/1eLamLMUPfybSnq5u78no4zKx/eaRlwuVEz56rQm63OeMTqYXv+sIjmLKUaRWuv+TNKS7t53SmhNHYK1BDfpvLlL+rekYpRtVk7Hddr+aMV1Cw/wtLbUYqnBUZKecve+S1PqalcVths0PX9tsy5SeQMFe2nb6xWXCXm07fHJks7PPI48p+g8O0dxLbVhDfgSx9ZKI3vW4cLdZdXUvqujHVHbeSZts4aG9rmlFcsWKk2fHsusj63ECAs9pcy5pYreyEPc/eS2sqem53J6XPdTrMn6gGIIV4qe1yUUWbpyHa+hXteL1NbrqmhI9lPMzNSecSonon61oqF3fepde6ViweaRknbP2L4utSRX6PRDLTvCIGlKe+MolT3HzKZkltEpE9Jakg41sy+6+6wOzw/TZ393z1u3U0zS8IzlL9ot2k3RKFlf0icLZZRdP9RpZO5YSX2ngSb/o2gw/sMi4czBGmowniipb4PR61lPKeuz3k5xMfBe9XhE0papgdjqMJrt5S8PUscxRB7Tnk+WdLINjQ4eY2Z9RweTDQvBwEmKJBhlTFD8ZgbaQQt6rYfLOkeY2YmK7+RJxRTdKxUdWY/13HC4qsdlqftx6Jp+DZY2VY8BEwqNl50lneiRkfAsM8sdfanj+F75u1VNa7I8czSzjzrWul9jZh/tEtxc22WbhaQ6fE3S1ywuo3CyYtpg7t+5vWIUcMG+7e5PWFxM+3ZFm6mfutpVVdWRxbpTsJc1TbjN4u1BmhSjWWkUM8eXNHTMaW+X5M662F/S6yW9pn1kz8z2zxnZU3RMtkahz3L3nXLee1hluwRiFrNtdtHQCFkvldsRquc88xrFTIWH3f1Bi6Uhf1TEWLmZNLtz90X2n6Sbezx3U8mytlZ8wZ+Q9OYB6nJD4fZtbc9dn1nGs4pekicVPRxPFO7/u0Rd9lU0Lh5U9HyO9vdyTdv94wu3f1einC0kXaJo/G8m6WZFJstHFNdByinjrkGeyyx7kqTrRvFzbe0f7fvIk5KeGOXveKbiBHy4pI0G2H5u4fa3Fb1frfs3VKjXJMX1XS4qsc18xUjYgYr1lG8s/hvFz7TyMaRP+Wtlvu66XvfLbl+hvqcpLnza/viekk7PLOM8RUrrUxRru16pNJOkRD0qH5frOg5VPQakY+jEdPt2SW8oPpdZh8rH95q+2/0U6+P2Uoz2fCNtv2TJ73fJVNbxqayJA+yrldsiihG5KxXnvG+kf5cq1r6sWqIuiysCw58ozpenK0b2c7e/Y5Dnurx+4HZVa59O/9p/d1nnO0VQeLpipPN7ius63luyHs+l72HdwmP3DLCPdD0u1nHMlLRf5uuul7RSh8cnK7+ten2n2yXru5wisDpecWkpS/vJfYpEHDllVG5H1HQsuk7SpHT7DYrkPzspZsacWfW7Hcv0/HWoo1dOkuTuF0u6uEJdKve6esWePTNbQdGT9lpFtqq3S/q1me3rFS/+XdKKxTvuvk/h7uQS5dQxwlBLL2UnHtOHqo4alHm/Onp+61J1ZG6CDWVr2kbRQGoZ+LjkMVrwP2q7sHkftWTCrEHlY0jV0cGklfVR0kBZG+v6TVTulXf37dJv9BWKKbKfkrSRmf1VkVDkCxll1PG7q2u0pGpdTpN0qZk9qtjH/jfVY13F5Rty1HF8r2PEZQ3FFLSXS5qrCHIeVkw3KrNmZ6aGpvq+XbGvlF37W7kt4jFFeksz21oxPV0qMTJvZq1j2DTFKPgsxXqosteWu9XMdve265Sa2QcVwX22Ku0qryHLsdcz66KOSxZJw4+rRbnZdPs5QDGTrJ86Rva8y+0yfqTIxnuVYo3ZgYrf/w6ev7ayjnZEHceiOmYqdNXINWq56p4HXbEulee511CHeyR9R9IxPrRGbdP02H1eIalByXrUklyhprUhqygufPi0OvwIvcI1VdJUt896RiZMDGdmhyoaRY8qFv++yt09NRhnuvtWY1SvgTNh1vDedayVqZSNsw4WqdGzMlzk2VMAAAJWSURBVCpmlldsuN4yaKdTWjuwlVKyKUkvcvcVem9Vj5E8Dg1Qly0UCSfO96FMaetLWtbdr+u5seo7vqdtKn+3ViFRTNq+jrW/Y94Wscjoeapi3dTAvz8zW10xi+UfGj6ldCnFvpqddKqJrGQW68J2rWBvV8WFpWeqxBT7kWb5Sfi6rvfLXQvY51yV05nX/ruboNQW8BLX3KyzHVHlWGRmNyuuLfmMmd2u6CC5rPWcV7zg9aIeqDXm5NcEZraGuz/Y5bmFenNHsB51JVeonI2r8PoqP8JOGb0mKYa3d3f3Ur2MCFUbjDXXpXImzCawitk4n2/M7JOKBvxWilGTKxQ9uFcopqQ912PzkahPLUHnWKrr+F5jfQZOFJO2r3RuSds879oihTW3pthXL+yzybgxaLA3kiw/CV8ryFroKY3SoEKqR+XfXdpuzNsRI93xvEgHai3Ph5Pf85ENT65Q+ntpwihlqkf7Yn9XJDQpO6UEDWQjkAmzCcZydLApzOybStdO85SSHvWoenyv4f3bE8X8TrFGrkyimPaGa/FcM0iCJdoiGDHWJwmfuy8yy5nq/N01wUgGjM+LQA0ABmWR8rh1wigeEBfVE8bzYnQQ6CWtF1pJkSDlSsVI6c1OowbA8wiBGgA8TzxfRweBTtoSxWyp2PezE8UAQNMRqAHA88TzbXQQyDGWiWIAYCQRqAEAgEVK0xLFAMBIWGQWHgIAACRTJJ0paX8SxQB4vmJEDQAAAAAaZrGxrgAAAAAAYDgCNQAAAABoGAI1AAAAAGgYAjUAAAAAaJj/DwngvEAjc8WbAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0xa3a4dd0>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Plot a bar chart of the lith codes for the rarer lithologies\n",
    "lithCounts[(lithCounts < 50)].plot.bar(rot=90,figsize=(15,5))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xn8VVW9//HXWxC1HAD9aggommiZv5wIKW9lYYg2oN00vJXgxfg9yib73cq63XDI5nJosChJHFLJMqm8IlHqo64DOKQCDjgQBALKoITpRT+/P9Y6sj2ec75n4/d8B77v5+NxHmfvtdfee629z9mfvdbeZx9FBGZmZs3aqqsLYGZmPYsDh5mZleLAYWZmpThwmJlZKQ4cZmZWigOHmZmV4sDRYpLmSzqiq8vRHUj6kqSfNZg+UdKfO7NM3Z2k4yQtkbRe0sFdXR4zcOB4RSQ9JunIqrSXHPwi4g0RcWM7yxkmKST1bVFRu4WI+FpEnAK9p84d4DvAJyJi+4i4q3pi/gw+kwPLekk3NFqYpCMl3SnpHzkgnVCY9s487SlJj0iaXGcZP8/7bp9C2ickzZP0rKSLq/JX9vX6wuu/CtNPkPQ/kjZIurHG+qZKekDSC5ImVk07QNIsSU9IetmP0vK6r5O0RtLjkn5Q+cxJ2lfStZJWSVqdl7Nf1fyn5fnWSZomaZuq6Z+W9Gjengsl7ZvT3yHpXklrJT0p6RpJgwvzbZOX91Re/merlnuKpEV5W10vafda+6KrOHD0Aj4492h7AvPbyfPeHFi2j4gx9TJJ2h/4BfCfwE7AQcAdedrWwDXAT/K0DwLfk3Rg1TL+BXhtjcUvA74KTGtQzv6Fcp5dSF8NnAd8o858fwU+DtxZY9r/AjOASXXm/RGwEhhEqu/b87IA+gMzgf2A3YDbgWsrM0o6CjgdGA0MA/YGzixMPyWv993A9sB7gCfy5AXAURHRH9gdeAi4sFCuM4DhpP37DuDzksbm5b4d+BowDhgIPApcUad+XSMi/NrMF/AYcGRV2kTgz7XyACOBecBTwArgezn9b0AA6/PrzaSg/mVgMemDfwmwU2G5J+VpTwL/VbWeM4Crgcvyuk7J674FWAssB34A9CssL0hfqIeAp4GzSQeIW/IyZhTzV9V5MXBoHv5wXtb+efwU4DeFcl3WoM4TgT+TzrLXkL4wRzfY/qcDD+fyLgCOa5D3DOCXeZs8DdwL7At8MW/fJcCYQv6dgIvytvo76aDYJ097LfDHvO2fAC4nHRSL+/w/gHuAdcBVwLZ1ylVzPwPb5O0SwD+Ah5v9DDbYBr8Azq4zbbe8rlcV0uYCJxbG+wJ3AW/MefepsZyvAhdXpQ3L+fu2U75TgBsbTP8zMLHOtH2AqJG+EDimMP5t4Cd1ljEwl3Pnwvb6WmH6aODxwn5bAoxuYrtvA3wdWFBI+3vV5+1s4Mo8/B3gh4Vpu+dyvbaZ/dwZL7c4Otf5wPkRsSPp4DMjp78tv1fOyG4hHUQnks5G9iad0fwAXjxz/BHwIdKZ1E7Ai83gbBwpePQnHdieB04DdiEdpEez6cyrYixwKDAK+DwwNa9jKHAAcGKdet0EHFGoyyOkM7vK+E015qlVZ4DDgAdyOb8FXCRJddb7MPBWUv3PBC6TNKhOXoD3ApcCA0gHwFmkA8Bg4CzS2XbFdGAj6YB0MDCGdGADEOlAsDvwetL2OaNqXSeQtudepAPtxDplmkiN/RwRz0bE9jnPgRFR6yy/4vLc3XJDdQuhyiiA3IWyXNJlkgYCRMQK0lntyZL6SHoz6Wy4eM3pNODmiLinwToaWSxpae7q2mUzl1HW+cB4Sa/KXUVHA9fXyfs2UmB4Mo+/gdTaqfgrsJuknYEh+XVA7vJ7VNKZkl48pkraQ9Ja4BnSicS3cvoA0menetlvqMyaXxTGIX0Hu4eujlw9+UU621tPOouvvDZQv8VxM+kAt0vVcoZRdUYGzAE+Xhjfj9Qs7wt8BbiiMO1VwHO8tMVxcztl/wxwTWE8gMML43cAXyiMfxc4r86yJgEz8/BC0gG2cva0GDikUK7LGtR5IrCoql4BvKbJ/XE3MK7OtDOA2YXx9+Z9V2lF7JDX1Z909v0ssF0h/4nAn+os+1jgrqp9/uHC+LeAH9eZt+5+LuyXl53ZF/IfDmyXt9UXgccptH6q8j6Xy7YvKUD9Cri8apusIAXMjcBHC9OGAovIrd565aJ2i2N7YET+7O5GOqGZVWPeVrQ4Xp8/yxtzmS8GVCPfEFIroNjCehgYWxjfOi9jGPCWPPz7/JkZBjxY3GaF+QYCXwBGFbZlUGiFAu8CHsvDo0kt2TfmffsT4IVi2br65RbHK3dsRPSvvHj5WXzRJNKX9n5JcyW9p0He3UkH3YrFbPri7U5qJgMQERtI3SZFS4oj+ULg7/KFuKdIfajVZ30rCsPP1BjfntpuAt4q6TVAH1LXzOGShpFaA3fXma+WxysDuV7UW6+kkyTdnS9AriWdkTU6k62uzxMR8XxhvLKuPUkHieWFZf8E2DWvd1dJV0r6e96Wl9VY7+OF4Q316kDj/dyuiPhLRDwTERsi4uukk5e31sn+DPDziHgwItaTPgPH5Dq9jrTfTgL6kc5+Py/p3Xne84CzImJdM+WqKuP6iJgXERsjtWw+AYyRtGPZZZWRz/5nAb8GXk3aRwOAb1blawNuAH4UEcVrCeuBYhkrw0+z6fPyrYhYGxGPkT4jx1SXIyJWk1qw1+brjeurllcZfjrnnwNMIQX2xaRg/zSwtLmat54DRyeKiIci4kTSAeibwNWSXk06+6i2jHQAq9iDdNa0gtTvPqQyQdJ2wM7Vq6savxC4HxgeqavsS7y0ObzZImIR6eD4KVJL52nSgXMyqfX1Qq3ZXsk6Je0J/JR0ENo5B+376Jg6LSG1OHYpnBTsGBGVroSvk8r/xrwtP/wK1ttoP2+OaFCWe6i/3Q8AHoiIWRHxQkQ8QDqbPjpPHw18O594VILiLZL+bTPLSINydpSBpLP7Stffk8DPKRzcc7fRDaQW8zlV888Hil1/BwIr8nIeILXgmv0c9yV973eMiDWk73D1sl+8CSIifhgRwyNiV1IA6Uv6fHcLDhydSNKHJbXlA+nanPw8sIrUFN27kP0K4DRJe0nannR2eFVEbCQ19d8r6S2S+pG6v9r7Eu5Ausi9Pp9dfqzDKpbcRDqIV65n3Fg1Xq1WncuoBNxVAJJOpoP6gCNiOelg8l1JO0raStJr890ukLblemBt7jf/3CtYXaP93FDuQz9cUj9J20r6HOms+i91Zvk56RrG3pJeReo++V2edhcwXOmWXEl6LekuoUo//L6kg9tB+QWpa+uaXJa+krYltTj75PJUbns9TNJ+eTvuDFxA6pJal6f3yfP2BbbK825dqGe/PF3A1nn6Vnma8rR+eXxb5VtmI+IJ0g0WH8vl6w9MqNQpt3hmAX+JiNNrbK9LgEmS9s8B5sukrq5Ka/gqUqtsB0lDgI9Wtqek9xfq3AZ8j9Sdubqw7C9LGpC/jx+tLDvX4YBctz1I1xrPzwGnW3Dg6FxjgfmS1pMv2kXEP/OH8BzgL7lrZBTptsZLSddFHgX+CXwSICLm5+ErSWcuT5PuyHm2wbr/A/i3nPenpA99R7qJdEC9uc74S9Spc9MiYgHpusstpLPz/0P9A+bmqHTZLCDd4XU16UYESIH6ENIdU78ndYVsrrr7uQk7kFqSa0j982NJd6E9CSDpQ5KKZ7HTSAes20hdIM+SWolExMPAv5MO6k+R9t+vSHeWERErI+Lxyisv8omIqHTZfJnUfXM6qQX2TE6DdHJwPemzd19eb/FGi4/k/BeSutmeIX1GK27IaW8hHUSfYdPNFXvm8Uo9nyG1Biren7fLKtI1mo2ki/wAxwFvIgXT4m9M9sh1vp50fepPeXstJnUhVXyCdAKxjPQ5/AWbbkceXKjzvaSTpOMK804hXUNZTNrW387rA9g2L2s96RbhW0h3TnYbyhdjrAfLZ6prSd1Qj3Z1ecxsy+YWRw8l6b1Ktxi+mnTf972ki2hmZi3lwNFzjSM1kZeRfoE6Ptx8NLNO4K4qMzMrxS0OMzMrZYt8+N0uu+wSw4YN6+pimJn1KHfccccTEdHWXr4tMnAMGzaMefPmdXUxzMx6FEmL28/lriozMyvJgcPMzEpx4DAzs1IcOMzMrBQHDjMzK8WBw8zMSnHgMDOzUhw4zMysFAcOMzMrZYv85bhtec6d/WDL13Hau/Zt+TrMtgRucZiZWSkOHGZmVooDh5mZleLAYWZmpThwmJlZKQ4cZmZWigOHmZmV4sBhZmaltCxwSNpP0t2F11OSPiNpoKTZkh7K7wNyfkm6QNIiSfdIOqSwrAk5/0OSJrSqzGZm1r6WBY6IeCAiDoqIg4BDgQ3ANcDpwJyIGA7MyeMARwPD82sycCGApIHAFOAwYCQwpRJszMys83VWV9Vo4OGIWAyMA6bn9OnAsXl4HHBJJLcC/SUNAo4CZkfE6ohYA8wGxnZSuc3MrEpnBY7xwBV5eLeIWA6Q33fN6YOBJYV5lua0eukvIWmypHmS5q1ataqDi29mZhUtDxyS+gHvA37ZXtYaadEg/aUJEVMjYkREjGhraytfUDMza0pntDiOBu6MiBV5fEXugiK/r8zpS4GhhfmGAMsapJuZWRfojMBxIpu6qQBmApU7oyYA1xbST8p3V40C1uWurFnAGEkD8kXxMTnNzMy6QEv/j0PSq4B3Af+3kPwNYIakScDfgONz+nXAMcAi0h1YJwNExGpJZwNzc76zImJ1K8ttZmb1tTRwRMQGYOeqtCdJd1lV5w3g1DrLmQZMa0UZzcysHP9y3MzMSnHgMDOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKwUBw4zMyvFgcPMzEpx4DAzs1IcOMzMrBQHDjMzK8WBw8zMSnHgMDOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKyUlgYOSf0lXS3pfkkLJb1Z0kBJsyU9lN8H5LySdIGkRZLukXRIYTkTcv6HJE1oZZnNzKyxVrc4zgeuj4jXAQcCC4HTgTkRMRyYk8cBjgaG59dk4EIASQOBKcBhwEhgSiXYmJlZ52tZ4JC0I/A24CKAiHguItYC44DpOdt04Ng8PA64JJJbgf6SBgFHAbMjYnVErAFmA2NbVW4zM2uslS2OvYFVwM8l3SXpZ5JeDewWEcsB8vuuOf9gYElh/qU5rV76S0iaLGmepHmrVq3q+NqYmRnQ2sDRFzgEuDAiDgb+waZuqVpUIy0apL80IWJqRIyIiBFtbW2bU14zM2tCKwPHUmBpRNyWx68mBZIVuQuK/L6ykH9oYf4hwLIG6WZm1gVaFjgi4nFgiaT9ctJoYAEwE6jcGTUBuDYPzwROyndXjQLW5a6sWcAYSQPyRfExOc3MzLpA3xYv/5PA5ZL6AY8AJ5OC1QxJk4C/AcfnvNcBxwCLgA05LxGxWtLZwNyc76yIWN3icpuZWR0tDRwRcTcwosak0TXyBnBqneVMA6Z1bOnMzGxz+JfjZmZWigOHmZmV4sBhZmaltPriuHWhc2c/2CnrOe1d+3bKesyse3CLw8zMSnHgMDOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKwUBw4zMyvFgcPMzEpx4DAzs1IcOMzMrBQHDjMzK8WBw8zMSnHgMDOzUhw4zMyslJb+H4ekx4CngeeBjRExQtJA4CpgGPAYcEJErJEk4HzgGGADMDEi7szLmQB8OS/2qxExvZXlNrP2+f9eeq/OaHG8IyIOiogRefx0YE5EDAfm5HGAo4Hh+TUZuBAgB5opwGHASGCKpAGdUG4zM6uhK7qqxgGVFsN04NhC+iWR3Ar0lzQIOAqYHRGrI2INMBsY29mFNjOzpNWBI4AbJN0haXJO2y0ilgPk911z+mBgSWHepTmtXvpLSJosaZ6keatWrergapiZWUWr/3P88IhYJmlXYLak+xvkVY20aJD+0oSIqcBUgBEjRrxsupmZdYyWtjgiYll+XwlcQ7pGsSJ3QZHfV+bsS4GhhdmHAMsapJuZWRdoWeCQ9GpJO1SGgTHAfcBMYELONgG4Ng/PBE5SMgpYl7uyZgFjJA3IF8XH5DQzM+sCreyq2g24Jt1lS1/gFxFxvaS5wAxJk4C/Acfn/NeRbsVdRLod92SAiFgt6Wxgbs53VkSsbmG5zcysgZYFjoh4BDiwRvqTwOga6QGcWmdZ04BpHV1GMzMrz78cNzOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKwUBw4zMyvFgcPMzEpx4DAzs1IcOMzMrBQHDjMzK6XV/8dhZmZZZ/xPe2f8R3tTLQ5JhzeTZmZmW75mu6q+32SamZlt4Rp2VUl6M/AWoE3SZwuTdgT6tLJgZmbWPbV3jaMfsH3Ot0Mh/SngA60qlJmZdV8NA0dE3ATcJOniiFjcSWUyM7NurNm7qraRNBUYVpwnIt7ZikKZmVn31Wzg+CXwY+BnwPNlViCpDzAP+HtEvEfSXsCVwEDgTuAjEfGcpG2AS4BDgSeBD0bEY3kZXwQm5XV/KiJmlSmDmZl1nGbvqtoYERdGxO0RcUfl1eS8nwYWFsa/CZwbEcOBNaSAQH5fExH7AOfmfEjaHxgPvAEYC/woByMzM+sCzQaO30r6uKRBkgZWXu3NJGkI8G5SSwVJAt4JXJ2zTAeOzcPj8jh5+uicfxxwZUQ8GxGPAouAkU2W28zMOlizXVUT8vvnCmkB7N3OfOcBn2fTHVk7A2sjYmMeXwoMzsODgSUAEbFR0rqcfzBwa2GZxXnMzKyTNRU4ImKvsguW9B5gZUTcIemISnKtxbczrdE8xfVNBiYD7LHHHmWLa2ZmTWoqcEg6qVZ6RFzSYLbDgfdJOgbYlvSjwfOA/pL65lbHEGBZzr8UGAosldQX2AlYXUivKM5TLMtUYCrAiBEjXhZYzMysYzR7jeNNhddbgTOA9zWaISK+GBFDImIY6eL2HyPiQ8Cf2PTjwQnAtXl4Jpu6xD6Q80dOHy9pm3xH1nDg9ibLbWZmHazZrqpPFscl7QRcupnr/AJwpaSvAncBF+X0i4BLJS0itTTG53XPlzQDWABsBE6NiFK3BJuZWcfZ3MeqbyCd+TclIm4EbszDj1DjrqiI+CdwfJ35zwHO2YxymplZB2v2Gsdv2XRBug/wemBGqwplZmbdV7Mtju8UhjcCiyNiaQvKY2Zm3VxTF8fzww7vJ/0eYwDwXCsLZWZm3Vez/wB4AulOpuOBE4DbJPmx6mZmvVCzXVX/CbwpIlYCSGoD/sCmR4eYmVkv0ezvOLaqBI3syRLzmpnZFqTZFsf1kmYBV+TxDwLXtaZIZmbWnbX3n+P7ALtFxOckvR/4F9Kzo24BLu+E8pmZWTfTXnfTecDTABHx64j4bEScRmptnNfqwpmZWffTXuAYFhH3VCdGxDzS38iamVkv017g2LbBtO06siBmZtYztHdxfK6kj0bET4uJkiYBzf51rNVw7uwHu7oIZmabpb3A8RngGkkfYlOgGAH0A45rZcHMzKx7ahg4ImIF8BZJ7wAOyMm/j4g/trxkZmbWLTX7fxx/Iv0Bk5mZ9XKb+38cZraZOuP61mnv2rfl67Dey48NMTOzUtziqMF3PJmZ1ecWh5mZldKywCFpW0m3S/qrpPmSzszpe0m6TdJDkq6S1C+nb5PHF+XpwwrL+mJOf0DSUa0qs5mZta+VLY5ngXdGxIHAQcBYSaOAbwLnRsRwYA0wKeefBKyJiH2Ac3M+JO0PjAfeAIwFfiSpTwvLbWZmDbQscESyPo9unV8BvJNNfwA1HTg2D4/L4+TpoyUpp18ZEc9GxKPAImBkq8ptZmaNtfQah6Q+ku4GVgKzgYeBtRGxMWdZCgzOw4OBJQB5+jpg52J6jXmK65osaZ6keatWrWpFdczMjBYHjoh4PiIOAoaQWgmvr5Utv6vOtHrp1euaGhEjImJEW1vb5hbZzMza0Sl3VUXEWuBGYBTQX1LlNuAhwLI8vBQYCpCn7wSsLqbXmMfMzDpZK++qapPUPw9vBxwJLCQ9uuQDOdsE4No8PDOPk6f/MSIip4/Pd13tBQwHbm9Vuc3MrLFW/gBwEDA93wG1FTAjIn4naQFwpaSvAncBF+X8FwGXSlpEammMB4iI+ZJmAAuAjcCpEfF8C8ttZmYNtCxw5H8OPLhG+iPUuCsqIv4JHF9nWecA53R0Gc3MrDz/ctzMzEpx4DAzs1IcOMzMrBQHDjMzK8WBw8zMSnHgMDOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKwUBw4zMyullU/HtV7i3NkPdnURzKwTucVhZmalOHCYmVkpDhxmZlaKA4eZmZXiwGFmZqU4cJiZWSktCxyShkr6k6SFkuZL+nROHyhptqSH8vuAnC5JF0haJOkeSYcUljUh539I0oRWldnMzNrXyhbHRuD/RcTrgVHAqZL2B04H5kTEcGBOHgc4GhieX5OBCyEFGmAKcBgwEphSCTZmZtb5WhY4ImJ5RNyZh58GFgKDgXHA9JxtOnBsHh4HXBLJrUB/SYOAo4DZEbE6ItYAs4GxrSq3mZk11inXOCQNAw4GbgN2i4jlkIILsGvONhhYUphtaU6rl169jsmS5kmat2rVqo6ugpmZZS0PHJK2B34FfCYinmqUtUZaNEh/aULE1IgYEREj2traNq+wZmbWrpYGDklbk4LG5RHx65y8IndBkd9X5vSlwNDC7EOAZQ3SzcysC7TyrioBFwELI+J7hUkzgcqdUROAawvpJ+W7q0YB63JX1ixgjKQB+aL4mJxmZmZdoJVPxz0c+Ahwr6S7c9qXgG8AMyRNAv4GHJ+nXQccAywCNgAnA0TEaklnA3NzvrMiYnULy21mZg20LHBExJ+pfX0CYHSN/AGcWmdZ04BpHVc6MzPbXP7luJmZleLAYWZmpThwmJlZKQ4cZmZWigOHmZmV0srbcc3MeoRzZz/Y1UXoUdziMDOzUhw4zMysFAcOMzMrxYHDzMxKceAwM7NSHDjMzKwUBw4zMyvFgcPMzEpx4DAzs1IcOMzMrBQHDjMzK8WBw8zMSnHgMDOzUloWOCRNk7RS0n2FtIGSZkt6KL8PyOmSdIGkRZLukXRIYZ4JOf9Dkia0qrxmZtacVrY4LgbGVqWdDsyJiOHAnDwOcDQwPL8mAxdCCjTAFOAwYCQwpRJszMysa7QscETEzcDqquRxwPQ8PB04tpB+SSS3Av0lDQKOAmZHxOqIWAPM5uXByMzMOlFn/5HTbhGxHCAilkvaNacPBpYU8i3NafXSX0bSZFJrhT322KODi229wZb0Zz5bUl2s++kuF8dVIy0apL88MWJqRIyIiBFtbW0dWjgzM9ukswPHitwFRX5fmdOXAkML+YYAyxqkm5lZF+nswDETqNwZNQG4tpB+Ur67ahSwLndpzQLGSBqQL4qPyWlmZtZFWnaNQ9IVwBHALpKWku6O+gYwQ9Ik4G/A8Tn7dcAxwCJgA3AyQESslnQ2MDfnOysiqi+4m5lZJ2pZ4IiIE+tMGl0jbwCn1lnONGBaBxbNzMxege5ycdzMzHqIzr4d18ysFN9a3P24xWFmZqU4cJiZWSkOHGZmVooDh5mZleLAYWZmpThwmJlZKQ4cZmZWigOHmZmV4sBhZmalOHCYmVkpDhxmZlaKA4eZmZXiwGFmZqU4cJiZWSkOHGZmVooDh5mZleLAYWZmpfSYwCFprKQHJC2SdHpXl8fMrLfqEYFDUh/gh8DRwP7AiZL279pSmZn1Tj0icAAjgUUR8UhEPAdcCYzr4jKZmfVKfbu6AE0aDCwpjC8FDitmkDQZmJxH10t64BWsbxfgiVcwf3expdQDXJfuaEupB2xBdfnsK6vLns1k6imBQzXS4iUjEVOBqR2yMmleRIzoiGV1pS2lHuC6dEdbSj3AdSmrp3RVLQWGFsaHAMu6qCxmZr1aTwkcc4HhkvaS1A8YD8zs4jKZmfVKPaKrKiI2SvoEMAvoA0yLiPktXGWHdHl1A1tKPcB16Y62lHqA61KKIqL9XGZmZllP6aoyM7NuwoHDzMxK6bWBo71HmEjaRtJVefptkoZ1fimb00RdJkpaJenu/DqlK8rZHknTJK2UdF+d6ZJ0Qa7nPZIO6ewyNquJuhwhaV1hn3yls8vYDElDJf1J0kJJ8yV9ukaeHrFfmqxLT9kv20q6XdJfc13OrJGndcewiOh1L9IF9oeBvYF+wF+B/avyfBz4cR4eD1zV1eV+BXWZCPygq8vaRF3eBhwC3Fdn+jHAf5N+1zMKuK2ry/wK6nIE8LuuLmcT9RgEHJKHdwAerPH56hH7pcm69JT9ImD7PLw1cBswqipPy45hvbXF0cwjTMYB0/Pw1cBoSbV+iNjVtpjHsUTEzcDqBlnGAZdEcivQX9KgzildOU3UpUeIiOURcWcefhpYSHqSQ1GP2C9N1qVHyNt6fR7dOr+q73Rq2TGstwaOWo8wqf4AvZgnIjYC64CdO6V05TRTF4B/zd0IV0saWmN6T9BsXXuKN+euhv+W9IauLkx7clfHwaSz26Iet18a1AV6yH6R1EfS3cBKYHZE1N0vHX0M662Bo91HmDSZpztoppy/BYZFxBuBP7DpLKSn6Sn7pBl3AntGxIHA94HfdHF5GpK0PfAr4DMR8VT15BqzdNv90k5desx+iYjnI+Ig0pM0Rko6oCpLy/ZLbw0czTzC5MU8kvoCO9E9ux7arUtEPBkRz+bRnwKHdlLZOtoW8+iZiHiq0tUQEdcBW0vapYuLVZOkrUkH2ssj4tc1svSY/dJeXXrSfqmIiLXAjcDYqkktO4b11sDRzCNMZgIT8vAHgD9GvsrUzbRbl6r+5veR+nZ7opnASfkunlHAuohY3tWF2hySXlPpb5Y0kvRdfLJrS/VyuYwXAQsj4nt1svWI/dJMXXrQfmmT1D8PbwccCdxfla1lx7Ae8ciRjhZ1HmEi6SxgXkTMJH3ALpW0iBSlx3ddietrsi6fkvQ+YCOpLhO7rMANSLqCdFfLLpKWAlNIF/2IiB8xfMjVAAADc0lEQVQD15Hu4FkEbABO7pqStq+JunwA+JikjcAzwPhuemJyOPAR4N7cnw7wJWAP6HH7pZm69JT9MgiYrvQnd1sBMyLid511DPMjR8zMrJTe2lVlZmabyYHDzMxKceAwM7NSHDjMzKwUBw4zMyvFgcN6BUk7F554+rikv+fhtZIW1JnnLElHNljmREk/aF2p2yfp2LJPcJX0B0kDWlUm2/I5cFivkH89f1B+RMOPgXPz8EHAC3Xm+UpE/KEzy7kZPg/8qOQ8l5KenGq2WRw4zKCPpJ/m/zW4If8SF0kXS/pAHn6TpP/JD7+7XdIOxQVIerekWyTtkue7IOd/pLKMnO9zkubmB06emdNeLen3edn3SfpgTv+GpAU573eqCy1pX+DZiHiiUN4Llf5z4hFJb1f6X5CFki4uzDoTOLFjN6H1Jr3yl+NmVYYDJ0bERyXNAP4VuKwyMT/K5SrggxExV9KOpF8VV6YfB3wWOCYi1uQnVgwC/gV4HelAfbWkMXldI0kPoJsp6W1AG7AsIt6dl7eTpIHAccDrIiIqj5eocjjpoXxFA4B3kh4t89uc5xRgrqSDIuLuXMZtJO0cEd3ucRrW/bnFYQaPRkTlERR3AMOqpu8HLI+IufDig/A25mnvAL4AvDsi1hTm+U1EvBARC4DdctqY/LqLdMB/HSmQ3AscKembkt4aEeuAp4B/Aj+T9H7SozyqDQJWVaX9Nj8i415gRUTcGxEvAPOr6rUS2L3RRjGrx4HDDJ4tDD/Py1viov7jqB8h/Zvcvg2WqcL71yvXWiJin4i4KCIeJD2x+F7g65K+kgPTSNKTXI8Frq+x7meAbeus94WqMrxQVa9tKbSazMpw4DBr3/3A7pLeBCBph/yYaoDFwPuBS9T+n/7MAv5d6f8gkDRY0q6Sdgc2RMRlwHeAQ3KenfKjvT9DuohfbSGwT9nK5Ke/vgZ4rOy8ZuBrHGbtiojn8gXr7+cL58+QHmNdmf6ApA8Bv5T03gbLuUHS64Fb8nWQ9cCHSQf/b0t6Afhf4GOkVsy1krYltVROq7HIm4HvSlLJJ7geCtxa6G4zK8VPxzXrwSSdT7qu0fRtw3memRExp3Ulsy2Zu6rMeravAa8qOc99Dhr2SrjFYWZmpbjFYWZmpThwmJlZKQ4cZmZWigOHmZmV4sBhZmal/H+A3QM6+kXspQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0xa6552b0>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.mlab as mlab\n",
    "import matplotlib.pyplot as plt\n",
    " \n",
    "# example data\n",
    "mu = np.mean(log_data['Thickness'].values) # mean of distribution\n",
    "sigma = np.std(log_data['Thickness'].values) # standard deviation of distribution\n",
    "x = log_data['Thickness'].values\n",
    "# the histogram of the data\n",
    "plt.hist(x, bins=[0,0.25,0.5,0.75,1.0,1.25,1.5,1.75,2,2.25,2.5,2.75,3.0], alpha=0.5)\n",
    "plt.xlabel('Thickness (m)')\n",
    "plt.ylabel('Count')\n",
    "mystring=\"Histogram with a mean of \"+ str(mu)\n",
    "plt.title(mystring)\n",
    " \n",
    "# Tweak spacing to prevent clipping of ylabel\n",
    "#plt.subplots_adjust(left=0.15)\n",
    "plt.show()\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import numpy as np\n",
    "# cmap = plt.get_cmap('viridis')\n",
    "# colors = cmap(np.linspace(0, 1, len(lithCounts.index)))\n",
    "# colors\n",
    "\n",
    "# for row in log_data.itertuples():\n",
    "#     boreid=row[3]\n",
    "#     for ind,value in enumerate(recs):  \n",
    "#         try:\n",
    "#             value.index(boreid)\n",
    "#             print(recs)\n",
    "#         except:\n",
    "#             continue\n",
    "#     #(row[3])\n",
    "\n",
    "\n",
    "\n",
    "# for ind, value in enumerate(recs):\n",
    "#     #Get the lat lon value\n",
    "#     lon=value[18]\n",
    "#     lat=value[17]\n",
    "#     #Get the Lithology unit\n",
    "#     value[]\n",
    "    \n",
    "#     #Now plot it\n",
    "#     plt.plot(lon,lat,\"|\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Exercise\n",
    "Go to http://www.bom.gov.au/water/groundwater/explorer/map.shtml and pick another River Region. Download the dataset in \"Shapefile\" format (this will downlaod the csv also). Once you have the data, follow the same routines as above and see what you can find out about the river region. Submit your jupyter notebook for review."
   ]
  },
  {
   "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}