{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "In the preprocessing section, we extracted the raw location records and stored them into MongoDB. Further, we removed outlier, and applied segmentation and interpolation to cater for the GPS outage and unavailability. We can now use these filtered location records for our first application of *Automatic bus-stop detection*. The bus-stop detection uses an unsupervised clustering based algorithm to detect stoppages on a route. We firstly describe the algorithm and then describe it's application for bus-stop detection." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## DBSCAN algorithm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Density-Based Spatial Clustering of Applications with Noise (DBSCAN) [1] is the data clustering algorithm which can detect clusters of arbitrary shapes in the data-set. \n", "\n", "### Steps to find the cluster in a data-set\n", "\n", "1. Find the points in the **{$\\varepsilon$} (eps)** neighborhood of every point\n", "2. Mark core points if more than **minPts** \n", "3. Identify connected neighbor points with the core points\n", "4. If the point is not in the neighborhood of core point mark it as a noise\n", "5. Repeat step 1-4 for unvisited data point" ] }, { "attachments": { "DBSCAN_Algo.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWEAAAFbCAYAAAAEHj8cAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAABmJLR0QA/wD/AP+gvaeTAAAAB3RJTUUH4QgKACwP54+/HQAAV+dJREFUeNrtnXeYFeX1xz+zld2l96YUEREEey+xF6TYYolGkxgTjSWJNdHErEZNNP40xq6JLRobogJi770LKiAWRFBAeoeFvfP745zxzt69/c7t5/M8++zs3rlz5075znnPe0ol0AO4FKgBDgLeIcwewElAPTAGeBPDMAwjUH4D7KfL4yNeuxwYocsv2KEyDMMIliq1hN/Wv2sjXn8RGA20Awba4TIMwwhehOcCXfXvdRGvPwdMQ1wVz9nhMgzDCBYH6AT8CZgKhIBPgaHAvYibYifgW2COWsY/cMkllxzvCfjq1avrGxoa1gC4rru0sbHxHju8hpE/nnjiiZ2bm5sPdxzn8dGjR9t8TgFbwkuBcxFXhGcJT9XfzyMTdWuA5sg3u677aUVFRQPAhAkT/vqTn/zkz/r/tXZoDSN/PP744+1CodCzjuO0A8549tln+xxwwAHL7cgUpiUcFA8BR9shNYz8M378+C7V1dXz1dBqDoVCfcaOHbvAjkzhUWGHwDBKjyOOOGKx67onApOAk0yAy4OH7BAYhmGYJWwYhmEibBiGYZgIG2WCCz1cuM+F+13oaUfEKEfMJ2zkU4T/6YKrPzfaETGKhSo7BEYRC68DtAe6AW18Ly22o2OYCBtG9kS3C5Lp2Q2o1pef1t8O8IoLdQ5Y0pBhImwYWRLdSDYCT+hyHbCHC+86sMyOpGEibBjJi25boDNSkySe6CaiDbC7Cx85UvvEMEyEDSOLohuNCmA7F9o68JkdccNE2CgXce0C7AK86sCKiNfqVWyzIbqxGOxCA2IVh+wMGSbCRikLcD3wAbApMM0VMfYs3a5Ibep80AeoUz9xkwsDgL7Aaw64duYME2GjVOirAgxSl3o0sLJA9q0zMmG3DmnXVQtcD5xlp83IF5YxZwTN58A4JDzs0QwE+FfA/cCvff/bAjgF2C6D/WsAjiXcyutHdsoMs4SNUqIncI/+pEtXYJQuHwo8AiwHLkHcHYcAJyMNCdLhHeBL3de/2ykzTISNksCVmN5tA9jUMmAe0AuYr39XAJX6eiWZTegtB36vy7NccMwvbJgIG8UuwLVIP8LKADa3ETgH2BKYrn+jVut+SHfw7wPa9QFAgwvvO+HPMYycYe2NjCAEuAKJguhSxF9jBfCOpTobucYm5ox0hbfShQddcRXckKIAtwF+hkyytSuQr9QeiZzoaGfXMBE2ioHtdeTTAYlgSMW1dRhwBBK+dkIBfac2wG6u+KINw0S4WJg4ceLeEydO3K/MvvbnhP2yn5OaPzUUY7kQqAR2cCUczjCyjvmEM2TSpEm/cl33Vv3zzNGjR99QLt/dhR7IRFkTsD6Ft9YARyKTeQ8Dqwv0K36LpTobWcaiIzIVItf9ISTLcZzty+wJvsCFx4F9U3xrE5KIkQsjI5PQsxapzna1G+aOKEyuA6YCnwD/V4bff00WXApbAFchYWpt0ryuLwTGA79MUbR3oGVGnpfq3NGFk1z4qRtMGJ5hmDvCCMwtsTepRTk46pKI5cK4EokRBrgFmJziLg3Qh6PuHscg9SLi0QD8FRikf98MPOl7/SjgRF0+z4Gr7cwbZgkbhUIq9SHaAzfpQ/uIGOssiLHcLknDYR7hScPPkhBgkDTpQb6/+0UI9GExXjOMjDCfsBEEq1JYdxvE1+oJ3/go69yEZMotQMpigsQUjwZmA+cTP6liHfBboD8SuZEMyyKWH/P9vb0+PFDr/Uo75YaJsFGsIjxdLed2SPpxLBF9MuJ/B/qs0C194hzphjhNt38d8GkK+/UM0Ix0+niSlq6SmSr6dcCr5K8mslHCItyAVK6aHWWdLvqatYcxgnBHLESSO7rEuN5i8ToShfG9iqJHvV7HK4DjgSH6/5HAAyls3wWei/HafOBUpOraDGCIC/MdmZQ0jIyoQGI9LwWGA2dGvD4UmazoBfzLDpcRg9WkFgq2KkUBRi3b04DTfZb3lsBdwN3Anogv2OO7gL/jUrXiXSQ6YriddiMoS/hIZPb5ecQ/d73v9e7ALOA9tQSMAmHixIl7AP8DVoRCocPGjh37Rb72xYFmV4br9Vn8GJfWXZN3IRzCtifiq52lVvF7Wf7a3V3o6YiVnMqXOAT4CfC4I8XvDbOE6UF4UqI24vUZagWfgfjMjMLhbGATYFhFRcUpBbA/q/LwmW8i/uNm4BX9/UIOBNhjuJtCXWNX/M2PIPUy7neht13GRgUwF/H5QutQnl8iaaV/R+IkjcLhTe/edl33jTIV4RlINbYTgdfy8PltSK3GREgfFJHLRhnjAJ2APyFZXyFkRnkocC+wGxIf+SYSkP9b/5svueSSsUhGEc8///xp++23382qCssbGxvH2+HNvkuioqJi1aGHHvpRvvfFleaeW0f8uxsShjZbLdRSxAVed5JsteTCXkjyyAQHnrar2KhCLp5z1RXhWcJT9fcbSD+uOqRpY0sFd5xFrus2AaxcuXId6h+rqKhYZYc2+4wePfq1AtqdaOf8HH2gAywGppSoIXOrC8OAfzpwe4KVX1HXiWH8cAEFhaUtlzHqGz044t//BAbq8uXEjgsuZjYDrtXlJqCdFfsxUrWEDSOIp/kGVxIcaiNE+MfANzqiisV2yPzDMiS2d2oRffUlOoJsA3xtAmyYJWzk0xrejfT6zF2HZLt53EHLtOFEdEdiiENI4Z1FOf7qfRG3y2SndaafYcTFCvgYQZLuXEBk4saWKb7/x0h9hx2RuPdcMxcJ4ayy1kiGibCRT1am+b7r1fr9Xq3Yx5EaEXW+dX4EPIhkyF2AxEh7LPQtL8jB9zwGSWy6lNZxwlulEjtsGOYTNgrBEt6g7gfPBfE7pE7EEiQscrlauHX6szsS/naurv+IincIeDnL39FRy7sKqQg3lJZRH17s8Cd2ORhmCRu5pj8yyZbpXMMu+rszsLkuR4qav8qZlyn3Epm1M0oG17cvK4Cvox0HV+LvDcMsYSM3uFK74SV9sE8CbktjMz2BXZH49P2ROhBeOcrbkUpqWyK1fSf63lerFujyHH3dy3Q/muNYyyNceNWahBomwkau2M43sto8zWH+FUgK/Xqk3OU6dUvMAL5UQY6sEbwJklbfFmmFlIvohI3ACMJtk86j9eRieyTi40u7NAxzRxi54H9qqc5BJtlSpRro6LNsq5CJr18jTT+7x3jfjoTbHu2Tw++7o/5uA2wVY50t3OxWljNMhA3jBzN2oQN7OFJD4kZksi0VmpCa1TOAe5Cwrx4+ge4c433vEa5n/HIOv/JkdUcsIHYmoNUdNpIaAgaFJWsYP+CKgO6Y4TW2F9IMdCpwJ9En3Rygg/5emuOvWaMPm0STge85LQvOG4ZZwkZWhLeDq52IHbEQP89wk68g4Wp3xBC6dmp13434jv20By5GWtMPzMLX3QI4nORqAlvssGEibGRdgIchk1OzXKmeBtILLpvJE9sgKcMOcGjEa4cAOwCDkSLq6XI40uJrV9//OiEFiY5HJhMT3Uep1h02fEyYMGHshAkT/jthwoRRJsKGEZuDfG6Bo9UadoEPEZ9tNpiOxOpCa7+sv8dcum6KfsDPkTrJZyM+XpDGuDW+5WSijCx2OA0mTpzY1XGchxzHOcFxnHGPPvpoRxNhw4hxvyBpxyF80RGO+EzfQcK6gmYREj1xJnBrxGt+4U9Uz8FRi9kvkiOAXxCOBV5JOOZ3LvAf4GOkUlwyldO82GG751Kguro65DvuzW3btnVL7TvaBWEEgiP+302Bbk6EIDqSzvxRlj56NdE7N7eJsRyNkxHf8W2Ea1KcB2yr1u+zSPcZvwA8DlxEam2VvNhhI0kOPvjgJaFQaIzjOLe6rjv6gAMOWF5q39GSNYwghXg9LdOJ/a/NcyVxYbMc7c6bwP36YHggyus1iB97FuJbBolPHorEOq9V9wpIGvUdAe3XFi7Mc2CNXTHJMXbs2Gf1QViSmCVs5JLptKx4hroBjgH6JHhvb+AsJGQtmbA3V0X4yhiW8iX6cz3wFOI2mQO85Xvdox1SxS2SHZBEklNIPhTPYocNs4SNvFnKrgsfIHUm6hEf7BVqlY5EJsFi1Vo4i3C/ugVIdl4mxscQXe6AWOhH07IWxLeItVrv+zuS3yBp1kMQn/EDSX5+dxd6WeywYZawkQ8hbkKy3JpV4Lwog7YJjILmGMvpEFLBXKeW78wY27xY3Rq3E71JqX/y78dq1SfLcIsdNswSNvIlxMtdsWZPUZFrixTeiRdlcB0Ss/st6TUM7azW6gZ1HQxFJuxq41jfM4G/xdnmjUjxoApEUC9HOo6vBP4P6TAdi1qs7rBhlrCRRxqBw5AkiKdJHGXwPRJ1MYnUawafiHTkuAUpEtQFiXxAf3dL8zvMQHzO3v7UILHFWyHJIomw2GHDRNjIGzW+5XSH5QORWN5tE6y3t/7uhkyKLUEmCdHfmTQGfRNpef8u4Uk9iD4ZGGVQwCWuZBne4gZby8UonpFhYFgBHyNpXBHQP6t74dUIUU7WgLgbmVhrBn5F68gLj5OQ9kiLkJTqpfr+XsjkWCjA+2kXxFc8NYn1K4FxhDPxdnJEzI0ywnzCRr6e/l8h0RC4EgK2R4rXYwXiV/WWOyMRFB0RX+0M37p3A08gKc6e3zlE9IiHDJ8tvJnC+s3AN0gCx0qit0oyzBI2S9jIiWXcE4m7TeWa3B6pWfEeIuQn6f+nqJVdDLTV7zEDeEEfToZZwoaRc2tgvgufEY7fTYaPkFTgDbQs0jO7iL76KsLF6Ie5UOuE/dWGibBh5FSIP1fXRJ8k33I0cJwu3wr8Ud0RbwW0S/3UxfBNDg/DIFeOhQmxibBh5IWPCGfTJaKHb7kn4vcNin2B3+ryNeS2ddIgF6ocqdJmlDheiFqDPvUjaUCKZvcluQ4ChpGpNRxCfLzrklj9QSQK4X3gsYB3ZSvZHRzyU+uhvyulLy1srQws4R7A+cCLwBikoInHTwhXktqH1t0LDCMbQrzOlRrEuxMO34rGPKTEZDZ4Cpkw85bzQT+gyoUPndQTVIwiEuEjkc6xzwPjI0T4XqSkX1+y1x3BMKIJ8XJXunLskKddmAn8jPyLXx/AUSEO2ZVRmu6IHsAy/bs24vW1+vt44D47XEaOhXgemTcLzYRCsT57AztYV47SFeG5SDk+iO6Hq1FxXmGHy8iDEM+gZb+4ZBgN3Il0ai4Vn2oPYGfXJtNL8RqnE+JXm4oMdz5FKkzdq+schwS/T4t8c2Nj4/5ISBCTJ08+e+TIkdcAOI6z8i9/+cvTdniNgMzRSsQ/3CHJa3oc4XoUF5A43KszMh/yHdLBoZD9r0uAt53s9Owz8kAVEuR+rlq7niXsz3u/P4ElXaHC63jLruvajK4RpKXQ7EpNhT1p7TKL5kKYhdT2XYOUlkzEWcB2uryS1FKPc01nYFcX3tImqkYJWMJBYWnLRrYt4s5I6ctEvtF6FdWZSAnMRPwN6TcHEhP8Ugq7Va/79JWKf65Yjghxk10ZJsImwkYuhbgviUtXpkof4FhkIvABUotC8AR8I3AmqRcF6o8UMlqEdPBYl8J7VwFvOqm9xyhAd4RhFJPVMNcV3/DAADf7LdIJIx36+u6lXmmI8MnA1rr8Fall/bUFdnNFiNfa1VGcWMiLUYxMQ9ojBUkndUXcTWqxybcRbl9/PKnXRV4aYzlZGoDdXWkeeooLe9nlYSJsGNm2hl2ka/PKADe7NzBIxfioFN73JlCny5shtYEjdpeuxM78uxm4A2mT9Eaa+16HJFrdBrzkSiSJYSJsGFkV4o1IxERQEQIzCfuCU0kQ2Ui4ats3tC7Mfr6K7NUxrOS1SN2L1zPc/z4+0d/Fhc5u/JRvo3Cu5cCwiTkj56QQMZEMv0bqoyxF2iAl6j3XDthSBbweicTwx+/WIDHLHueQvQzAEcgE3xzgBiRqwkUm75brzzJguSMdPYwCwSbmjGK3Ipa4UvJx6wA250VddAK2AZ6Ls26lWre9kN52p9E6gaJJXQy7IZNus6MYQceqG+NhpKh9ukwFfh9l++30x5tAdN3WwrzM6lKYCBtGJkL8jUZM9M9wUy8hlQOXInWNY1GPZJL20r+7IRZ5tMSQK/X1xbS0QLcFxhJOEumjQp6L0W+kMIdc8a8v8axlYJUjgj0AuEsfOr9wxOo3TIQNoxWfIiFbXTPYxgPAMzqEj5cEcRJwiC5vRARscAwRdmmdMNIWuIiWPuJVeTx2FchDzJ8WvsEVMT6dcMTFBUhInRHwwS8ZJk6c2PWhhx6qsdNaltawVww+05KrS6II8IFqpW6if3ePMGQ6IanPiSbC2iAdQLxi8ah1/JhazIVEtT7Q/A+HGXalmQjHE+Czge/r6uq+fuKJJ/qV/qlzK8C9AtxJ4O5qlzJoLYV3CLamwjDgDLV8z9f/+VsveX7gpcT3q3ZDMuJuAw5HMu1eAC5HoicWFehhnQxcphbw1XaVmTsiHj9T66JXKBQ6BLilxM/dKKSxJUj22FC7nMGBVa60O9qZYKJ/KqIsf0s4Y+92JE73NaJXX+sMnKJWpTfc3wO4Ry33gn/a64NtmnX3MBGOf6W47oOO4wwHFldWVj5bBufOP0xcapdyCyFe6Er5yiAeTB/rA31zwn3s/oWUd52HRDz8SEU2WrGgYwknT6xVl8QzRXhY29mVlbXrNTDyHif86KOP9t+wYcOio48+elXpnjK3EsnsGqQ3d39gAvAkOMtL/YLVxpc7A/Oc1iFfketuDWya5V26EokVDiHdmWdHuDJOJdxE916kX53XIGEL4AgkdnhcjO131e1/RLAZgqmy1BFr3zBLODaHH3741yUsPw4SxjSUcE3dNwnXvt0N3PfB+b7Er9krgfOA9S7s7ojrIZ4V2wB0yeL+dPO5Krr5RNgBLlQL0kU6Qz9CyzC1c5CJul2BL5AuN4uR0LchSI2Mq4H26gL5DflzCZglbCJczrhd1apqn+Bc7gTux+DMLuGDsY/+rkV8q+/HGeaFXPG77onE9maDm9Tl8DlSz8KPf6LuZVpnqvkrn52CRF98oAJcj4S8eee8Ny0bL2RCpW43FTdWlQttrGymuSPKTXwb9IbsneIbvwBnekkeERG829XiPECbgSZ6T3vEdZNro2MQcIC6EqJ16+gJHKxiGC3+thnxH++ElLgcF8A+1QP/UMF/Arg1hfe+6RRuFIeJsIlwoFJTrTfwQNIPI5wDTAWn5NJRXXBSnal3RfB2SPOaH4xMvr2DTMhlY0R6M9LME+TB4iK+/sm+9doi7qjPkESKdNgGuFSXm0itYtwnTm67h5g7wsi5vETz+6bLJkAbcN8Dp6SaQqYTKuXAfFfEa0iKb61R0apXq/VkJLW3Ri3VIIrhbEQiL36tf3+NxBFHuhCu1lGRV6sindZGn6uboyfwSorvbWv3qIlwKQtwNxXf9gFutJsMw923wSl7X54Dn7sywdQnhbdV+O6TSn7wvXO+iuDFyKRapjyFxBH3Bf4bQwB7+85rrFoViViNJJ90JvXC+DY5lwWsnnD+xbctuNsDuwQswB7qD3XNihE+Uks2WdYhERmvA9ciPtH91BJuS3CdLJoRF9QewI1Ih4/hvteXq7W8CngyTQHG54aYT+ojCruGzBIuKfENwu+bLPVyc7vvgrO4zK3hkCvF4PdEEieS4V398XhbH5obCS7rrQrY3rfcCUlv/ti3zh36ky9qXaixDs9FKsKTJk26zHXd84DnevXqddgOO+ywoTwPuVuB+GuHkHo/skyoFuFwPwLn23K+6B1Y58ok2+6k133iBeATpEZFUNmKG4GnEb+zx7TkvxJ76EPlRVrXNQ6StkiRIyO46zEwYkZHPPTQQ5V1dXXrPNEPhUJ7jB079vXSP7xuBRKQPwD4OxIXOoz8+9ZmgvNZuV/8riRF7FBgu1WHVGlri5TnTIYDET8vSK+5u7K4f1McaeNkBEROfMJHH310M+GUx/m1tbXTy+T4HgNchcxkP6hD2EKY3BgM7lYajVHOFsg8stduKF3WIjHQn6bwnm6+5VFqFXv39yYBj7hscq4YRRigbdu2hziOs++GDRu2Ovjgg5eUz33+A7UFtm8DxAp0y7oZpCM1cr8r8q8xiXCERg1wvC5fjEzyXRugENvkXMDkzCe8zz77rEP8VeXEA0ixmcEEk+3kUUkw8alat8B9B5xynmyZouLSvkj3fzlS2e1avTZm6kPfa520if58aSJcxpZwmdIBaRZ5E9HLHKbDwUhTyBtp2Y4mXTrJ8NVtKGNreCMyUbe+iL/G18CZwF9VkNcjE4ioayOoeiL1rkVVmQgXEVtlYZujdQSzCZI0EAQNKsSdy1iI1yLhZsWc5j0XCaXzRkn/VNfEHwk2YqJkHtiNjY1518AK30HtF8dlMSQgq6uMcPsCHbOw4Q/19zqSD2FKhhokhK1HGQvxEqR1fCFSR3qTYtmoQVz0LolHH320/8SJE7/Yfvvtl0+YMOHIfItwDyQ3frgOZ1oMPYA/6ck/14Q1aQGuJPUaBcnyH+Bs4FdIjdlOSPTF8UgscCZUAjuC27+MhXgO8FUaluGVwH1IoZ+gGYyEnd2Tpe2nStFHSFRXVx8DbAa0dRznd/kW4SORSk2TCNdq9ThKra32wBUmrkmzmVoufsshWnTEpsC+aQzvdkbiQgcDv0CaUB5Dy0D/DHSI4WUewjaN1Ooq7IJ0v2in5yEduiIRK7G2X6cPyUIQ4aK3hEOh0Mto5p/jOHltN+VZwsv070ihGIj4NaeT33TJYrKC2yDpyB57IwVZ7lZx9ugOXAP8DmhM4QN21Bt9J6Sdjj//P8iuCwOA7TThpNysYRcprp7sUP4Lwh2e04mBH4LU9b2O6KUl30bcT82kXvnMLOEojBkz5q2KiorBwPajRo36az73pQpx5nfVvyMrbS3U17/DgrSTZUtapsLurX/XI21svDChHoRjNzdJYfurIpb/o783A0Yg4VZzAvouvWUf3ffAKas0cwc2amrzniSOsZ0NnI6E/KXjU96RsCtpe1qHM36GdBOvJvk6wlvoqGZGFg5PvQsVTnFPYnLooYfOJriokYxEeBzi9+2JVGnaDimpeC+SivwHtZgjW7dw6aWX7hoKhdoBjB8/vscRRxxxIEBFRcXqiy+++HXKDrcDrcskvglsi4QMveP7/ydI14RhSMhZMnRUgf27WqpP6U35NZIphbo2LgrwS3UlXA5zbZkJ8RpXHpr/hwjOA8CaGKvPJ73KZlsjhXpAIhgmxVhvTQrbPIDw/M4NBN/duUKNilUYgYjwUmTSzd+/aqrPEj5Hb+xHo/hVeiF1SenUqVM92tk2FAotK9PjOYyWWXI/Qgp1L0MmP7+McB3ckMK2d9fz1Kzbui+OdZyN4eceKsQrYjyA+iOt4TcCp4Izt0TO6ZnAGJ9I3hPw9nckHHe7Bvg5Ei73QQbbHBzh6ngmS9eEiXBAIuwJgt8VERlTuDramxsbG8f7/jzw5Zdf/nf5Hkq3F627+h5MuCzhtmSWseRV/KpEJmqmRljb1yF+5klZ+oJt1CKO1dH5QuAgXT4fOKtETqz/XmjOwvZfBfZXI6i9/vwkQxF+Aol2coCJWTouljkXsAgbmQlwBeILjuRdtY6byLw32ctIVIS/GJLfqn4+R9fLjuBOiWLp+n1rX5fQyb0EcSVtQKKIuiO+2X2AFcBbaW63QR9aC4AT9TppVOH83He800my+FpGI1nF5ogCwhp9BiPCg2KIMMik2yqCqTvbgPgmC8E3G1EO061SC64ZuL9UG4wiE2eNwKH67//TB2SqnE+42tkVKuYDEB/8+2odn6oifWEa189OyETwS2SnOelypzAiNcwSNtxaWoakRTInwA9bXUBffLCG400Fx9VmoveU8pl2wHUlY7GT79/d09xcxyjLswh3M/bS0/voCOipFLbdCUlVrlSh/ymtI58yNd5+7MJY4J9OsNe4ibCRMluQeaZasbIpUFeKHZ3jqE+zK4kyd+pD8ak0N3U7cJJaui9Eef1DpJTAOlKrLaxGe9TloNiRcBH5rdVqN9LECvhkZgW3UyEqZ7rJsNetLZcv7IhA7oLU6002oWNzwsXXq5GJ2k+RJKhoZUTvIJyenqqluQy4HClpWQX8mWALu/u3VWc6YCKcTyJD0sqVjjLsLZ+Ozo5YqG+RXPnLExHf8S2I6+oIJPnip8Bxcd73Bal1hvbzHpJsU4VESgTZxul14H/IPNDJdvmbCOfLCu5By7YyuWAgkkxTiMLvdXQum3KYjrgj3iKcshyLET4LeBgtuzzXtdwkOyFZekHcm/5OzT8juO4uLpK4coqTnYw8E2Ej4TXoEDsaIpJBAQnn1kjnhEbEl1iIVKtroncZCfEKEtchfhyxnOci4YWPIL7k59Si9DgQyV49j5aFgDoihba29/1vONIs4EpiF4D6l2+5p1rGRoFhE3Pp0Z/k4iS3R/yGjt54d2fwmYN8Qr5FgT/Ytwe3HpwvykSIF7kSVrZDjIftqyq+/kmym6Ks1zvG8vlIIS1Xlz8Dfgn01df/hERDRLIayZY7EPFjB90l2TUpMEs4H1ZwNS3TQuOxme+mHJThB7+ABOGHdFvb6v/9xYLaIGFD+0YRgy6I/3GXHB2oLcupHKYjdSM+TkKweqprYPco60wAPkJKaT7o+38nn7vCW/bHiscLk7sBqcz2F4LP+DMRNks4L2xB/JnmgWolT9Xh5o5IV5IHM/zcpWrN9Ed8e4fo8P9cJBnkz0iNg0N84vys7/0XEC40fz5hX16lbmMo4ud7MsBjNUAeDO6H4DSX+oXhwGxXzk28kcqFeg5dPe6f+15brCOnaEJ6HFJs3isCdSVwFVK75X8Jdq3JJ+LbIVl4QSRwmAibCOfcCo7XBgq1Thv1Yn8YqSN8XoA7MEWFtlIFeaRav22QYkEdfet2jHivP3LB31V4K59V9rOARRiglwpxWXR0dmCmKw/HgTFWqfMJYpuI1yqReiN1SA0QL8HiU3U5RD6UT0lx98YiTQA8Yc+0sI+JsIlwzhmWwIXj99sOzsLnf4C0MqpB/HteW/Nmtbxf0nO6Eini4udfSFr5LKSmhce3erO3QUKisoHX0fltcFaXwXUyTY9ntImwq5EQtc+iuC9G+oS1M3BbwPvVP2KUYpawiXBRWcFdkULs8XgBydfvSPI1gv2civhzXwRujrGOv2btY4gPcQ3gVTaL1SVgBlICM5JFSLnGgYSbiAZFvVp3K5EZfK8u8fJSvlJ86c3VtA5j/Az4W4y3+kcrqba8Gq5GwEtErzNRg/iRlyINTSeYCJsIF5MAO2oFJ2IxkuWUDl3UEgLx645D6jkn4usAvuACUuuplgxD1TVTjSQqvKaWu1cOc0EpXzEOhFwZcexKy1oT8XhcLeA6xJUVKaLdkC43keLXTx++FTLi4Jwo2z6CcNGhl4B5dl+bCBcTm9LSj5oNlukN1ltvkKVFfsx2Iezz3JNw+U2vo/PH4Mwu5YtG60y8Iw8ehus5jnde1wA3Rvl/G6QfYV+k39zlEa93JOwm6xJj2zUxls0SDtpkc11n0qRJ1wJ7Oo7z91GjRj1sIpzZIa0iN3G5zWpFD0WaRRZ7QZy3kEmmaiRONkKf2FqzDn+NRI+cC87XpXb1ONDkSt3gC5AU5/PSGL1sSjgmeCc9pv4svanAeL1OY93w41TMK0kcTREp3k0mwskzefLkPZEmvLiuexsJXJMmwonZnODSPROxFsm+KgWmIdEWVUhWmZ/d9SLdQDjpZSNwbIleQwfo71q1iFMV4a+RYjyDkXmHDVHE8C4ZYXA24ha7VH/7rexUJvqq1cWxpbpJ7jARTtoSnqsPrhokrDAulqwR/3DWEzvUKGj2RQqi3EDyPsRCZ00UAQYJlWpDy6zDhSV8IV2vN+V8pBVVqjSpBf1TpI1VLI5A3GYDkJDFTI2PoTpqGR2hFVVqkQ8yjWjNqFGjvnIcZzfgjA0bNhySaH2zhOMzNIcPqrFqKW2KTOZMDmi7NUiKay8kbTobYWibICna75Nc2cUpSOLIBiTC4xNk8q4kceAuF+4HdiP93mwuidvdT0cmkJuRKIxMmK0Pxm5IpIe/NsbvET//ha74+o9zZD7DCAvx+3o/JMREOPY130mFK1dMUQumSYfyQbEX4psF8b96ySMD1V3wPVJgfH2a268F/q5W7dFIMkCiLg73IRNWy/XzPyj1+GEH1rvyfbNZ7vO/cixZisR/p8qeSPLOs/qwPgPp7DEripXsjaT3UhfIuaYZ6WEiHJutcvx5dyJWxRIkdjco/LPxS3zLJyP+SRC/VbqWdxufsLTVG/JqYk/meHwesY1yYGWWt+/qqCLd0cy56n7YRR/Qa2OMnO4DTvedt69MLkyEg76W+9I67Zcc3EAzs7Dd95HkgF7A0/q/XhGCkEk43HLgVuDnahXvgsSqvpCikJsIx6eS4Avw+KmIsRyNl9Vg6KBC/bBpholwkFpYSbjQTangnww6Bjheb57xapG+2Xr0zC7IxFoyhV4mI4XLd9OHSar+QRPh+OJ4IRL58AzR44iDYDaS2j4CqR+SKPKh2cmsNKthIhyTQRRu36weyMTLRxGuhVTYVX/XIRMvr0dZ53jEv4u6FpJpbX4tEl43l9S7LZSLCK9Ra7Yyhff0RSIRQOKN7yJ7Xbef159kR26GiXCgFvAOenEPyMOHt0XCjyqQyZVoYV31SNZUO2Qy69dpDk+f0fcupmUhHz+b+JaTbWS6HindmQ5lIcJaU2KVDuOTZT7hTMrPVMg9+iCTZO9E/D8nN4xpholwkNfTRcBlSMLAn0m9xXimHEW4DvBGxMcaSUfCcbXd1JJdlcZnTUZqB6yPI+IPqtW9mmBC5ToQ7i68sIwtYc8lkYoINwG/QyqgfekTv67APxE//GcEWzLVMBHOOfv6jsdWeRBhf3hYrPCu75DU093Uml2Vweclspq+0hs/IAOQq5DJwBVIC/fIz68At6Yc6g2Tnl94Ha1dPL0IZ3L20+PstjymjFBL2l95r7sK+kckjmAxS9hEOGdch0x8rFIrMdc8giQuVCIpovisna2QYPnlwD36kyvaIbGinZB27emEIrVBWvqAZHN1ifEQaJOhKBQLKwLazqeI/3YrJNMyUhRPR9KlNyJV1WapAN+gx/o9opc2NRHOkwg36A0/O8rrPX0Xz4rSPAzOBHC3RAqk54MmtXIjRekadUNk4gPOhAMJT+SdiJSmTER/pEVPDRIa9ykyiz4SeIPYGXW1ZXLPBRUrHCJ+CvMQ3z08SEV4U8Kuny2QhJ2TkfmBm4ifZFONpC/XIlmOKzACoQLx/V2KBO6fGeVEno0UVhle4seiJkvb9ay/VOlIOFa5GzIxV4PUBMhVPQv/MDbZ+rMH6AO9PeH6teOR1Ok74ryvrhxuOEdCAzfk4KMeVMGfQTgEcQoyGbsCmQD+hd7XexMuMhSL0UgCx3HIJLJZwgFawkciky/P681yfcQ6U5Cwo7dK/FhkwxIbjnS5rUImUVJxdXgV1QYicZsrkbbmu6oVdDbZz1R6XR/QHZEA/WT4GBglesPUPB//QmUVmRVpqtJzEi+z8hW9Z+sJzx9soGXnFX88/KIUzk8bE+FgRbgHUig62o0wD6ljsBrxJZ1mIpwSO/os7F1TEGEH+AfiCloBTNT/9/ONYDYhN+mi7+l1kmx947cQP3IVrWsOxKPcIiTSFeEGpNhRbyQDMlbyRg+9hjog7qDxUda5BUnWWYwUnN+K2GnPjyNzBG3UijYRDtAdMVeHj9DaJ1StN+E0WsaOmggnx2vIJNQGpG9cstQR7mfXHml5AxKo/50OL5MpibgtUkx8vwxcNFciE4cnpPC+OSkKcDmKcLrsQ7iB6N5x1ttOrWVH3xONdUhD2HokguUKpJt3NLx6xP+i+Lu+FJwlPA5pp90Tcbhvh5RwvBepqrSJPilbiUhjY+OIioqKBoD77ruvy/HHH78rQCgUWtvY2PiRiTAzkQmtKlILpl+DREEcgkxmzdX/N+vy+ySOJKgE/qCCvhviG0y1stYQpKg3wGFI4ZZsWUAmwsmxu2853kjIi6jpQGJXkr8u8OZJ7odZwgGK8FKkelKtzxL2fHkP683hID7KSEa4rtsNoF+/fl1c190FwHGcJUgcYpHg7qSCMz8LG28ivdCrR/THPyo5X8/TjojvdV6Cm6RJRThEepNBs5D06M4q/Nm88cpJhDOJLJhHuOnsI3HWmw+cou6LxQm2+QQyf1GlhpiJcA5xAtzWQ4TrDRQR7pk6xHJ16P1Gge5opbojOiD+2VOSuLk2QxJRPiDJAtMRdEV8lw5S0jCUzRMhYuCUxc2tfedSjcjpjkzIdlYBfjqPX2G1k1qlPCOOJVzu7OF7IA1NUYS3RML6liM+tWz6ypqRSIsfqaAuTuI9X+pPOgwHLtFr5BayU2Yz0iDwj8bKwSWRaujiGL1GUWs4ngjXIHM+6R7POsT3PCvGw9csYRPhwLjOJ8SpWhbHIlWu+iKJDQ8msErX6QV9AOKjfSfFz/uK7EZEVCLZVV4zSu/62IHg2i0lckmYCMdmTsRyO70G1+tItAo4SUcww1SE/+YbBe0AnAos0P/HSn1vq6PDrnqNXmYibCKcTQPsDaAPuLulcVPMQiIQvOVYnIIEu7tqwXbV5d+SeufdbDLE90DaCvEpVyLtbnKBTc7FphrYRkdbLyMT6qcRbl21Xo/fQRHv28cnwj9Rl0Z3HVE9EeOz+hGOmNqe1nUp2gFjXejpwAOmISbCQbE6DRG+G2muuJz4NXT39Q25O/iW69MYsu+v25hM8OULv9Xj0KDb9qq2dYjYhy2QiZ9lJsJpk8rk3AC1Yr3IiD2Q7EM3wjL1Xw8h/Z8/lPFLJBIiRHw31ed6PQ9BYtQjrd6zgJ2Bk11Y5cAkkw8T4SBIR9Bcwoku8ZirwgUS9L4pEnudakPPfQinlvdFsvCCZJlufzMV+531//7ojjP1tY1IJt8dBFfTwizh1uwEXKQPv5C6GLzR03/VAl6HRDW4SBRMFZIxt5GWxf9vQdKWv08wAnORSJwaokf2tPctdzXpMBEO0hIOir5Iz7UlwH+QwjcHqKUZrZB6Z6TEI0gt4RDir4sUt7YxloNkkf58hiSGLKFlpt/WvmtntFpUQc2SBxir7XZE4qNfB2d5AV5vzWqVLiF+l5QhhKOYPkXmLbwH/ypa1+N4LMFnxpqHqFSXxxX6EH4AuD/Gurcgvuh3gf+ZdJgIF6II/0KHj94Q/7EEN8ePVTBAkmYG6vvOjdivp5FuCl4qajZZhnSAjmS8PjA8YdgY4GcGVMTHrUUyPTeTobU7DJwNBXa9/Vsf1OuQuQF/zPdgHTm9hvjjd0ZcV/8juFrX+yOp9G8DhyNFoiYTTtw4LI4IzwL+4CTX9sowEc6qOyIWfstrHTIx16wX9doYgue3olGxHRZhuawHbs7zcVriE+C5wKsBbjsoS7iHCjBIBlh3Us8WzDaH+Fwwg30i3BP4u96b+6tAnx7wZ3dF3EoOMvHmdVfeTq/FjvoQI4HLwjARDhJnA7hNBFPS8jYkhGiZisFo/f8GxI8XyTgV7jaI33cA4rebHm+H1XLM9OFxkt58E0muR1xNxMMmmZuxDVLKsrNa13PirBfEufwG3FuQkot3gvNtAV5wbWIYAJ1992UPWkcmRJ6Ls3XkdI9azsnQhLi8KglnU9aqe2Ec0rXjS9MEE+F8WcM1GR7PZsLt5EGKnHvEmsBqBp7SG2qAb9gfa/LGK7iyKZJaHinsnfUGnpFAJAcipUxBwp1eSsK98CoSxN+L+HHRfg7UH9TqaowtKm4FOAFk5jmnUdhV/2arxenSsozkdOBRJETwkQTnbwfCbqxfpCDCrs/6rUJqjISQrMi+SL+6NnqNfaouikFIadOVJhMmwtlkNeFC6qmyB9KXbQUym+0NLx9QYdtI/Fz/SCszXlv0YYS7II+MEOHeSNREG6QX3Q1xtrNUXRy16mZIJsohRGxfYTLumWVJWIhryuBaG4tEICykZYy5S3RffDTm6Dmr1AdzP1p3x4llbCxUN81y4BvCSTIH6AMWJENvKVISsxKpxmcNRU2Esy7C6TJSRbSrCvLDvgs+2b5w/9abYRnxs/dmEi6sE1nScrBvqLt1gs9bCkxAJga7Io0hp2ThuL6iw+rOSFhb2YuwA7NdiZzpk8Fm5ujIZG8V4aOQWsOJaFYx3RYpBLUuwhI/TM/Xp8iEnWcQ9DKJMBHOhTsiXd4jnGWWrpAtInbs7zBkkmYFknJ6jg4dI7tXvKvDyk2TsLwhXCi+EsnIyoYIuyRf0L6cYoXrA9jGx4TrCn+dwvuWEj208C0kGaMNEqbo6INzCJIabZgIF6wlPB5JD11Bdgr5/BiZOe8JHKM3Xi0yCfhExHc4m/gTOn6eQmpFrCW1wvPZotxEuBuSqfkZ6UUcPItM4taQOKIhWWZHPEBvNmkwES4GEY68eIPmSySKwVXx9cRqd6LXAEj2hn4PqSkQojDCjspChF0ZeWyJTH7VIJmU/0lzc1Ps1jURLoXboh0y4bQtcDvp1d/NJvcSrlOxHonvbCC5sLJENBfQ9ywXS7geScjwJmNHFN9zxDARDpZjCLdoP6kARdiNGG7+TC3i1QV6POuR9liziV/cqJxF+E1kQrcH0RtxJkMXZEJuEZKV2Wy3solwsfIJ4QIpXxTB/nphb9Fu7kYkPvk/5K/7wvmE3Sdnk3zwfzmJ8DLCBZkiOQkJCXsdqScSi18Du+jy8oBGRkYOqbBD4OG8hfRuOxZp01QMxcXbIPVkt/X9b3tkJrsN+W035YUzOYQ7R5sItxThWHRAkmg66uisZ5L3cFWOjl8NcKAbFn/DRDgwIf4AnAfBWYq0OSp0IT4d+A3ShsgT4i8Rn7Fn3WfzRuwfZzR1G9IF5DlalvvsCIxC4pljbNftUwYXWzwRXoVUsEPdDPGqrN2CRLhMAk5A5jX2yvK+/xb4M/CGK0WADHNHZEWQV4P7JpIWWlugO+m3MD1r6TskXbcXwVXciqQaSQroh0wW/oHWEzXvE92v/heknkazPkS+873WCbga6AbuJeA0lqEI90OK6Z+H9JObQbimr6PXo6NuCi/l+SZkjsCr8zuS7FY46+vbnyG0ThgyzBIOTIhXIcHrTQW6g3cg8aUv0zLGdxESxB+tBkN34EYk2mL7ND+3O+Ekjy1pWeQ7mfeChGhFdjIZIQIMSJnHcrCER/hGBdshvd2uUDfE27RM+R4LXID42w9X14PnPvO7LLIduXAX0n1jPMnXDzHMEk5biFeoRbwrwVRYC5IZpJ7Lvw+wiS4fTnpRIPOAD1Q0XokQikTciMzmT4viLvlEHyBdKeFi4a5cR9U6pPfqTl+JZDk6vodbJH6h7aXie5T+Pd/3WiZxw3VI5ufMOOf1Q+AWR5JEDBPhghDiPRCf7LMqjBVIDG8hVpyaTjgKJF13RQiJwEimO3Jvtcy8gkZv6E80FiCxs+3AmV2KV5Irlv4Rer3s4Htpc6So+q7IxFy0lPPx6gpwkJKTYyIE+m0VyGfSvdD1YdAfyfo8jdip/EtNFwJSlwC39RD5nY3PxS3UHvHJVUdYJLfosVyJTJRdqeITrcxkIeB155iW5c/ZC6lxAZIZ9nqcoXl7eYg5X5ewBVyhD+nN9eHlRTKsReo1LEhxk+2QAvDeyOaBDEcQbdS94OnCWUSvR7HKKYwU95LAfMIpW8Q/BNh7xbQjH2TDVIAhXEM3U+oD/iLf5kCAUWvP0Z9tY6zTF4lnvg1Jny5l6gh3/GiDpCo/HUOA+yBdujtE/K+T7++VwIW6jYlIHeJMWKeGwxqkuM9ss4LNHVGIbAOc6htCnobM6G+nw8B5wGJk0um1AD7vt0jQ/qfAnyiujCivP5qry9HYWl03IP7NK0p42LnaFTfOGciEb6xaEZ2Qanq1iOvnUrV2T0YmiS9CJmRB/La3IK2Q9tbjnEnfv3v1Jx4mwibCecXfMNK72F+lZa+1U1WEv8vws2pVgD0Luz/ZaTuzOeKLfZXYkzEd1KJdFmc7ByIFhV5ESldOQ2JXIXqkBkgq9uHIZNy/S/nCUXfEq8BHCVbtQTgssgKJlPD+rlFD4DPf+mOREDVPwIOe1KxHms721BHLyyYDJsL5tGfeAPdnSHbdZKJ3wFgfgAB72/GiEOYgjTX3RNrdP0/yoXOdVUAXR3mtF+Kv9Ton/D7KOtup9VWh60aLC+2CxP06SNjVu0hdi0StiqapldegSTKlzGDEj5uIz/ThtIPvQbUSiZhYFWWE1dG3fIg+sK+lZdukTNjLty8nApeZDgQvwg1qicTyAfXTYXaTHTIA527gbnAdHW5vgvjmNlVr7rUAP+yviN90ngrwb/X/m6hVkoidgD/qcjQB7eJ7kMRKL96Z8GTkaUg7nG+jPDDWE+6MkUyL+Q3Ah+A0lfq15crE42bJr86lSEx1LeHGqG8j7qjIB9t4vX+31PM5HMlKvCuAXe+n59PVB+xMxyqoBUqF3niX6omLVkykK5J62tsOVysxdpE42+10iNjZN/wOimZ9ODbpufCfl2TYSUW2UsU0kk+RFkdfELsf3Ws+10tH4Lgo66xCfNYP6O9kRHUqOGtL/SpxRbyGk/pE+Pe07Ey9IcbIYpk+YP2haXMC2PWDgeuRCJe7kF5zf7T7PnhL+EgdVj+vT9TrI9b5BfmrxFUMQrwB3EmIz6wKCXLPFpPVMmlL4tC3nZDojQ8J1xJ4MYbVlcgX+zHSfPIU/XtejPVmpvD9vwHnuzK5SAboAzrbPKAP7HWIGytTvBrHjj70H8USNLIiwj0IF1iJrJFwEBKqcoIdqrhC/Dq4hwA/Irt1iFeqxZOILRAfroPEpXrnL5Mh/ySkkEwtmU/MrCF7dS0KzQqu0/ORo4+LmQgTj2pk0vcbwsWfQHrLba3XjRfdYpERWRDhub6hbWT20yBk4mYoEj50E+XRjjwdIX4O3MWEi5vkk/aE45fbE4y/1VWrek0A2/kAnI1lcmFsTWFPgDvA5UghnrnA73zXy8f6APd8wGucliJtBCTC4xAfXk+kMv92Krr3Inn+A5EiIesib8C//vWvm7mu2wbg3//+d/tf/vKXwwBCodD6xsbGL8rweE5FQrna5Xk/3tPzOpDEHXK7Ik1Ev9fhZsh3bXQCFurf5yPp2VOQSmihNPftszKIgvCeNn0JFyQqVBpUgFEDogct/cn+STizgrMkwksRf2atzxL2t1H/CnHQt6K5ufkAx3F6AGy99dZ9Xdc9CsBxnEUUR3eKoI2KZnDfRXywkdbPJsjk3bu0LLaSrWHpPUmuezZSsMW7yV7QG/MaHQU9rdvaw2fZ9SG1iZ+2SCzym+VyXWiRnmE5/tht9OE4NYX3rEIm3vfXa/PbOOsuM8nMzlAkKMqgdkTSt2BvWpaJrEeyoxpU6E4m+aymOiQaoRoJwl8Z8Pl/mHBRon8jkRLbIoXi0QfzLrr/2yMlDH+lN3s1EnVRQct6Gn5qdUTVHfEDjwAnVPJXgIwoc1mcfjThidMbSX0yvZrEYYWvOWYNZ8USNoJ/tn0Hbid1B3iWoJea2xGJpV2V5MaOAA7z/X1rgDtaEcWdsgGZhPsJMqF0BzhTwN1V//4CnBiV01wvFK6KcFjcloRrCA9TF8fiEhfg7jkWYJAIDI+Babw/kQCHSK1kqWEinHemqeB2Rvytd+uQ/tkUBJgIiznouhHNSGbVSOAdJITtc3DWgzsc6BYOI3M2kLBdktOs2/RNBLpvqTvjx/IAcUpdgCuRmOBscQRSye95JHrB4zEkGSSko5mgWe6kPw9gmDsib7dkLeIfriP9LKMavfGqkRqzQUenhBD/7szYFq6RgggPR8K9skF3wjHdISRyYVUONGI0MuF8oVOWcz1mCRczFYhvdUckkP6hNLbRpO/NBt8BM6SfnhGAAHci3PYpG6zWh3A9sILcNKLdCfilLvdBCjQZJsJFw16Eu9EelaYIZ4OFwHRwzMcXnABXIJEjThY/ZjXSY24bJJJhY26+WtRlw0S4KJhCuGfaWwWwP0vV8l1kpyYwhTocKcr+IbnpNjGb2IW2ssF9eg0PRWpHGCbCxYQzH9wtkWD499UqbsjwfO2rw9BXU7BMViFJEt/ZOQlUgCuRSce2SEH1b2hd77kSqXw3BLif3LUF6oAkXnyepgUbAqY4kkU3w862iXAxC/Eifiht6b6HREhUprmxE5BJOlTMn0yw/lq9Cb/Rim9GMOLbGSnc/yUS+dJWXQPRYri3VoEGqbe8HZKNuiCLu9gDuA7xHT9L66JcidgAvOuUeCihiXB5CvIKcKfojZgOyZaybFKB+KocEiPywB1INwuQbNNtkDoL0SqMfafno0Z/foREy0QrjN4PKR/wHpmFIw4m3Jdw2xTfuwZ428l+1IVhIpw3If4W3K5IAfhUuU+trnVEjwXdiHTH/byMCuTkg64RVqd/wrWbiulUFd/5SLeSsYQbv66NIZxXIRN8z6slmy4fyAOYTUmt+edStYCtSI+JcMnzMVLdrGOK75uHNIqMxIv1/UwSLYwscwZSeWwRLTuVdEUK49epEHrnag5SgXAWEsb2eJRt9iecwejvwFGJ+HeXpLB/q5FqaJUpWNTzgA+d4mokWxpmWYDbsmSNlHDrkBC2mkw2ojfPdHCsxGguz54kz+yr52+IWrnLkFBEb1h/bBL3n+err0e6VvRGMhPbAk8hMbq9kJb2t2fhq/RE/NMTrW2RiXA53so9kESOdM7DQmCa+JmNPIjwVki9BgeZaGungvoxUrvhPuCJGG8fCFysI9Er5Dz+wDDgb7q8knBZ1GREPVW8/oMVwC8d8XUb5o4oq2fgAnC/QMo8piK+M8BZZscvbwLcltiZcTcTvxwkwAGE2x0dEiHCqwk31Vyuf/ck+NC2jYg/24vUOQgTYRPhMuUzxDecqPj3MnU7WKJF/hlK2H9bgUQStEMm4eYl8f6PkaJJDpLQ4+drJHJiCNK4c4leHwsD3P91SEuzBcAYpNzoLXZaTYTL1Rp2wX0f8Q/XR1nBEi0KywruqhakRzfEZ4v+7kjrSbQKxH9chRRQfwOZ3KsgevZbE1LDeSBwdcACvAIJQVsHfOCKD7rCyU0KtBGFCjsEBSHEG5CMup10qFqt1spU4CUT4IIR4LOQ0EB/nPcCpAwoSFLOClr7+Ecj7/0NcIz+bw6x049/inRi2Z5wokcQzEcKs/9Q+MeBkAmwWcKGMArp9QcSkvQLrc9rFIYAb0U4dndr4HgknMtV90Et4Wp5y/Vceq4Jf7v7Lgk+qhpJfx6s2/4moK8wC/jUIiBMhI3Y9PctdzQBLjjWqehWIpEK/kzEOuAkpNh6jboo9kBaR4EkTPTQ++3+OJ+xDXCRfs59SFGgmZk/P/jUERE2TISNONyMpJh21RvRKCyqkeQLryGmG+FuGOn7u4mWzTaXAX9P4jP2U4saJKkjUwFuBt53slunwjARLhWcxcCRdhwKygXRAXgQGAT8F0kD/pH+VAAv6ar+NORXkaao0TLcHKSSnoNMzkW6BrxO3SGkfoSXHp1ORMw64B3H+sKZCBtGEXMsEj8LcCIt438394nwZP3dgKQkr42xvTFIp22AO2ld1+EVpGxkSIX/dhXsK2mZHk0MgT8Wiai4C7jLyU3nDcNE2DCyhieIFUhBnMkqck20zIZrRtKKE9E7xrKf730PAC+RYqckRHg74DjPleFY3K+JsGEUuSuiMzLJdjYS//sWErHwsroc/GGD/dVF8RGtky/8jAf6qtX6SIJdeBGZ3HOAF5LYZX+fwKV2BosHqx1hGK0FeCCwJS3j6Dure6Aa8ff+XFalQof/HdUiPoX0fLixjCQHKbKeDIMQP/KtTnChbYZZwoaRM/GtQSJUukd5uV4FGKR2RJWKYwXQRv/vXw6CjantPrc6Ej5nFBGVAW7rx4TjIg2j2AS4I5Iq3DHGKit0yF+JREp42W4hJAa3Hsmm+yALu1eNxCBXxXE1fGfWr7kjzB1hlJL7oZA4l3Do2rnAF1HWec0xX3BRYrUjjHIW32oXdkBq+BbyvdDfd79GK6G52AS4+EW4gdj1UXsjs8KGUWruh70IV0DzoiEKkbuROhTvAK9Hef1LO6PFSxWS034+EhIzhpbtsUciQelfIMHqf7JDZpSAAPdDCvJ4RsgJiCttMXAOqfVzywXv6k80VhG9y7NRRJbwkUgQ+iRgn4jXn0QqR80gnNNuGMUqvlWulIccEeF+2Ft/d9HXiokvrTJa8YtwD6TACFGE1kXqmh6H5LIbRrEKcAd1P0TLVHtJfy+mZeGdoPCage5PsCFs60ncSskoAnfEXMKFQiJzzXsixacvU2v5QTtkRhEKcF+1cGOFZN6LFNQ5SH/GkVyCRKKW8g5SEW8n3/92Bi4P6Kt9bS3qS0OExyG+3p7AY0gO+lC9MA/VC20lUdIxGxsbe+pTnhtvvLH+9NNP3xSgurp6w0UXXTTPDq+Rb/eDim+fJFY/Qi1lED/rxCTW/ynSE+4ioidJ9IoQYJLcl2Ro1s82SkCElyKxh7U+S9gbkv0HKVhdQfRi1D/xhnd77rnn5kgLFzZs2PA9cJUdXiOPAtweCT9riHhpEBJv+1XE/0Mt356QsWqgbIYUY38jyjoLkQSKTZFkjyV6H+2G1ArOJL15jiOFhIwix5I1jHJyPxyINNh0gWsJ+4JBfMZHq1g+QuKU4d8iRdhX6nIsQa1BojG+VhfHNfogWA2cSnr1fl3gJUcsdqMELGHDKBf3w1Cf8bFlhAgvRwr0JMu/kOghB+m40R64gXDTT48m4HNdrgQG6HID4gJMR4QXmACXDpYxZ5SKALdFSj/GEuCu6oJYjLT7eTLzj2QBcKm6GzomMRJsBqaplf0p6bcvsuQME2HDKDj3w5+QWPer1SqNFOAbgV+qW+AUgpnUGoDMmXjMiLNuN6SbxnAdgW6W5mcucwovmSQQXnzxxSrXdZ1yu35NhI1iFt9KV9rPb6tWaAOSYr97xKqb+sRyyyQ2vQMSGVSfYL1pwHRkUu9FZCLbT3tksvosJBHK3+7+e9JLsviiFM/lE088cfCqVasWT5o0ae6ECROGmQgbRvG4Hzb1CSI61P8sYvWPkWSj1UjoZTQ6IZE+2wAXA78GzvS9Xqnr+GkCLkAm8jZT4fZzDJKgsT+S9ORFM3ygn5Eqa4D5pXg+Q6HQafrQ6u04zonldC3bxJxRrO6HyOiH25A+bN9HEaoNiO82FlsgCRQ1SOsivwsBtYivRj73CeBW3zp9kFraqJvjBcIxw/7Y4fn6eidah8cly1elmqLsuu7TjuOMATaGQqHnzBI2jMIU3wpXCu9sS+vwMxeJb0/HUtyGcAW1zsDzyMTZ7T6R7qvL+0e8dxnh/m6LaJl1+jDa+ViXl2YgwBuQ7NWSZMyYMTe5rrvVxo0bNx87duyz5XRdW5ywUSwC3ID4attnYfN91RJuh8QPvxrxeoP+vyfwDBKK5qe3WubvItEX2eBzJ/7En2EibCJsBCq6w5EkiHeRtkFbE+7xlq1RYRWxs9CqkYm1fPhkQ8DzTuvaLkYJYD5ho1B5AEmuOBnxo36aJSNkBOK7/Zz4acAbyN+k2FwTYBNhwyjmUVosjkXKtLrAFcDb+v/uSPGdj2lZUyJffGWXg4mwYeQMV4pJXaUC+XmWrGCQGg6e4A9SEe4L/BOZqHtRl/MZkbDAkfoUhomwYeREgOuAXZEJrhuy/HHjkPjdVciEmyfMXqTE3kgkRmMca7QzUoWwHVJP4nOzgo1UsBA1o5AE2EvAaMjRR04HfgWcjZSdBCnAMw2p8+AgNSEuRupsR+MQJGyuH+LeCJIVTmblLg0TYcNISnx3dCUGdibSBiifrAH+ANzsc0N0Bk6Psf7XMZaD4Au7OswdYRi54BzCyRCHqADmm2eQ2sJ/VIs4Vkfj15HOGm2Bt2Ks0w4YiMT5rk/y89chbe4NE2HDyKoV3NPnCvBcBIXCW8CFSF2Il+Ks93Gc1+oQX3EXxF98TpKf/aVTGJEZhomwUcIC3BdJGX5Gh/LNFF6t3E9JPjpjNHAU0o/xn4iI9iRcPW0Q0m05UczvRko4RdkwEc4r7733XvWCBQt2dRxn+siRIxeWsQAPAIYRjgeemYWPqQfWkpsQMwc4CYms2Bt4WsX7a+A1pOHnBJJLupjtJNft2SgBbGIux8ybN++xUCj0cnNz87Tx48d3LyPR7e3CfS7c5UoNiK3IbkLGz5Csu+toWXg9CCqRspXHEK457BIOJ1sDfOf7/1VqId+TxLZDwCy7U8wSNrLHPvq7a01NzQigXMr2NSLduUGK8NyZ5c8bqb/7I4XcPwhw2wci9YYh3LUD4C9IKNvnSMW0tJ7TjljvhlnCRpa4AhlqPl9RUfFqGX1vf3WxFTn4PK8F/SKCT6BoiHB5eKxFoiW+z2Db1j+uzLAqavkYmruu4ziOW1bfWSakGpHkh6eQSbhsX9t9kciLRH7YSqRq2xySK0XZBvH/tlUXQ1C+/YVO7DC3guTxxx/ftrKysm7UqFFv2J1tImwUvhDvS+6y4VLhQmAXtWTPSEJUa4hfcS1d3nYys6JzyoQJE452HOcB1ZHfjR49+jq7ys0dYRSuALfLsQBvCfwXuAOJxEi0LsgEXqJ1f4fUnPhzwPfPymISYGVXnyG3u13lJsJGYdMrx583EuiATJwdmGDdB9QK/kh/YtEOseYBdkQ6agRF0fmCHce5FckC/MZ13WvsEk8Pi44wckXPLBgQbWjZTNPPx8CPxAjnkwTbekJ/ErFKRWcI8A2wIKDvsh74tthO6OjRo2f4RhGGibBRwK6IerVK4w1rf67C9g8VpXiZZZ2Q2NvuwN3A+CjrPIMUwAmyQaaL+I/7674GlVAxy1KUzR1hGNkS4BqkLXw8TlRLeSdgLyTl9yH9fzS2BXog/siD42z3K4JP/92o4h7UxFwzMNuuFBPhBqQeajT6Ea5wZRipCPBhSEv4t5Di6TEtQZ8gVSMVxwBGxVj/E8LdJhKFRnVCOnTsUqCHabaTnUgLo4jcET2A85FWLmOA632vn67WRC/EornFDpmRAr9GIg7q1MK9L8Z616qYzkEyzX6MFL15O8b63wOnqMAm8qWej9SoALiAwqrS5mIpyibCwJHAZOB5xLfmF+E7kYmPPsAldriMFHlC3QUbiB91sAEpcuPxGzUO4g3T1xB7Us5P+xjLEA6vylfizDwnue9glIEl7FkctVEudJC2LVfb4TJS5CYk9KuZ5DLRPNYS7lLRESn4M530Igiu1+t3FtK6yGMYMsm2Efgr+eliYf3jDKqAuUgsJUSfjT4OeJPgUjONEseVIuhXI3G1LyEdIpw0LE4H+JuOxNYi/eCWp7iNGUi6dCQH6f4B7J8HEV7spF/kxygxER4H/AmZnX4MqQI1FLgXSUM+FfgQmbw7xf/mK6+8st3atWsrAa6++urqc889tyNAXV1d8wUXXGBtusuXvyGTcgD76e/XkLCyVK9PL764DvETLw9oH6ciccQgRdjTpTMyXzI/yfVrkDmWj+wyMbyLfCnSsrvWZwlP1d8P6U/0MePatb9HIydGjx49wrvJ1q5dOx/pUGuUJ9FGTbsjkQ+pxNZuAG5HOla8HzF8bwscr9byfYSjJZLlOaSQfDPpJ0qMUCu7Ut0eicqS1iCTkJsAh/seUEYZYwV8jGy4IxpUPA9VC7YaeAW4AXFVzCLzCalf+Kztx4H/5OGr/gw4QpffBi5PsP4ApMi8R3sn9YeHUYKWsGEETbNap/erldgJWAJcg8QAzwPOJLP42KYYy343wUlI9t1dSYj+cGAwEqq5JMl9eAU4QEeRzySx/jdIOvVw4C4TYMNE2MgW/QgnAjUjhdW9tu8gPtFucdwAuwInIEVt/oVEMETysLorHODRKK//lHAXk2XA/+Ls76bAZbqtPYDfJ/k9v0Ky+ipIzs3SDFwEfONk5oc2SghLWzaCdkVUED37ciXwpArRq4R7sEXjZMRvujdSrSwa64EHkQpo66O8vibGcjQ6EnbNddbf3ZEJ6kQuu2ZSryFhscGGWcJGYKLbDmnZ1AH4o4pWbYzVb0ayLv2hag6SUtxWXQEbkbDJ7ipw36W5a/9VC3gdkowUj4+RKKEh+nsAEmJXjSSc3BrwYXPsyjFMhI2gOAvpRuGJy22JdfuHdQ9GEjF29Lkx/g1cCeyGJGykW9xmnbosknyWtOiEfJAKMEhX6KAxETbMHWEEhn9yqUmH9smwF3BahLuhu/7eFil/ma8au28i/ui1ahmbCBtmCRsFy836MO+AFOE5DCmwE63FfDUyWbYmwmWxTF0Q9yKTcn/Q/2+GlLXMNStIfnLOzyAkdG4RkrK9zkTYMBE2soojk1L/1G7KTxBu/3Mp8F7E6j8lHNt7o1qZ3ZFwtnn6/xG+9bsU2eE4mXDFts+I3a3DRNgwETYCpx9SDCrSteCnW4TA3hNlnefUAu5KsAkYVYjvtzmLx2BhjGUTYcNE2MgevrC0O9Ua/B54Icqq9yFREKuBiep6OBWpu3AZ4l9eR8ussiDYGonP3YikGc/M0qG4SS3ghbSs2BaJzcUYJsJGoPRBfLwzkcLpsZiLtIr3OAbJpuuEJEk8maX9OwBxl4D4pLMlwutIrmGoWcKGPZGNQOmf5vs8MWwmu6Uk3/O5It4vgONlImyYJWxk7IJoQIo2DVc3wxtpbOYWpP/cQoJvyOnnJaSu8EYkcsFE2DARNoqew4GRuvzTNEU4RPRQNj+1iKtiHjAtg/2dX0DHzkaghomwkTEfq2VZBbyLFFvvkIXP+a2KsIukRU+zQ2+YCBtljwNTXNgG2AKYhMQLd9e/gxTj3r4hfO8MRLiS7IanmSVs2MVg5FyIP3VgvANNDrgOLHCkxu47SBZcENyJFIF/Ham+liqdkay+BwknkhTAoTMMs4SN7CnMAmCBK8kZQ0i+nkQ0pqhLIl12RULoAMYQPX7ZRNgwETZKUowXAgsDEuN0+RQpLFSDNKw1S9gwETbKWoy3QJIzcsXXSJfwLuS+rb2JsGEibBSkGHdWMe6ao49eqj8FdCgMQ7CJOSMfCrTEkZq9rxNOnuiMpBR3tfvOMEvYMHIkxsCbLvQEPkKqsC0DfkXsWryGYU9kwwiYZsJlMDvqcN21+84wS9gwcmMRL3SlutrPgfsdeMiVaIau+tMdqCutr2wYwV8MDwFH2yE1soErPee6qcXcrcityRUOvGxn1TBL2Cgma2EN0nl5tispyJ3VSu5GOE16M+BQpGLaM2YJGybChpEdBWtGw92A6a5UWuuCpCf3AfZH4oNnFu5XMAzBJgiMUhDl9Q58R8uIilkq0iETYaMYLOEGHdrNjrJOJVK9ao4dLqPAORI4A3jdgfEAPtdFD/2pN+PHKLQncg/gfOBFYABwve/1DsC5SJrpGQm2ZRNzRsHjm+Dz/MnVediNdQ48a2fD8J7IRwKTkZqw+0S8vhy41Q6TUUJWxxoHZjvSa+5ppPTmdCRxxM3dbhhG2B3RA3hb/661Q2KUkSC7amgsB77IYWyyibDRQoTnEs7Xt1RRo5xFuQmZ4PsOWrkuuhNcNJGJsNFChMcBf0Ly9x8DtgOGAvfqOocDg4Etddj2A42NjRVDhw51AI477jjn/vvvrwSYNm2a29jYGLLDaxS764JwbHIFEgYXGZtsImxkLMJLkcm3Wp8lPNW3zvW0nKzzc8W0adP6AZxwwgk7T5s2zRPu74Bz7PAaJSTIIaLHJndFXHptUticRUcYWXkiW3SEUba40E7FuKuKczyhdR2ZCDeMwsmYe/rppxuampquBfo0Nzefd9hhh1lrc6OYrJmVyM8XcdKqPfq5MMJpOeI0TITzS1NT08lIGxoqKyurgIPs9BhFKsjR0qq94kNHAn8AHBdOduAOO2LlTSH5phbEWDaMYhfl9Q7M1djkjoTdgPvY0TEKRoRHjRr1EHC867rnhUKh0+3UGCXK7cC3amjcZIfDCJKH7BAYhmEUqSVsGIZhImwYhmGYCBuGYZgIG4ZhGCbChmEYJsKGYRiGibBhGIaJsGEYhmEibBiGYSJsGIZRrlxODxoz11ATYcMwjFRpZBwbmA+8zK2Zdew2ETYMw0hNgDsiJUkB9mABQzLZXJUdUSMeEydO7Oq67q8rKiq+HDVq1AN2RAwTYZbRyNNIzfMP6chMs4SNbHKH4ziXua57/6RJkw61w2EYAIwEtqAXO3MW680SNrJJ1xjLRgEwYcKEUUCfdevW3XP00UevtSOSM2s4BJlZwCbCRnJDpYqK34RCocuAL9asWXO/HZHCYeLEiYcBjwLU1dXtDPzCjkrxYSJsxOXQQw/9CBhlR6Ig6edb7m+Hw0TYMIwcUllZeWdzc/OujuP0Bc6zI2JYeyPDMIwUsegIwzAME2HDMAwTYcMwDMNE2DAMw0TYMAzDKBARbqBlPKJhGIYREInihHsA5wMvAmOA64Pegccee2yTysrKu4Fa4OTRo0fPsNNiGIZZwsKRwGRgErBPVp4CVVW/123v5rruRXZKDMMwEW5pCS/T5dp4K/br12+ndHYgFAp95i07jvNZJl+msbHx3EsvvXTrfB7QxsbGTRsbGy/L94ltbGz8b7734ZJLLjnt0ksv3TXPx6FrY2PjtXY+4JJLLjmpsbFx/3zuwzXXXFPX2Nh4q50PaGxsPKqxsXFMIhGeS7hy1roE61amsyNjxoy5FRjjOM7R77///hWZfCnHcarJ/2RjBVBTAA/Y2nzvgOu61aFQqDKf+1BTU+PY+QgPPB3Hyev5WLFihVMgx6IgzgdQlcgnPA74E9ATeCxbezJ69OiJNigxDKMcSSTCS4Fz9amxLoo5/Q80cmL58uXtlixZMgdg9erVa8eNG/dRrr/M22+/PWTGjBn7EHah5Jybb765YbfddtuMPFe1mjVr1q7kuZ7Ha6+9tvncuXMPAxblax+uuuqq2oMOOmhYvo9FIZyPF154YcCyZcvWAifnax/+9re/VR511FHb2fmAp59+uu9BBx10dTJV1NxYrojGxsbzfMsPXXvttUfn80vtvPPOf9xll12euvjiiz/M1z6cdtpp/YHfPProo+fn81gMGDDgISCv52OPPfY4C/igsbHxtXztw/nnn9+tqampETi93M/HvvvuezIwd8qUKU/nax/++Mc/1gM3AT8r9/Nx0EEHHQuWrGEYhpFXAhPhlStXfkJ+kjr8ySRfhEKh5b7XKoFNcrwPaxzHme57rR/QNw/H4gPf/3sDg3O9D47jzAKWRLzej+xPlP2wD01NTesdx/lE/1+l56Iv0D6P56MKGAJ0yPE+fOO67kLf/71j0TtX+9C5c+dmx3Gm+F7rAmyR5/ujCzAg18K1cePG+R9//HEoqO31AP4P6cBwZg6/R7zP7QD8Fbghj/twOnAI0nbm1Dztx0jgt8ChwGV5PBYgkTafk11/ebx92Aq4Bpnn2D1Px6IeaAR21OszX+fjFD0O5wJP5GkfhqprYm/gX3k6FnsAfwYOBs7OoXYFrk+/AfbT5fE5/CKJPrdvDkQ43j7U6+8+wL/ztB+O/t4M+Eeez8f5ej7652kftgJO0tcr83QsTkR8kfsBdXk8H3W+e+S0PO3D3khHkLbAA3k6FpcDI3T5hRwbw32BG4JyRySd1JGFp1s+PjfZfVijv48Frs7TfrjqkjkOeC+Px+Igvcg35nEf5gHT9PcNedqPgfowmA7ckcdj4XVmPh64L0/7MAPoBZwBPJOnY/EiMFpHRgPzISBBiXAqSR1Bkq/PTWUfjgPeBBbmaT96AnPUFXFSHo/FIBWfocBRvlFCLvehWh9E08j+XEGs/VgIvA98B7TL87VZo4K0Ik/78EvgYeDvek3k41g8B9wJfKvLOSeoIdks4Of6VJsOTM3R/vs/tw7oFvHZJwG76Y23KAf7MB2ZdNlP9+NoHQZ1QXyzE/NwLI5B/I/exNzzeTof7+rNPlJF8PU8nI9D1CIfDHyvD8dcH4uv1fps8Fli+bo/jtbPX5ijfYg8HxuBw9Ud0R54Kg/HYj9kvqQz8Jqen1xxErDb/wPAd8L7GgjWBAAAACV0RVh0ZGF0ZTpjcmVhdGUAMjAxNy0wOC0xMFQwMDo0NDoxNCswMDowMBBwQDwAAAAldEVYdGRhdGU6bW9kaWZ5ADIwMTctMDgtMTBUMDA6NDQ6MTQrMDA6MDBhLfiAAAAAAElFTkSuQmCC" } }, "cell_type": "markdown", "metadata": {}, "source": [ "