{ "cells": [ { "cell_type": "markdown", "id": "5e30f59a", "metadata": {}, "source": [ "# Demo of PG-RAD\n", "\n", "This is a demo." ] }, { "cell_type": "code", "execution_count": 1, "id": "415fdd25", "metadata": {}, "outputs": [], "source": [ "from matplotlib import pyplot as plt\n", "\n", "from pg_rad.dataloader import load_data\n", "from pg_rad.path import path_from_RT90\n", "from pg_rad.landscape import Landscape\n", "from pg_rad.objects import Source\n", "\n", "from pg_rad.isotopes import Isotope" ] }, { "cell_type": "markdown", "id": "8431d39b", "metadata": {}, "source": [ "## Loading a file" ] }, { "cell_type": "code", "execution_count": 2, "id": "5a0e470a", "metadata": {}, "outputs": [], "source": [ "FILENAME = \"B10_NaIR_MGS_ROI_CPS_IPL.CSV\"\n", "df = load_data(FILENAME)" ] }, { "cell_type": "markdown", "id": "baa7aba8", "metadata": {}, "source": [ "## Demo: Path regression" ] }, { "cell_type": "code", "execution_count": 3, "id": "2ec97553", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2026-01-27 20:41:53,431 - INFO: Piecewise regression reduced path from 105 to 4 segments.\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGwCAYAAABPSaTdAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAcMhJREFUeJzt3XdcVfUfx/HXZSoi4AYSV06cZWlUlqblzllp5sqsTC3NzGyoLTVbWpq2TBvqT00rd+YucZHknjkTpFRARFDg/P44efMKFwG53Au8n4/HeSDn+z3f8zmnm/fjOd9hMQzDQERERKSAcnN2ACIiIiKOpGRHRERECjQlOyIiIlKgKdkRERGRAk3JjoiIiBRoSnZERESkQFOyIyIiIgWah7MDcAVpaWmcOnWK4sWLY7FYnB2OiIiIZIFhGJw/f57g4GDc3Ow/v1GyA5w6dYqQkBBnhyEiIiI5cOLECcqXL2+3XMkOULx4ccC8WX5+fk6ORkRERLIiPj6ekJAQ6/e4PUp2wPrqys/PT8mOiIhIPnO9LijqoCwiIiIFmpIdERERKdCU7IiIiEiBpj47IiLiUKmpqVy+fNnZYUg+5Onpibu7+w23o2RHREQcwjAMoqOjiY2NdXYoko8FBAQQGBh4Q/PgKdlxhpQU8NCtF5GC7UqiU7ZsWXx8fDRpq2SLYRgkJiYSExMDQFBQUI7b0jeuI126BAcOwO7dsGeP+XP3bjh2DGJjwcvL2RGKiDhEamqqNdEpVaqUs8ORfKpo0aIAxMTEULZs2Ry/0lKykxvsJTUHD0JqasbHHDwItWvnbZwiInnkSh8dHx8fJ0ci+d2Vz9Dly5eV7DhNnz7w3Xfmq6ns2LNHyY6IFHh6dSU3Kjc+Q0p2blSxYtlPdMB88vPQQ7kfj4iIC/CuUoXq8fGZLs7odMWLw19/OTsKyQNKdm5UTp/O7N6du3GIiLiShATcL1xwdhQigCYVvHGhoTk7bs+e3I1DRERcQtOmTRkyZIj190qVKjFx4kSnxNKnTx86duzolHO7EiU7NyqnT3YOHDA7NouIiEvp06cPFosl3Xbo0CFnh2bX0aNHsVgsREZGOjsUl6Rk50aVKWNu2ZWSYo7IEhERl9OqVSuioqJstsqVKzs7LMkhJTu5Qa+yREQKFG9vbwIDA202d3f3DF8LDRkyhKZNm+b4XFfafP311ylTpgx+fn48/fTTXLrq6f/y5cu5++67CQgIoFSpUrRr147Dhw9by68kYrfccgsWiyVdPO+99x5BQUGUKlWKgQMHFrrlO5Ts5AZ1UhYRkRuwatUq9u7dy9q1a5k9ezYLFizg9ddft5ZfuHCB559/nm3btrFq1Src3Nzo1KkTaWlpAGzZsgWAX375haioKBYsWGA9ds2aNRw+fJg1a9Ywc+ZMZsyYwYwZM/L0+pxNo7Fyg5IdEZECZfHixfj6+lp/b926NfPmzXPY+by8vJg+fTo+Pj7Url2bN954g+HDh/Pmm2/i5uZGly5dbOpPnz6dMmXKsGfPHurUqUOZf7tTlCpVisDAQJu6JUqUYPLkybi7u1OzZk3atm3LqlWr6N+/v8Oux9Uo2ckNOX2NFRGRu3GIiEiuaNasGVOnTrX+XqxYMYeer379+jazTYeFhZGQkMCJEyeoWLEiBw8eZNSoUWzevJl//vnH+kTn+PHj1KlTJ9O2a9eubTPzcFBQEDt37nTMhbgop77Gmjp1KvXq1cPPzw8/Pz/CwsJYtmyZtbxp06bpesM//fTTNm0cP36ctm3b4uPjQ9myZRk+fDgpOZnk70bk9MnOkSPmOlkiIuJSihUrRtWqVa3blUUo3dzcMAzDpm5e9H9p3749Z8+e5fPPP2fz5s1s3rwZwKZfjz2enp42v1ssFmuyVFg49clO+fLlGT9+PNWqVcMwDGbOnEmHDh3Yvn07tf9NIPr3788bb7xhPebqzDc1NZW2bdsSGBjIxo0biYqKolevXnh6ejJ27Ni8u5ArI7L+/jv7x65ZYy45ISIiLq9MmTLs2rXLZl9kZGS6hCK7/vjjDy5evGhd+HLTpk34+voSEhLCmTNn2L9/P59//jlNmjQB4Ndff7U53uvfhaVT7a3HWMg59clO+/btadOmDdWqVaN69eq8/fbb+Pr6smnTJmsdHx8fm97wfn5+1rKff/6ZPXv28O2339KgQQNat27Nm2++yZQpU7KU7eaq22+3+TUVN9ZyL7PpxlruJdXerV6zJg+CExGR3HDfffexbds2vv76aw4ePMjo0aPTJT85cenSJfr168eePXtYunQpo0ePZtCgQbi5uVGiRAlKlSrFZ599xqFDh1i9ejXPP/+8zfFly5alaNGiLF++nNOnTxMXF3fDMRUkLjMaKzU1lTlz5nDhwgXCwsKs+7/77jtKly5NnTp1GDlyJImJiday8PBw6tatS7ly5az7WrZsSXx8PLsz6fybnJxMfHy8zXbDmjWz/nEBnajEUZqxlkeZTTPWUomjLKBT+uOU7IiI5BstW7bktdde48UXX+T222/n/Pnz9OrV64bbbd68OdWqVeOee+7hkUce4cEHH2TMmDGA+epszpw5REREUKdOHYYOHcq7775rc7yHhwcfffQRn376KcHBwXTo0OGGYypILMa1Lx/z2M6dOwkLCyMpKQlfX19mzZpFmzZtAPjss8+oWLEiwcHB7NixgxEjRtCoUSPrkLonn3ySY8eOsWLFCmt7iYmJFCtWjKVLl9K6desMzzlmzBibIX1XxMXF2Tw5ypbff4eGDVlAJ7oyH/Om/pdLWjDfj86nK51ZaHvs8eMQEpKz84qIuCAjOJi0fxcCddl1z11kIdA+ffoQGxvLDz/84OxQXFJSUhJHjhyhcuXKFClSxKYsPj4ef3//635/O300Vo0aNYiMjCQuLo758+fTu3dv1q1bR2hoKE8++aS1Xt26dQkKCqJ58+YcPnyYm2++OcfnHDlypM0jwPj4eEJuNNmoX59U/5I8FzcpXaIDYOCGhTSGMJEO/Ig7V3UOCw9XsiMiBUryn3/a/YISyWtOf43l5eVF1apVadiwIePGjaN+/fpMmjQpw7qNGzcGsK5PEhgYyOnTp23qXPn92nkGrubt7W0dAXZlu2Hu7myoM4CThGDvthq4cYIKbKCJbcHGjTd+fhEREcmQ05Oda6WlpZGcnJxh2ZUFzq4MAQwLC2Pnzp3ExMRY66xcuRI/Pz9Cczr3zQ2IqhR2/UpAFEG2O375BZz7NlFERJxkxowZeoXlYE59jTVy5Ehat25NhQoVOH/+PLNmzWLt2rWsWLGCw4cPW/vvlCpVih07djB06FDuuece6tWrB8ADDzxAaGgoPXv2ZMKECURHR/Pqq68ycOBAvL298/x6gsIqwXdZqEeU7Y7du2HlSnjgAYfEJSIiUpg59clOTEwMvXr1okaNGjRv3pytW7eyYsUK7r//fry8vPjll1944IEHqFmzJsOGDaNLly4sWrTIery7uzuLFy/G3d2dsLAwHnvsMXr16mUzL09eatK3KuU5ae2MnJ5BCMdpwob0RePGOTQ2ERGRwsqpT3a+/PJLu2UhISGsW7fuum1UrFiRpUuX5mZYOebu482kGp/Qdf9bWEjDSJdLWmjBL7adk69YuxY2bIAmTdKXiYiISI65XJ+d/K7zQ+7Mpys3YTuc0Z9zAMygT8bz7QA8/jgkJDg6RBERkUJFyU5u69aNzizkKJVYQ1Nm0Z01NOUfSvMMUzBwowffsYnG6Y89dAiGDMnzkEVERAoyJTu5rXZtaNECd9Joyjq6M4emrMODNCbxHO1YRBJFac8iDpHBXEFffgn/TpooIiIFx9q1a7FYLMTGxjo1jjFjxtCgQQOnnNtZ90DJjiOMHJnhbg9SmUM3GrKNfyhDG5byD6XSV+zfH64aTi8iUpilpprdGmfPNn86eq3LPn36YLFYsFgseHp6UrlyZV588UWSkpIce+ICpmnTpgxxkbcVSnYcoVkzaNQow6JiJLKYdlTkKAepTkd+IIlrhsmfPQuvvZYHgYqIuLYFC6BSJfOv1UcfNX9WquT4B+CtWrUiKiqKP//8kw8//JBPP/2U0aNHO/ak4jBKdhzBYrH7dAcgkNMspQ3+xPIbd9OLr0m7dvWYL7+EvXsdHKiIiOtasAC6doWTJ233//WXud+RCY+3tzeBgYGEhITQsWNHWrRowcqVK63laWlpjBs3jsqVK1O0aFHq16/P/PnzbdpYunQp1atXp2jRojRr1oyjR4/alGf0OmnixIlUqlTJZt/06dOpXbs23t7eBAUFMWjQIGtZbGwsTzzxBGXKlMHPz4/77ruPP/74w+b48ePHU65cOYoXL06/fv2u+4TqyqumJUuWUK9ePYoUKcIdd9xhs7r7mTNn6N69OzfddBM+Pj7UrVuX2bNnW8v79OnDunXrmDRpkvUp2dXXHxERwW233YaPjw933nkn+/fvzzSmG6Vkx1EefNDu0x2AUPbyAx3x5BLzeJiXGG9bITU104RJRCS/MQy4cCFrW3w8PPtsxpPLX9n33HNmvay0dyOT1O/atYuNGzfi5eVl3Tdu3Di+/vprpk2bxu7duxk6dCiPPfaYdcqUEydO0LlzZ9q3b09kZCRPPPEEL730UrbPPXXqVAYOHMiTTz7Jzp07+emnn6hataq1/KGHHiImJoZly5YRERHBrbfeSvPmzTl79iwAc+fOZcyYMYwdO5Zt27YRFBTEJ598kqVzDx8+nPfff5+tW7dSpkwZ2rdvz+XLlwFzcc6GDRuyZMkSdu3axZNPPknPnj3ZsmULAJMmTSIsLIz+/fsTFRVFVFSUzRqUr7zyCu+//z7btm3Dw8ODxx9/PNv3JlsMMeLi4gzAiIuLy92Gt283DE9PwzD/P8tw+5ZHrb9OYUD6Ohs25G5MIiJ54OLFi8aePXuMixcvWvclJGT616FDt4SErMfeu3dvw93d3ShWrJjh7e1tAIabm5sxf/58wzAMIykpyfDx8TE2btxoc1y/fv2M7t27G4ZhGCNHjjRCQ0NtykeMGGEAxrlz5wzDMIzRo0cb9evXt6nz4YcfGhUrVrT+HhwcbLzyyisZxrlhwwbDz8/PSEpKstl/8803G59++qlhGIYRFhZmPPPMMzbljRs3Tnfeq61Zs8YAjDlz5lj3nTlzxihatKjxv//9z+5xbdu2NYYNG2b9/d577zWee+65DNv+5ZdfrPuWLFliADaflatl9Fm6Iqvf33qy40gNGsDYsZlW6cEs3uIVAAbzMYtoZ1vhqafg4kUHBSgiIhlp1qwZkZGRbN68md69e9O3b1+6dOkCmItRJyYmcv/99+Pr62vdvv76aw4fPgzA3r17rYtXXxEWlrX1E6+IiYnh1KlTNG/ePMPyP/74g4SEBEqVKmUTx5EjR3IljqvrlSxZkho1arD33+4VqampvPnmm9StW5eSJUvi6+vLihUrOH78eJbavrLsE/y33mWMAwfmOHUG5ULh+edh2TJYvdpulZcZy1Eq8QX96cYc1nEvtxFhFu7ZAy++CB9/nEcBi4g4ho9P1udNXb8e2rS5fr2lS+Gee7J27uwoVqyY9XXR9OnTqV+/Pl9++SX9+vUj4d+LWLJkCTfddJPNcdlZl9HNzQ3jmvdrV14TARQtWjTT4xMSEggKCmLt2rXpygICArIcR068++67TJo0iYkTJ1K3bl2KFSvGkCFDuHTpUpaO9/T0tP7ZYjH7rKal2Vtq6cbpyY6jubnBzJlQooTdKhbgE56hJctJpBjtWMxRKv5XYfJk8/9oEZF8zGKBYsWytj3wAJQvbx5jr62QELNeVtqz105WuLm58fLLL/Pqq69y8eJFQkND8fb25vjx41StWtVmu9IvpVatWtb+K1ds2rTJ5vcyZcoQHR1tk/BERkZa/1y8eHEqVarEqlWrMozr1ltvJTo6Gg8Pj3RxlC5d2hrH5s2bM43DnqvrnTt3jgMHDlCrVi0AfvvtNzp06MBjjz1G/fr1qVKlCgcOHLA53svLi1RHzxOQRUp28kL58td9neVJCnN5mPpEcppA2rCUcwT8V6FvXzh92rFxioi4CHd3mDTJ/PO1icqV3ydONOvlhYceegh3d3emTJlC8eLFeeGFFxg6dCgzZ87k8OHD/P7773z88cfMnDkTgKeffpqDBw8yfPhw9u/fz6xZs5gxY4ZNm02bNuXvv/9mwoQJHD58mClTprBs2TKbOmPGjOH999/no48+4uDBg9bzALRo0YKwsDA6duzIzz//zNGjR9m4cSOvvPIK27ZtA+C5555j+vTpfPXVVxw4cIDRo0eze/fuLF3zG2+8wapVq9i1axd9+vShdOnSdOzYEYBq1aqxcuVKNm7cyN69e3nqqac4fc13VKVKldi8eTNHjx7ln3/+ceiTm+vKtEdPIeGwDspXu3TJMKpXv24vupMEG+U5boBh3MsaIwmv/8pbtDCMlBTHxSgikksy61SaHd9/bxjly9v+VRkSYu53lN69exsdOnRIt3/cuHFGmTJljISEBCMtLc2YOHGiUaNGDcPT09MoU6aM0bJlS2PdunXW+osWLTKqVq1qeHt7G02aNDGmT59u00HZMAxj6tSpRkhIiFGsWDGjV69exttvv23TQdkwDGPatGnW8wQFBRmDBw+2lsXHxxuDBw82goODDU9PTyMkJMTo0aOHcfz4cWudt99+2yhdurTh6+tr9O7d23jxxRez1EF50aJFRu3atQ0vLy+jUaNGxh9//GGtc+bMGaNDhw6Gr6+vUbZsWePVV181evXqZXPf9u/fb9xxxx1G0aJFDcA4cuSIte2r78H27dut5RnJjQ7KFsO4kQF5BUN8fDz+/v7ExcXh5+fnuBMtWAD/dnDLzA7qcje/ch4/evAt39Dzv1l4xowBTWwlIi4uKSmJI0eOULlyZYoUKXJDbaWmwoYNEBUFQUHQpEnePdEpjNauXUuzZs04d+6cw/v+ZEVmn6Wsfn/rNVZe6tQJstALvh47+Z4ueHCZ73iM13jzv8LXX4dffnFgkCIirsXdHZo2he7dzZ9KdCS7lOzkJYvFHFWVhf9T7+cXPuNJAN7mVb6gn1lgGOac6ddOKSoiIiIZUrKT1xo2hFGjslS1LzN4jTcAeJpprOABs+Dvv6FdOzh/3lFRiohIIdW0aVMMw3CJV1i5RcmOM7z8Mtx5Z5aqvs5oevI1qXjQlflEUt8s+OMP6NYNUlIcGKiIiEj+p2THGTw84NtvoXjx61a1AF/wBPexigSK05YlnKC8Wbh0qbk4jPqYi4iL0hgYuVG58RlSsuMslSubkwVmgReX+Z4u1GYXp7iJtiwhjn97nX/yiWZXFhGXc2WG3MTERCdHIvndlc/Q1bMuZ5eGnpOHQ8+vZRjQpw98/XWWqh+jAnewiWiCuJ+fWUJbPEkxOzyvWAF21k8REXGGqKgoYmNjKVu2LD4+PtZlAUSywjAMEhMTiYmJISAgwLqG1tWy+v2tZAcnJjsAFy5Ao0bmGlhZ8Du3cA/ruYAvfZnOl/Qz5+ApWRK2bTOfGImIuADDMIiOjiY2NtbZoUg+FhAQQGBgYIbJspKdbHBqsgOwdy/cdhtk8XHvUlrTnkWk4c7rjGLUlXl46tWD334DX18HBisikj2pqak2C1yKZJWnpyfumUzXomQnG5ye7IDZYblnzyxX/5QneZpPAZhJL3rxjVnQuDH88AMEBjogSBEREdehGZTzm8ceg0GDslz9KT5jBOMB6MeXrKaZWbB5M9x+O/z+uyOiFBERyXeU7LiSDz+ENm2yXH0sL9ON2aTgSWcWsJtQs+DkSbj7bpg/30GBioiI5B9KdlyJhwfMmQP162epuhsGX9GXu9lAHAG0Zhmn+Le3+sWL8NBD8Mor5ip6IiIihZSSHVdTvDgsXgzBwVmqXoRkfqQDNdjHCSrQjsUkUOy/CmPHmk+LzpxxUMAiIiKuTcmOKypfHhYtgqJFs1S9JOdYShvKEMN2buUR/kcKV/Ve//lnc00u9eMREZFCSMmOq7r1VvjyyyxXr8IRFtOOoiSylLYMYjI2w+yOHTPX41q4MNdDFRERcWVKdlxZ9+4wfHiWqzdiK7PpjoU0PuVpJvCibYXkZHj4YViyJJcDFRERcV1KdlzduHHQsmWWq3fgJyYyBICXeIc5PGJbISUFHn0UDh7MxSBFRERcl5IdV+fuDrNnQ4MGWT7kWT5mKB8A0JuZrKeJbYX4eOjQARIScjFQERER16RkJz8oUQJ++QXuuSfLh7zHC3Tmey7hTUd+YB81bCvs3QshIRqlJSIiBZ6SnfyiVClYuRKeeCJL1d0w+JbHuINwzlGSNizlNGVtK8XGQtmyMHEiXLqU6yGLiIi4AiU7+YmXF3z2GUyaBG7X/09XlCR+4kFu5hBHqMKD/EQi1wxnT0uDoUOhbl1YscJBgYuIiDiPU5OdqVOnUq9ePfz8/PDz8yMsLIxly5ZZy5OSkhg4cCClSpXC19eXLl26cPr0aZs2jh8/Ttu2bfHx8aFs2bIMHz6clJSUvL6UvGOxwLPPwtKl4O9/3epl+IeltKEkZ9hCYx5lFqkZ/Wc/cABatTJHf2ltWBERKUCcmuyUL1+e8ePHExERwbZt27jvvvvo0KEDu3fvBmDo0KEsWrSIefPmsW7dOk6dOkXnzp2tx6emptK2bVsuXbrExo0bmTlzJjNmzGDUqFHOuqS807IlbNsGdepct2p1DvITD+JNEj/SkaF8iN105r334PXXczVUERERpzJcTIkSJYwvvvjCiI2NNTw9PY158+ZZy/bu3WsARnh4uGEYhrF06VLDzc3NiI6OttaZOnWq4efnZyQnJ9s9R1JSkhEXF2fdTpw4YQBGXFyc4y7MURISDKN7d8Mwn8dkus2lq/XXDxhiv67FYhirVjn7ykRERDIVFxeXpe9vl+mzk5qaypw5c7hw4QJhYWFERERw+fJlWrRoYa1Ts2ZNKlSoQHh4OADh4eHUrVuXcuXKWeu0bNmS+Ph469OhjIwbNw5/f3/rFhIS4rgLc7RixeC778xOxtfxEPN5lxcAGMb7fE/njCsaBvTqBWfP5mKgIiIizuH0ZGfnzp34+vri7e3N008/zcKFCwkNDSU6OhovLy8CAgJs6pcrV47o6GgAoqOjbRKdK+VXyuwZOXIkcXFx1u3EiRO5e1F5zWKB556DLLy+G8b7PMMUDNx4jG8J546MK/71F3TqpP47IiKS7zk92alRowaRkZFs3ryZAQMG0Lt3b/bs2ePQc3p7e1s7RV/ZCoTRo80VzjNhASbxHO1YRBJFeZCfOMTNGVdevx7uv9+chFBERCSfcnqy4+XlRdWqVWnYsCHjxo2jfv36TJo0icDAQC5dukRsbKxN/dOnTxMYGAhAYGBgutFZV36/UqdQcXODb7+FKlUyreZBKnPoRkO28Q9laM0y/qFUxpVXrTIXJY2IcEDAIiIijuf0ZOdaaWlpJCcn07BhQzw9PVm1apW1bP/+/Rw/fpywsDAAwsLC2LlzJzExMdY6K1euxM/Pj9DQ0DyP3SWUKAGLFpmTEGaiGIksph0VOcohqtGBH7lIkYwrHz4MYWEwZYoDAhYREXEspyY7I0eOZP369Rw9epSdO3cycuRI1q5dS48ePfD396dfv348//zzrFmzhoiICPr27UtYWBh33GH2M3nggQcIDQ2lZ8+e/PHHH6xYsYJXX32VgQMH4u3t7cxLc67QUPMVVP36mVYL5DTLaE0A59jIXfRmJmlYMq58+TIMGmT2C1I/HhERyUecmuzExMTQq1cvatSoQfPmzdm6dSsrVqzg/vvvB+DDDz+kXbt2dOnShXvuuYfAwEAWLFhgPd7d3Z3Fixfj7u5OWFgYjz32GL169eKNN95w1iW5jtBQ2LoVPvrIfNpjRy328QMd8eQS83iYEbyTebtvvgnTpuVysCIiIo5jMQz9Mz0+Ph5/f3/i4uIKTmflq509CwMGwNy5dqt8x6M8xncATGYgA/nEfnseHmZfnmwsTCoiIpLbsvr97XJ9dsQBSpaE//0PXn7ZbpUezOItXgHgWT5iEe3st5eSAl27wvHjuR2piIhIrlOyU5i89RZ07Gi3+GXG8gSfk4Y73ZjDNhrab+vvv815eBITcz9OERGRXKRkpzCxWODzzyEoKONi4BOeoSXLSaQY7VjMUSrab+/33+Hxx9VhWUREXJqSncKmdGlzxXQ7nZY9SWEeD1GfSE4TSBuWco4A++39739aOFRERFyakp3CqEEDWLkSbropw+LiJLCEtpTnBHsJpRMLScbLfnuvvw5z5jgmVhERkRukZKewatgQIiOhbdsMi2/iFEtoS3HiWUdTHme6/Tl4APr2hS1bHBOriIjIDVCyU5iVLg0//QTvvWcOJ79GPXbyPV3w4DKz6MFrvGm/raQk6NwZzpxxYMAiIiLZp2SnsHNzg2HDYPVqyGDW6fv5hc/pD8BYXuFznrDf1l9/Qb9+6rAsIiIuRcmOmJo0gU8/zbCoDzMZhdkJeQBTWU5L++38+CNMneqICEVERHJEyY78p3dveO65DIvGMIZezCQVDx5iHpFksu7W889rlXQREXEZSnbE1nvvwX33pdttAT6nP/exigSK05YlnKB8xm0kJ0PTpuYQdxERESdTsiO2PDxg3jyoVi1dkReX+Z4u1GYXp7iJNiwlDjtrkSQkQPv2MGWKgwMWERHJnJIdSa9kSVi8GAIC0hUFEMdS2hDEKXZRl67M5zLpR3IBkJYGgwbBkCGQmurQkEVEROxRsiMZq14d5s8Hd/d0RRU4wRLaUowEfuF+nuQzMh1/NWkS9OwJly87LFwRERF7lOyIfc2bw0cfZVh0C5HM5WHcSWEGfXmT1zJva/Zsc+HQixcdEKiIiIh9SnYkcwMGQNeuGRa1YRlTGAjAaN5gJr0yb2vJEmjVCuLjcztKERERu5TsSOYsFvjsM6hQIcPip/iMEYwH4Am+YBXpR3LZWL8ewsJg69bcjlRERCRDSnbk+kqUgFmzzNmWMzCWl+nGbFLwpDML2EXtzNvbswfuuANeeslcZkJERMSBlOxI1tx1F7z9doZFbhjMoA9NWE88/rRhKacIyry9tDR45x245RYID3dAwCIiIiYlO5J1I0bA0KEZFnlziR/oSA32cYIKtGMxCRS7fpv79sHdd8MHH2hNLRERcQglO5J1Fgu8/745y7LFkq64JOdYShvKEMN2buVh5pJC+qHr6aSlmYuRjhyphEdERHKdkh3JHovFTEy+/x6KFk1XXIUjLKYdRUlkGW0YyJTM5+C52jvvmGtzpaXlasgiIlK4KdmRnOnUyRxZFRiYrqgRW5lNdyyk8RlP8Q4jst7uxx+bbavjsoiI5BIlO5Jzt90GmzZluI5WB35iEuYK6iMZz2y6Zb3dn34yh7rPnKllJkRE5IYp2ZEbU7EibNgA9eqlKxrMZIbyAQB9mMF6mmS93b//hj59oG5dc9mKlJRcClhERAobJTty48qVg7VrzckCr/EeL9CZ77mENx35gX3UyF7be/fCQw9B+fJmf57Nm9WJWUREskXJjuSOEiVg5Upo2dJmtxsG3/IYdxDOOUrSmmWcpmz22z992lyn6447oGpVeO012L07l4IXEZGCzGIY+mdyfHw8/v7+xMXF4efn5+xw8rfUVHjzTXMCwqtePf1NacII5zBVuZ0t/EJzfqchUQQRRBRN2IA7ORiFVbs2PPywudWsmYsXIiIiri6r399KdlCy4xDbt0PfvvDHH9ZdB6lKGOGcoTRFuEgS/w1dL88JJvEcnVmY83PWrQv33mu+Vitb9r+fgYEQEgIeHjdyRSIi4mKU7GSDkh0HuXQJxo83n/T8+5RnLC/xCmMB20kJLf8+1ZlP1xtLeOzx9jZHj911l7ndeSeULp375xERkTyjZCcblOw42LZt0Lo1qf+cpRJHOUl5rk12wEx4ynOSI1TO2Sut7KpRw0x6WrSA+++HMmUcf04REck1SnayQclOHti9m7VNXqPZuQXXrXorETRiC9U4aN2q8CdeXHZcfBaL+eSnVStza9RIr71ERFyckp1sULKTN2Z/GMWjz19nNXQ73EilIsdsEqArWyWO4kkuz8Pj5gZ+fuDvb7v5+ppbsWK2m4+P+aqsSBHz55U/Fy1qTpAYGJjhemIiIpJzWf3+1j9dJc8E3ZK1RGcE4/EgxSalSaA4R6jCEarwM7bD291JoTJHqMqhdIlQRY7hQQ5mYU5Lg9hYc8sNZcpAs2bQvLm5Vami5EdEJI/oyQ56spNXUlOhUiX46y8Dw8h6nx0DOE25DJ7pVOMQVUmkmN1zenKJyhzJ8IlQCCfypm9QRipUMJOfkBAoVQpKlvzvZ0CA+UTo2idF7llYQV5EpBDRa6xsULKTdxYsgK5dzT9f/cnL6WgsAzhFcIaJ0GFuthnefi0vkrmZw1TjYLqnQuU5iVvW12vPGx4eGb8q8/Y2X7lVqmQu21G/vrmpw7WIFHD5ItkZN24cCxYsYN++fRQtWpQ777yTd955hxo1/ltSoGnTpqxbt87muKeeeopp06ZZfz9+/DgDBgxgzZo1+Pr60rt3b8aNG4dHFjuYKtnJWwsWmCs/nDz5376Qon8z8eJTuTrsPA0Lf3GT3UToEt52jy3CRWsidO0WzKkMxpK5oKAgM+mpW9cceVa9uvmzTBm9QhORAiFfJDutWrWiW7du3H777aSkpPDyyy+za9cu9uzZQ7Fi5quJpk2bUr16dd544w3rcT4+PtaLSk1NpUGDBgQGBvLuu+8SFRVFr1696N+/P2PHjs1SHEp28l5qqrl+aFSU+Z3cpAm4R2yBCRPMVc8vO3DkFZCKGycIyTAR+pMqpOBp91gfLtg8Cbr6z4FEu34i5O9vJj41a8Itt0DDhubP4sWdHZmISLbki2TnWn///Tdly5Zl3bp13HPPPYCZ7DRo0ICJEydmeMyyZcto164dp06doly5cgBMmzaNESNG8Pfff+Pl5ZXumOTkZJKTk62/x8fHExISomTHVZw5A99/D7NmwTVP9fJCCu4co6K1T9DVidARKpOaSb9+X85n2FG6Ggcpw9+umwhZLGYC1LChOQS/UyfztZiIiAvLl8nOoUOHqFatGjt37qROnTqAmezs3r0bwzAIDAykffv2vPbaa/j4+AAwatQofvrpJyIjI63tHDlyhCpVqvD7779zyy23pDvPmDFjeP3119PtV7Ljgk6cgDlzzMTnqv/GznIZD45SKcMnQseoSBr2OxH7EZfuqCuJUSnOuF4i1LQpjB5t/hQRcUH5LtlJS0vjwQcfJDY2ll9//dW6/7PPPqNixYoEBwezY8cORowYQaNGjViwwJyc7sknn+TYsWOsWLHCekxiYiLFihVj6dKltG7dOt259GQnnzpwAObNg7lzYccOZ0eTTjJeHKFyhonQCUIwcLN7bADnMnwaVI2DlCA27y4iIw8+aL5evKovnYiIK8h3yc6AAQNYtmwZv/76K+XLl7dbb/Xq1TRv3pxDhw5x88035yjZuZb67ORD+/aZic/335uJj2t8jO1Kwps/qZJhOnOSkEyPLcU/GSZBVTmEP/G5El8qbmygif1V6D084OmnzSc9WlNMRFxEvppUcNCgQSxevJj169dnmugANG7cGMCa7AQGBrJlyxabOqdPnwYgMDDQMQGL89WsCa+9Zm7nz5tPfWJizO306f9+7t1rvv5KzcHEgrmoCMmEspdQ9qYrS6Qoh7k5w0QoimDOUJozlGYTYemOLUOM3SdCvlzIUmwL6MRzTLJJutKtQp+SApMnw9dfw1tvwcCB5izTIiL5gFOTHcMwGDx4MAsXLmTt2rVUrlz5usdc6ZsTFGTOxhsWFsbbb79NTEwMZcuWBWDlypX4+fkRGhrqsNjFhRQvbnastSchAbZsgd9+M7eNG80EyUX4cJG67KIuu9KVJVDM2kn62s7Spwnkb8ryN2XZyF3pjg0kyu4TIR8uAmai05X56WYU+oub6Mr89PMexcfDs8+arxKnT4dq1XLzVoiIOIRTX2M988wzzJo1ix9//NFmbh1/f3+KFi3K4cOHmTVrFm3atKFUqVLs2LGDoUOHUr58eevcO1eGngcHBzNhwgSio6Pp2bMnTzzxhIaeS8YuX4bNm2H5cnOLiHB2RDkST/F0CdCV7R8yn1AwmL+oykEiuI0LFCNHq9AXKQJvv21OmqTZnUXECfJFnx2LnYnNvvrqK/r06cOJEyd47LHH2LVrFxcuXCAkJIROnTrx6quv2lzUsWPHGDBgAGvXrqVYsWL07t2b8ePHa1JByZqYGPj5Z1i7FnbvNl9/xcWZm5Nff+VULP52Xm5V4xwls9XWGprSlEymAAgLg6++UgdmEclz+SLZcRVKdiRDhgGJif8lPufPw4UL6beEBEhKguRk25+xseZrs7//dvaV2DhLCQ5Sja/pxScMvG79WXSnO3Myr+TjA5MmQb9+mp1ZRPJMvuqgLOKSLBYoVszcgoNz1oZhwK5dsGoVrF5tPj1ycn+hkpyjMVu4SNEsJTveJF2/0cRE6N/ffC342WfmgqYiIi5CT3bQkx3JQykp8OefcPQonD1rzhZ99c+zZ83E4dqnRFf/TE6GS5duOJRU3KjEUf7iJjtzABmABX/O8Q4v0Z/Ps7Y4avny8O23cO+9NxyjiEhm9BorG5TsSL6TlmYmPMnJZqK0cyf88Ye5RUaaCVUWXBmNBdgkPBbSMIDKHOUIVQBozCam8TQN+CNrMfbrB2PHwr+jJEVEcpuSnWxQsiMFTny8OfHigQPmtn//f39OTLSpmtE8OyEcZyJDeJCf+IRneJW3OI8fbqTyHJN4ndEUJ+H6cfj7w5tvwoAB5sSEIiK5SMlONijZkUIjLQ2OHTOH21/ZVq687gzKpwhiKB8yl0cAuImT/046uCBra3rVrQtTppjL24uI5BIlO9mgZEcKNcOAb74x58w5cCDTqit4gGf4hD+5GYA2LGEyg6jM0eufx2KBESPMGZg1L4+I5IKsfn9rvneRws5igV69zFFjH30EpUrZrdqSn9lFHV7jDTy5xFLaEsoexjKSS3hmfh7DgPHjoVMnp49IE5HCRcmOiJg8PWHwYDh0CJ56ym61oiTxBqPZSV3uYxVJFOUVxtKASNaShRFYixbBnXeaI9JERPKAkh0RsRUQANOmwS+/QKVKdqvV4AC/0IJv6UFZTrOXUJqxlt7MIOY6y1Wwaxfcfjv8+muuhi4ikhElOyKSsebNzSHtA+1PPGgBejCLfdTkaaZiIY2v6U1N9vEZ/UnLrPvyP/+Y55hzndmZRURukJIdEbHP1xcmTzZngL7pJrvVShDLVJ4hnDAasJ1zlOQpPuNufuUP6tlv/9Il6N4dJkww+/SIiDiAkh0Rub777jMnLOzYMdNqjdnCVm7nQ4bgy3nCuZOGRDCM90igmP0DR4yAQYPy7cKrIuLalOyISNaUKgULFsDUqVCkiN1qHqQyhEnsoyZdmUcqHnzAMGqxl4V0tL/gxCefQOfO6SY9FBG5UUp2RCTrLBZ4+mlzMsJGjTKtehOnmMfDLKU1lfmTk4TQmYW0ZxFHqJTxQT/9BC1amGuEiYjkEiU7IpJ9oaEQHg5ffAGlS2datTXL2UUdXuEtPLnEEtpRm92MZ0TGc/OEh8Pdd8OJEw4KXkQKGyU7IpIzbm7mYp8HDpjz87jZ/+vEh4u8xWv8QX2asoaL+DCS8dzCdtaTwRISe/eac/Hs2ePACxCRwkLJjojcmBIlzJmXt2+H2rUzrVqLfazmPr6mJ2WIYQ+1uZf19GU6f3PNE6KTJ81XZZMmqeOyiNwQJTsikjvq1YONG6Ft20yrWYCefMs+avIU0wCYQV9qso8v6Gc7N8+FCzBkCNx1lzkRoYhIDijZEZHc4+cHP/4IL7xw3aolOcc0BhDOHdQnkrOUoj9f0IQN7KSObeXNm+HWW2HYMDhyxEHBi0hBpWRHRHKXuzu8+y5Mnw4eHtetfgeb2cZtfMBQipHARu7iFrYznAm2c/NcvgwffAA33wzt28Py5ZCW5sALEZGCQsmOiDhG376wdCkUL37dqh6kMpSJ7KUWnfmeVDx4j+GEsocfedC2smHA4sXQujXUqGHO8Jyc7KCLEJGCQMmOiDjO/ffDhg0QHJyl6iGc5Hu6spi2VOIIJ6hAR37kQX7kGBXSH3DokDkSrE4dWL8+l4MXkYJCyY6IOFb9+rBpk5mQZFFblrKb2oxkLJ5cYhEPEsoeJjCcy2TwauzQIbj3XjPxSUjIxeBFpCBQsiMijhcSYj7hueeeLB/iw0XG8gqRNOAe1pFIMUYwgVvYzgbuzvigyZPNUWGrV+dS4CJSECjZEZG8ERAAK1ZAp07ZOiyUvaylKTPoTWn+Zjd1uIcNPM6X/EOp9AccOQLNm8PAgVpnS0QAJTsikpeKFIF58+DJJ7N1mAXozdfsoyb9+QyAr3icGuxnOn1t5+a54pNP4JZbYOvWXAhcRPIzJTsikrfc3WHaNHj99WwfWoqzfMZT/Mad1GUHZylFP6ZzL+vYRQazNx84YC478dZbkJKSC8GLSH6kZEdE8p7FAqNGmUPTy5fP9uF3Ek4EDXmPYRQjgV9pwi1sZwTjuYCPbeWUFHjtNbO/0J9/5tIFiEh+omRHRJyndWvYvRueeSbbh3qSwjA+YA+hdGQhKXgygRGEsoefaJ/+gPBwuO02WLkyFwIXkfxEyY6IOJefH0yZYo7WqlEj24dX4AQL6cxPtKciRzlORTrwEx1ZyHFCbCufOwetWsF775mTE4pIoaBkR0Rcw913Q2QkvPFGlmZdvlZ7FrOb2rzEODy4zI90pBZ7eZcXbOfmSUuD4cPhjjtg//7ci19EXJaSHRFxHUWKmP1rTpyASZOgevVsHV6MRMbxMpE0oAnrSaQYL/Iut/I7v3GnbeUtW6BmTWjaFL78EmJjc+0yRMS1WAxDz3Lj4+Px9/cnLi4OPz8/Z4cjIlekpcGqVeZrrkWLsrXwpwHMpDcv8B5nKA1AP77gHUZQirPpD/D2NvsQtW0LbdpkeYkLEXGerH5/K9lByY5IvrBnDzz+OGzenK3DzlCSEbzDlzwBQCn+4T1eoDczM5qd5z8NGvyX+DRubA6ZFxGXktXvb73GEpH8ITQUfvvN7FxcpEiWDyvFWb6gP79yF3XYyRlK05cZ3Ms6dhNq/8DISHj7bbjrLrj5Znj/fbh48cavQ0TynFOTnXHjxnH77bdTvHhxypYtS8eOHdl/TYfBpKQkBg4cSKlSpfD19aVLly6cPn3aps7x48dp27YtPj4+lC1bluHDh5OiCcRECh53dxg2DP74w0xCsuEuNvI7tzKB4fhwgQ3cQwMiGclYEima+cHHjsELL0Dt2rB48Q1cgIg4g1OTnXXr1jFw4EA2bdrEypUruXz5Mg888AAXLlyw1hk6dCiLFi1i3rx5rFu3jlOnTtG5c2dreWpqKm3btuXSpUts3LiRmTNnMmPGDEaNGuWMSxKRvFC9OqxbB+PGgadnlg/zJIXhvMceQunAD6TgyXhGEsoeFtP2+g0cOQLt20PHjuafRSR/MFxITEyMARjr1q0zDMMwYmNjDU9PT2PevHnWOnv37jUAIzw83DAMw1i6dKnh5uZmREdHW+tMnTrV8PPzM5KTk7N03ri4OAMw4uLicvFqRCRPREQYRq1ahmHOnJOt7UfaGxU4at3Vie+N45TPehvNmxvGxx8bxp49hpGW5uw7IVLoZPX726X67MTFxQFQsmRJACIiIrh8+TItWrSw1qlZsyYVKlQgPDwcgPDwcOrWrUu5cuWsdVq2bEl8fDy7d+/O8DzJycnEx8fbbCKST916K0REwODB2T70QRaxh1Be5B08uMxCOlOLvbzP87Zz89izapV53tBQc9mLnj3h++8hOTkHFyIijuIyyU5aWhpDhgzhrrvuok6dOgBER0fj5eVFQECATd1y5coRHR1trXN1onOl/EpZRsaNG4e/v791CwkJybCeiOQTRYvCRx/B3Lng43P9+lcpRiLv8BK/cyt38SsX8OUF3uc2thHOHVlv6NQp+PZb6NoVKleGiRMhMTF71yEiDuEyyc7AgQPZtWsXc+bMcfi5Ro4cSVxcnHU7ceKEw88pInngoYdg40aoVCnbh9ZlF+u5hy95nJKcYQf1uZNwnuRTzlIie41FRcHQoWbS8+67kJCQ7XhEJPe4RLIzaNAgFi9ezJo1ayh/1QrIgYGBXLp0idhrZjY9ffo0gYGB1jrXjs668vuVOtfy9vbGz8/PZhORAqJ+fdi2De67L9uHumHwOF+xnxo8zpcAfM6T1GA/X9OTbE9KFhMDL75oJl/jxsFVgy9EJO9kKdm5tn9LVrasMAyDQYMGsXDhQlavXk3lypVtyhs2bIinpyerVq2y7tu/fz/Hjx8nLCwMgLCwMHbu3ElMTIy1zsqVK/Hz8yM0NJM5NESk4CpVClasMJ+q5OAfM6U5w5c8wXqaUJtd/EMZevM1zVjDXmpmP54zZ+Dll6FqVfj0U7h8OfttiEiOZWkGZTc3NyyWTOcatW3UYuHAgQNUqVIl03rPPPMMs2bN4scff6TGVasd+/v7U7SoOe/FgAEDWLp0KTNmzMDPz4/B/3ZC3LhxI2AOPW/QoAHBwcFMmDCB6OhoevbsyRNPPMHYsWOzFK9mUBYpwP7+Gz75xOxPc+hQtg+/hCcfMpTXGc1FfPDkEi/wHq/yFj7kcJLB6tVh7Fjo3Bmy8XeriNjK1eUi3Nzc+P77762jpDJjGAZt2rRh165d10127CVQX331FX369AHMSQWHDRvG7NmzSU5OpmXLlnzyySc2r6iOHTvGgAEDWLt2LcWKFaN3796MHz8eD48sjKZAyY5IoWAY5uKf334Lc+bAP/9k6/CjVGQwH7OY9gBU4ghTGEgbluU8psaNzY7Md2SjI7SIWOVqslO5cmW2bdtGqVKlsnTyOnXqsGzZsnwzyknJjkghc/kyLF8OCxbA0qVm35osMIAf6cCzfMQJKgDQhflMZAjl+StnsVgs8Pzz8NZb2VoGQ0S0EGi2KNkRKcTS0uD332HJEjPx2brVfAqUiQSK8Tqj+ZChpOKBL+d5g1EM5mM8SM1ZHLVqwddfw2235ex4kUJIC4GKiGSFm5uZYIweba6ofuCAudJ5Jny5wLu8yO/cShgbSaA4z/Mht7GNTTTOWRx795qvs8aMUQdmkVyWoyc7W7duZc2aNcTExJCWlmZT9sEHH+RacHlFT3ZExIZhwA8/wHPPwXXm4UrDwnQe50UmcI6SWEjjST5jHCMpQWzOzl+7Nnz8MTRrlrPjRQoJh73GGjt2LK+++io1atSgXLlyNp2MLRYLq1evznnUTqJkR0QydOECjB8PU6bAuXOZVv2b0rzIBGbQF4AyxPABz9OD78jxeKuHHoL33oMKFXLagkiB5rBkp1y5crzzzjvW0VIFgZIdEclUYqLZmfnnn831sE6dslt1HfcwgKnsxZznqxmr+YRnqMn+nJ27aFF46SUYPtz8s4hYOazPjpubG3fdddcNBScikq/4+MBjj5kdiE+ehH37zKc9nTqBp6dN1XtZTyQNGMdLFCWRNdxHPXbwKm9ykRyMtrp40exPFBoKP/543c7TIpJetpOdoUOHMmXKFEfEIiLi+iwWqFEDnnnGfNpz6JD5Zy8vaxUvLvMS77Cb2rRhCZfx4m1epQ67WE7LnJ336FHo2BFat4b9OXxKJFJIZfs1VlpaGm3btuXAgQOEhobiec2/ahYsWJCrAeYFvcYSkRv2118wYQJ89hkkJVl3G8BCOvEckziJOfdYV+YxkSHchP3XYZny9DQXGn3tNfD1zYXgRfInh73GevbZZ1mzZg3Vq1enVKlS+Pv722wiIoXSTTfBpEnmk57+/c0h7YAF6MxC9hDK87yPOynM5yFqso9JPEsK7tk/1+XLZmLVoAFERubmVYgUSNl+slO8eHHmzJlD27ZtHRVTntOTHRHJdXv3mot//vCDze4/qMfTTGMT5mLGDdjOpzxFI7bm7DxFi8JXX8Ejj9xgwCL5j8Oe7JQsWZKbb775hoITESnwatWChQvht9/g1lutu+uzg9+4i095khKcJZJbuINNPMMUYsnB0/GLF6FbN3PEVmoOZ28WKeCyneyMGTOG0aNHk5iY6Ih4REQKljvvhE2b4PXX4d/Fid0weJLP2UdNejETAzem8gw12M93PEqOxlu98w60b3/d+YBECqNsv8a65ZZbOHz4MIZhUKlSpXQdlH///fdcDTAv6DWWiOSJ33+HXr1g926b3WtoyjN8wj5qAXAfq/iEZ6jBgeyfo2JFc1V3raQuhUBWv789sttwx44dbyQuEZHC69ZbISICRo2C99+3vnZqxlr+oD7v8QJv8hqraU49djCCdxjJOIqSdJ2Gr3LsGDRpAm+/DS+8YO0oLVKYadVz9GRHRJxg50549llYu9Zm959UZhCTWYa5GOnNHGIKA2nJz9k/R6tWMHMmlC2bCwGLuB6tei4i4srq1oXVq2HuXAgJse6uwhGW0Jb5dCGYvzhMVVqxgkeYwymCsneO5cvNRUVnzdLMy1KoZSnZKVmyJP/880+WG61QoQLHjh3LcVAiIoWCxWIu9rl3L7z6Knh7m7uBLixgHzUZwoe4kcpcHqEm+/iIwaRm59+p//wDPXpAu3Zw/LhjrkPExWXpNZabmxszZ87M8qSB3bt3Z+fOnVSpUuWGA8wLeo0lIi7hzz/NmZF/+slm93Ya8DTT2EJjAG4lgmk8ze1sy177xYrBuHEwcKD68kiBkKurnrvl4H+KQ4cOKdkREcmJ5cvhuefgwH+jsVJx43P6M5JxxFICC2k8wye8xasEEJe99h94AGbPhpIlczlwkbyVq3120tLSsr3ll0RHRMTltGpldmB+6y3rExh30niaT9lHTR7jGwzcmMIgarKP2XTL3tw8P/8Mt99unkOkENBzTBERV+TlBa+8AitWQIkS1t3liOEberGK+6jOfk4TyKPM5gF+5iBVs97+n39CWBh8/70DghdxLUp2RERcWYsWsG0b1Kljs/s+1rCDerzJq3iTxC/cTx12MYbRJOGdtbYvXICuXc2kKiXFAcGLuAYlOyIirq5KFQgPNxOTq3hziVd5m93UpiXLuYQ3rzOGuuxkJS2y3v7YsXDXXXDwYC4HLuIaspzsnDp1ypFxiIhIZnx9zTl5Jk82X3Fd5Wb+ZBmtmctDBHGKQ1TjAVbSnVlEEZi19rdsgQYN4LPPNCePFDhZTnZq167NrFmzHBmLiIhkxmIxh41v2gTVqtkWAQ8xn33U5Dkm4kYqc+hOTfYxmYFZm5snMRGeego6dIC//nLMNYg4QZaTnbfffpunnnqKhx56iLNnzzoyJhERycwtt5hrbPXoka7Ij/NMZChbuZ3b2UI8/gxmMo3ZTAS3Zq39RYugZk344AO4fDmXgxfJe1lOdp555hl27NjBmTNnCA0NZdGiRY6MS0REMlO8OHzzDXz9tc1orStuZTvhhDGFZ/AnlghuoxFbGMxHxJGF+cQSEmDYMHPx0g0bHHABInknRwuBTp48maFDh1KrVi08PGwXTv/9999zLbi8okkFRSRfO33anITwf//LsDiacgzjfWZhPgkKJIqJDOFh5mLJ6jl69YJ339WiouJScnUG5asdO3aMvn37smvXLp566ql0yc7o0aNzFrETKdkRkQLhp5/gmWfs9rdZxX0MYCoHqQ7A/fzMJzxDVQ5nrf2AAHj7bbNfj7t7LgUtknMOSXY+//xzhg0bRosWLfj0008pU6ZMrgTrbEp2RKTAiI83k5E5czIsTsKbCbzIWF4mmSJ4k8TLjGUE7+DNpaydo2FDmDrVnIVZxIlydbkIgFatWjFixAgmT57MggULCkyiIyJSoPj5waxZ8M475uitaxQhmVG8yS7qcD8/k0wRRvMGddnJLzTP2jkiIqBxYxgwwEyuRFxclpOd1NRUduzYQa9evRwZj4iI3CiLBV58EZYuNV89ZaAqh1lBS+bwCIFEcZDq3M8v9OBboil3/XMYBkybBvXqwZo1uRu/SC7LcrKzcuVKypcv78hYREQkN7VqZU4WeM1SE1dYgEeYyz5qMpiPcCOVWfSgJvv4hAFZm5vn2DG47z6zg3RiYu7GL5JLtFyEiEhBVq0abN0KQ4bYreJPPB/xHJtpTEO2EUcAA/mEMML5nVuydp6PPjLn/9m0KXfiFslFSnZERAq6IkXgww9h5UoIDrZb7TYi2ExjJjMQP+LYSiNuZyvPMZF4il//PAcOmGtsjR8PaWm5eAEiN0bJjohIYdGiBezcCQ8/bLeKO2kM5BP2UZNuzCYNdz7iOWqyj7k8xHWH76alwciR5pITmm1fXIRTk53169fTvn17goODsVgs/PDDDzblffr0wWKx2GytWrWyqXP27Fl69OiBn58fAQEB9OvXj4SEhDy8ChGRfKRkSXNY+pw5EBRkt1oQ0czmUX7mfqpykCiCeYS5tGYZh6ly/fMsXmzOvrx1ay4GL5IzTk12Lly4QP369ZkyZYrdOq1atSIqKsq6zZ4926a8R48e7N69m5UrV7J48WLWr1/Pk08+6ejQRUTyL4sFHnkE9u2DoUMznSDwfn5hJ3UZzRi8SGYFrajDLt7iFZLxsnscYHZevvtuc9SWiBPlaLkIR7BYLCxcuJCOHTta9/Xp04fY2Nh0T3yu2Lt3L6GhoWzdupXbbrsNgOXLl9OmTRtOnjxJcCbvpq+mSQVFpFDbscOcefm33zKtdoBqDGQKv3A/ADXYxyc8w31kYej5s8+aC4tq5mXJRbk+qaCzrF27lrJly1KjRg0GDBjAmTNnrGXh4eEEBARYEx2AFi1a4ObmxubNm+22mZycTHx8vM0mIlJo1asH69fDV19BJhPGVucgP/MAs+lGOaLZT02as5qefM1prrNm1kcfQceO5gKjInnMpZOdVq1a8fXXX7Nq1Sreeecd1q1bR+vWrUlNTQUgOjqastcsSufh4UHJkiWJjo622+64cePw9/e3biEhIQ69DhERl+fmBn36wP795szIGcy+DObcPN34H/uoyUAmYyGNb+lJTfYxjadIy2xp0cWLoUkTOHnSIZcgYo9LJzvdunXjwQcfpG7dunTs2JHFixezdetW1q5de0Ptjhw5kri4OOt24sSJ3AlYRCS/K1ECPvkENm+Gq56aXyuAOCYzmM005lYiiKUEA5hGGOFsp4H99iMjzfl45s41Z2EWyQMunexcq0qVKpQuXZpDhw4BEBgYSExMjE2dlJQUzp49S2BgoN12vL298fPzs9lEROQqt99uThA4caI5T4+9amxjC434iMEUJ54tNOY2tjGUDziPb8YH/fOP2UG6a1c4fdox8YtcJV8lOydPnuTMmTME/TtcMiwsjNjYWCIiIqx1Vq9eTVpaGo0bN3ZWmCIiBYO7u7kMRGSkufCnvWqkMZjJ7KMmjzCHNNyZyFBqsZf5dLE/N8+CBRAaai5cqqc84kBOTXYSEhKIjIwkMjISgCNHjhAZGcnx48dJSEhg+PDhbNq0iaNHj7Jq1So6dOhA1apVadmyJQC1atWiVatW9O/fny1btvDbb78xaNAgunXrluWRWCIich01asCvv8LYseDpabdaMFHMoTvLacnNHOIvyvMQ82nLEv6kcsYHnT0LPXrAE0/ApUsOugAp9AwnWrNmjQGk23r37m0kJiYaDzzwgFGmTBnD09PTqFixotG/f38jOjrapo0zZ84Y3bt3N3x9fQ0/Pz+jb9++xvnz57MVR1xcnAEYcXFxuXl5IiIFT2SkYVSvbhjmsxi7WyJFjNd43fAiyQDDKEKi8RYvG0l42T+ubVvDuHjR2Vco+UhWv79dZp4dZ9I8OyIi2RAfD/37m52Mr2M/1XmGT1hNcwBqspepDKAp6zI+oHlz+PFHKFYsNyOWAqrAzLMjIiIuxs/PXG7i448zfa0FUIMD/EILvuNRynKafdSiGWvpxUxiyGBOn1WroFUrM6ESySVKdkREJPssFhg0CDZsgExGv4I5N8+jzGYfNRnAJ1hI4xt6UZN9fEb/9HPz/PqruWjpP/84Ln4pVJTsiIhIzjVubM7JU7fudauWIJZPGMgm7uAWfuccJXmKz7iL3/iDeraVt24119U6dsxBgUthomRHRERuTIUK5tOYVq2yVL0RW9lCIybyHMWJZxNhNCSCYbxnOzfP/v1w552wc6eDApfCQsmOiIjcOD8/WLTIXFA0CzxI5Tk+Yi+1eIi5pOLBBwwjlD0soNN/c/OcOmUuMbEmC4uNitihZEdERHKHhwdMmQLffmsuO5EFN3GKuTzCUlpTmT85SQhdWEB7FnGESmaluDhzlNbIkZqLR3JEyY6IiOSuHj1g927o0CHLh7RmObupzau8iSeXWEI7arObcbzEJTzBMEgdP4G1tQcy+92TrF0L/64JLXJdmmcHzbMjIuIQhmEOUR80yJwpOYv2UYMBTGUtzQCoxR66M4vPeIqThFjrlS9vMGmShc6dcz1yySey+v2tZAclOyIiDhUVBS1bZqujsQF8Rw+e5wP+puxVe/8bpm4hDSwW5s9XwlNYaVJBERFxDUFBZgfjhg2zfIgFeIzv2EMtipHAtYkOgIEbGAZDnkggNaXQ/7tdMqFkR0REHK9UKXN25DvvzNZhu6jLBXy5NtG5wsCNE+d82XD3SPMJkkgGlOyIiEje8PeHFSvMkVVZFEVQ1uptPgYNGkB4eA6Dk4JMyY6IiOQdX19YsgQefjhL1YPI2tOaIKIgJgaaNcvSAqVSuCjZERGRvOXtDbNnm6O0rqMJGyjPCbMzcoYMyhFFEzaYvyYnwyOPwNix5mgwEZTsiIiIM7i5wUcfwdtvZ1rNnTQm8RxABgmP2Wn5Mp6cIti26JVX4PHHNQmhAEp2RETEWSwWePllmDkTihSxW60zC5lPV27iL5v9N/EX5TnOWUrTmmXE4m974IwZ0LmzEh5RsiMiIk7WqxdERJgdjO3ozEKOUok1NGUW3VlDU45Rkd+4myBOsZs6dGIhyXjZHnilf9Dly469BnFpSnZERMT5QkNh82Z46SXziU8G3EmjKevozhyasg530qjACZbShuLEs5Zm9GEGadcOU//xR+jeXQlPIaZkR0REXIOXF4wbB+vWQcWKWT6sAX/wPV3w4DJz6M5IxqWv9P338NhjkJKSiwFLfqFkR0REXEuTJrBjh9nBOIvu5xemY9afwAgmMzB9pblzoV07iI/PrUgln1CyIyIirsfPD778EhYtgsDALB3Sk295m5cBeJaPWEjH9JVWrIC77oJjx3IxWHF1SnZERMR1tWsHu3ZBt25Zqj6ScTzFNAzceJRZbCQsfaVdu6BxY9i6NZeDFVelZEdERFxbqVLmJISffQbu7plWtQCTGUQ7FpFEUdqziP1UT1/x9Gm4915zvS4p8JTsiIhI/tC/PyxbZr7iyoQHqcyhG7ezhbOUojXLOE3Z9BUvXoQHH4Q//nBQwOIqlOyIiEj+cf/9sHHjdUdrFSORxbSjCoc5QhXasoQEiqWvmJgIHTvCmTOOiVdcgpIdERHJX2rXhk2boGbNTKuV5W+W04rS/E0Et/EI/yOFDF6DHT1qTjyoYekFlpIdERHJfwIDYfVqqJ5Bf5yrVOMQi2lHURJZSlue4RMyXB509Wp46iktLVFAKdkREZH8KSjITFJuvjnTao3Zwhy64UYqn/Mkb/NKxhWnTzeHpR8+7IBgxZmU7IiISP51002wZg1UrpxptQdZxGQGAfAabzGD3hlX3LYNbr3VnIBQCgwlOyIikr+FhMCGDXDLLZlWG8A0Xvp3KYn+fM4KHsi4Ynw8PPIIDB6sfjwFhJIdERHJ/266Cdavh/btM632Nq/Qg29JwZOuzGc7DexXnjzZnMxQ/XjyPSU7IiJSMPj6wsKFMGSI3SpuGEznce5jFQkUpw1LOUYF+21+/z107gxJSbkfr+QZJTsiIlJwuLvDhx/C22/breLFZRbQmbrsIJogWrOMs5Sw3+aSJeYTowsXHBCw5AUlOyIiUvCMHAmPPmq32J94ltKG8pxgL6F05AeS8Lbf3i+/QOvWSnjyKSU7IiJS8Fgs8PnnmXZaLs9fLKM1/sSygXvoxdekYbHf5oYNZh+etDQHBCyOpGRHREQKJh8f+OEHKFPGbpU67GYhnfDkEvN4mBd4L/M2Fy+GN9/M3TjF4Zya7Kxfv5727dsTHByMxWLhhx9+sCk3DINRo0YRFBRE0aJFadGiBQcPHrSpc/bsWXr06IGfnx8BAQH069ePhISEPLwKERFxWRUqmBMPlitnt0oz1jKDPgB8yPN8yJDM2xwzxkx6JN9warJz4cIF6tevz5QpUzIsnzBhAh999BHTpk1j8+bNFCtWjJYtW5J0Va/4Hj16sHv3blauXMnixYtZv349Tz75ZF5dgoiIuLo6dcy1tBo1slvlUWbzDi8CMIz3mUfXzNt87DE4dCg3oxRHMlwEYCxcuND6e1pamhEYGGi8++671n2xsbGGt7e3MXv2bMMwDGPPnj0GYGzdutVaZ9myZYbFYjH++usvu+dKSkoy4uLirNuJEycMwIiLi8v9CxMREdeQnGwYL7xgGJDhlgbGQD42wDC8uWis5267dQ0wjEqVDOPQIWdfVaEWFxeXpe9vl+2zc+TIEaKjo2nRooV1n7+/P40bNyY8PByA8PBwAgICuO2226x1WrRogZubG5s3b7bb9rhx4/D397duISEhjrsQERFxDV5e8O675lDyUqXSFVuASTxHRxaSTBE68CN7yWRl9aNH4e67Ydcuh4UsucNlk53o6GgAyl3znrVcuXLWsujoaMqWLWtT7uHhQcmSJa11MjJy5Eji4uKs24kTJ3I5ehERcVlt2sDvv2e4Yro7acziUcLYyDlK0pplRBFov63oaLjnHtiyxYEBy41y2WTHkby9vfHz87PZRESkEKlQwVxeok6ddEVFSeInHqQaBzhGJdqyhPP42m/r3Dlo3lwJjwtz2WQnMNDMpE+fPm2z//Tp09aywMBAYmJibMpTUlI4e/astY6IiEiGypWDtWvNVc6vUZozLKcVZTnNdm7lIeZxGQ/7bSUkmLMsR0U5Ll7JMZdNdipXrkxgYCCrVq2y7ouPj2fz5s2EhYUBEBYWRmxsLBEREdY6q1evJi0tjcaNG+d5zCIiks+UKgWrVkGtWumKqnCExbTDhwusoBVP8hlGZm3FxEDfvpp00AU5NdlJSEggMjKSyMhIwOyUHBkZyfHjx7FYLAwZMoS33nqLn376iZ07d9KrVy+Cg4Pp2LEjALVq1aJVq1b079+fLVu28NtvvzFo0CC6detGcHCw8y5MRETyj4AAcwHR4sXTFd3ONubyMG6kMoO+jGFM5m2tWGGuli6uJY9Gh2VozZo1BpBu6927t2EY5vDz1157zShXrpzh7e1tNG/e3Ni/f79NG2fOnDG6d+9u+Pr6Gn5+fkbfvn2N8+fPZyuOrA5dExGRAmzhQrvDzD+lv/XXz+mX+ZB0b2/D2LHD2VdTKGT1+9tiGEamT+UKg/j4ePz9/YmLi1NnZRGRwuzVV+2umP4ab/AWr+FOCj/xIG1YZr+doCBYvhzq1XNQoAJZ//522T47IiIiee71183FPjPwBqPozQxS8eAh5rGNhvbbiYoyh6SvW+egQCU7lOyIiIhc4e4O33wDvXqlK7IAn9Of+/mZRIrRliX8SWX7bcXFQcuWsGCB4+KVLFGyIyIicjUPD/jqKxg8OF2RJynMpysN2E4M5WjNMv4h/WzMVsnJ8NBDSnicTMmOiIjItdzcYNIkGDUqXZEf51lCWypwjAPU4EF+4iJF7LeVlgY9esAffzgwYMmMkh0REZGMWCxmH56ePdMVBRPFMloTwDnCuZMefEdqZl+pSUnmE574eAcGLPYo2REREcnM5MlQOX3fnFD28iMd8CKZhXRmCBMzn3Tw4EF44glzgLrkKSU7IiIimfHzMzstu6X/yryHDXyD+eRnMoN5n2GZtzVvniYddAIlOyIiItdz110wZkyGRQ8zjw8YCsBw3mMOj2Te1pAhMGtW7sYnmVKyIyIikhWvvmqufZWBoUxkCB8C0JuZrOVe++2kpZn9gObOdUSUkgElOyIiIllhscAXX8Czz2ZY/D7D6Mo8LuFNR35gF7Xtt5WWBo8+Ct9/76Bg5WpKdkRERLLKzQ0mToRx49IXYfANPbmbDcQRQGuW8ReZLEqdmmrO1rx0qePiFUDJjoiISPZYLPDSS/Dll+k6LRchmR/pQE32cpIQ2rCUODJZczElBR5+GCIjHRtzIadkR0REJCcefxzGjk23uyTnWEZrAoliB/XpwvdcwtN+OxcuQNu2cPKkA4Mt3JTsiIiI5NTw4Waico1KHGMJbfHlPKtoQT++zHwOnlOnoH17OH/eYaEWZkp2REREcsrNDWbOhAoV0hXdynbm0xV3UviWnrzKW5m3FRlp9uFJSXFMrIWYkh0REZEbUaqUOYzc2ztdUUt+5nP6AzCWV5jGU5m3tXQpPP20ZlnOZUp2REREblTjxjBnjrli+jX6MoPXMRcUHcgUfqJ95m19+aW5JpfkGiU7IiIiuaFjR5g9G9zd0xW9xpv04wvScKcbc9hMo8zbev11+Owzx8RZCCnZERERyS1du8J336Ubkm4BpjKA1izlIj60YzGHuDnztgYMgCVLHBdrIaJkR0REJDc98gh8/bU5H89VPElhLg9zKxH8QxlasZy/KW2/nbQ0eOwxOHrUsfEWAkp2REREcluPHvD+++l2+3KBJbSlEkc4TFXasZhEitpvJzbWnHTw0iXHxVoIKNkRERFxhCFD4Jln0u0O5DTLaUVJzrCFxnRjDimk7+djtXWrOZ+P5JiSHREREUewWGDSJGjdOl1RDQ6wiPYU4SKLeJDBfJz5pIMffWS+GpMcUbIjIiLiKB4e8L//Qf366YruJJzv6IGFNKYxgPG8lHlbffqYw9Il25TsiIiIOFLx4uZkgRnMstyZhUziOQBeZhzf8Jj9dgwDnngCPv7YUZEWWEp2REREHC04GJYvhxIl0hUNZjIv8C4AjzOdX2ieeVvPPgsTJjgiygJLyY6IiEheqFULFi+GIkXSFb3DCB5hDil40pkF/EG9zNsaMQKmTXNQoAWPkh0REZG8cued8O236Xa7YTCT3tzLWs7jRxuWcoLymbc1eDCsX++gQAsWJTsiIiJ5qUsXeP75dLu9ucQPdKQ2uzjFTbRmGbH4228nJcWcsfn4cQcGWzAo2REREclr48fDHXek2x1AHEtpQzB/sZs6dOQHkvGy387ff5trciUmOi7WAkDJjoiISF7z9DSHpJcsma6oAidYShuKE886mtKHGaRhyaCRf23fbg5LT011XLz5nJIdERERZ6hQAebMMefiuUZ9drCAznhwmTl05yXGZ97WvHnw3HPm8HRJR8mOiIiIs9x/P8yenWHC04JVTOdxAN7lRT5mUOZtTZkCb77piCjzPSU7IiIiztS1KyxcCF7p++b05Fve5mUAnmMSC+mYeVujR8PUqQ4IMn9TsiMiIuJs7dqZc/AUTb8C+kjG8RTTMHDjUWaxkbDM2xo4EH7+2UGB5k8uneyMGTMGi8Vis9WsWdNanpSUxMCBAylVqhS+vr506dKF06dPOzFiERGRHLr/fnOWZW9vm90WYDKDaM9PJFGU9ixiP9Xtt2MY0LMnxMQ4Nt58xKWTHYDatWsTFRVl3X799Vdr2dChQ1m0aBHz5s1j3bp1nDp1is6dOzsxWhERkRtwzz3w2WfpdnuQymy604jNnKUUrVlGNOXstxMTA48/rg7L/3L5ZMfDw4PAwEDrVrp0aQDi4uL48ssv+eCDD7jvvvto2LAhX331FRs3bmTTpk1OjlpERCSHevWCIUPS7S5GIotoz80c4ghVaMdiEihmv50lS7SkxL9cPtk5ePAgwcHBVKlShR49enD835kiIyIiuHz5Mi1atLDWrVmzJhUqVCA8PDzTNpOTk4mPj7fZREREXMa770Lz9AuCluVvltOK0vxNBLfxMHNJwd1+O88/D5GRjoszn3DpZKdx48bMmDGD5cuXM3XqVI4cOUKTJk04f/480dHReHl5ERAQYHNMuXLliI6OzrTdcePG4e/vb91CQkIceBUiIiLZ5OFhTjpYPX3fnKocZjHtKEoiy2jDAKZi92VVUhK0agWHDjk0XFfn0slO69ateeihh6hXrx4tW7Zk6dKlxMbGMnfu3Btqd+TIkcTFxVm3EydO5FLEIiIiuaRUKVixAoKD0xU1Zgtz6IYbqXxBf97iVfvtnD5tdn7+6y8HBuvaXDrZuVZAQADVq1fn0KFDBAYGcunSJWJjY23qnD59msDAwEzb8fb2xs/Pz2YTERFxOZUqmQlPiRLpih5kEZP/nWhwFG8yg9722zl6FFq2hLNnHROni8tXyU5CQgKHDx8mKCiIhg0b4unpyapVq6zl+/fv5/jx44SFXWcOAhERkfyiTh27c/AMYBojGQtAfz5nBQ/Yb2f3bujc2VwtvZBx6WTnhRdeYN26dRw9epSNGzfSqVMn3N3d6d69O/7+/vTr14/nn3+eNWvWEBERQd++fQkLC+OODFaSFRERybfuvBM+/TTDord5hcf4hhQ86cp8fucW++2sWwfvvOOgIF2XSyc7J0+epHv37tSoUYOHH36YUqVKsWnTJsqUKQPAhx9+SLt27ejSpQv33HMPgYGBLFiwwMlRi4iIOMBjj0G3bul2W4Av6UdzfiGB4rRlCUepaL+dN96AffscF6cLshiGZhyKj4/H39+fuLg49d8RERHXFRsL9epBBgNr4vCjCRvYST1qspffuIuSnMu4nTvvhA0bwM2ln3lcV1a/v/P3VYqIiBQmAQHwzTdgsaQr8ieepbShPCfYRy068CNJeKdvA2DjRnOV9EJCyY6IiEh+cu+9dvvdlOcvltEaf2L5lSb05BvSSJ8YAfDiixAR4cBAXYeSHRERkfzmhRdg6NAMi+qwm4V0wotk5vMQL/Bexm0kJUGnToViwVAlOyIiIvmNxQLvvQe9M55bpxlrmUEfAD7keT5kSMbtnDgBXbvCpUuOidNFKNkRERHJj9zc4IsvoEOHDIu7M4cJDAdgGO8zj64Zt7Nhg7mGVgGmZEdERCS/8vCAWbOgatUMi1/gPQbxMQZuPMa3rKdJxu1MmWLO1FxAKdkRERHJz3x84PPPMyyyABMZQkcWcglvOvAje6iVcTtDhhTY2ZWV7IiIiOR3TZvCU09lWOROGrN4lDA2EksJWrOMUwSlr7hvn92kKb9TsiMiIlIQvPMO3HRThkVFSeInHqQaBzhORdqyhPP4pq84ejTExzs40LynZEdERKQg8PeH774Dd/cMi0tzhuW0oiynieQWujKfy3jYVvr7b+jRo8CNzlKyIyIiUlDcey98+KHd4iocYQlt8eECP9OSJ/mMdGtGLV5srsF1+bJDQ81LSnZEREQKkkGDoG9fu8W3EcE8HsKdFGbQl9G8nr7SwoXmE54C0mFZyY6IiEhBYrHAJ59A48Z2q7RhGVMZAMCbjOJznkhfad48s9NzAVgvXMmOiIhIQVOkCHz1ld3+OwD9+YLXeAOAAUxlCW3SV5o+vUCM0FKyIyIiUhDVqgX9+2da5XVG05sZpOLBw8xlGw3TVxo6FPbvd1CQeUPJjoiISEE1Zgz4ZjDE/F8W4HP6cz8/k0gx2rKEP6lsWykxMd+P0FKyIyIiUlCVKwevvpppFU9S+J4uNGA7MZSjFcv5h1K2lSIiYMSIfNt/R8mOiIhIQfbCC9CsWaZVipPAUtpQgWMcpDoP8hOJFLWtNHEiPPsspKU5LlYHUbIjIiJSkLm7w//+B6GhmVYLIprltCKAc4RzJz34jtRr04TJk6F7d0hOdmDAuU/JjoiISEFXpgysXg01a2ZarRb7+IkH8SKZH+jEc0xKP+ng3LnQtm2+WlZCyY6IiEhhUK6cmfBUq5ZptSb8yrc8BsAUBvEeL6SvtGqV+Wrs9GlHRJrrlOyIiIgUFkFBsGaN3QVDr3iI+XzAUABe5F1m0y19pd9/h7vugsOHHRFprlKyIyIiUpjcdBN8840503ImhjKRIZjrbPVmJmtomr7S4cNw552wfbsDAs09SnZEREQKm2bNzFFa1/E+w+jKPC7jRScWsova6SvFxJgLkK5dm/tx5hIlOyIiIoXRm2/CLbdkWsUNg2/oyd1sII4AWrOMk2TwCuz8eWjdGpYvd1CwN0bJjoiISGHk7Q1z5kDJkplWK0IyP9KBmuzlJCG0YSlx+KWvmJQEDz4IP/zgmHhvgJIdERGRwqp6dViyBPz9M61WknMsozWBRLGTenRmAZfwTF/x8mXo2tVMolyIkh0REZHC7I47YMMGCA7OtFoljrGUNvhyntU053Gmp5+DByA1FR591JzI0EUo2RERESns6taFjRuhRo1Mq91CJPPpijspfMdjvMLbGVc0DOjVy1xTywUo2RERERGoWBF+/RUaNcq0Wkt+5nP6AzCOl5nK0xlXvHQJHnoIYmNzOdDsU7IjIiIiptKlzVmWW7XKtFpfZvA6owAYxGR+on3GFY8cMZ/wOHm1dCU7IiIi8p9ixeCnn+CxxzKt9hpv8gSfk4Y73ZjDJhpnXHHRIujQwamLhyrZEREREVuenjBzJgwbZreKBZjKANqwhIv40J5FHKRqxpUXLYLGjWH3bsfEex1KdkRERCQ9Nzd47z2YMMFuFQ9S+R+P0JBt/EMZWrOMGMpkXPmPP6BhQ5g4EdLSHBOzHUp2RERExL7hw+GTT+wW+3KBJbSlMn9ymKq0YzEX8Mm4cnIyDB0KLVvCX385KOD0lOyIiIhI5gYMgBkzzKc9GShHDMtoTUnOsJVGdGMOKbjbb++XX8zh7nPnOibeaxSYZGfKlClUqlSJIkWK0LhxY7Zs2eLskERERAqO3r1h9mxwzziJqcEBFtGeIlxkMe0ZxGRScGMt9zKbbqzlXlKvTjvOnYNHHoGePSEuzqGhF4hk53//+x/PP/88o0eP5vfff6d+/fq0bNmSmJgYZ4cmIiJScDz8MLz7rt3iOwlnFo9iIY1PeZpSnKEZa3mU2TRjLZU4ygI62R707bdQrx6sW+ewsAtEsvPBBx/Qv39/+vbtS2hoKNOmTcPHx4fp06c7OzQREZGCZcgQ6NjRbnEnfuBxvgQgngCbsr+4ia7MT5/wHD8O+/fnbpxXyffJzqVLl4iIiKBFixbWfW5ubrRo0YLw8PAMj0lOTiY+Pt5mExERkSywWOCrr6By5QyLU3FjBa0gg5WzjH/TjiFMtH2l1a4d9O/viGiBApDs/PPPP6SmplKuXDmb/eXKlSM6OjrDY8aNG4e/v791CwkJyYtQRURECoaAAJg/H3x90xVtoAknCcGciSc9AzdOUIENNDF3lCkDX3xhJlEOku+TnZwYOXIkcXFx1u3EiRPODklERCR/ufVWWLkS/P1tdkcRlKXDrfW+/BKueWCR2/J9slO6dGnc3d05ffq0zf7Tp08TGBiY4THe3t74+fnZbCIiIpJNd9wBmzaZkwX+K4ioLB0aRJT56qq9nXW1clG+T3a8vLxo2LAhq1atsu5LS0tj1apVhIWFOTEyERGRQqBmTdi4EV55BdzcaMIGynMCCxnPkmwhjRCO0+TmKPjggzwJMd8nOwDPP/88n3/+OTNnzmTv3r0MGDCACxcu0LdvX2eHJiIiUvB5ecFbb8H69bhXrsgkngNIl/Bc+X2i2/O4f/d1hn1+HKFAJDuPPPII7733HqNGjaJBgwZERkayfPnydJ2WRURExIHuugsiI+ncx5/5dOUmbJeEKM9J5tOVzq/VMRcGzSMWwzDSjw0rZOLj4/H39ycuLk79d0RERHLDggWkPvEUG87VJooggoiiCRtwb3Qb/PqrubL6Dcrq97fHDZ9JRERE5FqdO+N+xx00ffxxWDHH3OfjY86YnAuJTnYUiNdYIiIi4oKCg2HZMvj4YyhSBD78EKpVy/Mw9GRHREREHMdigUGDzCHmFSo4JQQlOyIiIuJ4FSs67dR6jSUiIiIFmpIdERERKdCU7IiIiEiBpmRHRERECjQlOyIiIlKgKdkRERGRAk1Dz4ErK2bEx8c7ORIRERHJqivf29db+UrJDnD+/HkAQkJCnByJiIiIZNf58+fx9/e3W66FQIG0tDROnTpF8eLFsVgsudZufHw8ISEhnDhxotAuMKp7oHsAugege1DYrx90DyD374FhGJw/f57g4GDc3Oz3zNGTHcDNzY3y5cs7rH0/P79C+8G+QvdA9wB0D0D3oLBfP+geQO7eg8ye6FyhDsoiIiJSoCnZERERkQJNyY4DeXt7M3r0aLy9vZ0ditPoHugegO4B6B4U9usH3QNw3j1QB2UREREp0PRkR0RERAo0JTsiIiJSoCnZERERkQJNyY6IiIgUaEp2HGjKlClUqlSJIkWK0LhxY7Zs2eLskBxizJgxWCwWm61mzZrW8qSkJAYOHEipUqXw9fWlS5cunD592okR37j169fTvn17goODsVgs/PDDDzblhmEwatQogoKCKFq0KC1atODgwYM2dc6ePUuPHj3w8/MjICCAfv36kZCQkIdXcWOudw/69OmT7nPRqlUrmzr5+R6MGzeO22+/neLFi1O2bFk6duzI/v37bepk5bN//Phx2rZti4+PD2XLlmX48OGkpKTk5aXkWFbuQdOmTdN9Dp5++mmbOvn5HkydOpV69epZJ8kLCwtj2bJl1vKC/hmA698DV/gMKNlxkP/97388//zzjB49mt9//5369evTsmVLYmJinB2aQ9SuXZuoqCjr9uuvv1rLhg4dyqJFi5g3bx7r1q3j1KlTdO7c2YnR3rgLFy5Qv359pkyZkmH5hAkT+Oijj5g2bRqbN2+mWLFitGzZkqSkJGudHj16sHv3blauXMnixYtZv349Tz75ZF5dwg273j0AaNWqlc3nYvbs2Tbl+fkerFu3joEDB7Jp0yZWrlzJ5cuXeeCBB7hw4YK1zvU++6mpqbRt25ZLly6xceNGZs6cyYwZMxg1apQzLinbsnIPAPr372/zOZgwYYK1LL/fg/LlyzN+/HgiIiLYtm0b9913Hx06dGD37t1Awf8MwPXvAbjAZ8AQh2jUqJExcOBA6++pqalGcHCwMW7cOCdG5RijR4826tevn2FZbGys4enpacybN8+6b+/evQZghIeH51GEjgUYCxcutP6elpZmBAYGGu+++651X2xsrOHt7W3Mnj3bMAzD2LNnjwEYW7dutdZZtmyZYbFYjL/++ivPYs8t194DwzCM3r17Gx06dLB7TEG7BzExMQZgrFu3zjCMrH32ly5dari5uRnR0dHWOlOnTjX8/PyM5OTkvL2AXHDtPTAMw7j33nuN5557zu4xBe0eGIZhlChRwvjiiy8K5Wfgiiv3wDBc4zOgJzsOcOnSJSIiImjRooV1n5ubGy1atCA8PNyJkTnOwYMHCQ4OpkqVKvTo0YPjx48DEBERweXLl23uRc2aNalQoUKBvRdHjhwhOjra5pr9/f1p3Lix9ZrDw8MJCAjgtttus9Zp0aIFbm5ubN68Oc9jdpS1a9dStmxZatSowYABAzhz5oy1rKDdg7i4OABKliwJZO2zHx4eTt26dSlXrpy1TsuWLYmPj7f5V3F+ce09uOK7776jdOnS1KlTh5EjR5KYmGgtK0j3IDU1lTlz5nDhwgXCwsIK5Wfg2ntwhbM/A1oI1AH++ecfUlNTbf7DAZQrV459+/Y5KSrHady4MTNmzKBGjRpERUXx+uuv06RJE3bt2kV0dDReXl4EBATYHFOuXDmio6OdE7CDXbmujP77XymLjo6mbNmyNuUeHh6ULFmywNyXVq1a0blzZypXrszhw4d5+eWXad26NeHh4bi7uxeoe5CWlsaQIUO46667qFOnDkCWPvvR0dEZfk6ulOUnGd0DgEcffZSKFSsSHBzMjh07GDFiBPv372fBggVAwbgHO3fuJCwsjKSkJHx9fVm4cCGhoaFERkYWms+AvXsArvEZULIjN6x169bWP9erV4/GjRtTsWJF5s6dS9GiRZ0YmThTt27drH+uW7cu9erV4+abb2bt2rU0b97ciZHlvoEDB7Jr1y6bvmqFjb17cHUfrLp16xIUFETz5s05fPgwN998c16H6RA1atQgMjKSuLg45s+fT+/evVm3bp2zw8pT9u5BaGioS3wG9BrLAUqXLo27u3u6HvenT58mMDDQSVHlnYCAAKpXr86hQ4cIDAzk0qVLxMbG2tQpyPfiynVl9t8/MDAwXWf1lJQUzp49W2DvS5UqVShdujSHDh0CCs49GDRoEIsXL2bNmjWUL1/euj8rn/3AwMAMPydXyvILe/cgI40bNwaw+Rzk93vg5eVF1apVadiwIePGjaN+/fpMmjSpUH0G7N2DjDjjM6BkxwG8vLxo2LAhq1atsu5LS0tj1apVNu8wC6qEhAQOHz5MUFAQDRs2xNPT0+Ze7N+/n+PHjxfYe1G5cmUCAwNtrjk+Pp7NmzdbrzksLIzY2FgiIiKsdVavXk1aWpr1L4KC5uTJk5w5c4agoCAg/98DwzAYNGgQCxcuZPXq1VSuXNmmPCuf/bCwMHbu3GmT9K1cuRI/Pz/rKwBXdr17kJHIyEgAm89Bfr4HGUlLSyM5OblQfAbsuXIPMuKUz0CudHOWdObMmWN4e3sbM2bMMPbs2WM8+eSTRkBAgE1v84Ji2LBhxtq1a40jR44Yv/32m9GiRQujdOnSRkxMjGEYhvH0008bFSpUMFavXm1s27bNCAsLM8LCwpwc9Y05f/68sX37dmP79u0GYHzwwQfG9u3bjWPHjhmGYRjjx483AgICjB9//NHYsWOH0aFDB6Ny5crGxYsXrW20atXKuOWWW4zNmzcbv/76q1GtWjWje/fuzrqkbMvsHpw/f9544YUXjPDwcOPIkSPGL7/8Ytx6661GtWrVjKSkJGsb+fkeDBgwwPD39zfWrl1rREVFWbfExERrnet99lNSUow6deoYDzzwgBEZGWksX77cKFOmjDFy5EhnXFK2Xe8eHDp0yHjjjTeMbdu2GUeOHDF+/PFHo0qVKsY999xjbSO/34OXXnrJWLdunXHkyBFjx44dxksvvWRYLBbj559/Ngyj4H8GDCPze+AqnwElOw708ccfGxUqVDC8vLyMRo0aGZs2bXJ2SA7xyCOPGEFBQYaXl5dx0003GY888ohx6NAha/nFixeNZ555xihRooTh4+NjdOrUyYiKinJixDduzZo1BpBu6927t2EY5vDz1157zShXrpzh7e1tNG/e3Ni/f79NG2fOnDG6d+9u+Pr6Gn5+fkbfvn2N8+fPO+Fqciaze5CYmGg88MADRpkyZQxPT0+jYsWKRv/+/dMl+/n5HmR07YDx1VdfWetk5bN/9OhRo3Xr1kbRokWN0qVLG8OGDTMuX76cx1eTM9e7B8ePHzfuueceo2TJkoa3t7dRtWpVY/jw4UZcXJxNO/n5Hjz++ONGxYoVDS8vL6NMmTJG8+bNrYmOYRT8z4BhZH4PXOUzYDEMw8idZ0QiIiIirkd9dkRERKRAU7IjIiIiBZqSHRERESnQlOyIiIhIgaZkR0RERAo0JTsiIiJSoCnZERERkQJNyY6IiIgUaEp2RKRAqVSpEhaLBYvFkm4Bxuxq2rSpta0r6/mISP6jZEdEXE5qaip33nknnTt3ttkfFxdHSEgIr7zySqbHv/HGG0RFReHv739DcSxYsIAtW7bcUBsi4nxKdkTE5bi7uzNjxgyWL1/Od999Z90/ePBgSpYsyejRozM9vnjx4gQGBmKxWG4ojpIlS1KmTJkbakNEnE/Jjoi4pOrVqzN+/HgGDx5MVFQUP/74I3PmzOHrr7/Gy8srW23NmDGDgIAAFi9eTI0aNfDx8aFr164kJiYyc+ZMKlWqRIkSJXj22WdJTU110BWJiLN4ODsAERF7Bg8ezMKFC+nZsyc7d+5k1KhR1K9fP0dtJSYm8tFHHzFnzhzOnz9P586d6dSpEwEBASxdupQ///yTLl26cNddd/HII4/k8pWIiDMp2RERl2WxWJg6dSq1atWibt26vPTSSzlu6/Lly0ydOpWbb74ZgK5du/LNN99w+vRpfH19CQ0NpVmzZqxZs0bJjkgBo9dYIuLSpk+fjo+PD0eOHOHkyZM5bsfHx8ea6ACUK1eOSpUq4evra7MvJibmhuIVEdejZEdEXNbGjRv58MMPWbx4MY0aNaJfv34YhpGjtjw9PW1+t1gsGe5LS0vLcbwi4pqU7IiIS0pMTKRPnz4MGDCAZs2a8eWXX7JlyxamTZvm7NBEJJ9RsiMiLmnkyJEYhsH48eMBc7LA9957jxdffJGjR486NzgRyVeU7IiIy1m3bh1Tpkzhq6++wsfHx7r/qaee4s4777yh11kiUvhYDP2NISIFSKVKlRgyZAhDhgzJlfaOHj1K5cqV2b59Ow0aNMiVNkUkb+nJjogUOCNGjMDX15e4uLgbaqd169bUrl07l6ISEWfRkx0RKVCOHTvG5cuXAahSpQpubjn/N91ff/3FxYsXAahQoUK2Z24WEdegZEdEREQKNL3GEhERkQJNyY6IiIgUaEp2REREpEBTsiMiIiIFmpIdERERKdCU7IiIiEiBpmRHRERECjQlOyIiIlKg/R9gBbWQG6TkkwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p1 = path_from_RT90(df, east_col = \"East\", north_col = \"North\", simplify_path = False)\n", "p2 = path_from_RT90(df, east_col = \"East\", north_col = \"North\", simplify_path = True)\n", "\n", "p1.plot(color='r', linestyle='-', linewidth = 10, label = \"Full path\")\n", "p2.plot(color='b', linestyle='-', marker = 'o', label = \"Reduced path\")\n", "\n", "plt.xlabel(\"X [m]\")\n", "plt.ylabel(\"Y [m]\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "da8620fa", "metadata": {}, "source": [ "## Making a landscape\n", "\n", "You can make a landscape with or without a path. Let's make an empty landscape." ] }, { "cell_type": "code", "execution_count": 4, "id": "24f1159d", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAG2CAYAAACEbnlbAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAK65JREFUeJzt3X90VPWd//HX5CckYSYGkxkihFCRH1FARA0jdV0lEjFaKHGXejgSXI4iDSgEKU1LQdhqWLr+YovSrSvRsyKWruiKBUlBwiqRH8EsEYSCi0mUTKJwMiEgSUju9w9P5tv5gDXRzEwSno9z5hzm8/nce9+ffJrm5Z1779gsy7IEAAAAn7BQFwAAANDVEJAAAAAMBCQAAAADAQkAAMBAQAIAADAQkAAAAAwEJAAAAAMBCQAAwEBAAgAAMBCQAAAADCENSI899phsNpvfa9iwYb7+c+fOKTc3V3379lVcXJyys7NVU1Pjt4/KykplZWUpJiZGSUlJWrhwoc6fPx/sqQAAgB4kItQFXH311frzn//sex8R8f9Lmj9/vt5++21t2LBBDodDc+bM0ZQpU/T+++9LklpaWpSVlSWXy6Vdu3apurpa06dPV2RkpJ544omgzwUAAPQMtlB+We1jjz2mN954Q2VlZRf0eb1eJSYmat26dbrnnnskSYcPH9bw4cNVUlKisWPHavPmzbrrrrt04sQJOZ1OSdKaNWu0aNEiffHFF4qKigrmdAAAQA8R8jNIR48eVXJysnr16iW3262CggKlpKSotLRUzc3NysjI8I0dNmyYUlJSfAGppKREI0aM8IUjScrMzNTs2bN18OBBjR49+qLHbGxsVGNjo+99a2urTp06pb59+8pmswVusgAAoNNYlqXTp08rOTlZYWGde9VQSANSenq6CgsLNXToUFVXV2vZsmW6+eab9dFHH8nj8SgqKkrx8fF+2zidTnk8HkmSx+PxC0dt/W1936SgoEDLli3r3MkAAICQqKqqUv/+/Tt1nyENSBMnTvT9e+TIkUpPT9fAgQP1hz/8Qb179w7YcfPz85WXl+d77/V6lZKSoqqqKtnt9oAdFwAAdJ76+noNGDBAffr06fR9h/wjtr8WHx+vIUOG6NixY7r99tvV1NSkuro6v7NINTU1crlckiSXy6U9e/b47aPtLre2MRcTHR2t6OjoC9rtdjsBCQCAbiYQl8d0qecgNTQ06JNPPlG/fv00ZswYRUZGatu2bb7+I0eOqLKyUm63W5LkdrtVXl6u2tpa35iioiLZ7XalpaUFvX4AANAzhPQM0qOPPqq7775bAwcO1IkTJ7R06VKFh4fr3nvvlcPh0MyZM5WXl6eEhATZ7XbNnTtXbrdbY8eOlSRNmDBBaWlpuu+++7Ry5Up5PB4tXrxYubm5Fz1DBAAA0B4hDUifffaZ7r33Xp08eVKJiYn64Q9/qA8++ECJiYmSpKefflphYWHKzs5WY2OjMjMz9dxzz/m2Dw8P16ZNmzR79my53W7FxsYqJydHy5cvD9WUAABADxDS5yB1FfX19XI4HPJ6vVyDBABANxHIv99d6hokAACAroCABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGLpMQFqxYoVsNpvmzZvnazt37pxyc3PVt29fxcXFKTs7WzU1NX7bVVZWKisrSzExMUpKStLChQt1/vz5IFcPAAB6ki4RkPbu3avf/e53GjlypF/7/Pnz9dZbb2nDhg0qLi7WiRMnNGXKFF9/S0uLsrKy1NTUpF27dumll15SYWGhlixZEuwpAACAHiTkAamhoUHTpk3T73//e1122WW+dq/Xq//4j//QU089pdtuu01jxozR2rVrtWvXLn3wwQeSpK1bt+rQoUP6z//8T1177bWaOHGi/vmf/1mrV69WU1NTqKYEAAC6uZAHpNzcXGVlZSkjI8OvvbS0VM3NzX7tw4YNU0pKikpKSiRJJSUlGjFihJxOp29MZmam6uvrdfDgwW88ZmNjo+rr6/1eAAAAbSJCefD169dr//792rt37wV9Ho9HUVFRio+P92t3Op3yeDy+MX8djtr62/q+SUFBgZYtW/Y9qwcAAD1VyM4gVVVV6ZFHHtErr7yiXr16BfXY+fn58nq9vldVVVVQjw8AALq2kAWk0tJS1dbW6rrrrlNERIQiIiJUXFysVatWKSIiQk6nU01NTaqrq/PbrqamRi6XS5LkcrkuuKut7X3bmIuJjo6W3W73ewEAALQJWUAaP368ysvLVVZW5ntdf/31mjZtmu/fkZGR2rZtm2+bI0eOqLKyUm63W5LkdrtVXl6u2tpa35iioiLZ7XalpaUFfU4AAKBnCNk1SH369NE111zj1xYbG6u+ffv62mfOnKm8vDwlJCTIbrdr7ty5crvdGjt2rCRpwoQJSktL03333aeVK1fK4/Fo8eLFys3NVXR0dNDnBAAAeoaQXqT9bZ5++mmFhYUpOztbjY2NyszM1HPPPefrDw8P16ZNmzR79my53W7FxsYqJydHy5cvD2HVAACgu7NZlmWFuohQq6+vl8PhkNfr5XokAAC6iUD+/Q75c5AAAAC6GgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgCGlAev755zVy5EjZ7XbZ7Xa53W5t3rzZ13/u3Dnl5uaqb9++iouLU3Z2tmpqavz2UVlZqaysLMXExCgpKUkLFy7U+fPngz0VAADQg4Q0IPXv318rVqxQaWmp9u3bp9tuu02TJk3SwYMHJUnz58/XW2+9pQ0bNqi4uFgnTpzQlClTfNu3tLQoKytLTU1N2rVrl1566SUVFhZqyZIloZoSAADoAWyWZVmhLuKvJSQk6De/+Y3uueceJSYmat26dbrnnnskSYcPH9bw4cNVUlKisWPHavPmzbrrrrt04sQJOZ1OSdKaNWu0aNEiffHFF4qKimrXMevr6+VwOOT1emW32wM2NwAA0HkC+fe7y1yD1NLSovXr1+vMmTNyu90qLS1Vc3OzMjIyfGOGDRumlJQUlZSUSJJKSko0YsQIXziSpMzMTNXX1/vOQl1MY2Oj6uvr/V4AAABtQh6QysvLFRcXp+joaD300EPauHGj0tLS5PF4FBUVpfj4eL/xTqdTHo9HkuTxePzCUVt/W983KSgokMPh8L0GDBjQuZMCAADdWsgD0tChQ1VWVqbdu3dr9uzZysnJ0aFDhwJ6zPz8fHm9Xt+rqqoqoMcDAADdS0SoC4iKitLgwYMlSWPGjNHevXv17LPPaurUqWpqalJdXZ3fWaSamhq5XC5Jksvl0p49e/z213aXW9uYi4mOjlZ0dHQnzwQAAPQU7QpIq1at6vCO77//fvXp06fD27W2tqqxsVFjxoxRZGSktm3bpuzsbEnSkSNHVFlZKbfbLUlyu916/PHHVVtbq6SkJElSUVGR7Ha70tLSOnxsAAAAqZ0Bad68eerfv7/Cw8PbtdOqqirddddd3xqQ8vPzNXHiRKWkpOj06dNat26dduzYoXfeeUcOh0MzZ85UXl6eEhISZLfbNXfuXLndbo0dO1aSNGHCBKWlpem+++7TypUr5fF4tHjxYuXm5nKGCAAAfGft/oht3759vrM036a9Z45qa2s1ffp0VVdXy+FwaOTIkXrnnXd0++23S5KefvpphYWFKTs7W42NjcrMzNRzzz3n2z48PFybNm3S7Nmz5Xa7FRsbq5ycHC1fvry90wIAALhAu56DtGzZMi1cuFAxMTHt2mlBQYFmz559wR1oXRXPQQIAoPsJ5N/vLvegyFAgIAEA0P10qQdFfvXVVzp79qzvfUVFhZ555hlt3bq1UwsDAAAIlQ4HpEmTJunll1+WJNXV1Sk9PV1PPvmkJk2apOeff77TCwQAAAi2Dgek/fv36+abb5Yk/fGPf5TT6VRFRYVefvnl7/Q4AAAAgK6mwwHp7NmzvrvUtm7dqilTpigsLExjx45VRUVFpxcIAAAQbB0OSIMHD9Ybb7yhqqoqvfPOO5owYYKkr2/Z5wJnAADQE3Q4IC1ZskSPPvqoUlNTlZ6e7nuq9datWzV69OhOLxAAACDYvtNt/h6PR9XV1Ro1apTCwr7OWHv27JHdbtewYcM6vchA4zZ/AAC6n0D+/e7Ql9U2Nzerd+/eKisru+Bs0Y033tiphQEAAIRKhz5ii4yMVEpKilpaWgJVDwAAQMh1+BqkX/7yl/rFL36hU6dOBaIeAACAkOvQR2yS9Nvf/lbHjh1TcnKyBg4cqNjYWL/+/fv3d1pxAAAAodDhgDR58uQAlAEAANB18GW14i42AAC6oy71ZbXS19/B9sILLyg/P993LdL+/fv1+eefd2pxAAAAodDhj9gOHDigjIwMORwOffrpp3rggQeUkJCg119/XZWVlb4vsgUAAOiuOnwGKS8vTzNmzNDRo0fVq1cvX/udd96pnTt3dmpxAAAAodDhgLR3717NmjXrgvYrrrhCHo+nU4oCAAAIpQ4HpOjoaNXX11/Q/pe//EWJiYmdUhQAAEAodTgg/ehHP9Ly5cvV3NwsSbLZbKqsrNSiRYuUnZ3d6QUCAAAEW4cD0pNPPqmGhgYlJSXpq6++0i233KLBgwerT58+evzxxwNRIwAAQFB1+C42h8OhoqIivf/++/rf//1fNTQ06LrrrlNGRkYg6gMAAAi6Dgekl19+WVOnTtW4ceM0btw4X3tTU5PWr1+v6dOnd2qBAAAAwdbhJ2mHh4erurpaSUlJfu0nT55UUlKSWlpaOrXAYOBJ2gAAdD9d6knalmXJZrNd0P7ZZ5/J4XB0SlEAAACh1O6P2EaPHi2bzSabzabx48crIuL/b9rS0qLjx4/rjjvuCEiRAAAAwdTugDR58mRJUllZmTIzMxUXF+fri4qKUmpqKrf5AwCAHqHdAWnp0qWSpNTUVE2dOtXva0YAAAB6kg5fg5STk6Nz587phRdeUH5+vk6dOiVJ2r9/vz7//PNOLxAAACDYOnyb/4EDB5SRkSGHw6FPP/1UDzzwgBISEvT666+rsrJSL7/8ciDqBAAACJoOn0GaP3++ZsyYoaNHj/p9zHbnnXdq586dnVocAABAKHT4DNK+ffv07//+7xe0X3HFFfJ4PJ1SFAAAQCh1+AxSdHS06uvrL2j/y1/+osTExE4pCgAAIJQ6HJB+9KMfafny5WpubpYk2Ww2VVZWatGiRdzmDwAAeoQOB6Qnn3xSDQ0NSkpK0ldffaVbbrlFgwcPVp8+ffT4448HokYAAICg6vA1SA6HQ0VFRXrvvfd04MABNTQ06LrrrlNGRkYg6gMAAAi6Dn9ZbU/El9UCAND9BPLvd4fPIEnS3r179e6776q2tlatra1+fU899VSnFAYAABAqHQ5ITzzxhBYvXqyhQ4fK6XTKZrP5+v763wAAAN1VhwPSs88+qxdffFEzZswIQDkAAACh1+G72MLCwjRu3LhA1AIAANAlfKevGlm9enUgagEAAOgSOvwR26OPPqqsrCxdeeWVSktLU2RkpF//66+/3mnFAQAAhEKHA9LDDz+sd999V7feeqv69u3LhdkAAKDH6XBAeumll/Rf//VfysrKCkQ9AAAAIdfha5ASEhJ05ZVXBqIWAACALqHDAemxxx7T0qVLdfbs2UDUAwAAEHId/oht1apV+uSTT+R0OpWamnrBRdr79+/vtOIAAABCocMBafLkyQEoAwAAoOvgy2rFl9UCANAdBfLvd4evQQIAAOjp2hWQEhIS9OWXX7Z7pykpKaqoqPjORQEAAIRSu65Bqqur0+bNm+VwONq105MnT6qlpeV7FQYAABAq7b5IOycnJ5B1AAAAdBntCkitra2BrgMAAKDL4CJtAAAAAwEJAADAQEACAAAwtDsgnThxIpB1AAAAdBntDkhXX3211q1bF8haAAAAuoR2B6THH39cs2bN0j/8wz/o1KlTgawJAAAgpNodkH7605/qwIEDOnnypNLS0vTWW28Fsi4AAICQafeDIiVp0KBB2r59u377299qypQpGj58uCIi/Hexf//+Ti0QAAAg2Dp8F1tFRYVef/11XXbZZZo0adIFr44oKCjQDTfcoD59+igpKUmTJ0/WkSNH/MacO3dOubm56tu3r+Li4pSdna2amhq/MZWVlcrKylJMTIySkpK0cOFCnT9/vqNTAwAAkNTBM0i///3vtWDBAmVkZOjgwYNKTEz8XgcvLi5Wbm6ubrjhBp0/f16/+MUvNGHCBB06dEixsbGSpPnz5+vtt9/Whg0b5HA4NGfOHE2ZMkXvv/++JKmlpUVZWVlyuVzatWuXqqurNX36dEVGRuqJJ574XvUBAIBLk82yLKs9A++44w7t2bNHzzzzjKZPnx6QYr744gslJSWpuLhYf/d3fyev16vExEStW7dO99xzjyTp8OHDGj58uEpKSjR27Fht3rxZd911l06cOCGn0ylJWrNmjRYtWqQvvvhCUVFR33rc+vp6ORwOeb1e2e32gMwNAAB0rkD+/W73R2wtLS06cOBAwMKRJHm9XklSQkKCJKm0tFTNzc3KyMjwjRk2bJhSUlJUUlIiSSopKdGIESN84UiSMjMzVV9fr4MHD170OI2Njaqvr/d7AQAAtGl3QCoqKlL//v0DVkhra6vmzZuncePG6ZprrpEkeTweRUVFKT4+3m+s0+mUx+PxjfnrcNTW39Z3MQUFBXI4HL7XgAEDOnk2AACgO+syXzWSm5urjz76SOvXrw/4sfLz8+X1en2vqqqqgB8TAAB0Hx26SDtQ5syZo02bNmnnzp1+Z6lcLpeamppUV1fndxappqZGLpfLN2bPnj1++2u7y61tjCk6OlrR0dGdPAsAANBThPQMkmVZmjNnjjZu3Kjt27dr0KBBfv1jxoxRZGSktm3b5ms7cuSIKisr5Xa7JUlut1vl5eWqra31jSkqKpLdbldaWlpwJgIAAHqUkJ5Bys3N1bp16/Tmm2+qT58+vmuGHA6HevfuLYfDoZkzZyovL08JCQmy2+2aO3eu3G63xo4dK0maMGGC0tLSdN9992nlypXyeDxavHixcnNzOUsEAAC+k3bf5h+Qg9tsF21fu3atZsyYIenrB0UuWLBAr776qhobG5WZmannnnvO7+OziooKzZ49Wzt27FBsbKxycnK0YsWKC57y/U24zR8AgO4nkH+/QxqQugoCEgAA3U+XeA4SAADApYKABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGEIakHbu3Km7775bycnJstlseuONN/z6LcvSkiVL1K9fP/Xu3VsZGRk6evSo35hTp05p2rRpstvtio+P18yZM9XQ0BDEWQAAgJ4mpAHpzJkzGjVqlFavXn3R/pUrV2rVqlVas2aNdu/erdjYWGVmZurcuXO+MdOmTdPBgwdVVFSkTZs2aefOnXrwwQeDNQUAANAD2SzLskJdhCTZbDZt3LhRkydPlvT12aPk5GQtWLBAjz76qCTJ6/XK6XSqsLBQP/nJT/Txxx8rLS1Ne/fu1fXXXy9J2rJli+6880599tlnSk5Obtex6+vr5XA45PV6ZbfbAzI/AADQuQL597vLXoN0/PhxeTweZWRk+NocDofS09NVUlIiSSopKVF8fLwvHElSRkaGwsLCtHv37m/cd2Njo+rr6/1eAAAAbbpsQPJ4PJIkp9Pp1+50On19Ho9HSUlJfv0RERFKSEjwjbmYgoICORwO32vAgAGdXD0AAOjOumxACqT8/Hx5vV7fq6qqKtQlAQCALqTLBiSXyyVJqqmp8Wuvqanx9blcLtXW1vr1nz9/XqdOnfKNuZjo6GjZ7Xa/FwAAQJsuG5AGDRokl8ulbdu2+drq6+u1e/duud1uSZLb7VZdXZ1KS0t9Y7Zv367W1lalp6cHvWYAANAzRITy4A0NDTp27Jjv/fHjx1VWVqaEhASlpKRo3rx5+vWvf62rrrpKgwYN0q9+9SslJyf77nQbPny47rjjDj3wwANas2aNmpubNWfOHP3kJz9p9x1sAAAAppAGpH379unWW2/1vc/Ly5Mk5eTkqLCwUD/72c905swZPfjgg6qrq9MPf/hDbdmyRb169fJt88orr2jOnDkaP368wsLClJ2drVWrVgV9LgAAoOfoMs9BCiWegwQAQPdzST4HCQAAIFQISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAQAAGAhIAAAABgISAACAgYAEAABgICABAAAYCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgKHHBKTVq1crNTVVvXr1Unp6uvbs2RPqkgAAQDfVIwLSa6+9pry8PC1dulT79+/XqFGjlJmZqdra2lCXBgAAuqEeEZCeeuopPfDAA7r//vuVlpamNWvWKCYmRi+++GKoSwMAAN1QRKgL+L6amppUWlqq/Px8X1tYWJgyMjJUUlJy0W0aGxvV2Njoe+/1eiVJ9fX1gS0WAAB0mra/25Zldfq+u31A+vLLL9XS0iKn0+nX7nQ6dfjw4YtuU1BQoGXLll3QPmDAgIDUCAAAAufkyZNyOBydus9uH5C+i/z8fOXl5fne19XVaeDAgaqsrOz0HzA6pr6+XgMGDFBVVZXsdnuoy7mksRZdB2vRdbAWXYvX61VKSooSEhI6fd/dPiBdfvnlCg8PV01NjV97TU2NXC7XRbeJjo5WdHT0Be0Oh4P/wXcRdrudtegiWIuug7XoOliLriUsrPMvqe72F2lHRUVpzJgx2rZtm6+ttbVV27Ztk9vtDmFlAACgu+r2Z5AkKS8vTzk5Obr++ut144036plnntGZM2d0//33h7o0AADQDfWIgDR16lR98cUXWrJkiTwej6699lpt2bLlggu3v0l0dLSWLl160Y/dEFysRdfBWnQdrEXXwVp0LYFcD5sViHvjAAAAurFufw0SAABAZyMgAQAAGAhIAAAABgISAACA4ZIPSKtXr1Zqaqp69eql9PR07dmzJ9Ql9Tg7d+7U3XffreTkZNlsNr3xxht+/ZZlacmSJerXr5969+6tjIwMHT161G/MqVOnNG3aNNntdsXHx2vmzJlqaGgI4ix6hoKCAt1www3q06ePkpKSNHnyZB05csRvzLlz55Sbm6u+ffsqLi5O2dnZFzyItbKyUllZWYqJiVFSUpIWLlyo8+fPB3Mq3d7zzz+vkSNH+h446Ha7tXnzZl8/6xA6K1askM1m07x583xtrEfwPPbYY7LZbH6vYcOG+fqDtRaXdEB67bXXlJeXp6VLl2r//v0aNWqUMjMzVVtbG+rSepQzZ85o1KhRWr169UX7V65cqVWrVmnNmjXavXu3YmNjlZmZqXPnzvnGTJs2TQcPHlRRUZE2bdqknTt36sEHHwzWFHqM4uJi5ebm6oMPPlBRUZGam5s1YcIEnTlzxjdm/vz5euutt7RhwwYVFxfrxIkTmjJliq+/paVFWVlZampq0q5du/TSSy+psLBQS5YsCcWUuq3+/ftrxYoVKi0t1b59+3Tbbbdp0qRJOnjwoCTWIVT27t2r3/3udxo5cqRfO+sRXFdffbWqq6t9r/fee8/XF7S1sC5hN954o5Wbm+t739LSYiUnJ1sFBQUhrKpnk2Rt3LjR9761tdVyuVzWb37zG19bXV2dFR0dbb366quWZVnWoUOHLEnW3r17fWM2b95s2Ww26/PPPw9a7T1RbW2tJckqLi62LOvrn31kZKS1YcMG35iPP/7YkmSVlJRYlmVZf/rTn6ywsDDL4/H4xjz//POW3W63GhsbgzuBHuayyy6zXnjhBdYhRE6fPm1dddVVVlFRkXXLLbdYjzzyiGVZ/F4E29KlS61Ro0ZdtC+Ya3HJnkFqampSaWmpMjIyfG1hYWHKyMhQSUlJCCu7tBw/flwej8dvHRwOh9LT033rUFJSovj4eF1//fW+MRkZGQoLC9Pu3buDXnNP4vV6Jcn3RY+lpaVqbm72W49hw4YpJSXFbz1GjBjh9yDWzMxM1dfX+85+oGNaWlq0fv16nTlzRm63m3UIkdzcXGVlZfn93CV+L0Lh6NGjSk5O1g9+8ANNmzZNlZWVkoK7Fj3iSdrfxZdffqmWlpYLnrbtdDp1+PDhEFV16fF4PJJ00XVo6/N4PEpKSvLrj4iIUEJCgm8MOq61tVXz5s3TuHHjdM0110j6+mcdFRWl+Ph4v7Hmelxsvdr60H7l5eVyu906d+6c4uLitHHjRqWlpamsrIx1CLL169dr//792rt37wV9/F4EV3p6ugoLCzV06FBVV1dr2bJluvnmm/XRRx8FdS0u2YAEXOpyc3P10Ucf+X22j+AaOnSoysrK5PV69cc//lE5OTkqLi4OdVmXnKqqKj3yyCMqKipSr169Ql3OJW/ixIm+f48cOVLp6ekaOHCg/vCHP6h3795Bq+OS/Yjt8ssvV3h4+AVXvtfU1MjlcoWoqktP28/6b62Dy+W64ML58+fP69SpU6zVdzRnzhxt2rRJ7777rvr37+9rd7lcampqUl1dnd94cz0utl5tfWi/qKgoDR48WGPGjFFBQYFGjRqlZ599lnUIstLSUtXW1uq6665TRESEIiIiVFxcrFWrVikiIkJOp5P1CKH4+HgNGTJEx44dC+rvxiUbkKKiojRmzBht27bN19ba2qpt27bJ7XaHsLJLy6BBg+RyufzWob6+Xrt37/atg9vtVl1dnUpLS31jtm/frtbWVqWnpwe95u7MsizNmTNHGzdu1Pbt2zVo0CC//jFjxigyMtJvPY4cOaLKykq/9SgvL/cLrUVFRbLb7UpLSwvORHqo1tZWNTY2sg5BNn78eJWXl6usrMz3uv766zVt2jTfv1mP0GloaNAnn3yifv36Bfd34ztdYt5DrF+/3oqOjrYKCwutQ4cOWQ8++KAVHx/vd+U7vr/Tp09bH374ofXhhx9akqynnnrK+vDDD62KigrLsixrxYoVVnx8vPXmm29aBw4csCZNmmQNGjTI+uqrr3z7uOOOO6zRo0dbu3fvtt577z3rqquusu69995QTanbmj17tuVwOKwdO3ZY1dXVvtfZs2d9Yx566CErJSXF2r59u7Vv3z7L7XZbbrfb13/+/HnrmmuusSZMmGCVlZVZW7ZssRITE638/PxQTKnb+vnPf24VFxdbx48ftw4cOGD9/Oc/t2w2m7V161bLsliHUPvru9gsi/UIpgULFlg7duywjh8/br3//vtWRkaGdfnll1u1tbWWZQVvLS7pgGRZlvVv//ZvVkpKihUVFWXdeOON1gcffBDqknqcd99915J0wSsnJ8eyrK9v9f/Vr35lOZ1OKzo62ho/frx15MgRv32cPHnSuvfee624uDjLbrdb999/v3X69OkQzKZ7u9g6SLLWrl3rG/PVV19ZP/3pT63LLrvMiomJsX784x9b1dXVfvv59NNPrYkTJ1q9e/e2Lr/8cmvBggVWc3NzkGfTvf3TP/2TNXDgQCsqKspKTEy0xo8f7wtHlsU6hJoZkFiP4Jk6darVr18/KyoqyrriiiusqVOnWseOHfP1B2stbJZlWd/r3BcAAEAPc8legwQAAPBNCEgAAAAGAhIAAICBgAQAAGAgIAEAABgISAAAAAYCEgAAgIGABAAAYCAgAegyUlNTZbPZZLPZLvgyyu6ssLDQN6958+aFuhwA7UBAAtCpWlpadNNNN2nKlCl+7V6vVwMGDNAvf/nLv7n98uXLVV1dLYfDEcgyVVhYqPj4+IAeo83UqVNVXV3NF2ED3QgBCUCnCg8PV2FhobZs2aJXXnnF1z537lwlJCRo6dKlf3P7Pn36yOVyyWazBbrUTtHS0qLW1ta/OaZ3795yuVyKiooKUlUAvi8CEoBON2TIEK1YsUJz585VdXW13nzzTa1fv14vv/xyh0NC25meTZs2aejQoYqJidE999yjs2fP6qWXXlJqaqouu+wyPfzww2ppafFt19jYqEcffVRXXHGFYmNjlZ6erh07dkiSduzYofvvv19er9f30ddjjz32rdv9dT3//d//rbS0NEVHR6uyslI7duzQjTfeqNjYWMXHx2vcuHGqqKj4vj9KACESEeoCAPRMc+fO1caNG3XfffepvLxcS5Ys0ahRo77Tvs6ePatVq1Zp/fr1On36tKZMmaIf//jHio+P15/+9Cf93//9n7KzszVu3DhNnTpVkjRnzhwdOnRI69evV3JysjZu3Kg77rhD5eXluummm/TMM89oyZIlOnLkiCQpLi7uW7e76qqrfPX8y7/8i1544QX17dtXCQkJuvbaa/XAAw/o1VdfVVNTk/bs2dNtzoIBuAgLAALk448/tiRZI0aMsJqbm791/MCBA62nn37ar23t2rWWJOvYsWO+tlmzZlkxMTHW6dOnfW2ZmZnWrFmzLMuyrIqKCis8PNz6/PPP/fY1fvx4Kz8/37dfh8Ph19/e7SRZZWVlvv6TJ09akqwdO3b8zfndcsst1iOPPPI3xwDoGjiDBCBgXnzxRcXExOj48eP67LPPlJqa+p32ExMToyuvvNL33ul0KjU11XfWp62ttrZWklReXq6WlhYNGTLEbz+NjY3q27fvNx6nvdtFRUVp5MiRvvcJCQmaMWOGMjMzdfvttysjI0P/+I//qH79+n2n+QIIPQISgIDYtWuXnn76aW3dulW//vWvNXPmTP35z3/+Th87RUZG+r232WwXbWu7WLqhoUHh4eEqLS1VeHi437i/DlWm9m7Xu3fvC+axdu1aPfzww9qyZYtee+01LV68WEVFRRo7dmz7JwqgyyAgAeh0Z8+e1YwZMzR79mzdeuutGjRokEaMGKE1a9Zo9uzZAT/+6NGj1dLSotraWt18880XHRMVFeV3UXd7t/u2444ePVr5+flyu91at24dAQnopriLDUCny8/Pl2VZWrFihaSvHwD5r//6r/rZz36mTz/9NODHHzJkiKZNm6bp06fr9ddf1/Hjx7Vnzx4VFBTo7bff9tXU0NCgbdu26csvv9TZs2fbtd3FHD9+XPn5+SopKVFFRYW2bt2qo0ePavjw4QGfK4DAICAB6FTFxcVavXq11q5dq5iYGF/7rFmzdNNNN2nmzJmyLCvgdaxdu1bTp0/XggULNHToUE2ePFl79+5VSkqKJOmmm27SQw89pKlTpyoxMVErV65s13YXExMTo8OHDys7O1tDhgzRgw8+qNzcXM2aNSvg8wQQGDYrGP9PBQDtkJqaqnnz5vXYr+P4+7//e1177bV65plnQl0KgG/BGSQAXcqiRYsUFxcnr9cb6lI6zSuvvKK4uDj9z//8T6hLAdBOnEEC0GVUVFSoublZkvSDH/xAYWE947/hTp8+rZqaGklSfHy8Lr/88hBXBODbEJAAAAAMPeM/zwAAADoRAQkAAMBAQAIAADAQkAAAAAwEJAAAAAMBCQAAwEBAAgAAMBCQAAAADP8PsN1mgrT5m78AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "landscape = Landscape()\n", "fig, ax = landscape.plot()\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "035b4f42", "metadata": {}, "source": [ "## Adding a source" ] }, { "cell_type": "code", "execution_count": 5, "id": "91019da5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[Source(name=Source 1, pos=(100, 100, 0), isotope=Cs137, A=100 MBq)]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cs137 = Isotope(name = \"Cs137\", E = 662.66, b = 0.851)\n", "my_source = Source(x = 100, y = 100, z = 0, activity = 100, isotope = cs137)\n", "landscape.add_sources(my_source)\n", "\n", "landscape.sources" ] }, { "cell_type": "code", "execution_count": 6, "id": "7913fe1e", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAG2CAYAAACEbnlbAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAL/xJREFUeJzt3Xt4VNW9//HP5Aq5zMRgLiAEqMolCoioYaSWKpGAoUKJLfXwSKAcRRpQLlKaFkE41VB6vNGieOqR4FMRS49oxYKkILFC5BJMiSAULCYomUThlwkBSUKyfn9MMzob1EQzmSS8X8+zH5m91t77u7PEfFz7MjZjjBEAAAC8ggJdAAAAQFtDQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsAhoQHrooYdks9l8ln79+nnbz549q6ysLHXp0kVRUVHKyMhQeXm5zz5KS0uVnp6uiIgIxcfHa968eTp37lxrnwoAAOhAQgJdwFVXXaW//e1v3s8hIZ+XNHv2bL3++utat26dHA6HZsyYofHjx2v79u2SpPr6eqWnpysxMVE7duxQWVmZJk2apNDQUD3yyCOtfi4AAKBjsAXyy2ofeughvfLKKyoqKjqvze12Ky4uTmvWrNEdd9whSTp48KD69++vgoICDR06VBs3btSYMWN0/PhxJSQkSJJWrlyp+fPn65NPPlFYWFhrng4AAOggAj6DdPjwYXXr1k2dOnWS0+lUTk6OkpKSVFhYqLq6OqWmpnr79uvXT0lJSd6AVFBQoAEDBnjDkSSlpaVp+vTp2r9/vwYPHnzBY9bU1Kimpsb7uaGhQSdPnlSXLl1ks9n8d7IAAKDFGGN06tQpdevWTUFBLXvXUEADUkpKinJzc9W3b1+VlZVp8eLFuummm/Tee+/J5XIpLCxMMTExPtskJCTI5XJJklwul084amxvbPsyOTk5Wrx4ccueDAAACIhjx46pe/fuLbrPgAak0aNHe/88cOBApaSkqGfPnvrTn/6kzp07++242dnZmjNnjvez2+1WUlKSjh07Jrvd7rfjAgCAllNVVaUePXooOjq6xfcd8EtsXxQTE6M+ffroyJEjuvXWW1VbW6vKykqfWaTy8nIlJiZKkhITE7Vr1y6ffTQ+5dbY50LCw8MVHh5+3nq73U5AAgCgnfHH7TFt6j1I1dXV+uCDD9S1a1cNGTJEoaGh2rJli7f90KFDKi0tldPplCQ5nU4VFxeroqLC2ycvL092u13JycmtXj8AAOgYAjqD9MADD+gHP/iBevbsqePHj2vRokUKDg7WnXfeKYfDoalTp2rOnDmKjY2V3W7XzJkz5XQ6NXToUEnSyJEjlZycrLvuukvLli2Ty+XSggULlJWVdcEZIgAAgKYIaED66KOPdOedd+rEiROKi4vTd7/7Xb3zzjuKi4uTJD3++OMKCgpSRkaGampqlJaWpqeeesq7fXBwsDZs2KDp06fL6XQqMjJSmZmZWrJkSaBOCQAAdAABfQ9SW1FVVSWHwyG32809SAAAtBP+/P3dpu5BAgAAaAsISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwaDMBaenSpbLZbJo1a5Z33dmzZ5WVlaUuXbooKipKGRkZKi8v99mutLRU6enpioiIUHx8vObNm6dz5861cvUAAKAjaRMBaffu3XrmmWc0cOBAn/WzZ8/Wa6+9pnXr1ik/P1/Hjx/X+PHjve319fVKT09XbW2tduzYodWrVys3N1cLFy5s7VMAAAAdSMADUnV1tSZOnKg//OEPuuSSS7zr3W63/vd//1ePPfaYbrnlFg0ZMkSrVq3Sjh079M4770iSNm/erAMHDuiPf/yjrrnmGo0ePVr/9V//pRUrVqi2tjZQpwQAANq5gAekrKwspaenKzU11Wd9YWGh6urqfNb369dPSUlJKigokCQVFBRowIABSkhI8PZJS0tTVVWV9u/f/6XHrKmpUVVVlc8CAADQKCSQB1+7dq327t2r3bt3n9fmcrkUFhammJgYn/UJCQlyuVzePl8MR43tjW1fJicnR4sXL/6W1QMAgI4qYDNIx44d0/33368XXnhBnTp1atVjZ2dny+12e5djx4616vEBAEDbFrCAVFhYqIqKCl177bUKCQlRSEiI8vPztXz5coWEhCghIUG1tbWqrKz02a68vFyJiYmSpMTExPOeamv83NjnQsLDw2W3230WAACARgELSCNGjFBxcbGKioq8y3XXXaeJEyd6/xwaGqotW7Z4tzl06JBKS0vldDolSU6nU8XFxaqoqPD2ycvLk91uV3JycqufEwAA6BgCdg9SdHS0rr76ap91kZGR6tKli3f91KlTNWfOHMXGxsput2vmzJlyOp0aOnSoJGnkyJFKTk7WXXfdpWXLlsnlcmnBggXKyspSeHh4q58TAADoGAJ6k/bXefzxxxUUFKSMjAzV1NQoLS1NTz31lLc9ODhYGzZs0PTp0+V0OhUZGanMzEwtWbIkgFUDAID2zmaMMYEuItCqqqrkcDjkdru5HwkAgHbCn7+/A/4eJAAAgLaGgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAAi4AGpKeffloDBw6U3W6X3W6X0+nUxo0bve1nz55VVlaWunTpoqioKGVkZKi8vNxnH6WlpUpPT1dERITi4+M1b948nTt3rrVPBQAAdCABDUjdu3fX0qVLVVhYqD179uiWW27R2LFjtX//fknS7Nmz9dprr2ndunXKz8/X8ePHNX78eO/29fX1Sk9PV21trXbs2KHVq1crNzdXCxcuDNQpAQCADsBmjDGBLuKLYmNj9dvf/lZ33HGH4uLitGbNGt1xxx2SpIMHD6p///4qKCjQ0KFDtXHjRo0ZM0bHjx9XQkKCJGnlypWaP3++PvnkE4WFhTXpmFVVVXI4HHK73bLb7X47NwAA0HL8+fu7zdyDVF9fr7Vr1+r06dNyOp0qLCxUXV2dUlNTvX369eunpKQkFRQUSJIKCgo0YMAAbziSpLS0NFVVVXlnoS6kpqZGVVVVPgsAAECjgAek4uJiRUVFKTw8XPfee6/Wr1+v5ORkuVwuhYWFKSYmxqd/QkKCXC6XJMnlcvmEo8b2xrYvk5OTI4fD4V169OjRsicFAADatYAHpL59+6qoqEg7d+7U9OnTlZmZqQMHDvj1mNnZ2XK73d7l2LFjfj0eAABoX0ICXUBYWJiuuOIKSdKQIUO0e/duPfnkk5owYYJqa2tVWVnpM4tUXl6uxMRESVJiYqJ27drls7/Gp9wa+1xIeHi4wsPDW/hMAABAR9GkgLR8+fJm73jKlCmKjo5u9nYNDQ2qqanRkCFDFBoaqi1btigjI0OSdOjQIZWWlsrpdEqSnE6nHn74YVVUVCg+Pl6SlJeXJ7vdruTk5GYfGwAAQGpiQJo1a5a6d++u4ODgJu302LFjGjNmzNcGpOzsbI0ePVpJSUk6deqU1qxZo23btumNN96Qw+HQ1KlTNWfOHMXGxsput2vmzJlyOp0aOnSoJGnkyJFKTk7WXXfdpWXLlsnlcmnBggXKyspihggAAHxjTb7EtmfPHu8szddp6sxRRUWFJk2apLKyMjkcDg0cOFBvvPGGbr31VknS448/rqCgIGVkZKimpkZpaWl66qmnvNsHBwdrw4YNmj59upxOpyIjI5WZmaklS5Y09bQAAADO06T3IC1evFjz5s1TREREk3aak5Oj6dOnn/cEWlvFe5AAAGh//Pn7u829KDIQCEgAALQ/bepFkZ999pnOnDnj/VxSUqInnnhCmzdvbtHCAAAAAqXZAWns2LF6/vnnJUmVlZVKSUnRo48+qrFjx+rpp59u8QIBAABaW7MD0t69e3XTTTdJkv785z8rISFBJSUlev7557/R6wAAAADammYHpDNnznifUtu8ebPGjx+voKAgDR06VCUlJS1eIAAAQGtrdkC64oor9Morr+jYsWN64403NHLkSEmeR/a5wRkAAHQEzQ5ICxcu1AMPPKBevXopJSXF+1brzZs3a/DgwS1eIAAAQGv7Ro/5u1wulZWVadCgQQoK8mSsXbt2yW63q1+/fi1epL/xmD8AAO2PP39/N+vLauvq6tS5c2cVFRWdN1t0ww03tGhhAAAAgdKsS2yhoaFKSkpSfX29v+oBAAAIuGbfg/SrX/1Kv/zlL3Xy5El/1AMAABBwzbrEJkm///3vdeTIEXXr1k09e/ZUZGSkT/vevXtbrDgAAIBAaHZAGjdunB/KAAAAaDv4slrxFBsAAO1Rm/qyWsnzHWzPPvussrOzvfci7d27Vx9//HGLFgcAABAIzb7Etm/fPqWmpsrhcOjDDz/U3XffrdjYWL388ssqLS31fpEtAABAe9XsGaQ5c+Zo8uTJOnz4sDp16uRdf9ttt+mtt95q0eIAAAACodkBaffu3Zo2bdp56y+77DK5XK4WKQoAACCQmh2QwsPDVVVVdd76f/7zn4qLi2uRogAAAAKp2QHp9ttv15IlS1RXVydJstlsKi0t1fz585WRkdHiBQIAALS2ZgekRx99VNXV1YqPj9dnn32m4cOH64orrlB0dLQefvhhf9QIAADQqpr9FJvD4VBeXp62b9+uf/zjH6qurta1116r1NRUf9QHAADQ6podkJ5//nlNmDBBw4YN07Bhw7zra2trtXbtWk2aNKlFCwQAAGhtzX6TdnBwsMrKyhQfH++z/sSJE4qPj1d9fX2LFtgaeJM2AADtT5t6k7YxRjab7bz1H330kRwOR4sUBQAAEEhNvsQ2ePBg2Ww22Ww2jRgxQiEhn29aX1+vo0ePatSoUX4pEgAAoDU1OSCNGzdOklRUVKS0tDRFRUV528LCwtSrVy8e8wcAAB1CkwPSokWLJEm9evXShAkTfL5mBAAAoCNp9j1ImZmZOnv2rJ599lllZ2fr5MmTkqS9e/fq448/bvECAQAAWluzH/Pft2+fUlNT5XA49OGHH+ruu+9WbGysXn75ZZWWlur555/3R50AAACtptkzSLNnz9bkyZN1+PBhn8tst912m956660WLQ4AACAQmj2DtGfPHv3P//zPeesvu+wyuVyuFikKAAAgkJo9gxQeHq6qqqrz1v/zn/9UXFxcixQFAAAQSM0OSLfffruWLFmiuro6SZLNZlNpaanmz5/PY/4AAKBDaHZAevTRR1VdXa34+Hh99tlnGj58uK644gpFR0fr4Ycf9keNAAAArarZ9yA5HA7l5eXp7bff1r59+1RdXa1rr71Wqamp/qgPAACg1TX7y2o7Ir6sFgCA9sefv7+bPYMkSbt379abb76piooKNTQ0+LQ99thjLVIYAABAoDQ7ID3yyCNasGCB+vbtq4SEBNlsNm/bF/8MAADQXjU7ID355JN67rnnNHnyZD+UAwAAEHjNfootKChIw4YN80ctAAAAbcI3+qqRFStW+KMWAACANqHZl9geeOABpaen6/LLL1dycrJCQ0N92l9++eUWKw4AACAQmh2Q7rvvPr355pu6+eab1aVLF27MBgAAHU6zA9Lq1av1f//3f0pPT/dHPQAAAAHX7HuQYmNjdfnll/ujFgAAgDah2QHpoYce0qJFi3TmzBl/1AMAABBwzb7Etnz5cn3wwQdKSEhQr169zrtJe+/evS1WHAAAQCA0OyCNGzfOD2UAAAC0HXxZrfiyWgAA2iN//v5u9j1IAAAAHV2TAlJsbKw+/fTTJu80KSlJJSUl37goAACAQGrSPUiVlZXauHGjHA5Hk3Z64sQJ1dfXf6vCAAAAAqXJN2lnZmb6sw4AAIA2o0kBqaGhwd91AAAAtBncpA0AAGBBQAIAALAgIAEAAFg0OSAdP37cn3UAAAC0GU0OSFdddZXWrFnjz1oAAADahCYHpIcffljTpk3Tj370I508edKfNQEAAARUkwPSz372M+3bt08nTpxQcnKyXnvtNX/WBQAAEDBNflGkJPXu3Vtbt27V73//e40fP179+/dXSIjvLvbu3duiBQIAALS2Zj/FVlJSopdfflmXXHKJxo4de97SHDk5Obr++usVHR2t+Ph4jRs3TocOHfLpc/bsWWVlZalLly6KiopSRkaGysvLffqUlpYqPT1dERERio+P17x583Tu3LnmnhoAAICkZs4g/eEPf9DcuXOVmpqq/fv3Ky4u7lsdPD8/X1lZWbr++ut17tw5/fKXv9TIkSN14MABRUZGSpJmz56t119/XevWrZPD4dCMGTM0fvx4bd++XZJUX1+v9PR0JSYmaseOHSorK9OkSZMUGhqqRx555FvVBwAALk42Y4xpSsdRo0Zp165deuKJJzRp0iS/FPPJJ58oPj5e+fn5+t73vie32624uDitWbNGd9xxhyTp4MGD6t+/vwoKCjR06FBt3LhRY8aM0fHjx5WQkCBJWrlypebPn69PPvlEYWFhX3vcqqoqORwOud1u2e12v5wbAABoWf78/d3kS2z19fXat2+f38KRJLndbklSbGysJKmwsFB1dXVKTU319unXr5+SkpJUUFAgSSooKNCAAQO84UiS0tLSVFVVpf3791/wODU1NaqqqvJZAAAAGjU5IOXl5al79+5+K6ShoUGzZs3SsGHDdPXVV0uSXC6XwsLCFBMT49M3ISFBLpfL2+eL4aixvbHtQnJycuRwOLxLjx49WvhsAABAe9ZmvmokKytL7733ntauXev3Y2VnZ8vtdnuXY8eO+f2YAACg/WjWTdr+MmPGDG3YsEFvvfWWzyxVYmKiamtrVVlZ6TOLVF5ersTERG+fXbt2+eyv8Sm3xj5W4eHhCg8Pb+GzAAAAHUVAZ5CMMZoxY4bWr1+vrVu3qnfv3j7tQ4YMUWhoqLZs2eJdd+jQIZWWlsrpdEqSnE6niouLVVFR4e2Tl5cnu92u5OTk1jkRAADQoQR0BikrK0tr1qzRq6++qujoaO89Qw6HQ507d5bD4dDUqVM1Z84cxcbGym63a+bMmXI6nRo6dKgkaeTIkUpOTtZdd92lZcuWyeVyacGCBcrKymKWCAAAfCNNfszfLwe32S64ftWqVZo8ebIkz4si586dqxdffFE1NTVKS0vTU0895XP5rKSkRNOnT9e2bdsUGRmpzMxMLV269Ly3fH8ZHvMHAKD98efv74AGpLaCgAQAQPvTJt6DBAAAcLEgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCW3TJ59I06dLSUlSeLiUmCilpUnbtwe6sm/mrbekH/xA6tZNstmkV14JdEUAgK8QEugCgAvKyJBqa6XVq6XvfEcqL5e2bJFOnPDvcWtrpbCwlt/v6dPSoEHST38qjR/f8vsHALQoZpDQ9lRWSn//u/Sb30g33yz17CndcIOUnS3dfvvn/UpLpbFjpagoyW6XfvxjT5BqNHmyNG6c775nzZK+//3PP3//+9KMGZ71l17qmaWSpP37pTFjPPuNjpZuukn64IPPt3v2Wal/f6lTJ6lfP+mpp776nEaPln79a+mHP2zmDwMAEAgEJLQ9UVGe5ZVXpJqaC/dpaPCEo5Mnpfx8KS9P+te/pAkTmn+81as9s0bbt0srV0offyx973ueS3tbt0qFhZ6Zn3PnPP1feEFauFB6+GHp/felRx6RHnzQsx8AQIfAJTa0PSEhUm6udPfdnsBy7bXS8OHST34iDRzo6bNli1RcLB09KvXo4Vn3/PPSVVdJu3dL11/f9ONdeaW0bNnnn3/5S8nhkNaulUJDPev69Pm8fdEi6dFHP79U1ru3dOCA9MwzUmbmNz5tAEDbwQwS2qaMDOn4cekvf5FGjZK2bfMEpdxcT/v773uCUWM4kqTkZCkmxtPWHEOG+H4uKvJcUmsMR190+rTnUtvUqZ/PdEVFeS6fffESHACgXWMGCW1Xp07Srbd6lgcflP7zPz2zN5MnN237oCDJGN91dXXn94uM9P3cufOX77O62vPPP/xBSknxbQsOblpdAIA2jxkktB/JyZ4ZHMlzg/SxY56l0YEDnhu8k5M9n+PipLIy330UFX39cQYO9NwkfqEwlZDgeVT/X/+SrrjCd+nd+5ucFQCgDWIGCW3DoUPSa69JJSWeR+3//ndp5kxp2DDPU2R79njuExo71tM/NVUaMECaOFF64gnPDdQ/+5nnXqXrrvP0ueUW6be/9dyb5HRKf/yj9N570uDBX13LjBnS737nuecpO9tzP9I773iepOvbV1q8WLrvPs/6UaM8N5Lv2SP9v/8nzZlz4X1WV0tHjnz++ehRT1iLjfW86wkA0LaYAMrPzzdjxowxXbt2NZLM+vXrfdobGhrMgw8+aBITE02nTp3MiBEjzD//+U+fPidOnDD/8R//YaKjo43D4TA//elPzalTp5pVh9vtNpKM2+3+tqeE5jp2zJjhw42RjAkKMiY01JiQEM9nyZjwcGMiIozp29eYBQuMOXPm821LSoy5/XZjIiONiY425kc/Msbl8t3/woXGJCQY43AYM3u2MTNmeI7XaPhwY+6///y6/vEPY0aO9Bw7OtqYm24y5oMPPm9/4QVjrrnGmLAwYy65xJjvfc+Yl1/+8vN8883Pz+mLS2ZmM39gAIBG/vz9bTPGepNG69m4caO2b9+uIUOGaPz48Vq/fr3GfeG9Nb/5zW+Uk5Oj1atXq3fv3nrwwQdVXFysAwcOqFOnTpKk0aNHq6ysTM8884zq6uo0ZcoUXX/99VqzZk2T66iqqpLD4ZDb7Zbdbm/p08SXqaryzMh8+unnj9Bb2WyemZxmjCcA4OLgz9/fAQ1IX2Sz2XwCkjFG3bp109y5c/XAAw9IktxutxISEpSbm6uf/OQnev/995WcnKzdu3frun9fVtm0aZNuu+02ffTRR+rWrVuTjk1ACpAFC6ScHM87jb7Ojh2ey2QAAPybP39/t9mbtI8ePSqXy6XU1FTvOofDoZSUFBUUFEiSCgoKFBMT4w1HkpSamqqgoCDt3LnzS/ddU1OjqqoqnwUB8MYbTQtHISGeFzYCANBK2mxAcrlckqSEhASf9QkJCd42l8ul+Ph4n/aQkBDFxsZ6+1xITk6OHA6Hd+nxxXfpoPXY7Z5H8b9OQ4PnXUMAALSSNhuQ/Ck7O1tut9u7HPvio+JoPbNmNW0GKTJSuusuv5cDAECjNhuQEhMTJUnlX/zy0X9/bmxLTExURUWFT/u5c+d08uRJb58LCQ8Pl91u91kQAGPGeL5ORPLcjG0VHOy5vPbcc57H4QEAaCVtNiD17t1biYmJ2rJli3ddVVWVdu7cKee/b9Z1Op2qrKxUYWGht8/WrVvV0NCgFOtbjtH22Gye7y9bu1a68UbfkBQR4fmus717pTvuCFyNAICLUkBfFFldXa0jX3h53tGjR1VUVKTY2FglJSVp1qxZ+vWvf60rr7zS+5h/t27dvE+69e/fX6NGjdLdd9+tlStXqq6uTjNmzNBPfvKTJj/BhgCz2aQJEzzLyZPSRx95Zo0uv1wKDw90dQCAi1RAA9KePXt08803ez/P+fdbiDMzM5Wbm6uf//znOn36tO655x5VVlbqu9/9rjZt2uR9B5IkvfDCC5oxY4ZGjBihoKAgZWRkaPny5a1+LmgBsbFcSgMAtAlt5j1IgcR7kAAAaH8uyvcgAQAABAoBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwIKABAAAYEFAAgAAsCAgAQAAWBCQAAAALAhIAAAAFgQkAAAACwISAACABQEJAADAgoAEAABgQUACAACwICABAABYEJAAAAAsCEgAAAAWHSYgrVixQr169VKnTp2UkpKiXbt2BbokAADQTnWIgPTSSy9pzpw5WrRokfbu3atBgwYpLS1NFRUVgS4NAAC0Qx0iID322GO6++67NWXKFCUnJ2vlypWKiIjQc889F+jSAABAOxQS6AK+rdraWhUWFio7O9u7LigoSKmpqSooKLjgNjU1NaqpqfF+drvdkqSqqir/FgsAAFpM4+9tY0yL77vdB6RPP/1U9fX1SkhI8FmfkJCggwcPXnCbnJwcLV68+Lz1PXr08EuNAADAf06cOCGHw9Gi+2z3AembyM7O1pw5c7yfKysr1bNnT5WWlrb4DxjNU1VVpR49eujYsWOy2+2BLueixli0HYxF28FYtC1ut1tJSUmKjY1t8X23+4B06aWXKjg4WOXl5T7ry8vLlZiYeMFtwsPDFR4eft56h8PBv/BthN1uZyzaCMai7WAs2g7Gom0JCmr5W6rb/U3aYWFhGjJkiLZs2eJd19DQoC1btsjpdAawMgAA0F61+xkkSZozZ44yMzN13XXX6YYbbtATTzyh06dPa8qUKYEuDQAAtEMdIiBNmDBBn3zyiRYuXCiXy6VrrrlGmzZtOu/G7S8THh6uRYsWXfCyG1oXY9F2MBZtB2PRdjAWbYs/x8Nm/PFsHAAAQDvW7u9BAgAAaGkEJAAAAAsCEgAAgAUBCQAAwOKiD0grVqxQr1691KlTJ6WkpGjXrl2BLqnDeeutt/SDH/xA3bp1k81m0yuvvOLTbozRwoUL1bVrV3Xu3Fmpqak6fPiwT5+TJ09q4sSJstvtiomJ0dSpU1VdXd2KZ9Ex5OTk6Prrr1d0dLTi4+M1btw4HTp0yKfP2bNnlZWVpS5duigqKkoZGRnnvYi1tLRU6enpioiIUHx8vObNm6dz58615qm0e08//bQGDhzofeGg0+nUxo0bve2MQ+AsXbpUNptNs2bN8q5jPFrPQw89JJvN5rP069fP295aY3FRB6SXXnpJc+bM0aJFi7R3714NGjRIaWlpqqioCHRpHcrp06c1aNAgrVix4oLty5Yt0/Lly7Vy5Urt3LlTkZGRSktL09mzZ719Jk6cqP379ysvL08bNmzQW2+9pXvuuae1TqHDyM/PV1ZWlt555x3l5eWprq5OI0eO1OnTp719Zs+erddee03r1q1Tfn6+jh8/rvHjx3vb6+vrlZ6ertraWu3YsUOrV69Wbm6uFi5cGIhTare6d++upUuXqrCwUHv27NEtt9yisWPHav/+/ZIYh0DZvXu3nnnmGQ0cONBnPePRuq666iqVlZV5l7ffftvb1mpjYS5iN9xwg8nKyvJ+rq+vN926dTM5OTkBrKpjk2TWr1/v/dzQ0GASExPNb3/7W++6yspKEx4ebl588UVjjDEHDhwwkszu3bu9fTZu3GhsNpv5+OOPW632jqiiosJIMvn5+cYYz88+NDTUrFu3ztvn/fffN5JMQUGBMcaYv/71ryYoKMi4XC5vn6efftrY7XZTU1PTuifQwVxyySXm2WefZRwC5NSpU+bKK680eXl5Zvjw4eb+++83xvD3orUtWrTIDBo06IJtrTkWF+0MUm1trQoLC5WamupdFxQUpNTUVBUUFASwsovL0aNH5XK5fMbB4XAoJSXFOw4FBQWKiYnRdddd5+2TmpqqoKAg7dy5s9Vr7kjcbrckeb/osbCwUHV1dT7j0a9fPyUlJfmMx4ABA3xexJqWlqaqqirv7Aeap76+XmvXrtXp06fldDoZhwDJyspSenq6z89d4u9FIBw+fFjdunXTd77zHU2cOFGlpaWSWncsOsSbtL+JTz/9VPX19ee9bTshIUEHDx4MUFUXH5fLJUkXHIfGNpfLpfj4eJ/2kJAQxcbGevug+RoaGjRr1iwNGzZMV199tSTPzzosLEwxMTE+fa3jcaHxamxD0xUXF8vpdOrs2bOKiorS+vXrlZycrKKiIsahla1du1Z79+7V7t27z2vj70XrSklJUW5urvr27auysjItXrxYN910k957771WHYuLNiABF7usrCy99957Ptf20br69u2roqIiud1u/fnPf1ZmZqby8/MDXdZF59ixY7r//vuVl5enTp06Bbqci97o0aO9fx44cKBSUlLUs2dP/elPf1Lnzp1brY6L9hLbpZdequDg4PPufC8vL1diYmKAqrr4NP6sv2ocEhMTz7tx/ty5czp58iRj9Q3NmDFDGzZs0Jtvvqnu3bt71ycmJqq2tlaVlZU+/a3jcaHxamxD04WFhemKK67QkCFDlJOTo0GDBunJJ59kHFpZYWGhKioqdO211yokJEQhISHKz8/X8uXLFRISooSEBMYjgGJiYtSnTx8dOXKkVf9uXLQBKSwsTEOGDNGWLVu86xoaGrRlyxY5nc4AVnZx6d27txITE33GoaqqSjt37vSOg9PpVGVlpQoLC719tm7dqoaGBqWkpLR6ze2ZMUYzZszQ+vXrtXXrVvXu3dunfciQIQoNDfUZj0OHDqm0tNRnPIqLi31Ca15enux2u5KTk1vnRDqohoYG1dTUMA6tbMSIESouLlZRUZF3ue666zRx4kTvnxmPwKmurtYHH3ygrl27tu7fjW90i3kHsXbtWhMeHm5yc3PNgQMHzD333GNiYmJ87nzHt3fq1Cnz7rvvmnfffddIMo899ph59913TUlJiTHGmKVLl5qYmBjz6quvmn379pmxY8ea3r17m88++8y7j1GjRpnBgwebnTt3mrfffttceeWV5s477wzUKbVb06dPNw6Hw2zbts2UlZV5lzNnznj73HvvvSYpKcls3brV7NmzxzidTuN0Or3t586dM1dffbUZOXKkKSoqMps2bTJxcXEmOzs7EKfUbv3iF78w+fn55ujRo2bfvn3mF7/4hbHZbGbz5s3GGMYh0L74FJsxjEdrmjt3rtm2bZs5evSo2b59u0lNTTWXXnqpqaioMMa03lhc1AHJGGN+97vfmaSkJBMWFmZuuOEG88477wS6pA7nzTffNJLOWzIzM40xnkf9H3zwQZOQkGDCw8PNiBEjzKFDh3z2ceLECXPnnXeaqKgoY7fbzZQpU8ypU6cCcDbt24XGQZJZtWqVt89nn31mfvazn5lLLrnEREREmB/+8IemrKzMZz8ffvihGT16tOncubO59NJLzdy5c01dXV0rn0379tOf/tT07NnThIWFmbi4ODNixAhvODKGcQg0a0BiPFrPhAkTTNeuXU1YWJi57LLLzIQJE8yRI0e87a01FjZjjPlWc18AAAAdzEV7DxIAAMCXISABAABYEJAAAAAsCEgAAAAWBCQAAAALAhIAAIAFAQkAAMCCgAQAAGBBQALQZvTq1Us2m002m+28L6Nsz3Jzc73nNWvWrECXA6AJCEgAWlR9fb1uvPFGjR8/3me92+1Wjx499Ktf/eort1+yZInKysrkcDj8WaZyc3MVExPj12M0mjBhgsrKyvgibKAdISABaFHBwcHKzc3Vpk2b9MILL3jXz5w5U7GxsVq0aNFXbh8dHa3ExETZbDZ/l9oi6uvr1dDQ8JV9OnfurMTERIWFhbVSVQC+LQISgBbXp08fLV26VDNnzlRZWZleffVVrV27Vs8//3yzQ0LjTM+GDRvUt29fRURE6I477tCZM2e0evVq9erVS5dcconuu+8+1dfXe7erqanRAw88oMsuu0yRkZFKSUnRtm3bJEnbtm3TlClT5Ha7vZe+Hnrooa/d7ov1/OUvf1FycrLCw8NVWlqqbdu26YYbblBkZKRiYmI0bNgwlZSUfNsfJYAACQl0AQA6ppkzZ2r9+vW66667VFxcrIULF2rQoEHfaF9nzpzR8uXLtXbtWp06dUrjx4/XD3/4Q8XExOivf/2r/vWvfykjI0PDhg3ThAkTJEkzZszQgQMHtHbtWnXr1k3r16/XqFGjVFxcrBtvvFFPPPGEFi5cqEOHDkmSoqKivna7K6+80lvPb37zGz377LPq0qWLYmNjdc011+juu+/Wiy++qNraWu3atavdzIIBuAADAH7y/vvvG0lmwIABpq6u7mv79+zZ0zz++OM+61atWmUkmSNHjnjXTZs2zURERJhTp05516WlpZlp06YZY4wpKSkxwcHB5uOPP/bZ14gRI0x2drZ3vw6Hw6e9qdtJMkVFRd72EydOGElm27ZtX3l+w4cPN/fff/9X9gHQNjCDBMBvnnvuOUVEROjo0aP66KOP1KtXr2+0n4iICF1++eXezwkJCerVq5d31qdxXUVFhSSpuLhY9fX16tOnj89+ampq1KVLly89TlO3CwsL08CBA72fY2NjNXnyZKWlpenWW29VamqqfvzjH6tr167f6HwBBB4BCYBf7NixQ48//rg2b96sX//615o6dar+9re/faPLTqGhoT6fbTbbBdc13ixdXV2t4OBgFRYWKjg42KffF0OVVVO369y583nnsWrVKt13333atGmTXnrpJS1YsEB5eXkaOnRo008UQJtBQALQ4s6cOaPJkydr+vTpuvnmm9W7d28NGDBAK1eu1PTp0/1+/MGDB6u+vl4VFRW66aabLtgnLCzM56bupm73dccdPHiwsrOz5XQ6tWbNGgIS0E7xFBuAFpednS1jjJYuXSrJ8wLI//7v/9bPf/5zffjhh34/fp8+fTRx4kRNmjRJL7/8so4ePapdu3YpJydHr7/+urem6upqbdmyRZ9++qnOnDnTpO0u5OjRo8rOzlZBQYFKSkq0efNmHT58WP379/f7uQLwDwISgBaVn5+vFStWaNWqVYqIiPCunzZtmm688UZNnTpVxhi/17Fq1SpNmjRJc+fOVd++fTVu3Djt3r1bSUlJkqQbb7xR9957ryZMmKC4uDgtW7asSdtdSEREhA4ePKiMjAz16dNH99xzj7KysjRt2jS/nycA/7CZ1vgvFQA0Qa9evTRr1qwO+3Uc3//+93XNNdfoiSeeCHQpAL4GM0gA2pT58+crKipKbrc70KW0mBdeeEFRUVH6+9//HuhSADQRM0gA2oySkhLV1dVJkr7zne8oKKhj/D/cqVOnVF5eLkmKiYnRpZdeGuCKAHwdAhIAAIBFx/jfMwAAgBZEQAIAALAgIAEAAFgQkAAAACwISAAAABYEJAAAAAsCEgAAgAUBCQAAwOL/A8FWFMInKZUgAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "landscape.plot()" ] }, { "cell_type": "code", "execution_count": 7, "id": "df4715c1", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAG2CAYAAACEbnlbAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAARYNJREFUeJzt3Xd4VGX+/vF7UkkhiUFIaKGogFGagDi42IgGFhQkuKgoyKIiBpSiq1EEddUg7iq6gli+UlREcSmC4IoIWIiUIEsVRJEgEEL5JaFIEpLz++PZDGQImEAyZzJ5v65rLmbOOTPzOZx15+Y5T3FYlmUJAAAALn52FwAAAOBtCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABubA1ITz/9tBwOR4lHixYtXPuPHz+u5ORk1apVS+Hh4UpKStK+fftKfEZGRoa6d++u0NBQ1alTR48++qhOnDjh6VMBAAA+JMDuAi677DJ9+eWXrtcBASdLGjFihD777DPNmjVLkZGRGjp0qHr37q3vvvtOklRYWKju3bsrNjZWK1as0N69e9W/f38FBgbqhRde8Pi5AAAA3+Cwc7Hap59+WnPnztW6detO25eTk6PatWtrxowZ6tOnjyTpxx9/1KWXXqq0tDRdddVVWrRokXr06KE9e/YoJiZGkjR58mQ99thj2r9/v4KCgjx5OgAAwEfY3oL0008/qV69eqpRo4acTqdSU1MVFxen9PR0FRQUKCEhwXVsixYtFBcX5wpIaWlpatmypSscSVJiYqKGDBmiTZs2qW3btqV+Z15envLy8lyvi4qKdOjQIdWqVUsOh6PyThYAAFQYy7J0+PBh1atXT35+FdtryNaA1LFjR02dOlXNmzfX3r179cwzz6hz587auHGjMjMzFRQUpKioqBLviYmJUWZmpiQpMzOzRDgq3l+870xSU1P1zDPPVOzJAAAAW+zatUsNGjSo0M+0NSB169bN9bxVq1bq2LGjGjVqpI8//lghISGV9r0pKSkaOXKk63VOTo7i4uK0a9cuRUREVNr3AgCAipObm6uGDRuqZs2aFf7Ztt9iO1VUVJSaNWum7du368Ybb1R+fr6ys7NLtCLt27dPsbGxkqTY2FitWrWqxGcUj3IrPqY0wcHBCg4OPm17REQEAQkAgCqmMrrHeNU8SEeOHNHPP/+sunXrql27dgoMDNSSJUtc+7du3aqMjAw5nU5JktPp1IYNG5SVleU6ZvHixYqIiFB8fLzH6wcAAL7B1hakRx55RDfffLMaNWqkPXv2aOzYsfL399cdd9yhyMhIDRo0SCNHjlR0dLQiIiI0bNgwOZ1OXXXVVZKkm266SfHx8br77rs1fvx4ZWZmavTo0UpOTi61hQgAAKAsbA1Iv/32m+644w4dPHhQtWvX1p/+9Cd9//33ql27tiTplVdekZ+fn5KSkpSXl6fExERNmjTJ9X5/f38tWLBAQ4YMkdPpVFhYmAYMGKBnn33WrlMCAAA+wNZ5kLxFbm6uIiMjlZOTQx8kAACqiMr8/faqPkgAAADegIAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADghoAEAADgxmsC0rhx4+RwODR8+HDXtuPHjys5OVm1atVSeHi4kpKStG/fvhLvy8jIUPfu3RUaGqo6dero0Ucf1YkTJzxcPQAA8CVeEZBWr16tN998U61atSqxfcSIEZo/f75mzZql5cuXa8+ePerdu7drf2Fhobp37678/HytWLFC06ZN09SpUzVmzBhPnwIAAPAhtgekI0eOqF+/fnr77bd1wQUXuLbn5OTo//7v//Tyyy/rhhtuULt27TRlyhStWLFC33//vSTpiy++0ObNm/X++++rTZs26tatm/7+979r4sSJys/Pt+uUAABAFWd7QEpOTlb37t2VkJBQYnt6eroKCgpKbG/RooXi4uKUlpYmSUpLS1PLli0VExPjOiYxMVG5ubnatGnTGb8zLy9Pubm5JR4AAADFAuz88pkzZ2rt2rVavXr1afsyMzMVFBSkqKioEttjYmKUmZnpOubUcFS8v3jfmaSmpuqZZ545z+oBAICvsq0FadeuXXr44Yf1wQcfqEaNGh797pSUFOXk5Lgeu3bt8uj3AwAA72ZbQEpPT1dWVpauuOIKBQQEKCAgQMuXL9drr72mgIAAxcTEKD8/X9nZ2SXet2/fPsXGxkqSYmNjTxvVVvy6+JjSBAcHKyIiosQDAACgmG0BqUuXLtqwYYPWrVvnerRv3179+vVzPQ8MDNSSJUtc79m6dasyMjLkdDolSU6nUxs2bFBWVpbrmMWLFysiIkLx8fEePycAAOAbbOuDVLNmTV1++eUltoWFhalWrVqu7YMGDdLIkSMVHR2tiIgIDRs2TE6nU1dddZUk6aabblJ8fLzuvvtujR8/XpmZmRo9erSSk5MVHBzs8XMCAAC+wdZO2n/klVdekZ+fn5KSkpSXl6fExERNmjTJtd/f318LFizQkCFD5HQ6FRYWpgEDBujZZ5+1sWoAAFDVOSzLsuwuwm65ubmKjIxUTk4O/ZEAAKgiKvP32/Z5kAAAALwNAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAQkAAMANAekU33wjFRbaXQUAALAbAekUPXpIjRtLs2fbXQkAALATAcnN7t1Snz6EJAAAqjMCkhvLMn8OH87tNgAAqisCUiksS9q1y/RJAgAA1Q8B6Sz27rW7AgAAYAcC0lnUrWt3BQAAwA4EpDNo2FDq3NnuKgAAgB0ISGcwZIjk7293FQAAwA4EJDchIebPCROknTttLQUAANiEgHSKBQtMx+w2baSsLDNxZG6u3VUBAABPIyCdonNnKTJSmj/fdNDeuFG64w7mQwIAoLohIJWiQQNp3jypRg1p4ULpkUfsrggAAHgSAekMOnSQpk83zydMkN5809ZyAACABxGQzuK226TnnjPPk5OlL7+0tx4AAOAZBKQ/8MQT0l13mX5It90mbd1qd0UAAKCyEZD+gMMhvf221KmTlJ1tRrYdPGh3VQAAoDIRkMqgRg1pzhypcWNp+3YpKUnKz7e7KgAAUFkISGVUp46ZJ6lmTWn5cjPTtmXZXRUAAKgMBKRyuOwy6eOPJT8/6d13pX/8w+6KAABAZSAglVPXrmbYvyQ99piZLwkAAPgWAtI5GDpUevBBc4vtzjulH36wuyIAAFCRCEjnwOGQXn1VuvFG6dgx6eabpT177K4KAABUFALSOQoIMP2RWrSQdu+WevY0YQkAAFR9BKTzEBVlRrbVqiWtWSMNGCAVFdldFQAAOF8EpPN00UVmjqTAQOmTT6SxY+2uCAAAnC8CUgXo3NnMti2Ztdvef9/eegAAwPkhIFWQAQOkxx83zwcNkr77zt56AADAuSMgVaDnn5duvdUsQ3LrrdKOHXZXBAAAzgUBqQL5+UnvvSddcYW0f79Z2DYnx+6qAABAeRGQKlhYmPTpp1K9etLmzdLtt0snTthdFQAAKA8CUiWoX9+EpJAQ6fPPpVGj7K4IAACUBwGpkrRrd3I022uvSZMm2VsPAAAoOwJSJerdW3rhBfP8oYekL76wtx4AAFA2BKRK9vjjUv/+UmGhdNtt0pYtdlcEAAD+iK0B6Y033lCrVq0UERGhiIgIOZ1OLVq0yLX/+PHjSk5OVq1atRQeHq6kpCTt27evxGdkZGSoe/fuCg0NVZ06dfToo4/qhBf1inY4pLfekv70Jyk314xsO3DA7qoAAMDZ2BqQGjRooHHjxik9PV1r1qzRDTfcoJ49e2rTpk2SpBEjRmj+/PmaNWuWli9frj179qh3796u9xcWFqp79+7Kz8/XihUrNG3aNE2dOlVjxoyx65RKFRxsliNp2lT65Rdz6y0vz+6qAADAmTgsy7LsLuJU0dHReumll9SnTx/Vrl1bM2bMUJ8+fSRJP/74oy699FKlpaXpqquu0qJFi9SjRw/t2bNHMTExkqTJkyfrscce0/79+xUUFFSm78zNzVVkZKRycnIUERFRaee2ebPkdJqWpAEDpClTTAsTAAAov8r8/faaPkiFhYWaOXOmjh49KqfTqfT0dBUUFCghIcF1TIsWLRQXF6e0tDRJUlpamlq2bOkKR5KUmJio3NxcVytUafLy8pSbm1vi4Qnx8dLHH0v+/tK0adKLL3rkawEAQDnZHpA2bNig8PBwBQcH64EHHtCcOXMUHx+vzMxMBQUFKSoqqsTxMTExyszMlCRlZmaWCEfF+4v3nUlqaqoiIyNdj4YNG1bsSZ1FYqIZ9i9JKSnS7Nke+2oAAFBGtgek5s2ba926dVq5cqWGDBmiAQMGaPPmzZX6nSkpKcrJyXE9du3aVanf5+7BB6WhQ83zu+6S0tM9+vUAAOAP2B6QgoKCdPHFF6tdu3ZKTU1V69at9eqrryo2Nlb5+fnKzs4ucfy+ffsUGxsrSYqNjT1tVFvx6+JjShMcHOwaOVf88LRXXjGtSb//Lt1yi7R7t8dLAAAAZxBQloNeK74nVA4DBw5UzZo1y/2+oqIi5eXlqV27dgoMDNSSJUuUlJQkSdq6dasyMjLkdDolSU6nU88//7yysrJUp04dSdLixYsVERGh+Pj4cn+3JwUESB99JHXqZDpv33KL9PXXZi03AABgrzKNYvPz81ODBg3k7+9fpg/dtWuXtm3bpqZNm571uJSUFHXr1k1xcXE6fPiwZsyYoRdffFH/+c9/dOONN2rIkCFauHChpk6dqoiICA0bNkyStGLFCkmmY3ebNm1Ur149jR8/XpmZmbr77rt177336oXiKazLwFOj2EqzY4d05ZVmbqTevaVZsyQ/29v1AADwfpX5+12mFiRJWrNmjauV5o+UteUoKytL/fv31969exUZGalWrVq5wpEkvfLKK/Lz81NSUpLy8vKUmJioSacsaubv768FCxZoyJAhcjqdCgsL04ABA/Tss8+W9bRs16SJNHeudMMNpsP26NEnlycBAAD2KFML0jPPPKNHH31UoaGhZfrQ1NRUDRky5LQRaN7KzhakYu+9Z5YkkaSpU808SQAA4Mwq8/fb6yaKtIM3BCTJtB49/7wUGCgtWSJ17mxbKQAAeD2vmijy999/17Fjx1yvd+7cqQkTJugLlqo/b88+K/XpIxUUSLfeKv38s90VAQBQPZU7IPXs2VPTp0+XJGVnZ6tjx4765z//qZ49e+qNN96o8AKrEz8/M8N2u3bSwYPSzTdLbrMcAAAADyh3QFq7dq06/+/ezyeffKKYmBjt3LlT06dPP6fpAFBSaKj06adS/frSli1S377SiRN2VwUAQPVS7oB07Ngx1yi1L774Qr1795afn5+uuuoq7dy5s8ILrI7q1ZPmzzdh6YsvpOHD7a4IAIDqpdwB6eKLL9bcuXO1a9cu/ec//9FNN90kyQzZt7ODs69p21b64APJ4ZAmTpRef93uigAAqD7KHZDGjBmjRx55RI0bN1bHjh1ds1p/8cUXatu2bYUXWJ316iWNG2eeP/yw9PnntpYDAEC1cU7D/DMzM7V37161bt1afv+b9nnVqlWKiIhQixYtKrzIyuYtw/xLY1nSoEHSlClSzZpSWpp02WV2VwUAgP28Zph/QUGBAgICdODAAbVt29YVjiTpyiuvrJLhyNs5HNLkydI110iHD0s9ekj799tdFQAAvq1cASkwMFBxcXEqLCysrHpQiqAgswzJRRdJv/5qbr0dP253VQAA+K5y90F68skn9cQTT+jQoUOVUQ/OoFYtacECKTJSWrFCuu8+c/sNAABUvHL3QWrbtq22b9+ugoICNWrUSGFhYSX2r127tkIL9ARv7oPk7ssvpa5dpcJC6bnnpCeftLsiAADsUZm/3wHlfUOvXr0qtACUT0KCGfI/ZIhZu615c7M8CQAAqDgsVquq1YJUbPhw6dVXpZAQaflyqUMHuysCAMCzvGYUW7Hs7Gy98847SklJcfVFWrt2rXbv3l2hxeHM/vlP6c9/ln7/XbrlFmnXLrsrAgDAd5Q7IK1fv17NmjXTiy++qH/84x/K/t9qqrNnz1ZKSkpF14cz8PeXPvxQuvxyKTPThKQjR+yuCgAA31DugDRy5Ejdc889+umnn1SjRg3X9j//+c/6+uuvK7Q4nF1EhFmzrXZtad066a67pKIiu6sCAKDqK3dAWr16tQYPHnza9vr16yszM7NCikLZNW4szZ0rBQdL8+ZJNOIBAHD+yh2QgoODlZube9r2bdu2qXbt2hVSFMqnUyfp3XfN8/HjzbIkAADg3JU7IN1yyy169tlnVVBQIElyOBzKyMjQY489pqSkpAovEGVz553SU0+Z54MHm5FthYXSsmWmr9KyZeY1AAD4Y+Ue5p+Tk6M+ffpozZo1Onz4sOrVq6fMzEw5nU4tXLjwtIkjq4KqOMy/NEVF0h13SB9/LIWHm8epdz0bNDBTA/TubV+NAABUlMr8/T7neZC+++47/fe//9WRI0d0xRVXKCEhoUIL8yRfCUiSGfbfqpW0ffvp+xwO8+cnnxCSAABVn1cFpOnTp6tv374KDg4usT0/P18zZ85U//79K7RAT/ClgFRYKDVsKO3dW/p+h8O0JO3YYaYKAACgqvKqiSIHDhyonJyc07YfPnxYAwcOrJCicO6++ebM4UgyC9zu2mU6c//0k/S/rmQAAOAU5V6LzbIsOYrv1Zzit99+U2RkZIUUhXN3tnB0qieeMI+AAKlpU6lZM/No3vzk87p1T96WAwCgOilzQGrbtq0cDoccDoe6dOmigICTby0sLNSOHTvUtWvXSikSZVe3btmOu+giac8e02dp2zbzcBcWdnpoKn6QhQEAvqzMAalXr16SpHXr1ikxMVHh4eGufUFBQWrcuDHD/L1A586mj9Hu3eZ2mrviPkhbt5rnu3efDEjbtpnt27aZPkpHj0o//GAe7mJiSgam4hDVtKmZtBIAgKqs3J20p02bpr59+5ZYZqSq86VO2pI0e7bUp495furVLc8otvx86ZdfSoam4sfZJkz38zOze5fW8tSggdkPAEBF8KpRbJKUnZ2tTz75RD///LMeffRRRUdHa+3atYqJiVH9+vUrtEBP8LWAJJmQ9PDD0m+/ndzWsKE0YcL5D/HPzS0ZmE4NUWdbMDckRLrkktNv1zVvLkVHn19NAIDqx6sC0vr165WQkKDIyEj9+uuv2rp1q5o2barRo0crIyND06dPr9ACPcEXA5JkhvwXj2qrW9fcfqvMof2WZVqX3EPTtm3Szz9LJ06c+b21apV+y+7ii02wAgDAnVcFpC5duqhdu3YaP368atasqf/+979q2rSpVqxYoTvvvFO//vprhRboCb4akLzJiROmX1NpLU+7d5/5fQ6HafkqraN4o0YVG/g8HSgBAOfHqwJSZGSk1q5dq4suuqhEQNq5c6eaN2+u48ePV2iBnkBAsteRI2bmb/dWp61bpVKm3HIJCjItTKW1PNWuXb4pCkq7JcnSLADg3Srz97vc8yAFBwcrNzf3tO3btm1T7dq1K6QoVC/h4VKbNuZxKsuSDhw4vZP4tm1mksv8fGnzZvNwFxlZekfxSy4x33eq4k7t7v9U2L3bbGdpFgCofsrdgnTvvffq4MGD+vjjjxUdHa3169fL399fvXr10jXXXKMJEyZUUqmVhxakqqewUMrIKL2/U0ZG6VMcFKtf/2Rguvhiadw46eDB0o9laRYA8F5edYstJydHffr00Zo1a3T48GHVq1dPmZmZcjqdWrhwocLCwiq0QE8gIPmW3383ncJLa3k6cODcPnPpUum66yq0TADAefKqW2yRkZFavHixvv32W61fv15HjhzRFVdcoYSEhAotDDhXISHS5Zebh7tDh0oGpsWLpVWr/vgzy7qECwDAN5zTPEi+hhak6mvZMun66//4OFqQAMD7eFULkiStXr1aS5cuVVZWloqKikrse/nllyukMMAT/mhplmIffGA6kUdFeaoyAICdyh2QXnjhBY0ePVrNmzdXTEyMHKeMpXaw9DuqGH9/M5S/Tx/TIdt9aZbi1++8I332mfT664xoA4DqoNy32GJiYvTiiy/qnnvuqaSSPI9bbDjb0izR0dL995upBSSpVy8TlKrgqjoA4FMq8/e73EuH+vn56eqrr67QIgC79e4t/fqr6Ws0Y4b5c8cOs/2666T166UnnpACAqS5c6X4eOmNNyS3O8wAAB9R7hak8ePHa8+ePVVyvqMzoQUJZbV+vXTffSdHvl19tfTWWyYwAQA8y6vmQSoqKlL37t21bds2xcfHKzAwsMT+2bNnV2iBnkBAQnkUFkoTJ5oWpaNHpcBA6cknpccfl4KD7a4OAKoPr7rF9tBDD2np0qVq1qyZatWqpcjIyBIPwNf5+0sPPWSWOPnzn6WCAunpp6W2baXvvrO7OgBARSh3C1LNmjU1c+ZMde/evbJq8jhakHCuLEv66CPTwTsry2wbMkRKTTXrwQEAKo9XtSBFR0froosuqtAigKrK4ZBuv13askUaONBse+MN0ydp7lxbSwMAnIdyB6Snn35aY8eO1bFjxyqjHqBKio6W3n1XWrJEuugiac8e6dZbpaQk8xwAULWU+xZb27Zt9fPPP8uyLDVu3Pi0Ttpr166t0AI9gVtsqEi//y49+6z00kumQ3dkpDR+vHTvvZJfuf9JAgA4E69aaqRXr14VWgDga0JCTB+kvn3NlABr1kiDB0vvv2+mBGjRwu4KAQB/hMVqRQsSKk9hofTaa9Lo0dKxY1JQkHn+2GPmOQDg3HlVJ20AZefvL40YIW3aJHXtKuXnS2PGSFdcIaWl2V0dAOBMyhSQoqOjdeDAgTJ/aFxcnHbu3HnORQG+pnFjaeFC6YMPpAsvNIHp6qulYcOkw4ftrg4A4K5MfZCys7O1aNGiMk8EefDgQRUWFp5XYYCvcTikO++UEhOlUaOkadPMordz50qTJkk332x3hQCAYmXqg+R3DkNvtm/frqZNm55TUZ5GHyTYYfFi6YEHpF9+Ma9vu830V4qNtbcuAKgqbO+DVFRUVO5HVQlHgF1uvFHasEF69FHTV2nWLOnSS6V33jEzdAMA7EMnbcBGoaFmjqTVq03H7exsMzXA9ddL27bZXR0AVF8EJMALtG0rrVwp/eMfZh6l5culVq2kF14wi+ECADyLgAR4iYAA03l740Zz+y0vT3rySaldOxOeAACeU+aAtIcFpQCPaNpU+s9/pOnTpVq1TD8lp1N6+GGmBAAATylzQLrssss0Y8aMyqwFwP84HNLdd0tbtkh33WU6bb/2mnTZZdJnn9ldHQD4vjIHpOeff16DBw/WbbfdpkOHDlVmTQD+p3Zt6b33pM8/N5NN7tol9egh3XGHtG+f3dUBgO8qc0B68MEHtX79eh08eFDx8fGaP39+ZdYF4BSJiaZv0qhRkp+fNHOmmRJgyhSmBACAynBOi9W+/vrrGjFihC699FIFBJScjHvt2rUVVpynMFEkqpI1a8xUAOvWmdc33CC9+aZ08cW2lgUAHmf7RJGn2rlzp2bPnq0LLrhAPXv2PO1RHqmpqerQoYNq1qypOnXqqFevXtq6dWuJY44fP67k5GTVqlVL4eHhSkpK0j63ewsZGRnq3r27QkNDVadOHT366KM6ceJEeU8NqBLat5dWrZJefFGqUUP66iupZUtp3DimBACAilKuFqS3335bo0aNUkJCgt58803Vrl37vL68a9euuv3229WhQwedOHFCTzzxhDZu3KjNmzcrLCxMkjRkyBB99tlnmjp1qiIjIzV06FD5+fnpu+++kyQVFhaqTZs2io2N1UsvvaS9e/eqf//+uu+++/TCCy+UqQ5akFBV/fyzNHiwtGSJed26tZmJu317e+sCAE+o1N9vq4wSExOtCy64wJo2bVpZ31JuWVlZliRr+fLllmVZVnZ2thUYGGjNmjXLdcyWLVssSVZaWpplWZa1cOFCy8/Pz8rMzHQd88Ybb1gRERFWXl5emb43JyfHkmTl5ORU4NkAnlFUZFlTplhWdLRlSZbl52dZI0ZY1uHDdlcGAJWrMn+/y3yLrbCwUOvXr1f//v0rNqGdIicnR5IUHR0tSUpPT1dBQYESEhJcx7Ro0UJxcXFKS0uTJKWlpally5aKiYlxHZOYmKjc3Fxt2rSp1O/Jy8tTbm5uiQdQVTkc0j33mCkB7rhDKiqSXnlFuvxyM/oNAFB+ZQ5IixcvVoMGDSqtkKKiIg0fPlxXX321Lr/8cklSZmamgoKCFBUVVeLYmJgYZWZmuo45NRwV7y/eV5rU1FRFRka6Hg0bNqzgswE8r04dacYMaeFCKS5O2rlT6tZN6tdP2r/f7uoAoGrxmqVGkpOTtXHjRs2cObPSvyslJUU5OTmux65duyr9OwFP6dZN2rRJGj7cTAkwY4bUooU0bRpTAgBAWXlFQBo6dKgWLFigpUuXlmilio2NVX5+vrKzs0scv2/fPsXGxrqOcR/VVvy6+Bh3wcHBioiIKPEAfEl4uLnN9v33ZtHbQ4fMbbibbjIduwEAZ2drQLIsS0OHDtWcOXP01VdfqUmTJiX2t2vXToGBgVpSPERH0tatW5WRkSGn0ylJcjqd2rBhg7KyslzHLF68WBEREYqPj/fMiQBeqkMHM29SaqqZEuDLL82UAC+9JDETBgCc2TlNFFlRHnzwQc2YMUPz5s1T8+bNXdsjIyMVEhIiyQzzX7hwoaZOnaqIiAgNGzZMkrRixQpJJ4f516tXT+PHj1dmZqbuvvtu3XvvvQzzB07x009mSoClS83rtm3NlABXXGFvXQBwrirz99vWgORwOErdPmXKFN1zzz2SzESRo0aN0ocffqi8vDwlJiZq0qRJJW6f7dy5U0OGDNGyZcsUFhamAQMGaNy4cafN8n0mBCRUF5Zllid55BHp//0/00dpxAjpmWek/009BgBVhs8GJG9BQEJ1s2+f9PDD0kcfmddNmkiTJ5s+SgBQVXjVUiMAqr6YGLPg7YIFUsOG0o4dZkHc/v2lAwfsrg4A7EdAAqqx7t3NlAAPPWQmnHzvPenSS6X332dKAADVGwEJqOZq1pRefVVascLMvn3ggHT33WY+pR077K4OAOxBQAIgSbrqKik9XXruOSk4WPrPf0xgevllpgQAUP0QkAC4BAVJTz4prV8vXXutdOyYNGqUCU/r1tldHQB4DgEJwGmaNZO++kp6+20pMtK0LLVvLz32mAlNAODrCEgASuXnJ917r7Rli9Snj1RYKI0fb2biPmVyewDwSQQkAGdVt640a5Y0b55Uv770yy9SQoJZ2+3gQburA4DKQUACUCa33CJt3iwlJ5spAaZNM1MCzJjBlAAAfA8BCUCZRURIr78uffutFB8v7d8v9etn5lPaudPu6gCg4hCQAJRbp07SDz9Izz5rRr4tWiRddpk0YYLpqwQAVR0BCcA5CQqSnnpK+u9/pc6dpaNHzcK3TqfZBgBVGQEJwHlp0UJatswsdhsRIa1ebaYESEmRfv/d7uoA4NwQkACcNz8/afBgMyVA795m5u1x46RWraSlS+2uDgDKj4AEoMLUqyf9+9/SnDnm+fbt0g03SIMGSYcO2V0dAJQdAQlAhevVy0wJMGSIef3uu2ZKgI8+YkoAAFUDAQlApYiMlCZNkr75xvRTysqSbr9duvlmKSPD7uoA4OwISAAq1Z/+ZBa6HTtWCgyUPvvMTAnwr38xJQAA70VAAlDpgoOlp582QalTJ+nIEemhh6Srr5Y2bLC7OgA4HQEJgMfEx5tbbpMmSTVrSitXSldcIY0eLR0/bnd1AHASAQmAR/n5mc7bmzdLPXuaKQGef15q3Vpavtzu6gDAICABsEWDBmY6gE8+kWJjpW3bpOuuk+6/X8rOtrs6ANUdAQmAbRwOKSnJTDB5//1m29tvmykBPvmEKQEA2IeABMB2UVHSm2+aW2zNmkmZmdJtt5n5lH77ze7qAFRHBCQAXuOaa8xCt6NHSwEB0qefmo7dEydKRUV2VwegOiEgAfAqNWpIf/+79MMP0lVXSYcPS0OHmvmUNm2yuzoA1QUBCYBXuvxy6dtvzYSS4eFSWprUtq00ZoyUl2d3dQB8HQEJgNfy9zetR5s3myVKCgpM61KbNiY8AUBlISAB8HoNG0rz5kkffyzFxEg//ih17iw98ICUk2N3dQB8EQEJQJXgcJiRbVu2SIMGmW1vvmmmBJg9297aAPgeAhKAKuWCC6R33pGWLpUuuUTau9fMpdS7t7Rnj93VAfAVBCQAVdJ115kpAZ54wkwJMGeOaU2aPJkpAQCcPwISgCorJMSs45aeLl15pZSba9Z5u/ZacysOAM4VAQlAldeqlbRihfTqq1JYmBnh1qaN9MwzTAkA4NwQkAD4BH9/6aGHzGSSf/6zlJ8vPf20mTtpxQq7qwNQ1RCQAPiURo2kBQukDz+U6tQxt9r+9CcpOdncggOAsiAgAfA5Dod0++0mHA0cKFmWNGmSWddt3jy7qwNQFRCQAPis6Gjp3XelL7+ULrpI2r1b6tVL6tPHTA8AAGdCQALg87p0kTZskB57zPRV+ve/zZQAb73FlAAASkdAAlAthIRI48ZJa9ZI7dubJUoGD5auv17auvXkcYWF0rJlpg/TsmXmNYDqh4AEoFpp00b6/nvp5Zel0FDp66/NNAF//7tZ661xYxOa7rzT/Nm4MUuZANWRw7Isy+4i7Jabm6vIyEjl5OQoIiLC7nIAeMivv5qJJT///MzHOBzmz08+McuZAPAelfn7TQsSgGqrcWNp4ULpvfckvzP8v2HxPyGHD+d2G1CdEJAAVGsOh9Sgwdk7a1uWtGuX9M03nqsLgL0ISACqvbIO+WdqAKD6ICABqPbq1q3Y4wBUfQQkANVe587mNltxh+zSNGxojgNQPRCQAFR7/v7Sq6+a52cKSU6nOQ5A9UBAAgCZIfyffCLVr19y+wUXmD8//lj64APP1wXAHgQkAPif3r3N3EhLl0ozZpg/9++X/vY3s/+vf2UkG1BdMFGkmCgSwNkVFUm33WZm1K5Vy8zEffHFdlcFgIkiAcBGfn5mMskOHaSDB6Xu3aVDh+yuCkBlIiABQBmEhkqffirFxUnbtklJSVJ+vt1VAagsBCQAKKPYWGnBAqlmTWnZMun++08uRQLAtxCQ4J327zeriMbFScHB5pcpMVH67ju7Kzs3X38t3XyzVK+eGUc+d67dFeEctWxpRrT5+0vTpkmpqXZXBKAyEJDgnZKSpB9+ML9A27aZexvXXWc6gFSmyrpncvSo1Lq1NHFi5Xw+PKprV+lf/zLPn3zSBCYAvoWABO+TnW3GUr/4onT99VKjRtKVV0opKdItt5w8LiND6tlTCg+XIiKkv/xF2rfv5P577pF69Sr52cOHm6BV7LrrpKFDzfYLLzStVJK0aZPUo4f53Jo1zRTKP/988n3vvCNdeqlUo4bUooU0adLZz6lbN+m556Rbby3nXwa81ZAh5n82ktS/vxnZBsB3EJDgfcLDzWPuXCkvr/RjiopMODp0SFq+XFq8WPrlF6lv3/J/37RpUlCQuX03ebK0e7d0zTXm1t5XX0np6WYCnBMnzPEffCCNGSM9/7y0ZYv0wgvSU0+Zz0G18o9/mDuneXkmu+/YYXdFACpKgN0FAKcJCJCmTpXuu88EliuukK69Vrr9dqlVK3PMkiXShg3mF6lhQ7Nt+nTpssuk1avNeOyyuuQSafz4k6+feEKKjJRmzpQCA822Zs1O7h87VvrnP82sgpLUpIm0ebP05pvSgAHnfNqoevz9zYSS11xj7gj36GFydlSU3ZUBOF+0IME7JSVJe/aYvkddu5ohQ1dcYYKTZFpuGjY8GY4kKT7e/DJt2VK+72rXruTrdevMLbXicHSqo0fNrbZBg062dIWHm9tnp96CQ7URHi7Nn2+WKNm82UwoWVBgd1UAzhcBCd6rRg3pxhvN7asVK0yforFjy/5+P7/Tx2CX9ssVFlbydUjImT/zyBHz59tvmyBV/Ni4kU4o1Vj9+iYkhYVJX34pJScz/B+o6ghIqDri400LjmQ6SO/aZR7FNm82Hbzj483r2rWlvXtLfsa6dX/8Pa1amU7ipYWpmBgzVP+XX8xaE6c+mjQ5l7OCj2jbVvrwQ5PL337b3IUFUHURkOAdtm41PV6HDZMGDzYh5403pPXrTT+jWbNMP6GePc3xCQlmQpp+/aS1a6VVq8xQomuvldq3N8fccIO0Zo3pm/TTT6b1aePGP65l6FApN9f0eVqzxrz3vfdMjZL0zDNm8pvXXjNTEGzYIE2ZIr388pk/88iRk61NkjmndevMSDz4jJtvPvk/g7/9TZozx956AJwHy0bLly+3evToYdWtW9eSZM2ZM6fE/qKiIuupp56yYmNjrRo1alhdunSxtm3bVuKYgwcPWnfeeadVs2ZNKzIy0vrrX/9qHT58uFx15OTkWJKsnJyc8z0llNeuXZZ17bWWJVmWn59lBQZaVkCAeS1ZVnCwZYWGWlbz5pY1erRlHTt28r07d1rWLbdYVliYZdWsaVm33WZZmZklP3/MGMuKibGsyEjLGjHCsoYONd9X7NprLevhh0+v67//taybbjLfXbOmZXXubFk//3xy/wcfWFabNpYVFGRZF1xgWddcY1mzZ5/5PJcuPXlOpz4GDCjnXxi8XVGRZT34oLm8ISGWtXq13RUBvqsyf78dlmXfnfJFixbpu+++U7t27dS7d2/NmTNHvU6Zt+bFF19Uamqqpk2bpiZNmuipp57Shg0btHnzZtWoUUOS1K1bN+3du1dvvvmmCgoKNHDgQHXo0EEzZswocx2VuRowziI3V2reXDpw4OQQencOh2nJKcf1BOx24oRpTfr8czMJ/MqVZlJ4ABWrUn+/KzxynSO5tSAVFRVZsbGx1ksvveTalp2dbQUHB1sffvihZVmWtXnzZkuStfqUf6ItWrTIcjgc1u7du8v83bQg2eTJJ02rUWktK+6PFSvsrhYol5wcy2rZ0vzPt2VL8xpAxarM32+v7YO0Y8cOZWZmKiEhwbUtMjJSHTt2VFpamiQpLS1NUVFRal/c50RSQkKC/Pz8tHLlyjN+dl5ennJzc0s8YIP//MdM+PhHAgLMhI1AFRIRYRa2jY013dRuv/3MDaUAvI/XBqTMzExJUkxMTIntMTExrn2ZmZmqU6dOif0BAQGKjo52HVOa1NRURUZGuh4NT51LB54TEWGG/PyRoiIz2QxQxcTFmam8QkKkRYukESPsrghAWXltQKpMKSkpysnJcT12nTpUHJ4zfHjZWpDCwqS77670coDK0KGD9P77pjvd66+bwY8AvJ/XBqTY2FhJ0r5TFx/93+vifbGxscrKyiqx/8SJEzp06JDrmNIEBwcrIiKixAM26NHDLCcimV8Pd/7+5vbau+9K0dGerQ2oQL17m7WXJdOKtGCBvfUA+GNeG5CaNGmi2NhYLVmyxLUtNzdXK1eulNPplCQ5nU5lZ2crPT3ddcxXX32loqIidezY0eM1o5wcDrN+2cyZUqdOJUNSaKj5VVm7VurTx74agQryyCPSvfeaRtPbby/bnKUA7GPrYrVHjhzR9u3bXa937NihdevWKTo6WnFxcRo+fLiee+45XXLJJa5h/vXq1XNNBXDppZeqa9euuu+++zR58mQVFBRo6NChuv3221WvXj2bzgrl4nBIffuax6FD0m+/mVajiy6SgoPtrg6oMA6HNGmSmSN0yRLTgLpypVmmBID3sXUepGXLlun6668/bfuAAQM0depUWZalsWPH6q233lJ2drb+9Kc/adKkSWp2ysrqhw4d0tChQzV//nz5+fkpKSlJr732msLL0amXeZAAeEp2tmkw3bLFLE/y9deMQQDOVWX+ftsakLwFAQmAJ+3YIXXsKO3fL91yizR7tulyB6B8KvP322v7IAGAr2rSRJo3z9xF/vRTs24bAO9CQAIAGzid0rRp5vnLL0uTJ9tbD4CSCEgAYJO+faXnnjPPhw41k8sD8A4EJACw0RNPSAMGSIWF0m23SRs32l0RAImABAC2cjikt96Srr1WOnxY6t5dOstKSQA8hIAEADYLCpL+/W/pkkukjAypZ0/p2DG7qwKqNwISAHiBWrWkzz4zq+qsWiX171+2pQoBVA4CEgB4iUsukebMkQIDTYvSk0/aXRFQfRGQAMCLXHONWZ9ZksaNO/kcgGcRkADAy9x1lzRmjHk+eLBZuw2AZxGQAMALPf20dMcd0okTUlKSWbsNgOcQkADACzkc5vZap05STo4Z/r9/v91VAdUHAQkAvFSNGtLcuVLTpmaB2169pOPH7a4KqB4ISADgxWrXNsP/IyOlFSukv/5Vsiy7qwJ8HwEJALxcixbS7NlSQID04YemfxKAykVAAoAq4IYbpDffNM+ffVZ67z176wF8HQEJAKqIv/5Vevxx83zQIOnrr+2tB/BlBCQAqEKef17q00cqKJBuvVX66Se7KwJ8EwEJAKoQPz9p+nTpyiulQ4fM8P+DB+2uCvA9BCQAqGJCQqR586S4ONOC1Lu3lJ9vd1WAbyEgAUAVFBtrhv/XrGn6It13H8P/gYpEQAKAKuryy6VZsyR/f3Pb7YUX7K4I8B0EJACowhITpddfN89Hj5Y++sjeegBfQUACgCrugQekkSPN8wEDzIzbAM4PAQkAfMD48dItt0h5eWbNtl9+sbsioGojIAGAD/D3l2bMkNq2lfbvN8P/s7PtrgqoughIAOAjwsKk+fOl+vWlH388OaEkgPIjIAGAD6lfX1qwwISlJUukBx9k+D9wLghIAOBj2rSRZs40s26/84700kt2VwRUPQQkAPBBPXpIr7xinj/2mDR7tr31AFUNAQkAfNRDD0lDh5rnd90lrV5tbz1AVUJAAgAf9sorUrdu0u+/SzffLGVk2F0RUDUQkADAhwUEmNm1W7aU9u0zw/9zc+2uCvB+BCQA8HE1a5qRbbGx0saNUt++0okTdlcFeDcCEgBUA3FxZo6kkBDp88+lhx9m+D9wNgQkAKgm2reXPvhAcjikSZOkV1+1uyLAexGQAKAaufVWs26bZBa4nT/f3noAb0VAAoBqZtQo6b77zC22O+6QfvjB7ooA70NAAoBqxuGQJk6UbrxROnrUTCq5e7fdVQHehYAEANVQYKD08cdSfLy0Z4+ZI+nIEburArwHAQkAqqmoKDP8v3Ztc5vtzjulwkK7qwK8AwEJAKqxJk2kTz+VgoNNh+1HHrG7IsA7EJAAoJq76ipp+nTzfMIEMwUAUN0RkAAA+stfpOefN8+HDTOTSQLVGQEJACBJSkmRBgyQiopMYNqwwe6KAPsQkAAAkszw/7fekq67Tjp82Cxsm5lpd1WAPQhIAACXoCDp3/+WmjWTdu2SbrlFOnbM7qoAzyMgAQBKiI6WPvtMqlVLWr1auvtuc9sNqE4ISACA01x8sTRnjmlRmj3b9E8CqhMCEgCgVJ07S//3f+b5+PHS22/bWw/gSQQkAMAZ3XWXNHasef7gg9KXX9pbD+ApBCQAwFmNHWuWITlxQurTR9q82e6KgMpHQAIAnJXDYW61XX21lJMj9eghZWXZXRVQuQhIAIA/VKOGNHeu1LSptGOH1KuXdPy43VUBlYeABAAokwsvNMP/o6KktDTpnnsY/g/fRUACAJRZixZm2H9AgPTRRyc7cAO+hoAEACiX6683S5JI0nPPSdOm2VsPUBkISACAchs48OTkkffdJy1fbm89QEUjIAEAzslzz0m33SYVFEi33ipt22Z3RUDFISABAM6Jn5+5vdaxo/T//p/Uvbt08KDdVQEVg4AEADhnISHSvHlSo0bS9u2mJSkvz+6qgPNHQAIAnJeYGDP8PyJC+uYb6d57Jcuyuyrg/BCQAADn7bLLpFmzJH9/6f33Tf8koCrzmYA0ceJENW7cWDVq1FDHjh21atUqu0sCgGrlppukiRPN8zFjpA8/lAoLpWXLzPNly8xroCrwiYD00UcfaeTIkRo7dqzWrl2r1q1bKzExUVksFgQAHjV4sDRqlHnev79Ut66ZN+nOO82fjRubiSYBb+ewrKp/p7hjx47q0KGDXn/9dUlSUVGRGjZsqGHDhunxxx//w/fn5uYqMjJSOTk5ioiIqOxyAcCnFRZKTqe0evXp+xwO8+cnn0i9e3u2Lvieyvz9DqjQT7NBfn6+0tPTlVI8Y5kkPz8/JSQkKC0trdT35OXlKe+UYRY5OTmSzF80AOD8FBZKv/1W+r7if5I/9JBpUfL391xd8D3Fv9uV0dZT5QPSgQMHVFhYqJiYmBLbY2Ji9OOPP5b6ntTUVD3zzDOnbW/YsGGl1AgAKGn3bik62u4q4CsOHjyoyMjICv3MKh+QzkVKSopGjhzpep2dna1GjRopIyOjwv+CUT65ublq2LChdu3axe1Om3EtvAfXwntwLbxLTk6O4uLiFF0JabvKB6QLL7xQ/v7+2rdvX4nt+/btU2xsbKnvCQ4OVnBw8GnbIyMj+R+8l4iIiOBaeAmuhffgWngProV38fOr+DFnVX4UW1BQkNq1a6clS5a4thUVFWnJkiVyOp02VgYAAKqqKt+CJEkjR47UgAED1L59e1155ZWaMGGCjh49qoEDB9pdGgAAqIJ8IiD17dtX+/fv15gxY5SZmak2bdro888/P63j9pkEBwdr7Nixpd52g2dxLbwH18J7cC28B9fCu1Tm9fCJeZAAAAAqUpXvgwQAAFDRCEgAAABuCEgAAABuCEgAAABuqn1Amjhxoho3bqwaNWqoY8eOWrVqld0l+Zyvv/5aN998s+rVqyeHw6G5c+eW2G9ZlsaMGaO6desqJCRECQkJ+umnn0occ+jQIfXr108RERGKiorSoEGDdOTIEQ+ehW9ITU1Vhw4dVLNmTdWpU0e9evXS1q1bSxxz/PhxJScnq1atWgoPD1dSUtJpE7FmZGSoe/fuCg0NVZ06dfToo4/qxIkTnjyVKu+NN95Qq1atXBMOOp1OLVq0yLWf62CfcePGyeFwaPjw4a5tXA/Pefrpp+VwOEo8WrRo4drvqWtRrQPSRx99pJEjR2rs2LFau3atWrdurcTERGVlZdldmk85evSoWrdurYkTJ5a6f/z48Xrttdc0efJkrVy5UmFhYUpMTNTx48ddx/Tr10+bNm3S4sWLtWDBAn399de6//77PXUKPmP58uVKTk7W999/r8WLF6ugoEA33XSTjh496jpmxIgRmj9/vmbNmqXly5drz5496n3KsuuFhYXq3r278vPztWLFCk2bNk1Tp07VmDFj7DilKqtBgwYaN26c0tPTtWbNGt1www3q2bOnNm3aJInrYJfVq1frzTffVKtWrUps53p41mWXXaa9e/e6Ht9++61rn8euhVWNXXnllVZycrLrdWFhoVWvXj0rNTXVxqp8myRrzpw5rtdFRUVWbGys9dJLL7m2ZWdnW8HBwdaHH35oWZZlbd682ZJkrV692nXMokWLLIfDYe3evdtjtfuirKwsS5K1fPlyy7LM331gYKA1a9Ys1zFbtmyxJFlpaWmWZVnWwoULLT8/PyszM9N1zBtvvGFFRERYeXl5nj0BH3PBBRdY77zzDtfBJocPH7YuueQSa/Hixda1115rPfzww5Zl8d+Fp40dO9Zq3bp1qfs8eS2qbQtSfn6+0tPTlZCQ4Nrm5+enhIQEpaWl2VhZ9bJjxw5lZmaWuA6RkZHq2LGj6zqkpaUpKipK7du3dx2TkJAgPz8/rVy50uM1+5KcnBxJci30mJ6eroKCghLXo0WLFoqLiytxPVq2bFliItbExETl5ua6Wj9QPoWFhZo5c6aOHj0qp9PJdbBJcnKyunfvXuLvXeK/Czv89NNPqlevnpo2bap+/fopIyNDkmevhU/MpH0uDhw4oMLCwtNm246JidGPP/5oU1XVT2ZmpiSVeh2K92VmZqpOnTol9gcEBCg6Otp1DMqvqKhIw4cP19VXX63LL79ckvm7DgoKUlRUVIlj3a9HadereB/KbsOGDXI6nTp+/LjCw8M1Z84cxcfHa926dVwHD5s5c6bWrl2r1atXn7aP/y48q2PHjpo6daqaN2+uvXv36plnnlHnzp21ceNGj16LahuQgOouOTlZGzduLHFvH57VvHlzrVu3Tjk5Ofrkk080YMAALV++3O6yqp1du3bp4Ycf1uLFi1WjRg27y6n2unXr5nreqlUrdezYUY0aNdLHH3+skJAQj9VRbW+xXXjhhfL39z+t5/u+ffsUGxtrU1XVT/Hf9dmuQ2xs7Gkd50+cOKFDhw5xrc7R0KFDtWDBAi1dulQNGjRwbY+NjVV+fr6ys7NLHO9+PUq7XsX7UHZBQUG6+OKL1a5dO6Wmpqp169Z69dVXuQ4elp6erqysLF1xxRUKCAhQQECAli9frtdee00BAQGKiYnhetgoKipKzZo10/bt2z3630a1DUhBQUFq166dlixZ4tpWVFSkJUuWyOl02lhZ9dKkSRPFxsaWuA65ublauXKl6zo4nU5lZ2crPT3ddcxXX32loqIidezY0eM1V2WWZWno0KGaM2eOvvrqKzVp0qTE/nbt2ikwMLDE9di6dasyMjJKXI8NGzaUCK2LFy9WRESE4uPjPXMiPqqoqEh5eXlcBw/r0qWLNmzYoHXr1rke7du3V79+/VzPuR72OXLkiH7++WfVrVvXs/9tnFMXcx8xc+ZMKzg42Jo6daq1efNm6/7777eioqJK9HzH+Tt8+LD1ww8/WD/88IMlyXr55ZetH374wdq5c6dlWZY1btw4Kyoqypo3b561fv16q2fPnlaTJk2s33//3fUZXbt2tdq2bWutXLnS+vbbb61LLrnEuuOOO+w6pSpryJAhVmRkpLVs2TJr7969rsexY8dcxzzwwANWXFyc9dVXX1lr1qyxnE6n5XQ6XftPnDhhXX755dZNN91krVu3zvr888+t2rVrWykpKXacUpX1+OOPW8uXL7d27NhhrV+/3nr88ccth8NhffHFF5ZlcR3sduooNsvienjSqFGjrGXLllk7duywvvvuOyshIcG68MILraysLMuyPHctqnVAsizL+te//mXFxcVZQUFB1pVXXml9//33dpfkc5YuXWpJOu0xYMAAy7LMUP+nnnrKiomJsYKDg60uXbpYW7duLfEZBw8etO644w4rPDzcioiIsAYOHGgdPnzYhrOp2kq7DpKsKVOmuI75/fffrQcffNC64IILrNDQUOvWW2+19u7dW+Jzfv31V6tbt25WSEiIdeGFF1qjRo2yCgoKPHw2Vdtf//pXq1GjRlZQUJBVu3Ztq0uXLq5wZFlcB7u5BySuh+f07dvXqlu3rhUUFGTVr1/f6tu3r7V9+3bXfk9dC4dlWdZ5tX0BAAD4mGrbBwkAAOBMCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAAABuCEgAvEbjxo3lcDjkcDhOW4yyKps6darrvIYPH253OQDKgIAEoEIVFhaqU6dO6t27d4ntOTk5atiwoZ588smzvv/ZZ5/V3r17FRkZWZllaurUqYqKiqrU7yjWt29f7d27l4WwgSqEgASgQvn7+2vq1Kn6/PPP9cEHH7i2Dxs2TNHR0Ro7duxZ31+zZk3FxsbK4XBUdqkVorCwUEVFRWc9JiQkRLGxsQoKCvJQVQDOFwEJQIVr1qyZxo0bp2HDhmnv3r2aN2+eZs6cqenTp5c7JBS39CxYsEDNmzdXaGio+vTpo2PHjmnatGlq3LixLrjgAj300EMqLCx0vS8vL0+PPPKI6tevr7CwMHXs2FHLli2TJC1btkwDBw5UTk6O69bX008//YfvO7WeTz/9VPHx8QoODlZGRoaWLVumK6+8UmFhYYqKitLVV1+tnTt3nu9fJQCbBNhdAADfNGzYMM2ZM0d33323NmzYoDFjxqh169bn9FnHjh3Ta6+9ppkzZ+rw4cPq3bu3br31VkVFRWnhwoX65ZdflJSUpKuvvlp9+/aVJA0dOlSbN2/WzJkzVa9ePc2ZM0ddu3bVhg0b1KlTJ02YMEFjxozR1q1bJUnh4eF/+L5LLrnEVc+LL76od955R7Vq1VJ0dLTatGmj++67Tx9++KHy8/O1atWqKtMKBqAUFgBUki1btliSrJYtW1oFBQV/eHyjRo2sV155pcS2KVOmWJKs7du3u7YNHjzYCg0NtQ4fPuzalpiYaA0ePNiyLMvauXOn5e/vb+3evbvEZ3Xp0sVKSUlxfW5kZGSJ/WV9nyRr3bp1rv0HDx60JFnLli076/lde+211sMPP3zWYwB4B1qQAFSad999V6GhodqxY4d+++03NW7c+Jw+JzQ0VBdddJHrdUxMjBo3buxq9SnelpWVJUnasGGDCgsL1axZsxKfk5eXp1q1ap3xe8r6vqCgILVq1cr1Ojo6Wvfcc48SExN14403KiEhQX/5y19Ut27dczpfAPYjIAGoFCtWrNArr7yiL774Qs8995wGDRqkL7/88pxuOwUGBpZ47XA4St1W3Fn6yJEj8vf3V3p6uvz9/Uscd2qoclfW94WEhJx2HlOmTNFDDz2kzz//XB999JFGjx6txYsX66qrrir7iQLwGgQkABXu2LFjuueeezRkyBBdf/31atKkiVq2bKnJkydryJAhlf79bdu2VWFhobKystS5c+dSjwkKCirRqbus7/uj723btq1SUlLkdDo1Y8YMAhJQRTGKDUCFS0lJkWVZGjdunCQzAeQ//vEP/e1vf9Ovv/5a6d/frFkz9evXT/3799fs2bO1Y8cOrVq1Sqmpqfrss89cNR05ckRLlizRgQMHdOzYsTK9rzQ7duxQSkqK0tLStHPnTn3xxRf66aefdOmll1b6uQKoHAQkABVq+fLlmjhxoqZMmaLQ0FDX9sGDB6tTp04aNGiQLMuq9DqmTJmi/v37a9SoUWrevLl69eql1atXKy4uTpLUqVMnPfDAA+rbt69q166t8ePHl+l9pQkNDdWPP/6opKQkNWvWTPfff7+Sk5M1ePDgSj9PAJXDYXni/6kAoAwaN26s4cOH++xyHNddd53atGmjCRMm2F0KgD9ACxIAr/LYY48pPDxcOTk5dpdSYT744AOFh4frm2++sbsUAGVECxIAr7Fz504VFBRIkpo2bSo/P9/4N9zhw4e1b98+SVJUVJQuvPBCmysC8EcISAAAAG58459nAAAAFYiABAAA4IaABAAA4IaABAAA4IaABAAA4IaABAAA4IaABAAA4IaABAAA4Ob/A0uxEJ3HIsPkAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "landscape.set_path(p2)\n", "landscape.plot()" ] }, { "cell_type": "markdown", "id": "17fec32f", "metadata": {}, "source": [] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.12.9" } }, "nbformat": 4, "nbformat_minor": 5 }