{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "bd0633c5",
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import seaborn as sb\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"id": "7fafa1f8",
"metadata": {},
"source": [
"**TODO: Update with correct information.**\n",
"\n",
"We performed 2 test experiments:\n",
"- test1 \n",
" - CTR+ (A01)\n",
" - CTR- (A02)\n",
" - dilution 1:10^5 (A03)\n",
"- test2 \n",
" - CTR+ (A01)\n",
" - CTR- (A02)\n",
" - CTR-scramble1 treat1 (A03)\n",
" - CTR-scramble2 treat2 (B01)\n",
" - hEvolvR treat1 (B02)\n",
" - hEvolvR treat2 (B03)"
]
},
{
"cell_type": "markdown",
"id": "3b78e06a-7048-414a-9c9f-62be0451bffe",
"metadata": {
"jp-MarkdownHeadingCollapsed": true,
"tags": []
},
"source": [
"## Load and normalize data"
]
},
{
"cell_type": "code",
"execution_count": 53,
"id": "f2697d3b-503f-4e30-86e0-d20dd759c33c",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Index(['Plate ID', 'Well Name', 'MEASUREMENT SET ID',\n",
" 'Cell: ObjectID (Arosio)', 'Cell: Big Clusters Count_Sum (Arosio)',\n",
" 'Cell: Big Clusters Area_Average (Arosio)',\n",
" 'Cell: Big Clusters Green Average Intensity_Average (Arosio)'],\n",
" dtype='object')\n",
"Index(['Plate ID', 'Well Name', 'MEASUREMENT SET ID',\n",
" 'Cell: ObjectID (Arosio)', 'Cell: Small Clusters Count_Sum (Arosio)',\n",
" 'Cell: Small Cluster Area_Average (Arosio)',\n",
" 'Cell: Small Clusters Green Average Intensity_Average (Arosio)'],\n",
" dtype='object')\n",
"Index(['Plate ID', 'Well Name', 'MEASUREMENT SET ID',\n",
" 'Cell: ObjectID (Arosio)', 'Cell: Cells Count_Sum (Arosio)',\n",
" 'Cell: Cells Area_Average (Arosio)',\n",
" 'Cell: Cells Green Average Intensity_Average (Arosio)'],\n",
" dtype='object')\n"
]
}
],
"source": [
"bc = pd.read_table(\"./220617/[ ID_23352 ] Arosio Big Cluster_SpikeTest3.txt\")\n",
"sc = pd.read_table(\"./220617/[ ID_23357 ] Arosio Small Cluster_SpikeTest3.txt\")\n",
"c = pd.read_table(\"./220617/[ ID_23356 ] Arosio Cells_SpikeTest3.txt\")\n",
"\n",
"print(bc.columns)\n",
"print(sc.columns)\n",
"print(c.columns)"
]
},
{
"cell_type": "code",
"execution_count": 55,
"id": "d6c897cd-2cf0-44ea-8ad7-25a640ade61b",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Plate ID | \n",
" Well Name | \n",
" MEASUREMENT SET ID | \n",
" Cell: ObjectID (Arosio) | \n",
" Cell: Big Clusters Count_Sum (Arosio) | \n",
" Cell: Big Clusters Area_Average (Arosio) | \n",
" Cell: Big Clusters Green Average Intensity_Average (Arosio) | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" 13098 | \n",
" A01 | \n",
" 23352 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 8812.0 | \n",
" 938.629333 | \n",
"
\n",
" \n",
" 1 | \n",
" 13098 | \n",
" A01 | \n",
" 23352 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 116821.0 | \n",
" 813.279724 | \n",
"
\n",
" \n",
" 2 | \n",
" 13098 | \n",
" A01 | \n",
" 23352 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 29865.0 | \n",
" 799.443054 | \n",
"
\n",
" \n",
" 3 | \n",
" 13098 | \n",
" A01 | \n",
" 23352 | \n",
" 2.0 | \n",
" 1.0 | \n",
" 8945.0 | \n",
" 899.276794 | \n",
"
\n",
" \n",
" 4 | \n",
" 13098 | \n",
" A01 | \n",
" 23352 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 44012.0 | \n",
" 809.692810 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Plate ID Well Name MEASUREMENT SET ID Cell: ObjectID (Arosio) \\\n",
"0 13098 A01 23352 1.0 \n",
"1 13098 A01 23352 1.0 \n",
"2 13098 A01 23352 1.0 \n",
"3 13098 A01 23352 2.0 \n",
"4 13098 A01 23352 1.0 \n",
"\n",
" Cell: Big Clusters Count_Sum (Arosio) \\\n",
"0 1.0 \n",
"1 1.0 \n",
"2 1.0 \n",
"3 1.0 \n",
"4 1.0 \n",
"\n",
" Cell: Big Clusters Area_Average (Arosio) \\\n",
"0 8812.0 \n",
"1 116821.0 \n",
"2 29865.0 \n",
"3 8945.0 \n",
"4 44012.0 \n",
"\n",
" Cell: Big Clusters Green Average Intensity_Average (Arosio) \n",
"0 938.629333 \n",
"1 813.279724 \n",
"2 799.443054 \n",
"3 899.276794 \n",
"4 809.692810 "
]
},
"execution_count": 55,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"bc.head()"
]
},
{
"cell_type": "code",
"execution_count": 56,
"id": "cc05424d-729a-438d-96a1-a127028a56ce",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"([13098], [23352], ['A01', 'A02', 'A03', 'B01', 'B02'])\n",
"([13098], [23357], ['A01', 'A02', 'A03', 'B01', 'B02'])\n",
"([13098], [23356], ['A01', 'A02', 'A03', 'B01', 'B02'])\n"
]
},
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" well_name | \n",
" object_ID | \n",
" count_sum | \n",
" area | \n",
" green | \n",
"
\n",
" \n",
" \n",
" \n",
" 0 | \n",
" A01 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 8812.0 | \n",
" 938.629333 | \n",
"
\n",
" \n",
" 1 | \n",
" A01 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 116821.0 | \n",
" 813.279724 | \n",
"
\n",
" \n",
" 2 | \n",
" A01 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 29865.0 | \n",
" 799.443054 | \n",
"
\n",
" \n",
" 3 | \n",
" A01 | \n",
" 2.0 | \n",
" 1.0 | \n",
" 8945.0 | \n",
" 899.276794 | \n",
"
\n",
" \n",
" 4 | \n",
" A01 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 44012.0 | \n",
" 809.692810 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" well_name object_ID count_sum area green\n",
"0 A01 1.0 1.0 8812.0 938.629333\n",
"1 A01 1.0 1.0 116821.0 813.279724\n",
"2 A01 1.0 1.0 29865.0 799.443054\n",
"3 A01 2.0 1.0 8945.0 899.276794\n",
"4 A01 1.0 1.0 44012.0 809.692810"
]
},
"execution_count": 56,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"def norm(df):\n",
" try:\n",
" df.columns = ['plate', 'well_name', 'set_id', 'object_ID', \n",
" 'count_sum', 'area', 'green']\n",
" plate = list(df[\"plate\"].unique())\n",
" set_id = list(df[\"set_id\"].unique())\n",
" well_name = list(df[\"well_name\"].unique())\n",
" df.drop(df.columns[[0, 2]], axis=1, inplace=True)\n",
" return plate, set_id, well_name\n",
" except:\n",
" raise Warning(\"DataFrame already normalized\", )\n",
"print(norm(bc))\n",
"print(norm(sc))\n",
"print(norm(c))\n",
"bc.head()"
]
},
{
"cell_type": "markdown",
"id": "77fe10a0-847f-4faf-8e3e-d3ec4698e6dc",
"metadata": {},
"source": [
"## The total counts"
]
},
{
"cell_type": "code",
"execution_count": 100,
"id": "9a952ea1-192c-4f4c-8ee2-ecaa432cb941",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" Single cells | \n",
" Small cluster | \n",
" Big cluster | \n",
"
\n",
" \n",
" well_name | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" A01 | \n",
" 64420 | \n",
" 267 | \n",
" 17 | \n",
"
\n",
" \n",
" A02 | \n",
" 31713 | \n",
" 1212 | \n",
" 4 | \n",
"
\n",
" \n",
" A03 | \n",
" 31339 | \n",
" 579 | \n",
" 8 | \n",
"
\n",
" \n",
" B01 | \n",
" 35907 | \n",
" 4291 | \n",
" 15 | \n",
"
\n",
" \n",
" B02 | \n",
" 63666 | \n",
" 4772 | \n",
" 16 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Single cells Small cluster Big cluster\n",
"well_name \n",
"A01 64420 267 17\n",
"A02 31713 1212 4\n",
"A03 31339 579 8\n",
"B01 35907 4291 15\n",
"B02 63666 4772 16"
]
},
"execution_count": 100,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"nbc = bc.groupby(\"well_name\")[\"area\"].count()\n",
"nsc = sc.groupby(\"well_name\")[\"area\"].count()\n",
"nc = c.groupby(\"well_name\")[\"area\"].count()\n",
"\n",
"pd.DataFrame(nc.rename(\"Single cells\")).merge(nsc.rename(\"Small cluster\"), \n",
" left_index=True, right_index=True).merge(nbc.rename(\"Big cluster\"), left_index=True, right_index=True,)"
]
},
{
"cell_type": "markdown",
"id": "7ddcb3b2",
"metadata": {},
"source": [
"Ecco i risultati con il primo metodo di clusterizzazione\n",
"Ho provato ad aumentare il size cellulare da 400µm2 ad 800µm2 e di conseguenza gli small clusters sono tra 800 e 8000, con i big cluster >8000.\n",
">Ho applicato il nuovo sizing ad entrambi gli esperimenti.\n"
]
},
{
"cell_type": "markdown",
"id": "ef868f1b-9fa5-4f1d-b1b0-1867c3b8c926",
"metadata": {},
"source": [
"## Big clusters"
]
},
{
"cell_type": "code",
"execution_count": 98,
"id": "d0eda200",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 98,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEHCAYAAABfkmooAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAnC0lEQVR4nO3de3xV9Znv8c8TbikCSkAQiR2cYi1UkRekauvloA1WO22cSanA6HgDPb3Y9DIzbXV62uq0tOOZekN0BhHRUgVLMyXHWhVQp+oRNQhBRS0ctRqgGJKAIA25PeePvXJlJ3sn7L3XZq/v+/XKKzu/ddlPFuw86/L7PT9zd0RERHqTF3YAIiKS/ZQsREQkISULERFJSMlCREQSUrIQEZGEBoYdQDqMHj3aJ0yYEHYYIiJHlA0bNux292PjLcvJZDFhwgQqKyvDDkNE5IhiZn/qaZluQ4mISEJKFiIikpCShYiIJJSTzyziaWpqorq6moaGhrBDSan8/HwKCwsZNGhQ2KGIREpNTQ3HHhv3WXBOikyyqK6uZvjw4UyYMAEzCzuclHB3amtrqa6u5sQTTww7HJHIqKqqoqysjIULFzJlypSww8mIyNyGamhoYNSoUTmTKADMjFGjRuXc1ZJINmtubmbBggW4OwsWLKC5uTnskDIiMskCyKlE0SYXfyeRbFZeXk59fT0AdXV1lJeXhxxRZkQqWYiIHI7a2lqWLFnSfjXf0NDAkiVLqKurCzmy9FOyEBFJ0rp162htbe3S1traytq1a0OKKHOULLp55513OOWUU+Iumz9/Plu2bOnzPjdt2sSjjz7a6zrLli3juuuuA+DHP/4x48ePZ+rUqZx00kmUlpb2631FUq2mpibsEEJVXFxMXl7XP5t5eXkUFxeHFFHmKFn0wZIlS5g8eXKft0smWXT37W9/m02bNrF161Zmz57N+eefH/kPqoSrqqqKWbNmsXnz5rBDCU1BQQHz588nPz8fiHVdnz9/PgUFBSFHln6RTxa33HILp5xyCqeccgq33XYbEOvtcOmllzJp0iRmzZrFgQMHAJgxY0Z7zaknnniCT3/600ybNo0vf/nL7N+/H4CXXnqJz3zmM5x22mmcfvrp7N27lx/+8IesXLmSqVOnsnLlyj7HOHv2bC644AIefPDB1PzSIn0U1R5A8ZSWlrYnh4KCAkpLS0OOKDMinSw2bNjAfffdxwsvvMD69eu55557qK+v58033+RrX/sar7/+OiNGjOCuu+7qst3u3bv5yU9+wtq1a3n55ZcpKirilltuobGxkdmzZ3P77bdTVVXF2rVrOeqoo7jpppuYPXs2mzZtYvbs2f2Kddq0abzxxhup+LVF+iyqPYDiGThwINdffz1mxg033MDAgdEYrhbpZPHss8/yd3/3dxx11FEMGzaM0tJSnnnmGU444QTOOussAC677DKeffbZLtutX7+eLVu2cNZZZzF16lTuv/9+/vSnP/Hmm28ybtw4PvWpTwEwYsSIlP1HcveU7Eekr6LcA6gnp512GqtWrYrMgDyI0Ajuvug+dqH7z+7OzJkzeeihh7q0v/LKK2mLaePGjRQVFaVt/yI96a0H0CWXXBJSVOGLUqkPiPiVxTnnnMNvf/tbDhw4wIcffsh//dd/cc455/Duu+/y/PPPA/Dggw9y9tlnd9nuzDPP5LnnnmPbtm0AfPjhh/zxj3/k5JNPZufOnbz00ksA7Nu3j+bmZoYPH86+ffv6HedvfvMbnnjiCebOndvvfYj0V5R7AEmHSCeLadOmceWVV3L66adzxhlnMH/+fEaOHMnJJ5/MokWLmDRpEvX19Xz1q19t38bMOPbYY1m2bBlz585lypQpfPrTn+aNN95g8ODBrFy5km984xucdtppzJw5k4aGBs477zy2bNnSpwfct956a3vX2eXLl/Pkk09G7kxGskOUewBJB8vFe+FFRUXefaa8119/nUmTJh3Wfk899VQqKiqyrmhfKn43kd40Nzdz2WWXsWPHDo4//niWL18emQe7PcnFqrNmtsHd497vTtuVhZktNbP3zezVbu3fMLM3zOw1M7u5U/v1ZrbNzN40s891ar8waNtmZt9PV7yJzJw5k1NPPTXrEoVIJkS1B1BPojjmJJ3/4suAO4EH2hrM7DzgYuA0dz9oZmOC9snAHOCTwPHAWjP7eLDZImAmUA28ZGYV7p7x4cxr1qxJyX7uu+8+br/99i5tZ511FosWLUrJ/kXSpa0HUK6dTfdV9zEnUbnKSttv6O5/MLMJ3Zq/Cvzc3Q8G67wftF8MrAja3zazbcDpwbJt7v4WgJmtCNY9YmtfXHXVVVx11VVhhyHSL1FPFBB/zEkUeoVl+gH3x4FzzOwFM/tvM/tU0D4eeK/TetVBW0/tIiIZF+UxJ5lOFgOBAuBM4J+Bhy1FEzKY2bVmVmlmlaqhJCLpoKqzmVMNlHvMi0ArMBrYDpzQab3CoK2n9kO4+2J3L3L3Il0qi0g6RHnMSaafyvwWOA94KniAPRjYDVQAD5rZLcQecJ8EvAgYcJKZnUgsScwB/j6VAX39W//Ert2pu4QcO7qARbf9e6/rXH311TzyyCOMGTOGV1999ZDl7s43v/lNHn30UYYOHcqyZcuYNm1aymIU6a9c7C7aF21jTtpuRUVpzEnakoWZPQTMAEabWTXwI2ApsDToTtsIXOGxgR6vmdnDxB5cNwNfd/eWYD/XAY8DA4Cl7v5aKuPctbuOt8fNSN0Odz6dcJUrr7yS6667jssvvzzu8t///vds3bqVrVu38sILL/DVr36VF154IXUxivRDVVUVZWVlLFy4MFI1kborLS2lvLycHTt2qOpsKrj7XHcf5+6D3L3Q3e9190Z3v8zdT3H3ae7+ZKf1f+ruH3P3k939953aH3X3jwfLfpqueDPp3HPP7fVMZPXq1Vx++eWYGWeeeSZ79uxh586dGYxQpCuVKO8Q1TEnkS73ka22b9/OCSd0PKopLCxk+/a4j2pEMkIlyruKYtVZJQsR6VWUu4v2JmrPbpQsstD48eN5772O4SXV1dWMH6/hJRKOKHcXlQ5KFlmopKSEBx54AHdn/fr1HH300YwbNy7ssCSiotxdVDpE48lML8aOLkiqB1Of9pfA3Llzefrpp9m9ezeFhYXceOONNDU1AfCVr3yFz3/+8zz66KNMnDiRoUOHct9996UsPpG+ausuunjxYhobGxk8eHBkuotKB5UozwG5/LtJdjh48CAXXnghLS0tDBgwgMcee4whQ4aEHZakWCglykUkd6xevZoBAwYAMGDAAFavXh1yRJJpShYi0qu23lCNjY0ANDY2qjdUBClZiEiv1BtKQMlCRBJQbygBJQsRSaCtN1R+fj5ApIrnSQclCxFJqLS0tD05RKl4nnSI/DiL67/9dfbW/jll+zt61HH87Nae59N+7733uPzyy9m1axdmxrXXXss3v/nNLuuoRLlkm7bieWVlZZEqnicdIv8vvrf2z3x/4h9Ttr+fb+t9+cCBA/nFL37BtGnT2LdvH9OnT2fmzJlMnjy5fR2VKJds1FY8L2o1kSRGt6EybNy4ce1XCcOHD2fSpEmHVJRViXLJVkoU0aVkEaJ33nmHjRs3csYZZ3RpV4lyEck2ShYh2b9/P1/60pe47bbbGDFiRNjhiIj0SskiBE1NTXzpS1/i0ksvjdurRCXKJVvV1NSEHYKERMkiw9ydefPmMWnSJL7zne/EXUclyiUbVVVVMWvWLDZv3hx2KBKCyPeGOnrUcQl7MPV1f7157rnn+OUvf8mpp57K1KlTAViwYAHvvvsuoBLlkp26z8G9fPlydZ+NmMj/a/c2JiIdzj77bBKVhTczFi3KbFzSs5qamsj3Aoo3B/cll1wSclSSSZFPFiK9qaqqoqysjIULFzJlypSww0m7ioqKQwoENjU1sWXLlvaTnIaGBhYtWsRTTz3FoEGDgFj9qJKSkozHK5mjZCHSgyjeelm7di0bX9lC69BudZ+GjT1k3c3v1gKQdyBWqlzJIrfl9v98kcMQ1VsvrUMLaJj8haTXz9/ySBqjkWyh3lAicbRN+NPQ0ADEbr1owh+JMiULkTg04Y9IV0oWInFowh+RriL/zOK6f7yOXbW7Ura/saPGcucv7uxxeUNDA+eeey4HDx6kubmZWbNmceONN3ZZ5+DBg1x++eVs2LCBUaNGsXLlSiZMmJCyGCWxtgl/2m5FacIfibrIJ4tdtbvYMX1H6na4offFQ4YM4cknn2TYsGE0NTVx9tlnc9FFF3HmmWe2r3PvvfcycuRItm3bxooVK/je977HypUrUxejJKW0tJTy8nJ27NihCX8k8nQbKsPMjGHDhgGx/utNTU2YWZd1Vq9ezRVXXAHArFmzWLduXcKBfJJ6bRP+mJkm/JHIU7IIQUtLC1OnTmXMmDHMnDmz1xLlAwcO5Oijj6a2tjaMUCOvbcKfKAzIE+mNkkUIBgwYwKZNm6iurubFF1/k1VdfDTsk6UXUS32IgJJFqI455hjOO+88HnvssS7tnUuUNzc3s3fvXkaNGhVGiCIigJJFxtXU1LBnzx4A/vKXv7BmzRo+8YlPdFmnpKSE+++/H4BVq1Zx/vnnH/JcQ0Qkk9L2xM7MlgJfAN5391O6LftH4N+BY919t8X+Et4OfB44AFzp7i8H614B/CDY9Cfufn8q4xw7amzCHkx93l8vdu7cyRVXXEFLSwutra1ccsklfOELX+CHP/whRUVFlJSUMG/ePP7hH/6BiRMnUlBQwIoVK1IXoIhIP6Sze8cy4E7ggc6NZnYCcAHwbqfmi4CTgq8zgLuBM8ysAPgRUAQ4sMHMKty9PlVB9jYmIh2mTJnCxo0bD2m/6aab2l/n5+fz61//OpNhiYj0Km3Jwt3/YGYT4iy6FfgusLpT28XAAx7rH7rezI4xs3HADGCNu9cBmNka4ELgoXTFLRJ1B/fV07R9S9Lr2756oPdJv+TIl9FnFmZ2MbDd3au6LRoPvNfp5+qgraf2ePu+1swqzaxS8wSLiKRWxkYZmdlQ4AZit6BSzt0XA4sBioqKNIJNpJ+GDB+Jj5+c/Pp730pjNJItMnll8THgRKDKzN4BCoGXzew4YDtwQqd1C4O2ntpFRCSDMpYs3P0Vdx/j7hPcfQKxW0rT3P3PQAVwucWcCex1953A48AFZjbSzEYSuyp5PFMxi4hITNqShZk9BDwPnGxm1WY2r5fVHwXeArYB9wBfAwgebP8r8FLwdVPbw24REcmcdPaGmptg+YROrx34eg/rLQWWpjS4Tr573XXs2fV+yvZ3zNgx3Hxn4u64LS0tFBUVMX78eB55pOu0lCpRLiLZJvJlNPfsep9Ld6VuPotfJbne7bffzqRJk/jggw8OWaYS5SLZr6amJlJ1w1TuIwTV1dX87ne/Y/78+XGXq0S5SHarqqpi1qxZbN68OexQMkbJIgTf+ta3uPnmmw+ZtrONSpSLZK/m5mYWLFiAu7NgwQKam5vDDikjlCwy7JFHHmHMmDFMnz497FBEpB/Ky8upr49VHKqrq6O8vDzkiDJDySLDnnvuOSoqKpgwYQJz5szhySef5LLLLuuyjkqUi2Sn2tra9nnZARoaGliyZAl1dbnfSVPJIsN+9rOfUV1dzTvvvMOKFSs4//zzWb58eZd1VKJcJDutW7eOlpaWLm0tLS2sXbs2pIgyJ/K9oY4ZOybpHkzJ7q8/VKJcJPsVFxdz9913d2lzd4qLi0OKKHMinyySGRORLjNmzGDGjBmASpSLHAncPbJX+boNJSKSpHXr1h2SLPLy8nQbSkQkqioqKg5JAk1NTYd0lW1qauKpp57i2WefBWK3qkpKSjIWZ6ZEKlnk4iWkBuuJpMfatWvZ+MoWWocWdF0w7NCpkze/GxsHlXcg1itKyeIIlp+fT21tLaNGjcqZhOHu1NbWkp+fH3YoIjmpdWgBDZO/kPT6+VseSbzSESoyyaKwsJDq6mpybRa9/Px8CgsLww4jp0WtBpBIPJFJFoMGDeLEE08MOww5wlRVVVFWVsbChQuZMmVK2OGIhEa9oUR6ENUaQCLxKFmI9CCqNYBE4lGyEIkjyjWAROJRshCJY926dbS2tnZpa21tjcTgK5F4lCxE4iguLj5kvpG8vLxI1AASiUfJQiSOgoIC5s+f3z6GJT8/n/nz51NQUJBgS5HcpGQhceXaeJT+KC0tbU8OBQUFlJaWhhyRSHiULOQQUZxfOJ6BAwdy/fXXY2bccMMNDBwYmWFJIodQspAuNLagq9NOO41Vq1ZpQJ5EnpKFdKGxBYdSqQ8RJQvpRGMLRKQnShbSTmMLRKQnST2xM7NS4N+AMYAFX+7uI9IYm2RYcXExS5Ys6dIWhbEF8Sa5SUauTnIjEk+y3TtuBr7o7q+nMxgJV9vYgrZbUVEZW9DjJDe9yOVJbvIO1PVpXobYsTgufQFJVkg2WexSooiG0tJSysvL2bFjR6TGFmiSm5h4V5Hvv/8+O3fu7DIro5kxbtw4xowZAxyX81efknyyqDSzlcBvgYNtje6urjI5pm1sQVlZmcYWRFBJSckhV0t1dXXMmTOnveMDwJAhQ7jrrrty/qpTOiT7gHsEcAC4APhi8JX8aZgcUTS2QDpT6ROBJK8s3P2qdAci2UVjC6SzqN6elA5JXVmY2cfNbJ2ZvRr8PMXMfpDe0EQkW6j0iSR7G+oe4HqgCcDdNwNz0hWUiGQf3Z6MtmRPD4a6+4tm1rmt16JBZraU2HON9939lKDtfxN73tEI/D/gKnffEyy7HpgHtABl7v540H4hcDswAFji7j9PMmYRSbGo3Z48uK+epu1bkl7f9tWTq92Ik72y2G1mHwMcwMxmATsTbLMMuLBb2xrgFHefAvyR2NUKZjaZ2JXKJ4Nt7jKzAWY2AFgEXARMBuYG64qISAYle2XxdWAx8Akz2w68DVzW2wbu/gczm9Ct7YlOP64HZgWvLwZWuPtB4G0z2wacHizb5u5vAZjZimDd5FO9iEg/DRk+Eh+f/PnpkL1vpTGacCXbG+otoNjMjgLy3H1fCt77amBl8Ho8seTRpjpoA3ivW/sZ8XZmZtcC1wJ89KMfTUF4mXfHHXewbdu2fm07ceJEysrKUhyRiEhMsr2hxprZvcAqd99nZpPNbF5/39TM/oXYM49f9Xcf3bn7YncvcveiqN1XFRFJt2RvQy0D7gP+Jfj5j8SuCu7t6xua2ZXEHnx/1jvqB2wHTui0WmHQRi/tOUdXBiKSrZJ9wD3a3R8GWgHcvZlYr6U+CXo2fRcocfcDnRZVAHPMbIiZnQicBLwIvAScZGYnmtlgYg/BK/r6viIicniSvbL40MxG0dEb6kxgb28bmNlDwAxgtJlVAz8i1vtpCLAm6Ia73t2/4u6vmdnDxB5cNwNfd/eWYD/XAY8T6zq71N1f69uvKCIihyvZZPEdYmf0HzOz54Bj6ejJFJe7z43T3ONtK3f/KfDTOO2PAo8mGackSXM4iCSmcu0dEiaLYKzD/wi+TiY28dGb7t6U5tgkjTSHg0jveiq73rlke9dS7ZDL5doTJgt3bzGzue5+K6BbQDlEcziI9Cxeufba2lrmzp3bPreHu1NXVxeJcu3JPuB+zszuNLNzzGxa21daIxMRyTJRnqc+2WcWU4PvNwbfjdjD7vNTHZCISLaK6jz1kHyyeIRYcmirJOjAB2Y21d03pSMwEZFsE9V56iH521DTga8A44Djgf8JfA64x8y+m6bYRESyTmlpaXtyiNJEUMkmi0Jgmrv/k7v/I7HkMQY4F7gyTbGJSJapqakJO4TQRXUiqGR/yzHAwU4/NwFj3f0vZnawh20kB1VXV/erLIkKHR75qqqqKCsrY+HChZGfAKltIqgo1aFLNln8CnjBzFYHP38ReDCoQqty4SI5rrm5mQULFuDuLFiwgOXLl3c5o+5vxeQj+SQiSokCki9R/q9m9nvgrKDpK+5eGby+NC2RSVYqLCzkjjvuCDsMybDy8nJqa2sB2L17N+Xl5VxyySUhRyWZlPTNtiA5VCZcUURySm1tLUuWLKGxsRGAxsZGlixZQnFxcfuD3iP16kCSl+wDbhGJqHXr1tHU1LW6T1NTUyQGokkHJQsR6dX06dNpaek6I0FLSwtFRUUhRSRhULIQkV5t2LCBvLyufyry8vKorNRd6ShRshCRXhUXFzNo0KAubYMGDYpEiQvpoGQhIr0qKCjgmmuuYfDgwQAMHjyYa665JhIlLqSDkoWIJFRaWsrw4cMBGD58eGRKXEiHaIxTl7gO7qunaXvyYyptXz25OguYJC+YElkiRlcWIpJQeXk5H374IQD79++nvLw85Igk03RlEWFDho/Ex09Ofv29b6UxGslWbYPyGhoaAGhoaDhkUJ7kPiULEXRLrje9zQ6nkh/RodtQItKr4uLiuOMs1HU2WnRlIYJuyfUmyrPDSQddWYhIQlGdHU46KFmISEJRnR1OOuhfXESSEsXZ4aSDrixEJGlKFB2iNh+5koWISB9VVVUxa9YsNm/eHHYoGaNkEZKonZWI5Iru85E3NzeHHVJGKFmEIIpnJSK5ory8nPr6egDq6uoiU/pED7gzrPtZyfLly0PrWZJ3oI78LY/0af2ojFoWiSfKpU+ULDIs3llJGCUT+jf69jiN2pVIi3LpEyWLDMqms5KSkhJKSkoy+p4iR7ri4mKWLFnSpS0qpU/0zKKTdD907u2sRESyX0FBAVdffXX7nB5mxrx583L+FhSkMVmY2VIze9/MXu3UVmBma8xsa/B9ZNBuZnaHmW0zs81mNq3TNlcE6281syvSFW8mHjqrIJtI7nH3sEPIiHReWSwDLuzW9n1gnbufBKwLfga4CDgp+LoWuBtiyQX4EXAGcDrwo7YEk0qZ6grXVpAtPz8fQAXZRI4wtbW1LF26tD1BuDtLly6lrq4u5MjSL23Jwt3/AHQ/ghcD9wev7wf+tlP7Ax6zHjjGzMYBnwPWuHudu9cDazg0AR228vJyamtrAdi9e3dau8KpIJvIkSvKt5Iz/cxirLvvDF7/GRgbvB4PvNdpveqgraf2Q5jZtWZWaWaVfXn2UFtbyz333ENjYyMAjY2N3HPPPWk7U1BBNpEjV5RvJYf2gNtj13Epu9nn7ovdvcjdi/pSv2bdunU0NTV1aWtqakrrmUJbQbYpU6ak7T1EJPWifCs506e1u8xsnLvvDG4zvR+0bwdO6LReYdC2HZjRrf3pVAY0ffr0uJeVRUVFqXybQ6ggW3bRAEVJVmlpKeXl5ezYsSNSt5IznSwqgCuAnwffV3dqv87MVhB7mL03SCiPAws6PdS+ALg+lQFt2LCBAQMG0NLS0t42YMAAKisr+eu//utUvtURpaamJjIJTQMUpS/abiWXlZVF6laypavbl5k9ROyqYDSwi1ivpt8CDwMfBf4EXOLudRbrtHwnsYfXB4Cr3L0y2M/VwA3Bbn/q7vcleu+ioiKvrKxMKs66ujrmzJnTPlAOYpeWK1asiMSlZTxVVVWUlZWxcOFC3SoT6UEunlCZ2QZ3j3tbJW3JIkx9SRYADz/8MIsXL6axsZHBgwdz7bXX5vzQ/Z40Nzdz6aWXsnPnTo4//vhQa1eJZLOoJQuN4CZ2D3L06NEAjB49OjL3IOOJakVNkb6IYuVoJQs67kECkboH2V1PtauiMOBIJFmaz0IiL8oDjkSSFdWrbyULOs4UgEidKXQX5QFHIsmI8tV3NO+3dJOOOSYqKir6dUZeXFwcWunwtgFHbR+GKA04EkmG5rOIiHh/wJuamtiyZUt7YbCGhgYWLVrEU089xaBBg4D+/QFfu3YtG1/bCMf0YaM9sW9hzjMR1QFHIsmI8nwWkUoWa9euZeMrW2gd2u1MedjYQ9bd/G6ssGBspG4//4AfA60zWhOu1ibv6fDvCkZ1wJFIMqJ89R25vwStQwtomPyFpNfvSwmIXNFWuyrX+pCLpEJUr77DP5WVrKREIRJfVCtHR+O3TKNcHMUpIr3r7er7jjvuYNu2bX3e58SJEykrK0tFeGmhK4vDEMVRnCISE7WTRF1Z9FP3UZyqoSQiQFZfHRwOXVn0U1RHcYpINEXuVPjgvnqatm9Jen3bV0/3SW56GsVZXFwciS50IhI9kUsW/bF7926+/OUvt//8wQcfsGfPni7rNDQ0MHfuXEaMGAHAn//8ZxiTySglXdSJQSSCyWLI8JH4+MnJr7/3Lerr69lXX8Po/NgAu48AHxkaZ+WDNTTV1ADQeHAAg8lPQcQSJk0EJRITuWTRX5PHDOIHRfuTXv+ap4ZyII3xSPqpE4NIBz3gFumBOjGIdFCyEIkjyqWoReLRNXWS3t0/gJ9UDkt6/YMtRsPeBnxr8nOc216D8f2JTlItyqWoReJRskjCyJEjGT36pD5t85GtW/mwqSFNEUm6RbkUtUg8ShZJGD16NHfccUeftikrK2P/9v20ntSHEuXbdVcwW0S5FLVIPPrrJNKD0tLS9uQQpVLUIvEoWYj0IKqlqEXi0f9+kV5oIiiRGF1ZiCSgRCGiZCEiIklQshARkYSULEREJKHIPeDOO1BH/pZH+rR+9/ksOlP5ahGJgkgli95G3+7fv59t27YxceJEhg3rXNbjuB63U/nqaNAJgUjEkkVJSQklJSWHtDc3N7dPblRfX8/ixYsT9qlX+epo0AmBSIyeWQCrVq2itrYWiJWiXrVqVcJtVL4693U/IWhubg47JJHQRD5Z1NbW8p//+Z/tP7s7ixcv7rUUtcpXR4NOCEQ6hJIszOzbZvaamb1qZg+ZWb6ZnWhmL5jZNjNbaWaDg3WHBD9vC5ZPSGUsFRUVtLS0dGlrbm6moqKix216K1/dxR7Iezov6S/2pOiXksOmEwKRrjJ+k93MxgNlwGR3/4uZPQzMAT4P3OruK8zsP4B5wN3B93p3n2hmc4B/A2anMJ4+b5NM+ep+lbIe38/tJOU0n4VIV2E9kR0IfMTMmoChwE7gfODvg+X3Az8mliwuDl4DrALuNDNz9+RnFerFOeecw9KlSw9pP/fcc3vcJpny1T09TJcjg+azEOkq47eh3H078O/Au8SSxF5gA7DH3dueIFbTMWfceOC9YNvmYP1R3fdrZteaWaWZVdbU1CQdz4YNG8jL63oY8vLyqKys7HU7la/ObW0nBPn5+QCaz0IiL+PJwsxGErtaOBE4HjgKuPBw9+vui929yN2L+tInvri4mEGDBnVpGzRoUMIzSJWvzn06IRDpEMYD7mLgbXevcfcmoBw4CzjGzNr+4hYC24PX24ETAILlRwO1qQqmoKCAa665hsGDBwMwePBgrrnmmqTOINvKV6v/fW7SCYFIhzCSxbvAmWY21GJPlz8LbAGeAmYF61wBrA5eVwQ/Eyx/MlXPK9qUlpYyevRoIDaFal/OIDWyN7fphEAkJoxnFi8Qe1D9MvBKEMNi4HvAd8xsG7FnEvcGm9wLjAravwN8P9Ux6QxSeqMTAhGwFJ+kZ4WioiJP9IA6HtUAEpEoM7MN7l4Ub1nkR3B3pkQhIhKfkoWIiCSkZBGSvowFEREJm5JFCKqqqpg1axabN28OOxQRkaQoWWSYyl6LyJFIySLDVPZaRI5EShYZpLLXInKkUrLIoKTnwRARyTJKFhlUXFwct8Ktyl6LSLZTssgglb0WkSOVkkWGqey1iByJlCwyTEULReRIpL9UIWgre61aVCJypNCVRUiUKETkSKJkISIiCSlZiIhIQkoWIiKSUE7OlGdmNcCfwo4DGA3sDjuILKFj0UHHooOORYdsOBZ/5e5xH6jmZLLIFmZW2dMUhVGjY9FBx6KDjkWHbD8Wug0lIiIJKVmIiEhCShbptTjsALKIjkUHHYsOOhYdsvpY6JmFiIgkpCsLERFJSMlCREQSUrI4DGb2t2bmZvaJTm1XmNnW4OuKTu0/NbP3zGx/ONGmV7LHwsyGmtnvzOwNM3vNzH4eXtTp0cf/F4+ZWVVwLP7DzAaEE3XqmVmLmW0Kfr+XzewznZZF6nPS12ORlZ8Td9dXP7+AlcAzwI3BzwXAW8H3kcHrkcGyM4FxwP6w4w7zWABDgfOCdQYH21wUdvwh/r8YEXw34DfAnLDjT+Fx2N/p9eeA/07ieOTk56SvxyIbPye6sugnMxsGnA3MA+YEzZ8D1rh7nbvXA2uACwHcfb277wwl2DTry7Fw9wPu/hSAuzcCLwOFIYSdFv34f/FBsM5AYn8UcrXHyQigPngdyc9JJwmPRTZ+TjSfRf9dDDzm7n80s1ozmw6MB97rtE510Jbr+nUszOwY4IvA7ZkKNAP6fCzM7HHgdOD3wKpMBptmHzGzTUA+sauF84P2KH5O+n0ssuVzoiuL/psLrAherwh+jqo+HwszGwg8BNzh7m+lMbZM6/OxcPfPEfsDMoSOPyK54C/uPtXdP0HsyuEBM7OwgwpJv45FNn1OdGXRD2ZWQOxDfaqZOTCA2O2D7wEzOq1aCDyd6fgy6TCOxWJgq7vflpFAM+Bw/l+4e4OZrSZ2ZbImE/Fmkrs/b2ajgWOB7UTsc9JZH49F1nxOdGXRP7OAX7r7X7n7BHc/AXib2D/8BWY20sxGAhcAj4cZaAb0+ViY2U+Ao4FvhRRzuvTpWJjZMDMbB+1nkH8DvBFW8OkU9AwbANQS+38Qtc9Ju2SPRdZ9TsLuJXAkfgFPEXsI1bmtDLgbuBrYFnxd1Wn5zcTuR7YG338c9u8RxrEgdubkwOvApuBrfti/R0jHYizwErAZeBVYCAwM+/dI4fFo6fRvXAX8TadlUfuc9OlYZOPnROU+REQkId2GEhGRhJQsREQkISULERFJSMlCREQSUrIQEZGElCxERCQhJQuRw2RmT5tZUfD6nWB0rkhOUbIQEZGElCxEAmb2z2ZWFry+1cyeDF6fb2a/MrMLzOz5YPKaXwflyPuy/wlm9rqZ3RNMaPOEmX0kWHaNmb0UTI7zGzMbGrQvM7O7zWy9mb1lZjPMbGmwn2Wd9n1YsYkkomQh0uEZ4JzgdREwzMwGBW2bgR8Axe4+DagEvtOP9zgJWOTunwT2AF8K2svd/VPufhqxEg/zOm0zEvg08G2gArgV+CSxgoVTg9teqYhNpEeqOivSYQMw3cxGAAeJTThTRCxZVACTgeeCytKDgef78R5vu/umTu83IXh9SlA47hhgGF0L6/0fd3czewXY5e6vAJjZa8H2hSmKTaRHShYiAXdvMrO3gSuB/0vsauI8YCKx6rFr3P1w5y052Ol1C/CR4PUy4G/dvcrMrqRr2eq2bVq7bd9K7DPckqLYRHqk21AiXT0D/BPwh+D1V4CNwHrgLDObCGBmR5nZx1P4vsOBncFtr0v7uG26YxNRshDp5hlis9Y97+67gAbgGXevIXbF8ZCZbSZ2m+cTKXzf/wW8ADxHH+e0yEBsIipRLiIiienKQkREEtIDbpEUM7NRwLo4iz7r7rWZjkckFXQbSkREEtJtKBERSUjJQkREElKyEBGRhJQsREQkof8PskuwwXGmL3MAAAAASUVORK5CYII=\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sb.boxenplot(y='green', \n",
" x='well_name', data=bc, hue='object_ID')"
]
},
{
"cell_type": "markdown",
"id": "54228a92-fc21-465a-b42b-24b616bdc079",
"metadata": {},
"source": [
"## Small clusters"
]
},
{
"cell_type": "code",
"execution_count": 99,
"id": "ce72872e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 99,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEHCAYAAABvHnsJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAn9UlEQVR4nO3df3wcd33n8ddHWkuy7Njy+kcS/8I2Mc0FLiagC1DSSwA9EkNoHIRDw4M7AlgJtDRugV4Tm9JcoKjQu8elhEBKqvgIHEd+uOIRXQtJvSK0JE1SZBI58Q8ixSGJnR92JNuJLevHaj/3x47slSXZ2rF2Z1Z6Px+PfWj2OzO7nxntzGfmO9/5jrk7IiIiYZRFHYCIiJQuJREREQlNSUREREJTEhERkdCUREREJLRE1AEU27x583zZsmVRhyEiUlK2bt36mrvPP7F8yiWRZcuW0dbWFnUYIiIlxcyeH61c1VkiIhKakoiIiISmJCIiIqEpiYiIxNj+/fujDuGklERERGKqvb2dtWvXsm3btqhDGZOSiIhIDKXTaRobG3F3GhsbSafTUYc0KiUREZEYam5u5sCBAwB0d3fT3NwccUSjUxIREYmZrq4umpqa6O3tBaC3t5empia6u7sjjmwkJRERkZhpbW0lk8kMK8tkMqRSqYgiGpuSiIhIzNTV1VFWNnz3XFZWRl1dXUQRjU1JREQkZpLJJA0NDVRVVQFQVVVFQ0MDyWQy4shGUhIREYmh+vr6Y0kjmUxSX18fcUSjUxIREYmhRCLBhg0bMDM2btxIIhHP/nLjGZWIiLBq1So2b97M/PkjemCPDZ2JiIjEWJwTCCiJiIjIaVASERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREIrWBIxs01mts/Mns4p+x9mtsvMtpnZT8ysJmfcBjPrNLPfmNllOeWrg7JOM7sxp3y5mT0elN9jZhWFWhYRERldIc9Evg+sPqFsC/A2dz8feAbYAGBm5wFXA28N5vmumZWbWTnwHeCDwHnAx4NpAb4J3OLu5wAHgHUFXBYRERlFwZKIu/8r0H1C2T+7+9DT5h8DFgfDa4C73b3P3Z8DOoELg1enu+92937gbmCNmRnwfmBzMP9dwJWFWhYRERldlNdEPgP8LBheBLyYM25PUDZW+VzgYE5CGioflZldZ2ZtZta2f//+CQpfREQiSSJm9mUgDfyoGN/n7ne4e62718a9R0wRkVJS9OeJmNmngA8DH3B3D4r3AktyJlsclDFGeRdQY2aJ4Gwkd3oRESmSop6JmNlq4M+BK9y9J2dUC3C1mVWa2XJgJfDvwK+AlUFLrAqyF99bguTzELA2mP8a4P5iLYeIiGQVsonvj4FHgd8xsz1mtg64DTgD2GJmT5rZ3wG4+3bgXmAH8ADweXcfDM4y/hh4ENgJ3BtMC3AD8EUz6yR7jeTOQi2LiIiMzo7XKE0NtbW13tbWFnUYIiIlxcy2unvtieW6Y11EREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEpFI6XHFIqVNSUQi097eztq1a9m2bVvUoYhISEoiEol0Ok1jYyPuTmNjI+l0OuqQRCQEJRGJRHNzMwcOHACgu7ub5ubmiCMSkTCURKTourq6aGpqore3F4De3l6ampro7u6OODIRyZeSiBRda2srmUxmWFkmkyGVSkUUkYiEpSQiRVdXV0dZ2fCfXllZGXV1dRFFJCJhKYlI0SWTSRoaGqiqqgKgqqqKhoYGkslkxJGJSL6URCQS9fX1x5JGMpmkvr4+4ohEJIyCJREz22Rm+8zs6ZyypJltMbOO4O+coNzM7FYz6zSzbWb2jpx5rgmm7zCza3LK32lmTwXz3GpmVqhlkYmXSCTYsGEDZsbGjRtJJBJRhyQiIRTyTOT7wOoTym4EWt19JdAavAf4ILAyeF0H3A7ZpAPcBLwLuBC4aSjxBNNcmzPfid8lMbdq1So2b97M+eefH3UoIhJSwZKIu/8rcGKbzTXAXcHwXcCVOeU/8KzHgBozOxu4DNji7t3ufgDYAqwOxs1y98fc3YEf5HyWlJD58+dHHYKInIZiXxM5091fDoZfAc4MhhcBL+ZMtycoO1n5nlHKR2Vm15lZm5m1qa8mEZGJE9mF9eAMwov0XXe4e6271+rIV0Rk4hQ7ibwaVEUR/N0XlO8FluRMtzgoO1n54lHKRUSkiIqdRFqAoRZW1wD355R/Mmil9W7gUFDt9SBwqZnNCS6oXwo8GIx73czeHbTK+mTOZ4mISJEUrF2lmf0YuASYZ2Z7yLay+gZwr5mtA54HPhZM/lPgQ0An0AN8GsDdu83sa8Cvgum+6u5DF+v/iGwLsOnAz4KXiIgUkWUvTUwdtbW13tbWFnUYIiIlxcy2unvtieW6Y11EREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCS0SJKImX3BzLab2dNm9mMzqzKz5Wb2uJl1mtk9ZlYRTFsZvO8Mxi/L+ZwNQflvzOyyKJZFRGQqK3oSMbNFwHqg1t3fBpQDVwPfBG5x93OAA8C6YJZ1wIGg/JZgOszsvGC+twKrge+aWXkxl0VEZKqLqjorAUw3swRQDbwMvB/YHIy/C7gyGF4TvCcY/wEzs6D8bnfvc/fngE7gwuKELyIiEEEScfe9wP8EXiCbPA4BW4GD7p4OJtsDLAqGFwEvBvOmg+nn5paPMs8wZnadmbWZWdv+/fsndoFERKawKKqz5pA9i1gOLARmkK2OKhh3v8Pda929dv78+YX8KhGRKSWK6qw64Dl33+/uA0Az8F6gJqjeAlgM7A2G9wJLAILxs4Gu3PJR5hERkSKIIom8ALzbzKqDaxsfAHYADwFrg2muAe4PhluC9wTjf+7uHpRfHbTeWg6sBP69SMsgIiJkL3AXlbs/bmabgV8DaeAJ4A7gn4C7zeyvgrI7g1nuBH5oZp1AN9kWWbj7djO7l2wCSgOfd/fBoi6MiMgUZ9mD+qmjtrbW29raog5DRKSkmNlWd689sXxc1VlmVm9mHWZ2yMxeN7M3zOz1iQ9TRERKyXirs/4G+H1331nIYEREpLSM98L6q0ogIiLFF/d728abRNqC/qs+HlRt1ZtZfUEjkykh7huISJTa29tZu3Yt27ZtizqUMY03icwCeoBLgd8PXh8uVFAyNZTCBiISlXQ6TWNjI+5OY2Mj6XT61DNFYFxJxN0/PcrrM4UOTiavUtlARKLS3NzMgQMHAOju7qa5uTniiEY33tZZbzGzVjN7Onh/vpn9RWFDk8msVDYQkSh0dXXR1NREb28vAL29vTQ1NdHd3R1xZCONtzrr74ENwACAu28juOlPJF+ltIGIRKG1tZVMJjOsLJPJkEqlIopobONNItXufmKXIqp/kFBKaQMRiUJdXR1lZcN3z2VlZdTV1UUU0djGm0ReM7M3Aw5gZmvJduMukrdS2kBEopBMJmloaKCqqgqAqqoqGhoaSCaTEUc20niTyOeB7wHnmtle4E+BPyxUUDK5ldIGIhKV+vr6Y9tEMpmkvj6ed1WMt3XWbnevA+YD57r7Re7+24JGJpNaqWwgIlFJJBJs2LABM2Pjxo0kEkXvL3dcxts660wzuxPY7O5vmNl5ZrbulDOKjGFoAwFivYGIRGnVqlVs3ryZ888/P+pQxjTe6qzvAw+SfRIhwDNkq7REQuvs7ATg2WefjTgSkfiK+9NYx5tE5rn7vUAGjj3rXM/ukND6+vq47bbbALj11lvp6+uLOCIRCWO8SeSImc3leOusdwOHChaVTHo333wzg4PZ45DBwUFuvvnmiCMSkTDGm0S+SPZxtG82s0eAHwDXFywqmdQ6Ozt5+OGHh5U9/PDD7N69O6KIRCSsUyYRMysHLg5evwt8FnhrcNe6SN42bdo0anlTU1ORIxGR03XKJBI8t/zj7p529+3u/rS7DxQhNpmkGhoa8ioXkfgab3XWI2Z2m5n9npm9Y+hV0Mhk0lqxYgUXXXTRsLKLLrqIFStWRBSRiIRl7n7qicweCgaHJjbA3f39hQqsUGpra72trS3qMKa8vr4+Vq9ezeDgIOXl5TzwwANUVlZGHZaIjMHMtrp77Ynl473D6x/JJhAL3jvwupm93d2fnJgQZSqprKzkIx/5CJs3b+ajH/2oEohIiRpvddY7gc8BZ5O94fCzwGXA35vZn+f7pWZWY2abzWyXme00s/eYWdLMtphZR/B3TjCtmdmtZtZpZttyq9HM7Jpg+g4zuybfOCQ66XSaRx55BMi2zNJDqURK03iTyGLgHe7+Z+7+JbJJZQHwn4FPhfjebwEPuPu5wCpgJ3Aj0OruK4HW4D3AB4GVwes64HYAM0sCNwHvAi4EbhpKPBJ/eiiVyOQw3iSyAMi9pXgAONPdj55QfkpmNpts8rkTwN373f0gsAa4K5jsLuDKYHgN8APPegyoMbOzyZ4JbXH3bnc/AGwBVucTi0RDD6USmTzGm0R+BDxuZjeZ2U3AI8D/NbMZwI48v3M5sB/432b2hJk1BZ9zprsPPaPkFeDMYHgR8GLO/HuCsrHKRzCz68yszcza9u/fn2e4MtH0UCqRyWO8XcF/jWxV0sHg9Tl3/6q7H3H3T+T5nQngHcDt7n4BcITjVVdD3+ccbwl22tz9DnevdffauHdmNhXooVQik8d4z0Rw9zZ3/1bwOp02snuAPe7+ePB+M9mk8mpQTUXwd18wfi+wJGf+xUHZWOUSc0MPpTLLNvYzMz2USqREjTuJTBR3fwV40cx+Jyj6ANkqsRZgqIXVNcD9wXAL8Mmglda7gUNBtdeDwKVmNie4oH5pUCYloL6+nmnTpgEwbdo0PZRKpERF9SSg64EfmVkFsBv4NNmEdm/wsKvngY8F0/4U+BDQCfQE0+Lu3Wb2NeBXwXRfdXddmS0RiUSCpUuX0tnZydKlS/VQKpESFcmWG9ygOOLOR7JnJSdO62Sf8T7a52wCRu/NT2Jv5syZw/6KSOkpenWWiIiMX9xblCqJiIjEVHt7O2vXrmXbtvg+eUNJREQkhtLpNI2Njbg7jY2Nse0aSElERCSGSqVrICUREZGYKaWugZRERERippS6BlISERGJmVLqGkhJREQkZoa6BqqqqgKgqqoqtl0DKYlIJFpaWujo6ADgtddeizgakfipr68/ljSSyWRsuwZSEpFIpFIpyAwyY3rVsRYoInJcIpFgw4YNmBkbN26MbddA8YxKpoQ3Lz0bgGf37DvFlCJT06pVq9i8eTNxfoSFzkRERGIszgkElEREROQ0KImIiEhoSiIiIhKakoiIiISmJCIiIqEpiYiISGhKIiIiEpqSiIiIhKYkIiIioSmJiIhIaEoiIiISWmRJxMzKzewJM/vH4P1yM3vczDrN7B4zqwjKK4P3ncH4ZTmfsSEo/42ZXRbRooiITFlRnon8CbAz5/03gVvc/RzgALAuKF8HHAjKbwmmw8zOA64G3gqsBr5rZuVFil1ECmj//v1RhyDjFEkSMbPFwOVAU/DegPcDm4NJ7gKuDIbXBO8Jxn8gmH4NcLe797n7c0AncGFRFkBECqa9vZ21a9eybdu2qEORcYjqTORvgT8Hhp5EPxc46O7p4P0eYFEwvAh4ESAYfyiY/lj5KPMMY2bXmVmbmbXpCEckvtLpNI2Njbg7jY2NpNPpU88kkSp6EjGzDwP73H1rsb7T3e9w91p3r4173/wiU1lzczNdXV1A9rHJzc3NEUckpxLFmch7gSvM7LfA3WSrsb4F1JjZ0JMWFwN7g+G9wBKAYPxsoCu3fJR5RKTEdHV10dTURH9/PwD9/f00NTXR3d0dcWRyMkVPIu6+wd0Xu/syshfGf+7unwAeAtYGk10D3B8MtwTvCcb/3N09KL86aL21HFgJ/HuRFkNEJlhraysDAwPDygYGBkilUhFFJOMRp/tEbgC+aGadZK953BmU3wnMDcq/CNwI4O7bgXuBHcADwOfdfbDoUUveWlpa6Ojo4NkXXubZF17m6NGjtLS0RB2WROyd73wng4PDN+HBwUFqa2sjikjGI9Ik4u6/cPcPB8O73f1Cdz/H3a9y976gvDd4f04wfnfO/F939ze7+++4+8+iWg7JTyqVgsEBVpw9hxVnz2F6RUJHm8LWrVspKxu+SyorK6OtrS2iiGQ84nQmMmWohRisWDiXb/zhFXzjD69gxcK5UYcjMVBXV8e0adOGlU2bNo26urqIIpLxUBIpMrWBFxldMpnk2muvpaKiAoCKigquvfZakslkxJHJySiJFJHawIucXH19PfPmzQNg3rx51NfXRxyRnIqSSBE1Nzdz4MABALq7u9UGXuQEiUSCq666CoA/+IM/IJFInGIOiZqSSJEMtYHv7e0FoLe3V23gRU6QTqe59957Abjnnnt0tl4ClESKpLW1lUwmM6wsk8moVZJIDp2tlx4lkSKpq6sbtfmiWp6IZOlsfXRxb82pJFIkyWSShoYGqqqqAKiqqqKhoUEtT0QCra2tI6qv0un0lD5bL4XWnEoiRVRfX8/MmTMBmDlz5pRseTJ0t/rul7q48fYWbry9hd0vddHR0aG71qe4urq6Ue9Yn6pn66XSmlNJpMiy3X4d/zvVpFIpSPexfMFMPN2HB8Ok+6b0Eadkr4GcuF24OwcPHowmoIiVyvUhJZEiam5u5siRIwAcOXIktj+KQlt+5iy+/l/eNey1/MxZUYclEdu0adOo5U1NTUWOJHqldH1ISaRISulHIdGI+wXUQmtoaMirfDIrpdacSiJF0traOqK+N64/Cim+UriAWmgrVqzgoosuGlZ20UUXsWLFiogiik4pteZUEimSurq6EfW9mUwmlj8KKa5SuYBaDDfddNOx4fLy8mHvp5JkMslnPvMZzAwAM2PdunWxbM2pJFIk7j5lL6bLyZXKBdRiqKysPNYB4/r166msrIw4oviI6/5DSaRIWltbR5yempmqs6Y4XSsbadq0aUyfPp2PfOQjUYcSma6uLjZt2jSsNeemTZti+btQEimSuro6ysvLh5WVl5erOmuK07WykY4ePUpfX1/UYURKF9ZlhKE71nOflaA71kXXymQ0pXRhXf0sF9GaNWv47ne/C2Qvpq5ZsybiiAqvpaVl2NFTR0cHnu7jy//n8WHT7X71dayrg/Xr1x8rq6ur44orriharFEY7VpZXOu+i6GlpeXYEXhLS8uk//+PZeig8zvf+Q7uHutuknQmUkQ/+clPjm0gmUyGn/zkJxFHVHipVIpndmxj4NVnGHj1GZbNcpYnK/CB3mGv5ckKls3yY9M9s2NbLE/dJ9po1VmDg4NTYtlHk0qlqK6ooLqiYsqugyH19fXHHhecTCZj202SzkSKpKuri+9973vDyu644w4uvfTSWB5dTKTlc6v46prlec3zl/c/V6Bo4uWcc84ZtXzlypVFjiQ+3jR3btQhxEIikThWpbVx48bYPqBLZyJF0tLSMuKIM51Oq9PBKW7z5s2jlt93331FjkTiqLy8nOrqas4///yoQxlT0ZOImS0xs4fMbIeZbTezPwnKk2a2xcw6gr9zgnIzs1vNrNPMtpnZO3I+65pg+g4zu6bYy5KPoZuGqqtncO5b30Z19YyII5I4UFcfxw318Px8VxfPd6ln5yFD+464iuJMJA18yd3PA94NfN7MzgNuBFrdfSXQGrwH+CCwMnhdB9wO2aQD3AS8C7gQuGko8cTR0AXCpcuX8+WvfoOly5cPK5epafbs2SPKzIyampriBxOxVCpFpr+fJbNmsWTWLDL9/VP+ukgpKHoScfeX3f3XwfAbwE5gEbAGuCuY7C7gymB4DfADz3oMqDGzs4HLgC3u3u3uB4AtwOriLUl+kskkCxcuHFa2cOHCSX89RE6utbUVgOoZ1VTPqAaydeFTdee5dPZsbrj4Ym64+GKWjpJgp5qjR49y9OjRqMM4qUiviZjZMuAC4HHgTHd/ORj1CnBmMLwIeDFntj1B2Vjlo33PdWbWZmZtUfSU2tLSwvr16zl06BDPP7ebr//ljTz/3G4OHTrE+vXrdco+hR27TubBi2wT3zjeDyDFNdTcOZPJxHofEdnlfjObCfwD8Kfu/npuvZ+7u5lNWGN5d78DuAOgtra26I3wU6kUO3f9hgWLlgJwpHeAeWcvAWDnrt8ApV+tdeL9IEM6Ojrwgd68W1s919WLvT78vpFck+UekkcffZSq6VUsWLIAgFdeeIXknOSUO0Mduh6S6evjm//yLwC8cPAgZUePTtn7RVKpFFWV048Nx3UdRJJEzGwa2QTyI3cf6m3uVTM7291fDqqr9gXle4ElObMvDsr2ApecUP6LQsZ9OhYsWsonrt8wovxH3/7rCKKZeKlUime2P8mbaoZ37bK0GsDIpPPrxuJNsw3opW/vUyPGPX8we/Qe141qvIZ2nM7x4xrHOXTo0JTbcaZSKTK9vSyeOZPMwAAAi2fMYM/hw7HegRbK0G9jyFAjgziuh6InEcuectwJ7HT3/5UzqgW4BvhG8Pf+nPI/NrO7yV5EPxQkmgeBxpyL6ZcCI/fSERjtLu30YGbUhLFv7wt0l5cdO+Iu5SPsN9WUc9P7qgv+PTc/1FPw7yiGVCrFwOAAyUVJetPZDhjnLJxD997uSb/jHG0bGUtHx9TrySCVSjGYznDW/MUAvLJ/T2x/E1GcibwX+K/AU2b2ZFC2kWzyuNfM1gHPAx8Lxv0U+BDQCfQAnwZw924z+xrwq2C6r7p7LLq4TKVSPL1zF7PnZy+kV83JXt7p6RsYMe3MeWcD8OJrr3No/0tA6R9hy+hG7QKGkbWrjk/6HWcqlWLXtm0sqqoC4CyAyspjZyFDFlZWgjtvPPMMAHuD3o4n07oYy1nzF9Nw1RcAaLrvloijGVvRk4i7PwyM1fD5A6NM78Dnx/isTcDoD2aO2Oz5C/ndj/1RXvP8273fLVA0EgepVIqndj7FzDNnApCYmyBBgqMDw1vfVC/Ins091529jnT41cPA5NtxLqqq4vo3Lctrnm8//9uCxBInQ1VZ6YHMseTx8r497D9QFssqLZtqnb3V1tZ6W1vbhHzWyS4m9w0MMnvBwlHmGtuhfS9ROa18zC4vinE0OtYynUpHRwfef2TENZFCeP7gIFYxI++uQYp1NH+y30VvuvdYEhmvw68epipRNeryxv0M5WTrYvDo0WNnIuO1t7eX8unTS3JdjNf69evZuWMXZyaHNzZ9tXsv/+G8c7n11lsjicvMtrp77Ynl8eyMpUSkUim2bd9JRc2C4SOqkySAI6NUX51MYvZ8BoFde7tGjOs/mG1nUOiNJJVK8czTv2bpzMFTT5xjyTRgGvhAfsscxtIZAAfp/e2vTjXpMS8czia3YuxkUqkU7TvaqZhbMXzELEiQoHegN6/PSyQTpEmz89Wdw8r7u/qBeJ+hpFIpdrW3s/CEfp/OBCgvz/v3srC8HPr7eX379mHlLwWPFI7zuhjLqNeHRju293heH1ISOU0VNQs4+5JPFPx7Xv7Fjwr+HUOWzhzkL2oPF+37iuGv2vI7+j9dFXMrOGvNWQX9jlfuf6Wgnz9RFiYSfK6msJ1J/N3BAwX9/EJKpVLs2L6L+bOz10drpi+A6TCYHv5QqnnB+P17DmX/HsreVqckEgOnU4XT158uyg6+7+A+Onq6x7xvYixxOFIpVafzu+gf6C/4Tr6/q5+Ok9xLczL5/i5OZ10MptMF38m/lE7zakdx1kVYJ6vaG/XM41RGOTMZUsztXkmE7JHAE0/tIFOd5w1eZWdAFaR7C1+FQ9UcXge2Pjv+HVNZT7axmpJIOKlUiie2PwE1ec4YtHIearZbMLPhMId5Yu8T+c13MPsnn99FKpVi15NPku+51fzg74mtribaWQBHjnDwySfzmm9oaypWNeeOp3eSnDl8LZ4xbR5MG3nmcSrJmdlWn6/8dniC7j6cXSolkSLLVCfpPe/DUYcxoap2/GOo+fa93scvfpvfNZG42/d6H0vDzFgDmUvy27jjruwX4Xo7OgtYN2bDytJ0Z6hTgPCSM8/iQ7WfLOh3/LTtBwX9/BMpiQT63jjAwN4dUYcxoeyNA5D3saPk6j3Ui3dMrhaMdsjG6GVOxuO0qjl70wXfyXe98QpvdLxWtKo9JREZYcGsSi5ZVoQquiJ6+LXKqEMoeV29vYy/PVxp6CL/2spUKsX2p3ZSU73g1BPnqC5LUl0Ng4OFPbMdimvvsyNbeZ7MwZ5wLUCVRAKVZ8zBF50XdRgTqvLQ7lDzvXC4vOitmQrthcPlvCXEfFWzq8isnGTVWXv1QNPTZYMJFs06N+owJtShN8J1+KEkEijr6Q59DSGushfW86vOCtsFedxvNnwLIZftYPhrCLF1kFDVWQNVVWyb6FgiNrnOt6OhJMLp7Th7+tNUnnizYQH0HdxHdUUiz7u0z8p72a644opQrTrWr1/PM9ufzHu+sFauXFnwO3dP53dxpOcINq2wF6F9wJlRnf+d+yzKf9lOZ10M9vSMuNlwor2UTlNeXZ33uqgh3LLNPiPJWxa+Pe/54mzv67tCzackwuntOLdt33nqCSeAUZwdZ1hjbYgdHR0c7TlC1bT8juZ7BzJMH2MH+ZYQO8Ewwv4uTnY/wJGeI5RX5He2Ntg/OGayKNb9AKezjexqby9ARCMVc/s42LOPh3bdXZTvKpaDPftYxNy851MSOQ0n23Ee6ekhMS2/i7npgT5mjHU0tWhurJ92N9ZO5sQd6muvvcZLL71EJpOhLHgQWSbov23hwoUsWHD8rK5Ub5Qc77ro6Oigp6eHROXwzTDdl6b6hN9Bqa6Lk20jPUeOUFmW38FFXyZD9YyRCXXWSb5rop1WzcWRHhLlFaee+DSkB/upnpH/Wdkiwu1j1AFjAYy249y3bx/9/dm+jsrKyshkshdrKysrmT9/PvPmzQNKd2eRj8997nPs2DGyOfXq1avZuHFjBBFFo6WlhVtvvZWB9AAVldkdS39fP55xLrvsMr785S9HHGHhnPTgIigbas6wePHiY9sHlO42crIz1J6eHqaV53fQOTDYN+JgY0gh1tFYHTAqiRTJ6tWr6ekZ+TCl6upqHnjggaLHE6UrrriCgwcPjiivqamJ9bOkC6Guru7YwUWuqfi7+NrXvsaWLVtGlN9www1cfvnlEURUHKMl1FdeeYV00Klkrrlz57JkyfEHvRYzoaoX34hdf/31fPOb3xy1fKr50pe+xFe+8pVRy6eaL3zhC/pdBB555JFRy7/97W9P6iQyWvXnWAdag4ODsbsuOsnaLsbX5ZdfzvTp04eVTZ8+fVJvHGO5+OKLSSaH91OWTCa5+OKLI4ooOvpdHDdW4pyKCfU973lPXuVRUhIponvuueek76eSH/7whyd9P5Xod5GlhHrcZz/72bzKo6QkUkQ1NTWsWrUKgAsuuICamppoA4rQGWeccezM433vex9nnHFGxBFFR7+L45RQs5LJJNdee+2wsuuuu27EGXwc6MJ6BO677z6uuuqqSGOIiwceeIDVq1dHHUYs6HeRdf3119Pe3s4FF1zAt771rajDiUw6neZDH/oQvb29VFVV8dOf/pREgW/aPBm1zgrEIYmIyMkpoWa1t7dz/fXXc9ttt3H++edHGouSSEBJRERKyf79+5k/f/6pJyywsZKIromIiMRYHBLIyZR8EjGz1Wb2GzPrNLMbo45HRGQqKekkYmblwHeADwLnAR83s8n1UBARkRgr6SQCXAh0uvtud+8H7gbWRByTiMiUUerdniwCXsx5vwd414kTmdl1wHXB28Nm9psixHYy84DXIo4hLrQujtO6OE7r4ri4rIs3jVZY6klkXNz9DuCOqOMYYmZto7VymIq0Lo7TujhO6+K4uK+LUq/O2gssyXm/OCgTEZEiKPUk8itgpZktN7MK4GpgavUlLiISoZKuznL3tJn9MfAgUA5scvftEYc1HrGpWosBrYvjtC6O07o4LtbrYsrdsS4iIhOn1KuzREQkQkoiIiISmpJIgZjZlWbmZnZuTtk1ZtYRvK7JKf+6mb1oZoejibZwxrsezKzazP7JzHaZ2XYz+0Z0URdGnr+JB8ysPVgXfxf0zjApmNmgmT0ZLN+vzex3c8ZNqW0k33URy+3E3fUqwAu4B/glcHPwPgnsDv7OCYbnBOPeDZwNHI467qjWA1ANvC+YpiKY54NRxx/hb2JW8NeAfwCujjr+CVwPh3OGLwP+ZRzrY1JuI/muizhuJzoTKQAzmwlcBKwj2+wYsj+QLe7e7e4HgC3AagB3f8zdX44k2ALKZz24e4+7PwTg2S5sfk32vp9JIcRv4vVgmgTZncVkbQEzCzgQDE+5beQEp1wXcdxOSrqJb4ytAR5w92fMrMvM3snoXbQsiiS64gm1HsysBvh9YDI91i7vdWFmD5LtH+5nwOZiBltg083sSaCK7NnF+4PyqbiNhF4XcdlOdCZSGB8n2xkkwd+PRxhLlPJeD2aWAH4M3OruuwsYW7HlvS7c/TKyO5ZKju9cJoOj7v52dz+X7JnGD8zMog4qIqHWRZy2E52JTDAzS5Ld4P+jmTnZmyAduAG4JGfSxcAvih1fsZzGergD6HD3vy1KoEVwOr8Jd+81s/vJnslsKUa8xeTuj5rZPGA+2S6LLskZPam3kRPluS5is53oTGTirQV+6O5vcvdl7r4EeI7sj+JSM5tjZnOAS8neaT9Z5b0ezOyvgNnAn0YUc6HktS7MbKaZnQ3HjjgvB3ZFFXwhBS3VyoEusr+DqbSNDDPedRG77STq1gmT7QU8RPYCWG7ZeuB24DNAZ/D6dM74vyFb55kJ/v73qJej2OuB7JGWAzuBJ4NXQ9TLEdG6OJNsv3DbgKeBbwOJqJdjAtfHYM7/uB24PGfclNlGwqyLOG4n6vZERERCU3WWiIiEpiQiIiKhKYmIiEhoSiIiIhKakoiIiISmJCIiIqEpiYgUiJn9wsxqg+HfBncji0wqSiIiIhKakojIKZjZfzOz9cHwLWb282D4/Wb2IzO71MweDR4qdF/Q7Xs+n7/MzHaa2d8HDxr6ZzObHoy71sx+FTy06B/MrDoo/76Z3W5mj5nZbjO7xMw2BZ/z/ZzPPq3YRE5FSUTk1H4J/F4wXAvMNLNpQdk24C+AOnd/B9AGfDHEd6wEvuPubwUOAh8Nypvd/T+5+yqyXV2sy5lnDvAe4AtAC3AL8FayHT2+Pag+m4jYRMakXnxFTm0r8E4zmwX0kX0QUC3ZJNICnAc8EvTgXQE8GuI7nnP3J3O+b1kw/Lagw70aYCbDOyT8f+7uZvYU8Kq7PwVgZtuD+RdPUGwiY1ISETkFdx8ws+eATwH/Rvbs433AOWR7493i7qf7zJi+nOFBYHow/H3gSndvN7NPMbx78KF5MifMnyG7bQ9OUGwiY1J1lsj4/BL4M+Bfg+HPAU8AjwHvNbNzAMxshpm9ZQK/9wzg5aD67BN5zlvo2ESURETG6ZdknzL4qLu/CvQCv3T3/WTPUH5sZtvIVhedO4Hf+xXgceAR8nymSBFiE1FX8CIiEp7OREREJDRdWBcpEjObC7SOMuoD7t5V7HhEJoKqs0REJDRVZ4mISGhKIiIiEpqSiIiIhKYkIiIiof1/DWKsqPU0QxgAAAAASUVORK5CYII=\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sb.boxenplot(y='green', x='well_name', data=sc, )"
]
},
{
"cell_type": "code",
"execution_count": 105,
"id": "d8151840",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" well_name | \n",
" object_ID | \n",
" count_sum | \n",
" area | \n",
" green | \n",
"
\n",
" \n",
" \n",
" \n",
" 622 | \n",
" A02 | \n",
" 3.0 | \n",
" 1.0 | \n",
" 3071.0 | \n",
" 6504.530273 | \n",
"
\n",
" \n",
" 623 | \n",
" A02 | \n",
" 4.0 | \n",
" 1.0 | \n",
" 2686.0 | \n",
" 6500.688477 | \n",
"
\n",
" \n",
" 1186 | \n",
" A02 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 1014.0 | \n",
" 6016.974609 | \n",
"
\n",
" \n",
" 1268 | \n",
" A02 | \n",
" 1.0 | \n",
" 1.0 | \n",
" 5082.0 | \n",
" 10465.741211 | \n",
"
\n",
" \n",
" 1346 | \n",
" A02 | \n",
" 2.0 | \n",
" 1.0 | \n",
" 2547.0 | \n",
" 7104.483887 | \n",
"
\n",
" \n",
" 9159 | \n",
" B02 | \n",
" 7.0 | \n",
" 1.0 | \n",
" 1078.0 | \n",
" 13170.574219 | \n",
"
\n",
" \n",
" 10156 | \n",
" B02 | \n",
" 3.0 | \n",
" 1.0 | \n",
" 983.0 | \n",
" 7135.930176 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" well_name object_ID count_sum area green\n",
"622 A02 3.0 1.0 3071.0 6504.530273\n",
"623 A02 4.0 1.0 2686.0 6500.688477\n",
"1186 A02 1.0 1.0 1014.0 6016.974609\n",
"1268 A02 1.0 1.0 5082.0 10465.741211\n",
"1346 A02 2.0 1.0 2547.0 7104.483887\n",
"9159 B02 7.0 1.0 1078.0 13170.574219\n",
"10156 B02 3.0 1.0 983.0 7135.930176"
]
},
"execution_count": 105,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"sc[sc[\"green\"]>6000]"
]
},
{
"cell_type": "code",
"execution_count": 106,
"id": "b47cefb7",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 106,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEHCAYAAABvHnsJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAj7ElEQVR4nO3df3RcZ33n8fdXGsmy5Bh7YsVx7Bgn4Dab0DghWsi22ZKANjikjUHr0tBSEohjyi9vS7tb7FNOCuz6QM/Zsg3lRxzjJimUJrhKo7aQNKMFygKhkcF2fha5ITR2jGNLdmzLUqSRvvvHXEkjaWRrrjX3h+bzOmeOZp65d+Y7j+693/vc+9znmrsjIiISRk3cAYiISHopiYiISGhKIiIiEpqSiIiIhKYkIiIioWXiDiBqS5Ys8VWrVsUdhohIquzateuIuzdPLq+6JLJq1Sq6urriDkNEJFXM7GelynU4S0REQlMSERGR0JREREQkNCURidXhw4fjDkFEzoKSiMRmz549rF+/nr1798YdioiEpCQiscjn82zduhV3Z+vWreTz+bhDEpEQlEQkFu3t7Rw9ehSA3t5e2tvbY45IRMJQEpHI9fT0sH37dgYGBgAYGBhg+/bt9Pb2xhyZiJRLSUQi19nZycjIyISykZERcrlcTBGJSFhKIhK51tZWamomLno1NTW0trbGFJGIhKUkIpHLZrNs2LCBhoYGABoaGtiwYQPZbDbmyESkXEoiEou2traxpJHNZmlra4s5IhEJQ0lEYpHJZNi8eTNmxpYtW8hkqm4sUJE5QWuuxGbNmjXs3LmT5uYpo0uLSEqoJSKxUgIRSTclERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEhGRxDl8+HDcIcgMKYmISKLs2bOH9evXs3fv3rhDkRlQEhGRxMjn82zduhV3Z+vWreTz+bhDkjNQEhGRxGhvb+fo0aMA9Pb20t7eHnNEciZKIiKSCD09PWzfvp2BgQEABgYG2L59O729vTFHJqejJCIiidDZ2cnIyMiEspGREXK5XEwRyUwoiYhIIrS2tlJTM3GTVFNTQ2tra0wRyUxULImY2YVm9i0ze9rMnjKz/xaUZ83sUTPrDv4uDsrNzO40s31mttfMXl/0WbcE03eb2S1F5VeZ2RPBPHeamVXq94hIZWWzWTZs2EBDQwMADQ0NbNiwgWw2G3NkcjqVbInkgT9w90uBq4EPmdmlwMeATndfDXQGrwFuAFYHj43AF6GQdIA7gDcCbwDuGE08wTS3F823toK/R0QqrK2tbSxpZLNZ2traYo5IzqRiScTdD7r7j4LnJ4BngOXAOuDeYLJ7gbcHz9cB93nBY8AiM1sGvBV41N173f0o8CiwNnhvobs/5u4O3Ff0WSKSQplMhs2bN2NmbNmyhUwmE3dIcgaR/IfMbBVwJfBDYKm7Hwze+jmwNHi+HHihaLb9QdnpyveXKC/1/RsptG5YuXLlWfwSEam0NWvWsHPnTpqbm+MORWag4ifWzWwB8LfA77n78eL3ghaEVzoGd9/m7i3u3qIFUyT5tJ6mR0WTiJnVUUggX3X30auGDgWHogj+vhSUHwAuLJp9RVB2uvIVJcpFRCQileydZcCXgWfc/c+K3uoARntY3QI8VFT+nqCX1tXAy8Fhr0eA681scXBC/XrgkeC942Z2dfBd7yn6LBERiUAlz4n8CvA7wBNmtjso2wJ8GnjAzG4Dfga8M3jvG8DbgH3AKeC9AO7ea2afAh4Ppvuku49ewvpB4B5gPvDN4CEiIhGxwmmJ6tHS0uJdXV1xhyEikipmtsvdWyaX64p1EREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0JREREQkNCUREREJTUlERERCUxIREZHQlERERCQ0JREREQlNSUREREJTEhERkdCUREREJDQlERERCU1JREREQlMSERGR0CqWRMxsh5m9ZGZPFpX9iZkdMLPdweNtRe9tNrN9ZvavZvbWovK1Qdk+M/tYUflFZvbDoPx+M6uv1G+Ryjl8+HDcIYjIWahkS+QeYG2J8s+6+xXB4xsAZnYpcDNwWTDPF8ys1sxqgc8DNwCXAu8KpgX4TPBZrwWOArdV8LdIBezZs4f169ezd+/euEMRkZAqlkTc/Z+B3hlOvg74G3d/xd1/CuwD3hA89rn7c+4+CPwNsM7MDHgzsDOY/17g7bMZv1RWPp9n69atuDtbt24ln8/HHZKIhBDHOZEPm9ne4HDX4qBsOfBC0TT7g7Lpys8Fjrl7flJ5SWa20cy6zKxLh0+Sob29naNHjwLQ29tLe3t7zBGJSBhRJ5EvAq8BrgAOAv87ii91923u3uLuLc3NzVF8pZxGT08P27dvZ2BgAICBgQG2b99Ob+9MG64ikhSRJhF3P+Tuw+4+AtxN4XAVwAHgwqJJVwRl05X3AIvMLDOpXFKgs7OTkZGRCWUjIyPkcrmYIhKRsCJNIma2rOjlO4DRnlsdwM1mNs/MLgJWA/8CPA6sDnpi1VM4+d7h7g58C1gfzH8L8FAUv0HOXmtrKzU1Exe9mpoaWltbY4pIRMKqZBffrwE/AH7RzPab2W3An5rZE2a2F7gO+H0Ad38KeAB4GngY+FDQYskDHwYeAZ4BHgimBfgj4KNmto/COZIvV+q3yOzKZrNs2LCBhoYGABoaGtiwYQPZbDbmyESkXFbYqa8eLS0t3tXVFXcYVS+fz/Pud7+bF198kQsuuICvfOUrZDKZM88oIrEws13u3jK5XFesSywymQybN2/GzNiyZYsSiEhKac2V2KxZs4adO3eiHnMi6aWWiMRKCUQk3ZREREQkNCUREZEES/ooG0oiIiIJlYZBSpVEREQSKC2DlCqJiIgkUFoGKVUSERFJmDQNUqokIiKSMGkapFRJREQkYdI0SOmMkoiZtZlZt5m9bGbHzeyEmR2vdHAiItUoTYOUzrQl8qfATe7+Kndf6O7nuPvCSgYmIlLN2traaGpqAqCpqYm2traYIyptpknkkLs/U9FIRESkJDOLO4RpzTSJdJnZ/Wb2ruDQVpuZJTMtiojMAe3t7fT19QFw8uTJ1HfxXQicAq4Hfj14/FqlghKR6pb0oT4qbc518XX395Z4vK/SwYlI9UnDUB+VNue6+JrZL5hZp5k9Gby+3Mz+uLKhiUi1SctQH5U257r4AncDm4EhAHffC9xcqaBEpDqlZaiPSpuLXXwb3f1fJpVV5y7CLKj2470ipaTpPEAU2traxpJGNptNfRffI2b2GsABzGw9cLBiUc1hOt4rUlqazgNEIZPJsHnzZsyMLVu2kMkk827mM00iHwLuAi4xswPA7wEfqFRQc5WO94pML03nAaKyZs0adu7cyeWXXx53KNOaae+s59y9FWgGLnH3a9z9+YpGNgfpeK/I9NJ0HiBKzc3NcYdwWjPtnbXUzL4M7HT3E2Z2qZndVuHY5hQd7xU5s7ScB5BxMz2cdQ/wCHBB8PonFA5pyQzpeK/ImaXlPICMm2kSWeLuDwAjAO6eB4YrFtUcpOO9IjOThvMAMm6mSaTPzM5lvHfW1cDLFYtqDtLxXpGZS/p5ABk30yTyUaADeI2ZfQ+4D/hIxaKao3S8V0TmmjMmETOrBd4UPH4ZeD9wWXDVupRBx3tFZK4xdz/zRGb/4u5viCCeimtpafGurq5YYzh8+LCa6yKSKma2y91bJpfPdFf4e2b2F8D9QN9oobv/aJbiqypKICIyV8w0iVwR/P1E8NconGR/82wHJCKi1np6zPTE+j8Afw/8Y/D4e+CvzeyKCsUlIlVK48uly0yTyFXA7wLLKFxw+H7grcDdZvY/KhSbiFQZjS+XPjNNIiuA17v7H7r7H1BIKucBvwrcWqHYRKTKaHy59JlpEjkPeKXo9RCw1N37J5WPMbMdZvbS6N0Qg7KsmT1qZt3B38VBuZnZnWa2z8z2mtnri+a5JZi+28xuKSq/ysyeCOa508ysjN8tIgmj8eXSaaZJ5KvAD83sDjO7A/gehXMiTcDT08xzD7B2UtnHgE53Xw10Bq8BbgBWB4+NwBehkHSAO4A3Am8A7hhNPME0txfNN/m7RCRFNL5caUm/id1Mh4L/FIWN+7Hg8bvu/kl373P3355mnn8GJu9CrAPuDZ7fC7y9qPw+L3gMWGRmyyicd3nU3Xvd/SjwKLA2eG+huz/mhQtd7iv6LBFJIY0vN1UaOhnMtCWCu3e5+58Hj7BX6y1199E7Iv4cWBo8Xw68UDTd/qDsdOX7S5SXZGYbzazLzLqSntWrjf4fMmp0fLn6+noA6uvrq3p8ubR0MphxEpltQQvizJfLz853bXP3FndvUd/z5EjDXpZEa926dQwPFwYIHx4eZt26dTFHFJ+0dDKIOokcCg5FEfx9KSg/AFxYNN2KoOx05StKlEtKpGUvS6L10EMPUVtbC0BtbS0PPfRQzBHFI02dDKJOIh3AaA+rW4CHisrfE/TSuhp4OTjs9QhwvZktDk6oXw88Erx33MyuDnplvafosyQF0rKXJdEZ3XAODg4CMDg4mNgNZ6WlqZNBxZKImX0N+AHwi2a2P7id7qeB/2Jm3UBr8BrgG8BzwD7gbuCDAO7eC3wKeDx4fDIoI5hmezDPvwHfrNRvkdmVpr0siU6aNpyVlqZOBhVLIu7+Lndf5u517r7C3b/s7j3u/hZ3X+3uraMJIeiV9SF3f427/1LxiXt33+Hurw0ef1lU3uXurwvm+bDPZDhiSQRtLKSUNG04Ky1NN7GL7cS6VC9tLKSUNG04o5CWm9gpiUjktLGQ6aRlwxmFTCbDb/zGbwDwm7/5m4m9iZ2SiMRCGwspJS0bzijk83keeOABAO6///7E9mBUEpFYaGMhpaRlwxmF9vZ2enp6ADhy5EhiezAqiUgstLGQUtT1uyBN3Z2VRCQW2ljIZOr6Pa6zs5OhoaEJZUNDQ4nswagkIpHTxkJKUdfvcVddddXY8C+jhoeHaWlpiSmi6SmJxKDaBx3s7OwsuYJU48ZCxqnr97hdu3aVrIuurrBj31aOkkjENOhgYWMxeY/T3atyYyHjRrt+j95frpq7fre2tlJXVzehrK6uLpHriJJIhDToYIG7oxtRSiltbW1jG89q7vqdzWa5/fbbJwyLf/vttycyoSqJREgnkws6OzunJJGamhodzhIymQwrV64EYMuWLVXd9butrY0lS5YAsGTJksQmVCWRiOhk8rjW1tax4b5HVeuxb5lqwYIFXHbZZVx++eVxhxKrTCbD5s2bMbNEJ1QlkYio58m4bDbL+973vrHWiJlx2223JbKpLvGYfD6gWq1Zs4adO3cmOqEqiUREPU9OT4Mwi5SW9LuxKolERPePHtfT08OOHTvGEoe7s2PHjqo8tFes2rt+S2lJXy6URCK0bt26sR5Z+Xy+au8fretEplLXbyklDcuFkkiEHnzwwbHzIiMjIzz44IMxRxQPXScykbp+SylpWS6URCLS09PDXXfdNaFs27ZtVXkIR9eJTKSu3xN1d3fT3d0ddxixS8tyoSQSkY6OjimHcPL5PB0dHTFFFB9dJzJOXb+llDQtF0oiEdGe97jW1tYpvbFGRkaq8nCWun5LKWlaLpREInLTTTdNSSRmxk033RRTRPFxd3XpDajrt5SSpuVCSSQi7j7lKu3Jr6tFZ2fnlBXEzBK5l1Vput+8lJKmSwKURCLS2dk5Ze/b3atyw1lq2JPa2tpE7mVFQfebl1LSckmAkkhE0nSTmUrTkN8TjY6RBBp0UMZNviTg7/7u7+INaBpKIhHZtWtXyR5JSbzJTBQ05LfI9Hp6eti2bduEsrvuuku9s6qZeiRNpCG/x41eVAYk+qIyiU5HR8eU5SCplwQoiURkuj2IY8eORRtIgmjI74K0XFQm0UnTJQFKIhHZsWNHyfLt27dHHEmyVPuQ32m6qEyic9NNN03pfJLJZBJ5SYCSSEQ2bNhQVrlUh87OzimHLap9MEopnCd8//vfP/bazNi4cWMiO58oiUTk4osv5pprrplQds0113DxxRfHFJEkQWtr65QkMjQ0VLXnygD6+/vp7++PO4zYrV+/fuxcYTabZf369TFHVJqSSITuuOOOsee1tbUTXkt10rkymU4mk2HVqlUAfOITn0hs5xMlkQjNmzeP5cuXA7Bp0ybmzZsXc0QSty996Usly7/whS9EHIkkURo6nyQztc1hzc3NNDU18Y53vCPuUCQBpjvGncRj3xKPpHc+URKJQWNjY9whJILuGTH94azRLr8iSafDWSIx+q3f+q2yykWSJpYkYmbPm9kTZrbbzLqCsqyZPWpm3cHfxUG5mdmdZrbPzPaa2euLPueWYPpuM7sljt8icjb27dtXslytNEmLOFsi17n7Fe4+OgLhx4BOd18NdAavAW4AVgePjcAXoZB0gDuANwJvAO4YTTwiadHa2jo23Peo+vr6qu7iK+PScKvgJB3OWgfcGzy/F3h7Ufl9XvAYsMjMlgFvBR519153Pwo8CqyNOOaypWGhkOhks1k2btw49jqTyST2ojKRUuJKIg78k5ntMrPRNWipux8Mnv8cWBo8Xw68UDTv/qBsuvIpzGyjmXWZWdfhw4dn6zeIzIq2trax1sh5552nEY0lVeJKIte4++spHKr6kJn9avGbXhjudtbun+ru29y9xd1bmpubZ+tjRWaFRjSWNIslibj7geDvS8CDFM5pHAoOUxH8fSmY/ABwYdHsK4Ky6cpFUicNF5VFZWRkZOxmTJJ8kScRM2sys3NGnwPXA08CHcBoD6tbgIeC5x3Ae4JeWlcDLweHvR4BrjezxcEJ9euDMpFUSvpFZSKlxNFuXgo8GIyXnwH+2t0fNrPHgQfM7DbgZ8A7g+m/AbwN2AecAt4L4O69ZvYp4PFguk+6u8bPFhGJUORJxN2fA9aUKO8B3lKi3IEPTfNZO4DSN+pIKI1OOk51MU499iStdAZPRCSh0rCjlaTrREREJGWUREREJDQlERERCU3nRCKm/u9SShqOfYuUopaIiIiEppaIxEatMpH0UxIREUmoNOxo6XCWiIiEppaISAKkYY9TpBS1RCR2HR0dcYcgIiEpiUjscrlc3CEkRrUn1OLfX+11kRZKIjHRCiKlVHtCzeVyzM9kmJ/JVH1dpIWSSEyqfQUpTqJHjhyJMRJJmgsXLuTChQvjDkNmSElEYpHL5Wial6FpXoajR4/GHU6slFALOjo66O7u5oXjx3nh+HG6u7vVYk8BJZEIaWMx0UVLF3LRUu1x5nI56hrqqGuoq+qEmsvlGBkYYEVjEysamxgZGKj6FnsaKIlEKJfLkalvIFPfUNUbi9E9zucOHee5Q8fp7++v+j3OV53/Kl51/qviDiN2K5oW8NHLL+ejl1/OiqYFcYcjM6AkErEFzctY0Lws7jBilcvl8KEBLloyn4uWzKchY1W7xzmaUI8dPMaxg8eqNqGO1sP+vpP82d69/NnevezvO1nVh7TS0lNNSSRCR44c4cThg5w4fLBqNxajLmpu5FPrL+VT6y/loubGuMOJTS6XYyA/QNPSJpqWNmF11ZlQc7kcI/39LG9owIfz+HCe5Q0NjPT3V2V9QKFO5tXPZ179/ETXga5Yj9DoIazG7Pmc6v05uVyOm266KeaoojW6x+lDA3x859MA/PTwKexYYY+z2uoDYMHSBVz57isB+PFXfhxzNPFZPn8+my6+eELZnc89F1M08RpdT8ABxlpkSVw/lEQi1pg9n0ve9l6e/cZfxh1KLAqHsvpZlZ2HD+cBWJWt5/ne/qpLqqMbioH8wFjyOHnoJN09yd1gzJaOjo4Je9fd3d0M9/dPSRr7+/up7e5m06ZNY2Wtra1zum6gsJ7kh4Y5L7scgJd6DyR2/VASicDoCtPf348NHuTZb/wlp3oO0n3iMJs2bZrTK0WpjUWwczWRF96rpo1FLpejP9/P/Ob55EcKCbWhuYH+w3M/oeZyOZ7ds4fl9fUAnA9QV4fn8xOmW15XB/k8J555BoADg4MAc7puRp2XXc5v3/BBAL76zS/EHM30lEQikMvl2PvUM9SfuwKAoeER6hYtZRjY+1Rh5ZirK0Uul+MnT+1h1eJaAF7dBDTVwMjQhOlWLa4BXmHwxScBeP7oMDC36kUJdaLl9fV8cOn5Zc3zhUM/r1A0yTHaQs0PDo8lj5d6DtB7ojaRLVQlkQootbEota2AwjakeIMxFzcWqxbXcsdbyuu++onOlysUTXxyuRx7nt5D/bmFvW8WQh115Icn7n3XnVtHnjzPHCrsYAz2pHfve/K6MKq7u5vhwcGyk8KBwcEph7dGzZV1J5fLMTQ4TPOiZYwMF0Z3XrJoGYePHUxkC9Xcp9u8zU0tLS3e1dVV0e/YtGkTe558htqFzWXNN3z8MGte9x+48847KxRZ5ZxuY+GDfaxaXN7+yvNH81h9E6tXry75fpI3GKeri1NDp8aTyAwN9gzSWNdYsi6SXA9QWBee3b2bZbW1Ff2eg8PDXHLFFXNi3enu7h5LIsUOHztIXX3thOUgyv+/me1y95bJ5WqJnIXTbjinbXtMz/EphzGKJXmDkcvl+MmTu3n1pAbHygaggbGT6DP16oUAfbzywu4p7/0saKQkuS52P7Wb2uykDecCqKFm7PzHTNUsrmGAAZ44+MSE8uHedBzyW1Zby/sXnFPR77jr5ImKfn4l5XI5nn7yWc5dWDi0t3BeM8yDkeGJ25Bzzym8f+jfjwHQc7zQiov7/6+WyFnYtGkTu598GluwZFY+73T85BGueN2lFd/Tmi4xnkl3dzf+St+UJFIJP3sZbN70rZTplJuEz6Yu+gb7piaRWTbcO0zTaVprpxNlXQz39UXSEqltiqYuwjrdTufQK3nOXVjeRcg9xw9SNy8TWQtVLZHTOJsVJESDI5wSJ1tnotyFqdCi+BErFwyX9T0X1gF14OXNFsrKBQDHGHj+8RnP8+8nCxuxcuvix0/9GBaVFx+NhUe5LY6yLYKTnOTHB8q8tuRY4U+5dfHs7t2UdxocRg/o+nBlF4zzAfr6OLZ7d1nzjZ6RiSKJ5HI5nnriGbILlk4oX5A5FzKMnf+YqcVNhc85+NPeCeW9Jw8B0bVQlEQINhZPPM1IY7a8GWvOKWwwyvznh9KY5Tiw699mfiKy5lRh4Sp3YWoYOcU1SwbKmifp2o83hJtxEYxcO7duXVvz7XADVZwP3IbNbjAx+3KIvcCz2umMSJgdTgjXglESCYw0Zhm49NfiDmNWNTz9D3GHkHoDLw/g3XPrkK+9bLC8/Pl6BgaYedsvHXoov6E52qJY1HheWfM11mRpbISRkcouT6NxHfi3nrLmO3bqJaD8nU4lkcBgz37sh38VdxizanCwD15T7gEIGKhp5P8dmVeBiOIzUBPymHwf2NNza++bvnCzHQO+NZtxJMCxkPPZcIblCy+ZzVBi9/KJ3jNPVIKSCIUm3MGDB8ue78SJE/T1TV0jBwcHCdthwcyory/dBbSpqYlzzimnl8siWltby/r+cqcfdeTIkZLD2w8NDZEvugrZ3U9bNzU144daMpkMdXV1JadbvHgxS5bMvEPDL1D+bwtbF93d3fSd6qv82pWHpsbwJ5MrOf2o7u5uTvX1UeldkleAxrM4sV6uvsGXefrA98ueL8n6BsNdm6XeWRUw+ZjpkSNHePHFFxkZmXpsvba2lmXLlo1tEJPcjXe2XHfddQyXONFaX1+f6NFKZ2q6Y+aTE+3Q0BCDwTAek5kZjY3joxtPlzSTvrzMRl1kMhnmzRtPQ3HXRdhzIlHuaJW7kzXqdHU4Xe8sJZGIXHvttSWTSE1NDd/+9rcjjydObW1tJe/suGTJEtrb22OIKD5aLsatXbuWU6dOTSlvbGzk4YcfjiGi+LzpTW8qmUjMjO985zsxRDR9EtH9RCJy6623llU+l3384x8vq3wu03Ix7iMf+UhZ5XPZzTffXFZ5nFKfRMxsrZn9q5ntM7OPxR3PdG699VbMJp6gNbOq3FhceeWVXHjhhRPKVq5cyZVXXhlTRPHRcjHuxhtvZP78+RPK5s+fz4033hhTRPH5wAc+UFZ5nFKdRMysFvg8cANwKfAuM7s03qim9/Wvf/20r6vJ3XffPeH1tm3bYookflouxt1///2nfV1N0lIXqU4iwBuAfe7+nLsPAn8DrIs5pmmdd955XHDBBQCsWLGC884rr5/5XNLY2Mj69esBeOc73znhJHK10XIxbtGiRaxZswYotFgXLVoUb0AxWrZsGeeeey4Azc3NLFtW3rAokRntCZDGB7Ae2F70+neAvygx3UagC+hauXKlx+1zn/tc3CEkxve///24Q0gMLRfjHnjggbhDSIzPfOYzcYfg7u5Al5fYDqe6d5aZrQfWuvuG4PXvAG909w9PN09cvbNERNJsrvbOOgAUn6FdEZSJiEgE0p5EHgdWm9lFZlYP3Ax0xByTiEjVSPWwJ+6eN7MPA48AtcAOd38q5rBERKpGqs+JhGFmh4GfxRzGEmDqJdvVSXUxTnUxTnUxLil18Wp3n3LP76pLIklgZl2lTlBVI9XFONXFONXFuKTXRdrPiYiISIyUREREJDQlkXhU7xgfU6kuxqkuxqkuxiW6LnROREREQlNLREREQlMSERGR0JREKsTM3m5mbmaXFJXdYmbdweOWovL/ZWYvmNnJeKKtnJnWg5k1mtk/mtmzZvaUmX06vqgro8xl4mEz2xPUxZeC2x7MCWY2bGa7g9/3IzP75aL3qmodKbcuErmelBqVUY9ZGWH4fuC7wCeC11ngueDv4uD54uC9q4FlwMm4446rHoBG4Lpgmvpgnhvijj/GZWJh8NeAvwVujjv+WayHk0XP3wp8Zwb1MSfXkXLrIonriVoiFWBmC4BrgNsojOcFhQXkUXfvdfejwKPAWgB3f8zdD8YSbAWVUw/ufsrdvwXghXvD/IjCgJpzQohl4ngwTYbCxmKu9oBZCBwNnlfdOjLJGesiietJqsfOSrB1wMPu/hMz6zGzq4DlwAtF0+wPyuayUPVgZouAXwf+PKpAI1B2XZjZIxRuvPZNYGeUwVbYfDPbDTRQaF28OSivxnUkdF0kZT1RS6Qy3kXhLosEf98VYyxxKrsezCwDfA24092fq2BsUSu7Ltz9rRQ2LPMY37jMBf3ufoW7X0KhpXGfTb7RfPUIVRdJWk/UEpllZpalsML/kpk5hdGFHfgj4NqiSVcA3446vqicRT1sA7rd/f9EEmgEzmaZcPcBM3uIQkvm0SjijZK7/8DMlgDNFO4FdG3R23N6HZmszLpIzHqilsjsWw/8lbu/2t1XufuFwE8pLBTXm9liM1sMXE9hCPu5qux6MLP/CbwK+L2YYq6UsurCzBaY2TIY2+O8EXg2ruArKeipVgv0UFgOqmkdmWCmdZG49STu3glz7QF8i8IJsOKyTcAXgfcB+4LHe4ve/1MKxzxHgr9/EvfviLoeKOxpOfAMsDt4bIj7d8RUF0sp3HBtL/Ak8DkgE/fvmMX6GC76H+8Bbix6r2rWkTB1kcT1RMOeiIhIaDqcJSIioSmJiIhIaEoiIiISmpKIiIiEpiQiIiKhKYmIiEhoSiIiFWJm3zazluD588HVyCJzipKIiIiEpiQicgZm9t/NbFPw/LNm9n+D5282s6+a2fVm9oPgpkJfD4Z9L+fzV5nZM2Z2d3CjoX8ys/nBe7eb2ePBTYv+1swag/J7zOyLZvaYmT1nZtea2Y7gc+4p+uyzik3kTJRERM7su8B/Dp63AAvMrC4o2wv8MdDq7q8HuoCPhviO1cDn3f0y4BjwX4Pydnf/j+6+hsJQF7cVzbMY+E/A7wMdwGeByygM9HhFcPhsNmITmZZG8RU5s13AVWa2EHiFwo2AWigkkQ7gUuB7wQje9cAPQnzHT919d9H3rQqevy4YcG8RsICJAxL+vbu7mT0BHHL3JwDM7Klg/hWzFJvItJRERM7A3YfM7KfArcD3KbQ+rgNeS2E03kfd/WzvGfNK0fNhYH7w/B7g7e6+x8xuZeLw4KPzjEyaf4TCuj08S7GJTEuHs0Rm5rvAHwL/HDz/XeDHwGPAr5jZawHMrMnMfmEWv/cc4GBw+Oy3y5y30rGJKImIzNB3Kdxl8AfufggYAL7r7ocptFC+ZmZ7KRwuumQWv/fjwA+B71HmPUUiiE1EQ8GLiEh4aomIiEhoOrEuEhEzOxfoLPHWW9y9J+p4RGaDDmeJiEhoOpwlIiKhKYmIiEhoSiIiIhKakoiIiIT2/wGlt4TRVcoYqwAAAABJRU5ErkJggg==\n",
"text/plain": [
"