2016-3-6FindingOptimalTimeWindow.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Finding the Optimal Time Window for Detecting RR Interval"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy\n",
    "import matplotlib.pyplot as plt\n",
    "# %matplotlib inline\n",
    "fps = 30.323427803"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function obtains a list of color component from txt file.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def readFromFile(fileName):\n",
    "    l =numpy.loadtxt(fileName)\n",
    "    return l[0], l[1], l[2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function runs through the color list, returns the peak that are local maxima for an interval of 2 times the frameThreshold. It returns the list of peaks. For exmple, if frameThreshold is 10, it returns only the maxima that are greater than the values of 10 frames before and those of 10 frames after."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def detectPeak(colorL, frameThreshold):\n",
    "    peakL = list()\n",
    "    colorlistLen = len(colorL)\n",
    "    for i in range(1, colorlistLen-1):\n",
    "        biggerThanPrevious = colorL[i] > colorL[i-1]\n",
    "        biggerThanAfter = colorL[i] > colorL[i+1]\n",
    "        #if it is a local maximum, then check if it is a RR peak by comparing its maxima interval to frameThreshold\n",
    "        if biggerThanPrevious and biggerThanAfter:\n",
    "            for k in range(2, frameThreshold+1):\n",
    "                if i-k >= 0:\n",
    "                    if not colorL[i] > colorL[i-k]:\n",
    "                        #this is not a rr peak\n",
    "                        break\n",
    "                if i+k < colorlistLen:\n",
    "                    if not colorL[i] > colorL[i+k]:\n",
    "                        #this is not a rr peak\n",
    "                        break\n",
    "            else:\n",
    "                #this IS a rr peak, store the array index\n",
    "                peakL.append(i)\n",
    "    return peakL"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function constructs rr Interval list from the peak list generated from detectPeak function"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def constructRRIntervalList(peakL):\n",
    "    return [ (peakL[i]-peakL[i-1])/fps for i in range(1, len(peakL))]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def constructTimeListFromRRList(rrL):\n",
    "    tL = list()\n",
    "    tL.append(rrL[0])\n",
    "    for i in range(1, len(rrL)):\n",
    "        tL.append(tL[i-1]+rrL[i])\n",
    "    return tL"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function graph the data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def plot(inVarList, deVarList, inVarName, deVarName, yBegin=0, yEnd=1.5):\n",
    "    plt.plot(inVarList, deVarList)\n",
    "    plt.xlabel(inVarName)\n",
    "    plt.ylabel(deVarName)\n",
    "    plt.title(inVarName+' vs '+deVarName)\n",
    "    plt.ylim(yBegin, yEnd)\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function defines the criteria for being a true RR interval."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def isTrueRRInterval(rrInterval):\n",
    "    if 0.6 <= rrInterval and rrInterval <= 1.2:\n",
    "        return True\n",
    "    else:\n",
    "        return False"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function detect the amount of errors in a list of RR interval, as judged by function isTrueRRInterval"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def getAmountOfError(rrL):\n",
    "    errorNum = 0\n",
    "    for rrInterval in rrL:\n",
    "        if (not isTrueRRInterval(rrInterval)):\n",
    "            errorNum += 1\n",
    "    return errorNum"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This function analyze the rgb txt file of a PPG video and return important information."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def analyzeVideoRGB(fileName):\n",
    "    print(\"Analyzing \"+fileName)\n",
    "    rL, gL, bL = readFromFile(fileName)\n",
    "    for frameThreshold in range(1, 30):\n",
    "        rRRL = constructRRIntervalList(detectPeak(rL, frameThreshold))\n",
    "        gRRL = constructRRIntervalList(detectPeak(gL, frameThreshold))\n",
    "        bRRL = constructRRIntervalList(detectPeak(bL, frameThreshold))\n",
    "        rError = getAmountOfError(rRRL)\n",
    "        gError = getAmountOfError(gRRL)\n",
    "        bError = getAmountOfError(bRRL)\n",
    "        print(\"For Frame Threshold of {:2d}: Red has {:2d} errors | Green has {:2d} errors | Blue has {:2d} errors\"\n",
    "              .format(frameThreshold*2, rError, gError, bError))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "videoList = [\"/home/mong/Desktop/map/map/videos/VID_20150123_172926_F.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173259_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173329_F.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173358_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173427_F.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173457_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173556_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173629_F.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173700_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173757_Q.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_173028_F.mp4.txt\",\n",
    "    \"/home/mong/Desktop/map/map/videos/VID_20150123_172857_Q.mp4.txt\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_172926_F.mp4.txt\n",
      "For Frame Threshold of  2: Red has 52 errors | Green has 11 errors | Blue has 46 errors\n",
      "For Frame Threshold of  4: Red has 28 errors | Green has  9 errors | Blue has 23 errors\n",
      "For Frame Threshold of  6: Red has  8 errors | Green has  4 errors | Blue has 16 errors\n",
      "For Frame Threshold of  8: Red has  6 errors | Green has  4 errors | Blue has  8 errors\n",
      "For Frame Threshold of 10: Red has  3 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 12: Red has  3 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 14: Red has  2 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 16: Red has  2 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 18: Red has  2 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 20: Red has  2 errors | Green has  2 errors | Blue has  3 errors\n",
      "For Frame Threshold of 22: Red has  2 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 24: Red has  2 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 26: Red has  3 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 28: Red has  3 errors | Green has  4 errors | Blue has  2 errors\n",
      "For Frame Threshold of 30: Red has  3 errors | Green has  4 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  4 errors | Green has  4 errors | Blue has  3 errors\n",
      "For Frame Threshold of 34: Red has  5 errors | Green has  7 errors | Blue has  6 errors\n",
      "For Frame Threshold of 36: Red has  9 errors | Green has  7 errors | Blue has  6 errors\n",
      "For Frame Threshold of 38: Red has 10 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 40: Red has  8 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 42: Red has  8 errors | Green has  5 errors | Blue has  8 errors\n",
      "For Frame Threshold of 44: Red has  6 errors | Green has  3 errors | Blue has  8 errors\n",
      "For Frame Threshold of 46: Red has  6 errors | Green has  3 errors | Blue has  8 errors\n",
      "For Frame Threshold of 48: Red has  6 errors | Green has  3 errors | Blue has  8 errors\n",
      "For Frame Threshold of 50: Red has  6 errors | Green has  3 errors | Blue has  8 errors\n",
      "For Frame Threshold of 52: Red has  6 errors | Green has  3 errors | Blue has  8 errors\n",
      "For Frame Threshold of 54: Red has  6 errors | Green has  3 errors | Blue has  7 errors\n",
      "For Frame Threshold of 56: Red has  6 errors | Green has  3 errors | Blue has  7 errors\n",
      "For Frame Threshold of 58: Red has  6 errors | Green has  3 errors | Blue has  7 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173259_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 40 errors | Green has 39 errors | Blue has 165 errors\n",
      "For Frame Threshold of  4: Red has 32 errors | Green has 32 errors | Blue has 85 errors\n",
      "For Frame Threshold of  6: Red has 26 errors | Green has 24 errors | Blue has 60 errors\n",
      "For Frame Threshold of  8: Red has 13 errors | Green has 13 errors | Blue has 35 errors\n",
      "For Frame Threshold of 10: Red has  3 errors | Green has  3 errors | Blue has 24 errors\n",
      "For Frame Threshold of 12: Red has  0 errors | Green has  0 errors | Blue has 20 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  0 errors | Blue has 17 errors\n",
      "For Frame Threshold of 16: Red has  0 errors | Green has  0 errors | Blue has 14 errors\n",
      "For Frame Threshold of 18: Red has  0 errors | Green has  0 errors | Blue has 13 errors\n",
      "For Frame Threshold of 20: Red has  0 errors | Green has  0 errors | Blue has 13 errors\n",
      "For Frame Threshold of 22: Red has  0 errors | Green has  0 errors | Blue has  9 errors\n",
      "For Frame Threshold of 24: Red has  0 errors | Green has  0 errors | Blue has  6 errors\n",
      "For Frame Threshold of 26: Red has  0 errors | Green has  0 errors | Blue has  4 errors\n",
      "For Frame Threshold of 28: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 30: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 40: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 42: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 44: Red has  1 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 46: Red has  2 errors | Green has  2 errors | Blue has  3 errors\n",
      "For Frame Threshold of 48: Red has  2 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 50: Red has  3 errors | Green has  3 errors | Blue has  4 errors\n",
      "For Frame Threshold of 52: Red has  5 errors | Green has  5 errors | Blue has  5 errors\n",
      "For Frame Threshold of 54: Red has  4 errors | Green has  5 errors | Blue has  5 errors\n",
      "For Frame Threshold of 56: Red has  6 errors | Green has  8 errors | Blue has  5 errors\n",
      "For Frame Threshold of 58: Red has  6 errors | Green has  9 errors | Blue has  7 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173329_F.mp4.txt\n",
      "For Frame Threshold of  2: Red has 101 errors | Green has 24 errors | Blue has 220 errors\n",
      "For Frame Threshold of  4: Red has 42 errors | Green has 14 errors | Blue has 129 errors\n",
      "For Frame Threshold of  6: Red has 28 errors | Green has  8 errors | Blue has 76 errors\n",
      "For Frame Threshold of  8: Red has 19 errors | Green has  8 errors | Blue has 44 errors\n",
      "For Frame Threshold of 10: Red has 16 errors | Green has  7 errors | Blue has 35 errors\n",
      "For Frame Threshold of 12: Red has 12 errors | Green has  5 errors | Blue has 22 errors\n",
      "For Frame Threshold of 14: Red has 13 errors | Green has  5 errors | Blue has 15 errors\n",
      "For Frame Threshold of 16: Red has 13 errors | Green has  2 errors | Blue has 13 errors\n",
      "For Frame Threshold of 18: Red has 11 errors | Green has  2 errors | Blue has 10 errors\n",
      "For Frame Threshold of 20: Red has  9 errors | Green has  2 errors | Blue has 10 errors\n",
      "For Frame Threshold of 22: Red has  9 errors | Green has  2 errors | Blue has  9 errors\n",
      "For Frame Threshold of 24: Red has  9 errors | Green has  2 errors | Blue has  8 errors\n",
      "For Frame Threshold of 26: Red has  8 errors | Green has  2 errors | Blue has  5 errors\n",
      "For Frame Threshold of 28: Red has  8 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 30: Red has  8 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  8 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  8 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  8 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 40: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 42: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 44: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 46: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 48: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 50: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 52: Red has  8 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 54: Red has  8 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 58: Red has  7 errors | Green has  0 errors | Blue has  1 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173358_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 36 errors | Green has 37 errors | Blue has 123 errors\n",
      "For Frame Threshold of  4: Red has 30 errors | Green has 32 errors | Blue has 91 errors\n",
      "For Frame Threshold of  6: Red has 22 errors | Green has 23 errors | Blue has 64 errors\n",
      "For Frame Threshold of  8: Red has 11 errors | Green has 11 errors | Blue has 51 errors\n",
      "For Frame Threshold of 10: Red has  6 errors | Green has  4 errors | Blue has 28 errors\n",
      "For Frame Threshold of 12: Red has  2 errors | Green has  1 errors | Blue has 13 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  1 errors | Blue has  7 errors\n",
      "For Frame Threshold of 16: Red has  1 errors | Green has  1 errors | Blue has  5 errors\n",
      "For Frame Threshold of 18: Red has  1 errors | Green has  1 errors | Blue has  2 errors\n",
      "For Frame Threshold of 20: Red has  1 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 22: Red has  1 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 24: Red has  1 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 26: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 28: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 30: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 40: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 42: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 44: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 46: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 48: Red has  1 errors | Green has  2 errors | Blue has  1 errors\n",
      "For Frame Threshold of 50: Red has  1 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 52: Red has  1 errors | Green has  2 errors | Blue has  4 errors\n",
      "For Frame Threshold of 54: Red has  5 errors | Green has  5 errors | Blue has  5 errors\n",
      "For Frame Threshold of 56: Red has  9 errors | Green has 11 errors | Blue has  9 errors\n",
      "For Frame Threshold of 58: Red has  9 errors | Green has 10 errors | Blue has  7 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173427_F.mp4.txt\n",
      "For Frame Threshold of  2: Red has 32 errors | Green has 36 errors | Blue has 126 errors\n",
      "For Frame Threshold of  4: Red has 15 errors | Green has 22 errors | Blue has 60 errors\n",
      "For Frame Threshold of  6: Red has 10 errors | Green has 18 errors | Blue has 40 errors\n",
      "For Frame Threshold of  8: Red has  3 errors | Green has 12 errors | Blue has 20 errors\n",
      "For Frame Threshold of 10: Red has  3 errors | Green has  6 errors | Blue has 15 errors\n",
      "For Frame Threshold of 12: Red has  2 errors | Green has  6 errors | Blue has 15 errors\n",
      "For Frame Threshold of 14: Red has  2 errors | Green has  2 errors | Blue has 10 errors\n",
      "For Frame Threshold of 16: Red has  2 errors | Green has  2 errors | Blue has  6 errors\n",
      "For Frame Threshold of 18: Red has  2 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 20: Red has  0 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 22: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 24: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 26: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 28: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 30: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 32: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  0 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 40: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 42: Red has  1 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 44: Red has  2 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 46: Red has  3 errors | Green has  5 errors | Blue has  4 errors\n",
      "For Frame Threshold of 48: Red has  4 errors | Green has  5 errors | Blue has  4 errors\n",
      "For Frame Threshold of 50: Red has  4 errors | Green has  6 errors | Blue has  6 errors\n",
      "For Frame Threshold of 52: Red has  7 errors | Green has  7 errors | Blue has  7 errors\n",
      "For Frame Threshold of 54: Red has  8 errors | Green has  7 errors | Blue has  8 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  8 errors | Blue has  8 errors\n",
      "For Frame Threshold of 58: Red has  6 errors | Green has  8 errors | Blue has  8 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173457_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 37 errors | Green has 37 errors | Blue has 158 errors\n",
      "For Frame Threshold of  4: Red has 27 errors | Green has 29 errors | Blue has 107 errors\n",
      "For Frame Threshold of  6: Red has 19 errors | Green has 17 errors | Blue has 78 errors\n",
      "For Frame Threshold of  8: Red has 11 errors | Green has  7 errors | Blue has 56 errors\n",
      "For Frame Threshold of 10: Red has  0 errors | Green has  0 errors | Blue has 31 errors\n",
      "For Frame Threshold of 12: Red has  0 errors | Green has  0 errors | Blue has 17 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  0 errors | Blue has  8 errors\n",
      "For Frame Threshold of 16: Red has  0 errors | Green has  0 errors | Blue has  6 errors\n",
      "For Frame Threshold of 18: Red has  1 errors | Green has  0 errors | Blue has  4 errors\n",
      "For Frame Threshold of 20: Red has  1 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 22: Red has  1 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 24: Red has  1 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 26: Red has  1 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 28: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 30: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 40: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 42: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 44: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 46: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 48: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 50: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 52: Red has  1 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 54: Red has  3 errors | Green has  3 errors | Blue has  4 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  7 errors | Blue has  8 errors\n",
      "For Frame Threshold of 58: Red has 10 errors | Green has  9 errors | Blue has 10 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173556_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 33 errors | Green has 35 errors | Blue has 131 errors\n",
      "For Frame Threshold of  4: Red has 23 errors | Green has 24 errors | Blue has 100 errors\n",
      "For Frame Threshold of  6: Red has 19 errors | Green has 19 errors | Blue has 76 errors\n",
      "For Frame Threshold of  8: Red has 10 errors | Green has 10 errors | Blue has 53 errors\n",
      "For Frame Threshold of 10: Red has  4 errors | Green has  0 errors | Blue has 32 errors\n",
      "For Frame Threshold of 12: Red has  0 errors | Green has  0 errors | Blue has 11 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 16: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 18: Red has  1 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 20: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 22: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 24: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 26: Red has  1 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 28: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 30: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 40: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 42: Red has  1 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 44: Red has  1 errors | Green has  1 errors | Blue has  2 errors\n",
      "For Frame Threshold of 46: Red has  1 errors | Green has  1 errors | Blue has  2 errors\n",
      "For Frame Threshold of 48: Red has  1 errors | Green has  2 errors | Blue has  3 errors\n",
      "For Frame Threshold of 50: Red has  4 errors | Green has  4 errors | Blue has  4 errors\n",
      "For Frame Threshold of 52: Red has  5 errors | Green has  5 errors | Blue has  7 errors\n",
      "For Frame Threshold of 54: Red has  7 errors | Green has  7 errors | Blue has  7 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  7 errors | Blue has  9 errors\n",
      "For Frame Threshold of 58: Red has  7 errors | Green has  7 errors | Blue has  8 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173629_F.mp4.txt\n",
      "For Frame Threshold of  2: Red has 191 errors | Green has 229 errors | Blue has 237 errors\n",
      "For Frame Threshold of  4: Red has 91 errors | Green has 134 errors | Blue has 140 errors\n",
      "For Frame Threshold of  6: Red has 50 errors | Green has 98 errors | Blue has 97 errors\n",
      "For Frame Threshold of  8: Red has 35 errors | Green has 76 errors | Blue has 64 errors\n",
      "For Frame Threshold of 10: Red has 25 errors | Green has 61 errors | Blue has 43 errors\n",
      "For Frame Threshold of 12: Red has 17 errors | Green has 43 errors | Blue has 33 errors\n",
      "For Frame Threshold of 14: Red has 10 errors | Green has 33 errors | Blue has 28 errors\n",
      "For Frame Threshold of 16: Red has 10 errors | Green has 27 errors | Blue has 21 errors\n",
      "For Frame Threshold of 18: Red has 11 errors | Green has 21 errors | Blue has 18 errors\n",
      "For Frame Threshold of 20: Red has 11 errors | Green has 15 errors | Blue has 12 errors\n",
      "For Frame Threshold of 22: Red has 10 errors | Green has 13 errors | Blue has  8 errors\n",
      "For Frame Threshold of 24: Red has 10 errors | Green has 11 errors | Blue has  4 errors\n",
      "For Frame Threshold of 26: Red has  9 errors | Green has  9 errors | Blue has  2 errors\n",
      "For Frame Threshold of 28: Red has  9 errors | Green has  3 errors | Blue has  2 errors\n",
      "For Frame Threshold of 30: Red has  9 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 32: Red has  9 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 34: Red has  4 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 36: Red has  4 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 38: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 40: Red has  3 errors | Green has  1 errors | Blue has  5 errors\n",
      "For Frame Threshold of 42: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 44: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 46: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 48: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 50: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 52: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 54: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 56: Red has  3 errors | Green has  1 errors | Blue has  4 errors\n",
      "For Frame Threshold of 58: Red has  3 errors | Green has  1 errors | Blue has  3 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173700_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 37 errors | Green has 35 errors | Blue has 135 errors\n",
      "For Frame Threshold of  4: Red has 28 errors | Green has 31 errors | Blue has 92 errors\n",
      "For Frame Threshold of  6: Red has 17 errors | Green has 21 errors | Blue has 71 errors\n",
      "For Frame Threshold of  8: Red has  9 errors | Green has  8 errors | Blue has 33 errors\n",
      "For Frame Threshold of 10: Red has  2 errors | Green has  2 errors | Blue has 13 errors\n",
      "For Frame Threshold of 12: Red has  0 errors | Green has  0 errors | Blue has  9 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  0 errors | Blue has  6 errors\n",
      "For Frame Threshold of 16: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 18: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 20: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 22: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 24: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 26: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 28: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 30: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  0 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 36: Red has  0 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 38: Red has  0 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 40: Red has  0 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 42: Red has  0 errors | Green has  0 errors | Blue has  0 errors\n",
      "For Frame Threshold of 44: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 46: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 48: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 50: Red has  0 errors | Green has  0 errors | Blue has  3 errors\n",
      "For Frame Threshold of 52: Red has  2 errors | Green has  3 errors | Blue has  6 errors\n",
      "For Frame Threshold of 54: Red has  6 errors | Green has  6 errors | Blue has  6 errors\n",
      "For Frame Threshold of 56: Red has  8 errors | Green has 10 errors | Blue has  7 errors\n",
      "For Frame Threshold of 58: Red has  9 errors | Green has 10 errors | Blue has  8 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173757_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 36 errors | Green has 41 errors | Blue has 103 errors\n",
      "For Frame Threshold of  4: Red has 31 errors | Green has 31 errors | Blue has 76 errors\n",
      "For Frame Threshold of  6: Red has 14 errors | Green has 16 errors | Blue has 70 errors\n",
      "For Frame Threshold of  8: Red has  3 errors | Green has  5 errors | Blue has 65 errors\n",
      "For Frame Threshold of 10: Red has  0 errors | Green has  3 errors | Blue has 57 errors\n",
      "For Frame Threshold of 12: Red has  0 errors | Green has  1 errors | Blue has 39 errors\n",
      "For Frame Threshold of 14: Red has  0 errors | Green has  1 errors | Blue has 27 errors\n",
      "For Frame Threshold of 16: Red has  0 errors | Green has  0 errors | Blue has 20 errors\n",
      "For Frame Threshold of 18: Red has  0 errors | Green has  0 errors | Blue has 12 errors\n",
      "For Frame Threshold of 20: Red has  0 errors | Green has  0 errors | Blue has  9 errors\n",
      "For Frame Threshold of 22: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 24: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 26: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 28: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 30: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 32: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 40: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 42: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 44: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 46: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 48: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 50: Red has  0 errors | Green has  1 errors | Blue has  2 errors\n",
      "For Frame Threshold of 52: Red has  1 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 54: Red has  4 errors | Green has  5 errors | Blue has  5 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  8 errors | Blue has  8 errors\n",
      "For Frame Threshold of 58: Red has  9 errors | Green has  9 errors | Blue has  8 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_173028_F.mp4.txt\n",
      "For Frame Threshold of  2: Red has 40 errors | Green has 25 errors | Blue has 48 errors\n",
      "For Frame Threshold of  4: Red has 17 errors | Green has  6 errors | Blue has 35 errors\n",
      "For Frame Threshold of  6: Red has  7 errors | Green has  6 errors | Blue has 22 errors\n",
      "For Frame Threshold of  8: Red has  5 errors | Green has  6 errors | Blue has 12 errors\n",
      "For Frame Threshold of 10: Red has  5 errors | Green has  6 errors | Blue has 12 errors\n",
      "For Frame Threshold of 12: Red has  5 errors | Green has  6 errors | Blue has 12 errors\n",
      "For Frame Threshold of 14: Red has  5 errors | Green has  6 errors | Blue has 12 errors\n",
      "For Frame Threshold of 16: Red has  6 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 18: Red has  6 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 20: Red has  6 errors | Green has  7 errors | Blue has 11 errors\n",
      "For Frame Threshold of 22: Red has  6 errors | Green has  7 errors | Blue has 11 errors\n",
      "For Frame Threshold of 24: Red has  6 errors | Green has  7 errors | Blue has 10 errors\n",
      "For Frame Threshold of 26: Red has  6 errors | Green has  7 errors | Blue has 10 errors\n",
      "For Frame Threshold of 28: Red has  6 errors | Green has  7 errors | Blue has 10 errors\n",
      "For Frame Threshold of 30: Red has  4 errors | Green has  7 errors | Blue has 10 errors\n",
      "For Frame Threshold of 32: Red has  5 errors | Green has  5 errors | Blue has 11 errors\n",
      "For Frame Threshold of 34: Red has  7 errors | Green has  8 errors | Blue has 13 errors\n",
      "For Frame Threshold of 36: Red has  7 errors | Green has  7 errors | Blue has 13 errors\n",
      "For Frame Threshold of 38: Red has  6 errors | Green has  8 errors | Blue has 13 errors\n",
      "For Frame Threshold of 40: Red has  5 errors | Green has  7 errors | Blue has 12 errors\n",
      "For Frame Threshold of 42: Red has  5 errors | Green has  6 errors | Blue has 11 errors\n",
      "For Frame Threshold of 44: Red has  5 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 46: Red has  5 errors | Green has  6 errors | Blue has 10 errors\n",
      "For Frame Threshold of 48: Red has  5 errors | Green has  6 errors | Blue has  9 errors\n",
      "For Frame Threshold of 50: Red has  5 errors | Green has  6 errors | Blue has  9 errors\n",
      "For Frame Threshold of 52: Red has  5 errors | Green has  6 errors | Blue has  9 errors\n",
      "For Frame Threshold of 54: Red has  5 errors | Green has  6 errors | Blue has  9 errors\n",
      "For Frame Threshold of 56: Red has  5 errors | Green has  6 errors | Blue has  9 errors\n",
      "For Frame Threshold of 58: Red has  4 errors | Green has  6 errors | Blue has  9 errors\n",
      "\n",
      "\n",
      "Analyzing /home/mong/Desktop/map/map/videos/VID_20150123_172857_Q.mp4.txt\n",
      "For Frame Threshold of  2: Red has 32 errors | Green has 36 errors | Blue has 127 errors\n",
      "For Frame Threshold of  4: Red has 15 errors | Green has 22 errors | Blue has 60 errors\n",
      "For Frame Threshold of  6: Red has 10 errors | Green has 18 errors | Blue has 38 errors\n",
      "For Frame Threshold of  8: Red has  3 errors | Green has 12 errors | Blue has 20 errors\n",
      "For Frame Threshold of 10: Red has  3 errors | Green has  6 errors | Blue has 15 errors\n",
      "For Frame Threshold of 12: Red has  2 errors | Green has  6 errors | Blue has 15 errors\n",
      "For Frame Threshold of 14: Red has  2 errors | Green has  2 errors | Blue has 10 errors\n",
      "For Frame Threshold of 16: Red has  2 errors | Green has  2 errors | Blue has  6 errors\n",
      "For Frame Threshold of 18: Red has  2 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 20: Red has  0 errors | Green has  2 errors | Blue has  2 errors\n",
      "For Frame Threshold of 22: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 24: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 26: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 28: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 30: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 32: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 34: Red has  0 errors | Green has  0 errors | Blue has  1 errors\n",
      "For Frame Threshold of 36: Red has  0 errors | Green has  1 errors | Blue has  1 errors\n",
      "For Frame Threshold of 38: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 40: Red has  0 errors | Green has  0 errors | Blue has  2 errors\n",
      "For Frame Threshold of 42: Red has  1 errors | Green has  1 errors | Blue has  3 errors\n",
      "For Frame Threshold of 44: Red has  2 errors | Green has  3 errors | Blue has  4 errors\n",
      "For Frame Threshold of 46: Red has  3 errors | Green has  5 errors | Blue has  4 errors\n",
      "For Frame Threshold of 48: Red has  4 errors | Green has  6 errors | Blue has  5 errors\n",
      "For Frame Threshold of 50: Red has  6 errors | Green has  6 errors | Blue has  7 errors\n",
      "For Frame Threshold of 52: Red has  7 errors | Green has  7 errors | Blue has  7 errors\n",
      "For Frame Threshold of 54: Red has  8 errors | Green has  7 errors | Blue has  8 errors\n",
      "For Frame Threshold of 56: Red has  7 errors | Green has  8 errors | Blue has  8 errors\n",
      "For Frame Threshold of 58: Red has  6 errors | Green has  8 errors | Blue has  8 errors\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for video in videoList:\n",
    "    analyzeVideoRGB(video)\n",
    "    print(\"\\n\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#This is Debug"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl8VPW5/z9PAiEQFlkDSSAMixDZBQXFJQERtV5xqdYa\nVNy1ltpaW9GKQPPzVlrr7d7bW2tdYnvr1bbaRQWlqS2IuLHIDgkBEhIiiwIiS/L8/njm61nmnJlJ\nMktm8rxfr7wyc2bOme+Zk3w/51m/xMxQFEVRFC8ykj0ARVEUpe2iIqEoiqL4oiKhKIqi+KIioSiK\noviiIqEoiqL4oiKhKIqi+KIiobRbiGggEX1CRJTssShKW0VFQmk3EFEVEU0zz5l5FzN3Zy0WUhRf\nVCQUJUUhosxkj0FJf1QklHYBET0DYBCAvwRdTPcRUSERNRFRRvA9/yCiMiJaTkSHiOglIupFROVE\n9DERvU1Eg2zHHElES4hoHxFtJKKrfT77GiJ6x7XtG0T05+DjS4hofXBcu4joXp/j3EhE/yaix4no\nIwALYvT1KIovKhJKu4CZbwCwE8ClQRfTY+Yl11u/BKAUQB6AYQBWAPgNgJ4ANiE4MRNRFwBLAJQD\n6APgWgA/J6KRHh//FwCnEtFQ27YvA3gu+PgJALcxc3cAowEsC3MqkwFsA9APwCMRTltRWo2KhNLe\niBSk/i0z72DmQwBeAbCdmf/BzE0A/g/AhOD7LgVQxczPsLAGwB8BhFgTzHwUwEsQYQARDQcwAsDL\nwbccBzCKiLox88fMvDrM+GqY+RfM3MTMx6I7ZUVpOSoSiuKk3vb4qMfzrsHHhQCmENH+4M8BANcB\n6O9z3N8jKBLB9/3ZNslfBeALAKqDLq8pYca3K/pTUZTW0yHZA1CUBBLLLKZdACqYeWaU718KoC8R\njYO4pr7++aCY3wNweTAQPRfA85D4iReaiaUkFLUklPZEHYAhrm0trZH4KyTOMJuIOhBRRyKa5BOT\nADOfhLirfgCJbywFgOB+1xFRd2ZuBHAIQGMLx6QoMUdFQmlPPApgftA9ZDKI7HfmUd+lM/NhABdC\nrILa4M+jALLC7PZ7ANMBPB+McRiuB1BFRAcB3A5xRylKm4DiWUdERAUAngGQC6AJwK+Z+SdE1BPA\nHyB+3R0ArmHmj237DQKwHsACZn48bgNUFEVRwhJvS+IkgHuZeRSAswDcHTTH5wF4nZlHQNL9HnDt\n90MAf4/z2BRFUZQIxFUkmLnOpPMFzfONAAoAzALwdPBtTwO43OxDRLMAVEIsCUVRFCWJJCwmQUSD\nAYwHsBJALjPXAyIkEHcUiKgrgG8DWISWBxQVRVGUGJEQkQhO/i8AuCdoUbgDISaItwDAfzHzp2bX\nRIxPURRF8SbudRJE1AEiEM8y80vBzfVElMvM9UTUH8De4PbJAK4iou9D0gQbiegoM//CdUzNFVcU\nRWkBzNysm+9EWBJPAtjAzD+2bXsZwJzg4xshLQvAzOcx8xBmHgLgRwD+0y0QBmZO258FCxYkfQx6\nfnp+7fH80vncmFt2bx1XS4KIpkKapa0jog8gbqYHASwG8DwR3QygGsA18RyHoiiK0jLiKhLMvByA\nX8/7CyLsuyj2I1IURVGag1Zct0GKi4uTPYS4oueX2qTz+aXzubWUuFZcxwsi4lQct6IoSjIhInAb\nDFwriqIoKYqKhKIoiuKLioSiKIrii4qEoiiK4ouKhKIoiuKLioSiKIrii4qEoiiK4ouKhKIoiuKL\nioSiKIrii4qEoiiK4kvc15NQFEVpT1RXVeGp+fPRVFODjPx8zCkrQ2EgkOxhtRjt3aQoihIjqquq\n8NMZM7Bo+3bkADgCYMHQoZi7dGmbEArt3aQoipJEnpo//3OBAIAcAIu2b8dT8+cnc1itQkVCURQl\nRjTV1HwuEIYcAE21tckYTkzQmISiJJB081crTjLy83EEcAjFEQAZeXlJGlHr0ZiEoiSItu6vVlpP\nW7/GLYlJqEgoSoJYNHs27nvuuZC7zMdKS7GgvDxZw1JizLYtVfjCyPkY2LEW4y/Kw9wftR1rsSUi\noe4mRUkQ6eivVkLp2CmAI3nlOFwAXHk/UGjTh1R0N6pIKEqCSEd/tRLKzp3AoEFAv35AXZ213dMV\ntXJlm3FF+aHZTYqSIOaUlWHB0KE4Enxu/NVzysqSOSwlxuzaBQwcCPTvD9TXW9tTNT1WLQlFSRCF\ngQDmLl2Kc8bMR6/jtRh5fh6+/T+WuyEVXRFKKMaSyMlxWhKp6m5UkVCUGBNush9YGMCGE+UomQZc\ncIflr05VV0R7JJKY79oFjBgBdOwIrF5t7Zey7kZmTrkfGbaitD12VFbyN4cO5cMAM8CHAf7m0KG8\no7KSmZnr6pj79GG+9VbmX/3K2m9haenn+7Bt34WlpUk6k9RmR2UlLywt5YeLi3lhaenn33+0r4c7\nbrjry8x82WXMf/yj/Mya1bx9401w7mzefNvcHdrCj4qE0laJNNm/+y7z+PHM8+YxP/KItd/DxcWO\nfczPwyUlSTqT1CXSZNyayToaMT/9dOZVq5hXrGCePDl0bJN7lfI5KOH7r4henGJFS0RCA9eKEkMi\n+Z2rqoDCQqBvX6ChwXqPcUXYSQlXRBskUoC4NQHkaOIKu3cDBQUSuLbHJAAgf2AA6z4rx9bcZbhl\ncXlKuBJVJBQlhkSa7LdtA4YPB/r0AT76yHqPZj7FjkgTeWsCyJGu77FjwIEDkv6amysiwba6361b\ngbw8YMAA4PDhZpxUElGRUJQYMqesDA8M9J/st20Dhg0LtSRM5tPMfqW4pGsJ5owr1aB1C4k0kbfG\naosk5jU1IgKZmUCXLkDnzsC+fdb+q1cD48cDXbuqSChKu6QwEMCk7y7FeJTi2rwSPFbqnOztImG3\nJMy+dd3KUfDlZRg4LTVcEW2ROWVleHiI/0Q+p6wMX+/TMqutMBDAV5csxaTMUlyYXYIHZjqvr3E1\nGYYNE+vBYESie3fg4MHQ41dXVWHR7NlYUFKCRbNno7qqqtnnH2s0BVZRYkwTAtiGcpw7E1jwpPO1\nrVvF3dTU5LQkAHFL7N4NjB0LvP9+4sabbhQGArj62aUYP3U+RvcO9k+ypakWBgLI/tJSjP/5fEwr\nqkXe6c7XI3FKrwBqupRjeDFQcpOz7YZbJIYPl2t+1lnyfPVqYO5ceV91tfO4bTUNWkVCUWJMXZ24\nHD77zLn9yBFg/36ZRD77DNi7V4SBgu3WGhqkACsvD1i2LPHjTieMUA84DVjg0Tvx0OEAaruUY9Sd\nwNe+5nwtUh2EEYKCAnEv2XGLxKmnWpYEM/DBB2JJbNoEVFY69/ULqD82f35SG0CqSChKjKmrAwKB\nUJGorJTtGRnir+7WTdo29O8vr5t2Dj16eLsilOjZuRMYOVKCyF7U1ACjRoV+z9HczRshyM+Xx3Z2\n7RIXk2H4cOAvf5HHGzZILKKgABg8GFixwrlvW63I1piEosSYmhpgyJBQkTDxCMOQIc67SSMSp5yi\nIhGJSL77nTuBcePEcvOittZbJKJJjzXXKRpLYvRo4F//AhobgZdeAi6+WLYPGADs2ePct62mQasl\noSgxZtcu4NJL5bedrVudIhEISN3E2Wdb+4UTCe3tJERzt29E4uWXvY+xZw9QVCRuHzvR3M2b65Sf\nH1kkxo4FPv5YBGnzZuDNN2X7gAGhNRRzysrwnX+vxCPVrgWLkpwGrZaEosSYnTvFF+1lSQwfbj03\nImGwi8THHzv3NRPjfc89h0UVFbjvuefw0xkz2kT2S6KJ5m7fuJtOnpTaBTtHjwKffiqWnFuMo7mb\n373bsiTc7ia3SBDJtSwrA+67D5g6VbabQjt7DUVhIIAJC5diclYpLs4JzYxLFnEVCSIqIKJlRLSe\niNYR0deC23sS0RIi2kxErxFRj+D2M4joA9vP5fEcn6LEmuPHJbXVKyaxebNTJPzcTT16yMRin0BS\ntc10PIj2br+wEOjVKzQusWeP3Mn37BkqEtEUNe7aZcUkamqs63T8uNREmBiTITMTuPpq4Ac/kHgU\nIDGprKzQm4F9+wMYdnE51p2yDAvK20YadLwtiZMA7mXmUQDOAnA3EY0EMA/A68w8AsAyAA8E378O\nwERmngDgYgC/IiK1dpSUoaZGJqCcHLljNTQ1SWbLhAnWNj9LokMHKcKyF1u11aBma2lJXUA0d/um\nXXfPnqFxCSMSXm49U9Q4tUspzkUJvntF6N28uU7dusm1MseorRWByMyM7tz79w+NS3z4IXD++c4C\nvGQT15gEM9cBqAs+PkxEGwEUAJgF4Pzg254GUAFgHjPb7706A2iK5/gUJdaYyWnQIBGMI0dEMLZu\nlVYcvXtb7/UTCcCawLp1k+cp22Y6DC2tC5hTVoYFK1c697P57g8fFoHu3VtEwm1J1NZKmrFf7Gdg\nYQAbT5aj9wDg6cecdRDMzutkrImePUNdTZEwcYmiImvb+vXAzTdbbrJOnaI/XrxI2F06EQ0GMB7A\nSgC5zFwPfC4k/WzvO5OIPgSwBsCdzKxCoaQMO3fKBNK1KzBpElBRAZw4AbzwAjBxovO9AwdKQdUv\nfynv2bNHJh0gdAJLx95OkVxoflaGqXo+o0MpzqMSlF3pvNvftUtEmshbJMJZEoDUq3TvLj/ueMbB\ng2I9GPEuKJBr/u67wJ/+1DyRcFsSTU0iEqNHA9nZoZ+dLBKS3UREXQG8AOCeoEXBrrd8/pyZVwEY\nTUQjADxDRK8w83H3MRcuXPj54+LiYhQXF8dj6IoSQrgso+pqyYEHgJkzgb//HXjxReCVV4AnXdXX\nHTsCzz8PPPSQ+KZ79bLuHN21EsYNcuGU+cjYW4tzrsrD3B+kdnZTOBdaJCsjp1sA9d3LgVzghjLn\n3b6x5gD5Tt3uJmNJmNhPU5MVKwDEMsjPF6vBPVHbrQhAspdWrAAeeUSef+Mb0Z9/ICDJDIbqahGu\nU06Rv4NYiERFRQUqKipadYy4iwQRdYAIxLPM/FJwcz0R5TJzPRH1B7DXvR8zbyaiwwBGAwhpUmAX\nCUWJJWFFIMLktWMHMHmyHOdLX5LgdI8eEqDu1Sv0s66+WtwUzz/vnHxOOSX0DrgwEEBlRjnqCXju\ncWsiTFXCudAiVR+bmpOOHUNFwC4S/fo515kGRCRGjhSLoEsXcU917269bkRi714JRttxi8SMGcBl\nl1nP7a9FYtQoq9AOAN55x4pZZWeHJj60BPcN9KJFi5p9jES4m54EsIGZf2zb9jKAOcHHNwJ4CRCX\nFBFlBh8XAhgBYEcCxqi0MZLV6CxSqmkkF0lVlWVJBALAH/4AvPqqt0AYuneXgKV9gsnLk8nMzief\nyM+wYZLCmeqEc6FFCtQbkejVKzTIa1x+gHyP7uDwnj2yHfB2OZnYgtfdfHW18zrNnAk88wzw6KPy\nvDnuptGjxb1keOstK0U2VpZELIh3CuxUAKUApgVTWt8noosALAYwg4g2A5gOIPgV4xwAa4jofQAv\nAriLmX1qJpV0JZk1AZFEINLktWOHJRIAcM01wJQp4T+zRw+5m7UX2g0aJJOdHdMcsGvX9BAJ40L7\n4tBSTO9YgtmjrdhCNOtyDB3q7U6yC/WAAaFia9xNgLdIVFdL+qzXRO2udQHEYpw2TR67XwvHyJHA\n9u0SjwJEJEwjwLYUk4irSDDzcmbOZObxzDyBmU9n5leZeT8zX8DMI5j5QmY+GHx/OTOPDr5vEjP/\nJdJnKOlHvGsCwlkpkUQg3OTV2Ch3oc11AxlXx2mnWdu8RGLzZmDECHGReIlEW2wzHYnCQADHBpWj\nx6XLkDvVqguIFKg3guknEiZU49VfqbZWxANomUjYxdyQmyu/x46N/tyzs+VYU6eK22nzZuCMM+S1\ndmNJKEpLiGdNQCQrJdIdbLjJa/duWSeiuWmLPXrIb3sqpJ9InHqqt0ikckX29u3AmWc619cwVsZZ\n2aWY1qEED17kzGDaskW+i0giMWyYMzh89Kj8GPefl0js2CEikZUVGpPwsiQAuV4nTzoD4NHwve+J\nWFx+OXDllXJtAfkbikVMIhaoSCgJJ9IdbzwbnUWyUuaUlWFegf8drEm/nJghxVbf/7I1eW3aJC6E\n5tK1q/x2i4S799O6dcCYMWJ5uCt123JFdrjrfeyY1AqMGiU1JXb65gawhcvRoWQZLvyqZWUwi0h4\nWRJHj8pz86eSlyeuPPN97dkjqaemPXtzLInGRhGgIUO8zzPaIjo7l14q/ZyeegqYN8/a3pbcTdrg\nT0ko0RRQzSkrwz2vrsSP9/k3Omtps7tIVkphIIBx85di/B3zcfqAWpw2LXRBmp69A6jNKcf+bOD/\nHrfaMGzc6Jzoo6VvX/ltz7DJzxe3SGOjNfmsWSN3nsuXhzaWa6sV2ZGud3W1BHtNXMZOVZVM1vn5\nzuDz3r1yl9+rV2gWmAksm++MSMRkyxZx5ZiWGga3SBw7JoHwvLxQkaiulmtl7vZjyfXXO5+ru0lp\nt0Rzx1sYCKDjVUtxbtdSXNM/tNFZa1wr0Vgpx08GcCi3HEcneffPMdkv7k6eLbUk+vVz9mkCZJLI\nzbVcTh9/LKmcw4Z5N5Zrq22mIxbMBetKcnJCLYnKSglOu9tqG1cTEJoqatbssDNihLjqmpqA4mJL\nlAGpyravELhnj3zvmZkiRPaJ+v33ZcGgRKDuJqXdEu0d75GjAQy9qByHxodO1K1xrcwpK8ODg8JX\nLtfUSLW0O7/esHu33N26K2Zbakn4MXKk1cp67VpxNWVmeotENBXZ4dw+8Qp6R7repqbBTySGDAkV\nCRO0BkLvuDdscCYAAPL8738X1w7g/JzBg53LiO7fb8UrOnVyxiTeey+0aj5eqLtJabdE24Nozx6g\npESKzNy0xrVSGAjg7EeXYvx18zF5UC2GnxvqTtq9WyaDp57yPkZNjUzUzJYlceyYtTRlrJg0Cfjr\nX2WhmtWrZX0EQD7bHa8wsZLTh81HP67F9OvyMPf/RVcECCBuaytHut6mOC2cJeFeztVuSbhF4sMP\ngXPOcR7n/POBBQuAW2+V9ij2lOTBg539sw4ccIqE25KYO7cZJ98KPGs0krSeiIqEklAiNWcz7Nol\nGS8//nHoMVrb7O74yQC2UznOnwEseCL09d27gauuAhYvdq5BbX+9oECyWYwu3X67TE7hiuaayze+\nIa6oKVPkLtZMbl6WBAB06hzAgT7laDgJ/PknzmaC4SqYEXwcj7WVI13vnTtl0aWuXb1FYto0cQ/Z\nLYk1a+T7BrxF4s47ncc57zzZxys91d1kcf9+6fcESLzCWJP19VIRbeoY4k23blI4aWhpM8RYoO4m\nJaGY1MYbx5bifCrBtSNCWzEzy+QxaZL805486TzGnLIy3NnV37USyXVSXS1uIT930o4d4sfu0AE4\ndCj0deNuKiyUY9XWSp77iy+27Dvxo29f6fd0113A738PXHSRbDeB3MZG5/u3b/cvMAtnfcUz6G2u\n96xBkg320BdKQ1aQi9bdxAz813/JcqAlJfIeu0g0NorLz+1uIvKvXygokJiEOYbdkrjoIuC550Rg\nHnxQLBG78MYTt4stmdlrakkoCacwEEDfs8tRmwN0DzibswEywWVlyZ1c796SzWI3EgoDAXw4ZCnO\n3zEfEwbUomCS5TKK5o5r2za5I1y3LnRsx4+LCAQC8tn79jmzjgBxN112mYzxhRek7cYFF8gaELHm\nppvEjdXUZFURZ2XJ2OrqrK6x5ryGDpX37tvnzOePZH3Fsw15YSAAOrUc/94JLPi6dzO+rCwRgePH\nrcem3qFDBznXt94CfvQj4I035G8DcIrE22/L+buvVzhMjGfnTvm+GhosITjzTBHoW26R73bLlph8\nHVGRlydiaEhm9ppaEkpS2L1bmpnZM0sMptUz4L0wCzNQVR3AhGvKEbjBGdiO5o5r40bxU3tZEjt2\nyKRhJmL3HbkZe36+WBtr1wJ//jNwySXN/gqi5ic/AX72M+e2QYOcAVfAaUm4+xnNKSvDV3p4W19z\nyspwX//4tiHfv18mYfu1bGx0NuMzrdNXrpTf3bqJGyo7W9xuTzwhBWemgSLgFImXX7aC081h2DDJ\nfgIk9jNmjPXavfeKq+/MM2UMicLduyuZ2WsqEkpSqKmRIO/ekP6/zgZtbrMbkMm9Y0eZdJrjVgFE\nYDZtklTIvXtDU0+3brXaLhhLwmvsBQXibiouFleTcQUlCuPqsmOqgb3ErTAQwN7JS3F251LcPMyZ\nVlwYCCDvrqUYj1LcMjw+ayvX1or70H4td+1y1h0UFUkW0llniSDYi9bOOw/47W/l+7ZjRIJZ3H1X\nXdX8sdkb7a1aJYJg6NIFeO01+exE4v67T+Z6IupuUpLCwYNyJ+5nSdhFwl6LAMgkP2KEBBjNHaAh\nkltlzx6ZWAYOlBYKhw9bC8gAzrYLXiJx9Kjs06ePPC8vlwCze13jeOMnEsOG+Ytbw74ABk0vR944\nYMH/c772yaEAtqEcd90pd8+xxKz9PH688+7Y3Qdp8mTg61+Xx+vXO0Xi/vslNjV9uvPYRiQ2bZLP\nOf305o9vwgTghz8UcTBddu1ceGHzj9la3JbE5wstjZiPPidrMfXKPMx9LDHZTWpJKEnh0CHxNzc0\nhN7N210Q+fmh6Z6bN0sNgdeqY5HuuOwFb15rDdgtiQEDQiuba2rkH9hkPGVlJS7jxY5bJJidTe+8\nRGLnTpkQvV4z/n93r6JYYNZ+HjjQOfHZv2tA4i8TJ4pYrF3rFImxY6XtumlhYjAi8Y9/iIC4M9Gi\nYeZMcTNNmybpxi05Rqzp3VsEy5651a2HVPrXDFmGWxaHFnnGCxUJJSkcOiR34126hPbOsVsSY8ZI\ndokdM9F7iYTJprlkgGTTfHuG03ViL3jLzQ11d9ktiZEj5f12jKsp2RQWOhsA7t8vk1uvXt7upqNH\nZdIpKvIWCZPR5VXAFTFbLMLr5nqa9hgGd7O8gQNlGdAJE6waiUgYkaioCHVFRUvfviKyy5cDv/hF\ny44RazIyRFjtomqE3KyolyjU3aQknBMnxHWQnS3/oA0NVm464AxcT5gAfOtbzv03bZK7vl69QkUC\nEKE4XliOmg7AV25yZtPYLYnc3PCWRFER8PTTztdN0DrZBAIykRq2bJFxE3m7m8xE3aePvyVxzjne\nBVxhey9FkU1m6kqKimScph/Vxo3AueeGjsW4//wa6dnp1EkEcNkycRm1hrPPbt3+scZ0AjYGg2ku\nuH9/YkVCLQkl4Rw+LG4DInH5uO/m7YHroUNFCOwTWzh3k2H/fqlQdrfbtlsS7s+2p78C8r6NG53u\nMDPhJZshQ2TSMLUSGzZY5+W3WtugQd4C8skn0ieooCDU3RQpWyyabDLzneXkiDAbcfvgA+8YgnEp\nRSMSHYK3uRdd1LylQ1MBt0vRWBLduzsL7eKNioSScA4dsu4WjSVhaGyU4LK5W8/IEEEwboqjR+X1\nQEBEwitFFZDt48eHFwm3JbFjh3xuVpY1tsxM53t2724bk1GXLiJyxrNjLyLzShu2i4T7OzNtLrxa\nQUTKFosmf9/uPjztNBE005PK67s0QhVNdieR1Ee4Lb50YMgQp3vONC9MtLtJRUJJOHaRyMtzTuR1\ndTKR2Rfu6d7daiO9bZtVYNWtm4iGWf7RwCwWxtixzqB3fb1U9ZqJyW1JeC0oY6wJg7vVdDKZOlUC\ntoBT/AIBETy7BWQW0vGyMsy+XiIRKT8/mvx9u/U1ZYr045o+HbjxRu8gsVm4J9oFfM48s2VrObR1\nLrgA+OMfrY4DyYpJqEgoCeejj6yq1rFjJZPFYNYQsJOTY63EZo8pZGR4Lxpz6JBUPw8d6hSgZcsk\nuGkmn379nFaMO9sGkMlzwwbruf2uONnceivw0ENSRGd3N3XvLhO+/dxMi4uuXUVU7W2od+yQycfd\nGhuInC0WTf6+/TubMwf43e9kXehHH4UnCxa0nTbZyeTccyWG9Oyz8lzdTUq7oabGcidNmCC+aYOJ\nN9ix9/XZsMH5uldcwrR7Lix03lG/8YYzz969DKi9BbXh4otlMlu1StIk25JITJ8uk+33vy9Cac8G\nmjhRsnUMRiRMYNu+VOi+fZb15o5JmGyx6b0kW+wB1zKi5vXrikpxHpXghjGhhXg7d1qWxMCBYrn8\n4Af+55WZ2fwlYNORjAzg4YclIG+q0wcP1uymlCZZrXxTDbtIjBljFUJlZckEEk4k3nsPuOEG6zUv\n94kRid695bfx17/xhrNQzL1gzbZtIgp2Lr9c7sjPPtsKEieyPUMk7r5bvq9LLnG6XGbMEFfUFVfI\ncyMSgJU1YybufftErP1WQxs0OIBNjeUIjAf+497QXluFgQByxpej4RSg71jn64cOiatwwABrW0sW\nZmqvFBfL329trfwtd+4sloR93e54o5ZEjEjlhegTjV0kunSRO34TyPRa3a1LF6dITJpkvebVtuPA\nASul9rLLpN34f/6nbLcvCuQWCS93E5G0pbb/UzZ3sft4MmKEtLB+8knn9kmTpPHgI4/IJPPxx1ZV\n+LBhcq4GuyXhJRINDSJAp57qnyhQUyPWi7s63jQdbEvfWSrRoYN8d1u2WFl3GpNIUdryQvRtDbtI\nAJKFZFxOXqu7GUuirk4mdXvMwksk7KuLfec7Mlk+9ZSY7vZAqV0kjh2TcfkZfoMHt92JbtIkydSy\nc+aZIgQPPSQCEghY4x8+3FsksrK8K65N9pNXC3KDn0h4ufCU5tG5s7hZ7SKRyJiEupuaQTh3Ultd\niL4t4haJCRPE3z9lisQI3HfzOTky2f3616FtE9w9bgCnSPTtG9rjyNC5syUSH34ok5lJf/Wiqcnp\nNmnL5OTI97JokSQG2MVv+HDpmArIOW3ZIq6o6mrvXlpbt4YXCebwIuG+nkrz6NxZelmZa9i9u8Yk\n2iSRKktbu1pae8JdtTx+PPCnP8lkc911oemMJk11/nzg6qudrw0YAKxY4dxmX10sHHZL4q23gDPO\nCP/+detk8k0VBgyQCWXtWuCLX7S229tjrFolGTS5uZIptm2b3KXa12QwlkR2dqggA+LGy8oSt1Jd\nnXM1v63V5dMKAAAgAElEQVRbQ5cTVZqHEYmbbpLn6m5qo0RyJ80pK8PdPv36FYumJplI7No5ZYr8\n0b/+uvcawnfeKe2ab7kldEEZL0vCvrpYOIxIHD4s6zX8x3+Ef//o0f7uqLZK167yfdurl8eMkbTZ\n0lJpTnjffbK9UyepgH7rLecxNm+2Ggd6WRLGMuzSRY5hn8DU3dR6jEiY7LVEu5tUJKIkkjupMBDA\n9lFLcUWgFFf1jU9P/nRg7165Y3UXy23YIC4fs/qanaIi/3bNfjGJaC2Jo0elWnfYMGDWrKhPI2Xw\nanHRubOsz/D738u533KL9dq55zpTZwGxRMaM8RcJu2Xorvb2KlBUmkfnznLjY65hot1NKhJREk1l\n6ZZtAdzw3XI0jHSulqZYuOMRrcXLkrCnd4bDWBK/+pU0EWyrgenW4NcH6cUXJWB9ww1W/yNAAt7v\nvms9/+QTsURMTMKrOaB7JUETl9ixI7XiOG0dM9V06yaJHO41zuNFGv5bxIc5ZWW4p7e/O6m+XipZ\nJ070XhZTEfbsia4nT7T07St3VfasnGjvXrOz5Z+tslJaXKQjRgDcbbc7dfK2tiZNkpRaU4C4Zo24\n2TIzw69TYURi0CCrn9SSJW1nfYZUxswnJlaXmSmxMdOqJt6oSERJYSCAE5csxQW9vd1Ja9dK11Gv\nldTaI35rDHz0kbWqWyzIyJDiNvOdHz8u1oqX28qNmUDHj3feTacTxi1hlgiNRH6+LA1r2pl88IFk\nnwHerj3AaUmMHWut//Huu871qJWW4ZUs0KOHfwfkWJOm/xrxYW9DAOffWo6KCuDFcudra9bIP0iP\nHpJz/+mnof+Y7aUiO1wm2P79gaiCys3BTF6DBolVMHCgTHTR0KcP8IUvxHY8bYk772z++RlrorBQ\nUpPNRN+zp4iwafVusFsS48YBf/ubuJmWLAG++tXYnEd75plnQhMmjFsvmpuh1qKWhI1IK2xt2SKd\nGd3tpwGxJMaOFdPaazGb9lSRHS4TzBRuxZKCAuuaNDdQWl8PzJsX2/G0JTp3lnhCczjjDKu77Ntv\nW2s+EIml4V7S1b7+x7hx8r9w++0ygY0d26rhKwCuvz40jdjPqosHKSsS9sZssSDSJG7cGFOnisvE\n3b7AiAQgKu8WiVStyI4knF6EywSzF7rFilGjJDMKaH7KZUaG+szd3HSTWANf/rL8n02caL1WUCDZ\nTIamJvm/MCJh6i3Ky6X2RYkPeXkqEhH57ndje7xIk3hlpfyDdO4sKm6/mzpxQnLJR42S5+7XgdSs\nyI4knH4CEi4TbM+e2DfIGzfO8oNrymXrKSgA/vUv4H//F7jrLmfWV0GB8297715Jyezc2dr2z39K\nHUY0achKy1BLIgreey+2x4s0iZuqU0AmIdOQDpDHhYVWDOLUU50rSgHRpdC2NcIJZzgB8Vtj4Mbv\nlmHFCkmzjCV2kdi0SUUiFgwcKALgXl/cbUl4tU4fOLBtrAOezgwYEBrQPnQodG2VWBBXkSCiAiJa\nRkTriWgdEX0tuL0nES0hos1E9BoR9Qhuv4CI3iWiNUT0DhGV+B3bvhBMtIRznUSaxO1ujMmTxVdr\nsLuaAOliahcRQFJov5UXviK7Ja6deBJOOMMJiFljYM64UhRnlODqYZIJdvRYADk5VpAzVgwdKi7A\nhgYJuEZqr6FER9++oa64/HynSNiD1kri8HI3Pf888PWvx/6z4p3ddBLAvcy8moi6AniPiJYAuAnA\n68z8fSK6H8ADAOYBaABwKTPXEdEoAK8B8CyLsi+aEg2Rei/NKSvD15esxI8abK8PHYq5wUl8yxZL\nCKZMkTYOBrdIjBgB/M//OD+/MBDAkHuW4qwF8zE4uxYTv5CHubbspkjjSwbh+lFFU4E+aHo59vUC\nuJOsMfDf/w2cf37sx5mZKcLw4IOybkSsA+OKRUEBsHSp9VxFIjl4uZvs3o5YEldLgpnrmHl18PFh\nABshk/4sAGbp8qcBXB58zxpmrgs+Xg8gm4g8kxkbG73bGvsRKeZQGAgg97aluGRAKWZ2Dq2D2LTJ\nugDGkmhqkuerV3tbEvY1hgHg8KcBDJxWjl2DQyuy22JgO9zSlNG4zxoaJDPGZB7985/xEQlAFmd5\n4glZN0KJH9G4m5T44+Vu2rxZblBjTcJiEkQ0GMB4ACsB5DJzPSBCAiAklElEXwTwPjOfcL8GSGn6\noUPRf340geOOWQGMv7ocb3VwTuKNjVJUZFIB+/WTDI433hChePttZ9FQ795SnBWSBlstx/CqWm2L\nge3CQABfXbIUkzJLMb1jCeZ/odRheUVa23jvXsmM2blTBDOeIjF3rqydPG1afI6vCAUFIgwGtSSS\nQ79+0kfr5ElrW0paEoagq+kFAPcELQrXPbbzedDV9D0At/sds7kiEc2d7+HD4nNtbHQee/Nm8c/a\n3RjTp8s6vcuXyz+Oe9GXkSNlPzvV1TJpeolEWw1sd84JYG+PcpyYsgyzvlXuWNv47tdEQC7qUoJv\nXxDa0LChQRrnZWeLQGRmxq+Las+ekrKpxJfcXCvtFVCRSBYdOkghqLkRPXlSMjDjsXZH3CuuiagD\nRCCeZeaXgpvriSiXmeuJqD+Avbb3FwD4I4DrmXmH33E/+2whFi+WP9ri4mIUFxeHHcecsjI8+K+V\n+M+d3jEHQIRhyBDJVNqxQzpfAt7B0DvukPWFf/Yz4NZbQz+vqEgybux3zjt2SJrsiRPSWC472zm+\nBStXOmMSrvElg+3bJTDcr5+1roMho0MAn+SWI2cKUHxt6NrHDQ0irpMnS7zg/PO1JiHVIRJr7dpr\ngccfV3dTMsnPF5HOzxcrIj/fmYoMABUVFaioqGjdBzFzXH8APAPgcde2xQDuDz6+H8CjwcenAFgN\n4PIIx+QpU5j//W9uFi88X8mTe5XyuSjhh79cyjsqKx2vX3cd87PPMs+axfzCC9b2u+9mfuwx57Eq\nK5kB5oEDmQ8cCP2sP/+ZefJk5qYmeX7yJHNWFvNnnzH378+8e3foPmtWV3JRh1K+oFMJf/M/QseX\nDJ56Sr6Xu+5i/vnPna+9+irztGnMd9zB/ItfOF9ramLOzmY+fJh52TLmwYOZ//nPxI1biR/19cz3\n38982mnyN33yZLJH1D75yleYf/ADefzEE8yzZ0feR6b85s3h8U6BnQqgFMA0IvqAiN4noouCIjGD\niDYDmA7g0eAudwMYCuBh2/s928E1190EACebJHC8PW8Zbl0c2sr78GE57siRssiHYdWqUEvCFITd\neKPEJ9xceqlYC3/4g5jnvXtLP6FOneSxl8vp2PEAsseU49D4ZbjqgbbRatwsXellSWzaJBZTnz6h\ny14eOSJFWDk5QEmJdAY977zEjVuJH/36AbfdJmnoeXmhKwkqiWH6dImLAuK1MDHTWBNXdxMzLwfg\n9yd0gcf7HwHwSDTH9ptow7F7t8QPzBoEbjO5vl6Oe9ZZwM9/LgvSLF8uVbzubpZmGcvSUu/PyswE\nfvITYM4c8Rd+/LHVdM5v7CbwdOSI91rDyWDrVlmM5+BBq/WFYeNGaSPd2ChuKTt794qrSUlPevSQ\n3xqPSB7FxdJC5cgRuWG7+OL4fE7KVlx7NdGLhBEJryZlzDLpFRXJ6lwrVwJlZcCMGcC99zpXUjOc\nPClWhx/nnit30NdfL8/NP1YkkejbN/SuPVls2SJFhH6WxMiR3pZEXV1oMF9JH8zfcjSLOynxoVcv\nafz3619LBwrT0j3WpGyr8Nzc5k+ku3dLIZzXamY1NRL0MRlMw4YB3/uerOB1xRXex4tkZhMBb74p\nonTggBVU8ppUAZmQL7lEmge2BUuC2eqFdPSov7upqSl0vBrQTG+MVazWYnK59Vbgyiul627//vH5\njJQViX79QvsjRcLubnJbEhs2SLWu4Uc/Al56SQSiNRk5554bum3QIGeuuWHzZrFadu5M3KpT4air\nkwysU04JtW4OHBAzNy9PLDp3BbyKRPtALYnkcsUVEjM1zUXjQcq6m4qKpNK5OdjdTW5Lwi0S550H\n/PCH8UnZHDzYWuLR0NRkVUx27ix37m4S3dvJ3q/K7W4yriYiERC1JNon6brsaypxxhnRrzzYElLW\nkjjjDJm0/va36FbeamyUO94BA7zdTevXxy87wM3gwVIzYaemRloud+8ud++ffeZ8PRm9newi0bOn\nWDfHjkl8xriaAHGfffSRuKeMqO7aFbpQipJeuNvOKOlJyloSHTtKv/vbb5fitEiYzKWsLO/AtduS\niCeBQKglsXmzFQT3Eolk9HYy6a+ApLMWFUkzQ0CC/Ga8nTvLnYzd5aSWhKKkBykrEoD4+/PyZIGU\n++6zllz0orramrTclkRjY2gn13gyYID49O0upU2brOZcXu6mZPR2cq/ydu65sm4xIK6+0aOt10aM\ncLZHV5FQlPQgpUUCkNzgN96Q+MGTT/q/z35X3LOnuE2OBJslbdokmQGJWkkrI0NcTpWV1jZ7B0cv\nSyIZvZ3cInHTTcBPfwpcdhnw1ltOd1JRkSUSR45IXciAAXEbmqIoCSLlRWLSJJmwgPCTvMn3B8Rv\nPmSINam9+64cJ5G4FyayWxJeIjGnrAz39Q/fdTWWNDVJgZy9YdjEicDs2cCKFcAXv+hcq7qoSFxQ\ngIhfIOBc9lJRlNQk5f+NJ02y3Ezhlu5zt9GdNg14/XV5/K9/xX5JzUicdpq1uh4zsG6dlcaWnR3q\nbioMBJB9zVJM6VSKa/qHrncRa6qqJCDdrZtz+2OPScLAb3/r3G4XPdMUUFGU1Cdls5sMeXnAK69I\ndpLxl3vhFokLLwQeekgmwj/9CVi4MO5DdVBUBLz6qjzes0eqt02Lg86dLUuivt6qXN5eFcCEa8px\n8CCwoDy+41u9Ghg/3vs1LwthwAArRVZFQlHSh5S3JADgooukktqvAK2pKdS/fuGFwOWXA+XlIhaJ\nLgoqKrJ6Ib33nrhyTPpodrYsZDRxosRK3n9frI333gNmzkxMNXY4kfAiJ8eK8ahIKEr6kPKWhCEn\nx18kamutGgRDp07AokWJGZsX48aJe+yOO8RSsNdomHUmLrtMmut997vA/feLhTFpUvPX924Jq1cD\nN98c/fvtIrFiBXDddfEZl6IoiSVtRKJrV2uScrN5c3yW9WsNWVkymc6cKbEJEx8BxKq54QZZkjM7\nW85tzRrggQek8jkRlsSaNSJk0WJEorpaalDOOit+Y1MUJXGklUj4WRIrV4auB9EWKCiQlNLycuda\nC127Ak8/bT3/5S9FVObMkedHjkgBoWmyFmv27xcrZ/Dg6Pfp0kXW97j/fmk4pmsMKEp60C5E4s03\nga98JbHjiZZp0+QnHHfc4Xzeu7e4nOx1CNVVVXhq/nw01dQgIz8fc8rKWpz5tHNn81NYO3cW4Vq1\nyrlgk6IoqU3aiESXLpI22tTknNzmzZOsp9/9LnljizWmoZ4RiVj3ddqzp/mFcCboPm5c6Dq7iqKk\nLhHvFUlo8w0WMjJkcvr0U+f2xYulUtisE5EOmIZ6hlj3daqtbXm1dKL6XymKkhgiikRw8ey/J2As\nrcbL5VRUBHzzm8kZT7xwt+aOdV+n+vqWL2By9dUt209RlLZJtO6m94noDGZ+J66jaSVeIvHpp/Ht\ntZ4M3CJh+jrZhaI1fZ0+/tjZciNatHW0oqQf0YYmJwN4i4i2E9FaIlpHRGvjObCW4CUSR4+mn0i4\n3U1zysowf0js+jp98omzpkRRlPZLtJbEzLiOIka0F0uiXz+rWhuQvk4zf7YU51w9H92O1GLadXmY\n+/9ant2kIqEoiiEqkWDm6ngPJBa4RYJZRCLdsm0CAVl/286RowEMml6O118H/vbfoY35msMnn7Ru\nf0VR0oe06N1kcIvEiRNS1NUhbRJ9hREjpIrcTlWViEfv3sC+fa07vloSiqIY0kok3P2b0tHVBEgl\n9N69znRfFQlFUeJBWomE25JIV5HIzJRFk7ZutbZVVcm2WIjE4cMiuIqiKGklEm5L4pNPRDjSkaIi\nWajIYFaDi4VIHD8uXXIVRVHSylufk+N0wezalfh1IhLFtGnAww8Dv/89cMUVwI4d4oZyp8e2hOPH\npaGgoihKWolEly5SLWyorgYKC5M3nnhyww2yKl9lJXDbbVJg17VrbCyJEydUJBRFEdLO3WS3JHbu\ntJYETTe6dpWV9b7xDXlu2mjEyt2kIqEoCpBmItGli3Phobq6ljeqSxVM91VT16AioShKLEk7kbBb\nEvX1Up3cHjBtmlorEswiEvFa0EhRlNQi7UUiNzd540kUZ5wB3HijPG6tSDQ2inWiK8spigKkWeDa\nrLNsaC8isWqV9bhPn9aJhAatFUWxk7aWBHP7EQk7ffpIG/GWtu3WeISiKHbSViRMUV26FtP50a2b\nuIoOHGjZ/ioSiqLYiatIEFEBES0jovXBNSi+Ftzek4iWENFmInqNiHoEt/cKvv8QEf2kuZ+XnQ0c\nOyaPjRVhsn/aC0RSeV1VZW2rrqrCotmzsaCkBItmz0a1/UUXKhKKotiJd0ziJIB7mXk1EXUF8B4R\nLQFwE4DXmfn7RHQ/gAcAzAPwGYCHAIwO/jSL7Gzgs8/kcXt0NRmMSEycKALx0xkzPl8D+wiABStX\nYu7SpZ7rTahIKIpiJ66WBDPXMfPq4OPDADYCKAAwC8DTwbc9DeDy4Hs+ZeYVAI615PM6dVKRAJyW\nxFPz538uEIAscbpo+3Y8NX++576a/qooip2ExSSIaDCA8QBWAshl5npAhARATKoZvNxN7ZFAQNp1\nAEBTTQ3cDV1zADTV1nruq9lNiqLYSYhIBF1NLwC4J2hRuHNvWpiL40TdTcKIEcD69fI4Iz8fR1yv\nHwGQYarvXKi7SVEUO3GvkyCiDhCBeJaZzaKb9USUy8z1RNQfwN7mHnfhwoWfPy4uLkZxcTE6dgRO\nnpSCsL17gdNOi8UZpB5TpwIffAAcPAjMKSvDgpUrnTGJoUMxt6zMc191NylK+lBRUYGKiopWHYO4\npQn10X4A0TMAPmLme23bFgPYz8yLg4Hrnsw8z/b6jQAmMfNcn2Oy37i7dJFW2bfeClxyCTB7dkxP\nJ2W45hrpgvv22xK8vmzUfHQ/WouzZuXh7v8q8wxaA8CKFcB998lvRVHSCyICMzcr5zOulgQRTQVQ\nCmAdEX0AcSs9CGAxgOeJ6GYA1QCuse1TBaAbgCwimgXgQmbeFO1nmrjEoUNW07v2yFNPSYuOhgag\nT78AtnA5ho0GvrwAKPTWBwBAUxOQkVbVM4qitIa4igQzLwfg1wXoAp99wkxhkTFxicOH27dIdOkC\nlJQA//63LGs6ZIisObF/f/j9GhtVJBRFsUir3k2AlQbb3i0JADj7bGD5conPTJgg30skkVBLQlEU\nO2k3HRhLQkVCAtgrVgBvvAFceCHQq1fk5n9NTdoBVlEUi7QWifbWt8nN5MnA2rXAm29KO/HevdWS\nUBSleaTddNCxoxSEqSUhcYmiInE3nXpq9JaEioSiKIa0mw6ysiTX/+hRmSTbO6+9BmzfLi6kaCwJ\nDVwrimIn7QLXWVmSAtvYCHRIu7NrPr16yY95HI27SWMSiqIY0u6eMStLVqfr2LH9tQmPhLqbFEVp\nLmk3HdhFQnGigWtFUZpL2k0HHTuKSGiTulCicTdpTEJRFDtpNx1kZUm1tVoSoXTvDnz8cfj3aExC\nURQ7aSkSn36qIuFFp06SHhwOdTcpimIn7aYDE5NQd1MomZkiAo2N/u9RkVAUxU7aTQcdO6q7yQ8i\nq47ED41JKIpiJ+2mA81uCk+nTuFFQmMSiqLYSVuRUHeTN36WxIYNsqqfupsURbGTdtOBZjeFx1Sk\n2zl+HBg1CnjlFRUJRVGcpF3jCs1uCo+XJbF2rfyuq5PfKhKKohjSbjrQYrrweMUkKivl97ZtGpNQ\nFMVJ2omEupvC4+VuqqyUlh319epuUhTFSdpNB5rdFB4vd1NVFTBpkjT/U5FQFMVO2k0Hmt0UHi+R\nqKwUkdi/X0VCURQnaTcdqCURHq+YhLEk9u/XYjpFUZyk3XTQsaOsSqcLDnnjjkk0NgK7dgETJlju\nJg1cK4piSDuRyMoCmHWi88MtErt3A337AgMGAAcOqCWhKIqTtJsOTCxCRcKbyZOB3/3Oel5VBQwZ\nIt9bdjZw8KCKhKIoFmnnlDEioe4mb26/HRg+HPjsMxGFykogEJDXevcGGhqsNbEVRVHS7p7RBKzV\nkvCmXz9g2jTgO9+R58aSAEQcGhr0u1MUxSLtRELdTZH5zW+A554DVq+WKmsjEr17Ax99pO4mRVEs\n0m46UJGITK9ewMKFwH33AevWAWPGWNsbGlQkFEWxSLvpQGMS0XHrrcA//gGsXw8UFck2E5NQkVAU\nxZB2U6laEtHRoQPw4otATY0U2AFiSezfr9+doigWaScSGriOnssvdz7v3Vt+qyWhKIoh7aYDtSRa\njkl9VZFQFMWQdtOBikTL6d9ffqtIKIpiSLvpQAPXLWfoUPmtAqsoiiFtRUInuuZTWCi/T55M7jgU\nRWk7xFUkiKiAiJYR0XoiWkdEXwtu70lES4hoMxG9RkQ9bPs8QERbiWgjEV3Y3M/UwHXLMd/dzp3J\nHYeiKG2HeFsSJwHcy8yjAJwF4G4iGglgHoDXmXkEgGUAHgAAIjoNwDUAigBcDOAXRETN+UC1JFpP\nVVWyR6AoSlshriLBzHXMvDr4+DCAjQAKAMwC8HTwbU8DMMmYlwH4X2Y+ycw7AGwFcGZzPtPEIjQm\n0TJKSoBzzkn2KBRFaSskbColosEAxgNYCSCXmesBERIi6hd8Wz6At2y71QS3NeNzxJpQS6JlLFuW\n7BEoitKWSIhIEFFXAC8AuIeZDxMRu97ifh6RhQsXfv64uLgYxcXFnz/v2FFFQlEUpaKiAhUVFa06\nBjE3e35u3gcQdQDwVwCvMPOPg9s2Aihm5noi6g/gH8xcRETzADAzLw6+71UAC5j5bdcxOdy4e/UC\nHnkEuOuuOJ2UoihKCkJEYOZmxXkTkQL7JIANRiCCvAxgTvDxjQBesm2/loiyiCgAYBiAVc39wKws\njUkoiqLEgrhOpUQ0FUApgHVE9AHErfQggMUAnieimwFUQzKawMwbiOh5ABsAnADwlbAmgw8ak1AU\nRYkNcRUJZl4OwG+6vsBnn+8B+F5rPldFQlEUJTakXcU1oIFrRVGUWJGWIqGWhKIoSmxIW5HQwLWi\nKErrSVuRUEtCURSl9ahIKIqiKL6kpUho4FpRFCU2xL3iOh5EqrhesgQ4/XSgT58EDkpRFKWN05KK\n67QUCUVRFCWUttqWQ1EURUlRVCQURVEUX1QkFEVRFF9UJBRFURRfVCQURVEUX1QkFEVRFF9UJBRF\nURRfVCQURVEUX1QkFEVRFF9UJBRFURRfVCQURVEUX1QkFEVRFF9UJBRFURRfVCQURVEUX1QkFEVR\nFF9UJBRFURRfVCQURVEUX1QkFEVRFF9UJBRFURRfVCQURVEUX1QkFEVRFF9UJBRFURRfVCQURVEU\nX1QkFEVRFF9UJBRFURRfVCQURVEUX1QkFEVRFF/iKhJE9BsiqieitbZtY4loBRGtIaKXiKhrcHtH\nInqSiNYS0QdEdH48x6YoiqJEJt6WxG8BzHRtewLAt5l5HIA/Afh2cPttAJiZxwK4EMAP4zy2NktF\nRUWyhxBX9PxSm3Q+v3Q+t5YSV5Fg5n8DOODaPDy4HQBeB3Bl8PFpAJYF92sAcJCIJsVzfG2VdP9D\n1fNLbdL5/NL53FpKMmIS64nosuDjawAMDD5eA+AyIsokogCAibbXFEVRlCSQDJG4GcDdRPQOgBwA\nx4PbnwRQA+AdAI8DWA6gMQnjUxRFUYIQM8f3A4gKAfwlGGtwvzYcwLPMPMXjteUAbmHmTR6vxXfQ\niqIoaQozU3Pe3yFeA7FBwR95QtSXmRuIKAPAQwD+O7i9M0S0PiWiGQBOeAkE0PyTVBRFUVpGXEWC\niH4HoBhAbyLaCWABgG5EdDcABvBHZn4q+PZ+AF4jokaI2+n6eI5NURRFiUzc3U2KoihK6pJyFddE\ndBERbSKiLUR0f7LHE2uIaEew0PADIlqV7PG0Fp+Cyp5EtISINhPRa0TUI5ljbCk+57aAiHYT0fvB\nn4uSOcbWQEQFRLSMiNYT0Toi+lpwe7pcP/f5zQ1uT/lrSESdiOjt4DyyjogWBLc3+9qllCURjGNs\nATAdQC0kE+pav9hFKkJElQAmMrO7viQlIaJzABwG8IxJXiCixQD2MfP3g0Lfk5nnJXOcLcHn3BYA\nOMTMjyd1cDGAiPoD6M/Mq4OdEd4DMAvATUiP6+d3fl9CGlxDIuoSjPFmQrJFvwbgKjTz2qWaJXEm\ngK3MXM3MJwD8L+SiphOE1LsuvvgUVM4C8HTw8dMALk/ooGKEz7kBtkSNVIaZ65h5dfDxYQAbARQg\nfa6f1/nlB19O+WvIzJ8GH3aCxJ8ZLbh2qTYZ5QPYZXu+G9ZFTRcYwFIieoeIbkv2YOJEP2auB+Qf\nFZK0kE58lYhWE9ETqeqKcUNEgwGMB7ASQG66XT/b+b0d3JTy15CIMojoAwB1AJYy8ztowbVLNZFo\nD0xl5tMBXAIpOjwn2QNKAKnj84zMLwAMYebxkH/OlHZZAEDQFfMCgHuCd9zu65XS18/j/NLiGjJz\nEzNPgFh/ZxLRKLTg2qWaSNQAGGR7XhDcljYw857g7wZIA8QzkzuiuFBPRLnA537hvUkeT8xg5ga2\nAn2/BnBGMsfTWoioA2QCfZaZXwpuTpvr53V+6XYNmfkTABUALkILrl2qicQ7AIYRUSERZQG4FsDL\nSR5TzCCiLrbW6TmQbrgfJndUMcFRUAm5ZnOCj28E8JJ7hxTCXSza3/balUj96/ckgA3M/GPbtnS6\nfiHnlw7XkIj6GDdZsFB5BiTm0uxrl1LZTYCkwAL4MUTgfsPMjyZ5SDEj2NjwTxATsAOA51L9/OwF\nlVm6bYgAAAHuSURBVADqIQWVfwbwf5AGjtUArmHmg8kaY0vxObcSiG+7CcAOAHcYH3CqQURTAbwJ\nYB3kb5IBPAhgFYDnkfrXz+/8rkOKX0MiGgMJTGcEf/7AzI8QUS8089qlnEgoiqIoiSPV3E2KoihK\nAlGRUBRFUXxRkVAURVF8UZFQFEVRfFGRUBRFUXxRkVAURVF8UZFQlCggoh5EdFfw8QAiej7ZY1KU\nRKB1EooSBcEGcH9h5jFJHoqiJJRErHGtKOnA9wAMIaL3AWwDUMTMY4joRki75RwAwwD8EEAWZPnd\nzwBcwswHiWgIgJ8D6APgUwC3MfOWJJyHojQLdTcpSnTMA7A92KH3W3B2zxwFEYozATwC4HDwfSsB\n3BB8z/8A+CoznxHc/5eJGriitAa1JBSl9fwjuMDLp0R0EMBfg9vXARgTbNZ4NoD/IyLTDLBjEsap\nKM1GRUJRWs8x22O2PW+C/I9lADgQtC4UJaVQd5OiRMchAN2Cj5u1tCUzHwJQRURfNNuIaGwMx6Yo\ncUNFQlGigJn3A1hORGsBfB/+K3r5bZ8N4JbgkpgfArgsDsNUlJijKbCKoiiKL2pJKIqiKL6oSCiK\noii+qEgoiqIovqhIKIqiKL6oSCiKoii+qEgoiqIovqhIKIqiKL6oSCiKoii+/H8JOwGA+hHabwAA\nAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9c6854dd30>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "rL, gL, bL = readFromFile(\"/home/mong/Desktop/map/map/videos/VID_20150123_172926_F.mp4.txt\")\n",
    "def constructTimeList(numOfElement, deltaTime):\n",
    "    l = list()\n",
    "    l.append(0)\n",
    "    for i in range(1, numOfElement):\n",
    "        l.append(l[i-1]+deltaTime)\n",
    "    return l\n",
    "tL = constructTimeList(len(rL), 1/fps)\n",
    "pL = detectPeak(rL, 10)\n",
    "plt.plot(tL, rL)\n",
    "plt.xlabel(\"time\")\n",
    "plt.ylabel(\"r\")\n",
    "plt.title(\"time\"+' vs '+\"r\")\n",
    "plt.ylim(199, 204)\n",
    "plt.plot([p/30.32 for p in pL],[rL[p] for p in pL],'ro')\n",
    "plt.show()\n",
    "# plot(tL, rL, \"time\", \"R color\", yBegin=199, yEnd=204)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}