Skip to content

Commit fb40ee8

Browse files
committed
review analysis and figures
1 parent 47f1944 commit fb40ee8

16 files changed

+4796
-349
lines changed

data_exploration/dataClassification.ipynb

+7-14
Original file line numberDiff line numberDiff line change
@@ -46,7 +46,7 @@
4646
"pd.options.mode.chained_assignment = None\n",
4747
"\n",
4848
"import networkx as nx\n",
49-
"import community.community_louvain as community_louvain # pip install community"
49+
"# import community.community_louvain as community_louvain # pip install community"
5050
]
5151
},
5252
{
@@ -1473,17 +1473,10 @@
14731473
" ax=ax1,dodge=True, size=5)\n",
14741474
"ax1.set_xlabel('# Parameters')\n",
14751475
"ax1.set_ylabel('$AIC$')\n",
1476-
"ax1.grid(False)p.ma.log2(probs)).sum()\n",
1477-
"\n",
1478-
" return info\n",
1479-
"\n",
1480-
"mutual_info(traindf.mc_mean,traindf.wd_mean)\n",
1481-
"/tmp/ipykernel_136646/267784704.py:13: RuntimeWarning: invalid value encountered in true_divide\n",
1482-
" probs = joint_prob / (np.reshape(probs_x, [-1, 1]) * probs_y)\n",
1483-
"Out[23]:\n",
1484-
"8.442297294007197\n",
1485-
"\n",
1486-
"ax1.scatter(np.arange(0,len(aic_table),1),aic_min,c='red', marker='o', facecolor=None,s=100)\n",
1476+
"ax1.grid(False)\n",
1477+
"\n",
1478+
"\n",
1479+
"# ​ax1.scatter(np.arange(0,len(aic_table),1),aic_min,c='red', marker='o', facecolor=None,s=100)\n",
14871480
"ax1.plot(np.arange(0,len(aic_table),1),aic_min,c='red',alpha=0.2)\n",
14881481
"\n",
14891482
"ax1.set_xticklabels(np.arange(1,len(aic_table)+1,1),rotation=90)\n",
@@ -1847,9 +1840,9 @@
18471840
],
18481841
"metadata": {
18491842
"kernelspec": {
1850-
"display_name": "pymac",
1843+
"display_name": "vscode38kernel",
18511844
"language": "python",
1852-
"name": "pymac"
1845+
"name": "pywork38"
18531846
},
18541847
"language_info": {
18551848
"codemirror_mode": {

data_exploration/dataInterpolation.ipynb

+24-15
Original file line numberDiff line numberDiff line change
@@ -64,19 +64,19 @@
6464
},
6565
{
6666
"cell_type": "code",
67-
"execution_count": 54,
67+
"execution_count": 2,
6868
"metadata": {},
6969
"outputs": [
7070
{
7171
"name": "stdout",
7272
"output_type": "stream",
7373
"text": [
74-
"['OctoberDesert1020_0.h5', 'OctoberDesert1020_1.h5', 'OctoberDesert1020_2.h5', 'OctoberDesert1020_3.h5', 'OctoberDesert1020_4.h5']\n"
74+
"['OctoberDesert1020_0.h5', 'OctoberDesert1020_1.h5', 'OctoberDesert1020_2.h5', 'OctoberDesert1020_3.h5']\n"
7575
]
7676
}
7777
],
7878
"source": [
79-
"folder_path = '../../data/BagsOctober/Run05/h5/'\n",
79+
"folder_path = '../../data/BagsOctober/Run01/h5/'\n",
8080
"\n",
8181
"file_list = [file for file in os.listdir(folder_path) if file.endswith('.h5')]\n",
8282
"sorted_file_list = sorted(file_list) # Sort the file list if needed\n",
@@ -111,7 +111,7 @@
111111
},
112112
{
113113
"cell_type": "code",
114-
"execution_count": 57,
114+
"execution_count": 3,
115115
"metadata": {},
116116
"outputs": [],
117117
"source": [
@@ -130,7 +130,7 @@
130130
},
131131
{
132132
"cell_type": "code",
133-
"execution_count": 58,
133+
"execution_count": 4,
134134
"metadata": {},
135135
"outputs": [],
136136
"source": [
@@ -162,7 +162,7 @@
162162
},
163163
{
164164
"cell_type": "code",
165-
"execution_count": 59,
165+
"execution_count": 5,
166166
"metadata": {
167167
"scrolled": true
168168
},
@@ -180,7 +180,7 @@
180180
},
181181
{
182182
"cell_type": "code",
183-
"execution_count": 60,
183+
"execution_count": 6,
184184
"metadata": {},
185185
"outputs": [],
186186
"source": [
@@ -199,15 +199,15 @@
199199
},
200200
{
201201
"cell_type": "code",
202-
"execution_count": 61,
202+
"execution_count": 7,
203203
"metadata": {},
204204
"outputs": [
205205
{
206206
"name": "stdout",
207207
"output_type": "stream",
208208
"text": [
209-
"dt Odor: 0.005000114440917969\n",
210-
"dt Wind: 0.02499699592590332\n",
209+
"dt Odor: 0.0049953460693359375\n",
210+
"dt Wind: 0.024998903274536133\n",
211211
"dt Gps: 0.20000004768371582\n",
212212
"dt Gps Vel: 0.20000004768371582\n",
213213
"dt Imu: 0.003999948501586914\n"
@@ -242,7 +242,7 @@
242242
},
243243
{
244244
"cell_type": "code",
245-
"execution_count": 62,
245+
"execution_count": 8,
246246
"metadata": {},
247247
"outputs": [],
248248
"source": [
@@ -333,7 +333,7 @@
333333
},
334334
{
335335
"cell_type": "code",
336-
"execution_count": 63,
336+
"execution_count": 9,
337337
"metadata": {},
338338
"outputs": [],
339339
"source": [
@@ -362,13 +362,22 @@
362362
},
363363
{
364364
"cell_type": "code",
365-
"execution_count": 64,
365+
"execution_count": 10,
366366
"metadata": {},
367367
"outputs": [],
368368
"source": [
369369
"# odor.to_hdf('~/Documents/Myfiles/DataAnalysis/data/Sprints/HighRes/Forest/ForestDataJuly_S1.h5',key='odor', mode='w')"
370370
]
371371
},
372+
{
373+
"cell_type": "code",
374+
"execution_count": 12,
375+
"metadata": {},
376+
"outputs": [],
377+
"source": [
378+
"odor['time']=odor.master_time-odor.master_time[0]"
379+
]
380+
},
372381
{
373382
"cell_type": "markdown",
374383
"metadata": {},
@@ -1071,9 +1080,9 @@
10711080
],
10721081
"metadata": {
10731082
"kernelspec": {
1074-
"display_name": "Python 3 (ipykernel)",
1083+
"display_name": "vscode38kernel",
10751084
"language": "python",
1076-
"name": "python3"
1085+
"name": "pywork38"
10771086
},
10781087
"language_info": {
10791088
"codemirror_mode": {

data_exploration/dataKalmanFiltering.ipynb

+6-6
Large diffs are not rendered by default.

data_exploration/figure/Supplemental/klmsupplemental.ipynb

+3-3
Large diffs are not rendered by default.

data_exploration/figure/Supplemental/lookbacktime.ipynb

+3-3
Original file line numberDiff line numberDiff line change
@@ -81,7 +81,7 @@
8181
"outputs": [
8282
{
8383
"data": {
84-
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAClCAYAAACdkQIdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAOKklEQVR4nO3df0jU9x8H8KfntIlzhOHMmxOx1SU1+8MgGDmamj/vOgtL0GJNqj8aDRzEdNSVM9ruj/2RThkNZjqDNRnT2UQl+sM5yjYJrDmVmU1DXaWIm/2463x///Cr6Kzu4/nyfny/zwcInbzveF7wefq5z30+n5efUkqBiEiAztMBiOh/BwuFiMSwUIhIDAuFiMSwUIhIDAuFiMQ4LRSr1YrExEQYDAb09vY+dY3D4UBxcTGSk5Oxfft21NbWigclIu/ntFCSkpJw/vx5vPrqq89c09DQgIGBAbS0tODChQsoKyvDnTt3RIMSkfdzWiibN29GRETEc9c0NjZi9+7d0Ol0CA0NRXJyMpqamsRCEpFveEHiRYaHh6HX62cfR0REYGRkRPPzp6amMDk5iYCAAPj5+UlEIqJFUkrBbrcjODgYOp1rh1dFCmWpJicnn3l8hojca926dQgJCXHpuSKFEhERgaGhIcTFxQFYuMfiTEBAAIDpNxIYGCgRyS1u3ryJjRs3ejrGovlibmZefjabDb29vbPboytECiUtLQ21tbVISUnB+Pg4Ll26hPPnz2t+/szHnMDAQKxYsUIiktv4Wt4Zvpibmd1jKYcdnH5QOnXqFN566y2MjIzg3XffRWZmJgDg4MGDuHHjBgDAbDYjMjISKSkp2LNnD9577z289tprLociIt/kdA/l2LFjOHbs2ILff/nll7P/9vf3R3FxsWwyIvI5PFOWiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMRouqdsf38/CgsLMT4+jpUrV8JqtSI6OnremtHRURQVFWF4eBhPnjzBli1bcOzYMbzwglfcWJ+I3EDTHsqJEyeQm5uL5uZm5ObmwmKxLFjzxRdfYM2aNWhoaMAPP/yA3377DS0tLeKBich7OS2U0dFRdHV1wWg0AgCMRiO6urowNjY2b52fnx8mJycxNTUFm80Gu92O8PDw5UlNRF7J6eeR4eFhhIeHw9/fH8D0DalfeeUVDA8PIzQ0dHbd4cOHceTIEWzduhUPHz5EXl4e4uPjFxXm5s2bi4zveR0dHZ6O4BJfzM3M3k/sAEdTUxMMBgOqqqowOTmJgwcPoqmpCWlpaZpfY+PGjT41x6Sjo2PRpekNfDE3My+/x48fL/mPutOPPBEREfjrr7/gcDgAAA6HA3fv3l0wQL2mpgY7duyATqdDSEgIEhMT0d7evqRwRORbnBbKqlWrEBsbi4sXLwIALl68iNjY2HkfdwAgMjISra2tAKZHGl65cgVr165dhshE5K00fctz8uRJ1NTUIDU1FTU1NbNDveZOD/zoo4/Q0dEBk8mErKwsREdHY8+ePcuXnIi8jqZjKGvWrEFtbe2C38+dHhgVFYXKykq5ZETkc3imLBGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRhNhdLf34+cnBykpqYiJycHt2/ffuq6xsZGmEwmGI1GmEwm3L9/XzIrEXk5TXdsmxn0ZTabUV9fD4vFgurq6nlrbty4gc8//xxVVVUICwvD33//jcDAwGUJTUTeSWzQ17lz55Cfn4+wsDAAQEhIiE+NxCCipXNaKM8b9DVXX18fBgcHkZeXh507d6KiogJKqeVJTUReSWzQl8PhQE9PDyorK2Gz2XDgwAHo9XpkZWVpfg1ODnQfX8zNzN7PaaHMHfTl7+//zEFfer0eaWlpCAwMRGBgIJKSktDZ2bmoQuHkQPfwxdzMvPzcMjlQ66Avo9GItrY2KKVgt9tx9epVrF+/fknhiMi3iA36yszMxKpVq5CRkYGsrCy8/vrryM7OXr7kROR1xAZ96XQ6FBUVoaioSC4dEfkUnilLRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJEJwcCwK1bt7Bp0yZYrVapjETkIzQVyszkwObmZuTm5sJisTx1ncPhwIkTJ5CcnCwakoh8g9jkQAA4e/Ystm3bhujoaPGgROT9nN6k+nmTA+eO0uju7kZbWxuqq6tRUVHhUhgO+nIfX8zNzN5PZHKg3W7H8ePH8cknn8wWjys46Ms9fDE3My8/iUFfIpMD7927h4GBARw6dAgAMDExAaUU/vnnH5SUlCwpIBH5DqeFMndyoNlsfurkQL1ej/b29tnHZWVlePDgAT788MPlSU1EXklsciARkdjkwLmOHDmytFRE5JN4piwRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYTTdY6u/vR2FhIcbHx7Fy5UpYrdYFozLKy8vR2NgInU6HgIAAFBQUICEhYTkyE5GX0lQoM4O+zGYz6uvrYbFYUF1dPW9NXFwc8vPzERQUhO7ubuzduxdtbW148cUXlyU4EXkfsUFfCQkJCAoKAgAYDAYopTA+Pi6fmIi8ltNCed6gr2epq6tDVFQUVq9eLZeUiLyeyKCvua5du4YzZ87gq6++WvRzOTnQfXwxNzN7P5FBXzOuX7+Oo0ePoqKiAjExMYsOw8mB7uGLuZl5+UlMDnT6kWfuoC8ATx30BQCdnZ0oKChAaWkpNmzYsKRQROSbxAZ9FRcX49GjR7BYLDCbzTCbzejp6Vm+5ETkdcQGfX333XdyqYjIJ/FMWSISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISo6lQ+vv7kZOTg9TUVOTk5OD27dsL1jgcDhQXFyM5ORnbt29/6h3eiOh/m6ZCmZkc2NzcjNzcXFgslgVrGhoaMDAwgJaWFly4cAFlZWW4c+eOeGAi8l5O7yk7MzmwsrISwPTkwJKSEoyNjc27831jYyN2794NnU6H0NBQJCcno6mpCQcOHHAaQikFALDZbK6+D495/PixpyO4xBdzM/Pymtn+ZrZHVzgtlOdNDpxbKMPDw9Dr9bOPIyIiMDIyoimE3W4HAPT29i4qvDfwxeFkgG/mZmb3sNvtLs8kF58c6Irg4GCsW7cOAQEB8PPz83Qcov9LSinY7XYEBwe7/BpikwMjIiIwNDSEuLg4AAv3WJ5Hp9MhJCTEhfhEJMnVPZMZYpMD09LSUFtbi6mpKYyNjeHSpUtITU1dUjgi8i1+SsMRmL6+PhQWFmJiYgIvv/wyrFYrYmJicPDgQbz//vt444034HA48PHHH+Pnn38GMD1VMCcnZ9nfABF5D02FQkSkBc+UJSIxLBQiEsNCISIxLBQiEuPWQvHFiwy1ZC4vL0dmZiZMJhN27dqFn376yf1B/0VL7hm3bt3Cpk2bYLVa3RfwKbRmbmxshMlkgtFohMlkwv37990bdA4tmUdHR3Ho0CGYTCakp6fj5MmTePLkifvD/pfVakViYiIMBsMzz053eTtUbrRv3z5VV1enlFKqrq5O7du3b8Ga77//XuXn5yuHw6FGR0dVQkKCGhwcdGfMebRkbm1tVQ8ePFBKKfX777+r+Ph49fDhQ7fm/DctuZVS6smTJ2rv3r3qgw8+UJ9++qk7Iy6gJXNnZ6dKT09Xd+/eVUopNTExoR49euTWnHNpyXzq1KnZ/1ubzaays7PVjz/+6Nacc/3yyy9qaGhIvf3226qnp+epa1zdDt22hzJzkaHRaAQwfZFhV1cXxsbG5q171kWGnqA1c0JCAoKCggAABoMBSimMj4+7O+4srbkB4OzZs9i2bRuio6PdnHI+rZnPnTuH/Px8hIWFAQBCQkKwYsUKt+cFtGf28/PD5OQkpqamYLPZYLfbER4e7onIAIDNmzcvONP931zdDt1WKM+7yPDf61y9yFCa1sxz1dXVISoqCqtXr3ZXzAW05u7u7kZbWxv279/vgZTzac3c19eHwcFB5OXlYefOnaioqFjS1bFLoTXz4cOH0d/fj61bt87+xMfHeyKyZq5uhzwoK+jatWs4c+YMPvvsM09Hccput+P48eMoLi6e3SB8gcPhQE9PDyorK/H111+jtbUV9fX1no71XE1NTTAYDGhra0Nrayt+/fVXj+11Lze3FcrciwwBOL3IcMbw8LDH/tprzQwA169fx9GjR1FeXo6YmBh3R51HS+579+5hYGAAhw4dQmJiIqqqqvDtt9/i+PHjXpsZAPR6PdLS0hAYGIiXXnoJSUlJ6Ozs9ERkzZlramqwY8eO2YtgExMT0d7e7onImrm6HbqtUHzxIkOtmTs7O1FQUIDS0lJs2LDBE1Hn0ZJbr9ejvb0dly9fxuXLl/HOO+9gz549KCkp8drMwPRxira2ttlL7a9evYr169d7IrLmzJGRkWhtbQUwfROjK1euYO3atW7Puxgub4eih4+d+OOPP1R2drZKSUlR2dnZqq+vTyml1IEDB1RnZ6dSavpbB4vFopKSklRSUpL65ptv3BlxAS2Zd+3apbZs2aJ27Ngx+9Pd3e3J2Jpyz1VaWurxb3m0ZHY4HOr06dMqLS1NZWRkqNOnTyuHw+HVmf/880+1f/9+ZTQaVXp6ujp58qSy2+0ey1xSUqISEhJUbGysevPNN1VGRsaCzK5uh7w4kIjE8KAsEYlhoRCRGBYKEYlhoRCRGBYKEYlhoRCRGBYKEYlhoRCRmP8A35KePnUHdxEAAAAASUVORK5CYII=\n",
84+
"image/png": "iVBORw0KGgoAAAANSUhEUgAAARQAAAClCAYAAACdkQIdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAOKklEQVR4nO3df0jU9x8H8KfntIlzhOHMmxOx1SU1+8MgGDmamj/vOgtL0GJNqj8aDRzEdNSVM9ruj/2RThkNZjqDNRnT2UQl+sM5yjYJrDmVmU1DXaWIm/2463x///Cr6Kzu4/nyfny/zwcInbzveF7wefq5z30+n5efUkqBiEiAztMBiOh/BwuFiMSwUIhIDAuFiMSwUIhIDAuFiMQ4LRSr1YrExEQYDAb09vY+dY3D4UBxcTGSk5Oxfft21NbWigclIu/ntFCSkpJw/vx5vPrqq89c09DQgIGBAbS0tODChQsoKyvDnTt3RIMSkfdzWiibN29GRETEc9c0NjZi9+7d0Ol0CA0NRXJyMpqamsRCEpFveEHiRYaHh6HX62cfR0REYGRkRPPzp6amMDk5iYCAAPj5+UlEIqJFUkrBbrcjODgYOp1rh1dFCmWpJicnn3l8hojca926dQgJCXHpuSKFEhERgaGhIcTFxQFYuMfiTEBAAIDpNxIYGCgRyS1u3ryJjRs3ejrGovlibmZefjabDb29vbPboytECiUtLQ21tbVISUnB+Pg4Ll26hPPnz2t+/szHnMDAQKxYsUIiktv4Wt4Zvpibmd1jKYcdnH5QOnXqFN566y2MjIzg3XffRWZmJgDg4MGDuHHjBgDAbDYjMjISKSkp2LNnD9577z289tprLociIt/kdA/l2LFjOHbs2ILff/nll7P/9vf3R3FxsWwyIvI5PFOWiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMSwUIhIDAuFiMRouqdsf38/CgsLMT4+jpUrV8JqtSI6OnremtHRURQVFWF4eBhPnjzBli1bcOzYMbzwglfcWJ+I3EDTHsqJEyeQm5uL5uZm5ObmwmKxLFjzxRdfYM2aNWhoaMAPP/yA3377DS0tLeKBich7OS2U0dFRdHV1wWg0AgCMRiO6urowNjY2b52fnx8mJycxNTUFm80Gu92O8PDw5UlNRF7J6eeR4eFhhIeHw9/fH8D0DalfeeUVDA8PIzQ0dHbd4cOHceTIEWzduhUPHz5EXl4e4uPjFxXm5s2bi4zveR0dHZ6O4BJfzM3M3k/sAEdTUxMMBgOqqqowOTmJgwcPoqmpCWlpaZpfY+PGjT41x6Sjo2PRpekNfDE3My+/x48fL/mPutOPPBEREfjrr7/gcDgAAA6HA3fv3l0wQL2mpgY7duyATqdDSEgIEhMT0d7evqRwRORbnBbKqlWrEBsbi4sXLwIALl68iNjY2HkfdwAgMjISra2tAKZHGl65cgVr165dhshE5K00fctz8uRJ1NTUIDU1FTU1NbNDveZOD/zoo4/Q0dEBk8mErKwsREdHY8+ePcuXnIi8jqZjKGvWrEFtbe2C38+dHhgVFYXKykq5ZETkc3imLBGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRgWChGJYaEQkRhNhdLf34+cnBykpqYiJycHt2/ffuq6xsZGmEwmGI1GmEwm3L9/XzIrEXk5TXdsmxn0ZTabUV9fD4vFgurq6nlrbty4gc8//xxVVVUICwvD33//jcDAwGUJTUTeSWzQ17lz55Cfn4+wsDAAQEhIiE+NxCCipXNaKM8b9DVXX18fBgcHkZeXh507d6KiogJKqeVJTUReSWzQl8PhQE9PDyorK2Gz2XDgwAHo9XpkZWVpfg1ODnQfX8zNzN7PaaHMHfTl7+//zEFfer0eaWlpCAwMRGBgIJKSktDZ2bmoQuHkQPfwxdzMvPzcMjlQ66Avo9GItrY2KKVgt9tx9epVrF+/fknhiMi3iA36yszMxKpVq5CRkYGsrCy8/vrryM7OXr7kROR1xAZ96XQ6FBUVoaioSC4dEfkUnilLRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJYKEQkhoVCRGJEJwcCwK1bt7Bp0yZYrVapjETkIzQVyszkwObmZuTm5sJisTx1ncPhwIkTJ5CcnCwakoh8g9jkQAA4e/Ystm3bhujoaPGgROT9nN6k+nmTA+eO0uju7kZbWxuqq6tRUVHhUhgO+nIfX8zNzN5PZHKg3W7H8ePH8cknn8wWjys46Ms9fDE3My8/iUFfIpMD7927h4GBARw6dAgAMDExAaUU/vnnH5SUlCwpIBH5DqeFMndyoNlsfurkQL1ej/b29tnHZWVlePDgAT788MPlSU1EXklsciARkdjkwLmOHDmytFRE5JN4piwRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYFgoRiWGhEJEYTTdY6u/vR2FhIcbHx7Fy5UpYrdYFozLKy8vR2NgInU6HgIAAFBQUICEhYTkyE5GX0lQoM4O+zGYz6uvrYbFYUF1dPW9NXFwc8vPzERQUhO7ubuzduxdtbW148cUXlyU4EXkfsUFfCQkJCAoKAgAYDAYopTA+Pi6fmIi8ltNCed6gr2epq6tDVFQUVq9eLZeUiLyeyKCvua5du4YzZ87gq6++WvRzOTnQfXwxNzN7P5FBXzOuX7+Oo0ePoqKiAjExMYsOw8mB7uGLuZl5+UlMDnT6kWfuoC8ATx30BQCdnZ0oKChAaWkpNmzYsKRQROSbxAZ9FRcX49GjR7BYLDCbzTCbzejp6Vm+5ETkdcQGfX333XdyqYjIJ/FMWSISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISw0IhIjEsFCISo6lQ+vv7kZOTg9TUVOTk5OD27dsL1jgcDhQXFyM5ORnbt29/6h3eiOh/m6ZCmZkc2NzcjNzcXFgslgVrGhoaMDAwgJaWFly4cAFlZWW4c+eOeGAi8l5O7yk7MzmwsrISwPTkwJKSEoyNjc27831jYyN2794NnU6H0NBQJCcno6mpCQcOHHAaQikFALDZbK6+D495/PixpyO4xBdzM/Pymtn+ZrZHVzgtlOdNDpxbKMPDw9Dr9bOPIyIiMDIyoimE3W4HAPT29i4qvDfwxeFkgG/mZmb3sNvtLs8kF58c6Irg4GCsW7cOAQEB8PPz83Qcov9LSinY7XYEBwe7/BpikwMjIiIwNDSEuLg4AAv3WJ5Hp9MhJCTEhfhEJMnVPZMZYpMD09LSUFtbi6mpKYyNjeHSpUtITU1dUjgi8i1+SsMRmL6+PhQWFmJiYgIvv/wyrFYrYmJicPDgQbz//vt444034HA48PHHH+Pnn38GMD1VMCcnZ9nfABF5D02FQkSkBc+UJSIxLBQiEsNCISIxLBQiEuPWQvHFiwy1ZC4vL0dmZiZMJhN27dqFn376yf1B/0VL7hm3bt3Cpk2bYLVa3RfwKbRmbmxshMlkgtFohMlkwv37990bdA4tmUdHR3Ho0CGYTCakp6fj5MmTePLkifvD/pfVakViYiIMBsMzz053eTtUbrRv3z5VV1enlFKqrq5O7du3b8Ga77//XuXn5yuHw6FGR0dVQkKCGhwcdGfMebRkbm1tVQ8ePFBKKfX777+r+Ph49fDhQ7fm/DctuZVS6smTJ2rv3r3qgw8+UJ9++qk7Iy6gJXNnZ6dKT09Xd+/eVUopNTExoR49euTWnHNpyXzq1KnZ/1ubzaays7PVjz/+6Nacc/3yyy9qaGhIvf3226qnp+epa1zdDt22hzJzkaHRaAQwfZFhV1cXxsbG5q171kWGnqA1c0JCAoKCggAABoMBSimMj4+7O+4srbkB4OzZs9i2bRuio6PdnHI+rZnPnTuH/Px8hIWFAQBCQkKwYsUKt+cFtGf28/PD5OQkpqamYLPZYLfbER4e7onIAIDNmzcvONP931zdDt1WKM+7yPDf61y9yFCa1sxz1dXVISoqCqtXr3ZXzAW05u7u7kZbWxv279/vgZTzac3c19eHwcFB5OXlYefOnaioqFjS1bFLoTXz4cOH0d/fj61bt87+xMfHeyKyZq5uhzwoK+jatWs4c+YMPvvsM09Hccput+P48eMoLi6e3SB8gcPhQE9PDyorK/H111+jtbUV9fX1no71XE1NTTAYDGhra0Nrayt+/fVXj+11Lze3FcrciwwBOL3IcMbw8LDH/tprzQwA169fx9GjR1FeXo6YmBh3R51HS+579+5hYGAAhw4dQmJiIqqqqvDtt9/i+PHjXpsZAPR6PdLS0hAYGIiXXnoJSUlJ6Ozs9ERkzZlramqwY8eO2YtgExMT0d7e7onImrm6HbqtUHzxIkOtmTs7O1FQUIDS0lJs2LDBE1Hn0ZJbr9ejvb0dly9fxuXLl/HOO+9gz549KCkp8drMwPRxira2ttlL7a9evYr169d7IrLmzJGRkWhtbQUwfROjK1euYO3atW7Puxgub4eih4+d+OOPP1R2drZKSUlR2dnZqq+vTyml1IEDB1RnZ6dSavpbB4vFopKSklRSUpL65ptv3BlxAS2Zd+3apbZs2aJ27Ngx+9Pd3e3J2Jpyz1VaWurxb3m0ZHY4HOr06dMqLS1NZWRkqNOnTyuHw+HVmf/880+1f/9+ZTQaVXp6ujp58qSy2+0ey1xSUqISEhJUbGysevPNN1VGRsaCzK5uh7w4kIjE8KAsEYlhoRCRGBYKEYlhoRCRGBYKEYlhoRCRGBYKEYlhoRCRmP8A35KePnUHdxEAAAAASUVORK5CYII=",
8585
"text/plain": [
8686
"<Figure size 595.276x841.89 with 1 Axes>"
8787
]
@@ -151,9 +151,9 @@
151151
],
152152
"metadata": {
153153
"kernelspec": {
154-
"display_name": "Python 3 (ipykernel)",
154+
"display_name": "vscode38kernel",
155155
"language": "python",
156-
"name": "python3"
156+
"name": "pywork38"
157157
},
158158
"language_info": {
159159
"codemirror_mode": {

0 commit comments

Comments
 (0)