A1_20230523_localmax.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "c3e47375-5973-4818-86a2-b65d238763ea",
   "metadata": {},
   "source": [
    "# env: pm-minmax"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "f81badeb",
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt \n",
    "import numpy as np\n",
    "from PIL import Image\n",
    "from scipy.ndimage import maximum_filter"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "c2fea336",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x15bf98d0b80>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAEpCAYAAAAzlkD8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABG6UlEQVR4nO3de4yc1X0//vfcbzs7a3u9F8fGcYIhwQZKcOtASbgkuHFTFCCtSNKmRlGjEC4FuVESg77B9KewkCiIVBC3UOQaBep8JXBKBcU4IjaNqBPj4C8OBGqKbdbY67V3vde5z5zfH5RN1nveH5hl/WB23y9pJTxnnpnnOec5zx5mn/d8Qs45BxEREZGAhN/rHRAREZGZRYsPERERCZQWHyIiIhIoLT5EREQkUFp8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERAKlxYeIiIgESosPERERCVT0RL3wj370I3z/+9/HoUOHsGTJEtx99934xCc+8bbb1et1HDx4ENlsFqFQ6ETtnoiIiEwh5xyGh4cxb948hMNv89mGOwE2btzoYrGYu//++91LL73kbrzxRpfJZNz+/fvfdtvu7m4HQD/60Y9+9KMf/bwPf7q7u9/2d33IuakvLLd8+XJ87GMfw7p168Ye++hHP4rLL78cXV1d5raDg4NoaWnB0qv+DyLx5IT20Xn805DiKWXv4x86pZduc2nby7Rteep/aNuHYkXaFoF/H4/WeVfvq8yibU8OnEXbfnn4FNo2+HrO+3j6QIRuk+zn+xgfqtM2ph41xms2XxkX2vh+VBaWaNup8/xj/fE5e+k256VfpW0LosO0LWZ8MDdU9/fxr0sL6DbbBz9M2355cCFtK+/L0ram1/07mezjYxkfqtG2kHG1qCV5h5Sa/f1RaOXbjC6q0ra2Bcdo2yc7+Hj+YeY17+PzowN0m1mRCm07VovRtpcrHbTtV0P+sd55dD7d5ujrLbQtdZB/iJ066h+0+CA/B8L8FEA1xces0MrndLHVvx/ulALdZun8g7Ttotmv0LY/TO6jbWclEt7HS46P80jd/7sFAHaX+fx7engJbfv5G4u9jw/tbaHbpA/w/s0c5uMZG2n82l1L8HGuZPz7UasUsfv//n8YGBhALuf//fOWKf+zS7lcxs6dO/Htb3973OMrVqzAs88+O+H5pVIJpdLvfpkMD795sY/Ek97FR8TokHDK3yHRjP9kA4BkE++CTJoPdDbG29jio2gsPtIVviCIV/kFLjLMjy2cmth/ABBJ8PeKxPk+RmNTu/iIxI0+TPL9qBkXPzbWySbeh5kM749slO+jtfio1/3bpWL8fIvX4rQtkjbGOekfZwCIxP07aY1lNDa5xUfI6JBq3N/H9nzmi4+IMacTxlinm/z70WSMczbC2yo13pYqG2Nd9+9jJN/4fAaASIK/F5vT1jlgfWLuyDn15n7wDcNkTrs0P6liGT4nUsa1u4n8LgCAZrKPJce3CZH5DACZMr9+JJxx7SZz2pzPRv9GzDnd+LXbms9149oN4B3dMjHlN5wePXoUtVoN7e3t4x5vb29HT0/PhOd3dXUhl8uN/SxYwP/PUERERN7/Tlja5fiVj3POuxpas2YNBgcHx366u7tP1C6JiIjISWDK/+zS2tqKSCQy4VOO3t7eCZ+GAEAikUCC/A1OREREpp8pX3zE43Gce+652LJlC6644oqxx7ds2YLPfe5z7/r1a9Y9CeRvxJ3pQbrNvBi/ca09wm+ESoZ41xWd/2/mRcf/Nvjr/Adp24sD/Ma1Ywf4TT1N+/3vlz3A//6X6uU3XUUK/G/wjvwNsJLh/VTOGn9TNM7MSJTvf1PMfzPq/Hg/3aYlzMc5YvzpsmLc/9BdbfE+/t8FPpa7jnyAthVfN24qPch3MvuGf8ziA3wso4P8hl5n3C9UzfC/b1fI/VPWTW1I8HGenx2gbZ1xPt87ov62duOm0gHj7/1Djv+P04t5Pp6/OdbpffzIgRa6TdM+4/60g7yvkv3+61Fs2JjPxokfMiZnqYX3Fbv8JZK87+ckRmmbde2eH+XHBvjHrOj4NoeNe3teKvKbhP+rdxFtY9fu5v3G/Sr7+f1YyX5+U2xkiLfVU/7xLM/i99uUsv59dA18O8YJ+Z6P1atX48tf/jKWLVuG8847D/fddx9ef/11XHPNNSfi7UREROR95IQsPq666ir09fXh7//+73Ho0CEsXboUTzzxBBYu5HFBERERmRlO2DecXnvttbj22mtP1MuLiIjI+5Rqu4iIiEigtPgQERGRQGnxISIiIoE6Yfd8vFsu4o9m1Yyv4p2d88eyPpI5TLdZEOujbXMjvHtiIR45PFLzR7b2VWbTbX4zPI+27T80h7ZlXuf7mO32x++aXufR0uhAnrahwqNo9Wza+3g1yfvJitM6o60pw+vqnNbkr+3SEuHHNduIWfLwKHCwlqJtL5f8UcoXBnn88uhBHpvOHOb/n9B0gMfvMvtGvI+H8jxOG6ry16tn+TGHSGQPACoZfwbPms+ZFn6ezkvxOO2HEny+zw37j9soZWJ6ucTn7b48n7dvHG3xPp7ex8+45n08Tps+zMcz2sf7kXFJPpbVNJ/TxjeUo06+JqElzffv1DSvy/UBox5POmTNXL98nZ8Fu42aTDuGeJDijSMttC2z39/HmUPGOL/Br2ORYzyWbAr753R01Lh2x/373siXuOuTDxEREQmUFh8iIiISKC0+REREJFBafIiIiEigtPgQERGRQJ20aZdyUwgRT9GpWgtPJsxrGvI+3hobptssiBrpDqN7huu8UM/Bmj/58avRD9Ntdvf60xEAEHudF6/K7rfugPfvY+z1I3QbV+b9G4ryu5/DMX9fhWpJuo2lmuV3nlt3x8+K+e/4nhflRahqRoG4UeP2/RdLPLmye9hfbOqVgxMrO78l9bqRdNjLxzn7mj/RAgDhoyQV4vhBm+dAgu/jZJJNlRw/rvYMn5uLUzwF8aEoLyKYDftTN0dqvCLWa1UjpTbKz4EXDvMkTGifP2GQOcTHpWkfTzNEBvmcCOVJOqzO+941+a9hAIA2nngyameinvHP6XlNPLnUTgoBAsDcCL8Gx0J8H1+v+udLd5Uf8wt5nnb59SHeFjnAr39srM00Yi/vDzdqnAMRfh0Lh/znfrWN90c95t+m3kBlOX3yISIiIoHS4kNEREQCpcWHiIiIBEqLDxEREQmUFh8iIiISKC0+REREJFAnb9S2xSGcnBhFimV4vGpx1h+/+3CcF5qyVByPou2v8sjhL/On+h/v+yDdZuRQE22b9QZtQqaH90fif/z94Uq8CJUr8KJtSBvxOxLXqmT5KVaaxWNZoRSP2nam/ZFqAMhF/HGzuhGZPWwUiHujOou2vZTnUcpdR0hbN3+vzEEjZnmAj0ukj0fJaz2Nn/vhFl7gzkV4lrKa5n1cavGPdWQ2PxfnNw3QtvYYb7PkScS4v87HpbvMC8Tt6vNHqgFgtDdD21oO+Psj220UiDvE4+JumJ8D9aL/NUPZLH89Ep0HeMwSAKpGsj7S5I9wz4rzSPUHYpOLyI/UeT/21/zX7l8XFtFt/usob8sf5NfuHBlnAMh2++d09Ai/vrkBI2pbM8q6JfnAuFTc+3ilyZrP/sdrvNsn0CcfIiIiEigtPkRERCRQWnyIiIhIoLT4EBERkUBp8SEiIiKB0uJDREREAjXlUdu1a9fitttuG/dYe3s7enp6GnqdWsrBeaK2HS28gmdrzN8WAc9kVYy41sE6j0ntq7TStlfy/sqlew/zyF76AB+K7IEqbYsf5LGsOolluQKvfujqvEOizTyaV23yV96tZHgflmfx90o382hpNmbEgYmBOo8JD5AqxADwWqmNtr04yCsR97/R4n0828v7o3m/EbPsnWT8ruo/d8IZHgOFEdlzCR61rRhRWzbWTRk+lgvSPGbZEubxzJJRWrWfZEGtmOWuYV619ODRFtqW3s/ndHO3f1zih3nlWjdonANknAHAkXhx2Kh0WkvzrxIo5Xj/WnO6ucl/3ZlDKlEDQMUZlcWNtoMVf3wUAHYVF3of/+UgPwf2v8Gv98ke3h/N+xq/doeM6rQ1ow2kYjMAhCdx7S618POjONd/jagXjLjvcU7I93wsWbIEP/vZz8b+HTG+G0BERERmlhOy+IhGo+jo6DgRLy0iIiLvcyfkno89e/Zg3rx5WLRoEb7whS/gtddeo88tlUoYGhoa9yMiIiLT15QvPpYvX44HH3wQmzdvxv3334+enh6cf/756Ovr8z6/q6sLuVxu7GfBAv73VREREXn/m/LFx8qVK/H5z38eZ555Jj796U/j8ccfBwBs2LDB+/w1a9ZgcHBw7Ke7u3uqd0lEREROIie8sFwmk8GZZ56JPXv2eNsTiQQSCf/dtiIiIjL9nPDFR6lUwm9/+1t84hOfaGi7Wq4Kl5oYU5qfHaDbtMf8kcNwiMd/uqvNtK23xuNJvxg6jbb96pA/yhU6wCtnpg/xiFqij0cwQyVe1bZe9reF4jyGFiLVaQEASb5IrDT7X7PUbMT5EnxccmkeKYuEeF+N1PxRyr4wrzy5v8RjdL8e5H8G/J+eubQt3udPeKV7+DFHh/lYhvI8klov+6uFAgDC/v2wzgG08Uh4qZWfw8U5RrXTWf59nJ3hkdlkmB9X3fjQ9nWjEnFPtcX7+J4Cj1TTCsUAcJBXC00dNWLro/6qzeFB/lUCzohSssq1ABBpJud+mo9leTY/LlahGABqGV6Nujnp38ewMZ/7anzeFh2PA+8r8zn9/0h0eucbfK5He/h8aTpgXLv7jWs3ic3WjUg1jErroTgfz/oc/nuu3Ozvx3KWj7PL+uemi/Jo8fGm/M8u3/jGN7Bt2zbs3bsXv/zlL/Hnf/7nGBoawqpVq6b6rUREROR9aMo/+Thw4AC++MUv4ujRo5g7dy4+/vGPY/v27Vi40P9pgIiIiMwsU7742Lhx41S/pIiIiEwjqu0iIiIigdLiQ0RERAKlxYeIiIgE6oRHbScr0VxCJD0x6hMP8yjPcM0fNbJiV5b/N3oKbdvVN5+2DRzxx8MyfTy6lOrnxxUuGPElI2IXIgX9QlE+7CEjflebw6PH5Rb/a5Z46hEuzWN52Tg/rlKN7//Rir/v3yi10G3252fTtpeP+CsUA0D9MI8jZsl35SX7+TFHhnmc1uWNSsQV49yJkzgiexxALcfPgXw73644h0cOIxl/NC8T4/Fiy/+UeTS27owYbsk/1i/0f4Bu038oR9uajCrFiUEei4yRWLUz5jOMAp00TgsglPPHLKutfD5blWuLPIkNZPi5mIj42yLGVyEcqfJ9PFDn8/bVPD8/dvb4r92lHl7duumINc7GnB4xxrNGtjO+7iDcxKtRh1r4eVrJ8K9JGO3wX0+LrXw+J7P+46pFjOM9jj75EBERkUBp8SEiIiKB0uJDREREAqXFh4iIiARKiw8REREJ1EmbdonHaojEJt4dbRUVO1T23+2bj/KiQL1lfjf1bwZ4QamDffzO4niPPxGQNApNRUq8LeR4m1XsLczujg8bxd5a+XEV23kKYrTd/16VLN/3cNJILpX5cQ3GeMpkuOrfrjfPx/nIML+DPN/L21J9vB9jo/7jjg8a6Y4qv2s+ZI2zUXAsFPOfi/V2nhQozeH9W242isc18dRCPOY/tmSEF48brPLzrVLnaYwDxRba1j3ij18dOMJjWbGj/DKZOGYUFRvg53eo4u+PUJQfF2C0GQm2eos/CTPZca4l+DHHUnw8qySFdKDA+/5gqIW2sbkOAHv6ecJxZMCfakn28v5N9vFjjg++82Jq45AxCzfza5WVeLLSiIVOY6xJocBqll+Pkp7fzQAQJvPc+9x3/EwRERGRKaDFh4iIiARKiw8REREJlBYfIiIiEigtPkRERCRQWnyIiIhIoE7aqG0iVkUkNjFWNFjmkaEoKVD0eoHHCo8UeEGm1/t5BKx+jMe8EqP+6FKkYsROKzymiDpvc0aBMKT8++g8/fqW8lwebyy28O2qJJFaSxkRYp7mw0iR9+/RMB+z4ZI/Vj04xCOztREjStnL29KH+bGxQoHmOFuRaiMeHUrzglgsOl2ZxedRoZUfc4knseHifP+jUX8E71iJ73s8zGN7+So/7w8bserD/f4ia+4IP99iw/xEjRb5eIaqxniSOe2y/Dw153qE72Nprr+Pi3Os4nH89epJfszW/80OFPzn3CtVXgTOctSY0+U8/3qFaI+/LdHH3ys+Yhyzee3mTY5cnxHh1+C6USCu0szbCrP4yJT9UwIuw+dfNOw/sBB53EeffIiIiEigtPgQERGRQGnxISIiIoHS4kNEREQCpcWHiIiIBEqLDxEREQlUw1HbZ555Bt///vexc+dOHDp0CJs2bcLll18+1u6cw2233Yb77rsPx44dw/Lly3HvvfdiyZIlDb1PqRJFpDJx9w6NklwQgKMkNluo8Ija4AiPHFb6eVt8wKhoOuJ/PFrk0bt61Ii2pXlszBnRKxf172OlmQ97uYnH70ot/JgrTf5jCxsFH+tV/nqjo7zvS2Vj/0nEzuX5ccX7jbZB2kQr1wKgEbt6nL9XqMmI2GV5JNXFjHFp8Z8flQzfjwpPMqOWNI7ZaCqTMesd5m/WN2pEiOv8mPN5PidYRD4xifkMAGFexBW1hFE9uonE4Gfxc8CKbVYzfE4USKS2aMQva3z6wcX4QFfL/Lwadv4XHTEi95USPy5nROQjo8acPtb4VyFERycXkXcp/run3kSu63X+elYlYis6Xc4ZVYqNr0NgqmT+WfPyeA1/8jE6Ooqzzz4b99xzj7f9e9/7Hu666y7cc8892LFjBzo6OnDppZdieHi40bcSERGRaajhTz5WrlyJlStXetucc7j77rtxyy234MorrwQAbNiwAe3t7Xj44Yfxta997d3trYiIiLzvTek9H3v37kVPTw9WrFgx9lgikcCFF16IZ5991rtNqVTC0NDQuB8RERGZvqZ08dHT0wMAaG9vH/d4e3v7WNvxurq6kMvlxn4WLFgwlbskIiIiJ5kTknYJHVe0wzk34bG3rFmzBoODg2M/3d3dJ2KXRERE5CQxpYXlOjo6ALz5CUhnZ+fY4729vRM+DXlLIpFAIsHvThcREZHpZUoXH4sWLUJHRwe2bNmCc845BwBQLpexbds23HnnnQ29Vj4fRxgTFyWFIo8uORLzqY7wbUJl/uGPFaeND1iVLv3RpWqCbxNO85gUQkac1oip1ZIkaps2KmAacdo6T/waO2JUtR3g41JL8GhbvcJ3JFL073+UVBoGgGieNiE2wvc/UjIiauTtakljnMGPy9qubkRtSzkWs+T9UWkyTiorT1vj25WHSOQ3ZszNqHEOjPJLlzWnY0Ps/KCbIFLmx1yLG/2YNcY67B9rK3LvwrytmjLmdI5cF3mS2YxUh4tGbDNiREudUZWXvVeJj2WkzPfDuj7Hh/zHFh82IrPGuFRyxjEX+TnArs/1GH+vIpnPAFDJWOcHbUI9So67wl9vlHxFRd24lh6v4cXHyMgIXn311bF/7927F7t27cLs2bNxyimn4KabbsLtt9+OxYsXY/Hixbj99tuRTqfxpS99qdG3EhERkWmo4cXHc889h4svvnjs36tXrwYArFq1Cv/yL/+Cb37zmygUCrj22mvHvmTsqaeeQjabnbq9FhERkfethhcfF110EZzxjW6hUAhr167F2rVr381+iYiIyDSl2i4iIiISKC0+REREJFBafIiIiEigpjRqO5Vqg3G48sQ4mpHmQ5jELCM1vk3EiI1ZbSGjyCGN31nVCkkcDgCKRvzV4sjo1o0EYN1Iw9WMr2Nh1Wtjw3zfq0acL37MiFIa48kirjGjrqFVndasXGtV48w0HqNj0TsAKBuxzZr1mmTMakmjirJ1VTCOOVIwXtOR/WePAwgbUcqwMSyRMm8Ll8hrGsdljlnciOEafVwmVWiN7jD3sWbE+Nl4Wu9F+wlAyKjkGxviJw+L71rXUjMiX+DbhUtGG7lWOeMyW01N7hpsVY9m12ErNl01zqlq5h3v1jgh8ks1XOD77kgE2hWMwTz+9d/xM0VERESmgBYfIiIiEigtPkRERCRQWnyIiIhIoLT4EBERkUBp8SEiIiKBOmmjtvGjUUSSE3fPGRm7cNUfGbKiXDDawkak04qpseqvdaMCphUftaKPLDYG8NisFSmz4nxWP4ZJvDFiRN6iRqVI67isvmKVZq34pdW/ViViK/YWJnHEchN/L4tVUdiOgvoft84Bq+9jI8YJYsRfEfJvZ8W+Q9brme/Fm9j5aFaHNvreirhOZh+tqDs7pwB7PFmbNTet9zI72Djm2LB/O2s/LOGqUYXWKsxMVIyIqzUu1nlqjQurYMy+IuHt9qPWeNFgAHxO16zIetV/YDWjiu+E13jHzxQRERGZAlp8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERAJ10qZdUn1AxHuXuZUY8d92bN2RbqZWjN6x74An72X1tnHHdN0oXmUlP0Ik/TPVBZnefE2yj1ZdNis9Y9zJzu4SB4AwOQccSVsAAIxzwCqmZxUOY3e5T7ZwWMhK/1ipCiJaNBon8XoATxoBMI+NsdJE1nlfS/K2ySTArDGb/HaNd7JZ6NJ4OZZcsc6pWHFyJ4F1bWEJFDZn3w4t4Ak7FcLOKytRZqVnrPcyscST8XvCSu1ZbeZukOtwot+6IPkfthIyx9MnHyIiIhIoLT5EREQkUFp8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERALVcNT2mWeewfe//33s3LkThw4dwqZNm3D55ZePtV999dXYsGHDuG2WL1+O7du3N/Q+yWN1RGITM0D1SOPxKhbBBYBy8+RiuJVmo6hR1N9WT1oViHhTuGCsEY2MXWyIHJvxXlZcMpbn2yUG/dlHKzJrxSUtlQzvDxbDrZMxAeyCf9U0b6tYReLY201yuW/FIq0iYLFhso0xJ6y4ZLTA89FW3JONddUo5sWKFQJ2BLpu5CJZjLGStc6PyUWIXazxCKk110POGLNRviOxEVJs0YjVp/r55AxXJhmDJ9eCSuadFyMbt12at1nnFYtiW+dUNTXJGLnRxs5vViD17USN67N1XWdfk2AVVmUFIWtlq4rrca//jp/5v0ZHR3H22Wfjnnvuoc/5zGc+g0OHDo39PPHEE42+jYiIiExTDX/ysXLlSqxcudJ8TiKRQEdHx6R3SkRERKavE3LPx9atW9HW1obTTjsNX/3qV9Hb20ufWyqVMDQ0NO5HREREpq8pX3ysXLkSDz30EJ5++mn84Ac/wI4dO3DJJZegVPL/gbGrqwu5XG7sZ8GCBVO9SyIiInISmfLaLlddddXYfy9duhTLli3DwoUL8fjjj+PKK6+c8Pw1a9Zg9erVY/8eGhrSAkRERGQaO+GF5To7O7Fw4ULs2bPH255IJJBIGLcZi4iIyLRywhcffX196O7uRmdnZ0PbxYbqiHqitlZl2GKLP/9jRSmtCpjlFh4bqmWMyGHKn1EKRfg29QqPm9WMmGhsgG/HIpgsegcAyWO8LT7Ms1exUX8WNFyaXJ62mualIq0iksVZ/v4oZ/k5UM5ZUVveH9UU34960j/WrNIwAIQrvM2KUoaMyCSL0UWNqqWJAT5mkQaidL8vVCdxzyI/f0tkPgNvE482xqWS8+9HtcXKMhsxxTTfrl4zIr8FciEzKtdaJXRjo3yz+LB//1NH+ThHR3h+O1I0jjnBL9C1JBlPI8XK5jNgx3orGf6a5Vn+c9j6KgSXNK5jZSMebcVm6/7tIkbFaRadB/g4A0Asz+dtbMTfFjGu3VUyltXKO7/eN7z4GBkZwauvvjr2771792LXrl2YPXs2Zs+ejbVr1+Lzn/88Ojs7sW/fPtx8881obW3FFVdc0ehbiYiIyDTU8OLjueeew8UXXzz277fu11i1ahXWrVuH3bt348EHH8TAwAA6Oztx8cUX4yc/+Qmy2ezU7bWIiIi8bzW8+LjooovgjG/a27x587vaIREREZneVNtFREREAqXFh4iIiARKiw8REREJ1AmP2k5WYqCMaHTi2qjczIOWrHJfadbkopS1WTxSFsvwkpvJpD+mVi7z7q5HeRSqUonTNhjJR1ad1IrlZQ7y3GakwON34VGyXdWIXsWM/kjwiF01ZrQl/WNdzRhRW6NCcWW2sf9GBJNxxlCGRvhxhepGDNeqQksitakjfCzjfTzrFzKidC5pXE7IfWL1GM/Fhoxzu2KMpxmPzvr3P9ZsnPcRI34e49eIQpFfq+pF0ldGRV4rZpkY4J3FIrWJo/zEscY5VOTnTsi4xqHu/z6ncjP/voNawpi3LfytSq3G1xo0+8csmuJjGYny/igX+Di7Ud4f7KsQrOq01jgnjzX+VQgAEBklO2Lc2xkm52+4apSiPv657/iZIiIiIlNAiw8REREJlBYfIiIiEigtPkRERCRQWnyIiIhIoLT4EBERkUCdtFHbcKmKcNUTAWrhsaZKkz+WVW7m71OZwyNIiSyP33XOGuIvSpRZVUcA/cO8DGO1YlSzHOZRtASpUNv0Bo9DRYf5MYdGjXKLA/4cYKgpTTdxUaMir1FRsZbgedUKqV5bnGvEaVt5dDBhRDDDRtS2WvWPWWWAxwrtqra0CfFBvh+ZHv9YW1VLw8M86xca5fFMl+XnsEv7Y5ahmhF1N2KWZtXSOfzcSbf5O7Ilw48rHuGvN1T0HxcAO1dNKt5G8/yYk/28r9K9xnXssP+Yw0NGprPMzw8rIo8Iv1a5uL+tZlQotr4moTiXx07dbH6Ny5DreiLGj7laM+L9VaPyrhGdjhb8bdY4J/v5uZgwIvKRY/wCEiJfh+CscW72R+TDtXde1VaffIiIiEigtPgQERGRQGnxISIiIoHS4kNEREQCpcWHiIiIBOqkTbu4cAjOc+d0Jc3XS6UW/93Dpdn8ruhYbnKJljlJfvdwte6/+/lwvoluU2aFpgBERvgxxwdpE1Lkzuh4P7+zP3yUv6ArGmkXdnd8mqc76kmeBqhmeKqpOIv3R7HVf6d4pYXfhZ1q4ceVjBtFtEL8rvRS2L//FaMWXXSU3xkfHzKSDkeMAohkrMN9/Nx2BSPRQluAUIUnP+pR/93x5Rw/79l8BoDSHGNOt/A5PafJn/BIx4wEWJi/V6HCz1MY50d0wH+NSBzjL5c8xvcj3s/PYTanXYn3k1VULAQ+p51Z9NHfVmg15rORUkMr3//W2SO0LZc0rmNEf54XQKyTZBtgX7tZgi0xaBSP6+UJpUivce3OG8mmuP86zGcfUCfjXA/x8T+ePvkQERGRQGnxISIiIoHS4kNEREQCpcWHiIiIBEqLDxEREQmUFh8iIiISqIaitl1dXXj00Ufx8ssvI5VK4fzzz8edd96J008/few5zjncdtttuO+++3Ds2DEsX74c9957L5YsWdLQjtVTUdSjE3ev1Nx4LCs0h0eyFrQO0LbTc718Bw1HS/6qV84oMlQf5pG9hBHBtGJZiX5/fDB8zF8EDgDcCI+omSIkYsUeB1Br4lHb0Q4jajvHiGDO9cdO03N5NHqWUVQsEeUx1mNG/K5U9O9/uGAUCTRqFaaP8KhwsocfG41ZGrFpV+BtIc+c/F0jH5da0r9dMWdE52fxt6rn+LjMaeaxwnlN/v6Ih/nrlev8mI+CV7irjPLzOzXi76t0rzGfj/HYd+SoEZ0mY10f4nM9lOSx6VCOV+qsZfgxF+aSc6CVbmIWfWyfw69jpzTzzHI05O/jYo2Pcz/4XHdFfo2zCgWmj/j3I9PN57M1zvV+fswhq0gciVVbhT/rJDZdN4rsHa+hTz62bduG6667Dtu3b8eWLVtQrVaxYsUKjI7+rrO+973v4a677sI999yDHTt2oKOjA5deeimGh/mJIiIiIjNHQ598PPnkk+P+vX79erS1tWHnzp345Cc/Cecc7r77btxyyy248sorAQAbNmxAe3s7Hn74YXzta1+buj0XERGR96V3dc/H4OCbH2HOnj0bALB371709PRgxYoVY89JJBK48MIL8eyzz3pfo1QqYWhoaNyPiIiITF+TXnw457B69WpccMEFWLp0KQCgp6cHANDe3j7uue3t7WNtx+vq6kIulxv7WbBgwWR3SURERN4HJr34uP766/HCCy/gX//1Xye0hY678cw5N+Gxt6xZswaDg4NjP93d3ZPdJREREXkfmFRhuRtuuAGPPfYYnnnmGcyfP3/s8Y6ODgBvfgLS2dk59nhvb++ET0PekkgkkEjwO6tFRERkemlo8eGcww033IBNmzZh69atWLRo0bj2RYsWoaOjA1u2bME555wDACiXy9i2bRvuvPPOhnasMCeOaGxibKvQxqNL1Tn+aOkHWnm1PytOGw3xeGPYqFiZr/rjZkN5Xg0ykp9c5dqmg0Y1zn5/ZMuN8ihizYjfheM8/hrKZr2P142qtqXZVpzWilQbVR9b/bHZBbMG+DYRK2bJo2PxKD8/6nX/eRofMqK2eX5OxQf5PoaH+HjWjw34G4xYrKvw9wpleLS03pymbcW5/v/BKFuVa9t4/zbP4XFEK2aZjfpj93UjBj9U5efp4CiPYEYG+OU12ecf69RR3vfxXj433QjvD1f0H3PIms8JHpl1Kf4/i2ycASDf5j/3S618nJvm8HP71JajtC0V4RHdBIlVvzrMM78jxrU7eoyPc7qHz+lUr38fw0NGVem8UZG3xvvRhfl1h/1Vot7Ez+0iuXZXK3wfjtfQ4uO6667Dww8/jH/7t39DNpsdu48jl8shlUohFArhpptuwu23347Fixdj8eLFuP3225FOp/GlL32pkbcSERGRaaqhxce6desAABdddNG4x9evX4+rr74aAPDNb34ThUIB11577diXjD311FPIkv8zFhERkZml4T+7vJ1QKIS1a9di7dq1k90nERERmcZU20VEREQCpcWHiIiIBEqLDxEREQnUpL7nIwj5tggi8Ykxx9Jsft9JS6s/isYqWQLA7JhRQZBUPwSA1wuzadvhEf/NtYVjPLqU7uPrwMwhHl+KjvBIGfrJcdf4cYUzPC6JOt8OuSbvw5VW/nr5Nh5jLZAKxQDg2niV4o4W/9fztyb5OFtx2uEKjw4OjBiVLvv9UcUo3w2kjvL+jR8xNswb0TwSv2PxOgAIN/E4LebyUrPlObw/CiQ6XTTmc3QOP64P5PicTkd5/DwV8bf1l/kx9xeNCPEIj6Smj/I5zaoUxwb4uR0a4udA3ahEDBKzNCsUt/DKtaUO/1wHgOEP8NfMd5J48TweIV44i8emO5P8HLCi04dL/mMbKfO5XjYqFKcHjcq1vfzaHSn553vIiNO6uhFljRnR6SSPCrtZ/v4odvDzfrTdf82slU9QVVsRERGRd0uLDxEREQmUFh8iIiISKC0+REREJFBafIiIiEigtPgQERGRQJ20UdtiKxDxpYPm8RhSZ7M/ZrkgxeNalnyNx6tGSeVaABgc8cearCqXiQGjoukQj1eFCjxqGyIRu3rVqFpqxO9CzTxiV2nPeR8fncf7Kd/OI2rluXwfZ+V4pcs5JFIbNmLT4RBfg+crPL5WKfG+ig37X5NVMwWAVC+PWaJqVKy0xpNEakNxPi6hHK/DVMnxOO3IJMa6PI/HYtub+TgvyPA53RTh/Viq+8fzSIGf24f6/Oc2AESO8mO2qlGzOR3OG9F5a5zTVkSenDtzeGy62sr7wxrnQjvfDdfpv3ZbFadPMcaZVacF7Kht1fnn5tEhHreOGJVrE/1GXDzP522sj0SnrXE25q2ptYU2VWb7z518W+Ox6Xrx7UuwvEWffIiIiEigtPgQERGRQGnxISIiIoHS4kNEREQCpcWHiIiIBOqkTbuUOisIpyYWqWmfNUy36Uz50y7WXdEVxwvhHC7xu/73DfLCcpUBf9ol3c/vwE4M8DRGdNS4+7nE0wJI+O+MDkeM4j8Rvh6tzOPHXGjzF2XKt/PXK7fwY8608SJabU28ENVkjBjF44byRkGmfr5d/Jh/rBNGcgl1405xo6ifVTSKtsV5iqfaxouKDX+Qp13ynfz8Lrb59z+V4+k1qyCkZajK97Gv5L+z/+AQP+baUT7OmSP8mFN9jc9pcz4n+X6EjGKRLus/5nIbT7RYKbXRD/BjLrXza9W8Vv94dqb9120AyEb5+ZEM82RQd5EneV4bmON9vHSMzyOr8GdslM/b2AjvDziynVEgDmHe92jiiaeKMdYjH/CfV6PzGk8j1gvG9e04+uRDREREAqXFh4iIiARKiw8REREJlBYfIiIiEigtPkRERCRQWnyIiIhIoBqK2nZ1deHRRx/Fyy+/jFQqhfPPPx933nknTj/99LHnXH311diwYcO47ZYvX47t27c3tGPpOQVE0hPjY7NTvNhUNOyP+dSMNdbREo8gvT7Mo6X9A3y72DF/lNUsNDVsRCmtCKYRmayTqK2l1szjfEUSpwWAkU7/MY+cYhxXGy8Alk3xiF1LokDb4iRWPVDmMbRDwzxSnR/i8bvkMX5exYf8YxYp8bF0ESNGRwrEAYBLG3HgpP/8qJPHAaDQwV+vONuITuf4sYVm+8e6rZnHptNRHqWsk+JgADBc5edp95A/gpkf5dtER4xxHjSKihX5uU9jllEeg3dsGwDOiFnWsv5jK8zl58Boh1FscR6PU6bn8oh8Z8Yfqc3F+Hy24rS9ZT5v94/wa/fQKCn8Ocz7PsIvR4iPGjFnY0q7uP/Xr8vx3y2I8nEptfKI+Whn40Ufi3P5cSVm+8esljc66jgNffKxbds2XHfdddi+fTu2bNmCarWKFStWYHR0/An3mc98BocOHRr7eeKJJxp5GxEREZnGGvrk48knnxz37/Xr16OtrQ07d+7EJz/5ybHHE4kEOjo6pmYPRUREZFp5V/d8DA6++XeE2bPHf8S1detWtLW14bTTTsNXv/pV9Pb20tcolUoYGhoa9yMiIiLT16QXH845rF69GhdccAGWLl069vjKlSvx0EMP4emnn8YPfvAD7NixA5dccglKJf/ffLu6upDL5cZ+FixYMNldEhERkfeBSdd2uf766/HCCy/gF7/4xbjHr7rqqrH/Xrp0KZYtW4aFCxfi8ccfx5VXXjnhddasWYPVq1eP/XtoaEgLEBERkWlsUouPG264AY899hieeeYZzJ8/33xuZ2cnFi5ciD179njbE4kEEgl+p7mIiIhMLw0tPpxzuOGGG7Bp0yZs3boVixYtettt+vr60N3djc7OzoZ2LJ0oI5KYGAFqivF4JovfvVFoodscyvNqlr1DPPJUNypdJgdZRVMjkmUUmq0neGO12ahoSqJ5lWYeu6pk+HuNzONt+Q7/e9UzPJbX0sQjdrkEj2yFwSOHQ2V/3KyvwKOII0bl2kg/jyNa8bsoObRwxYijGpHq6hx+LloR3VraP8WLs/jUL7Tyv8YW2mgTym28gmdTk3/exiL8/Biq8HE5WszQtsES325whLQZ8zk2zPs3mudjZseq/X1cb+L7XjdilvWkMTfb/fM938Zfr9hqxPtbePx1TpP1VQj+61+1zve9r8zP+54ij9oeNuLzpX7/NSJe4ONsVa615rQzxqxKItCW8ix+7S5njXh0h1GJeJZ//90cXmF5VtY/zrUw//18vIbu+bjuuuvw4x//GA8//DCy2Sx6enrQ09ODQuHNK+3IyAi+8Y1v4L/+67+wb98+bN26FZdddhlaW1txxRVXNPJWIiIiMk019MnHunXrAAAXXXTRuMfXr1+Pq6++GpFIBLt378aDDz6IgYEBdHZ24uKLL8ZPfvITZLN8JSoiIiIzR8N/drGkUils3rz5Xe2QiIiITG+q7SIiIiKB0uJDREREAqXFh4iIiARq0l8ydqLFwjVEPTG8Yo1HH/NVfwypUOXb9AzyG2ELJJIFAHEjfhchMUtnVCYNGQUw6zG+RrRes5r2R9jKzUbEzqhaWpxDm1CZ5T+ARAvPoybjPLLnjHKQR4o8fjdS9p8D/UM8mlke5vG1WMmKWdImqpoy1vuOn6cuzPej7omkv6WUZecA38Ya53KLUaU4wWOzLGZZqvJLUKXGI5hFY7tjwzxWXe0jMcshK2ZJmxCu8vvgakmrpKn/2KoZflzW61WM84pFpwttfN8rc3lsOp3hccqWJI/PM6M1Pv8Gy0bcusCvA8Mj/NodKvv7MZo3KkeHjL5KG9fnsPErlrxdPWrMzRb+XpUm4zow29j/2f6xTqT49TlG5nOIPO6jTz5EREQkUFp8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERAKlxYeIiIgE6qSN2oZDDmFPvGm4zCsBFir+qOKQUbW0MMxfLzLIuydiVEBkqSxSdBcAUIsbMTojflc3RrCS9r9mzYpmzuKvV5rLo5Ro8se1YjG+TdWIUlqVSSs1Iw5c9p8D5REe54v2WZVrJxePrrKkX8iocmlFM42mqjGeLFJbbuGvV2k2qi+n+HhGInb5BZ/+UR6LtZRLRkR3kM/p6LC//635bFUttWKR1SQf61rMGFDCqlpaamm8amnFiE1HMjxmmUnyaqdlY07XnX9c+owKxSw6D9iR6lqenx+xUX8/hnm62Jx/7DoLAGUSdQcAVpzbGdf0cs76agW+XTVrROTJnLau3SzqXqsZvyOOo08+REREJFBafIiIiEigtPgQERGRQGnxISIiIoHS4kNEREQCpcWHiIiIBOqkjdrmKzFEPLHJQolHr8olf2TSil1FBnibVeUwwpNoVCVjVCaNW5UR+WvWjcgei16Vszw6WM3wNhfjca1o3B+xsqrTjhSMSGSUR7ZKVsySRKdDJd6JVmQ2wlOFqPHdB+r+466kjf6NTu4cqPFUMt1Ha5zrcd5mVa6lGXMAo3n/jtSN86NWsDKHvEOiwzzeGCXRaWucrYrC1dTkqtqyatR1I3JvRSlLLcacbvKf4C5pxab5pKgaUff+At/JCKl4OlLkE6lU5DH4yihvC4/ycyBEDrvOX86+dhvbkeLFAIAqmbfm3Iwa5xsZZ+Btxjrq365a5TvPahfXyqpqKyIiIicpLT5EREQkUFp8iIiISKC0+BAREZFAafEhIiIigWoo7bJu3TqsW7cO+/btAwAsWbIE3/nOd7By5UoAgHMOt912G+677z4cO3YMy5cvx7333oslS5Y0vGMDxzIIFyfeDuyKxu3DJGEQyfM1VqTIXy5SMt7K6DlHiopZhYuqRnEi625qq63a5L8z2rpjup4y7laO87Za2T8u+crk1rdFz9iPMRISLNUSNQqHRUpWKoHvhnUO1HKsZXLF46z3qhoJGvaatbRxZ7xxfsBIoNRKxp34ZG6GylaizCr2ZmxnjrX/cSvxZM0xkNQKAFR4vTQ6LrQgIYBqhu+kmbhI+LcLGfPZUqrwc6BojGel7N/OSsTVhvmBhWrGdcC41rLtaknjumilCpv4e1nXDzZv68lJzs0wbwsbKbU6uUaXyZwFeGHHev6dF5ds6DfD/Pnzcccdd+C5557Dc889h0suuQSf+9zn8OKLLwIAvve97+Guu+7CPffcgx07dqCjowOXXnophoeHG3kbERERmcYaWnxcdtll+NM//VOcdtppOO200/Dd734XTU1N2L59O5xzuPvuu3HLLbfgyiuvxNKlS7Fhwwbk83k8/PDDJ2r/RURE5H1m0vd81Go1bNy4EaOjozjvvPOwd+9e9PT0YMWKFWPPSSQSuPDCC/Hss8/S1ymVShgaGhr3IyIiItNXw4uP3bt3o6mpCYlEAtdccw02bdqEM844Az09PQCA9vb2cc9vb28fa/Pp6upCLpcb+1mwYEGjuyQiIiLvIw0vPk4//XTs2rUL27dvx9e//nWsWrUKL7300lh76LgbsJxzEx77fWvWrMHg4ODYT3d3d6O7JCIiIu8jDdd2icfjOPXUUwEAy5Ytw44dO/DDH/4Q3/rWtwAAPT096OzsHHt+b2/vhE9Dfl8ikUAiYRXJEBERkenkXReWc86hVCph0aJF6OjowJYtW3DOOecAAMrlMrZt24Y777yz8RceiQG1iTGr6IhRUIpE7OoRHv8xY5bW50JGoqhO1lI1I0bnjJiUVTjMis26BGkzaoNZQkaBPkf62IzDGbE8UoPqze2MmCWLTFrRTGssLWWjmBc7rypG0SiziB2vp2hG82j/W+NitFnC1niS6KMViQyb+2jtB29j54fVv9Z1wFlxa6NAnyMJ0rrxTQJWPBrGNY7Fel2VH1i5bBTw5O8EWNF6soshI9JpzVtzThtdVWkm0WMzbm3MsarxNQnGVxew42bXUgBmnBbGPtatIo1M0SjGyeZm8Z1fTBvao5tvvhkrV67EggULMDw8jI0bN2Lr1q148sknEQqFcNNNN+H222/H4sWLsXjxYtx+++1Ip9P40pe+1MjbiIiIyDTW0OLj8OHD+PKXv4xDhw4hl8vhrLPOwpNPPolLL70UAPDNb34ThUIB11577diXjD311FPIZrMnZOdFRETk/aehxccDDzxgtodCIaxduxZr1659N/skIiIi05hqu4iIiEigtPgQERGRQGnxISIiIoF611HbEyV+JIJIwpM7M9JV4Yr/8agRhaoaMVbrvawKiLSSKIu+AnBWLC/U+HuZbUacz4q/WpGyKKkma0UirWibFaOz4plmBVLCqhRZmjO5vqeVZo33qmV5m9X3VkybxoiNcpvRvBFvnEScFuDVo61qvdbrWeeOVUmUxtaNLqwa8WgzajuZ6qSxScYsrXg0qyhsVPu24uyTja2zyLJVWdxZRcyta6ZVdZz0vXXeOCPGChLdBezLM9h1vWTlra2LDt8uXDYqu5OvqJjMXK+R3wPe13jHzxQRERGZAlp8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERAKlxYeIiIgE6uSN2g4BEV91WKtoHllKWRUrrShoJcffzIp51VPkRY2lXiTNc4U1IyYFozIlSCQubMShrCq/ESuCSWLOZlxykjHcGqkabO0Hjb4CqDZNLk5bT1ulVclrTrJSsjNKq1pVaKMkThkpWtvw/bBizpESb6OsrnfGHIvy/a808ddk54EV0a5ljHFmcx1AJMHb6mzeWufAZCqTAogO+SOYYWO8rEqz5jhb12f6FQTGyxnRdCsba85NMi6ROL/oxOL8xK8bfRWL8f0o5P2/mJzxeta4xAYaj9MCvJq2dX6EyWHVzJLHx73GO3+qiIiIyLunxYeIiIgESosPERERCZQWHyIiIhIoLT5EREQkUFp8iIiISKBO2qhtYtAh4omzWlVoXdgfJzJjdGYEk8ekQkaUK57y5z1DRnXaaJTHvMphPkyVUX5woYo/YhcbnlxFU6v6ZDTvPzYrMjuZ2PSb2zUes6xbcT6rOqYRl4wa8egoidgZ6VEzslcp8xM/bMVmyVhHC3w/4kNGdV0jwRgfNfqRjGc5Y1TrtSrGpnhb3YjWs2rUdSNOG87wcZ5sBDOa9L9fmcQv306sn1c0ZfF56xyIjk5qN+xK4JM4NFqFGEA9YVSTTfFxiZAKtc1ZnjHPxEmGH0Cxyk/UoVF+AM6qUEvE+vmFMTZixHCN8YyUyLV7El+FEC5bF/XjnvuOnykiIiIyBbT4EBERkUBp8SEiIiKB0uJDREREAqXFh4iIiASqobTLunXrsG7dOuzbtw8AsGTJEnznO9/BypUrAQBXX301NmzYMG6b5cuXY/v27Q3vWKTsEPFFA4ykQ2mW//Fqmr9PJcvvmA5n+R3OqQyvoJOI+bcjYZy3xQoQAQCMonMsuZLoa7xAHADEjDRDtOhvC1f5NrU4349qcnKd5Uj4p5YwkhhGaiVhjHMqwTsrHPafV2Xjzvj8iBXJ4U0J4w54llpI9fHzno3lm69nxZe4esI/ni7M973cxM+BmpFgs+Z0vcU/1sksr6LVlOJtETLOAFCq8LEulPwHYCUgooM80ZI4ZmxHzoFowRhnkoAA7KJ+xuUZrrnxNKKVnnFG4jDdxMcsk/TP6dkpnnapGYUda8ZBh1mBSfBCgZFR/l7WOMeH+XvFRoy2vP8cZulRAHDkVKxW3vn1oaHFx/z583HHHXfg1FNPBQBs2LABn/vc5/D8889jyZIlAIDPfOYzWL9+/dg28fjkomMiIiIyPTW0+LjsssvG/fu73/0u1q1bh+3bt48tPhKJBDo6OqZuD0VERGRamfQ9H7VaDRs3bsTo6CjOO++8sce3bt2KtrY2nHbaafjqV7+K3t5e83VKpRKGhobG/YiIiMj01fDiY/fu3WhqakIikcA111yDTZs24YwzzgAArFy5Eg899BCefvpp/OAHP8COHTtwySWXoFTif3/r6upCLpcb+1mwYMHkj0ZEREROeg1/vfrpp5+OXbt2YWBgAI888ghWrVqFbdu24YwzzsBVV1019rylS5di2bJlWLhwIR5//HFceeWV3tdbs2YNVq9ePfbvoaEhLUBERESmsYYXH/F4fOyG02XLlmHHjh344Q9/iH/6p3+a8NzOzk4sXLgQe/bsoa+XSCSQSBh3+YuIiMi08q4Lyznn6J9V+vr60N3djc7OzoZfN1JyiNYnxoPMmBeL/6SMuJNRUKolx6NXczK8LRHxx/lGKzz5c2SYVESDHb+LjPL4XeqwfzurcFhykEelIkUjnpn392Mtwf+yF8rwfa/H+DFbhcOqpHCYFcvL5HjFvJY0r74VNWKWlbr/uAslvvP1EZ45jBuFw+IDtAnJY/59TB7j/REb5hFiK35nqZDiiFY0s5oyothNRoQ7y4+NRWpbmvh8TpPoPADUrei/EbUtkVh15BjfJj5oFA4zYpYsVs0ilgDMaLc1p615WyFFBOvGbyGruGd6Fp+bWSMePT874H08bBT+rJL5DACDRV48rkqKewJAdMjfxq7bAJA6wvcxMWTM6RGjcGLVfx7UEnzfy8287Z1qaPFx8803Y+XKlViwYAGGh4exceNGbN26FU8++SRGRkawdu1afP7zn0dnZyf27duHm2++Ga2trbjiiive9Y6KiIjI9NDQ4uPw4cP48pe/jEOHDiGXy+Gss87Ck08+iUsvvRSFQgG7d+/Ggw8+iIGBAXR2duLiiy/GT37yE2Sz2RO1/yIiIvI+09Di44EHHqBtqVQKmzdvftc7JCIiItObaruIiIhIoLT4EBERkUC967TLVHP/W0yuVvEnEGpGIbVaidxNbRTKqhf4ney1PL9jugreFo747yyuVow79PM86VA3CkDBODbWH7WycVe3URjIWUWDqiTtEuHjZd0JPplxBvhY1wu8eJw5zo63wUi7sLvja3kjKVAw2opWXxkF2MiYVSv87vdQderTLlWS/KiVjeMqWf1hzWlrrMl1JWScA5NMu9Ss8SRFxUJFfkm2+mMyczpkzWcr7WIUA6wbleDYnK4ZU2yy87ZW522VsL+w3GTTLtZ+1I357oqsP6w+tMbZmtOTSLtE+Nxk1+63fm87X1HY4/fJvZNnBejAgQP6kjEREZH3qe7ubsyfP998zkm3+KjX6zh48CCy2SxCodDYN552d3ejubn5vd6995z6Yzz1x3jqj/HUH+OpP8ZTf/zOVPSFcw7Dw8OYN28ewsYnZMBJ+GeXcDjsXTE1NzfP+JPj96k/xlN/jKf+GE/9MZ76Yzz1x++8277I5XLv6Hm64VREREQCpcWHiIiIBOqkX3wkEgnceuutKj73v9Qf46k/xlN/jKf+GE/9MZ7643eC7ouT7oZTERERmd5O+k8+REREZHrR4kNEREQCpcWHiIiIBEqLDxEREQnUSb34+NGPfoRFixYhmUzi3HPPxX/+53++17sUmGeeeQaXXXYZ5s2bh1AohJ/+9Kfj2p1zWLt2LebNm4dUKoWLLroIL7744nuzsydYV1cX/vAP/xDZbBZtbW24/PLL8corr4x7zkzqj3Xr1uGss84a+zKg8847D//xH/8x1j6T+uJ4XV1dCIVCuOmmm8Yem2n9sXbtWoRCoXE/HR0dY+0zrT8A4I033sBf/dVfYc6cOUin0/iDP/gD7Ny5c6x9JvXJBz/4wQnnRygUwnXXXQcgwL5wJ6mNGze6WCzm7r//fvfSSy+5G2+80WUyGbd///73etcC8cQTT7hbbrnFPfLIIw6A27Rp07j2O+64w2WzWffII4+43bt3u6uuusp1dna6oaGh92aHT6A/+ZM/cevXr3e/+c1v3K5du9xnP/tZd8opp7iRkZGx58yk/njsscfc448/7l555RX3yiuvuJtvvtnFYjH3m9/8xjk3s/ri9/3qV79yH/zgB91ZZ53lbrzxxrHHZ1p/3HrrrW7JkiXu0KFDYz+9vb1j7TOtP/r7+93ChQvd1Vdf7X75y1+6vXv3up/97Gfu1VdfHXvOTOqT3t7ecefGli1bHAD385//3DkXXF+ctIuPP/qjP3LXXHPNuMc+8pGPuG9/+9vv0R69d45ffNTrddfR0eHuuOOOsceKxaLL5XLuH//xH9+DPQxWb2+vA+C2bdvmnFN/OOfcrFmz3D//8z/P2L4YHh52ixcvdlu2bHEXXnjh2OJjJvbHrbfe6s4++2xv20zsj29961vuggsuoO0zsU9+34033ug+/OEPu3q9HmhfnJR/dimXy9i5cydWrFgx7vEVK1bg2WeffY/26uSxd+9e9PT0jOufRCKBCy+8cEb0z+DgIABg9uzZAGZ2f9RqNWzcuBGjo6M477zzZmxfXHfddfjsZz+LT3/60+Men6n9sWfPHsybNw+LFi3CF77wBbz22msAZmZ/PPbYY1i2bBn+4i/+Am1tbTjnnHNw//33j7XPxD55S7lcxo9//GN85StfQSgUCrQvTsrFx9GjR1Gr1dDe3j7u8fb2dvT09LxHe3XyeKsPZmL/OOewevVqXHDBBVi6dCmAmdkfu3fvRlNTExKJBK655hps2rQJZ5xxxozsi40bN+LXv/41urq6JrTNxP5Yvnw5HnzwQWzevBn3338/enp6cP7556Ovr29G9sdrr72GdevWYfHixdi8eTOuueYa/O3f/i0efPBBADPzHHnLT3/6UwwMDODqq68GEGxfnHRVbX9fKBQa92/n3ITHZrKZ2D/XX389XnjhBfziF7+Y0DaT+uP000/Hrl27MDAwgEceeQSrVq3Ctm3bxtpnSl90d3fjxhtvxFNPPYVkMkmfN1P6AwBWrlw59t9nnnkmzjvvPHz4wx/Ghg0b8PGPfxzAzOqPer2OZcuW4fbbbwcAnHPOOXjxxRexbt06/PVf//XY82ZSn7zlgQcewMqVKzFv3rxxjwfRFyflJx+tra2IRCITVlq9vb0TVmQz0Vt3rs+0/rnhhhvw2GOP4ec//znmz58/9vhM7I94PI5TTz0Vy5YtQ1dXF84++2z88Ic/nHF9sXPnTvT29uLcc89FNBpFNBrFtm3b8A//8A+IRqNjxzxT+sMnk8ngzDPPxJ49e2bc+QEAnZ2dOOOMM8Y99tGPfhSvv/46gJl5/QCA/fv342c/+xn+5m/+ZuyxIPvipFx8xONxnHvuudiyZcu4x7ds2YLzzz//Pdqrk8eiRYvQ0dExrn/K5TK2bds2LfvHOYfrr78ejz76KJ5++mksWrRoXPtM6w8f5xxKpdKM64tPfepT2L17N3bt2jX2s2zZMvzlX/4ldu3ahQ996EMzqj98SqUSfvvb36Kzs3PGnR8A8Md//McTovn//d//jYULFwKYudeP9evXo62tDZ/97GfHHgu0L6b09tUp9FbU9oEHHnAvvfSSu+mmm1wmk3H79u17r3ctEMPDw+755593zz//vAPg7rrrLvf888+PRY3vuOMOl8vl3KOPPup2797tvvjFL07baNjXv/51l8vl3NatW8dFxPL5/NhzZlJ/rFmzxj3zzDNu79697oUXXnA333yzC4fD7qmnnnLOzay+8Pn9tItzM68//u7v/s5t3brVvfbaa2779u3uz/7sz1w2mx27ds60/vjVr37lotGo++53v+v27NnjHnroIZdOp92Pf/zjsefMtD6p1WrulFNOcd/61rcmtAXVFyft4sM55+699163cOFCF4/H3cc+9rGxaOVM8POf/9wBmPCzatUq59yb8bBbb73VdXR0uEQi4T75yU+63bt3v7c7fYL4+gGAW79+/dhzZlJ/fOUrXxmbF3PnznWf+tSnxhYezs2svvA5fvEx0/rjre9liMVibt68ee7KK690L7744lj7TOsP55z793//d7d06VKXSCTcRz7yEXffffeNa59pfbJ582YHwL3yyisT2oLqi5Bzzk3tZykiIiIi3El5z4eIiIhMX1p8iIiISKC0+BAREZFAafEhIiIigdLiQ0RERAKlxYeIiIgESosPERERCZQWHyIiIhIoLT5EREQkUFp8iIiISKC0+BAREZFAafEhIiIigfr/AZx5IyQItUwBAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "img = np.asarray(Image.open(r\"C:\\Users\\lociu\\Desktop\\PriyaM\\Data\\20230522_sample_minmax\\priya_samplefile.tif\"))\n",
    "plt.imshow(img)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "a4514aa8",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[(0, 25), (0, 60), (1, 42), (2, 6), (2, 24), (3, 41), (9, 72), (22, 5), (23, 40), (24, 58), (25, 21), (27, 72), (32, 13), (36, 7), (36, 25), (36, 43), (36, 61), (36, 72)]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAEpCAYAAAAzlkD8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAdW0lEQVR4nO3de3CU1f3H8c9C4DHgZivFZHeHGFMNKlcVbAilErBE87MMF9vx0loyTh2RQGGoA4J/kOmMLOrIaIeaVnQojNj4B9eOSIiDCXUoGigZI1gah6CxJKZSsxtSuuFyfn+0bF2TIBs2J9ns+zVzZrrnnN398t3EfubJ8zzrMsYYAQAAWDKgtwsAAADJhfABAACsInwAAACrCB8AAMAqwgcAALCK8AEAAKwifAAAAKsIHwAAwCrCBwAAsIrwAQAArErpqRd+6aWX9Nxzz6mxsVGjR4/WCy+8oO9///vf+LwLFy7o5MmTcrvdcrlcPVUeAACII2OMWltb5ff7NWDANxzbMD2grKzMDBo0yKxfv94cPXrULF682AwdOtR88skn3/jchoYGI4nBYDAYDEYCjoaGhm/8/3qXMfH/Yrnc3FzdfvvtKi0tjczdcsstmj17tgKBwCWfGwwG9a1vfUtT9H9K0aB4lwbEzba/1Xa5NmfkWIuVALhcXf3e8jt75c7prN7VLrW0tMjj8Vxyb9z/7NLe3q5Dhw7pySefjJovKCjQ/v37O+wPh8MKh8ORx62trf8tbJBSXIQP9F1p7q4PK/KzC/RNXf3e8jsbB/89lHE5p0zE/YTTL774QufPn1dGRkbUfEZGhpqamjrsDwQC8ng8kZGZmRnvkgAAQB/SY1e7fD35GGM6TUMrVqxQMBiMjIaGhp4qCQAA9AFx/7PL8OHDNXDgwA5HOZqbmzscDZEkx3HkOE68ywAAAH1U3MPH4MGDNWHCBFVUVGjOnDmR+YqKCs2aNSvebwf0mrv9t/Z2CYiT8pM1Xa7xOfcvfJ59Q4/c52Pp0qV6+OGHNXHiROXl5enll1/Wp59+qvnz5/fE2wEAgATSI+Hj/vvv16lTp/SrX/1KjY2NGjNmjHbt2qWsrKyeeDsAAJBAeuwOpwsWLNCCBQt66uUBAECC4rtdAACAVYQPAABgFeEDAABY1WPnfPR3XJqXHPickwOfJWAXRz4AAIBVhA8AAGAV4QMAAFhF+AAAAFYRPgAAgFVc7dJNnB2fHPic+5eurl7icwbs4sgHAACwivABAACsInwAAACrCB8AAMAqwgcAALCK8AEAAKziUlsASYNLaoG+gSMfAADAKsIHAACwivABAACsInwAAACrCB8AAMAqwgcAALAq7uGjpKRELpcrani93ni/DQAASFA9cp+P0aNH6+233448HjhwYE+8DQAASEA9Ej5SUlI42gEAADrVI+d81NXVye/3Kzs7Ww888ICOHz/e5d5wOKxQKBQ1AABA/xX38JGbm6tNmzapvLxc69evV1NTkyZPnqxTp051uj8QCMjj8URGZmZmvEsCAAB9iMsYY3ryDdra2nTDDTdo2bJlWrp0aYf1cDiscDgceRwKhZSZmal8zVKKa1BPlgYAAOLknDmrSu1QMBhUWlraJff2+BfLDR06VGPHjlVdXV2n647jyHGcni4DAAD0ET1+n49wOKyPPvpIPp+vp98KAAAkgLiHjyeeeEJVVVWqr6/Xe++9px/96EcKhUKaN29evN8KAAAkoLj/2eWzzz7Tgw8+qC+++ELXXnutJk2apAMHDigrKyvebwUAABJQ3MNHWVlZvF8SAAD0I3y3CwAAsIrwAQAArCJ8AAAAqwgfAADAKsIHAACwivABAACsInwAAACrCB8AAMAqwgcAALCK8AEAAKwifAAAAKsIHwAAwCrCBwAAsIrwAQAArCJ8AAAAqwgfAADAqpTeLgAAACSO8pM1nc6HWi/ompGX9xoc+QAAAFYRPgAAgFWEDwAAYBXhAwAAWEX4AAAAVhE+AACAVTGHj3379mnmzJny+/1yuVzavn171LoxRiUlJfL7/UpNTVV+fr6OHDkSr3oBAEAvutt/a6djzsixl/0aMYePtrY2jR8/XuvWret0/dlnn9XatWu1bt06VVdXy+v1asaMGWptbY31rQAAQD8U803GCgsLVVhY2OmaMUYvvPCCnnrqKc2dO1eStHHjRmVkZOj111/XY489dmXVAgCAhBfXcz7q6+vV1NSkgoKCyJzjOJo6dar279/f6XPC4bBCoVDUAAAA/Vdcw0dTU5MkKSMjI2o+IyMjsvZ1gUBAHo8nMjIzM+NZEgAA6GN65GoXl8sV9dgY02HuohUrVigYDEZGQ0NDT5QEAAD6iLh+sZzX65X0nyMgPp8vMt/c3NzhaMhFjuPIcZx4lgEAAPqwuB75yM7OltfrVUVFRWSuvb1dVVVVmjx5cjzfCgAAJKiYj3ycPn1aH3/8ceRxfX29ampqNGzYMF133XVasmSJVq9erZycHOXk5Gj16tUaMmSIHnroobgWDgAAElPM4ePgwYOaNm1a5PHSpUslSfPmzdPvf/97LVu2TGfOnNGCBQv05ZdfKjc3V3v27JHb7Y5f1QAAIGG5jDGmt4v4qlAoJI/Ho3zNUoprUG+XAwAALsM5c1aV2qFgMKi0tLRL7uW7XQAAgFWEDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYRfgAAABWET4AAIBVhA8AAGAV4QMAAFhF+AAAAFYRPgAAgFWEDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYRfgAAABWET4AAIBVhA8AAGAV4QMAAFgVc/jYt2+fZs6cKb/fL5fLpe3bt0etFxUVyeVyRY1JkybFq14AAJDgYg4fbW1tGj9+vNatW9flnnvuuUeNjY2RsWvXrisqEgAA9B8psT6hsLBQhYWFl9zjOI68Xm+3iwIAAP1Xj5zzUVlZqfT0dI0cOVKPPvqompubu9wbDocVCoWiBgAA6L/iHj4KCwu1efNm7d27V88//7yqq6s1ffp0hcPhTvcHAgF5PJ7IyMzMjHdJAACgD3EZY0y3n+xyadu2bZo9e3aXexobG5WVlaWysjLNnTu3w3o4HI4KJqFQSJmZmcrXLKW4BnW3NAAAYNE5c1aV2qFgMKi0tLRL7o35nI9Y+Xw+ZWVlqa6urtN1x3HkOE5PlwEAAPqIHr/Px6lTp9TQ0CCfz9fTbwUAABJAzEc+Tp8+rY8//jjyuL6+XjU1NRo2bJiGDRumkpIS3XffffL5fDpx4oRWrlyp4cOHa86cOXEtHAAAJKaYw8fBgwc1bdq0yOOlS5dKkubNm6fS0lLV1tZq06ZNamlpkc/n07Rp0/TGG2/I7XbHr2oAAJCwYg4f+fn5utQ5quXl5VdUEAAA6N/4bhcAAGAV4QMAAFhF+AAAAFYRPgAAgFWEDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYRfgAAABWET4AAIBVhA8AAGAV4QMAAFhF+AAAAFYRPgAAgFWEDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYRfgAAABWxRQ+AoGA7rjjDrndbqWnp2v27Nk6duxY1B5jjEpKSuT3+5Wamqr8/HwdOXIkrkUDAIDEFVP4qKqqUnFxsQ4cOKCKigqdO3dOBQUFamtri+x59tlntXbtWq1bt07V1dXyer2aMWOGWltb4148AABIPC5jjOnuk//xj38oPT1dVVVVuvPOO2WMkd/v15IlS7R8+XJJUjgcVkZGhp555hk99thj3/iaoVBIHo9H+ZqlFNeg7pYGAAAsOmfOqlI7FAwGlZaWdsm9V3TORzAYlCQNGzZMklRfX6+mpiYVFBRE9jiOo6lTp2r//v2dvkY4HFYoFIoaAACg/+p2+DDGaOnSpZoyZYrGjBkjSWpqapIkZWRkRO3NyMiIrH1dIBCQx+OJjMzMzO6WBAAAEkC3w8fChQv1wQcf6A9/+EOHNZfLFfXYGNNh7qIVK1YoGAxGRkNDQ3dLAgAACSClO09atGiRdu7cqX379mnEiBGRea/XK+k/R0B8Pl9kvrm5ucPRkIscx5HjON0pAwAAJKCYjnwYY7Rw4UJt3bpVe/fuVXZ2dtR6dna2vF6vKioqInPt7e2qqqrS5MmT41MxAABIaDEd+SguLtbrr7+uHTt2yO12R87j8Hg8Sk1Nlcvl0pIlS7R69Wrl5OQoJydHq1ev1pAhQ/TQQw/1yD8AAAAklpjCR2lpqSQpPz8/an7Dhg0qKiqSJC1btkxnzpzRggUL9OWXXyo3N1d79uyR2+2OS8EAACCxXdF9PnoC9/kAACDxWLvPBwAAQKwIHwAAwCrCBwAAsKpb9/lIFuUna7pcu9t/q7U6AADoTzjyAQAArCJ8AAAAqwgfAADAKsIHAACwivABAACsInwAAACruNT2EricFkg8XCIP9H0c+QAAAFYRPgAAgFWEDwAAYBXhAwAAWEX4AAAAVnG1C4B+hStakgNXNSU2jnwAAACrCB8AAMAqwgcAALCK8AEAAKwifAAAAKsIHwAAwKqYLrUNBALaunWr/vrXvyo1NVWTJ0/WM888o5tuuimyp6ioSBs3box6Xm5urg4cOBCfioFL6OryOy69A/oXfqcTW0xHPqqqqlRcXKwDBw6ooqJC586dU0FBgdra2qL23XPPPWpsbIyMXbt2xbVoAACQuGI68rF79+6oxxs2bFB6eroOHTqkO++8MzLvOI68Xm98KgQAAP3KFZ3zEQwGJUnDhg2Lmq+srFR6erpGjhypRx99VM3NzV2+RjgcVigUihoAAKD/6nb4MMZo6dKlmjJlisaMGROZLyws1ObNm7V37149//zzqq6u1vTp0xUOhzt9nUAgII/HExmZmZndLQkAACQAlzHGdOeJxcXFevPNN/Xuu+9qxIgRXe5rbGxUVlaWysrKNHfu3A7r4XA4KpiEQiFlZmYqX7OU4hrUndKQxDjhFAB6xzlzVpXaoWAwqLS0tEvu7dYXyy1atEg7d+7Uvn37Lhk8JMnn8ykrK0t1dXWdrjuOI8dxulMGAABIQDGFD2OMFi1apG3btqmyslLZ2dnf+JxTp06poaFBPp+v20UCl4sjHADQ98V0zkdxcbFee+01vf7663K73WpqalJTU5POnDkjSTp9+rSeeOIJ/fnPf9aJEydUWVmpmTNnavjw4ZozZ06P/AMAAEBiienIR2lpqSQpPz8/an7Dhg0qKirSwIEDVVtbq02bNqmlpUU+n0/Tpk3TG2+8IbfbHbeiAQBA4or5zy6XkpqaqvLy8isqCAAA9G98twsAALCK8AEAAKwifAAAAKu6dZ8PAACQnLq6mWOo9YKuGXl5r8GRDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYxaW2AADgsnX17eHnzFlJxy/rNTjyAQAArCJ8AAAAqwgfAADAKsIHAACwivABAACsInwAAACrCB8AAMAqwgcAALCK8AEAAKwifAAAAKsIHwAAwCrCBwAAsCqm8FFaWqpx48YpLS1NaWlpysvL01tvvRVZN8aopKREfr9fqampys/P15EjR+JeNAAASFwxhY8RI0ZozZo1OnjwoA4ePKjp06dr1qxZkYDx7LPPau3atVq3bp2qq6vl9Xo1Y8YMtba29kjxAAAg8biMMeZKXmDYsGF67rnn9Mgjj8jv92vJkiVavny5JCkcDisjI0PPPPOMHnvssct6vVAoJI/Ho3zNUopr0JWUBgAALDlnzqpSOxQMBpWWlnbJvd0+5+P8+fMqKytTW1ub8vLyVF9fr6amJhUUFET2OI6jqVOnav/+/V2+TjgcVigUihoAAKD/ijl81NbW6uqrr5bjOJo/f762bdumUaNGqampSZKUkZERtT8jIyOy1plAICCPxxMZmZmZsZYEAAASSMzh46abblJNTY0OHDigxx9/XPPmzdPRo0cj6y6XK2q/MabD3FetWLFCwWAwMhoaGmItCQAAJJCUWJ8wePBg3XjjjZKkiRMnqrq6Wi+++GLkPI+mpib5fL7I/ubm5g5HQ77KcRw5jhNrGQAAIEFd8X0+jDEKh8PKzs6W1+tVRUVFZK29vV1VVVWaPHnylb4NAADoJ2I68rFy5UoVFhYqMzNTra2tKisrU2VlpXbv3i2Xy6UlS5Zo9erVysnJUU5OjlavXq0hQ4booYce6qn6AQBAgokpfHz++ed6+OGH1djYKI/Ho3Hjxmn37t2aMWOGJGnZsmU6c+aMFixYoC+//FK5ubnas2eP3G53jxQPAAASzxXf5yPeuM8HAACJx8p9PgAAALqD8AEAAKwifAAAAKtivs8HEl/5yZou1+7232qtDgBAcuLIBwAAsIrwAQAArCJ8AAAAqwgfAADAKsIHAACwivABAACs4lLbJMTltACA3sSRDwAAYBXhAwAAWEX4AAAAVhE+AACAVYQPAABgFeEDAABYRfgAAABWET4AAIBVhA8AAGAV4QMAAFhF+AAAAFYRPgAAgFUxhY/S0lKNGzdOaWlpSktLU15ent56663IelFRkVwuV9SYNGlS3IsGAACJK6ZvtR0xYoTWrFmjG2+8UZK0ceNGzZo1S4cPH9bo0aMlSffcc482bNgQec7gwYPjWC4AAEh0MYWPmTNnRj1++umnVVpaqgMHDkTCh+M48nq98asQAAD0K90+5+P8+fMqKytTW1ub8vLyIvOVlZVKT0/XyJEj9eijj6q5ufmSrxMOhxUKhaIGAADov2IOH7W1tbr66qvlOI7mz5+vbdu2adSoUZKkwsJCbd68WXv37tXzzz+v6upqTZ8+XeFwuMvXCwQC8ng8kZGZmdn9fw0AAOjzXMYYE8sT2tvb9emnn6qlpUVbtmzRK6+8oqqqqkgA+arGxkZlZWWprKxMc+fO7fT1wuFwVDgJhULKzMxUvmYpxTUoxn8OAADoDefMWVVqh4LBoNLS0i65N6ZzPqT/nEB68YTTiRMnqrq6Wi+++KJ+97vfddjr8/mUlZWlurq6Ll/PcRw5jhNrGQAAIEHFHD6+zhjT5Z9VTp06pYaGBvl8vit9m4RSfrKm0/m7/bdarQM9q6vPWeKz7k/4nJMDn/Pl66pXodYLumbk5b1GTOFj5cqVKiwsVGZmplpbW1VWVqbKykrt3r1bp0+fVklJie677z75fD6dOHFCK1eu1PDhwzVnzpxY3gYAAPRjMYWPzz//XA8//LAaGxvl8Xg0btw47d69WzNmzNCZM2dUW1urTZs2qaWlRT6fT9OmTdMbb7wht9vdU/UDAIAEE1P4ePXVV7tcS01NVXl5+RUXBAAA+je+2wUAAFhF+AAAAFZd8dUu8XbxtiPndFaK6Q4kfUeo9UKn8+fMWcuVoCd19TlLfNb9CZ9zcuBzvnxd9Sp0+j/zl3P7sJhvMtbTPvvsM+5yCgBAgmpoaNCIESMuuafPhY8LFy7o5MmTcrvdcrlckTueNjQ0fOMd05IB/YhGP6LRj2j0Ixr9iEY//icevTDGqLW1VX6/XwMGXPqsjj73Z5cBAwZ0mpjS0tKS/ofjq+hHNPoRjX5Eox/R6Ec0+vE/V9oLj8dzWfs44RQAAFhF+AAAAFb1+fDhOI5WrVrFl8/9F/2IRj+i0Y9o9CMa/YhGP/7Hdi/63AmnAACgf+vzRz4AAED/QvgAAABWET4AAIBVhA8AAGBVnw4fL730krKzs3XVVVdpwoQJ+tOf/tTbJVmzb98+zZw5U36/Xy6XS9u3b49aN8aopKREfr9fqampys/P15EjR3qn2B4WCAR0xx13yO12Kz09XbNnz9axY8ei9iRTP0pLSzVu3LjIzYDy8vL01ltvRdaTqRdfFwgE5HK5tGTJkshcsvWjpKRELpcrani93sh6svVDkv7+97/rpz/9qb797W9ryJAhuvXWW3Xo0KHIejL15Prrr+/w8+FyuVRcXCzJYi9MH1VWVmYGDRpk1q9fb44ePWoWL15shg4daj755JPeLs2KXbt2maeeesps2bLFSDLbtm2LWl+zZo1xu91my5Ytpra21tx///3G5/OZUCjUOwX3oLvvvtts2LDBfPjhh6ampsbce++95rrrrjOnT5+O7EmmfuzcudO8+eab5tixY+bYsWNm5cqVZtCgQebDDz80xiRXL77q/fffN9dff70ZN26cWbx4cWQ+2fqxatUqM3r0aNPY2BgZzc3NkfVk68c///lPk5WVZYqKisx7771n6uvrzdtvv20+/vjjyJ5k6klzc3PUz0ZFRYWRZN555x1jjL1e9Nnw8d3vftfMnz8/au7mm282Tz75ZC9V1Hu+Hj4uXLhgvF6vWbNmTWTu3//+t/F4POa3v/1tL1RoV3Nzs5FkqqqqjDH0wxhjrrnmGvPKK68kbS9aW1tNTk6OqaioMFOnTo2Ej2Tsx6pVq8z48eM7XUvGfixfvtxMmTKly/Vk7MlXLV682Nxwww3mwoULVnvRJ//s0t7erkOHDqmgoCBqvqCgQPv37++lqvqO+vp6NTU1RfXHcRxNnTo1KfoTDAYlScOGDZOU3P04f/68ysrK1NbWpry8vKTtRXFxse6991794Ac/iJpP1n7U1dXJ7/crOztbDzzwgI4fPy4pOfuxc+dOTZw4UT/+8Y+Vnp6u2267TevXr4+sJ2NPLmpvb9drr72mRx55RC6Xy2ov+mT4+OKLL3T+/HllZGREzWdkZKipqamXquo7LvYgGftjjNHSpUs1ZcoUjRkzRlJy9qO2tlZXX321HMfR/PnztW3bNo0aNSope1FWVqa//OUvCgQCHdaSsR+5ubnatGmTysvLtX79ejU1NWny5Mk6depUUvbj+PHjKi0tVU5OjsrLyzV//nz94he/0KZNmyQl58/IRdu3b1dLS4uKiook2e1Fn/tW269yuVxRj40xHeaSWTL2Z+HChfrggw/07rvvdlhLpn7cdNNNqqmpUUtLi7Zs2aJ58+apqqoqsp4svWhoaNDixYu1Z88eXXXVVV3uS5Z+SFJhYWHkf48dO1Z5eXm64YYbtHHjRk2aNElScvXjwoULmjhxolavXi1Juu2223TkyBGVlpbqZz/7WWRfMvXkoldffVWFhYXy+/1R8zZ60SePfAwfPlwDBw7skLSam5s7JLJkdPHM9WTrz6JFi7Rz50698847GjFiRGQ+GfsxePBg3XjjjZo4caICgYDGjx+vF198Mel6cejQITU3N2vChAlKSUlRSkqKqqqq9Otf/1opKSmRf3Oy9KMzQ4cO1dixY1VXV5d0Px+S5PP5NGrUqKi5W265RZ9++qmk5PzvhyR98sknevvtt/Xzn/88MmezF30yfAwePFgTJkxQRUVF1HxFRYUmT57cS1X1HdnZ2fJ6vVH9aW9vV1VVVb/sjzFGCxcu1NatW7V3715lZ2dHrSdbPzpjjFE4HE66Xtx1112qra1VTU1NZEycOFE/+clPVFNTo+985ztJ1Y/OhMNhffTRR/L5fEn38yFJ3/ve9zpcmv+3v/1NWVlZkpL3vx8bNmxQenq67r333sic1V7E9fTVOLp4qe2rr75qjh49apYsWWKGDh1qTpw40dulWdHa2moOHz5sDh8+bCSZtWvXmsOHD0cuNV6zZo3xeDxm69atpra21jz44IP99tKwxx9/3Hg8HlNZWRl1idi//vWvyJ5k6seKFSvMvn37TH19vfnggw/MypUrzYABA8yePXuMMcnVi8589WoXY5KvH7/85S9NZWWlOX78uDlw4ID54Q9/aNxud+S/ncnWj/fff9+kpKSYp59+2tTV1ZnNmzebIUOGmNdeey2yJ9l6cv78eXPdddeZ5cuXd1iz1Ys+Gz6MMeY3v/mNycrKMoMHDza333575NLKZPDOO+8YSR3GvHnzjDH/uTxs1apVxuv1GsdxzJ133mlqa2t7t+ge0lkfJJkNGzZE9iRTPx555JHI78W1115r7rrrrkjwMCa5etGZr4ePZOvHxfsyDBo0yPj9fjN37lxz5MiRyHqy9cMYY/74xz+aMWPGGMdxzM0332xefvnlqPVk60l5ebmRZI4dO9ZhzVYvXMYYE99jKQAAAF3rk+d8AACA/ovwAQAArCJ8AAAAqwgfAADAKsIHAACwivABAACsInwAAACrCB8AAMAqwgcAALCK8AEAAKwifAAAAKsIHwAAwKr/BzllsvFTGbYbAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "mx = (img == maximum_filter(img, size=3))\n",
    "\n",
    "plt.imshow(mx)\n",
    "\n",
    "index = np.where(mx)\n",
    "coord = list(zip(index[0], index[1]))\n",
    "print(coord)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "2b27c330-c802-4396-a05a-de9bbc200337",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[False, False, False, ..., False, False, False],\n",
       "       [False, False, False, ..., False, False, False],\n",
       "       [False, False, False, ..., False, False, False],\n",
       "       ...,\n",
       "       [False, False, False, ..., False, False, False],\n",
       "       [False, False, False, ..., False, False, False],\n",
       "       [False, False, False, ..., False, False,  True]])"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mx"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "0031553b-f326-479b-92fa-0cb00d39b186",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "18"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(coord)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "68e146b8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.26597595"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "img[(2,6)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "69d20df2-bd08-44a5-9735-957de9a2fb18",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[1., 1., 1.],\n",
       "       [1., 1., 1.],\n",
       "       [1., 1., 1.]])"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.ones((3,3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "f978f402",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.19295397"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "img[(36,7)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "736731b3-55e7-4057-b35f-b11ec604b517",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mx[(0,25)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3b09a28b-2b32-4266-a17d-2161761058db",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}