2024-heraklion-scientific-p.../notebooks/walker/Step_2_plotting/Step_2_plotting.ipynb

167 lines
158 KiB
Plaintext
Raw Permalink Normal View History

2024-08-27 14:52:41 +02:00
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"ExecuteTime": {
"end_time": "2022-08-18T09:50:40.616906Z",
"start_time": "2022-08-18T11:50:40.181358+02:00"
},
"execution": {
"iopub.execute_input": "2022-08-20T06:27:54.689Z",
"iopub.status.busy": "2022-08-20T06:27:54.685Z",
"iopub.status.idle": "2022-08-20T06:27:55.297Z",
"shell.execute_reply": "2022-08-20T06:27:55.319Z"
},
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"%matplotlib inline\n",
"\n",
"from walker import Walker"
]
},
{
"cell_type": "markdown",
"metadata": {
"pycharm": {
"name": "#%% md\n"
}
},
"source": [
"# Simulate a trajectory"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"execution": {
"iopub.execute_input": "2022-08-20T06:27:55.307Z",
"iopub.status.busy": "2022-08-20T06:27:55.303Z",
"iopub.status.idle": "2022-08-20T06:27:56.162Z",
"shell.execute_reply": "2022-08-20T06:27:56.165Z"
},
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"# Create a Walker instance\n",
"walker = Walker(sigma_i=3, sigma_j=4, size=200, map_type='hills')\n",
"\n",
"# Sample a next step 1000 times\n",
"i, j = 100, 50\n",
"trajectory = []\n",
"for _ in range(1000):\n",
" i, j = walker.sample_next_step(i, j)\n",
" trajectory.append((i, j))"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"execution": {
"iopub.execute_input": "2022-08-20T06:27:56.177Z",
"iopub.status.busy": "2022-08-20T06:27:56.172Z",
"iopub.status.idle": "2022-08-20T06:27:56.321Z",
"shell.execute_reply": "2022-08-20T06:27:56.344Z"
},
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAawAAAGkCAYAAABtmxHBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOxddXwVx9p+Zs85EQLBIYTg7q4tlFIDqlChVC7121u3W7/1Url1pXZbqtSgSh0plEJxdw8EDwSJnd35/tgd3z05obQ0/fb5/ZI5MztrM7s77/O+77xDKKUUIUKECBEixF8c1pG+gBAhQoQIESIZhANWiBAhQoSoEAgHrBAhQoQIUSEQDlghQoQIEaJCIBywQoQIESJEhUA4YIUIESJEiAqBcMAKESJEiBAVAuGAFSJEiBAhKgTCAStEiBAhQlQIhANWiBAhQoSoEKiwA9ZLL72EJk2aIC0tDd26dcPUqVOP9CX9LXDfffeBEKL8ZWVl8e2UUtx3333Izs5Geno6BgwYgCVLlhzBK66Y+Pnnn3HqqaciOzsbhBB89tlnyvZk2rm4uBjXXnstatWqhYyMDJx22mnIzc39E++iYqGsNr/ooouMZ793795KnbDNjywq5ID14Ycf4oYbbsBdd92FefPmoV+/fhg8eDA2btx4pC/tb4F27dohLy+P/y1atIhve/zxx/HUU0/hhRdewKxZs5CVlYUTTjgB+/btO4JXXPFw4MABdOrUCS+88ILv9mTa+YYbbsD48eMxduxYTJs2Dfv378cpp5wC27b/rNuoUCirzQFg0KBByrM/YcIEZXvY5kcYtAKiZ8+e9Morr1TKWrduTW+//fYjdEV/H9x77720U6dOvtscx6FZWVn00Ucf5WVFRUW0atWqdPTo0X/SFf79AICOHz+e55Np5z179tBYLEbHjh3L62zevJlalkW//fbbP+3aKyr0NqeU0pEjR9LTTz89cJ+wzY88KhzDKikpwZw5c3DiiScq5SeeeCKmT59+hK7q74VVq1YhOzsbTZo0wbnnnou1a9cCANatW4etW7cqbZ+amopjjjkmbPvDiGTaec6cOSgtLVXqZGdno3379mFf/A5MnjwZderUQcuWLXH55Zdj+/btfFvY5kceFW7A2rlzJ2zbRt26dZXyunXrYuvWrUfoqv4+6NWrF95++2189913eO2117B161b07dsXu3bt4u0btv0fi2TaeevWrUhJSUH16tUD64QoHwYPHoz33nsPEydOxJNPPolZs2Zh4MCBKC4uBhC2+V8B0SN9AYcKQoiSp5QaZSHKj8GDB/PfHTp0QJ8+fdCsWTOMGTOGG6DDtv9zcCjtHPbFoWP48OH8d/v27dG9e3c0atQIX3/9NYYNGxa4X9jmfx4qHMOqVasWIpGIIdFs377dkEhD/H5kZGSgQ4cOWLVqFfcWDNv+j0Uy7ZyVlYWSkhLk5+cH1gnx+1CvXj00atQIq1atAhC2+V8BFW7ASklJQbdu3fDDDz8o5T/88AP69u17hK7q74vi4mIsW7YM9erVQ5MmTZCVlaW0fUlJCaZMmRK2/WFEMu3crVs3xGIxpU5eXh4WL14c9sVhwq5du7Bp0ybUq1cPQNjmfwkcUZePQ8TYsWNpLBajb7zxBl26dCm94YYbaEZGBl2/fv2RvrQKj5tvvplOnjyZrl27ls6YMYOecsoptEqVKrxtH330UVq1alU6btw4umjRIjpixAhar149WlBQcISvvGJh3759dN68eXTevHkUAH3qqafovHnz6IYNGyilybXzlVdeSXNycuiPP/5I586dSwcOHEg7depE4/H4kbqtvzQStfm+ffvozTffTKdPn07XrVtHJ02aRPv06UPr168ftvlfCBVywKKU0hdffJE2atSIpqSk0K5du9IpU6Yc6Uv6W2D48OG0Xr16NBaL0ezsbDps2DC6ZMkSvt1xHHrvvffSrKwsmpqaSvv3708XLVp0BK+4YmLSpEkUgPE3cuRISmly7VxYWEivueYaWqNGDZqenk5POeUUunHjxiNwNxUDidr84MGD9MQTT6S1a9emsViMNmzYkI4cOdJoz7DNjywIpZQeGW4XIkSIECFCJI8KZ8MKESJEiBD/PxEOWCFChAgRokIgHLBChAgRIkSFQDhghQgRIkSICoFwwAoRIkSIEBUC4YAVIkSIECEqBCrsgFVcXIz77ruPB6YM8ccjbPM/H2Gb//kI2/yviwo7D6ugoABVq1bF3r17kZmZeaQv5/8Fwjb/8xG2+Z+PsM3/ujiiDCtc5j5EiBAhQiSLIzZghcvchwgRIkSI8uCIqQR79eqFrl274uWXX+Zlbdq0wRlnnIFHHnkk4b6O42DFihVo27YtNm3aFNL2PwkFBQVo0KBB2OZ/IsI2//MRtvmfC0op9u3bh+zsbFhWYg51RAaskpISVKpUCR9//DGGDh3Ky6+//nrMnz8fU6ZMUeoXFxcrBtDNmzejbdu2f9r1hggRIkSIPxabNm1CTk5OwjpHZMXh8i5z/8gjj+D+++83yvv2vhUkPR0A4MTckdlJIbBj7uqfNOqmDk+hpDTi5SOE/+appabwFhSl8sKi+iKjAUM/cbyUAnC0MlvLe6llU6Meq8u3sfPxY1Bvu5Sn2nmYfKJfK8v73B/1VlMVbePlWRvJ7eeV6W3O29rrGycibQvoFxqlanmU8rrUokodeHlEWN67BYvybYSwlN1fQId5nUxZPG8AlJU53s56arM8QOLub+JtY3neL3Eo+1i2KON1tLoWfw4o364/M3rK+z6BTMr6Vn++ledfexeciN7/aupEfbaxvNdfRrmXh0XFNqL2pbhGquSVZ1a7VdYHvJy/i0S0sc36i+W9uj79ZrH+0FIS9969UrU8EqeiD3kdqu4Td7zziHJiO941USUFT9l98AcUBth7y/qYtSN/f0XKf0fdSpT1sZc6fLv4psrfT/dmvUsharnc18b31ALskiIse/tBVKlSxbwHDUdkwGJIdgnwO+64AzfddBPPM8pOMirB8gYsmuK2BIkRkJh3PK9x4eXZ3ZKo2pAkCt7YolO1VH85CFDmgMXeK+kDogxe8Bmw2MvCUravTc1t2r7GIOdIDzr/eGnX5vOg6x8xnmoPPM9HxIPN25G1NR98vI8CK4+JbbxujA1C4rgAAK8cUSr28QYmEnXUPEstkVpEK9MGKv2RE98AIn47Fi9z81657ZWzD2Oc8EGMfQjZhw9soOIDGng50euwPuSDm5oSW9TRP2LQBqwgQco9kJoSo6/Fh40LAbwP1JS/TxbEl4UJEFodRPy3w6KiL/m7pw5QvgIHK+P3rA1UerkDIWTwFEpeDFzeKeIEpBRq3VK2zUv5YCS2s4GJ8LpMsNQGMD4wOsrgBcjvsf4+awOZDDbI6EKJV0644CEGLCZo8nfcq2MxoTLC3mPC+4xoA5PxnEiCLetTKj9nxd4xfb79Oo7IgFXeZe5TU1ORmpr6Z11eiBAhQoT4C+KIDFjyMveyDeuHH37A6aefnvRx7JgFmmp5vz0JISapozxmxVWBrI6uppJUGDBUGZ7U78OwqCYQGBomg2ERU31jqXUNlRBLIwSUl6mqP3BpX1NtONRQiZgMy5Rq9Hs11USmSkioB7wy1gc8722XUlHmtXFMPwY1UyaRRxmj8tQpLO/doMXKLQrLEr/lOgy6cM7bgRLOqBzH8VKipt71sLanliV+E1VaZWoV1udcdUgg1JasDutbJhwzadYrp4pKUGNshroIwdDUvjD6GoZK0OhrH9VgsOpPVfOy/mTvGSyY6l12aZb2ECcSyDVGxRUJkipXvyauGmRsifUTbxvK6R1nUHp7cTUcy1Pxvuh14yzvsRjGtCyLM3j2fBOtTxnjYvV8oV2TzrD49UQtRb3vpm7e0VWDUh8bKj/9ufBRHfupBP3IYRCOmErwpptuwoUXXoju3bujT58+ePXVV7Fx40ZceeWVR+qSQoQIESLEXxhHbMAaPnw4du3ahQceeAB5eXlo3749JkyYgEaNGiV9DDuVgKa4Q7UtsSfOrHRGFWTcVxiCauilus7Vz9Drgeo/dEHQoSb
"text/plain": [
"<Figure size 480x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"walker.plot_trajectory(trajectory)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAGBCAYAAAC0IMw4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACMo0lEQVR4nO29e5hUxbnv/63Vt7kPA8NlkGFQGBWVQCJiDCreAEN23GRrYthBQcWYPMlOomdvzyESQaOiZsfEnX0Sf3uTIJGoZx/jMdFEUY8gQlT0qKh4yWgEEUS5zEzPtae7V/3+WKtW11q9LtUzq3suvJ/naR26q6u+Xauq3nV533oZ55yDIAiCIIqANtgCCIIgiJELGRmCIAiiaJCRIQiCIIoGGRmCIAiiaJCRIQiCIIoGGRmCIAiiaJCRIQiCIIoGGRmCIAiiaEQHW8BIQ9d17N+/H9XV1WCMDbYcgiCIAcM5R0dHByZOnAhNK+zahIxMyOzfvx+NjY2DLYMgCCJ09u7di0mTJhX0HTIyIVNdXQ3AOBg1NTWDrIYgCGLgJJNJNDY2WutbIZCRCRlxi6ympoaMDEEQI4r+PAKgB/8EQRBE0SAjQxAEQRQNMjIEQRBE0SAjQxAEQRQNMjIEQRBE0SAjI/HSSy9h0aJFqKurQ2VlJebMmYP777+/NI1nM0CqG9Cz3mU4B/p6gHTK+NurTCZtlNN177p03Wgvm/ZvL50C+nq92ytUe6ZPTTsfStp1Re19CtqzpvZMgPZe40XaC9Du056q9r5e//kFGOMu1R08v/p6jPHspz2dKp12Mb+CtIcMuTCbbNmyBQsXLkQ8HsfXv/511NbW4uGHH8Y3vvEN7N69Gz/84Q+L0zDXjcEhFsxsGojGgVgCYOY5AOfG+329AMwBlOkD4uVARDqEetYoo2dyZWIJIJoAhOuhGNiZVK49LQrEywAtkqtLtCcWnWwfECsDIjGprn5qZykP7T25RSdjtheN+2uPRIFYOSBHIbtqN9sTdem6sRDmaS+zt5fpM9oU2jXRD/3R3mt85qWdc2MhSUvaMwPV3pvTqZljxnaczfYs7WlzzARpjxnt5WmXFsxM2hxXTu09uQUzmzbGZyxRJO26OWaCtKfthjHTZ2iPxOzHOd1boPY+Q7Obdtmoi75yzi837fGy/PmVpz1gbcim89eGIsE49zOhRweZTAYnnngiPvroIzz//PP47Gc/CwDo6OjAGWecgXfffRdvvfUWmpubA+tKJpOora1Fe3u7f5yM20RyIiZCOuV9xheJGQMl25cbjE6YlhtMaZ+zpmgciMSBjDSRnGgRQ5eYcO4N5gyNvGD2V3ssYfztd3YcTRiTKp3KTSRP7RnTcHhoj5cZ/w/SHosbi4OKdtnI9kt71KgrTO3pPu+rQaYZfQXurz2WMLTJi72X9mwmd4KQ1x4z2wvQHo0D0ZiadrFIe2o3T6z8tEei5hhV0Q7zSsJPe9wo46Vdixjthaa9LGB+sfyTGBeU1zW3n1RQ6RHKM888g/fffx//+I//aBkYwIje/9GPfoRMJoP169eH26jvIm2S7jXPPnxuKWTThuHwGkSAOWB7gL5u/8vyTJ/9LNNLd19PgHZxyyTg9ouq9r6e4FsKmZQ54VS0ey3Spva+HjXtfb3q2r0WC6G9L0h7Rl170G0vod33dqNujJcg7eJ2j98Y1c0rHK9FGsjdTu3r9tee6VPULsaon/beYO3ZjHmcVbQH9Lu4reen3Ta/wtAeNL+40edFhIwMjFtlALBgwYK8z8R7zz77bCklEQRBjAjomQyAlpYWAHC9HVZXV4f6+nqrjJNUKoVUKneWk0wmiyOSIAhiGEJXMgDa29sBALW1ta6f19TUWGWcrF27FrW1tdaLdmAmCILIQUZmgKxcuRLt7e3Wa+/evYMtiSAIYshAt8uQu4LxuloRnhVuJBIJJBKJwhtlmvHye1AIZnh8hOL/J+oKwz9exeWRmb/P58GkMprR5JDTrlAX00yHhYCDyIKLqI0FBmgM0EMYNMraFXWptAeEc5wtTYGdOvC2gHC1a5oZx+KnnUGtTwOqAeyu1UWArmSQexbj9tyltbUVhw4dUnJfLghNA8qqcq6PeZ9HAHDDe4RpcB1QmmZMJuFh4jVYmFkX173LiPf1rFEncxkalmEMaE9FO4sY7+tZ4//MT7sevnav9grSnvHRzozjw83FQkm75q9d9Tjr+gD7XRjZEmrXIkZ7vMTaeZB2pjbeLe1+7UnaXbNLMqmvgrRzNe1+/c6YEU+TqCxqrAwZGQDz5s0DADz55JN5n4n3RJlQYWY8SVl1LujLNpFMxISxBrAYjLrdrdc52Zk0keQyQG7QOQcjYNTJdXMisFx5MZFsdUmT3U+75tDOzYlkFMqfMNoAtcuT3U07F9od7YWmndujqp3anScIoj0V7cKYaj595dRuM7Iu2uXjLE5I/LS7HWfniQDz0V5Ivw9YO5Pq8hijNu2SocqbX0Ha3caoZDiEdueJgJ/2oPmlrF0ao1Fz3ZEDV4sEBWPCCMY84YQTsG/fPrzwwguYNWsWAHsw5q5du3D88ccH1jWQoKVAn3YAULr+hXQ254N1WR5YWXCbWsTfZ1+1HkBRu0p7ipRaO4uEdBsRatpVNKmWK3W/HxXaGXxjwFTbUyVebhiXAhjIukbPZABEo1GsW7cOCxcuxFlnnYUlS5agpqYGDz/8MD744APccsstSgZm4KicUQzGOcEQPA8ZgpKIwYQGhDrFvXJxQkbG5Nxzz8W2bduwevVq/Nd//Rf6+vpw8skn48c//jG+8Y1vDLY8giCIYQkZGYk5c+bg8ccfH2wZBEEQIwZ68E8QBEEUDTIyQwWVeAQAyvdTw/QYKbL3iUuDoRQpsGCJ6gm3qpJT8kcfqp01FDu11NoV6wkl3kwdMjJDAZGUKNOXc8t1Q7g5Wn+7wCT3Us8gK9MtUtdhc+N1bQ+GAfRrT7hoaiFp5yrafWIEbG0UoN1rkmpmrA7gEwchtPvE8xSsHWrahWurV1084Dhbru5B2jX/2IyiaHe4mvdHu3D3DqvfxfwaqHZN0u4ZJybFzoS1NqR7jfWmRMaGjMxgIrYI7+3MbfMu/Nud8QaA3YUxb8I4gv/k8nmTh0sutHJciTmAnQF0Vl2OIDJnDIfuot0tliAvnqef2m1xJZJ2ZyyB7jLZ3bQ7DZJXHIQ82a1FoFTaTUMlaxfxI77axXF20e4XiwTAFswra7eVKbV2TU27FjHGpN4P7c75JbSL41yIds1Fu1t8TV6MlLPfnWtDP7Vn00BPR3AW0BAgIzNY6LpxkL3iYsSg8dvexDZhHMF/bnU5B6OtLmnCOAPocoWMNtwWX1ftPr79lvaofRFwrYsVpt31no4cYKqiXQvQznPavc4ILe0+7RWsnXuPh4K0RxS0a5JGn/aUtDtPbtzqYgX0u9cYhbRbA4KPs5925/wK0u43RuUTMD/tPAvXk0pX7drAtad7jZPcIhoa8i4bLDwnZVEbDaka1XpULsdV6gqxn5SrCkuXYoNhdkNoC4bi7RSl9hTLhFWXah+olFOqKswxqtrvYZUp7m0zupIhCIIgigYZGYIgCKJokJEhCIIgigYZmaGO0j3jMO8Hl/p+N2kPv73hrD2k9pSfhZH2YkNGZrDQIsZuqH6BjlrEdFrxOUzCS8gv8ZBwZ/TK01FIXc7t5T3bjMAzt0sh7RWq3SveQGgSZb0L5VxlS6kditr9yljafWJBCtIuXNaDtAfoYhFj3VTq9wDtrADtvn0l6XHN7SK15xdrZRQKv9/D0K4paGca5ZMZsTBmbLddVm3kdpCRXYSt/CiOye7MU+GVyEl2mfRK5GTLU8E9gtYc8SBeyZDk3Ca23C5O7Swc7ZpDu2uwoCMPjK92nnOVdQuS89LunOzK2lGYdq8EWjbtLrEgAOw5bMRxdtFuxbJwuMYP+Wl3G6MF9buPdhGoK2t3Luqydt++krS7JXlz1e4xRoO0u84vp3YpVkzW7tqeonbf+cWMhIllVUCkuE7G5MI82DAGxMuAaMzwWc9m3f3jhZuhHBWfX8g+2blHXXIci5cPvWh
"text/plain": [
"<Figure size 400x400 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"walker.plot_trajectory_hexbin(trajectory)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"hide_input": false,
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.11"
},
"nteract": {
"version": "0.28.0"
},
"toc": {
"nav_menu": {
"height": "12px",
"width": "252px"
},
"navigate_menu": true,
"number_sections": true,
"sideBar": true,
"threshold": 4,
"toc_cell": false,
"toc_section_display": "block",
"toc_window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 1
}