{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# An example of polygon plotting with folium \n", "We are going to demonstrate polygon plotting in this example with the help of folium" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import geopandas as gpd\n", "import folium\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We make use of nybb dataset" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "scrolled": true }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
BoroCodeBoroNameShape_LengShape_Areageometry
05Staten Island330470.0103321.623820e+09(POLYGON ((970217.0223999023 145643.3322143555...
14Queens896344.0477633.045213e+09(POLYGON ((1029606.076599121 156073.8142089844...
23Brooklyn741080.5231661.937479e+09(POLYGON ((1021176.479003906 151374.7969970703...
31Manhattan359299.0964716.364715e+08(POLYGON ((981219.0557861328 188655.3157958984...
42Bronx464392.9918241.186925e+09(POLYGON ((1012821.805786133 229228.2645874023...
\n", "
" ], "text/plain": [ " BoroCode BoroName Shape_Leng Shape_Area \\\n", "0 5 Staten Island 330470.010332 1.623820e+09 \n", "1 4 Queens 896344.047763 3.045213e+09 \n", "2 3 Brooklyn 741080.523166 1.937479e+09 \n", "3 1 Manhattan 359299.096471 6.364715e+08 \n", "4 2 Bronx 464392.991824 1.186925e+09 \n", "\n", " geometry \n", "0 (POLYGON ((970217.0223999023 145643.3322143555... \n", "1 (POLYGON ((1029606.076599121 156073.8142089844... \n", "2 (POLYGON ((1021176.479003906 151374.7969970703... \n", "3 (POLYGON ((981219.0557861328 188655.3157958984... \n", "4 (POLYGON ((1012821.805786133 229228.2645874023... " ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "path = gpd.datasets.get_path('nybb')\n", "df = gpd.read_file(path)\n", "df.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot from the original dataset" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAFpCAYAAAB3UOSMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd8W9Xd+PHPkWTJe+8VZ++9EyABAmGPssJKOhmFFrr70D6llPJ7aKFAmQUKLZRVdgIBQggZhOzpDMc7cRzvPWWt8/tD144dy5bsyLZsn/frlRfSuedeHUG4X90zvkdIKVEURVGU9nQD3QBFURTF96jgoCiKonSigoOiKIrSiQoOiqIoSicqOCiKoiidqOCgKIqidKKCg6IoitKJ2+AghEgRQmwUQmQIIY4IIe7TymcIIXYIIQ4IIfYIIeZp5UII8bQQIkcIkS6EmNXuWquEENnan1XtymcLIQ5p5zwthBBaeaQQYr1Wf70QIsL7/woURVGUM3ny5GADfiGlnAgsAO4RQkwC/go8JKWcAfxBew9wKTBW+3MH8AI4b/TAg8B8YB7wYLub/Qta3dbzLtHKfwtskFKOBTZo7xVFUZQ+5jY4SCmLpZT7tNf1QAaQBEggVKsWBhRpr68GXpdOO4BwIUQCsBxYL6WsklJWA+uBS7RjoVLK7dK5XPt14Jp213pNe/1au3JFURSlDxl6UlkIkQbMBHYC9wPrhBCP4wwyi7RqScDJdqcVamXdlRe6KAeIk1IWgzNICSFi3bUxOjpapqWl9eRrKYqiDBt79+6tkFLGuKvncXAQQgQDHwD3SynrhBB/Bn4mpfxACHEj8AqwDBAuTpe9KPeYEOIOnN1SpKamsmfPnp6criiKMmwIIU54Us+j2UpCCD+cgeFNKeWHWvEqoPX1ezjHEcD5yz+l3enJOLucuitPdlEOUKp1O6H9s8xV+6SUL0kp50gp58TEuA2IiqIoihuezFYSOJ8KMqSUT7Q7VAQs0V5fAGRrr9cAK7VZSwuAWq1raB1wsRAiQhuIvhhYpx2rF0Is0D5rJbC63bVaZzWtaleuKIqi9CFPupUWA7cDh4QQB7SyB4AfAX8XQhgAM1q3DvAZcBmQAzQB3wOQUlYJIR4Gdmv1/iSlrNJe3w38GwgAPtf+ADwKvCuE+AFQANzQi++oKIqi9JAYavs5zJkzR6oxB0VRFNeEEHullHPc1VMrpBVFUZROVHBQFEVROlHBQVEURelEBQdFURSlExUcFEVRlE5UcFAURVE6UcFBURRF6UQFB0VRFKUTFRwUZRhptthZm16M2Wof6KYoPq5HKbsVRRncXtt+nEc/P0ZUkJHbFozg9oUjiA42DXSzFB+kgoOiDCO5ZQ0AVDZa+PuGbF7YnMuyibFcPCmeS6fGYzLoB7iFiq9QwUFRhpETlU0d3ltsDj47VMJnh0p4cI0fl09L4JZ5qUxJChugFiq+QgUHRRkmpJRkFNd1eby22cpbOwt4a2cBU5JCuXFOCldNTyQ80OiVz39ifRajooO4ZmaS+8rKgFPBQVGGiSNFddS32Dyqe/hUHYdPHeHPn2Zw0eQ4VsxNYfHoaHQ6Vxs3uldSa2ZLVjl3njeqV+cr/U8FB0UZJrbnVvb4HIvdwdr0YtamFzM6JogVc1O5emYisSH+Hl9DSklBVRMWm4PXt59gf0E1dy8dzczUiB63R+k/aj8HRRkmbv3nDr7N6XmAOJNeJ1g6LoYb56Zw4YRYDPruZ8TvOV7Fba/sxGx1dChfMCqSp2+e2aNAo5w9T/dzUE8OijIM5JU3sCu/yn1FD9gdkg3HythwrIzIICPXzEjixrnJTIgPdVn/SFFdp8AAsCOvise+yCQswI+fLhtLqL+fV9qneIcKDooyDPx3z0msdu/3ElQ1Wnj123xe/TafyYmhrFw4giunJxJodN5aDpys4d09J7s8/729hQBcODGOhaOjvN4+pffUCmlFGQZOVTf3+WccKarjNx8cYs6fv+Ln7x5gS1Y5a9OLOFLU9QypVttzK/q8fUrPqCcHRRnipJTsPu6dLiVPNFnsfLjvFB/uO+XxOf/cms/tC9OICVGrtX2FenJQlCEuq7SB0rqWgW5Gt5osdlpsKt+TL1HBQVGGuG9zfL/LRggIUQPSPkUFB0UZ4jZllQ90E9ySEjYPgnYOJyo4KMoQ1thiGxRPDgBv7DjB4VO1DLW1V4OVGpBWlCFsR14ldsfguNnuyq/iime2Mio6iFvmp3LNzCSVTnwAqScHRRnCBmNXTV5FI39em8GiR7/m3rf2sS23AodDUlZvZs/xKhyDJNgNdurJQVGGKCklGzLKBroZvWaxOfg0vZhP04tJjghg8eho1hwsIjrEyMoFaVw2LYGk8ICBbuaQpXIrKcoQlVNWz7Intgx0M/qMEDBnRATLJ8ezfHI8KZGBA92kQUHlVlKUYW4wPzV4QkrYfbya3cer+fPaDKYmhXHtzCQumhSnAoUXuB1zEEKkCCE2CiEyhBBHhBD3tTv2EyFEplb+13bl/yOEyNGOLW9XfolWliOE+G278pFCiJ1CiGwhxH+FEEat3KS9z9GOp3nriyvKULfmYNFAN6FfHTpVy58+Pcq5f93I8ie38I/NuZyobBzoZg1abruVhBAJQIKUcp8QIgTYC1wDxAG/Ay6XUrYIIWKllGVCiEnA28A8IBH4ChinXS4LuAgoBHYDN0spjwoh3gU+lFK+I4T4B3BQSvmCEOLHwDQp5V1CiBXAtVLKm7prr+pWUpSh36XUE9OSw7h0SgJXzUhUYxR43q3k9slBSlkspdynva4HMoAk4G7gUSlli3as9Rn2auAdKWWLlDIfyMEZKOYBOVLKPCmlBXgHuFoIIYALgPe181/DGXxar/Wa9vp94EKtvqIo3fh4//B6auhOemEtf/niGLf9c6daQ9EDPZrKqnXrzAR24nwaOFfr7tkshJirVUsC2ufoLdTKuiqPAmqklLYzyjtcSzteq9VXFKULUkpWH/Q86d1wcffS0QghOFZSh83uoLTOPNBN8mkeBwchRDDwAXC/lLIO52B2BLAA+BXwrvar3tUve9mLctwca9+2O4QQe4QQe8rLB9+8bkXxpq05FZys6vsU3YNJamQg189KZvfxKh5Zm4HNITlUWDvQzfJpHgUHIYQfzsDwppTyQ624EOc4gZRS7gIcQLRWntLu9GSgqJvyCiBcCGE4o5z252jHw4BOuYellC9JKedIKefExMR48pUUZcgKMRm45/zRLo/NGxnJvJGR/dyigTcpIZS1h4q58cXt7Myr4vZXdvLYukxe3358oJvmszyZrSSAV4AMKeUT7Q59jHOsACHEOMCI80a/BlihzTQaCYwFduEcgB6rzUwyAiuANdLZCbgRuF677ipgtfZ6jfYe7fjXUnUaKopLzRY7b+w4wVcZZVw+NZGfXDCm7ViIycCs1HB25VexK79q2AWI62Yn09hiQ0qw2B3sPl5NZmk9O/OcvzUHS4qR/uTJOofFwO3AISHEAa3sAeBV4FUhxGHAAqzSbtxHtNlHRwEbcI+U0g4ghLgXWAfogVellEe06/0GeEcI8WdgP85ghPbP/wghcnA+Maw4q2+rKEOU2Wrn2ue/5VhJfdv7ey4YQ3SwiYziOhpabHyaXtxWf9+J/tv8Z6BNTw7joklxfHG4pEO5Ua/je4vTVGDogtvgIKXciuu+f4DbujjnEeARF+WfAZ+5KM/DOZvpzHIzcIO7NirKcPfunpNtgQEgLMCPf317nJmp4QSZDFQ3WrhpbgprDhTx3t5CUiODyKsY+msALp+WwM+WOWfS78rvGBDvXjqamakRmK12gkynb4UWm4OaJguxof792lZfo1ZIK8og12Kz8+LmvA5lr3ybzzUzkvj9R4c5VeMcnBYCfrV8PCsXpJJRUj+kg0OAn55mq527l4wmLMCPZzZk89/dBR3qhAf68c7uApaMi+kQHIwG3bAPDKCCg6IMep8fKmkLAK1qmqz8e9vxDmVSwvMbc/nk3sWc/7fN/djC/mEy6GixOQD46YVjySypY1JCKFc+u5UjRXWd6j/0yVFCTAaCjAaSIzqm27DaHfzvx4f5xcXjh+2+1io4KMogtynT8xxKyREBHK9s6sPWDBwhnEHh1a351DRbePyG6Xx+uMRlYFg4Kgq7lIT6G7hgYmyn4356HX++ZgrN1uG7r7UKDooyyBXVeL6Ya8n4GPafrOnD1gwcs9XB+qOl/OW6aYyLC0YIwZNfZXWooxPw4u1zuGhSHFJKKhoshHaxd7VBryNEP3y3vBm+31xRhgCb3UFhtWdPAn56wc1zU/nqaGkft2rgZBTX8dKWXKoaLXx2qJicsoYOx+emRXLRpDgAhBAedRlll9ZjHoZPEOrJQVEGsc8Ol1BU69mTw6VTEiitM3O0uHM3y1BysLCWhz456nKK6nnjerZItrbJyti4EG81bVBRwUFRBimb3cFTZ3SbdGXZxDjuXzaWVf/a1cet6l9GvQ6L3dGpvKsAOCo6qEfXDwt03eU0HKhuJUUZpN7cWUBeuWfTUW+dn8Jnh4qHXM6lv14/jZ9eOJYQk2e/c786YwOk7NJ6yupVAj5XVHBQlEGotsnK4+syPar7u8smUtts4/EvPXvKGCzOHRvN8snx/PyicTx/26y28ruXus4rBXD8jM1/xsaFEBXUcdyhoqGFMpWxVQUHRRmMXv02n/oWW7d14kJNTEoI5YY5yTzyWUY/taz/fJNdwbInNlPbZOXcsTHcvmAEd5w3isWjo13Wn54Szgu3ng4inx0q5tWt+eh1HRNARAebiAkx8Wl60bDeSU6NOSjKIFNc28zzm3K6rWMy6Hjgsokkhvnzq/fTKa9v6afW9a9TNc28t/ckPzx3FL++ZDzNVjsRgUYCjXqaLB1nGP3y4nHEhJgoqGwi0KQnKTyAkV2MQQgheH37CbJK6/nk3nOG5Z7U6slBUQaZJ9dnYbV3nyzOISXRwSY2HCtn/RCeugpQrM3WCvH3IzbEHz+9jkQX24EGGvUcKapj5as72ZpdwfSUcCYmhHZ5XZNB53Kl+XChgoOiDCKZJfW8u6fQbb3lk+MJMup5cUtuP7RqYNU0WTuVTWp3058QH8K6+89jfHwo23IrOF7Z1JZuxOGQrDtSwp7jnbPUtm4O8MrWfP6z40TfNN6HqeCgKIPIq1vzPaq3Ym4qD35ylKG8+0lEoB/j4oIJP2O6aUOLjbWHnOnJQ/wNvLxyDuPjQwg2GbhggjNVxtHiOurNVnQ6wdLxMYyN7biWocVmx6pNkV00OopxscH98I18ixpzUJRB4lhJHR/td7839KVT4skoruPgEE2T0eqJG2dw/oTOeZGCTQbGxgZzrKSeJoudUzXNbWMGYQFGEsL8+TangjqzjRB/P0wGPSaDvsM1jHodUcFGooNNPH7DdJfdVEOdenJQFB8npWT1gVPc/soulwu+2tMJuGvJaJ75OrufWjdw0rpZ0PbwNVMA5w5vB7QgaXdIAox6VsxNpabJyi0v7+Afm3NptnROjSEl1JttNFtsRAYZ++YL+DgVHBTFxz31VTb3vXPAoxlH185MJv1ULXXm7qe5DgWPrTvW5bG5aZGcOzaaxDB/LpkcD4BeJwg2GZic6ByPOFHZxLu7T+Lv1/k2qNMJfnnxeBotdpdjGsOBCg6K4uP89F1txNiRUa/jjvNG8tzX3U9zHSrWHy2ltpsb95JxMfzkwrEdnjCklLy/9/SAfl5FI7e9spOims4rx0P8DSSFBxAfNjw3/lFjDori4/Q6z37D3TI/lf/sOEHJMFndO2dEJEEmPVJKhOgcQH9wzkjOzL1X32LjiyMd95I+UFBDZYOl07jCqZpmkiKG31hDK/XkoCg+bkSU+wVYgUY9V0xP4J1dJ/uhRb4ht7yBZzfmUFjdjMNFBlYhRKfVz0FGA3PTIjqUTU8JZ1x859lI54yJZnxcyJBOcd4dFRwUxceN8yBl9G0LRvDathPYXNwkh6IAPz3+fnre2llAWb0Znc6zrje9TvDItVNp/6CxbGIcVrvsFGCEEPzxqsksHB3lzaYPGio4KIqPy69wn99n2cRY1p3RXTKUNVvtzB8ZyUNXTe60ytnh6Hyjby8uxJ/x7QKu0aDDbLXjomcKvU4Q5GHG16FGBQdF8XHu8iilRAZQ22TFYut+mutQ8/6+QhaPjSbQePrmvSWrnEkPfsHNL+/ocve2sEA//veKSW3vn1ifRYPZ5nLcYjgbniFRUQaJrNJ69hd0v5jt9gUj+NCDxXFDjZSwv6CGJe12dztvXAwZf7oEq11iNHT923f2iAjCAvyobbZS1Wghu6yh23UTw5F6clAUH/bB3u7zKJkMOpZPimdTZnk/tahvTIgPIbGHU0ZD/Q0up/kKIboNDAD+fnq2/fYCwgKcqTfsjuH11OUJFRwUxYdtzur+pn/NjCQ+PVRMcxddKINFRYOF0T3MX3TDnBQWdbF3gycCjXqWjo9hXFwwS8Z1TsMx3KluJUXxUWarneyyhm7r3DQ3mR++vrefWtR36sxW/P307iu2k1/R2OUaB08IIXjs+ulUN1kIMPbss4cD9eSgKD7qaHEd9m5m3SRHBGA06KlqtPRjq/qGxebAw9mobRpaejeILKVs2+PCaNARFzo8V0C7o4KDovioN9zsIRATYiLPg2mug0FSeAAGfc9uR6W9XAkuhOCiSXG9Onc4UcFBUXxQSa2ZtenF3dZxZg4d/EnhooNNLJ8c7/b7nulEZRM5brrdlN5zGxyEEClCiI1CiAwhxBEhxH1nHP+lEEIKIaK190II8bQQIkcIkS6EmNWu7iohRLb2Z1W78tlCiEPaOU8L7VlRCBEphFiv1V8vhOi47l1Rhqinvsqixc26hcrGFqKDTf3Uor5z3ewkXt9+vFfnPr8pBzmUdzQaQJ48OdiAX0gpJwILgHuEEJPAGTiAi4CCdvUvBcZqf+4AXtDqRgIPAvOBecCD7W72L2h1W8+7RCv/LbBBSjkW2KC9V5QhrbC6qUPm0K4U1ZiJGgJ7DdQ0Wnud9uPj/adIL6z1cosU8CA4SCmLpZT7tNf1QAaQpB1+Evg10P6/7NXA69JpBxAuhEgAlgPrpZRVUspqYD1wiXYsVEq5XTp/ArwOXNPuWq9pr19rV64oQ9IXh0t44KPDHt0s7Q7ZaQezwajF1vtpuA4JH+0/xa78qm5TZig916OprEKINGAmsFMIcRVwSkp58IwZA0lA+9SQhVpZd+WFLsoB4qSUxeAMUkIINRlZGbI+2l/Iz/57sEfneJjN26cdLa5j0egotuVW9vjcUH8Dv7t8In56HaV1ZgKMekL9/dyfqLjl8V8tIUQw8AFwP86upt8Bf3BV1UWZ7EW5x4QQdwgh9ggh9pSXD+6VosrwZLM7+PtXPd/a02aXBA7yOfpZpQ3MSu3dcKLNIdlfUMM/v8njoic287N3Dni5dcOXR8FBCOGHMzC8KaX8EBgNjAQOCiGOA8nAPiFEPM5f/intTk8GityUJ7soByjVup3Q/lnmqn1SypeklHOklHNiYmJcVVEUn5Zb3kh1L7ajrDfb2lJADGZ1ZisJvdhxrdlqZ+WrO/nz2gzqzDbyKxq7XRuieM6T2UoCeAXIkFI+ASClPCSljJVSpkkp03De4GdJKUuANcBKbdbSAqBW6xpaB1wshIjQBqIvBtZpx+qFEAu0z1oJrNY+fg3QOqtpVbtyRRlS9DrBxAT3+zacqd5sJcR/8Cc6+PpYWa/WHkgJZuvpWV15FY3klavprd7gyd+qxcDtwCEhROsz2wNSys+6qP8ZcBmQAzQB3wOQUlYJIR4Gdmv1/iSlrNJe3w38GwgAPtf+ADwKvCuE+AHOGVE3ePi9FGVQOVZSx468KvcVz9BstRPQw7QTvshsdaDzUsrsmubBv/bDF7gNDlLKrbgeF2hfJ63dawnc00W9V4FXXZTvAaa4KK8ELnTXRkUZzE5WNfGb99N7dW5FQ8ugn7Hk76fjN5eM5/Ap70xJfWRtBu/csQCTQaf2aDgLQ2Cug6IMXnaH5FfvH6TR0rvpnGarAz/D4L4Bnjc2hnPGRvPmzgL3lT1w4GQN5z++ifMf38TeE9VeueZwpIKDogygF7fk9qo7qdVAhgVvrc5euTCNEH8/DC72ZuipsbHBJIb5U1xr5nhlE9e9sI0/rD7shVYOPyo4KMoAKa9v4ekNPZ++2l54oN+Azc45f3wML6+cw/3Lxno8Y+rMzXnGxAazaHQUwSYDUUFnH2xyyxu4cW4Kwdq+z5dNjedXy8ef9XWHIxUcFGWAvLg5t8NMm96w2iVNveySOlvv7S3kzv/sobHFxssr53RbN9CoJybExMqFaVw1PZHWoYDlk+PQ6QSF1U2U17d49LkJYf7cPC+Vh6+Zwg/PGcn8kZGEajO2HBJOVTdz99LRTIgP4Y9XTSZELYrrlcE/B05RBqHC6iZe3959Sm5PRAYZqenF+ghvcUh4+Zt8Fo+J5vKpCaw95DqzapPFTpPFzitb84kONnH/heN4dmM21850JkPILW/EYu86UIb6G7h5Xio3zEkhOSKg08ZADodkS3Y5d/xnL2sPFfPAZRN5766FKjCcBRUcFGUAPLcxt9uboafiQk00WWxeaNHZWZtezAOXT6S22UphdRMFVU101dsVYNRx64JUgv0NjIl1ru2YkRJOiMlAfcvp7xIW4MdTK2awYGQULTY7YQF+Xc4+0ukES8fH8vvLJ/KH1Ud4YVMus1IjmJSogkNvqeCgKP2spNbMBx5kXfWEXid8Yirrp+nF/P6KSbzxw/kANFls7D1RTU5ZAw99crRD3QvGxxIdbOIH54xsKwsL8OP2hSN4flNuW9kfrpjE+eOd6dQ83cZz5cI0UiID+cvnx/jDamcCw/PGxXDXklEEGtXtrifUmIOi9LPHv8z0ylMDwLbcSq6YluCVa52NZqudD/edDniBRgPnjo3hqumJneruyKuiuLa5U/mKuakd3s8bGdmrtpw/Ppa1Pz2Xd+9cyKLRUTy9IZsL/7aZNQeLsDskDS0D/6Q1GKjgoCj96ODJGj7Y552nBoB/fpPPDbNTmJQQetbXCjTqmZ4SzrKJsSwd3/McZVuyOie9LK41c2ZPUGZpPe/sOtmpblJEQIfvsft476f46nUCnU7w60sm8OLts4kP8+enb+9n2RObeeWb/F5fdzhRz1mK0k+klDz86VG8uXFZQ4uNX39wkAevnMRNL+3o1TWSIwJYPDqa78xKYm5aJBJwSEm92cbL3+Tx390nqWq0uL3OzfNSO5WNjw/hT1dN5qP9pzhSVEeQycCIqECqmzpfT68TfPKTc/j6WBlbssq9Npi8fHI8F0+K4zcfpPPunkLyK1TuJU+IobbF3pw5c+SePXsGuhmK0snnh4q5+819fXLtJ26YDgJ+/m7P9oO4cEIsT66Y0e0eCGvTi/nJ2/u6HGAGmD8yknfuWODT6SrsDkleeQOpUYE+MU4zUIQQe6WU3c89Rj05KEq/sNkdPP5lZp9df016Ec/dMotfLTdTb7aSXdbAhgyXGe7bjIsL5vnbZmEy6JFS8sXhEi6ZEt/pBn/5tASqGiezI7+KmSnhTEsO57Xtx1mb7py2qtcJ/njVZJ8ODOBs59i4nme+Ha5UcFCUfvDe3kJyyxv77Prbcpy7qEkp+cfmPIx6wazUcPYV1LisbzToeHnlnLZf0Kdqmrn7zX28cOssLp3aeYD79oVp3L4wre19ZJAf189OJresAaNBx0QvjHl4QkrJ+3sLOX9CrNfSdyiuqeCgKH3MbLX3ape3nrDYHdQ0W3lfmyJrsUsMeh1xoSZK6zqvPH7oqsmMiApqe1/XbOOG2cks83BPhTGxIYyJDeH88bH9tnezlJLnN+Xy2LpMxsQG88zNM/stKA1HaraSovSxN3acoKTO3OefY7bYOV7Z1PZ+V34VUsJFk+Lacg0BrJibwk1zUjqcOykxlMdumI6f/vQt4WRVE2/tLODB1Ye7XWin0519d1KzxU5tsxVXY6ANLTbK6s08sT6Lx9Y5u+Zyyhr4zvPbeG5jDttzK6k3qz0cvE09OShKH6o3W3luY06/fFZMiJGR0UHkV5zuviqrb2H90VLuPX8MoQEGdELw/cUjXd7QH1+XyfLJ8UxNDkNKyS3/3MHJKud6hKzSBv79/bleHchtaLERbDJgttr5z47j5JY18pfrp3WoU17fwv99nsGH+051Or/ZauexdZnavg0wb2QUc0dEMDkplLK6FsIC/IgMMjJ7RAQGffe/g81WO2/tLKCsvoU7zxtFRJDRa99zsFLBQVH60Etb8nq1N3RPnTMmmhB/P2wO14vr/vVtPtsfuLDbWUkrF45oG1Q+XtnUFhgAtudVcs1z2/j+4jQmxIcyOTG0V08MFpsDo0HHvW/t47NDxTx+w3Q2ZpazKbOMerONey8YQ0pkYFv96GAjJ6uaurkitNic33lLVrnLtRYT4kO4fnYyyyfHd7h2q2aLnXvf2seGY84B/B15lbz9owUer8oeqlS3kqL0kbI6M//shwVXyybG8fTNM9meV9nhht5eo8VOQWX3N9nYUH9iQpyDvOmFnQeyM4rr+NX76Vz57FYuf2YrOWU9Wy+wNr2Y77zwLZ+mF7E9txKHdE69/eRgEfVmZ7fVZ2ck7hNCMD05vEefc6ZjJfX8eW0G+wpcb/wTYNTz3cVpbe8PnKzhfS8uVBysVHBQlD7y5FfZNFv7Np32Ly8ex8srZxMZZORoUR1Gg+v/pZeMi2FKUpjH161zsw9zRnEd1/9jW4+29my02Dh8qo5739pPZReL6rbmVHR4X1TTzNfHup+S66mxsZ2nsdodEodD0nxG2vO3dxZ0KhtuVHBQlD6QV97Au3s6p4jwpl8tH8+9F4ylocXGI2uP0mJzMCYmuFO9+FB/nr55JuBMbf3V0VIKKptcDv62mpToPpDUNFm54pmt2qpv9zOWJieGYuoieLW6blZy2+uHPz3K0sc3kVfhnSnAedrK6GaLndI6M5UNLWSV1vPQJ0c6JQc8WlzHVxmlXvncwUqNOShKH/jb+qw+3aEtLtSZ1VRKyS/fO8i6I65vZHNGRPDGD+fj76cBqLtiAAAgAElEQVSntsnKY18e440dzr2aJyWEMn9UJDNTI1gyNoawwNPjETNTwvni/nO58pmtWO3df49XtuZz9YxEpnXT/SOl5NZ/7mwbH+jKpswygk0GloyPIb2wBoub+j3x3p5CJsSH8N1/7aaw2nX3W3tv7Szg8qkJXpmNNRip4KAoXnb4VG3b6uG+ctv8Ec5cROnFXQYGo17Hw9dMwd9Pj83u4O4397Itt7Lt+NHiOo4W1/Gvb4+TEObPWz9awMho59oHnU4wIT6UeSMj+Tan0uX12yuuNTMtufvjF0yIdTnrqL2PDxTx8YEiJieGkllS7/Zze2JzVjmbn+g8YH2mUH8DdWYbmaX1lDe0EBfq79V2DBYqOCiKF7XY7Pzli2N9+hlhAX6smJdKvdnGH9cc6bLezy8e17ZI7OVv8jsEhjMV15q5/oVtvP6DeUxu16X02PXTeXfPSQqrm9lfUO1ylfekhFCWjo9pm5raXrPFzt++zOTD/ac8St7X6khRncd1vc1sc/CbSybw3UVpw3rGkgoOiuIldodkW04l32RXuK/cS356waPfmUpMiIlHPz/W5Q13eko4Pzp3FAB7T1Tx1FdZbq9d2WjhoU+O8u6dC9vKEsMDuH/ZuLb3FQ0t7M6vorLRwpy0CFqsDmx2Bz95az/P3jKr0zVzyxt4c2dBnw/Me5PF5mDviWruWjJqoJsyoFRwUBQv0Ql45uu+S5MxOiaI526dxYR459NAZUPntBitlk+OQ68TZJXW8+CaI277+gEC/PT89tIJVDVaiOxiEVh0sKlD7iWHQ/KdF7ZxsLCGJosNo6HjeVOSwrhoUhxrDxX36RiMt32VUcraQ8VcMa3zZkXDhQoOiuIlXx8r6zLRnTc8et20tsAA8IcrJ7F0fCzphTV8tP8UZfUtBJsMXDolnu8vHklueQM3vridGjeL8OJCTVw7M5m4UBP/b20GN85NYcHIKJIiAtC7GIxtsth4+NOj/HjpGDZllXOwsIbvLx5JeKDrgPLUTTNYNimOJ9dnMTM1HJtdsuZg0dn9y+gHv//4MFOTwjrkoHJFSsnfvsxiQkLIkAomaj8HRfECKSXXPr+NAyf7Ljh88+vzSYkMpKKhheIaM1OTT48NVDa08MT6LG6dP4JJiaGcrGpixUs7OFXjflbO+eNj2Jh5eqDWTy9Iiwrigx8v6rSiut5s5VRNM1c9+y06AQLBxZPj+PuKmR5/jy8OF3PXG32zr4W3TUwI5eN7FnWbNqT99/nRuSP5zSUT3KbrGEhqPwdF6SeNLTbSC2v7NDCA88lk1aI0bHbJiOiOaSCigk08cu1UAErrzFz17FaP03acOVBttUtqmq0dppEePlXLj9/cR3Ftc6eprZ8dKmZ0TDB3Lx3dIXGfK2ar3aMuLl+RUVzH8xtzuXvpaPz9XAeIZ9vlznr5m3xyyhp45pZZnQbnBxvfDW+KMkgEmQx82A/pFloXZcWH+RPq70dts5VNmWWcqDw9g2htejHX/2Nbj/I5ubpZVzVayClraFvcVtVooaCqyeWaB6td8sT6LK59/tsuU1S0Mhl0zEyJIKCLG60v+nB/YZeBoazOzNEzZlZtzCzn2ue+7XF6EV+jgoOinKV/fZvPe3v7PjjsO1HdljrbZnfwvX/t4rv/2s2Ff9tMdmk9T2/I5p639nWZX6knwgP8+GjfqbbFYp7sm3D4VB3feX4bj607RovN9ewkIQSpUYH4+w2eW8/JquYunwr3nqh2uX1qdlkDVz+7ldUHTnm0etwXuf0vJIRIEUJsFEJkCCGOCCHu08ofE0IcE0KkCyE+EkKEtzvnf4QQOUKITCHE8nbll2hlOUKI37YrHymE2CmEyBZC/FcIYdTKTdr7HO14mje/vKKcLZvdwT825/bLZ0WHmNq6bT4+UNQ2+G1zSC56cgtPrHc/XdVTlY0WVh88xdFi569io16HwcOVws9tzOXyp7ey+3iVy+NSSmyDaOYSwLNfu067fu64GEK66D5qtNi5750D/HHNEWz2wdOV1sqT8G0DfiGlnAgsAO4RQkwC1gNTpJTTgCzgfwC0YyuAycAlwPNCCL0QQg88B1wKTAJu1uoC/AV4Uko5FqgGfqCV/wCollKOAZ7U6inKgHM4JDllDfx723GXO631hV9ePL4tOHyT7X6l79kyW53ptb84XMLL3+T16NycsgZu+Md2fvtBOrX9kLK8r319rJRDhZ2TDAabDG53z3tt+wlu/edOyuv75++Jt7gNDlLKYinlPu11PZABJEkpv5RStm4PtQNoXTx/NfCOlLJFSpkP5ADztD85Uso8KaUFeAe4WjgTyF8AvK+d/xpwTbtrvaa9fh+4UPj6LubKsJBX0cihUzU808UvSm+7YXYyV05PZPfxKhpabCwaHdUv/fbf+9du7npjL89uzOnVr/13dp9kyeMbeXuXM5+TlJIjRXVtKboHC4eEH7+1l7L6jjv6ma3OJH7u7Myv4spntvYoi+1A61HHn9atMxPYecah7wOfa6+TgPbpKAu1sq7Ko4CadoGmtbzDtbTjtVp9RRkwLTY7Y2KD+Sarglo3qa29Yd7ISB65dipv7Szghn9sRwA3zU3l8/vOxejDUyZb1TRZOwzOhvr7te0bMZicrGrmu6/u7rAl6XMbc7pNS9JeSZ2Z617Yxps7TwyKcQiP/2YJIYKBD4D7pZR17cp/h7Pr6c3WIheny16Ud3etM9t2hxBijxBiT3l53z9uK8PbtzkV/Hd3AZ/2cXI9cCaB+/uKGew+XsWDaw6zaHQUgUY9UkqOldRjGQR92VdNT+S+ZWOB0wPS05M931vClxwtruMX7x5su7lfMCG2R+e32Bz87qPD3PfOAcw+nlLEo+AghPDDGRjelFJ+2K58FXAFcKs8HQoLgfa7lycDRd2UVwDhQgjDGeUdrqUdDwM6jXJJKV+SUs6RUs6JiYnx5CspSq8U1zYTG+JPXkVjv9yYf3/FJCICjfxxzRFmpkTw0so5HCmq47oXtnHXG3v7/PO9YVJiKH66jreaqKDB9+TQ6sujpW37gputvfs7sOZgETe9uN2jRYoDxZPZSgJ4BciQUj7RrvwS4DfAVVLK9vsPrgFWaDONRgJjgV3AbmCsNjPJiHPQeo0WVDYC12vnrwJWt7vWKu319cDXcjA8jylD1smqZmwOSZ6L7KTeNntEBNfPSualLXmcrG7iyRUzWHOgiCue2dqnaTq8ISzAj3vPH8Ofr5nCneeN6pTd9MrpibQfPQw06rvM5+SLXtqSx9++zEQn6HXX3sHCWq5+9ts+XzzZW27TZwghzgG+AQ4BrWHyAeBpwAS0drjtkFLepZ3zO5zjEDac3VCfa+WXAU8BeuBVKeUjWvkonAPUkcB+4DYpZYsQwh/4D85xjipghZSy22kTKn2G0lcqGlrILm2gqKaZX3+Q3qeJ5AKNej6/71zsDsklT33DTy4Yw0WT47jqmW8HRVcSwH9+MI9zx3b9JF/R0MLr208wMd65fefzm3I5NIgGbMH532nR6Ci+yji7rUz/cMUkvrc4jf6Yb+Np+gyVW0lRPGB3SHblVyGRPPFlFntOdL8S+Gw9fM0UbpidzO2v7MRic/DenQu57h/bB9XNc3RMEOvuPw+gy1xDUkp+8d5Bt5sA+bLoYBMh/gbyz3I70xVzU3j4miluU5CcLU+Dg+9PdVAUH1DdZCHE38DvPjrc54HhnDHR3DY/lb98cYwDJ2v4y/XTWH2waFAFBoDc8kZe+iYPs83hcnaO3SF56JOjgzowgPMJKNhk4Gx/9L+z+yQ/eG1Ph9lQA0kFB0XxQFSQkdUHTp31r0N3kiMCePKmGeSWOxfY3b1kNAmhATz+ZWaffm5f+esXmfzl82MU1zrXAtgdkpJaM2arnfve2c+/tx0f2AZ6SUZxHXNHRDI16exmYW3JKue6F7ZRUNnkvnIfU91KiuLG/oJqAo0Glj+1pU8/J9hk4MMfLyI5IoAVL+0gLMCPV1bN5Y+fHOGtnQV9+tl9zWjQcd7YaK6cnohDSv7+VTbHfeAG6G0h/gZabI4OGW17IzLIyN9XzOh2zKa3VMpuRfGSOrONH77Wtz84hIBnbpnJuLgQ/vfjw5TVtfDKqrkcKaptW108mFlsDr7KKDvrgVtfV2+2ERNiOutUGVWNFla9uos/XjWZlQvTvNO4HlLdSorixoT4EJc7onnTLfNSOX98LHtPVPGfHSe447xRRAYZ+cPqIwyxh/shr6KhxSsr1x0S/rD6CA99cmRAtlhVwUFR3HhyfRZlfZg0LSUygF9cPB6z1c6v3k/n8qkJfHdRGm/uPDHoBqEVkBKvTjf+17fHueXlHZTUus/h5E0qOChKN7JL63ln90n3FXvJTy948bY5RAYZ+esXmYQH+PG3G6dT3tDCY+sG5yC04n0786u44plv2O5hHidvUMFBUbpgsTlYe6hv8yfde/5YJiWGsimzjI/2F/L3FTMxGXT86v30QZe5VOlbFQ0WbntlJ5uz+id/nAoOitIFP73g62N9N4A6PSWce84fjcMheXNnAW/8cD4pkYG8ubOALf10A1AGF7tD9tvTgwoOitKFjOJ6jpyxP7C3+OkFj18/DYNex9HiOu49fwyTE8Morm3mv33YjaUMfpUN/bNpkJrKqigumK12tuVWuMwZ7w33nD+GsXHOnEJTtIVTjS02Hl+XhU4n0Alc7k2sKH05OaI9FRwUxYWimmaqmyx9stdxbIiJu5eObttLuTWXzmPrMvlgXyGBRj06IXCoOayKC57sPOcNKjgoigsOCS29zNXvzv99Zyomg3PDHp2WkKe2ycqn6c5tTJosvr0JjDKwGlr6Z6KCGnNQFBdigk3ctXQ0iWH+Xr3u9bOTuXCic0N6IUTb4rqc8oY+C0bK0FLb1D+J+VRwUBQX6sxWtmSVU+TlhUffW5zWqayioYW739hLvYe/CJdN7NnWlMrQUt9io7kfni5VcFAUF4RwJsLzpohAPyYnds7amV5Y4/EgY4jJQGOL6nYa7vpj3EEFB0VxYVNmOY98luHVayZHBLosb13UpBOQGhnItOQwgs7YVrPVZVMT2J7Xf6tkFd/UH3tPqwFpRXHh0inxvL79uFevOSY22GX59txKzhkTzS3zU/kqo5TqRguLF6YRF2riiyMl7MirYnpKOFFBRjZlDe2spop7I6ICmZUa0eefo4KDorgQGWQkLMDPq9e0dpGM7cufLeGb7HJ+9PoezNqg9MZM59PEsomxvP79uewrqOG5jTlY7Wp663B3z9IxBHTxZOlNKjgoigtCCGqbvTcrZFRMEAtGRbk8Vlpn5idv728LDO19lVFGk8VOTIiJBy6byCcHi9hXUOO1dimDS0KYP1fPTOyXz1JjDoriQnFtM5UNFq9dL8Tfjyunuf6f+v8+y6Cmm+mJRTXNrD5QxEOfHOXCiXF9vreE4rvuPG8UJkPfPzWACg6K0klhdROVDRYaLd5bbJRZUkeQqfP/1GV1ZnLKG7o9t/12ml8eKWHBqEivtUsZPGJDTKyYl9pvn6eCg6KcITkikDGxwS67eXrLbHVwvLKxU3lsqD/z0lx3N7lysLC2y+4pZWi7c8lo/P3656kBVHBQlE6e+DKTn797wOvX/d1Hh11uPH/uuOgeXWf90VLuXDIK1bs0fMSEmLh1fv89NYAakFaUTrLLGvpkQ5Wd+VXsyq/inLEdg8GISNfrH7qSXlhLbbOVn144Fr0QVDa2sK+ghvRCtaXoUPXDc0b261MDqOCgKB1IKXnhttlsy63glpd3ev36gS7GHU60G1Pw1InKJp76KhuAEH8D542NYfHoaF7YnHvWbVQGllGvY1RMEBPiQ5iQEMrCUVFMTeq8sr6vqeCgKO1UNlrYkVfJzryqPrl+bIjJZXlSeECvV73Wm22sPVTMzXNTWDgqSq2gHoQC/PRcPDmOq6Ynsmh0dL+sY3BHBQdFaSc62MQV0xLJLes8eHy2Ao16EsMCOpXvK6hmcmIo4+NDzmpb0h35VXx8zyIeXH2ErNIGjpXUqQ2DfNyCUZHcNDeFiybFez2X19nyrdYoio9wN720N8bHh6A7YxRZSsmn6UXkV/S8a+lM+RWNFFY3Mz4+lJzyBiYlhnL4VN9sc6qcnUWjo7h/2TjmjfTdackqOCjKGRpbbBwp8v7g7sSE0A7vzVY732SXeyUwtFqbXozF5lBBwcf94cpJTIgPdV9xALmdyiqESBFCbBRCZAghjggh7tPKI4UQ64UQ2do/I7RyIYR4WgiRI4RIF0LManetVVr9bCHEqnbls4UQh7RznhbCuT1WV5+hKH1pU2Y5eeXe71aaGB/S4X1OWT0Prjns1c/4JL2IOWnqfxNftyHD9xMoerLOwQb8Qko5EVgA3COEmAT8FtggpRwLbNDeA1wKjNX+3AG8AM4bPfAgMB+YBzzY7mb/gla39bxLtPKuPkNRvK7FZudYSR2jYoL65PrTU8LbXltsDt7cWUBRjXc3iz9Z1czR4nqvXlPxvm25FQPdBLfcditJKYuBYu11vRAiA0gCrgaWatVeAzYBv9HKX5dSSmCHECJcCJGg1V0vpawCEEKsBy4RQmwCQqWU27Xy14FrgM+7+QxF8bqcsgYuf3prn1w71N/QYaOfPcereHvXSa9e/zuzkqlstHCyqonrZiXxwb5TXru+4l3pJ2txOGSnMShf0qMV0kKINGAmsBOI0wJHawBp3bswCWj/t75QK+uuvNBFOd18hqJ4XW/WG3jqvHExHRLmfZJe7NXrx4b6k1fRyCcHi9hfUM2Y2GBSIjvPjFJ8Q32LjXwX6VR8iccD0kKIYOAD4H4pZZ02LOCyqosy2Ytyjwkh7sDZLUVqav8uMVeGjr4MDpdOSWh73WKzsza9yGvXXj45DqtdsjHT2Y+9eEw0f/ki02vXV/rGocJaRse43gDKF3j05CCE8MMZGN6UUn6oFZdq3UVo/2wdYSkEUtqdngwUuSlPdlHe3Wd0IKV8SUo5R0o5JyYmxpOvpCidzEmLwKj3frqxIKOeCyeefug9WlRHndl7GV+DTX6U1pmR2k+qd/ecZFZqePcnKQPuYKFv78vhyWwlAbwCZEgpn2h3aA3QOuNoFbC6XflKbdbSAqBW6xJaB1wshIjQBqIvBtZpx+qFEAu0z1p5xrVcfYaieN3a9GIsXezWdjaunZXUIS/OkSLvTTPVCYgI9ONYyelBaKtdkllSrxLz+bgjPj7d2JNupcXA7cAhIURrqsoHgEeBd4UQPwAKgBu0Y58BlwE5QBPwPQApZZUQ4mFgt1bvT62D08DdwL+BAJwD0Z9r5V19hqJ43adnMQ5gMuhoOSPj6lXTE0mLDuKuJaM6lK8+4L2BYoNOx8HCGiIC/ahotzlRo8XOjORwDvj4r9PhrLfpUvqLJ7OVtuJ6XADgQhf1JXBPF9d6FXjVRfkeYIqL8kpXn6Eo3lZWZ6aiwbNppULAnBERTEsOZ0pSKFMSwxgVE8zbuwqIDjaRHBFASmRgpz2oT1Y1kRDmz6lq790ULHYHu49XMyUptENwAJA9G7pT+pnR4Ns7JqgV0ooCIMBPL7Dau7+hJoUH8PTNM5k9ovNCs9sWjOjyvLzyBhwSPj9cQmpUIEW15rNucnuuVkQfLKxl/shIdub3TRJB5eyYrfaBbkK3fDt0KUo/Wb2/yG1gAGfaA1eBoT27Q1JvtpJf3kiTttVooNHAe3tPklFc16/pl3fmVzE5MZTwQD/3lZV+5evBQT05KArO1cuBRj1NFtf/w353URrLJsZ16CpqstjYllPJjNRwPjlYRGZJPbuPV3GisgmjXmB1SG6dP4L7LhxLXKiJxLAAHluXyZ3njcKgE9j6KWXqkaI6YkJMjIkNJqfM+wkFld45c4zK16jgoCjA3LQIblswgpe25HUoD/E38D+XTmRiQggzUsJpXd9TVmfm7V0nePKrHIKMehrPCCqtN/5/bzvOknExLBoTxQUTYkkvrOXFLXmMigkiq7T/btTl9S3UNFlIDg+g0McHQocLs9WOlJJu1owNKNWtpCiAEII7zhvVthlP66Dz69+fx3dmJXUIDA6H5N6397elpzgzMJzpkc8yyC5t4OtjZUQE+jE6NrhfA0MrvU4Q4GInOmVgOCQedWUOFBUclCEpr7yB5zflYO3BuoXoYBPv3rmQu5aM5o5zR/HfOxfi76fH30/f4dddWX0LS8fFUFDl2S/wnLIGnlyfxcbMMorrzF6drdQTZqsDk4/PkBlu6szWgW5Cl1S3kjLkVDS08Ojnx9h/soYQk4HbF6Z5fG5adBC/Xj6+LSHamXsw1DVbeHpDFltzPN+KUwjYlV9FfYv3VkX3VmldCxPigzlWosYefEGzm6fOgaR+RihDzrbcStZnlPL3m2Zww5wU9yecobtMmd9kV9JktVNQ5XkeJinxicAAzrGHrNIG5vvwDmTDSU+ebPubenJQhpyrpieSGhnIjBTv5heSUnLBhFhe3JLr1ev2N4d0TnGdlBBKSa2ZqiaL+5OUPtFfM9Z6QwUHZUjydmAAeGVrPicqm0gv9P4WogOhsLqpQ86n3kqLCiQ2xB8Hkj3Hq73QsuHDl9c6qOCgKMDhU7WcqmkmIcyfqUlhnaYXVjVaOFZSz8kedCf5OiEEsaEmyup7txtdYrg/AX56cssbOV7ZxLTk/lvcN1T48loHFRyUYa+2ycotL+9oS6N9yeR4fnPpBEZGn94u9JvschJC/fl4X2FXlxl0aput1J6y9irFxqzUcI6V1FNkOZ0GpLJRdU/1VIvVd4ODGpBWhoxPDhbx+48PUd2Dm5TdIbnvv/s77K/wxZESjhWfzlVUZ7ay70Q1yZEBLB4b7dU2+4LsHq6anjcykn0FNZ1Wk5fUmlEzZXumxea73UrqP6UyZExLDmN/QQ1v7y7w+JzH1mWyKbO8U3lsqKnt9ZoDRYyICuLB1UfYnOX7G8P3VFWjhTGxQe4r4gwMu7p4yrA7JEnhamvSnmhWYw6K0vdGRAXxyb3neJyoes3BIv6xuePMo2CTgZULRzB7hHOqZ0FlEyermjhaVIvZh/uHz1aIv/vEfPNGRnQZGFpFBZs44eHiQAVqmtQiOEXpF92tUWgvs6SeX79/sO19dLCRe88fw3Wzk9tulHaH5FhJHQ4J3/Rg0dtg1NjNOgwBzB4R4dHOZXq1/VyP9KQLtL+p4KAMS/kVjZitDvQ6wfWzkvnF8nGEBxjbNmCpbbLw6BeZvLOrYFhsmZNV2sDctAh2H69mXFwwFfUt6HU6ooKNNLTY2HPCsymq9V7cG3s48JXFka6o4KAMS8snx/HB3QsxWx0sHtNxkPlEZSMldWZMBt2wCAytDhfVkRweQF2zjbgwf8rrWzrsTe2Jmmbf7SbxRTofzcgKakBaGaaEEMweEcmi0VEdyvcVVPPo58eICjKxKbNsgFo3MJotdhosNkrqzGQU13fadtQTpXVm5o2MIDLQ2ActHHqcuyr7JhUclGGtfRruN3ee4NaXd+KQkl+8e4DjlUNnwZunznaAVErYle+c9qu4V+3DqUtUt5KiAB/sK+R3Hx0GYN2R0gFuzeBn0An8dM7d8JSuGfS++/vcd1umKP1oYkIoPtz9O+jsK6ghNSpwoJvh86Yk+m7KERUcFAXILW8gMUx1hXhThBp3cOtcH15xr4KDogDLJsYhhLM7RPEO9STWvSlJoaRE+u7TlRpzUBQgyGTgnDHRNFns7MqvoqTO7P4kpVu1alprtybGh7qvNIBUcFAUoMli49yxMTRb7dilZG168UA3adA7qaXREDCs1ot4KiLIt7vdVLeSogB55Y1Y7Q4unBDLdbOSmJd2ehvNYJOhQ0I5X+4n9iXNVjvxof5MVfs8uOTrYzIqOCgKMCY2mF3Hq4gIMjIjJYInV8wgJsSZmXX55Hjqzc4ukvBAP6YmqZudp0rqzBTVqER8rkQH+3ZwUN1KigL4++mxaFlXI4OMlNWbuWlOCueOjaa6yUpxbTNxof78cvl4fvr2/gFu7eBS0WAhItCPah/OQDoQooNN7isNIBUcFEXTfnLNmgNF3DQ3hZgQE00WO8snx7Epq5yrn93aq7QSw11ieIAKDmfy8dlcbruVhBCvCiHKhBCH25XNEELsEEIcEELsEULM08qFEOJpIUSOECJdCDGr3TmrhBDZ2p9V7cpnCyEOaec8LbR8BkKISCHEeq3+eiFEhHe/uqJ0NHdkJM9syObB1YepN9tIiQxErxMYDTpO1TSzaHQUyybGDXQzB6Ugk36gm+Bzskt7ltSwv3ny5PBv4Fng9XZlfwUeklJ+LoS4THu/FLgUGKv9mQ+8AMwXQkQCDwJzcE5c2CuEWCOlrNbq3AHsAD4DLgE+B34LbJBSPiqE+K32/jdn9W0VpRs3zklBStmWb6nFZsdql3x1tJSDhTWkRQVxx3mj+GBfIVa7mn/jqalJYZgtQ3ejpN7y9b9Dbp8cpJRbgDO3f5JA6yTdMKBIe3018Lp02gGECyESgOXAeilllRYQ1gOXaMdCpZTbpTM94evANe2u9Zr2+rV25YrSZ4QQNLTYMFvtPL8xl0v/voWaJgtldS388ZMjXPHMVp//n9r3SNJP1bqtdfWMRPz0nvW1zBsZSVyob/fZuxMa4H73vYHU2zGH+4F1QojHcQaYRVp5EnCyXb1Cray78kIX5QBxUspiACllsRAitpdtVRSPldSayato4LmNOTx2/XSarXb++MlRjAYdUkKTxXf3/PVVgUbPbjPnjIlm9YGibussHhPFyoVp7M6v4oIJsYT6+2FzOFh9oIi9Hm5I5CvWHDjF7QtGDHQzutTb4HA38DMp5QdCiBuBV4BluB5ikb0o7xEhxB04u6ZITU3t6emKAsDBkzVc98I2zp8QyyWT47nq2a1MTHA+IFuG8P7RfWVUdBDRwUbyKxo9qm93SAw6ga1dJlejXofF7vx3n68CS/AAACAASURBVBIZwB3njeZHr+1pKxvMdh+v5uDJGqanhA90U1zq7TqHVcCH2uv3gHna60IgpV29ZJxdTt2VJ7soByjVup3Q/tnlzitSypeklHOklHNiYmJ69YWU4a3JYsNPr+O3l05gQ0Yp/7v6CBUNFr7Jrhjopg1KY2KDSYsOIshkoKy+xaNzjhbXERfqD8CIqEDGx4UQ4n/69+vs1AjWphcNicDQav1R300P39vgUAQs0V5fAGRrr9cAK7VZSwuAWq1raB1wsRAiQpt1dDGwTjtWL4RYoM1SWgmsbnet1llNq9qVK4rXlda18N7ek+SWN6C2IDh74+ND+O6iNLJKGzw+50hRHZMSnU9qUoLV7uiwFmDR6Gh2Hx9cXUfu7Cvw3e/jyVTWt4HtwHghRKEQ4gfAj4C/CSEOAv8PrUsH52yjPCAHeBn4MYCUsgp4GNit/fmTVgbOLqp/aufk4pypBPAocJEQIhu4SHuvKF5ntTuw2BzMSAnn7V0n3Z+guLXxWBn7TlSzeEyU+8qa9MIapmmrzwuqmhgdG4zJ7/QtakRU4JBbbZ1f0eizW4W6HXOQUt7cxaHZLupK4J4urvMq8KqL8j3AFBfllcCF7tqnKGejsqGFv36RSYBRz878MyflKa7Eh/oTHWLk8Km6Lus0Wew8tSGbu5aMIthkoKHF5va6VrskLPD0DJ5Ao570wtOznA6dqmVKUtigG3juTnGtGbPVQYDR99aBqNxKyrBVUmumuNZMQVUT/952nIzirm92ymkldWYigzybRvrhvlM8f8tMJsSHeFRfAK1bagT4dbxh7iuoHnJJDyfEh2Ay+OZt2DdbpSh9zGp38Gl6Ef/ZfgKDh3PrldN25FUyf2Sk23pl9S18fLCIP109hSXjYlg6PoaXV87p8oaYXdZAWlQQADVNVi6fltB27ODJWmb46MyenjAadFwwIZa3f7SAz+87F52PbjClcispw5KfXsfmrHI1G6kbUUFGWmwOl11CFpuD8EA/Avz0NFu7X/thttqZNzKSeSPnYbE5+NHre2jpYmrw4VO1TE0OI6+ikWMldXxw9yIOFdYyOiYIg17H+PgQIoOMVDUOrvxWI6OD+PHS0cxMjSAtKhCD3vd/l6vgoPS7zw4V8/SGbIwGHXPTIjHoBIvGRLNkXN9PQ7bZHXx5tJQLJsQyMzVCBYduRAYZeejqyfzsvwcores8HXVDRhlXTEvgYzcL16a0S3H+wb5CNmeVd1n3WEk9P79oPIFGA8kRAUQFm9jy6/OpaGghyGigtM7MK6vmcNOLO3xmSuvEhNAuuyQDjXoeuGwiK+amDIqA0J4KDkq/kFKyOaucV7bmd7ghtw44vrglj6tnJPLY9dMx9mEfrEGvY8m4GNYcLCKnzLcTnw207LIGHluXyYyUcNb9//bOOzyqKn38nzMzmUx6771RQocAoSkgVlCsqNhwce26rq5dt/tby1fddXV1d9UVuyxiW0XEjlJDDZ0AIZX03jNzfn/MTZgkk2QISSblfJ4nT2bOPffed87cue895217O/rjN1skFXVNhHqbuiyrmhR80t6QU1bb5TlrG81E+Lnxl0vHtWlvcWmNDfQgFg8ePH8UT67ZPyBSmYwJ92b+qCBe+u5Ih23/uGYyc0cOzuQOSjkMEirqmjDqdZilxNPVwP78SkK9TU4vNVjfZOaTnbmtniv1TWYMekFqfAB+7kYOFVRRWd/MDwcL2ZXTdX6dT3bm4eaiZ/nsOBKCPPtkLTanrJYIXzeySmr5Iv1Erx9/qLEjq5wFo0M6RC638N3BIp5bMoG3Nx1ne1a53WNsPlrC2cnWbLbLZsaxdm8BGYUn4x/8PYwEehoZEeKFn7vRoUzWy2fHcV1qDK/8cIRXfjji1LQmq7blMCXGjz8tHsMrPxwl18bddqDXbOgKMVB9bHtKSkqKTEtLc7YYvUZ2aS07s8u5b+UuvN1cuPfsEVw8KZxb394OwIobp7ZmEW2htrGZRz/aw9gIH26cGdsnN1kpJbtzKvj9Z3vZ0clNoaeYXHR8csdsRjro4eII/03LRicED69OZ3ZSIN8dLGSIXfp9hqtBxw0zY5kQ6csL3xzmYLtU0x/cnEpMgAepf/nG7v7zRgbxnxuntb4vrm5g2/EyEoM9ifRzw9Vwem6c1Q3N/PvHo7z8wxGnpjlJCPLgpjnxPPJReuu19c19Z5IQ5Ok0mewhhNgmpUzprp+aOQxgtmeVce2rm3HXfKCLqxv4KaOINzYc41BBNXqdYPX2XOYkBeLt5kJ5bRPBXq64uej5cs8JjHodopf1QmV9ExaL5H+783lyzQGH/NdPFT93I2nHS0kM9kTfC4pt7d4T3L9qd+v7bw90molFYYeGZgv/Xn+U3184hscXJXP3+zvaGIT/suYAq2+bSXKYN/u0tfcpMX7cd84IYgM8CPMxtTleoKcr544J7TX5PF0N/PrsEVw/I4YnvtjP6u25vXbsU+FIUQ0/HS5m3sjg1mvM2zSwM692hZo5DFD25lVw9b82UVnf9c3X3ainttGMyUXH+Ehfnr1iAnnldfi4uxDp546na9f632KRHCyoYm9eJau2ZXPwRBVxgR5cOTWKCyeEoxMCg06wK6eCD7ZmUdtoZmtmqV0DZW+xYHQwvz57BAlBnphcTj846I53t/P57vxekGx4IwS8fM1kthwr4/Wfj7XZ9sq1k6mqb+b+VbuZmRDAK9dNcdqN8cCJSh5ctbvbZcy+wGjQ8fiiZB7/eA8GnWDrowucvvTbHjVzGMQcKarm+te2dKsY4GQK6fomC1uOlXLTijQSQzz5y6Xj8HQ1YLZIa2CRTtDYbGk19mYUVlHbaOaxj/e0iUIFKMsqZ3tWOc+tO0SErxuNZgvHimqo6ad13YRgT8aE+3Tf0QFyy+vQ9fb0aZgiJfz6g13899YZrNmTT37FSSP0U18e5Nv7zsTLZGBqrH+fKIZtx8uYFOXb7TLpqFBvPr5jFmv3nuCPn+0jr6JzY3lvY9AJzGYLQliD+AaKR1VPUMphgJFdWsu1r26mpId+3AcLqiiva+T+/+7i7ORQYgLcifJzx9Nk4LNdeeSX1/Hv9ceoazJ3amRsoaCyoU9nCJ1R1IvnNOgE/9vdtaulwnHqmsw88fl+nrhkLKu25bQa9Y+X1PD814f59YKkbuMeAEprGvnDZ3t5fslEuzf7ndnlPP3lAV5cOhl/7ck7s7gGP3cX4h1YwxdCcN7YMFLjA3h3SxbvbMpqYyjuKxKDPTlSVMOcpCBS4/1Zt6+AawdwzYauUMphAGG2SK59bXObJ7KeUFDZwNq9BazdW4BOQJiPG56uBg4WVLXJj9+VYnAmgV6upGWWkhLbfQRud3yyM1cZnnuZjUdL8DIZ+PvSSRj1uzlUUM2LSycR7G1CCOFQcZ9V27L5ZGceEb5uPHDeqA7b4wI92JpZyiX/+Jm3l08nyt+dstpGjpfWOqQcWvB1N3L73ESunxHLs18dZMWGzD7Nujs52o8x4d7cMS+RAE8jy1ekcfGkiG6Xdwcig0/iIczHO3I5XtK1H/ipYpG0eWIaDNNcHzeXXlEMUko+3OYc4+RQZ2e21UPtr1dN6tH+O7LKiQ/yYFai/VxJPm4uTIvz5+eMEpb8cyOJwZ7szqngljPjmdeDuAFPVwO/u3AMt89N5JpXN51SKnFHCfQ0ctboYOYkWYM59+RWEOTpSnlt46BUDoMrZG8Isy+vksc/2eNsMQYEx0scqxzWHfvyKzu4XSpOn2AvV16+dvIpu6BuzSzleEkNhwuqWLPnBPcsGNGpcgC4cHw4YM1cuv5wMRV1TXx9msVx/NxdeOSC0Syd3nsVIwM8jNw2N4GfHpzfqhjAGhn+7JIJRPq599q5+pPBp86GIBuPlLB8xVZVn1ijpRrY6fLxDjVr6G1CvU386/opjI90PAHeN/sLWLv3BCvTcnjh6kkEe7myYHQIC8eFdbnf9PiOtSCySns2szZbJM0WCwadjrkjg5k7Mpi4AA9e+eFIj+17YT4mls+O48qpUXgNYpfVzlDKwcm0JCJTiuEkveG+2hKLoYCbz4hnaqw/v/9072kZZSdE+vDOL1NPeYnkrU3H+f5gEb7uLhRW1nPRhHCmx/l3CN5sz8YjJR3aSmoaKatpxM/DyIaMYnbnVnDN9Ogub86NzVbvofYznV+eEc/V06O57rXNpxTI6ePmwu1zE1g2K/a0A/gGMko5OJn9+ZV9Ekg2mNl+3Fp4PSHYs8drtTuyy0/bsD9U+NePR6ltbObD22bynw3HeG39sR45I0yO8evR93HZ5EhmxAdw6eTIVlfqo0U1xAV6dOmWeqSoo11ASvjxcBGLJ0bwyc48PkjLpqnZwl1nJXV6nK5ydXm6Gljxi2mM//1X3X6OAA8jS6dHc+f8xCGtFFpQNgcn81OGygranqr6ZiZE+VJQWU9NDxWnvRvLcObtTVlc8MJ64gM9+N/ds5ka6+fwvhdOCOfpy8azO6eC+m7cVLNLaylsl4RPCKhpNONlMuBtMvDlnhNc//qWLl1e6xrNfHfQfiT7v9cfRUrJBK22w6HC6tMqtemIHWNchA9f33sm950zclgoBlAzB6dT0EU2y+FKZX0TwGnlpCmq6v/4jIFOaU0jD36YzoRIH35/0RgOF1TzlzX7Katt6nSf2+cmcP+5IxFCcPGkiC6fwj/Zmcsjq9OZHh/A68umtrb/dLiY97dm8+8fj+Jm1FNa08i1qdF4aLMQs0Xy6a5cJkX5ERtoLfTjZtRz9/wkHv0ovUPw5Z7cSs7763revzmV7w4W8kV6PrfPTWB0mHebfrnldfzli/1kFFZz79kjOMdOyo6Cynr+/Pl+AC6dHMFD542ioLKBa1/bTEXdyXG59cyEARfp3NeomYOTCfNxc7YIAw5Hgqi6o7x2cBWD6U925VRw6csb2Ha8jA9vm8mvzkqyW5ltSoxfq2KAzpdnahubeWDVLn71/k5qGs3EBLT1zqnSZn91TWZKaxoRApZOi8FikUgp0esE80eFtCqGFs4bG4pbJzETBwuqeG9LFsdLajDqdfxkkwbebJF8f7CQRS+s53+78zlUUMWPh4t4a2MmTTau3FX1TSz7z1ZKaxrxNhl49ILRBHubGBfpw+vLTlari/Z355wxId2M6tBDzRycjLqJdSS7tJacstrTcgE84YTI7sGElPBBWjYf7cjloztmcvmUSP7fF/tZs8ca8WzQCe47Z0S3RuP9+ZXc+e52jhRZ3Y/djXpumhPPvrxKwn1N+LobeeLisfzuwmRMLnpMBj16AXq9DotF8vtP93LH/ESCvdp6qB0rqmHtvny8TQaKqzt+l94mA3kVdTxz+QQu/sfPHC22LiPWNjaz5J8bW1PIAwR4upJRWM31M2Jx0Qru1DeZue3t7a1Feu6cn0iATXrtKTH+/O2qSdz69jbOGh3cut9wYvh94gFGTj+E9A82LBKe++rQaR3Dz33ouRb2BY1mCw9oGWtfvnYK7/5yOjfMiGH9g/OYmdB5DAJY60hf9vKGVsUAVqVjsUhiA93xdbcuw/i6Gwn2MuFtcmF7VhkL//4TR4uq0ekEiyaEI9pVcKhvMrM7p5wD+VU0WewHbc4dGcyfFo8lMdgTKSHU2zoD/+FgURvFADB/ZDDv3pRKbaOZN34+xvcHC7n17W2t9r4IXzeunxFr5xxBXDghHH/34bWc1IKaOTiZhHZTaYWV1TtyuWxKZJdBUl0xGCNSncXevEpue2cbLnodcxID+cPisd3u88nOXO5ftbtD/YS6JjPPrTvE81dObG0zWyT3fLCTZTNjmBzth14nsGgG5KmdRMJvzypj/eHiTmMQPt2Vx+HCap6+bDzzRgYxOcZqnG7v/eSiFzx4/ih0OsE/vsvgKzvG5/vPHdnBfdpskQgB2zJLefC8kZgtslfSxw8m1MzByYzoxYI2Q41fvLGVbw8U9MgTJdhr8FbgcgZ7civZkVXOi99ldPA2ak9lfROPfrSn08I6/u0Mt/VNZoK9XEnPqWB7Vhn/WTaVxODOr/uNR0vIKavrNjitsq4Jg14wf1Qw6bnWzMKpcQGYXKy3NS9XAyNCvPD3MJJbXsf3B9vWrg7ycuXlpZO5eFJEh2NbpMSo1/HbRcksfOEnvu/Ec2ooo5SDk1EG6c5paLZwy1vb2J5VdsoVvpJ7KeX3cMPP3djh5t6enVnlncbmuBp0HQr5eLgaeHxRMstmxZEaH0BwNxHwbi76ViN2V/z2wmRGh3lz7phQjmlLWz7uLlyi3eyrG5tbZzBr95ywm1dse3YZzXbaXfQ6ahrNPLvuEBV1TV3WyB6qKOXgZOpUZHSXNJkl17y6mfTcUytFOiJkYJVmHEj4ubtwZUqU3W0psX4YujG+2jPOTo31488Xj+XHB+YxLe70kiZmFjuWWytUUzLB3iY8XA1UaS7QqVraDSmtiqbZbOGdzcc77F9U1YBer+PAiY75tywWyb0f7OSwVuu6YBgGVKqFWSeTpwzS3VLfZOGWt7bz+d2zHc675OtuJMDD2OO8OUMZV4OeZbNiyS6rZUO7FBV78yqxWGSXkcup8f48c/l4CirrcTMamJ0Y2Gv1vn86XMxfvz5MlH/3M+qfMopbA+EeWzi61SZw/tgw/uJ9gGaLhTAfE3/63742RnNbymsaeXh1Ok9cMrZNvqi6JjPzRgXj6qJn45ESiqqH33WkZg5OpMls4e/fHXa2GIOC4uoG7n5vxyntc+GE8D6SZvCxcFxYqx3mRGU9S/65kbvPSiKunUNETlkd727J4lhxTae2HiEEV6REcef8JJbPjjstxZBRWE1D88nZ8+ykQFbeMoNrU2NIe2wB3/1mLtelxjA+0gdDO4X1xoZM3tqYSX2TGYNe1+p266IXGA06bpwVh1lK3tzUcdYA1lxRFXVNpOdWcMUrG6ltPLmU5W7UMzbch1kJAdx8RhzuRj0ZhcMr6l7NHJzI4YJqskvVzMFRNh8rJaOwqktjpi13zU/kyz0nhuV6sS1GvY7lc+L4PP1kIsKq+mZ+899d/HZRMg98uJtymyjpxz62po6/fW6C3UI8p0uT2dK6NPXH/+3DRSf48yVjySuv57NdeTx43igWT7TaDQI9XfnDRWNotkie+HwfKzYeR68TmC2SoqqGTuuMv3ZDCrGBHqRlliEAWzXnbTJwbWoM3+wvZM2eE5hcdDx12fg2RYr25lVy4Ys/tTmmi17HQ+d3PR51jWYMejEk4iIG/ycYxFhUibJT5oFVu7vN79NCgKcr/74+Zdi5ILbH282FQjtBgTlldfztm8O8ceM0fNvFhZydHMI3+wvZdLRjZtSe0PJUfuBEJRe/9HObWck3Bwr55ZtpXPbyBrteQTqddSZw5/wkfN1duGNeYus2ewkEhRAkhXjhotcR5GXk9xeNYXzkSQeF8ZG+fLIzj4yiapbNjGXLowtalVELvu4uXNeuvOebGzO57rXN3P/fXZTYBOZJKdl4pIRb39rG1Ce+ZvZT33LQjh1jsNHtzEEI8TqwCCiUUo61ab8LuBNoBj6XUj6gtT8MLAfMwN1SyrVa+3nA3wA98KqU8kmtPQ54H/AHtgPXSSkbhRCuwJvAFKAEuFJKmdkbH3qgkN3D3PTDme1Z5fx4qMhunhx7jIv0YVZiID8eKuq+8xDlr1dO5K73ttvdtjevktXbc7jvnJE8/vHJYlPr9hUwPtKHMeHedvez5a1Nx5kQ6dNljYfPduWRnlvBntxKGpstSGlNyHdEW6ppCVz71YIk3IzWmUB+RR1hPm4cLarGz91IkJcrL1w1iTV78vEyGaiqb+bjnbmcMSKo0/MmBnuRGOzFdakxrNqWw57cCmYkBPDZrnzumJdIss3nM1ska/bk02S2MCsxkMcWjebiSeE89vFe9udXUttoZv3hYiL93LjoxZ8J8XbF282FvXmVbXJ5VTc08/y6Q7xy3ZRux24g48jM4Q3gPNsGIcQ8YDEwXko5Bvg/rT0ZuAoYo+3zDyGEXgihB14CzgeSgau1vgBPAc9LKZOAMqyKBe1/mZQyEXhe6zeksE3spXCcX3+w85QSFj62cDRBwzTuIT7QA0+TocvkekeLajq4CptcdKy8ZYZDRWwMOsHbm45TUt3AoU4q7y0aH853B4qI8HVjza/moNMJKuqaWpMsuugFZ4wI4pJJkWSX1vLgh7uY9eS3/HioiLpGM1/utab1mJMUyPhIX6rqrTORrhIsltU08sCqXVz28gYeXp3OpZMiuHxKFNuzynls0eg2ikFKyf3/3cWd7+7gwVXprNyajatBz+RoPwLaufbmldeRW17H9qxyvj9YZFeGtftODHobRbfKQUr5I1Darvk24EkpZYPWp2UuuBh4X0rZIKU8BmQA07S/DCnlUSllI9aZwmJhtSDNB1Zp+68ALrY51grt9SrgLNFdopdBRoSfinHoCTWNZt7fko3FwZoEI0K8+P2FY/pYqoHJGSOCcDd2nWK6uLqBExVtbV8LRoc4XHRp0fgwtmaWkVteh6+bfWXi4WrguSUTeHHpJHRCkFFYzS1vpbXe5EeEePGHi6zf0TNrD7JuXyEWaQ2IW5mWzbbjZYD16f6zXXmtx53QyWwlt7yOm95Mo6KuCYNO8NW+Ao6V1DIu0od7zx7RIb7otZ+OsVqrHNhotrS681Y1NHdIq+/IZSclrEzL7r7jAKanBukRwBwhxBNAPfAbKeVWIALYZNMvR2sDyG7XPh0IAMqllM12+ke07COlbBZCVGj9h0wBhAlRvrgb9aoKXA94/utDHCyo5KWlk7tNDgewcHwY1Q3jeOzjPTSZh4+t540NmUyJ8ePiieF8vNN6U10wOpiHLxhNRmE1d727gxEhXswfFcL6w8UkhViXYGzrPVTUNrWm2g716ehK7GVyYem0aEaGenVZ6+Bk2U9JVkkN1Q3NLBofRnpuBRG+bq2eU7MTAyiorCcpxJM75iVyqKCKGH9rEsbXfz7W6n7r5+7CDTNj7Z4rv7yOHVlluBr0RPq5sSQlkhazdHult/FICf/vi/1t2kaGWJ0eaht6/tv8149HGR3mxSWTInt8DGfSU+VgAPyAVGAqsFIIEQ/Y+5VK7M9QZBf96WZbG4QQNwM3A0RH917h8L7G2+TCg+eN4nef7nW2KIOSLcfKTqn/lVOjmRDly9o9BWQUVXO4oMq6pGInQnaw4elqoLax2e5T7SMfpfP28unkVdSz5Vhp65PwgtEhbHt8QevS0Zf3nNFmv4raJl7/+Rgvf3+ERrMFL1cDmx89q41XTwvXz4zB1aAnu7SWKP+us+kKYV1CMlskx4prmD8quI29YsnUaJZMPfk7nhx9UlGdqDi5hHP/uaMwWySvrj/KdTNi2iimcF83bpub0Oo1FB/k2amXW0ZRdZtxc3PRt9YSOV2D/EMfpiMQdlN0DHR6qhxygNXS6nKwRQhhAQK1dtvQy0igZQ5or70Y8BVCGLTZg23/lmPlCCEMgA8dl7cAkFL+C/gXQEpKyqB6LLxhZizrDxfz9f7uq1Ep2lJc3UB+RT3hvo4vz40K9WZU6Mm15mazhaPFNTyz9iCbjpa0LnMMJuaPCua5JROY8/R3duWvqm/m5rfSeHv5dDxNBjxcDXi5GiiobKCkpoE9uRWs21fIkaJq4gI98HVzobCqgYMFVW3W0z1NBtw6WWpy0eloNlvanL+8thFPV4PdiGuDXsfZY0LbuLU6gm2K+1AfV9KOl/Lnz/fz2k/HeOem6cRrN/VwXzfuP7dzt9OymkaeW3eIRy4Y3SHr6rljQojWalKcrq2qodnCIx+lsyA5ZNAlg+yptB9jtRV8L4QYARix3ug/Bd4VQjwHhANJwBass4AkzTMpF6vReqmUUgohvgMux2qHuAH4RDvHp9r7jdr2b+Xp1AIcwNw1P1Ephx5yrLjmlJRDewx6HSNCvPj39SkAFFbW80V6Pm9syCSzZHB4k206WkJaZhnv35zKda9tobSmkXkjg4gL9OStTZk0mSUFlQ0s+vtPBHu7EhvgQVltY4fU1mAdz844UVlPbnmd3Tobe/MqGRfp02rkrapv4qIXf6airolIPze8TAbmjgzm+hkxbWYep6IYymoa+Wy39dnRx82FmQmBZJbUYNAJ8ivq2ZVT3qocuqK4uoFl/9nC1dOi2Xi0mEc/Tm/dFuHrxvzRIWzPKmNytB/rD5/+KnZto5miqoZBpxy6/WaEEO9hvUGPFELkCCGWA68D8UKIPWg3dWllL7AS2Ad8CdwhpTRrs4I7gbXAfmCl1hfgQeBeIUQGVpvCa1r7a0CA1n4v8FDvfOSBx/hIH8LtrOUquseen/vpEOxtYtmsOL769Zn8esGIQREjUdto5qY301ixIZPXl00lNd6fZ5dM5LGFo9n40Fmcq1Uxa2i2kF1ax/rDxXYVQ3dISYfMpi20FNtp4Zv9hWSV1lJR18TevEo2HS3lyTUHeGR1ut39u8NikRwrqWFkqBc6Yb2Jm1z0jAr15pYz4wHszmqktFabazZbqKhr4oOtWSx5ZSMzEwKZNzKYP362j/LaJnzcXPjynjl8cEsqr64/yurtOUDv1SL3Mg0uxQAOzByklFd3sunaTvo/ATxhp/0L4As77UexejO1b68HruhOvqGAEIK5o4J5d3OWs0UZVFwxJZLZPaz30B1Gg45fLUgiJdaP29/ZPijcjlem5ZCeW8mbv5iGXgjmP/s9hVUN/GnxGKbG+rfWSj4dNh4p4dp2wWEAo8O8aWg2t677nzU6mFGhXh2S2n2x5wT3OWCXaMFskTRbLDz0YTpTYvxaldq4iJNBbXfNT6KoqoHpcQGtbSXVDfzhs33syC6jpsFanrSFMB8Tv5wTz3WvbW6dHS5JiWRUqDffHyxkd04FNQ3N1DeZWTA6mPom82nPIAbbrAFU+owBw8B/Ph1YLBgdzF8uHdfnT/azEgP55I5ZXPPqZnIHNuba7AAAIABJREFUQZLE/fmVLF+xlbd+MQ0hBLWNZu77727+syyFqbF+bM08NSN+ew4VVCGl7OAhFhPgzvGSWkZoXj5eJhdeXDqJRX//ifomi9Zm4PGFyQ4phh1ZZTywajcWKTFbJJkltXykuZp6GPXcOjehta/JRc/Tl09ofd9stnD1vzdxqMD+U/+cpED+tzuvjeJ6a9NxGpotvKM9oB0pqmHU4186MiRdkhzmTYCnkYq6JoddgwcKKn3GAKCsppFPduZ131EBwKhQL164elK3qaV7i9hAD96/OZVAz8FRLnJ3TgUPf5TOihuntT5h379qN3eflXTayvRwYbVdu4SrQd+qGFpIDPbisztn87erJvLS0slseGg+S6baTxV+rLiGu97bwQdbs7BYJMtXpHG4sJojRTUdbD+zkwI7JAxsodls4cCJKuaNCu70MzSbrQrHlvomC29uPN6h/XQ5WlzNlmOlXQbrDVTUzGEA8J8NmZ0WT1G0xcvVwMvXTrHrTtmXRPm78/iiZH71/s5+Pa8j3H1WEhuPFLeZFXyRfoL6Jgt/u2oij328hw1HSvgi/QQzEwJOe4mktKaR+M4zVrQhKcSLpJCuEyW+vek4z607RGlNI1+k51PdYGblzalsyypDSoj2t9aj9jIZKK5uINCzcw+i1TtyeWDVbh5flMzEKF92ZnesA7JufwE/9FE6FaNBx9Jp0cxKDCTK3426RjPZZXXdFlAaiCjl4GQ2Hinh9Z+OOVuMQYGnq4E3l0/r9Kmxr7lgXBgr07LJL6/nqIMFafqD/PI6bjkjga2ZaW3avz1QSGZxDX9fOomLXvyZrNIah+thdMXmY6WkdFL7uSfEB3mQGOxJXnkdOWV1/Ol/+9gzKaJNHeoWOluSqms083l6PkeKqpkU7cuXe/IJ9zWx006Qcl+5K89KDOC5JRM7jPEkmziNwYRSDk6koq6Jm99KU7MGB3npmslO/aG56HX8/erJPLBqFyYXPfvyT93jpy/44VARri4nl9g8XQ2t19TR4hpe/DaDxxeOxt3VQFkvFD/6OaO4TWZUW7YdL2XzsVIunxzZbTnQFmYmBDIlxo9NR0u54fUtABRU1lNYVU+wV/fHyCuv4/WfjrXaDWYlBlDdYOaL9BOOf6hucHPRMzHKlykxfoT7ulFW28jBE1X8nFGMEDA11p+nLx/vUC6qwYJSDk6isKqehz9MH5RBV87gljPjObOL7Jv9hb+HkXvPHsnil37qvnM/UVXfTLOWEiQp2BN3o55dORWt29fsOYGbi57nrpzIGz+f/ix1w5ESCivrO9z89+RWcMPrW6luaOadTVl8f/9ch+MYXA165iQG8otZcazYmElpTWO3JXSllLy96TjPrjvUph7Fzxk9i2oO9HRlYpQvE6N8CPYyEepjItzXRLS/BzpBv9m4BgpKOTiJ7cfL+OZAx9z1io54mQz8esEIZ4vRypGiapKCvZgW58+bGzMdSsTWl9Q1mUmND2BshA+v/3ysjWJoYfOxUqobmsnvpcJHlfXNBNtk8y6pbuDSf2xoTUWSW17He1uyuH5GrMPH1OkEv70wmbvPSqSh2dLlElhhVT13vrODLZl2kyY4zIgQT+aNCubiiRGMCvVyKE/XcEEpByfhMQj9np2Bq0HHc0smDhg3wJLqBqL93fnTxWP57Sd7nK4YAPQ6wYnKehaOC+Pr/QUcbVcv+epp0fx2UTJuRn1rHILRoOPh80exalsOKTF+XDcjlld+OMKqbTkOnXPDkWISg09GI6/entshR9XmY6WnpBxa8HXvaLxtNltYn1HM4YIqEoI8eXh1OoWn6AFk1OsYG+HNjIQAkoK9mJ0U2KVxe7ij7lBOwqAbXlPUniAEvLR0MguSQ5wtSituRj0f7cjlrU297/bYU8wWyZNrDvDpzjyevnw8OiEI8Xalsq6ZsRE+TIzy5ck1+wnzdWNMuDdPXDKWSyZF4G40cOOsuNbj3LMgyWHl8Mr3R7hiSlRrYZ7yuo62jIldFP85Fb7ZX8AjH6VTYKeaXVcY9TrGRfowJymQ1PgARod64+M+dGwCfY1SDk6isGp41zV2hN+cM3JAKQaw+sMfKqhqoxjCfUyU1ja2Bns5i335lVzz6mauS42hyWIhwteNHw4V8tSXB1r7pD22oNOnZX8PI3GBHl3mV2ohr6KeHdllzEywRqi3rwIX5OXaaUyDo1TUNfHZrjxe++mYQ4rBRS+IC/Rg3qhgzkgKIjnMG79B6EI6UFDKwUkMhmhbZzJvZBC320TBDhRyy+paA5p0Aq6fEcvFkyK4+KWfe/U8F04IJ9LPjRh/d/Q6QZNZ0mS2oNMJCivrEUB1g5nc8lrKa5sYHeZNSqwfJoOeZovkREUdX+8vJC2zjEBPI6NCvckpqyW3rK5T5eBuNHDrmfE8+GH3+Y9GhniRapOuIjbgpHvxtanR3DU/CZ9OCv90Rn2Tmf35leSU1fFFej7fHCjsUKHOFqNex4LkYC6fEsnIUG+CPF0xGtSMvLdQysEJSClZf2jI1CzqddyNeu6YlzggjIM/HCpi45ESgrxcSY33J8DTyItLJ2N1XBEkBnvyu09O1l4O8XY95eUPe5gtFtIySzmQX0lqfACp8QGMCfdu4zEjpaTZIhHY96RZNiuOoqoGvtp3gkMnqojyd2dffiUTojpf7vFxc+xJ+/Z5Cehsoq1Hhnrx1ysn4u1mYN7IYLvfnZSSHw8X886m45TVNpIc5s2oMG/qm8xU1DWxM7u808R+LRh0gsnRflw4MZxF48LUzKAPUcrBCWSV1rLxNIuIDFWMBh3/WTbVoSCrTUdLiPJ3J+I0UnZ3RnZpLX/4bC9f77fvUXbf2SNYkByClJJfLRiBp8nAvrxKcsvrekU52Prof6fdMN2NejxcDTSZLTSbJY1mC2aL5PGFo1lmYzuwJcjLlWumd0yU1xmOFLeZHufPwnFhHdq7KmhTWtPIox+ls2bPyc91KnmeAjyMXD4lkhtmxp5WinaF4yjl4AScvTY9kHl84WgmRjtmyFyZls05yaG9rhyklNzx7nZ223EJbeHZdYd4dt0hwn1MjIv0IczHjdFh3q038r6gttFst6Tstqxyls3qnXPMTAjgjQ2ZHdq9TQbmJAVx1bQo5iSdWrzJdwcLeWDV7lPOL+Rh1HPmyCDOGxvGOcmO17RW9A5KOTiBDUfUkpI9HrlgFJdPicLFQU+uESFejAztOm9PT6hrMnMgv6r7jlgNs3kVznUu2Jl9eplWbTk7OYQIXzfCfU0snhhBtL87cYEehPqYTqkwD0BtYzNPfL6/NdOpo/i6u3DnvESunBo1pCKOBxtKOTiBLcdOL3BnKHLHvARuPuPUDNC3ntk3Bmt3o4HRYV52g8mchZfJgK+7C9mlVkcGg04Q5muisLKBpmZJfZO5V56shRCsu/eM005suCu7nHs+2OmQ51MLJhcd16XGcNdZSXgrpeB0lHJwAhmFvVNdaqjw1GXjuHRypLPFaMP548L6XTm46AV+7ka83VwwWyTR/u6MjfAmKdiLhePD0AtBVX0zbkY9RoMOs2aM1vVyTYvTUQxNZgv/+O4IL3x72OE4kCAvV5bNjGXptGhlYB5AKOXQz5TWNJLRS6UHhwL3nj2CK6dGO1uMDtx6ZgILRoeQUVjNxiPFfJ6eT3H16SetM7noiAv0ZGKUD5F+7oR6m4j0cyPE29S6dKMTdOqpZRvENdBKmO7KLufBD3d3qP7WGdH+7tx8RjxXpES2Rm4rBg5KOfQze/MqkAMjsNbp3DQ7jrvm28/uORBIDPYkMdiT88aGcu85I3n0o3T+tzu/2/28TAZCvU2EaH/B3q6E+7oxKcqXuEAPXA26IZXEra7RzHPrDvLaT8ccSicSF+jBbXMTuGRSxCnbMRT9h1IO/cz+AZLm2dksGB3CIxeMHhCxDI7g4+bCNdNj8HQ14GbUU1LdiFlK3Fz0hPmYSAz2JCHIkyg/a72BloLyFinR68Sg+Zynys8ZxTy8Op2s0tpu+06L8+e61BgWjgvr9aUwRe+jlEM/89EOVQ50VKgXf1g8ZtDdIKbF+TMjIaDbfmaLbP1sukFaHby+yYwQ1lTaUkpKahrxMhnIK6+nvLaRJrPk/S1ZrNbqOnfF7MRAfn32CCZH+w5ZJTkUUcqhH8ksrhn2M4ekYE/euWk6AYMwG6aja/wDzRZwquzJreCmFWlU1DUxOcaX4qpGDhZUEeXv1uot1R06AeMiffntotFMiem9qnGK/kMph35kZZqdmoXDiCumRPL05eN7/enxx0NFJAZ7qshZB2g2W7q1dzzx+X5OaHUfbAvnOKIYZiUGsGxmHNPi/PEw6oeUbWW4oZRDP/LVvgJni+A0/Nxd+szGMDXWv8sEbQqrUvjuYBFf7ytgQXIIZ9vJdttktvDsV4dOObWLQSe4IiWK285MIDrAfo1nxeBDKYd+orSmcdjGN8QFenDLGfF95sPuZtS31hVQdCSjsJpff7CT9Fxr3MZHO3O5aEI4T1wylsLKBspqGzEadDz0YTo7s8sdPu7IEC9unBXLwvFhKpJ5CKKUQz/h6WpAJxgQlcP6k/GRPnx420wMg3wdfrBQUdfE2j0nyC6rpabBzLbjpR2C+RqbLazalsPHO3Jp7sEFOSLEkzvnJ7FwXNigt68oOkcph36itKZx2CkGgPvOGal82fuJY8U13LRiK0eKHEtZcSqKwdfdhcsmR3LJpAjGRvj0VETFIEIph37iu4P2Uz8PVXzdXbjljARmJwY6W5RhwQ+Hirj97W3U2Mna2lOMeh0zEgJYOj2aM0cEqayowwylHPoBs0WyertjtXmHCr9dlDzg8iUNNeqbzORX1PPN/gKeXHOgR0tE7dEJGBvhw9Jp0cxMCFQG5mFMt8pBCPE6sAgolFKObbftN8AzQJCUslhYXVH+BlwA1ALLpJTbtb43AI9pu/5ZSrlCa58CvAG4AV8Av5JSSiGEP/ABEAtkAkuklL2Xm7iXMFsk724+Tm2jGS+TC9Pi/PBwNeBuNOBq0HGsuIZ7V+4aVvEN8YEenDV6YNV+Hmr8Ny2b33+697RmCiHerlyZEsWMhEAq65sYEeKFv4cRb5NBBaspHJo5vAG8CLxp2yiEiALOBmyTtZ8PJGl/04GXgenajf53QAoggW1CiE+1m/3LwM3AJqzK4TxgDfAQ8I2U8kkhxEPa+wd79jH7hpqGZi588SeO2qzx6oS1ZONVU6PYdryMAyeqHM5OORQYE+7Nx3fMGqRxwYODb/YXcP+q3T3eP9TbxL1nj+DSyREqDkHRKd0qBynlj0KIWDubngceAD6xaVsMvCmllMAmIYSvECIMmAusk1KWAggh1gHnCSG+B7yllBu19jeBi7Eqh8XafgArgO8ZIMqhyWyhuLqBj3fktVEMYPVGamy28ObG406Szjksnx2n1Q5GGaD7mLc29ezaCvF25Q6tiI7Kgqrojh7ZHIQQFwG5Uspd7aafEYBtGHCO1tZVe46ddoAQKWU+gJQyXwgR3BNZ+4InPt/Pe1uyaDSrwCuwRsU+csFo5dbYTxw5xZTvfu4u3H1WEldPi1ZGZYXDnLJyEEK4A48C59jbbKdN9qD9VGW6GevSFNHRfVsboLCynoMnqmhQEbmANTr2/10yTimGfsRaJa37VBbeJgPLZ8dz4+xYVVlNccr0ZOaQAMQBLbOGSGC7EGIa1if/KJu+kUCe1j63Xfv3Wnuknf4ABUKIMG3WEAZ06gsqpfwX8C+AlJSUPlvgr28yc+d7O1SZTxuumhZFTICHs8UYVnRnJwjxdmX57DiunhatIpcVPeaUlYOUMh1oXeIRQmQCKZq30qfAnUKI97EapCu0m/ta4P8JIfy03c4BHpZSlgohqoQQqcBm4Hrg71qfT4EbgCe1/7a2jX5DSkna8TKi/d25d+VOpRiwptyeEuPHiBAvzh8b6mxxhh1mi/1Za5CXK/csSOLyKaqymuL0ccSV9T2sT/2BQogc4HdSytc66f4FVjfWDKyurDcCaErgT8BWrd8fW4zTwG2cdGVdo/2BVSmsFEIsx+oRdcUpfbJe4vl1h3jh2wxnnLpfcTXouHRyJAvHhVFR18S7W463ycjZwpQYP1beMkMtIzmRirqmNu+9TAZunBXHL+fEqZmCotcQcojVrExJSZFpaWm9cqwdWWVc8o8NvXKsgYyXycDX956Jn7sRo8G6ZGGxSJav2Mp3B4sAiPJ34zfnjCQ1PoAQb5MzxR32jPvdWqoamokNcOe6GbEsSYlUSkHhMEKIbVLKlO76qQjpLvj3+qPOFqHPMep1PH3ZeLxMhlbFAKDTCf58yTjmPvMdQgheuXYKY8JVTh1n02S2gIDHFo5m2cxYFaeg6DOUcuiCRxcm8+2BQuqbhq5n0h8Xj2HeqGC7Lo4Rvm5E+LqxZGqUUgwDhNoGM1sfXaBcUhV9jnrs6IIIXzcunhjRfcdBzJhwn05vNIcLqnAzGvjlnPh+lkrRGT7uLkoxKPoFpRy6YVzk0H5iLqqu79AmpSSrpJYb39jK/eeOUBHPCsUwRP3qu+GMpCBni9Cn2HN5LKlp5PZ3txHp58a8kQMmMF2hUPQjSjl0g9Gg48wRQ09BGPU6Nj18FuPbzYwamy08uGo3e/MqeXxRssrOqVAMU5RBuhtCvE08c8V4znn+R8prm7rfYYDjZTLwm3NGsnhiOL7u1prOD324Gw9XAxYp2ZNbwdbMMpakRCojtEIxjFHKwQGCvUyMDPFi8yCPjtYJ+Oj2WSQGe7ZpX3+4mNzyk7l6Qr1NPHLB6P4WT6FQDCDUspKDtL+hDjbGhHuz2o5iALhhZkzrayHg+Ssnts4qFArF8EQpBwfxdhu8EajjI314/+ZUJkb52t1+3pgwwDqzeOrS8cxICOhP8RQKxQBELSs5yJr0fGeL0CPmjQzi5WundOkbH+5r4tkrJuDt5sLZyaq8p0KhUMrBIY4WVZNZUutsMU4JT1cDz1w+nnPHhKLrJkmeQa/jsimRXfZRKBTDC6UcHGB/fpWzRXCYxRPDuW1uAiOCvbpVCgqFQtEZSjk4wO6ccmeL0C0BHkb+b8kEFbSmUCh6BWWQdoCrp/Vt6dGeohMwIdKHCZE+fHLnLKUYFApFr6FmDg7gaRp4wxTgYeTla6cwNdYPKVFLSAqFoldRMwcHcDXo+P2Fyfg40Z010NPITM3FNMLXjY9un8W0OH+EEEoxKBSKXmfgPRIPQLxMLkyM9iPE27VDica+xtWg475zRrB0egyerga2HCvF3agnOsC9X+VQKBTDC6UcOsFikeh0gve2ZPHpzjz25FVQVd/crzIY9TpW/GIaqfEng9Kmxfn3qwwKhWJ4opSDHSwWSUZRNWaL5NmvDlJc3egUOZ66fFwbxaBQKBT9hVIOdtibV8kvVmyluLoBKZ0jw7KZsVwySQWmKRQK56CUQzsKq+pZ8s+N1DWZnSZDarw/jy9Kdtr5FQqFQimHdryzKatPFYNRryM1IYDkMG+mx/uTGOTJpS9voKiqAYBwHxMvXD0JvfJAUigUTkQph3ZsO17WZ8dODvPm1RtSCPd1a9P+54vH8uhH6ZTUNPJ/V0wg2MvUZzIoFAqFIyjlYENWSS27+ihVhkEneOXaKR0UA8C5Y0IZF+FDaU0jYyNU9TWFQuF8lHLQaGg2s/TVTX3mrnrHvMQuYxPCfd3sKg6FQqFwBko5aBj1Okr6wGXVRS+4NjWGu+Yn9vqxFQqFoq9QykFDCIGk9/1W/++KCSyeGNHrx1UoFIq+ROVWsiHMp3eXdS4YF8pFE8J79ZgKhULRH6iZgw3XTI/mz5/vP6V9InzdGB3mTbPFwp7cSi6dHEF8oAcjQ72YGOWLEMolVaFQDD66VQ5CiNeBRUChlHKs1vYMcCHQCBwBbpRSlmvbHgaWA2bgbinlWq39POBvgB54VUr5pNYeB7wP+APbgeuklI1CCFfgTWAKUAJcKaXM7KXPbZerp0Xz9JcHaTRbuuzn6WrgzJFBLJsZy9TYk7mOpJRKGSgUiiGBI8tKbwDntWtbB4yVUo4HDgEPAwghkoGrgDHaPv8QQuiFEHrgJeB8IBm4WusL8BTwvJQyCSjDqljQ/pdJKROB57V+fYqHq4E/Lh6Dp2vnOnN6nD8/PjCPl5ZObqMYAKUYFArFkKFb5SCl/BEobdf2lZSyxedzE9CSBGgx8L6UskFKeQzIAKZpfxlSyqNSykasM4XFwno3nQ+s0vZfAVxsc6wV2utVwFmiH+6+V02L5trUmDZtXiYDXq4G5iQF8s5N0/H3MPa1GAqFQuFUesPm8AvgA+11BFZl0UKO1gaQ3a59OhAAlNsoGtv+ES37SCmbhRAVWv/iXpC5S+5ZkMSRomq2Hy/jpWsmkxofoJaMFArFsOK0lIMQ4lGgGXinpclON4n9GYrson9Xx7Inx83AzQDR0adf79nkouff16e0P8dpH1ehUCgGCz12ZRVC3IDVUH2NlK2JrXOAKJtukUBeF+3FgK8QwtCuvc2xtO0+tFveakFK+S8pZYqUMiUoKKinH0mhUCgUGj1SDprn0YPARVLKWptNnwJXCSFcNS+kJGALsBVIEkLECSGMWI3Wn2pK5Tvgcm3/G4BPbI51g/b6cuBbGyWkUCgUij7EEVfW94C5QKAQIgf4HVbvJFdgnbbcsklKeauUcq8QYiWwD+ty0x1SSrN2nDuBtVhdWV+XUu7VTvEg8L4Q4s/ADuA1rf014C0hRAbWGcNVvfB5FQqFQuEAYqg9jKekpMi0tDRni6FQKBQDEiHENillSnf9VPoMhUKhUHRAKQeFQqFQdEApB4VCoVB0QCkHhUKhUHRAKQeFQqFQdEApB4VCoVB0QCkHhUKhUHRAKQeFQqFQdEApB4VCoVB0QCkHhUKhUHRgyKXPEEIUAcf78ZSB9EONiR6iZOsZA1W2gSoXKNl6ijNki5FSdpu+esgph/5GCJHmSJ4SZ6Bk6xkDVbaBKhco2XrKQJZNLSspFAqFogNKOSgUCoWiA0o5nD7/crYAXaBk6xkDVbaBKhco2XrKgJVN2RwUCoVC0QE1c1AoFApFR6SUw/YP+BWwB9gL3KO1PQMcAHYDHwG+Nv0fBjKAg8C5Nu3naW0ZwEM27XHAZuAw8AFg1NpdtfcZ2vZYR2Sz2fYbQAKB2nsBvKAdbzcw2abvDdr5DwM32LRPAdK1fV7g5CzSH1in9V8H+DkqG3CXNg57gaf7e9w6+T4nApuAnUAaMK2/xgx4HSgE9tgco7O+/SFPyzkqsZbxzXBArms0eXYDG4AJffH92VwjFUC5I2Nms30qYAYu78MxaylXXOqobFjLK+/Eej3+0Mfj1ua31Sv3x9482GD6A8ZivZG4Y62l/TWQBJwDGLQ+TwFPaa+TgV3aFxUHHMFaD1uvvY4HjFqfZG2flcBV2utXgNu017cDr2ivrwI+cEQ2bVsU1lrcxzmpHC4A1mgXciqw2ebCPar999Net1z0W4AZ2j5rgPO19qdbLljgoZbP78C4zdNeu2r9gvtz3LqQ6yubz3YB8H1/jRlwBjCZtjeTzvr2hzwt5zgD682z1gG5Ztqc73wbuXrt+6PtNXIlkO3ImNnI8S3wBZpy6KMxE8CtWBWLI9+nL7APiG73e+ircWv9bfXaPdLZN2ln/QFXAK/avH8ceKBdn0uAd7TXDwMP22xbq11kM4C1Nu0Pa38Ca3BLi6Jp7deyr/baoPUTjsgGrAImAJmcVA7/BK626X8QCAOuBv5p0/5PrS0MOGDT3tqvZV/tdRhw0JFx0y7sBXbGuV/GrQu51gJX2nzOd/tzzIBY2t5M7PbtD3lsz6HJ1WDTr8vvXWv3A3Lbfy+n+/3ZuUZ+AI50N2ba+3uAO4A3OKkc+mTMtPdH2h2js/1uB/5sZwz7ctxa+/XG33C2OewBzhBCBAgh3LE+IUS16/MLrE8NABFYn2hayNHaOmsPAMqllM3t2tscS9teofXvUjYhxEVYf5y72sl5qrJFaK/btwOESCnzNdnygeB25+ps3EYAc4QQm4UQPwghpvZQtp6OW2dy3QM8I4TIBv4P6w+qv8fMls769oc87Y/VdIqfYTnd/x568v21P9YJwMXmvV3ZhBARWB/gXmknZ1+O2QmsN+guZcP6e/ATQnwvhNgmhLi+G9l6Y9xs9zltDN13GZpIKfcLIZ7Cuk5YjXV61vLFIIR4VHv/TkuTvcNg36gvu+jf1bG6k+1RrMte7enseKfa3i1dyGbA+mSZinUNeKUQIr6Lc/XquHUh123Ar6WUHwohlgCvAQu6OFavj5mD9Ic8Pf4MQoh5WJXD7G6O1ZPvz167I/wVeFBKaRaizSH6eswcwYDVvnEW4AZsFEJs6qFsTrkmh/PMASnla1LKyVLKM7Aamg4DCCFuABYB10htvoZVK9vOLCKBvC7aiwFfIYShXXubY2nbfbTzdyVbJtZ1xV1CiEzteNuFEKE9kC1He92+HaBACBGmyRaG1aDqyLjlAKullS2ABWvemH4bt07kugFYrXX5LzCt/bH6Y8xs6Kxvf8jT/lgujnwGIcR44FVgsZSypBt5e/L9tT9WKNZZTQudyZYCvK/9Hi4H/iGEuLgL2XpjzEKxeYjsZr8vpZQ1Uspi4Eesy8F9OW62+5w+vbU+NRj/OGkkisbqoeSH1ZNgHxDUru8Y2hp/jmI1Lhm013GcNDCN0fb5L20NTLdrr++grYFppSOytdueyUmbw0LaGjO3aO3+wDHtc/lpr/21bVu1vi2GuQu09mdoa2B72sFxuxX4o9Y+Aut0V/TnuHUi135grtZ+FrCtP8eMjjYHu337Q55257iEtgbpzvaJxuoNM7PdWPfa90fHayTLkTFrJ88btDVI98WYpWL1PnLk+xwNfKONkzvWZc+xfTxuR1EG6V5TDuuxKoJdwFlaWwbWG9tO7e8Vm/6PYjVIHUTzctDaLwB/e6xrAAAA/0lEQVQOadsetWmPx+odkaF98S2ePCbtfYa2Pd4R2dptz6StK+tL2vnTgRSbfr/QzpMB3GjTnqJdsEeAFzlp2A3QLurD2n9/B8fNCLytHXM7ML+/x60TuWYD27S2zcCU/hoz4D0gH+tTcA7WZZnO+vaHPC3nqNJkanZArleBMk7+HtL64vuzuUYqsT4Vdztm7b77N2jrytrbY3YEq4ttsaOyAfdjvR730Nbluy/Grc1vqzf+VIS0QqFQKDowrG0OCoVCobCPUg4KhUKh6IBSDgqFQqHogFIOCoVCoeiAUg4KhUKh6IBSDgqFQqHogFIOCoVCoeiAUg4KhUKh6MD/BwJIUhrXbUSyAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "df.plot(figsize=(6, 6))\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "One thing to notice is that the values of the geometry do not directly represent the values of latitude of longitude in geographic coordinate system\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'init': 'epsg:2263'}\n" ] } ], "source": [ "print(df.crs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As folium(i.e. leaflet.js) by default takes input of values of latitude and longitude, we need to project the geometry first" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "{'init': 'epsg:4326', 'no_defs': True}\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
BoroCodeBoroNameShape_LengShape_Areageometry
05Staten Island330470.0103321.623820e+09(POLYGON ((-74.05050806403247 40.5664220341607...
14Queens896344.0477633.045213e+09(POLYGON ((-73.83668274106707 40.5949466970158...
23Brooklyn741080.5231661.937479e+09(POLYGON ((-73.86706149472118 40.5820879767934...
31Manhattan359299.0964716.364715e+08(POLYGON ((-74.01092841268031 40.6844914725429...
42Bronx464392.9918241.186925e+09(POLYGON ((-73.89680883223774 40.7958084451597...
\n", "
" ], "text/plain": [ " BoroCode BoroName Shape_Leng Shape_Area \\\n", "0 5 Staten Island 330470.010332 1.623820e+09 \n", "1 4 Queens 896344.047763 3.045213e+09 \n", "2 3 Brooklyn 741080.523166 1.937479e+09 \n", "3 1 Manhattan 359299.096471 6.364715e+08 \n", "4 2 Bronx 464392.991824 1.186925e+09 \n", "\n", " geometry \n", "0 (POLYGON ((-74.05050806403247 40.5664220341607... \n", "1 (POLYGON ((-73.83668274106707 40.5949466970158... \n", "2 (POLYGON ((-73.86706149472118 40.5820879767934... \n", "3 (POLYGON ((-74.01092841268031 40.6844914725429... \n", "4 (POLYGON ((-73.89680883223774 40.7958084451597... " ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df = df.to_crs(epsg=4326)\n", "print(df.crs)\n", "df.head()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEfCAYAAACtRRYAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XdAnPX9wPH3c4tj771JyN4hO5oYs9Q46qrGWfeqVWu11l9ra7WttraO1r1nTNwaZ2L2BhIyCJCEvTccHMet5/fHERLggCPcwQHf119w94wvJHye577P5/v5SLIsIwiCIIwMisEegCAIgjBwRNAXBEEYQUTQFwRBGEFE0BcEQRhBRNAXBEEYQUTQFwRBGEFE0BcEQRhBRNAXBEEYQUTQFwRBGEFUgz2AzkJCQuSEhITBHoYgCMKQkpaWVi3Lcmhv27ld0E9ISCA1NXWwhyEIgjCkSJJU4Mh2YnpHEARhBBFBXxAEYQQRQV8QBGEEEUFfEARhBBFBXxAEYQQRQV8QBGEEEUFfEARhBBFBXxAEYQQRQV8QBGEEEUFfEIawnAodB4rqsVrlwR6KMES4XRkGQRAcY7XK3PzOPopqW4jw03L+5EhWTY1kWkwACoU02MMT3JQI+oIwRLWarZTVGwAobzTw5o483tyRh7dGybxRwSwZF86isaFEB3gO8kgFdyKCviAMUWkFdZjtTOs0Gy1sOFrJhqOVAIyP9OPccWGcMy6MabEBKJ30KaCgppnN2VWsnhOHWilmiocKEfQFYYj64Ui5Q9sdLWvkaFkj/910nAAvNYvGhLJkXBhLxoXhq1Wf8fn/+UM2v1qQgFISU0lDiQj6gjAEmSxWvjtc1uf96vUmvjxQypcHStGoFMxJDGLRmFAWjw1jVKg3kgMBXJZl3t6ZT2p+HTkVOhSSREKwN1fNjmVuUjBatfJMfiRhgEiy7NhTf0mSlEAqUCLL8ipJkhKBNUAQkA5cJ8uysdM+GuAVIAWwAr+RZXlzT+dJSUmRRT19QejZ5/uLuf/jDKceMyHYi8Vjw1g0NpR5PQTvnAodq1/bQ3VTa5f3LpoaxXXz4okK8BTPEgaYJElpsiyn9LpdH4L+A9iCt19b0F8LfCbL8hpJkl4GMmRZfqnTPncDKbIs/0qSpDDgO2CWLMvW7s4jgr4g9MxssTL37xupbjL2vvEZ0qoVzE0KZsXECM4dH0aYrxaA45U61uwt4vXteT3uf+HUKF64errLxid05WjQd2h6R5KkGOAC4EngAcn2GXAJsLptk3eAPwMvddp1ArARQJblSkmS6rFdOPY6cl5BELqq1LW6NOADGExWNmdXsTm7CrA9DF40JhSFRK8BH2BLdiUVjQbC/bQuHafQd44+cn8WeAjbFA1AMFAvy7K57ftiINrOfhnAxZIkqdqmg2YCsZ03kiTpNkmSUiVJSq2qqurTDyAII83XGaUDfs6jZY28vOUEL24+4dD2jQYzz2885uJRCWei16AvSdIqoFKW5bTTX7azqb15ojexXRBSsV04dgLmzhvJsvyqLMspsiynhIb22tdXEEa0z/eXDPYQHBIT6DXYQxDscOROfwFwkSRJ+dge3C7BFsADJEk6OT0UA3S5/ZBl2SzL8v2yLE+TZfliIAAQl39BOEOHSxrIKtcN9jAccibZRYLr9Rr0ZVl+RJblGFmWE4CrgJ9lWb4G2ARc3rbZDcCXnfeVJMlLkiTvtq+XAWZZljOdNXhBGGnWphYN9hAcdrC4gbs/SGfN3kIqGg2DPRyhTX/y9B8G1kiS9ASwH3gDQJKki7Bl7PwJCAN+kCTJCpQA1/VzvIIwYhnN1kGZz++P9YfKWH/Idsc/IdKvPRtoYpQfkiTRaDDh66FyaH2A4BwOp2wOFJGyKQj2fXuojLs+SB/sYThFhJ+WRWNCKWs0kF3eyDljw1g+MZyUhCD8+rFKeCRzasqmIAiDbyhN7fSmvNHAx6f9PGv2FbFmXxEqhcTEaH/mJAYxLymYmQmB4iLgZCLoC8IQkFfdzJac4Z/ObLbKZBTVk1FUz6tbc1EpJKbGBpASH8jsxCDmjQrGSyPCVn+I354gDAGvbs3FzWZiB4TZKpNWUEdaQR2vtF0EpsT4s2B0CPOSgpkRHyhq/fSRCPqC4OYqGg18mlY82MNwC2arTHphPemF9bzw83HOSg7hvZvnDPawhhRRBFsQ3Nzr23IxWrotVzVihfl68PjFkwB4b1c+ZouVrSNgCqy/xJ2+ILix5lYzJotMiI+mS70dT7WSyTH+FNQ0U9HYteLlcDc1NoCGFhMXPL8NncEMkkS93kiQt4ZJ0f6DPTy3Je70BcENWawyn6QV89aOPCL8tXxyx3wSgk+VNZgS40+or4a9ebXoW81MjPIbxNEOPKVC4vazk9iQWcGR0kYKa/X88YvDPPNjDi/8bFv0bxafjuwSd/qC4IYe+uQgn6bb5vE1SgUXTY3i4fPGoVJIGExWdp6o5ucsWztEXasFzxH2MPO2s5NISQhif2F9h9ej/LX8bsU4DCYLHqqO97R6o1lk/iCCviC4nR+PlLcHfACjxcrevFq+OViGv6ea9YdKiQ304tLp0YwK8+HRzw9hMFsGccQD68u7FzAlxp/86mY+2FPQ4b3HLpyAyWKludWMVu3R4T2LnX7CI5FYkSsIbqTFaGHJM5spa3CsVs3ti5LwVCt5dsPwrGOoUkiYrTJeGiWrpkSiVir446oJXPv6HlIL6rpsH+ilJtLfk98uH8O548M7vCfLMh/sKWTZhPBhWedfrMgVhCFoT16NwwEfYFSoDz9lVrhwRIMr3E/LhVOjeGtHHvNGBfOL6TE8uT6zS8D391SzfEI4rWYrKqVESnxQl2NJksSVKbHoDKaBGr5bEkFfENxIUV2Lw9sGe2uYPyqYx7484sIRDa6S+hY2HK3gf6tncO74MPbl13bp3BUb5MkXdy0g2Mc2nWMwWbpdsKVRKQj28cBksaKQJJSKkVfoTQR9QXATrWYLb+/ovRXhSY9dNJG1qcW0mIb3fP7xyiYe/vQgV86K5ccj5V1WJq+eHd8e8IFeV+iaLVYMJgtGs7XDfiOFCPqC4Cb+/WMOJ6qae93Ox0PFw+eNQ6OUhmVLQm+NkmZjxwtZTbORl7pp1TgzPrBPxzdZZHxHcBE3kacvCG7gp8wKXtma69C2l86Ixluj5M9fDb9+RGclh/DjA4v484UTCPLWtL/eU7n9T9I6Vh8trNFT09xxsZrJYsXalr3jqRlZ6a2diaAvCIOsttnIQ59k9LiNRmn7U7397CQunRHD7z87RPkw7EalUkgcq9Bx44JEXrs+hRAfDxaODuHR88fb3T46wJMFo0Pav6/StWIwWwj27jhtYzRbqdMbMZgsmEb4oi0xvSMIg8hqlbn3o/3U6XvOKFk1NZLEYG/mjgrmlndSMZqHZ+DalF3FzhM1bHv4HGbGB/LOTbNIDvMlzU56JsBzV01jZnwgJfUtBHqp2XG8Gj9PFWPCfTts5+2hwttDxd0fpHO0vJFXr0thdJjPQPxIbkfc6QvCIHptWy7bj1f3up2/p5p5o4K54c29VDcN7zo7rWYr+dV6ACZG+aNRKZgY3bXMRLifB1EBnmzKruTCF7aTml/HJdOjWTIunO7WH3lqlORWNfPHLw7T3Gp26c/hrkTQF4RBklfdzDM/5vS6XYiPhitTYvj1R/vRG4d3ps5JnS9sDXpT+xy/r4eKv/1iMnv+sJSoAE+KaluobTZyuLQBsK283Xmihn35tV2Oa227GOzKreH/vjg8IlfpiukdQRgErWYLf/jsUK8lkyUJ/nnFFB765FCfFm0NRb4eKt69eTZatZIof88O77WYLNQ226qMRgZoOX9yBGBbZTs9LgBvjZIv9pdw3dx4fLXqDvP8J1msMrTF+H9fOZVLZ8S49gdyU+JOXxAGkNli5bP0Yi7+7w525db0uv1tZyexJbuaQyUNAzC6wZUU5sPkaH/GR/rh79UxpXJMuC8PrxwHQE5FE4dLGgFbII8J9GLVlChyKpq49vU9vLe7gFY7tYgUkm16RyHBnKRg1/9AbkoEfUEYQN8fKeeBtRlklet63XZsuC+rJkfy/u6CXrcdDjKK6nl7Z363799+dhKxQZ7cMC+e2Ym2MgsqpaKtfr5tzj+juIH//nyMiobW9hTNkyRJ4qEV45AkaUSXYhBBXxAGUKS/Y4W+NEoFf75oAnd/uB/zCJp3rtebun0Iq1BIPHXZFH67Yiya08oml9a38NhXp0pRVDS2csHz2zhS2tjlGFVNrYT7ejC2U3bPSCKCviAMoIhOc9XdeWBZMm9sz6OwVu/iEbmXzLJG0gvru81Qmj8qBL9Oq2m9NMouVTPHRfqSGOrdZf/4YC8WjwvjhyPlzhv0ECOCviAMoCh/LX7anvMnksN8CPfTsuFo5QCNyj0EeWswWay8tPlEe5aNIwK8NPx39fQOr10xMxagyxSPWqngyUsmMd/Og96RQgR9QRhAXxwoodHQc374tXPjWbOvqMdthqPaZiPzR4Xw3FXTCPM9dee+KauSFzYe41hF989BRof5csm0qPbvD5bU09hiQmGniqYkSV0+LYwkIugLwgBpNVt4+vvsHreJC/JiTIQPe/K65piPBJ+kFeHVqTbO6DAfxkf6UanrflGav6eaB5aNbf/+/d2FlDU4XqZ6JBFBXxAGyBf7S3rNtX945The2eJY4TV3t2xCeO8bdZIY4k1rpxITsUFeLJ0Qbjf3/nRxwV5senAxaqXt7n6kLGTrK7E4SxAGyDcHy3p8f0ZcAEoFbM6uGqARuZb3GVSzfHDF2F7r4fckIdiLh1eOo7iuhQWjRu68fU9E0BeEAdBqtnRbNOykXy9J5vmfh099/M537I7ILG1kXETXOju9MZgsZJfrmBobwC1nJfV5/5FETO8IwgB4a0d+r9MNSWHeHCnpmls+FCkVEpEOpqeerqbJeEbn06qVTI0NOKN9RxoR9AXBxQpqmnnBgQ5XBqOl11o8Q4GfVsWN8xP4YE/fVxKvTS0akUXQBpII+oLgQlarzIPrMrq0/7OntMFAmO/Q79l6/uRIPtxTeEbTO8cqm3j6h6xuV+UK/SeCviC4gCzL7Dhezd++Pcq+/J7n8k8qqG4mNsjLxSNzvWOVTYT24+L15vY8dp6ooWiErUYeKCLoC4ILfJpewjWv7+H17XkO71OrNxHh51htHneWVlDHiol9T9c86Zo58cxLCibSX8tPmRU0jdBmJ64igr4gOFmVrpXHvz7S+4ad1OuNhPkN/ekdgPwaPVNj/M9o3/2FdTz+TSbz//Ezt76byscjcHWyK4mgLwhOZrJYe11IZE9ts5EgL40LRjTwfsqs4MKpUb1vaEdGcQNv78xvX4H7w+GRWxzNFUTQFwQnSy+sO6Mqjk2tZnx7KcY2lJgsznkYe7yqadg2gh8MIugLghNtzq7kN2sOcCZZh+UNhn6tRnUnT182hWBv53xqqW02cveH6aQX1lFaL+rp9Nfwua0QhEG2L7+WO99PP+M8c4tVRmmnKuRQMy8pmCtSYvjFizudcrxgbw0/ZVbwU2YFGpWCK2bGcPPCRJJCfZxy/JFG3OkLghPUNLVyx3tptJjOvMiXn6e6T3Xkne3q2XE8ev54rp8Xj7+nY6WHNcqOIUSS4NELxiNJEjPjA50yrlBfD66aZauPbzRbUUgSIov/zIk7fUHoJ1mW+b8vDlPTfGYlBE5SShINLYPXu/WjvYUArJwYwZrb5nLZSzu7LR2hVEj4aVVMiPIjyNuDrLJGjlU2sXB0CJOi/ZFlmW3Hei8cp1ZKjI/0Y0qMrSG6j4eK0noDB4vrKWswkFFcT1a5jjsXjyK3upl5ScH8esloVEpxv3qmRNAXhH76NL2E75yQYRIT6Em9fvAbdn9/pJxpcQH89eJJ/HZdht1tLFaZOr2JHcdrAPjgljn8+qP9/OH88e3bNLd2vWCE+GhYOj68/ZPAykkReKmVSJKtuUlnO49Xc+Pb+/jDZ4d4YfV0zhkbZnc7wXGSuy13TklJkVNTUwd7GILgkOqmVpb+e4tTgvWvFiSgkCTe6MOCLleJ8NOy6cHFtJgsVOoMFNW2kF3eSHWTkbd35nfYdlSoNz/dvwizVe7QsPy7Q2Xc+UF6+/eJId58eud8gvr4gPdwSQN//SaTKl0ro8N8SEkIZPHYMMaM4Obm9kiSlCbLckpv24nPSIJwhk7W1XHW3fmmrEounhaFr8fgfwAvbzTw1s48grw1jIvwY9mEcO5Zksz9y8YQ4NVxvr+6ycjOEzUdAj7AykkRTIg8VSZ5+YTwPgd8gEnR/nxwyxy+u+8sHr94ErtO1LD8P1u54c29bM2pornVjKEfz1JGGoeDviRJSkmS9kuS9E3b94mSJO2RJOmYJEkfS5LU5V9TkiS1JEnvSJJ0SJKko5IkPeLMwQvCYHpx83GnNjzJr9Hz/MZjvHZDCqp+ZvFoVAounxnDv6+cyl8vnsgj541j1ZRIPFSO3+edqGzu8prBZCHUx6PD+BpaTPzl6yPojR3LJUiSxMPnjWN2QhCS1L/mMCqlAg+Vkgh/LW/eOIuPb5tLVICW69/cy8wnfmJdWvEZH3ukcXh6R5KkB4AUwE+W5VWSJK0FPpNleY0kSS8DGbIsv9Rpn9XARbIsXyVJkheQCSyWZTm/u/OI6R1hKEjNr+WXr+52SRngWxYmkhzuw39+OkZ5Y8/tFTtbOj6M+5aOITncBw/VqZx/q1XGIstUN7Xy0CcH2XasusfjRAd4sv7ehQTYWSFsNFvJqdCRVa7DU63ER6siJtCTpBDvbufbm1vNNBpMZ1RjvyffHy5j14karpodx/jIvjdfGU4cnd5xKOhLkhQDvAM8CTwAXAhUARGyLJslSZoH/FmW5RWd9rsaWA38AvAHdgFzZVnutuuzCPqCu9MZTJz//DaKal2zUEitlNj428UcKm5gT14NBpOFn7Mqqe6lwci9S0bzwHJbc/DKRgMnqpqZNyq4y3Ymi5UnvsmkocVEcrgvCcHevLkjr72zl0apYO0d85gmmpIMKY4GfUcnD58FHgJOPjkJBuplWT75ea4YiLaz3yfAxUAZ4AXc31PAF4Sh4LEvj7gs4IOtfMGPR8qZkxjMvWv2Y7HKzEoIpKHF1G1pg9sXJbUHfIBP0ot5Z2c+P963CP9Oc/BqpYK/XDypw2spCYE0t5rJr2lGo1QOaMDfeLSCjKJ67l4yusOnE8E1eg36kiStAiplWU6TJGnxyZftbGrvf+NswAJEAYHANkmSNsiynNvpHLcBtwHExcU5PnpBGGDfHy7js/0lLj+P3mjhw72F7dNH+/LrSIkPRKNSkFZQ16FBybVz43h4xbgO+y8dH85lM2I6BHyrVaa0oYVWs5VRnVazhreVdHbmKleLVUaW5W5z6k0WK+/szOfJb48iy7bG8Q8sH8N5kyLRG834ah1bICb0jSN3+guAiyRJOh/QAn7Y7vwDJElStd3txwCldvZdDXwvy7IJqJQkaQe25wIdgr4sy68Cr4JteudMfxhBcKXKRgOPfHZoQM7lp1VhtnScT09tm36ZGOnL8kkR5FfrWTQmlEumd/2QHRvoxWvbcrlj0Sg0KgWVjQZ+8eJOShtaUEgSj5w3jpsXJjo15/1kGQmj2UpxnZ53dxVw39LkDs8FZFlm67FqXtuaS3mjgeOVTe3v5VY3c8+H+0kKycEqy8xJDGZKrD+JId4khfhQ2tCCt0ZFVIDWoQuC1SqzO6+GVpOVs5JDxIKuNr0GfVmWHwEeAWi7039QluVrJElaB1wOrAFuAL60s3shsESSpPexTe/MxXbBEIQhxWqVefCTg9QNwOKpUF8Pzpscyfu77feYPVKmY1pcIP/55bRuj+GpUXLh1ChOxvQfMisoaStWZpFlnlh/lG8PlbFsQgSxQZ4sHB1i96FtT2RZpqqplTBfLS9sPMb/Nh/n0hkxKCWJn7MqKalvIdxPy52LR7XvI0kSOoOJ7ce7f5CcW23LGsqv0fNxatda+iqFxLnjw5iTGMxlM2PslowwmCzc+X4am9oyhs5KDuGla2fi4wbpsIOtP7+Bh4E1kiQ9AewH3gCQJOkiIEWW5T8B/wPeAg5jmxJ6S5blg/0bsiAMvLd35rM1x3npmfZolApuX5TE7YtG0WK09LhIy5G578QQ7/av9+TWdHk/vbCe9MJ6AHw8VDy4fAw3zE9w+O5/67Fqbnp7H/NHBVNQo8dgsvLhnsIO27y3K5/bzk7qUEguKqB/GTxmq8wPRyr4MbOCy1Ni7G6jVSu5aFpUe9Dfdqyad3bmc/c5o/t17uGgT0FfluXNwOa2r3Oxzdl33uYr4Ku2r5uAK/o7SEEYTFnljfzjuyyXniM2yJNXr0thfKQfh4obqG5qZcHoELbmVHVpML5yYgQPrbQ9tNUbzegMZsJ8PXoM1jGBPffebWo18+evM1l/qIy7zxnN4rFhvY55YpQfQd6aHtM/Q309aDVb8NKoKG8w8PKWE3x7qKzXYzvCU63sUvANbM8Kth+r5ttDHUtjPLfxGGcnhzL5DDt6DRfis44g9KDVbOG+NQcwWlzXxEMhwb+vnMb4SD++PFDC/R/br8cvSfDi6hksnRCOWqlgS04Vj399hBNVzYT7eZASH8S4CF/GRPhyVnIIXppTf953LEoiNsiTRz8/3ONY9uXX8au39/HT/YsYHdbzQ9273k+nqq27VXeMFpnvDpWzMDmETdmVXUo49IfeaOHdXflMjg7g8W8yaWwxIUlQ02S0W+3UaLby6rZcnvvlNBTDoIT1mRJBXxB68O8fc8gq17n0HKumRDErIYg9uTU8uC6j2wYsj180kfMmRwLw6tYT/O3bU58+KhpbWX+ojPVtd9Fjwn3495XTmBRtu6sN8NJw1aw4Pksvac/H744s2xZTGc3WLqUVwHYnvfNEDaG+HmiUih4viEfLGvntugx8PVQuuXCe/jvoSbC3BoPJQr3eSEOLiUAnNXgZikTBNUGww2SxsjWnilveTcWVfyJjwn1Yd/t89CYz5z23rds6PvctTea+pWMA+CStmAe7qX55OrVS4n+rZ7B8YkT7a7Jsq45ZqTOQVabjaFkjmWWNmCxWksN8aTVbmJsUjLeHikVjQrt08iqtb2Hls1tpNlpcshrZlf72i8msnjN8U8KdvThLEEaUFpOFP315xGUB38dDxX1Lk7l2bjxatZI9R2q6DfhB3hruPmc0ZouVT9OL+b2DaaOxgV4khXrTYrTgqbEFb0mSCPLWtBdS65zuabJYufGtvaTm15H+x2VdjhkV4Ml5kyJZm9Y1q8bdPbE+k2mxAUyIGuHlGsSdviB09eC6DD5xYRGvv1w0kRvmJ3R47WhZI7lVzezLr+XDvYVIwOUzY/jVgkTC/Ty458P9bHEgg+i8SREoFRIBXmpS8+u4bEYMyyeGEx/sbXf70voW/vTlYW5akMj3R8r5Yn8Jv1sxluvmJdjdHmzljp/5MZvJbdNHL2/NHRLNy8N8PXj/ljm9lmWWZZkP9xZSpWvljkWjhkTvYqfW3hlIIugLg0lnMHGgqJ7r3tjr0vOcDPp6o5kdx2tYNiG8w/vFdXrAlnWTXljHg+syyK3qWvWyM61aQbC3R3tO/kmLx4by9q9OJduV1rewO7eGgho9FY0G1uwr6rDtM1dMJdjHo9fzWa0y+TXNLHlmS6/buotgbw2f3TWfCH9tt6mv3x8u54730wCID/biyUsmszA5ZCCH2Weinr4gnAEvjYpXtuT2vmE/PfNjNrIs46lWtgf8/OpmqnStyLJMTKAXvh5qXt+Wy/Vv7HUo4AMYTNYuAR/gYHEDza2nSh9/vr+EB9Zm8NzGYx0CPthKIC98ahNPf59Fpa7nKp8ytlo+vtqhM1Nc02zkt2szug34sizz7Iac9u8LavRc+8Yefrcug5qmnrOVhgIR9AWhjdUqc9u7qT2uFnWWRoOZ9MK69tz6Z37MZvG/NjPryQ3884dsvjtUxopnt/LE+qM0tZp7OVrvRof68OC6jPaHr1N6yVVvMVl4cfMJFj61iX98l0VdN/1/lQqJ2CAvJkUNrdz3tMI6jpQ22H0vq1xnN2NrXVoxi/+1mbd25NFqHrpNW0TQF4Q2X2WUsjGrckDOFeGnZVyE7YFiWkEtL/x8vP29Fzef4M4P0vtcS78n6YV1/JRZ0d7oxNE5aqPZystbTrDgqZ95cn0mld2MyTDEgqAsw6OfH8ZqJwNpbLhvtxdFncHMX77O5PzntpFWMDQLBougL4x4u3Nr2JBZwaOfD0wxNX9PNW/9ahaeaiU1Ta3oDP2/k++N2SqjViq458P9nP30Jla/trtP++uNFl7blseCp37mkc8OkVNx6k7YbLF2mDoaKg4U1fPkt0e7BH6FQuKiqVE97nuiqpnLXtrFb9dmUD3EpnxE0BdGNJ3BREKwN3/79ijNRtffrUb5a/nsrvlYrDJz/r6RjOJ6Fo8NY+8fziW5lxWw/dVisrAlp4rCWn23dfl7Y7LIfLS3kCfXH21/rU7ffZ1/d/fG9jx+/dH+Dj12Nx6t4Okfsh3a/9P0Ypb9ewvv7S7A7MJV284kgr4wYukMJn44UsFfv8lsr+zoSsHeGt6/ZQ71ehNXvrKLxGBv5o8KoUFvqzo5EBU8neHsMaE8ffmU9u9DfT1YOr73Wj3uav2hMm54c2/71NekaP8+pZ/W6U388YvDXPD8dlLz3X/KRwR9YURqajVTUKMnLsiLYjvZLs6mVEi8sHo63h4q7vogjQunRPHGjSn884dsFv1rEw8MgWkCf081cUFe3LtkdHvTlZMWJod22V6tHDr1bfbk1XLru6k06E1nPN2WXaHj8pd38fAnB6nX99zacjCJPH1hRMou11Hd1MqBonr+6eBH+f741xVTWTUlkitf2YUkSbx302xuey+V3bnuf2cItoqW6+6Y117LpzOrVWZjViWfpRcT4KUmOsCT3KrmAeky5kwapYLFY0NpMVl6bR7fk3A/D+5fOoYrU2IHrLibWJwlCN2o1BnILteRXa7juQ3H0Ln4IeTvVozl5oWJ/Pqj/Ryr0PHJHfN4Yv1Rvjhgr9mc+4oJ9OTb35yFp1qJSiFfrcgHAAAgAElEQVR1W8r57R15/OWbTJfWLHK1EB8P/LSqfk/7TY8L4Jkrpjq1DWV3xOIsQeiG0WxlS3YVT6w/6vKAf/PCRO5clMRv12awN6+WN2+cxcasyiEX8AGK61q4/d00Supa7Ab0umYjd76fxp+/HtoBH6C6qRWjxUpCsBfB/ajIub+wnvOe28ZrW3PdpkCduNMXRgyzxUqLycKPRyp45LNDLq2RD7Y7/DsXJfHatjxe3HyC129IQZbhmtd3D9lsF7BNgZw9JpQVE8OZER+ISiHx/eFyXt+e12t9/aFGo1SQGOLNsUpdtyWvHTU52p+/Xzq52ymy/hLTO4LQSUOLiRc3H+eD3YVOWeXaHUmCpy6bwpUpsby69QTPbzzOh7fOIdTXgwtf2E51k/s+5OsrlUJCIUkuv4AOthAfD6c8aJckuGFeAr9bMRZvJ/frFdM7gtCJv6eaFqPFpQEfaH+A9+ORcp76Pps/rZrAmHBf7ng/fVgFfLAt+hruAR9s0z32WjP2lSzb+i2veHYr2465tudyd0TQF0aMLw+U8O6uApee467Fo/j1ktGkF9bx0KcH+d/q6VyREsMfvzhMRlG9S88tuJYzL27FdS1c98Zebn03lYIa168ROZ0I+sKI0NRq5rkNx1x6jl+mxPK7FWPJrtBx+3tpPPvLaaycFMkLPx9nnQtr8wtD10+ZFSz7z1ae23DMbh0gVxBBXxgRDhbVu3TV7dLx4Tz5i0m0mq28uOkE7940m8Vjw/hifwn//imn9wMII5bRbOU/G3IobXD9IkEQQV8YAXKrmvjXj65bgDUxyo/nr56GSqmgXm/iiV9MYnykH2kFtfz9u6O9H0AQgJI6EfQFwSlqm42MdlExM1+titdvSMFLo0KWZSL8tfhp1Xx3qIxrX99LuJ8Wb437t9oTBl/RAAX9odPuRhDOUGGtniOljU4/rkal4MVrZhDp74ksy5gsMhqVRHOrmad/yKbFZOFgsf1GHYLQ2UDVXhJ3+sKwZrJYiQ3yYtGYrgXB+uvR88dzVluhMUmS0Khsf07fHCylbIDmZ4XhY6BaMYqgLwxrEhDopaG+xblliwO91Fw7N77L61nljTz86SEMJsfS+y6fGYNWLf4MBSitd16ntJ6I6R1hWGs2WnhvVz6bnNwGcWKUP0o71RO/bqupI0kQ5e+Jr1ZFeaOBeju18qfHBXCgqN7hC4QwvOUNQE8HEEFfGOb8PdUsGB3CO05elDUm3Nfu6zXNRn61IIHEYG9+OlpBnd7InMQgRoX6YLRY2ZdXi5+nGoPJQmm9gZIBqOUvDA0Lk0MG5Dwi6AvDnp+n2qnHC/fzwNpNzaq/XzqZP315hD99daT9tcMljUAFKoXEpTOimT8qhANF9eRWDexKTMF9JYV687sVYwfkXCLoC8OaLMtklTk3cycpxIf7l42x+97/Nh3nvd32P1WYrTLphfV8vr8ECYm7zhnFe7sKqGkeXvV4hL576rIpqJ1Q28cR4gmSMGyZLVbyqpupdXLv2fyaZvztfHqobTayKbvnZwfHK5swWWxFyl7dmssvZ8U6dWzC0HP9vHhmJQQN2PlE0BeGLbNVRpIk9E6uqlnWYLDbADvIW8PF06IdPo7eaKGm2UhSiJczhycMITGBnjy8ctyAnlMEfWHYOljcwDn/2szr2/OcfuxrXt9DUa2+y+vzkoL7dJyP9xWRGOLDnYuSWDExnIlRfvg4uc664L7+cekUp9fV74343yUMW1Ni/En9v6Xc9UE6e/Oc24C81Wyl1dw11bL2DObnN2ZVsjGrEq1aQUKwN+eOD2NSlD//+jHb7jmEoSvQS01ckBejQn1YmBwyYBk7pxNBXxi28qqbqWkyku+C/Ge1UiIm0LPL6w0tJpZPCCerXEehnU8CPTGYrGSV68gq17Ehs4Lnr57Gb9cedHnTF8E1NCoFsxOCWJgcwqyEQMZGuMenuMEfgSC4yPhIP5pazVS5YHn7hEg/tOquhdS+PljKpuxKQn08+nX8ZqMFncHCF3cv4NkNOVQ3tZKaX4fZTZprC/b5aVVcMj2aZRPCmZUQZPf/yGATQV8Y1vbk1uCKNtCdsy2MZitZZY2sP1iGVYbShv4vqX99Wy6PXzyRbw+V9bspt+B6Ny9M5OGV49prMLkr9x6dIPRDXbORtalFLjn27MSOQT+9oJbb3t3n1OCcVa4jvbCecRF+zjuo4DL51c1uH/BB3OkLw9S61CKKavVUNDp/asdDpWD+6FMP4CoaDTz57VHKdc5bZKVRKkgK9eZoWSOzEgLJdPICM8H5tuRUoTOY8NU6dwW4s7n/ZUkQzsB7uwt4/ufjHHBBM/LlEyM6PJB7e2c+h0qcG5TvOmcU1U2tfHmglBAfDyZF+yF1re8muBGzVWbXiZrBHkavRNAXhh2LVSa7XOey418/71RJ5aZWM+93U3bhTPhpVVw3N54vD5RS3WQkwk/LjhPVHC5pdMmzCcG5doqgLwgDr8lg5vzJkS459oy4gA4PcXcer0ZncF5K5YQoP4rq9O1ldssbDVQ0GvqdDSQMjO3Hqwd7CL0SQV8YdjbnVPL5/hKnHtNbo2RchC9P/mJy+2uyLDv9PEazleZOefl51XriRamGIeF4ZdOAdcA6U+JBrjDs9GX1rVatYGKUP0kh3iSEeJMU4o2PVsWmrCrC/DwI9fEgPtiLmfGBSKdNqpc3GAj383D6St/0wnpm2ym+ZbJYUUiI1M0hQG+00LdiHANLBH1h2LHXpaozhQS3nJXE3eeMtlsx82TvW3tyq5poNJhJLaglMcTb6aWR99op5pZR1EBckBc+Hkoyy1z3vELoP4ubX5kdDvqSJCmBVKBEluVVkiQlAmuAICAduE6WZWOnfa4BfnfaS1OAGbIsH+j3yAXBDqtVZmtOVY/bzE4MotVs5d5zk7ssi7dYZSxWmZwKHUfLGtEZzKgUkF/Twuo5cYwO88HbQ8W/fszGYpWZEOVHakGdK3+kdifLOkyLDaC0voVKnXtPI4xUeqNlsIfQo77c6f8GOAqcXCnyFPAfWZbXSJL0MnAz8NLpO8iy/AHwAYAkSZOBL0XAF1xJoZB4+6bZXPbSzi7vxQR6cvXsOFZMDCch2BtVW9OK6qZWXtlygi8PlDI+0o99+bV2/3C3Havi2aumMSHSjweXj+W+jw9QrzcRF+TV5zo7/XGgqB6tWsGshED25Q/MBUdwXIvJvWslOfQgV5KkGOAC4PW27yVgCfBJ2ybvAJf0cpirgY/ObJiC4LiZ8YE8uPxUZ6uEYC+umhXLJ3fM58b5CYwO820P+EW1em5+Zx9v78ynUtfKlpyqbu/UjlU28ddvMsmpaOKHIxXMSgiiStc6oAH/JIPJitki46UWuRjupsXo3pVRHb3TfxZ4CDjZDToYqJdl+eQlrRjorXvEL4GL+zxCYURLK6jlpc0nuGxGDOf1IQ3zniXJLJsQQXphHWclh+CtUeHvqUahOPUw9liFjuwKHfpWMyaLY/OwGUUN/O6TDMwWGa1aMahTLAeK6xkV6sPxyqZBG4PQVX2Le7e/7DXoS5K0CqiUZTlNkqTFJ1+2s2m3fzWSJM0B9LIsH+7m/duA2wDi4uJ6G5IwQpQ1tHDH++lcMi2K5jOYJx0b4cvYCF+778myzPeHyyiua+FYpeOll1tMFg4WN/R5LK4gy1DRYGB2QhDphbWI0vvuYTjM6S8ALpIk6XxAi21O/1kgQJIkVdvdfgxQ2sMxrqKHqR1Zll8FXgVISUlx70ffwoCJ9PfkLxdNZGZ8IOF+Wqcf/5yxYax+fY/TjzuQdK1m9ubXEhWgJdTHgww3uSCNZCaLe199ew36siw/AjwC0Han/6Asy9dIkrQOuBxbBs8NwJf29pckSQFcAZztpDELI4grVtbqW808/X0WOZW21MvhoLTeQExg/xdwRfpriQrwRCmBRZZJK3B+7aLhzmAa4kG/Bw8DayRJegLYD7wBIEnSRUCKLMt/atvubKBYluXcfo1UEByw60QNb+7IQ5ZlJkT6cfnMWML9PfBQnWpmsf5QGVaGRp2Uvsgu1zErIZDU/Lru51q7EeSlISnUm/TCOsraegFMjwtw/iBHAINp6E/vtJNleTOwue3rXGC2nW2+Ar7qtM/cfoxRGGGKavX8lFnBJdOjCfLWOLxfToWOe9fsp6rt4eqGo5WkFtTxwS1z2rfZl19Lpc6Al0aJp1pBi5vflfVFQ4uJffl1zEkMYk8fVgrPTgziUElDl/UGlS4oSz0StLj5nL7I9xLcjp9WzfGqJp7+Pguzg/OjGUX1XPnKrvaAf5JSIbWXT2g1W9h0tIKiGj0vb8kdVgH/dNUO1n7x1aqYHO3H3rxau4GqpL4FX61YtN9XLcPpTl8QBoK/l5q/nVbYrDe7c2u4+e197Rk+kmR7SDsp2p+bFiS0b7cpqxIZWJNa7OQRu5dKXSv+nioaWrp/XuGpVuDnqe61D0BSiLd4ONxHouCaILiQ1Spz67upNBsteKgUnDcpgnuWJBPsrcFXq0KlVFDRaOCVLSf4OqN02N7dn05nMDMh0o8mQyMW2TZfX6s3olZKBHlpiAv2IqtcR0ldS6/HGgrt/9xNddMQz9MXBHemUEi8f/McMssaSYkPJCHEG3XbalurVWZvXg1NBjPr0oqdWvfe3Z38fTQaTBTVtjAnMYicCh2NreY+lW4YSb8zZzG6+YIJcRkXhrypsQFcPTuO5HDf9oBfWt/CdW/uwWSR2XrMuY1OhorUgjqK61poMVnYk1dLnd7U54eMhbV6UhIC8fMU94cOc/O2luJfUhhWzBYrT357lDV7i2gxWSirN5Bb7fiK2+Gmv6tD9UYLqfl1TI3xF3P7DnL3OX1xpy8MK2arTH51c3sGxUgO+M6UWdZIkFfXvgNCV4kh3oM9hB6JoC8MKzqDmUAvx3P7BceYLDLxwaJloyPOm+Sa/szOIoK+MKyE+nqwbEI4aqVEgghSTnWyHLXQvUAvNRdMEUFfEAaUl4eKf10xldFhPkT5O79Q20hVIVbo9mp0mA9atbL3DQeRCPrCsCLLMkdKG7hwShS/P28cq6ZGdVhVOibcp/1rHw8V4yP97B1GsKOwVs/kaH9mJwYO9lDcliuqwTqbCPrCsCJJEvpWC2WNBuKDvXlg2Rj+cP54ogM8SQj2Yv6oEMA2DfTQyrHoDL03URdOOVTSwPE+9B8YaaIDPAd7CL0SKZvCsJMU6k1eVTPRAZ6olbbqk1sfOofS+haajWbGhPsSE+jJ2tQiih1YlSp0VNtsJDrAk5J68bvrLNTXY7CH0CsR9IVhZ3ZiEOsPlpFRXI9GqWBuUhBKhYS/l5oALzUJwd5szq4krUA0FT9Tkf5aEfTtqGl27xIMIIK+MAzFBHpx29lJgG26R5Zl9ubVEh3oyXeHytiTV8tlM6JZMTGCt3fmD+5gh6CZcYFu3xJwsJS39SJwZyLoC8OSJElYrTJv7shj5aQIPksvZlduDeG+Wur0Ru54P32whzhkma1WMst6rs4JcOfiUXy0t5B6fc/PTcL9PLhkWjRKhcSJqiaOlukorNU7a7gDSqVw8xoMiKAvDFPPbzxGcZ2e+aNCeODjA6ycFEFaQR178x1vLiJ05OOhIjnMx+EesKE+HlitHXt4KSQ4/aU/rZpAhc7AJ6nF6FrNxAR6Mj7Sl8kx/qw/WObM4Q+Iz/eXcPuiJEaH+Q72ULolgr4w7ORWNTEm3JfN2ZWsTS1GIdGnTlJCVxOj/IgN9KKqqZXDRb3f5YOtz8H4SD/25NXipVGSHOZDXnUzzUYLFqtMkLeGSH8tj3+T2b5PblUzuVVDNzvIbJV5fVse/7hsymAPpVsiZVMYVlrNFnIqdOzNq+VYZRPQ8c5SODPjIvy49exECmocD8h78mpJSbDl9LeYLPhq1YT4euCtsS1eOm9SBBuzKl0y3sH0s5v/TCLoC8OC2WIlr7qZdanFKBUK1qUWjchyymdiVGjvBcI+TS/md+sOcsP8BIeP29BiIsTHlsIoy+CpUdJqstLY9u+iUkhY5eF3Ra7UtVLkxs8kRNAXhrwWo4XyRgNbsiv5vy8Oc8f7aehaRcB3VGGtnoWjQ3rdLre6mc3Zlfxp1XjC/RzLR28ymNsDv7+nmqrTyg5vOFrJhVOjzmzQbq7VjRupiDl9YcjbX1THuzsL0KoVSBJYxHxOn5gsMrlVTYyP9ONoL1k5aQX1jIvw43+rZ/BpejGSJDEtJoA/fnnYbqDbm1/LjLgAfsysYHduDX84bxx//to2h19S30KkvxYPlcKtg6QjfLUqJkX5s2pqJCsmRrRf6NyRCPrCGWtqNbMlu4pAbzWjQ31QKiSCvDVI0sCkrRXX6YkO8ESW4fsj5QNyzqHsqlmxbD9ebXcVcmmDgQlR/pTU6dunX7rj46EiJSGIlIQg8qqbuentfd0G7dT8Oh5YNoYfMyuoaDRw1ew4NCol+TXNLBwdQn51Mw8uH8vfvjvKUJrp8VQrWTUlkitSYhkb4Yu/59DpNSCCvtBn2eU6vsoo4d1dBV3mzROCvbhnSTKXzYh2efD39VCzKbsSi9V2pyXm8Hu2MauSe89N5rkNx6i2091pw9EKrp4dy0d7i3o8zqIxoe1ff7yviLweGtW0mCzMTAjk0J+XYzRb0aqVrJ4T1/6+yWJFrVSgVkrtnwAGm1at4Lq58Xx3uNzuBfKqWbE8vHIcgd5Ds2+DJLvZ5TUlJUVOTU0d7GE4nSzLlDca8NWqsVhl/D3V5FTo8PFQEeUGRZqyy3V8d7iMhhYTEhItJjO+WjVTYvxpNVmpamqlrtnIrtwaDjrQNu/auXH8YnoM02MDULhgwUp+dTM+WhXlDQaufnW3mMN3kKdayYLRwWw4aj/DJMpfy31Lx/DK1hOc6CZ18sHlY7hnSTIAhTV6bnsvlaxyHWC7+AZ5a/DTqokJ9MRPq+bKWbHMjO+9Mmd5g4FN2ZW8sT2P422ZV4NFqZC4Zk4cEvDFgVIaWk4tMPvm1wuZFO0/eIPrhiRJabIsp/S6nQj6riXLMl8cKCG/Ws8rW08gIfHQyrHMTQrmsS+PUNbYwrrb5xNhp+77He+lYbZa+dcVUwlwQTcog8lCQ4uJJ9Yf5dtDZS6ZC3/nptkd7gz76/VtuWSWNbI1p4rmVgsqpSTu8Pso1NeD8ydFcN28BK5/Yw+lp5UOUCslPr9rAbIMF/53e4f94oK8CPP1YOWkCG45K6n9dYPJQlGtnuhAT7w0zpk8OF7ZxH9/PsZPmRU0D2LJh+QwH66bF8+fvzrSnvq75XeLiQ92v5aIIui7Ab3RzEOfHOSbg2VolAqMbSsZfbUqWowWzG3/i+YlBTM2whe1UmJUqA9zk4KJDvTk3V0FqJUS189LcOh8FqvMgaJ69uXXEuilJtxPy5zEYFuqnNmCh0pJemEdCkkiNb+WrHIdPx4p73UO90ydlRzCtXPjOTs5FE9N/xtLfH+4TJRPcKJLZ0Rzx6JRXPfGng4NUmKDPPn+N2dz8zv72J1bS4SflmevmsbcpOABH2O93sjbO/NZu6+ow8VpIM2MD2R0qA8fpxahVkpkPLbcaRc3Z3I06LvfyIeJolo9t72X1p4NYTxt6XrnO9NduTXsyq1p/37p+DDKGgx8dNtc/LRqrFYZSbLVkzFbrO1t6xoNJlpNVtamFvHtoTKyy3XtF5KT4oO9OGdsGJU6A5Ik8VNmBcYBypQYFerD8gnhTpnbP/3jteAcn6WX0KA38Yfzx3P/xwfa72SLalv427dH+e/qGby3q4Br58a7pGRwWUMLvlo1Ph7dh6EALw33LR3DHYtG8cORct7akc+Bonqnj6UnxXV6zh0fBoBaqcBgsjKU2zCLO30ns1hl9uTVcPcH6dT1UmiqJz4eKkaF+bB4TCheGiWTY/yZGR/Ih3sKKW8w8MWBEqp0rUiS5LYpijfOT+CeJaOdkr5mtcqseHZr+ypbwXmWjAvj7nNG88LPx9icXQXAtNgAbjkrkQsmR9La9gC2J0W1eh794jD/Wz0dX23XTJaGFhNXvryLP66awMJk25qAvOpmdhyv5tq58Q6PVZZltuRUsSmrknVpxQNS7XPRmFCCvDU0tpiYFhtAuL+WK1NiXX7evhJ3+oOkoKaZa17f0+/0s6ZWMxlF9WScdlejUkhd7uTdOc/NQ61w2orLH46Ui4DvIj9nVeKlUfL69Sl8caAUP62KZRPCkWXbp0tHer5+tLeQrTlVPP51Jv+8YmqX9/091Xh7KLn2jT3celYiF0yJIjW/lvw+lHUA23gWjw1j8dgwfn1uMu/tKuDrjFJye8gg6q9zxoYyKdqfKTEB1DS3ctWru1k4OsQtEjDOhAj6TlRcp+emt/e5LA53CfhuTt9qIcy3/z1DZVnmv5uOO2FEgj2R/lqun5eASqng8pkx7a/3NitXpWvFS6OkqdXMR3sL+c25ydy3NLnb7X+1IJH0wv28ti2P17blARDsreGxCyeiVva9OECIjwf3LxvD6jlx3LfmAPuL6jCYnDd1OTXGn5sWJrJiYkT7hU+jVPDWjbOGRC/c7oig7yRfZ5Ty128yqdR1zX8eqZyVx7ztWDVHSh2r7DgS+GlVTnv4vmJiOP+4dEqf/q2+PFDCrhM1rNlXxKvXzSTYx4MFo0O4Z8noHp/fLJ8YjlopYbKcunmpaTZSWKtnVKhPt/vZI8syJovtWVeojwcf3TaXvOpmHv38EPvyazucoy9Gh/mwakokl82IITbIq8v7wT4eBLvxaltHiKDvBDkVOn790f7BHobbadAbsVhllP3M0395ywknjWjoC/LW8J9fTqNa18orW0+QU3HmU15PXzaFK2f1fW76tW25ZJXpWDUlkmajmeXxEUyLDej133nXiRq7wTitoI5RoT4cLWvk799lcWVKDKum9FyTx2ix4qHqOO2UGOLNh7fO5ViFjrs/THf4d+OtUbJ0QjjXzo0nJT5wwFaUDxYR9J3gmyHY7GEg5NXoKahpxipDVID2jNLcmlrNHTKbRrraZiM3vLmXGXEB3Ld0DEqFxCtbTpBe6HhGy9TYACSgVt97P9fTs8VOevKSSaiVCiZE2RYo1TS18tT3WTx24US8e8jE+epAqd3X39yexyXTomlqNbM1p4pjFTpWTozoct7TdQ74p0sO92V2YlCvQT86wJO7zhnFyokRQ/7uvS9E0HeCIyW9r1AdiSoaDCT18WN7Z4U1end+Vj1o0gvrueuDdOKDvbh5YSK/PjeZt3fksyWnqtt9fD1UPH35FM6bHIksyz02NjdbrKxNLebfP+Vw/bx47j331Fz9y1ty+TGzgikx/niqbWs/Vk6MaA/4siyz43gNsUGeHRYxnT0mlK3HqruUgMgq13HtG3v4x6WTuTIlhg1HK3lzRx63nT2qy5j25NXy8b4ipsT4c+3ceLsPmb8/XN5eSuLeJaOZFO3PNwfL+Crj1EVHIcFbv5rFmHD37XDlKiLoO4ErcpiHg1q9EVmW+/Vxubyx+8AkQEGNnj99eYQALzXXzInj1rOS+DmrkrWpRTR1Kk3xn19OY+mEcMCWBRMT2HXOGmBzdiV/+/Zo+51yUqd6+95tbRPNFhkDFqbGBPDAsrFYrTKlDS3EBHoxIz6gyye7lZPCeXNHnt26P3vzann088MkhHjZVrHvL20P+gaThdT8Op7bmMO+/DoA8muaqdebWDYhnKmxAe3HeW93AY99eRirbEsZfmD5WACWTQhnQpQfT3+fhVWGG+YnjMiADyLoO4VIJbSvStfKp+klHTJC+mog8rCHg3q9if9tOsEX+0v56Na53LcsmY/3FvHWjjyCfTy4Zk5ce8Dvzslc+62nfVqYGuPPBZMjOVTcQHK4D1q1kn+dlpJ5oqqJnHIdccG2C8gnacWcPzmyQ0CVZZnaJiOvbctFwlbXpvPakmBvDbeencjMuECOVTR1SOW8+Z197Dh+aopPpZBIDvPlkunR7eVLzBYrj3+Tybu7CgCYEOnHI+ePa99HkiTuWDSKMF8Pfv/ZIaeWBhlqRNB3gtARNB/YV3/4/BDjInzPuEBV0BCtZDhYSupbuOh/27l+bjwmq8x7t8whMdi716J3a/YW8o/vs6jvtKAwo7iBtalF/HJWx8qYb2zP49Lp0QR6aXjsqyOE+HowKyGI+5aOsXv8p3/I4khZI0fLdHYXE9a3mNiaU42ExMz4QBJDvNvPtS+vrsO2F06N4pkrp1JUq2f2kxsI8fGgTm9sX+muVSt4/urpXeb9zRYr4yL8mBztz+yEoB5/H8OZ6JzlBNPjAnrfaIQymq1c+tJO3tieh8HU97t2dyxs5e7q9Sae//k4L20+wQsbj/Ua8DdkVvD7zw51CfgndX6mYjBZqNeb2JVbQ0FNM1/ds5BZPQTRQyUNNBrMHC5p7Hb1uMVqK0xY32Lk+vkJ5LR9elYrFSydENa+XVyQFxdNs2X2PL/xGHqjhcJafXvAv2xGDB/dOpfRYR2fJdlKmUj4alUsHR/G/Kd+xmQZ2o1bzpQI+k4QOJQLcQwAo9nKX7/J5In1mV3mmXsTHeBJ9BBd+ThQTnYMs+es5N6nMbYds//wN9zPg6tnx3LBlMgOr/tq1fz+vHFcPC2a6XGBdivEni7MV2t3Hr+za+bYynFHB3gyNtyH2mZbdtFDK05N01w4NZJ5ScEcr2yymzXno1WxL7+2y+sKhYRSIbE5u5Knf8imXm+irH5wCrgNNjG94wzDO63Xad7fXcihkkbevWl2nzoNTYnx7zHTZKT7+6WTyS5vsrueQaXs/T/n+ZMjyalowiLLeGmULBwdwjnjwvq8YMqe45VN/PmrIw49zK9pOpVC+tDKcZTWtxDkrSE+2Iu5SUHszq0l0EuDwWTh8W8yabHzyfGnI+VMifEns7SRCVF+Hd7bX1jHmzvy2z+5FNbq259FjCSi4Fo/6Qwmlv9nK2WDVPZ1KJoY5cfnd3oAwgcAACAASURBVC1Ao3Lsg+bWnCquf3Ovi0c1dCSFehPq48GevFN3tLcvSsJHo+KZn3I6bDsp2o//u2ACMYGe3WbrONOJqiaiAzzbUyktVpmjZY3szaslMcQbtVLB90fK2H6smkaDuf1uHmy1/J+8ZDLzRwd3GevzG4/x0d5CfnpgEc/8mM1bO/K7nNtTrWT+qGAOldRT02zi/qXJ7c1ewPaJM6dCR2GtnuzyRsJ8tfxiRrRblkk+E6Lg2gDZk1srAn4fHSlt5LP0Yq6aHdf7xtjyux89fzxPfnvUxSMbGlZNiepQiA/glS25LBkXxhOXTOJv3x5tz3o6XNLIVa/uJtJfy98uncw5Y8PsHfKMybJMdoWOcRG2u+q/f3uU0noDf71kEjkVOr47XM7Tl03hpoWJ7fssTA7BbLGy/Xg1N761jxAfDdVNRkwW29oBe9N5N8xP4Lq58UjA7tyu0zerZ8eRW93EjhPV+GnVXDw1iqs7/f966vss3tie1+G1YB8Plk8IRwa7K4rNFis5FU2E+w398gsniaDfT3UOrGoUuvr9Z4eIDfJiwegQh7a/9ewkKhoNvN7pj3akmRYbQKivBxnFXVfg/pxVSWOLibdunMXN76S2Pz8J9tbgoVLwn59yHAr66YV1TIvpuc3ltmNVhPh48N3hctbuK2LXI0uQJInyRgOZZY1c/vJOZBnOHRfWZR2L2WJFqbBVy7xvaTI5FTrKGwykF9aTUVxPi8nS5e775HSgLMt8csc89hfW893hMtamFnHxtCjGhPuyJaeK285K4uazkuxOH14+M4ZxEb68uPlEe1/fO95Pw0+rYkpMAKvnxBHup8VDpSCrXMfm7Ep+zqpEb7QgSfDRrXMHpZGMs4mg309FdhonC46558N01t0xj9Fhji2SuW1REgeLG9hr50HdSKBSSPxyVixfHSjtNtMmtaAOtUqBn1bVHvRrmo3MiA/khaunO3SeN7bnsXJiBLMTg9CqlXYDaKS/lov+u4PJ0f68foNtRsFgslB62sPRxBBv/nnFVFrNFr7YX8K6tGIeWDaG0WE+rN1XzG+WJnPf0jF8lVHK39s+xaXm16HpofzC3rxaDhTVU99i4t5zkzl7TCi+Hio+SStm/b0LO7QVbdCbeOjTDDZlVxEX5MWHt8zhihTbg+lZT2xob8PYaDBzqKSBuz7oviubLMMfPjvEt785y6FS0+5MBP1+ig4YuiVWB1ud3sRfvs7kuaumO5SPH+ar5dXrZ/Krt/exvw+1ZoaLMF8P5iUF883BjjVsYgI90Rst1OmNyLLt4enp1V6TQr25a/Eoh4PVvUuSeeyrwyQEezMqzH7K7KhQH969aTYz4gKxtjU2eXHzCfRG24VmWmwA/7x8CkHeGp5cn8l3h8sprmvh64xS/r+9+45vszobPv47GpYs2Zb3XvHKcIaTOAlJCBmMUAghQEkJu0B5gFJ4C4WOp6WFlrfreYCn0EAD9AUedimjjAIpJBBKBiGTJCSxEzt27MR7L43z/iHZeMi25CXLOt/PJx9s6b6lc7B96b7PuK4ZSZauO5V2m50dx6u7hkinJYa5zbkjpeT373/N8apmmtps1LdaMeq03HlONs3tNvLTI3vMEdU2d7DuyW1dBduPVzV3Lbh4Z195n7q7nlRmO1bVzAvbT3Bjt6Eqf+Rx0BdCaIGdwEkp5SohxCTgZSAS2AVcI6XsM9YhhJgJ/AUIAxzAPCnlhBkEX5kbzx8/OExVkxrmGYotR6s447cfsfUnKzwaMw03BfHyzWfw18+KeHFHMSdrW/GzMgNDVlbfxnc2bGX9VXOpbtrfFdB+8q0prJqZyFNbjmG1SxZmRPHUdfmU1LSQlxLBjORvNsbtLanj0yOVGPQasmJDWDGl7y7dyfGh3LEiu8d5vQkhyO9cm++AKXGhZMeGsCgziu3HakiLMnXdwV02J5nUKDPhwXrOnx5PS4e9q1rWp0eqeH7bCcCZ1vjBNdPdZmaVUrLx4GkKK5uJCTUQaQriXNf6/d5J3uwOyXc2bO2RcC0lIrirtkPBMHbQ//qdg1jtDv7jrAy/zcbpzZX+ncAhnMEb4PfAw1LKl4UQTwA3Ao93P0EIoQOex/mBsFcIEQVMqGKn4aYgXr91MZc+/m8V+Ieow+agoKLJ44kyg07LrcsyuXVZJla7g1P1bZTUtlBa08oXRTW8tqvUb5O0hRl1PHjJDP68qaArqHd3uqGdm5/byR8vn0lDq/OKN921ge36Rem0WO18criSoipndtPDpxs5UFZPZWM7//q6oscE8GVzkt0GfYCFmVFUNLZR32Il25VSob88ShqNID48mHtXTqGgoonvLcnoEYinJIQxJeGb5ZOW4G+uyLuv35+XHklSRDC5v3yf83PjuWVZZtcEsUaj4d07ltBmtSOlc6dubK9CJieqW6hubmd2akSfvTPdU0hXNAzvmvN3//yaXcW1/OWauX4Z+D0K+kKIZOBC4EHgLuHs6QrgStchzwK/olfQB84D9kkp9wJIKSdkjtzUKBP3r57O91/sf0xQGdjB8gYWDGGSTK/VkBJpcha8yHT+cf96zXQa22wcKKtn09cVvLmnzG8Kqze02SiububRdbN54J2DbDlaRXqUicvzU3hvfzkHyhqobu7gxmd3MjkulHCTnr2ldWwrrKbMy1VkA5UqPFjeQG6ipevquKiqmauf3k6UOYjoEAMmg44FkyK5tNuSR4tJz9z0CI/fv8Pm4PltxV3fXzI7CY0QRJiCeHNPGXZJj3kIo17b7xDV23vLeHjjEd74/mJ+/c7BHhu0Ls5LJNpsoHN5+kgU5Pnw4GnqW6095hD8hadX+o8A9wKdM25RQJ2UsnN7ZSmQ5Oa8HEAKIT4AYoCXpZR/6H2QEOJm4GaA1FTPlvGNN+dPjyfRYvT6D09x8maz1mA6g0NnLdV7z5/CU1uOs35zAe228b/1/r8+PMLz205wzcI0LpuTTGZMCJEhQSSGG4kwBfHklmP8u6Da7Z2AN74sruXwqUYmx/edSP/qZD25id8M73x40DkmX9pt4cLbe8v44MApnrw23+vJzaqmdqSEdpuDYL2WDruDWSkWDDotz94wn0vXf054r98JKSWVje1UNXXQarXR0GrjRE0Lmw5X0Nxu49kb5vOH97/mhe3O4aK7zs1h3fxU9p+s44ZndpIZG8LUhNARSZAoBIQMUDtgPBu01UKIVUCFlPJLIcSyzofdHOruhloHnAnMA1qAj1wbCD7qcaKUG4AN4Nyc5Xnzxw+tRrAkO4ZXdpb4uil+xRSk5bErZ/c7zDASzAYdd56TzZrZifz09f18Xjj+bzhPNbTxxw8Os2xyDH9aN5u7X93LxoOnMQdp+eVF04gyG3rkhx+qHUU1boP+osyoHmPrN52ZQXF1S1dA7WTQaTwaSuuwOXhvfzkddgdWu4PHPi4gJy60a3x93fyUrgRpOXGhvHHbIuK6pXfYfaKWO1/ew4maFrev//pti7jntb091vCnR5v5qqyeX7x5AIDf/fMQCzOiWJkbx+FTjRRVu38tT8SEGPxyaAc8u9JfDKwWQlwAGHGO6T8ChAshdK6r/WTA3W9gKfCJlLIKQAjxHjAH+MjNsX6trqVjwAIWSl8RJj3/e+OCIWfg9FZalJnnb1zAI/86wp8+9o9C65sPV3LZ+s958tp8TEFa3tpTxs/e+Ir/uSKPoxVNHCof3lDF33aWcI1rUrW7qBBDj8lUjUbw64unk5tooaS2BaNOy/IpMcxMdp9ssKqpnVd3lpATG8o50+K46bmdPVI2A10rdiLNQTxw8fQez2V3S81c1dROSW1rvwEfnInY9pf2LGZ0x0u7EeKbhHFfFNV25eMfjhlJFlqtdupaOvxyw5ZXaRhcV/o/cq3e+Rvw924TufuklOt7HR+BM8CfCXQA7+Oc/H23v/fwtzQMnX76+n5e2nFi8AMVwFl45n9vnN81UTfW7nplD6/vPtkjKIwHa/ISyY4L5Y8fHO7xeKLFyFVnpJERY+bOl/ag0wp+vWY6d7+6d1jvZwnW8+XPzxmwNKG31m8u4Oktx6l2pVi4flE6l89NprimxTkBG2ok1Kgj1KijodWGXifc/h7YHZJTDW2s+K/NzEy2cMGMBO5/+6Db95wSHzrs4a7+pEQGs2pmIlMTwogOCaKpzUZlUzvnTYsfVwWUxiINw4+Bl4UQvwF2A0+73ng1kC+lvE9KWSuEeAj4Aufwz3sDBXx/ZHdIfv7m/q7ybMrg8tMi2HBtvk9z5f/nhVOZlRLOQxuPjKtJ3g8OnObbc1P6fBiV1TuHe1bPSuT2FZk8tqmQ1BHIpVPfauXzwmrO6qeoSHO7Db1W43GeJIDvLppEmFHPrhO1vL7rJM98XsT24zW8dsvCAWvodpJS0tBq443dpfzty1K+NT2e8vq2PncK3Q0n4AsBQVpNn/kevVbwh2/PZE1ekt8O5bjjVdCXUm4GNru+PgbMd3PMP4B/dPv+eZzLNiecNqudDZ8eUwHfC3NSw3n2hvke/fGPpqgQA0tzYhACfvPOITrGSW71VqudfxdWIaUzEOWlhrOnpI4OV0D6x94yVs1M4NN7lhEdYkAjGPY+hee3FbsN+ofKG7jhmS+obeng1qVZ3HlOtpuz+woO0nL1GWlcnp9MWV0r247VuHLYNDI7deDVPbXNHTy2qYA3d5/sulM4croRnUbjNqumJ1IjTSSGG7EE64kKMRAdYiA8WE90qIEp8aFdieAa2qwUVjRx+FQjZoOOWcnhEzILp39OP48Tr3xRwkO9shoq/QvWa9lwbb7PA36npIhgDpU3ctUZqRwoc2aCHA8aWq38zxV5vLOvnI0HT/d5fsvRKg6daiQrRnYF/O53BhoBOo3G4w+yw6f7XiVvPlzBDc980fX6j3x0hPz0CI9zJYFzP8ULN53B0YpGbHY54NxNQ5uVxz4u4MXtJ/rUXLDaJVa7ZwE/OsTA1IRQzsqOYf6kSNKjzFhMnq0MCzPqmZ0aMegHk78bH399fircw18mBbJjQ3j0ytlEj6OJr4KKJpZkR/PuvnJ2jqN8PkcrmkgMD+aGxZP4oqimR56dS2cn8cvVuRh0GuparAgBc1IjuGpBKlJCc4eNRZnRRJj0XPnkdrcBvbfi6ha2FlazMPObfRKPby7scQchJRwoq/cq6INzVVvv8XopJQfKGtj0dQX7XUtDn99eTGXj4IVWOgXrtSRHBDMrJZzcxDBy4kJJsBjJGIEaABOdCvrD4G0VqEA1OzWcZ66f7/EV11gJ0ml44O2DnOq2QzNIp+kaSvGVHcdr2FlUw4wkC/esnIwpyLmUMViv4/2vyln4249otdp5YHUuB+5f2W8++AtnJnB4o2dj3Xe/uodP713eNaGb2Cu9sSVYz8V57rbieKe+1cr6TQX85dNjXY996OZuprtQg46pCWFkxYWQlxxOXmo4GdHmEZ18DiQq6A9DSY3KsDmYvJRwXrhpwbgsVLHx4Gk6VyVGhwRx81kZAPzf974e0fc5e0osCeFG4kKN6HUa7A5Jh82BTiOobbFidzjosDs3HrVZ7Rh0GuZNiiTEoKPD5qCh1cZXZfVsK6ympqWDGUkWOmwOiqtb+hT/7u7WZZm88kWJR1XHzpkW1yOIzkq28Mbuk2g1gusXpfOdeSnEhXmXXNBqd1Bc3UxVUwfbj9Ww5Wglu07UDjgHkR0bwtlT48iOdW5IS4kIJi3KOeaujIzx95foJ9qsdt7b37dGp/KNjBgzv710hs8DvpSSXSdqeWlHCbXNHWTEmEmNMhMfZuT+i6czJT6U2DADBp2WK5/cBjiX6bVZHV4NOfTnVEMbBZVNaDWC3EQLuYlhzE1z5oexORzYHRK7Q9JucxAXZnS7WQqcxb0PnWpgX2k9ze02TEE62qz2fudIdBpBcNDgO2VjQw38aOXkHo9dszCdeEsw2XEh/ZZNrGpq55UvSvjo0GlaOuzO/6+RZjpsDprbbRRWNrGzeOB18dEhBhZnRXFGRhRLc2L63GEoI0+VSxyiz45WcfXT233djHFrfnokT12fT5hx8CGdd/eVU17fyk1LMka8HVVN7Vz31x395lvRCJia4AzCnYG4rtXKVyfr2dBtCGKszE4N543bFo/Ia7VZ7cz7zb9o7GcYMjokiHOmxnHNwrQeKRcG09hm5clPj/HUZ8e7KnR5Iz3KxPIpsZw7LY4Fk6LcVqxSvKfKJY6yIx5MkAWqS+ck8bMLpmLwcG13c4dt1Nbsv/Zl6YAJthzSmYDrQFkDz20t7ve4sbK/tJ66lo4RSeRl1Gu5+7wcHnzvENOTLOSnRRBvCcYSrGdqQiiT40K9Ghdvs9p5bmsR6zcX9lvExR0hnENcF85MYHFmNDGh/pvCYCJQQX+I1NBOXwadhrvPy+HmszK9Om9tfsrgBw3RvM6c7+PI1IQwEi1GDp9uJEinwaDTEhNqoKqxHbNBy8m61hHL3njF/FQum5tMqAd3XP2x2R289mUpj/zraI9J78HEhhq4ZHYS1yxMG5Oi7IpnVNAfgnabnYPDzHky0UyKNvP/rp9H2jjbzJKXEs7ty7PYdqyaouqWHvnbR5IpSEukOYiYUAMRpiDiwpz/DTXqsTucud/To8wkRwQTF2YcsyGNgdIRD8Zmd/DOvnL+9PFRjlX2n4a5u3CTnlUzE1g9K4n8tIgB6+wqvqGC/hB8drRqSGOZE9XCjKhxtwa/k1YjekxSFlU189aeMh7/pIA2q2dLM/VaQbgpCEuwnjCjjtRIE+nRZtKiTKRHmcmOC/XbNLvuWO0O3tx9skcB8YHoNILlU2K5fG4yyybHepWyQRl7E+c3dQx9OciKhEBy0axEHlo7y2+W1KVHm7nznGyEgIrGNhwSWjvsOKQkWK8lNsxIckQwyeHBRIcaMBt0rjS6zvN1rivXiTgm3WFz8PquUv68uWDQ5chCwOLMaFbmxnHRrES/LCYSqFTQ91Jrh71PTvFAdf2idH58/hS/Cfjd/WBFlkeB2+GQfj9E0dxuY/PhStKiTKRGmfi6vNG15FPLnpI6dFoN1U3tvLH7ZFe64/5kxJhZm5/CpbOT+pQrVPyDCvpeen136bjKyugLM5Mt3L48i/Ny40f0dWubO9BoxIhW0eqPp1fq/h7wKxrauGT95x5t0HLHqNewNCeG9GgzZ2XHsCgzakLe5QQSFfS9IKVk/aZCXzfDp+LDjDx93bxRySNu0Gs4erqJWSnuC3Mo3zhZ18rWwmrarHauWpDqNhCX1rbwq38c8DrgG/UaVubGs2pmIkuyo4c8EayMTyroe+HI6aYhXzH5O61GsGJKLLcuyxy1whGmIJ0K+AOw2R0cLG/gmX8X8caek11ZNZ/4pJCfXTCV5IhgdhXXMj3JwocHT/PM50Ue5xEKN+m5eFYiK3PjmZ0agVGvUVf0E5QK+l7wdLPRRPTzC6fy3cWTfN2MgLK/tJ7HPyngeFULNruD0tpWtznlS2tbue2FXV6/fmyogTOzozkrO4bzp8erK/oAoYK+F7YeG/8FtUdaiEHHedPiuG5huq+bEjCa2m384f2veWH7CezDrZDSjRDOfQtLsqJZOjmW2Snhfj9noXhPBX0PtVntAVkD98XvLei3+LUycupbrByvbuZYZRMPbTxCae3whxHNQVriLUZmJoezNCeG5ZNjCQvWqWGbABfwQV9KyZ6SOl7cfgKjXkuESc/c9EgcDsmkaDNhwXpe+aKEF3cUB1Qq5VCDjluWZTJjgGpHyvA5HJKHNh5hw6fHhlWyccGkSM6fHs/irGjCTXpa2u2kRJrQiIm5p0AZuoAP+s9vK+YXbx1w+1xGjNnj7ecTxbLJMdS3WrljRTbLp8T6ujkT3gs7TvDYpoIhnRsdYuCyuUmszU/pm/7YfXZmRQnMoC+l5M6X99Bhc/Dx4Yp+jwu0gL9ufgoPrpnh62YElI8PDVw1qjchYMXkWK6Yn8qyyTF+uTFO8a2ADPpbj1XzyZHKgN9k1V1+WgS/WTNDTeyNMU+zX6ZGmliTl8jaeSkqY6UyLAET9DuLxby0o4T73voK2wiuivBXWo3AoNMQaQ7iv9fOUsUsfGCg/+URJj0XzUrk4rwk5qSGq7F5ZURM+KAvpeTxTwrZW1JHbbOVHUU1vm7SqDLonLspE8KNHDnVyJajVW4/4FbNTOAXq6YRYtBhCtKqgOIjjW19q1rNTLZw1YJULs5LUmvnlRE34YP+QxuP8OjHQ5so8ze3Lsvku4vSeyTC2l9az5VPbusqmXfpnCSW5sSwMldtxhkPKlw1eHPiQlg1M5ELZiSQFeu+Jq2ijIQJHfTrW6y8tKPE180YE/+xNINbzsokLLjnj3RGsoVfrJrGvX/fx7UL07h/da66qh9HZiZbuPu8HJbmxKifizImJvTUv8Wk58lr5074sepLZydx97mTsZj0bgOHViNYNjmGX16kAv54c99F01g2OVb9XJQxM6GDPsDs1AiyJ/jtcmFlk9tqRTa7g/e/KufPmwp4eG3ehP/w80cGnRpiU8bWhA/6AGdPnbibjEKNOjJjQ3D0mqztsDl4bmsxt7+4mwcunk6EWVU2UhRlgo/pd7p5SSafHa1ib2m9r5syYoK0GhZnRfG9JRksyoqmprmD+lYrBp2GDpuD57cV89Rnx/n+8kzOzI72dXMVRRknAiLoW0x6/nJNPmf89iNfN2VELMqM4qG1ecRbvlml8+7+cn7x5lc9jluaE8P/OSdnrJunKMo4FhDDO+AcBokdpeIfYyXcpOe+VdNYf9WcHgEfnJO5SeHBXd/PnxTJE1fPVdv0FUXpISCu9MG5gkXnxxOZs5ItPLpuDqlR7rfgmw06LpmdxJNbjrFufio/+dYUtQ5fUZQ+Aibobz5cSVl9m6+b4bWUyGC+vyyLy+YmD3rV/qOVk7l9RZYK9oqi9Ctggv4nRyp93QSPZUSbuX5xOrmJYWTFhmIJ9iwpF6ACvqIoAwqIoG93SA6WN/i6GYOKMgfx0wumsiYvEZ0ai1cUZRQERNDXagSXz01mb0mdr5vSR0aMmVuWZlLd1MG6+SmEm9R6ekVRRk9ABH2AtCgTUxPCOOTjK/4Ei5FTDW1ICbcvz+IHZ2epXZmKooyZgAn6Z2ZFc/vyLH74yp5h1SIdiiCdhpvOnMTa/BTSo80cOd3Iscomzp+eMKbtUBRFmbBBX0qJzSHpsDm4+unt1DR3UFzdMubtiA4x8Mx35zG9W4HxnLhQcuJUEVNFUcbehA36hZVNPPpxAaW1rew+4Zux/LgwAy9+74y+RasVRVF8ZEIG/Y8OneaHr+yhwU1VorGSEhnMa7csIi7MOPjBiqIoY2TCBf02q53fvHtoVAN+TlwIS7JjSI82k5sYRnqUmcuf+JzCymbAuSLnhZsWqICvKMq4M+GC/omaFo5XNY/a69+zcjK3LcvsU/Ti3TuW8PjmQvaV1vGr1bkkWIL7eQVFURTfmVBBv6CiiZuf2zlqr796VqLbgA/OnbA/PFdltFQUZXybUEH/YHkDRaOwQifcpOfqBWn84OwsVdZOURS/5nHQF0JogZ3ASSnlKiHEJOBlIBLYBVwjpezodU46cAg47Hpom5TylhFot1v6UciiGWrU8fqti8hQK3AURZkAvLnSvxNnAA9zff974GEp5ctCiCeAG4HH3ZxXKKXMG14zPXNGRhShRh2NQ5jENeo1OBxgdTjIigkh3mIkI9rM2nkpKuArijJheBT0hRDJwIXAg8BdwjnGsQK40nXIs8CvcB/0x0yEOYi7zs3h/rcPDnpspDmI/LQILpmdxNz0CGJCDLRa7dQ0d5Ac4T5nvaIoir/z9Er/EeBeoHMbaRRQJ6XsvKQuBZL6OXeSEGI30AD8XEq5pfcBQoibgZsBUlNTPWySe+vmp1Lb3MGjmwqQsu/zwXotd52bw3WL0gnS9cxkaQrSYQqaUNMciqIoPQwa4YQQq4AKKeWXQohlnQ+7OdRNiKUcSJVSVgsh5gJvCiFypZQ9sp5JKTcAGwDy8/PdvY7HjHotd503mbpWK89tLUarEazJSyIp3IhWo+HM7CjmpkUO5y0URVH8lieXtYuB1UKICwAjzjH9R4BwIYTOdbWfDJT1PlFK2Q60u77+UghRCOTgnBAeVfevzuXSOcnotYLcRMvgJyiKogSAQSt1SCl/KqVMllKmA1cAH0sprwI2Ad92HXYd8Fbvc4UQMa5VPwghMoBs4NgItX1AQgjyUsJVwFcURelmOOWZfoxzUrcA5xj/0wBCiNVCiAdcx5wF7BNC7AVeA26RUtYMp8GKoijK0AnpbrbTh/Lz8+XOnaM++qMoijKhCCG+lFLmD3acKsSqKIoSQFTQVxRFCSAq6CuKogQQFfQVRVECiAr6iqIoAUQFfUVRlACigr6iKEoAUUFfURQlgIy7zVlCiEqg2NftGIJooMrXjRglqm/+Z6L2C1Tf+pMmpYwZ7KBxF/T9lRBipye74fyR6pv/maj9AtW34VLDO4qiKAFEBX1FUZQAooL+yNng6waMItU3/zNR+wWqb8OixvQVRVECiLrSVxRFCSAq6HtBCPGKEGKP61+REGJPr+dThRBNQogf9XP+C0KIw0KIr4QQfxVC6Mem5YMbgb7dLoQoEEJIIUT02LR6cCPQr0lCiO1CiKOu1woam5YPrr++CSHmd3t8rxDikn7OXyGE2OX6fXxWCOFJ+dQxMQJ9O9vVtz1CiM+EEFlj24P+jUDftnQ7rkwI8aZXDZBSqn9D+Af8N3Bfr8f+DvwN+FE/51yAs6i8AF4CbvV1P0awb7OBdKAIiPZ1H0awX68CV7i+fsIffmaACdC5vk4AKjq/73a8BigBclzfPwDc6Ot+jETfXM8dAaa6vr4NeMbX/RipvvU6/+/Atd68p7rSHwIhhADW4gzcnY+twVn/90B/50kp35MuwA6cBeXHlWH0bbeUsmjUGzhEQ+mX65wVOEt9efDDBgAAAudJREFUAjwLrBndlnqvd9+klC1SSpvraSPgbuIuCmiXUh5xfb8RuGy02+qtIfYN1+Nhrq8tQNlotnMohtG3zvNDcf5+enWlr4L+0CwBTkspjwIIIcw4awbf78nJrmGda4D3R62FQzesvo1jQ+lXFFDX7Q+xFEga1VYOTY++AQghFgghDgD7cdamtvU6pwrQCyE6NwJ9G0gZk9Z6Zyh9A7gJeE8IUYrzb+13Y9Ja7wy1b50uAT6SUjZ486bjZgxvvBBC/AuId/PUf0op33J9vY5uV4w4A8fDUsom54f3oNYDn0optwyrsV4ao76NuVHsl7snxnS52xD7hpRyO5ArhJgKPCuE+KeUsq3b81IIcQXwsBDCAHwIDBRgRtxo9c3lh8AFUsrtQoh7gIdwfhCMiVHuW6d1wFNeN87XY1r+9g/nB+VpILnbY1twjmUXAXVADXB7P+f/EuftmMbXfRnpvrmOL2KcjekPtV84g34V34yzLgQ+8HV/Buubm2M2AfmDvM55wKu+7s9I9A2IAQq7fZ8KHPR1f0by54bzLrQaMHr73upK33vnAF9LKUs7H5BSLun8WgjxK6BJSvlY7xOFEDcBK4GzpZSOMWirt4bct3FuSP2SUkohxCacQx8vA9cBbzG+9OmbEGISUCKltAkh0oDJOD/cehBCxEopK1xX+j8GHhyjNntqqH2rBSxCiBzpnLM4Fzg0Rm321JB/bi6XA+/I/u8C+qXG9L13Bb1uyQYihHhPCJHo+vYJIA7Y6lpudd9oNHAYhtw3IcQdrvHTZGCfEML7287RM5yf2Y+Bu4QQBTivrp4ehfYNh7u+nQnsdS0FfAO4TUpZBX36do8Q4hCwD3hbSvnxWDXaQ0Pqm3SOg38P+LsQYi/OMf17xrDdnhjOz62/8z2iduQqiqIEEHWlryiKEkBU0FcURQkgKugriqIEEBX0FUVRAogK+oqiKAFEBX1FUZQAooK+oihKAFFBX1EUJYD8f+5eyHS0r6QkAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "df.plot(figsize=(6, 6))\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Initialize folium map object" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
" ], "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "m = folium.Map(location=[40.70, -73.94], zoom_start=10, tiles='CartoDB positron')\n", "m" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Overlay the boundaries of boroughs on map with borough name as popup" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
" ], "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "for _, r in df.iterrows():\n", " #without simplifying the representation of each borough, the map might not be displayed \n", " #sim_geo = gpd.GeoSeries(r['geometry'])\n", " sim_geo = gpd.GeoSeries(r['geometry']).simplify(tolerance=0.001)\n", " geo_j = sim_geo.to_json()\n", " geo_j = folium.GeoJson(data=geo_j,\n", " style_function=lambda x: {'fillColor': 'orange'})\n", " folium.Popup(r['BoroName']).add_to(geo_j)\n", " geo_j.add_to(m)\n", "m" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Add marker showing the area and length of each borough" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
BoroCodeBoroNameShape_LengShape_Areageometrylatlon
05Staten Island330470.0103321.623820e+09(POLYGON ((-74.05050806403247 40.5664220341607...40.580858-74.153369
14Queens896344.0477633.045213e+09(POLYGON ((-73.83668274106707 40.5949466970158...40.707604-73.818485
23Brooklyn741080.5231661.937479e+09(POLYGON ((-73.86706149472118 40.5820879767934...40.644734-73.947677
31Manhattan359299.0964716.364715e+08(POLYGON ((-74.01092841268031 40.6844914725429...40.777276-73.967159
42Bronx464392.9918241.186925e+09(POLYGON ((-73.89680883223774 40.7958084451597...40.852627-73.866524
\n", "
" ], "text/plain": [ " BoroCode BoroName Shape_Leng Shape_Area \\\n", "0 5 Staten Island 330470.010332 1.623820e+09 \n", "1 4 Queens 896344.047763 3.045213e+09 \n", "2 3 Brooklyn 741080.523166 1.937479e+09 \n", "3 1 Manhattan 359299.096471 6.364715e+08 \n", "4 2 Bronx 464392.991824 1.186925e+09 \n", "\n", " geometry lat lon \n", "0 (POLYGON ((-74.05050806403247 40.5664220341607... 40.580858 -74.153369 \n", "1 (POLYGON ((-73.83668274106707 40.5949466970158... 40.707604 -73.818485 \n", "2 (POLYGON ((-73.86706149472118 40.5820879767934... 40.644734 -73.947677 \n", "3 (POLYGON ((-74.01092841268031 40.6844914725429... 40.777276 -73.967159 \n", "4 (POLYGON ((-73.89680883223774 40.7958084451597... 40.852627 -73.866524 " ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df['lat'] = df.centroid.y\n", "df['lon'] = df.centroid.x\n", "df.head()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
" ], "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "for _, r in df.iterrows():\n", " folium.Marker(location=[r['lat'], r['lon']], popup='length: {}
area: {}'.format(r['Shape_Leng'], r['Shape_Area'])).add_to(m)\n", " \n", "m" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.8" } }, "nbformat": 4, "nbformat_minor": 2 }