501 lines
1.0 MiB
Plaintext
501 lines
1.0 MiB
Plaintext
|
{
|
|||
|
|
"cells": [
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 3,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"import numpy as np\n",
|
|||
|
|
"import matplotlib.pyplot as plt\n",
|
|||
|
|
"%matplotlib inline\n",
|
|||
|
|
"from scipy.io import loadmat\n",
|
|||
|
|
"import pathlib\n",
|
|||
|
|
"from scipy.interpolate import interp1d"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 5,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"m = loadmat('sargolini2006/all_data/10073-17010302_eeg.mat')"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 8,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"data_path = pathlib.Path('sargolini2006/all_data/')"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 15,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"data = {}\n",
|
|||
|
|
"for fname in data_path.iterdir():\n",
|
|||
|
|
" if not fname.is_file():\n",
|
|||
|
|
" continue\n",
|
|||
|
|
" try:\n",
|
|||
|
|
" action, ftype = fname.stem.split('_')\n",
|
|||
|
|
" except Exception as e:\n",
|
|||
|
|
" print(fname)\n",
|
|||
|
|
" raise e\n",
|
|||
|
|
" if ftype == 'EGF':\n",
|
|||
|
|
" continue\n",
|
|||
|
|
" if action not in data:\n",
|
|||
|
|
" data[action] = {}\n",
|
|||
|
|
" data[action][ftype] = loadmat(fname)\n"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 25,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"action = data['10704-08070402']\n",
|
|||
|
|
"x, y, t = map(action['POS'].get, ['posx', 'posy', 'post'])"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 54,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"import pandas as pd"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 77,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"analog, cells = [], []\n",
|
|||
|
|
"for k, v in data.items():\n",
|
|||
|
|
" if 'POS' not in v:\n",
|
|||
|
|
" continue\n",
|
|||
|
|
" if set(['EEG', 'EG2']).intersection(set(v.keys())) == set():\n",
|
|||
|
|
" continue\n",
|
|||
|
|
" analog.append({\n",
|
|||
|
|
" 'action': k, \n",
|
|||
|
|
" 'eeg': None if 'EEG' not in v else v['EEG']['EEG'], \n",
|
|||
|
|
" 'eeg2': None if 'EG2' not in v else v['EG2']['EEG'],\n",
|
|||
|
|
" 'x': v['POS']['posx'],\n",
|
|||
|
|
" 'y': v['POS']['posy'],\n",
|
|||
|
|
" 't': v['POS']['post']\n",
|
|||
|
|
" })\n",
|
|||
|
|
" for kk, vv in v.items():\n",
|
|||
|
|
" if kk.startswith('T'):\n",
|
|||
|
|
" cells.append({\n",
|
|||
|
|
" 'action': k,\n",
|
|||
|
|
" 'channel': kk[1],\n",
|
|||
|
|
" 'unit': kk[3:],\n",
|
|||
|
|
" 'spikes': vv['cellTS']\n",
|
|||
|
|
" })"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 78,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"analog = pd.DataFrame(analog)\n",
|
|||
|
|
"cells = pd.DataFrame(cells)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 79,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"text/html": [
|
|||
|
|
"<div>\n",
|
|||
|
|
"<style scoped>\n",
|
|||
|
|
" .dataframe tbody tr th:only-of-type {\n",
|
|||
|
|
" vertical-align: middle;\n",
|
|||
|
|
" }\n",
|
|||
|
|
"\n",
|
|||
|
|
" .dataframe tbody tr th {\n",
|
|||
|
|
" vertical-align: top;\n",
|
|||
|
|
" }\n",
|
|||
|
|
"\n",
|
|||
|
|
" .dataframe thead th {\n",
|
|||
|
|
" text-align: right;\n",
|
|||
|
|
" }\n",
|
|||
|
|
"</style>\n",
|
|||
|
|
"<table border=\"1\" class=\"dataframe\">\n",
|
|||
|
|
" <thead>\n",
|
|||
|
|
" <tr style=\"text-align: right;\">\n",
|
|||
|
|
" <th></th>\n",
|
|||
|
|
" <th>action</th>\n",
|
|||
|
|
" <th>channel</th>\n",
|
|||
|
|
" <th>unit</th>\n",
|
|||
|
|
" <th>spikes</th>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" </thead>\n",
|
|||
|
|
" <tbody>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>0</th>\n",
|
|||
|
|
" <td>10704-08070402</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>3</td>\n",
|
|||
|
|
" <td>[[1.9145208333333332], [2.144375], [2.22252083...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>1</th>\n",
|
|||
|
|
" <td>10704-08070402</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>[[0.4719583333333333], [1.0331770833333334], [...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>2</th>\n",
|
|||
|
|
" <td>10704-08070402</td>\n",
|
|||
|
|
" <td>4</td>\n",
|
|||
|
|
" <td>5</td>\n",
|
|||
|
|
" <td>[[0.23117708333333334], [0.2521979166666667], ...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>3</th>\n",
|
|||
|
|
" <td>10704-08070402</td>\n",
|
|||
|
|
" <td>4</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>[[2.9638541666666667], [4.54359375], [5.202333...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>4</th>\n",
|
|||
|
|
" <td>10704-08070402</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>[[4.531625], [4.648020833333334], [5.679010416...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>...</th>\n",
|
|||
|
|
" <td>...</td>\n",
|
|||
|
|
" <td>...</td>\n",
|
|||
|
|
" <td>...</td>\n",
|
|||
|
|
" <td>...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>453</th>\n",
|
|||
|
|
" <td>11084-09030503</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>[[0.98903125], [1.9578645833333332], [2.223593...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>454</th>\n",
|
|||
|
|
" <td>10884-08070405</td>\n",
|
|||
|
|
" <td>4</td>\n",
|
|||
|
|
" <td>2</td>\n",
|
|||
|
|
" <td>[[2.1585], [2.2601458333333335], [2.3705625], ...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>455</th>\n",
|
|||
|
|
" <td>11138-08040501</td>\n",
|
|||
|
|
" <td>7</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>[[3.645072916666667], [3.6803020833333333], [3...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>456</th>\n",
|
|||
|
|
" <td>11016-02020502</td>\n",
|
|||
|
|
" <td>5</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>[[3.2897708333333333], [3.314229166666667], [3...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" <tr>\n",
|
|||
|
|
" <th>457</th>\n",
|
|||
|
|
" <td>11016-02020502</td>\n",
|
|||
|
|
" <td>7</td>\n",
|
|||
|
|
" <td>1</td>\n",
|
|||
|
|
" <td>[[0.44884375], [0.53609375], [0.56032291666666...</td>\n",
|
|||
|
|
" </tr>\n",
|
|||
|
|
" </tbody>\n",
|
|||
|
|
"</table>\n",
|
|||
|
|
"<p>458 rows × 4 columns</p>\n",
|
|||
|
|
"</div>"
|
|||
|
|
],
|
|||
|
|
"text/plain": [
|
|||
|
|
" action channel unit \\\n",
|
|||
|
|
"0 10704-08070402 1 3 \n",
|
|||
|
|
"1 10704-08070402 1 2 \n",
|
|||
|
|
"2 10704-08070402 4 5 \n",
|
|||
|
|
"3 10704-08070402 4 2 \n",
|
|||
|
|
"4 10704-08070402 2 1 \n",
|
|||
|
|
".. ... ... ... \n",
|
|||
|
|
"453 11084-09030503 2 2 \n",
|
|||
|
|
"454 10884-08070405 4 2 \n",
|
|||
|
|
"455 11138-08040501 7 1 \n",
|
|||
|
|
"456 11016-02020502 5 1 \n",
|
|||
|
|
"457 11016-02020502 7 1 \n",
|
|||
|
|
"\n",
|
|||
|
|
" spikes \n",
|
|||
|
|
"0 [[1.9145208333333332], [2.144375], [2.22252083... \n",
|
|||
|
|
"1 [[0.4719583333333333], [1.0331770833333334], [... \n",
|
|||
|
|
"2 [[0.23117708333333334], [0.2521979166666667], ... \n",
|
|||
|
|
"3 [[2.9638541666666667], [4.54359375], [5.202333... \n",
|
|||
|
|
"4 [[4.531625], [4.648020833333334], [5.679010416... \n",
|
|||
|
|
".. ... \n",
|
|||
|
|
"453 [[0.98903125], [1.9578645833333332], [2.223593... \n",
|
|||
|
|
"454 [[2.1585], [2.2601458333333335], [2.3705625], ... \n",
|
|||
|
|
"455 [[3.645072916666667], [3.6803020833333333], [3... \n",
|
|||
|
|
"456 [[3.2897708333333333], [3.314229166666667], [3... \n",
|
|||
|
|
"457 [[0.44884375], [0.53609375], [0.56032291666666... \n",
|
|||
|
|
"\n",
|
|||
|
|
"[458 rows x 4 columns]"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"execution_count": 79,
|
|||
|
|
"metadata": {},
|
|||
|
|
"output_type": "execute_result"
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"source": [
|
|||
|
|
"cells"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 87,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"l = analog.apply(lambda x: len(x.eeg), axis=1)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 90,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"text/plain": [
|
|||
|
|
"150000"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"execution_count": 90,
|
|||
|
|
"metadata": {},
|
|||
|
|
"output_type": "execute_result"
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"source": [
|
|||
|
|
"min(l)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 89,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"text/plain": [
|
|||
|
|
"[<matplotlib.lines.Line2D at 0x7f2b55804780>]"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"execution_count": 89,
|
|||
|
|
"metadata": {},
|
|||
|
|
"output_type": "execute_result"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAD4CAYAAADy46FuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3df4wc533f8ff39k5UHbmWZDKCTNEmYzM1qBSlbVpWUTdg7QSS1SCyAMOhmtqC61QqSqF267aRjAJS0BBIAstqDNsEpJCRXBhiBNuI2UBtYCtGlPxhyZTM2iJV2YQlRyQoibFk2a1iaWf22z/mmdu5uZnd493s3Dyczws43O2ze8PZ2ePznef5Pj/M3REREakzt94nICIi3aZAISIiEylQiIjIRAoUIiIykQKFiIhMNL/eJ9C0jRs3+tatW9f7NEREovLoo4/+rbtvqnrunAsUW7du5ciRI+t9GiIiUTGzH9Y9p64nERGZSIFCREQmUqAQEZGJFChERGQiBQoREZloaqAws4Nm9ryZPV4ou93MTpnZ0fB1TeG5W83shJk9aWZXhbItZvYNMztuZsfM7GOrPZaIiLRrJcNj7wE+C3yhVH6nu3+qWGBmO4A9wOXAG4Cvm9kvAgnwCXd/zMxeCzxqZl9z9+Nneyx3T8/qHYqIyJpMbVG4+0PACys83rXAIXd/xd2fAk4AV7j7aXd/LBzvp8ATwObVHGuF5yGy7n76syFfPXpqvU+jk9ydrzx2kr97Na77viNPv8CTz/50vU+jdWvJUdxsZt8JXVMXhbLNwDOF15ykFBDMbCvwNuDhtR6rcMwbzeyImR05c+bMqt+QSJP+5+PP8rFDR3n2pZ+t96l0ztM/epn/cP//5utPPLfep3JW/sufPs4fPvi99T6N1q02UOwH3gzsBE4Dd6zkl8zsAuDLwMfd/SdrOVaRu9/l7rvcfdemTZUz0EVa90oyAuDV8F3GXkmylkRs1+bVZBTdOTdhVYHC3Z9z99TdR8DdjLuETgFbCi+9LJRhZgtkQeKL7v6VtRxLJAZpmlUoyah/Fcs0SZrtrJmO4tphMxk5SWTn3IRVBQozu7Tw8DogHxF1GNhjZhvMbBuwHXjEzAw4ADzh7p9ey7FWc74i6yGvUGKrDNuQX5PYKt105L38PKeOejKz+4DdwEYzOwncBuw2s52AA08DNwG4+zEzux84TjbSaa+7p2b2buBDwHfN7Gg49Cfd/QHgD87mWI28a5EWJJFWhm3IW1lpZK2tYTpimMZ1zk2YGijc/fqK4gMTXr8P2Fcq+2vAal7/obM5lkgsFu+aUwWKsvyaDCO7Nn1tUWhmtsiM5JWhchTLpZF2yylHISKNShe7V/pXsUwTa7ecWhQi0qhYK8M2jFsUcbW2ktGol12JChQiMxJr90obYg2ialGISKOGiwnbuO6a25Dkc0wiujt3d4apM4ysFdQEBQqRGVGOol6MLYr8VPv4eSpQiMxIjJVhW2LMUeSj12JqBTVFgUJkRpSjqBdjEO3z56lAITIjMVaGbYmxW67Pn6cChciMjBO28XSvtGUY4czsPk+gVKAQmZE+34FOE3OOIo0ouDVFgUJkRvrcpz1NjEE01hVvm6BAITIjMVaGbYkyRxHpHhpNUKAQmZHFFoVyFMvEGETHLYr+fZ4KFCIzMlzc4S6eyrAti4nhiIJoHiBGDqOefaYKFCIz0uc+7Wli3P2v+Dn27TNVoBCZkRgrw7bkOYqYKtzijOy+faYKFCIzoh3u6sUYRNMlLYp4usyaoEAhMiNJhHMF2pKm8QXRYusnpgDXBAUKkRlJlMyuFWOLoph4j2lGeRMUKERmRMnsennXTUx7O6RqUYhI0xLlKGrFOGs9UY5CRJoW43pGbUkizFGoRSEijUsiHALalvhbFPGcdxMUKERmpM9rA00zjHA5jGIyO6aWUBMUKERmJL/r7NsImZWIcsKdchQi0jTlKOopRxEXBQqRGdE8inrKUcRFgUJkRmKsDNsS5zLjo8LP8Zx3ExQoRGZkGGFl2JbFbUUj6pYr5pqGES2P3gQFCpEZGS8K2K9KZSWUo4jL1EBhZgfN7Hkze7xQdruZnTKzo+HrmsJzt5rZCTN70syuCmVbzOwbZnbczI6Z2ccKr7/YzL5mZt8P3y8K5WZmnwnH+o6Zvb3Zty4yO+6uJTwmiPHaKEcx2T3A1RXld7r7zvD1AICZ7QD2AJeH3/m8mQ2ABPiEu+8ArgT2htcC3AI86O7bgQfDY4D3AdvD143A/lW8P5F10ee7z5WIcVHAJTmKiFpCTZgaKNz9IeCFFR7vWuCQu7/i7k8BJ4Ar3P20uz8WjvdT4Algc+F37g0/3wu8v1D+Bc98E7jQzC5d4XmIrKs+332uRIz7T/f5M11LjuLm0CV0MO8uIqv8nym85iTjgACAmW0F3gY8HIoucffT4edngUtWeqzCMW80syNmduTMmTOrfDsizVGLYrK8oo1p/+lUO9ydtf3Am4GdwGngjpX8kpldAHwZ+Li7/6T8vLs7cNafgLvf5e673H3Xpk2bzvbXRRpXTNL27e5zJYoJ/tTjuD7DJS2KeFpCTVhVoHD359w9dfcRcDdwRXjqFLCl8NLLQhlmtkAWJL7o7l8pvOa5vEspfH9+2rFEuq5YkWjU03JLthWNpL8/HWmtp7NSyhVcB+Qjog4De8xsg5ltI0tEP2JmBhwAnnD3T5cOdxi4Ifx8A/DVQvmHw+inK4GXCl1UIp2mrqfJYlw3qc9boc5Pe4GZ3QfsBjaa2UngNmC3me0k6yZ6GrgJwN2Pmdn9wHGykU573T01s3cDHwK+a2ZHw6E/GUZL/R5wv5l9FPgh8MHw/APANWQJ8ZeBj6z97Yq0o8+Jz5WIMZCmPe5OnBoo3P36iuIDE16/D9hXKvtrwGpe/yPgvRXlDuyddn4iXRRjRdimJV1zkVyfpS2KOFpBTdHMbJEZyJd4GMxZNF0rbUpSZzCX3TvGEkiT0WjxnPu2dLwChcgM5JXf+fNzvUt8rkQycs6fz6qfWNZNSgvnHEtwa4oChcgM5N0UGxYG0XSttCkdORsWBos/xyBJx+fct89UgUJkBvLKb8P8XDQVYZuS0YgN4e48lko3HfniOStHISJrlhQCheZRLLe00o0jUCQj57zIgltTFChEZiAPDucvDKKpCNvi7gxT5/y8GyeSHE4yGrEwmGPO4jnnpihQiMxAMUcxVKBYIr8c4/7+OFpcSerMzxnzgzm1KERk7ZSjqJcHhhhzFPMDY37OlKMQkbVLSoHCI1n4rg3FIFp83HXJyBnMzYW5MXGcc1MUKERmIF28a45rCGgbxkE0rhxFOgpdT3PWu89TgUJkBvKZuxsW4upeaUNSujaxVLrDNJuZPZib08xsEVm78czsfk7QmiTPUeTXZhhJf//SFkUc59wUBQqRGcgDw/n5XXPP7kAnSSO9NsnImR/MMT9QjkJEGlDOUcQyBLQNi11PkbW2lKMQkUbF2g/fhsVRT5Fdm2zUk2nUk4g0IykNAe1bxTLJ8nkUcbS2knQUWhT9W5ZFgUJkBsY5Cg2PLYv12qSFFkUs59wUBQqRGUjztZ4i23OhDXm3XH5tYplHkYychcEcC0pmi0gTims9QTx3zW1IS9cmlkpXLQoRaVR5mYpYKsM2lPM3scxJSEbFHEW/Pk8FCpEZKC9T0bc70EnyRHBsw2Pzfb77uA+6AoXIDCz2w2sJj2WWTbiL5Nok+TwK5ShEpAnpaIQZ4x3RlMxeVB71FMu6SWk+M1s5ChFpQn73OZizxceSWb7MeBxBNM9RDJSjEJEm5CNk5ufi6l5pQx40Y9t/evyZqkUhIg0Yps582OQG4qkM25B3w+X7T8dQ6eb7fM/PGYOBRbPibVMUKERmIB2NFu8+88eSyYPm4lDTCAJFfoqDOeUoRKQh2SzebIQMxJOwbUNeyS4u2R1Boj8fDpvtma0chYg0QDmKesUWRSwrsaZLWkFqUYhIA7JRT8pRVMm74WJKDOef3yDkKPr2eU4NFGZ20MyeN7PHC2W3m9kpMzsavq4pPHe
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"source": [
|
|||
|
|
"plt.plot(l)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 35,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"text/plain": [
|
|||
|
|
"((30000, 1), (30000, 1), (30000, 1))"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"execution_count": 35,
|
|||
|
|
"metadata": {},
|
|||
|
|
"output_type": "execute_result"
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"source": [
|
|||
|
|
"t.shape, x.shape, y.shape"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 36,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"x, y, t = x.ravel(), y.ravel(), t.ravel()"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 37,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": [
|
|||
|
|
"sx, sy = interp1d(t, x), interp1d(t, y)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": 49,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hVVdbGf/uc229IQgktlBBC7wRCR4qiNMGCBbuOztjHMpaxjnUcO/Y6lm8URBEBARvSCSX0FjoJIQQCJKTcfvb3x7k3uTUJCFPwvs9zH8g5u5199n732mutvY6QUhJHHHHEEcfZCeU/3YA44ogjjjjOHOIkH0ccccRxFiNO8nHEEUccZzHiJB9HHHHEcRYjTvJxxBFHHGcx4iQfRxxxxHEWI07yccQRRxxnMeIkH8fvBkKI8qCfJoRwBP19lRCiqxDiByFEsRAi6gESIcRkIcQaf55CIcQ8IcRg/70rhBC5QohSIcRhIcSnQojEf+9TxhFHKOIkH8fvBlLKhMAPyAPGB137F+ABvgJuipZfCHEv8BrwHNAEaAW8DUzwJ1kGDJJSJgHpgAF45kw+Uxxx1AbDf7oBccTx3wIpZS6QK4TICL8nhEgCngJukFLOCLo12/9DSpkfls0HRJQVRxz/TsRJPo446oYBgAX4tqZEftXN90AiUAlcdOabFkccsREn+TjiqBsaAsVSSm9NiaSUS4EkIUQqcDOw79/QtjjiiIm4Tj6OOOqGo0AjIUSdBCMpZQEwH5h6RlsVRxy1IE7yccRRN6wAXMDEk8hjANqemebEEUfdECf5OOLwQ+iwACb/3xYhhBlASlkKPA68JYSYKISwCSGMQojRQoh/+NNfJYRo5f9/a+BZ4Jf/yMPEEYcfcZKPI45qtAYcwBb/3w4gN3BTSvkycC/wKHAEyAfuAGb6k3QGlgshKtDdKXPR9fJxxPEfg4h/NCSOOOKI4+xFXJKPI4444jiLESf5OOKII46zGHGSjyOOOOI4ixEn+TjiiCOOsxj/VSdeGzVqJNPS0v7TzYgjjjji+J9CTk5OsZQyJdq9/yqST0tLY82aNf/pZsQRRxxx/E9BCLE/1r24uiaOOOKI4yxGnOTjiCOOOM5ixEk+jjjiiOMsRpzk44gjjjjOYsRJPo444ojjLEac5OOII444zmLEST6OOOKI4yxGnOTjiCOOOM5i/FcdhoojjjMKnw/mzYOcHP3/qgqZmTB6tP7/f2c75syB6dP1vydNgnHj9P/Pmwfr1kH37vrfGzdCr14wciQ8/zwsWwaDBsGjj4LJBG43PPNM9fWHH4ZffoE1a2D7djh8GAYPrk5fU5sCdffq9e/vk1PB/2Kb/wP4r4on36dPH/k/f+L13z3w/tcGek3tDb83ahT8+COsXAmLFkFeHtjtkJICQ4fWTlzh9Z5/PmRnQ0VF9XWLBZo3hxdf1P/++ms4eFC/dumletsCRBvetyfT9wEyXroUdu/WnyUY3bpBeTkUFoLTqZcjJWgaGI16XZpG1WytXx+xaxekpurpAzAYkEYjOBwhxXvrJXLizXdJnj8HVWrg9UJuLjRoALfdBs8+C9u2gcej192qFWzZovdvtGcMfvbwBSnw3qLdi9WHdVl4w+ucMgVWrdLfp90O/frBDz9U5fNpElUR0d/HWQYhRI6Usk/Ue3GS9+N0kGWASFaujDnwqhAugZ0MYZ1KfdGecdQo/e+pU2HzZhACLrwQHn8cl6KSf6QM1+zvMW3aQHF6J3yaRuO9uXi6dsNqVEnauQ21dy/M48dhsRgRopYJFd5es7maYIWA++/XCdbp1J/FZNL/H0ZYVahfHw4d0tOFP1uY5Ovr1h159dUY3K6IYiSgEUN3qSggJdJmw53Zl8PTZ1Hm1Th4tIKuN0yi4db1GJwOvBYrh9p3J2fiNTTcvZ2CNh3Z0GUAZT5JZXklU+4ZQ0JlWVWx4T0lY1wPvi/C/vYhUJER16OVIWu4F6tN2GyQno7cvh28XjSTmYrMLPL+9Q3trr8M06pshNOpvzvQFyWbDWk2IzwefdHyj0OpaUiTGVfjJmz9y5OUn3sByVYjHa67FPOaVYjKyup6ExL0cTx3buhiMWUKruXZen+rBgw+L2r1skelycId4x9gYXomw/bk0KVoN9sat+Hmoem0yNtJw8H9sFw4rm5z41SFpdMxr08BcZKvDeXl0LatvrUF/aUMHgzz5uGb8z2uqdNwejT2jxxLbuZQih1eAIxI0lYvpvGurZR27IrB52PAY3eiBhGJz2bn4Dsfo144nmSbEatR1SdA06Zw/Hh1G/yEVYnC4RMuDpe5OFzm5ITDSz2LAaMqEEKgCIEiQBGCRgvm0/neW1BdofXtee09DEaVels34+3ZE+38CzAbFJImjkXNzganE2EwgKIgPR59cgahxGjjwTF38cjCT2hSfhSjz4MUAoEAqSHQCUECDqOF9c3a81nWBHodz2NLk3R+aNkLi8VEfbuJZJuJxjaVMfnr6bpkLukL56F6vVV16QQrEAKErBthBd+f3XEwq7JG8ef575FcUozqcSNtNoTTBb7qejQhUMLKP1lUGC3ceeEDLMjIYsSuVbwx6x/YPdVStA/wKSqq5sNptLCueQeuvewp7l72JXctn/qb6g5HrEWhJpI/2fqj1eFF8E7/S7h15TcYYnBHTXVJwKUaWdOiM59kjuf1WS9i90YuvA6ThaNNW5JSdACjywFC1Do+fMBrgyfTL38LvQ7mYvG4QAgkEiElTqOFba068fXT7/Ggkk/9HVuqhZ0xY3Tho7xc3901bcrxLj1wbdqKy+tjX7M2ZPc9j029BiNUAzaTSs+W9endKpnuLZIxGxTKyyupl9YSEWVeo6pndMf9uyf5nUVl7D5STov6Ntw+DafHR6LFSOuGNupJL9Juh7ABpAGV7Tpi3ZmL4h/uGoLsVt245vKnAfjsq8erBpPTYEKREovPHVFOdstuZLfqxuamGSxMz4w66XXCGsK9Y+9h6L71dCnazZYmbVmYnommqCiaj+G7VzNu+xJAMLfDQB5d8BGtSosi2220YPB5MWreKrL5pPdY3p35fMjErIkQToYUfULBo6gYfV6cRjPrmnfgxosf5Y7s6Zy7cyXNy4qxeFyYNZ1060o2dSF5/SdQapFoT4XkwuFD8NrgyWxp0pabVs+kf95GgqdpeB2BReGm1TMZmLfxtJJ8LJyO56ypLAmUmmwkuSvrvOuIhgqDmY3N2tEvf3PUHZQEPIqKSfPVuVwvMK/9AC7YkY2R6LxWYbSwL7kZaSWFWL0ufeeR3hZ270atrFbjxWLFknr1efeK+/mpbRZ7jrtQNF/VrqFdcR7jty+JbOMjj+hqwrruuE8Bv2uS9/g02j0yL+b9Kd+9EPXFBAgkfAA6FQO3XfRXFKnx5sy/Y9FCpdJY23DQJZic1E6AZGDepsjtt1AoM9kw+TxYvO4qwrz+0if4ZPqT9M/fhCqrFxyQhA+RaNKXWzFwKKE+LU8cqRMBnCxRhKd3qkZUzYdBalXXToV4fgth1UbytalHoiGEIDzOqh1NLPiAV4dcjdHnPe2SfDQEdkVO1YjN5666frr6MHCttnLr8t4kgGoAn7dG9VJd3f/CWSxW/Zo/bfC8cSsGDJq3TnVV1TNsGMenz8Q4eCDW3TsRXg/EEDYqO3bBdmA/ory8+obRqNtTWreGu+6CCRN+E+HXRPJnvQulUVV4emLXmPd7FO6IeS/aQDFpXroU7eaxnz+okkxrggj6WXwe+udv5pg1MWo6g9RIdpVj97pQkdg9TnodzOXOFV/Rp2AbBr9kLQAVGfXlRSMeo+al5Ykjtbb1dMHs82AIUuucDjXByaK2vB5FpcJgjimxBbdFQyf4/cnNSDtWgN3jRKFaZRWrDAEM2rMWm/SgCRG08wj9nQ4E2nnbhffx5/H3k5fYBK2GtLX9Au2P9kx1eafBZWmIqM8pYhB8ANpJjABB3domiCQ9Qx3mcXg9cuFCfK1bk5C7FYPXgwoh9oFgHD14BBlM8ID0eJD79ukOBZdcAuedp9sFzgDOekk+AK9PY/u
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hUVdrAf+dOTe8VAiGBJCT03gm9C/besay9rd+669rW9VN37Yq9rA1ERVAEVEqkl4QAIZUkQHohkD597vfHTYZMZhKC4O5+Or/nyZOZueeee257z3ve9z3vEbIs48GDBw8efptI/+kGePDgwYOHXw+PkPfgwYOH3zAeIe/BgwcPv2E8Qt6DBw8efsN4hLwHDx48/IbxCHkPHjx4+A3jEfIePHjw8BvGI+Q9/G4QQjR3+LMLIQwdvl8thBgkhPhBCHFCCOF2AokQ4iohRHrbPpVCiPVCiElt264XQmQIIRqFEGVCiOeFEOp/71l68OCMR8h7+N0gy7Jv+x9QAizq8NtngAVYCdzsbn8hxAPAy8AzQATQB1gGLG4r4g3cB4QCY4EZwEO/4il58HBGhGfGq4ffI0KIY8BSWZY3utnWHzgiy7Lo8FsAUA7cKMvylz08xgPANFmWF52fVnvwcPZ4NHkPHnrGeEAPfHMW+0wBsn+d5njw0DM89kIPHnpGCHBClmVrTwoLIW4CRgFLf9VWefBwBjyavAcPPaMOCO2JI1UIsQT4X2CeLMsnfvWWefDQDR4h78FDz9gFmIAl3RUSQswF3kVx6mb9OxrmwUN3eIS8Bw9tCAU9oG37rhdC6ABkWW4AHgPeEEIsEUJ4CyE0Qoh5Qojn28pPBz4DLpZlee9/6DQ8eHDCY5P34OE0fYGjHb4bgONALIAsyy8IIaqAR1GEeROQAfy9rfxfgQBgnRCOwJxtsizP+9Vb7sFDF3hCKD148ODhN4zHXOPBgwcPv2E8Qt6DBw8efsN4hLwHDx48/IbxCHkPHjx4+A3zXxVdExoaKsfGxv6nm+HBgwcP/6/IyMg4IctymLtt/1VCPjY2lvT09P90Mzx48ODh/xVCiONdbfOYazx48ODhN4xHyHvw4MHDbxiPkPfgwYOH3zAeIe/BgwcPv2E8Qt6DBw8efsN4hLwHDx48/IbxCHkPHjx4+A3jEfIePHjw8Bvmv2oylAcPv0tsNli/HjIzYfhwmDcPVKrzU9+QIcpvhw6dn7r/P3C+r+f/czxC3sN/D2d6Of+bXt7z1RabDebMgT17oKUFfHxg7Fj44Ycu6zNb7WjVXQzCbTbss2fDnj2IlhZkSSknZBl0OoiKQrz0EixceHbtdXe+0PU16O76/NJr13m/2bPhxx9d29T5eo4ZA/fc8/vq6DrgEfIeuqTRaKGkrpUmoxWDxYrBbEenlvDTq/HVq/HXa5TPOjVq1Tla/s4k7Mxm5WXNzVU+63TQqxe8+GKPBJYsy6zaX85ne47TZLTSarah00gEe2sJ0quYULCX+LIjGAYNIeCiCxjaNwQvrVKnxWyhefV36LIOIkaMoGnKNLwWL8IrMwOVoRW7RktraDjp9z5G6fhptNqhtdVE5K4thOTnUBQzgD1JY2m2yNQ0GbFYZYxWG3a7zPSivTyzdQfeFqPS0OZmDNt28PgNz3DCYCWluojsiHjS4kZil84smKYX7uW1bTvxaatP2O2nNxqNyEePYrv8cuzjJyD9+CNCpaK83kBhTTONRgsj+gQRE+ztfGssVozTZ6LP2IdkaEXWaCAxCSksFPbtc7lfsiwj5s49fS91OmxRUVQ9/gxlY6eScNPl+B/cj2RoBZ0OER3t/j7abLB2LXz5JdjtGDIPYi8qxstiwqLTI2s0SFYrapMBs1bP4d5JvDV0Aa9s3eE4f5qbkdPSYOs2hNWiPDeJiXDRRTBy5O9C4P9XrQw1atQo+b8xd43VZudki5maJhO1zSZqm0w0GiyoJIFaJaFp/68SqCUJtUo4PmvVEv56DYHeGgK8NHhrVXRYGs4tsizTaLBS2WigssGI1SajUQkk0f4Hou2/JJ3+rX27EKBWCXRqFVq1hK7tT6uW0Kokl+NbbXbK6w18e6CC5XtLqGgwnvU18tWpCfDS4CXJjMzdzaCqQvTYiT5ZhV6rpnruBZxKnYVaq0YIgcliw2ixY7TYMFnt9N6xiQv/8Ud0xlZHnS0aPXdf8DBpcSNZ+9F9DKw9SseWy4BJreVI/8F8//xHDIz2Z3TuHgLzDyONGIFqwXwEIK9fR866rbxS5+siLLVmAz+9fxe9mmoRsh2jSktm74Fcd9lT2CUVkt3GxysfY3hFPnqLCaNGx/HAKPqeqsDHanJui0pDeu9kbrjkcT766kmnfTKjEx11duSe7Z9x/47lLudVEhBBaGvDGffvzN07lnPf9s9R0f173X5tN/cf0205aOs41jzncr52hNNx2usEeO3b508LWk5fn6Lg3sTWV7rdlt47mesue4peob5MiPHjgT9fQ/ixI20lFDpfp47fjSoNGwaMZ3He1m7LOb6r1RAaCsuWwYIFTqOCE5On893havYdO0mjwYoQoFUp77ZakkiO9md4n0CG9g7ERy2wrv0e1cEDiBEj/iMdhxAiQ5blUe62/W40+byqRpqNVoJ8tJgsdkxWG6G+OqIDvVBJArtd5mBZPaWnDPjqVOjUKtSSwGi1c9+KTE61Wn71Nkp2G6nFGWetvf2ax+hp+WaTlVaDSRGI5Xl4dxAIALYd69kdM5jrLv+by/5qq5lPVr6JpoOAB9BbTCTXFAPQv66Uzl2jAPRWM3FF2RR+9AWT0r8luDwXrc2CWaVhT3QSQgiGVRYwyGLitU7C0re5ngNvXIsK2VG3j83M8LJcUosz2Nx/DKnFGQyvyHcIJR+LkcQTxxGy3bUtNgujynL4pL1TsFkc+wyvyHfU2RFVF0pWVFMdWrv1jPuDInwuHdWbtYcqyY6Ix6jROQlRd+gtRm5KXw1wxmcgpboIfaf7KQCpU0fSfr+ELKO3uJbX2ywMqCtF3XZenbe1n2MaI7nun0sJ79SpnwmdzUJq8ZmVREedVityVRVcfAkmH19UVgtqkxHZ25vC8AH87ZIniQnzI8RHi4xiJrPaZExWG99nVQKglm189IWzElCeOIQ++7ah12vPovW/Hr8LIV9S18rcl7f9p5vRLe40xp5qbz2pO7U4g0FVhcwp2EXsqQr0VvMZj3G2bXIIxE4CAUAty4wqz3URVGqrmX2vX0egqdllH6Nai8Zm5eZ9q10EQ0e8LEYW5G5lXGkW6jahqbdZmFCahVVIaNsEso/FyJjSbJ49/A1rZl3Op89dA+AiSLxsZhbmbiUtbqQi4DoJLKmTgO+IzmZhbOlh1zotRkedHa+dTbg3c2k6na/eYmLp3m+4ed9q0nsn89r4y7CqFSFittn5bE8Jkt2GJNs54R2IqukEOjfCtB0VMOn4ISZV5SMHBmJJnUbz4gupmzqberOdvKomtuTVsDmvhtzwfghlTNTleQMYNTpywuOUNun0eJkMLmXUdiuypAK7zWWb3mIkpboQgAFuOvXOuOv0A8yt7op2X4dsR9fc6KhPtLQw7Phh7tv+GRaVhlMDBjKqTxAj6kvxGz+awIsvoMFsJ/NoHfpnnmZkabZThxybk8GG4TM5ERPHqAtnMmjpFXyXXc3T3+fyt8WDmDso8qzaeK78LoR8ZICemyb244MdR39xHb+2lu1OY+xOe+uujYOqClHJdmxCkBMRzw3p3zK8sgAvixHB6ZfjTMfoqk1vh9Zwy8kol/LuBGJHtDYLKdWFTse6e9dKAk3NLsNpGTBLapbuW+1od1cIYFpxhotWLABNJ4GssVtZ8sOnzNn0haOMu/rm5+8gvOUU/xq+AIQA2dVkIHf63v65qzoXF+xgdkMxz8y8lRXhg7FLKrIj4rAJ4eicukJCZnxpFgATSg5x+6F1ZC+5BtvOXeyMTOSNsRfzwaqn2zpkI/a2VnS+rnaUuGnH7yYTVFej+WIFAStX0jJ0DMffXoEAJsSHMLJvEFH648irXNtkR2BRqdHYrFh
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hUVfrHP3d6Jr33QkhIQg+B0HtTmmUVQRSxd9e+a1/Fsm5TF3XVtSuiIApIVTqEUBIIJUAS0kjvdZLp9/fHTSaZZAIBsfzY+TxPniS3nHvmzr3f8573vOc9giiKOHHixImTyxPZb10BJ06cOHHyy+EUeSdOnDi5jHGKvBMnTpxcxjhF3okTJ04uY5wi78SJEyeXMU6Rd+LEiZPLGKfIO3HixMlljFPknfzPIAhCc6cfqyAIrZ3+X9Tl2G2CIIiCICg6bRMEQXhIEIQTgiDoBEEoFgRhlSAIg9r2/0UQBFOX60T/2p/TiZPOKM5/iBMnlweiKLq1/y0IQgFwhyiKW7se1yb4SgdFvAXMBu4EUgA5cE3btuNtx3wjiuJNl7bmTpxcPE6Rd+KkE4IgeAIvAIuB1E7bY4H7gdGiKB7sdMryX7eGTpxcGE53jRMn9rwK/Aco77J9KlDcReAdMVcQhFpBEDIFQbj3F6mhEycXgFPknThpQxCE4cBYYJmD3b5A2XmKWAkkAP5ILp3nBUFYeEkr6cTJBeIUeSdOAEEQZMC7wB9FUTQ7OKQGCD5XGaIonhRFsVQURYsoivuQfPjXXfraOnHSe5wi78SJhAcwHPhGEIRy4FDb9mJBEMYD24CwNmu/t4iAcGmr6cTJheEUeSdOJBqAEGBo28+stu1JwAFRFHOQLP0VgiBMEgRBJQiCRhCEBYIg/BlAEISrBEHwbgu1TAYeAtb++h/FiZMOnNE1TpwAorSwgm2wVRAETdufFZ3cNw+1/bwD9AHqgL3AS237FwAfA2qgGHhdFMXPfvnaO3HSM4Jz0RAnTpw4uXxxumucOHHi5DLGKfJOnDhxchnjFHknTpw4uYxxirwTJ06cXMb8rqJr/Pz8xKioqN+6Gk6cOHHy/4r09PRqURT9He37XYl8VFQUaWlpv3U1nDhx4uT/FYIgFPa0z+muceLEiZPLGKfIO3HixMlljFPknThx4uQyxinyTpw4cXIZ4xR5J06cOLmMcYq8EydOnFzGOEXeiRMnTi5jnCLvxIkTJ5cxv6vJUE6cOLlEWCywaRMcOQKJiXDllSCX//7KvJR1AVi/HlatAqMRSkvhzBnw84M5c8DFBZKSftt6/wY4Rd7Jb8PFCMalEBmjEV5+GVJSYOxYePZZUKnOfQ34RcRNFEUE4cJXB2zUmxBF8NAoHJ5vNZkxTZ+O4tAhZC0tWBQKGiL70rh7H95ernhpVQ5KPQcWiySejz4qCafBAK6ukJwMDz0Ex47Z35fO93j0aElYT5ywP6a9zFWrpGtcf70kxHDue93aCkuWwLp1YDJJ5bi5SXURRdi9W9qGtPYigFBRAZmZ0j9KJfTvDwcP2n/vlzG/q0VDhg8fLv5u0hr8nqyW34Bmg5mCah2VTXpajBZajRYsVhE3jQI3tQJ3jRIPjQI3jQIPjRKtSt57wbJYYOZMOHAAdDpJMEaOhC1ber7HF3NOJ0RR5POd2VwzewTurU227U0u7ix+5Xvc3N3w1sh4+O8PEpFzHIW+FbSuWJNHYLWKyNMOIbS0IGq16IYmse/dryhtMqIzmDFZRKwmMxEHdxGQe4rSPvFkDR1LixV0BgvNBjM6g5lmg5niulaaDWZUchlGixWFTEApl6FSyKTfcgGFXEZ1s4EWo6V397MLU84cZNm6v+Fq0nd8fuCkfx/mLnmTUD83HpwSS3KUD37uagwmC7lVOnKrmlErZCT38SHMW2t338W9KWDQ2y1YKwoCKBQIJpMktCNHwg8/IIaGQl2d/bEAMhmGkDB2f7uNYQ8twTctFaxWAOnZGThQsroPHer4jpOT4f774dtvobAQcd8+W5mdyzfL5GC19N5q9fAAX1+p/E8+kbbdeqsk/u3bXFx6W9pvjiAI6aIoOlx/2CnyjrBYYMYM2LcPUa9HVGtoSBzBrmVfUtVipklvQi6ToZALqLASdXA3/rknaYgfROXYySjlMkJSd+CTnYlp8FC48kq83DV4uihxUfYshqIo0mQwU1avp6yhlfIGSWB9XFXIZQIyQUAuk14ImSAgE0AmCAhtv23bZAIqhQy1QoZaIW/73fa3UoZKLkMmk+pgNFsprNGx8Xg5Xx86S1mD3mHdzodCJuCqVrTdGwGtSmoMXNVyIn1diQlww0OjRCkXCEvZxpSXHkbV2mI736B24e27lrJMG++wfEfCpVNqeHDekxwdMo7BYZ4kRXozPMqHQA8NWpUchUxAEARaTRa+SC3E5eUXeWjf193E599jFvDG+JscXkMvVwKgsZi6XXd7TDIAMquFz1c+T2JpFhqTAb1SzZGQOBbPfwmr7Nc3DB5MWcEje5d3G3AzyhT8Z9R1mGVyMgP7sjM66bz1c3RP2um6SrlOqWFb3xHMPb2nx9XLRUAvUyII9ve0fZ+I/UChiICIaFeeo7LFc+zrqR624wVB6gV0wTJqNKUefpQ2m9k+aAL7E0ajUEnPsChCsKeGxAhvhkV4Ex+gRdy4EevhI2iSEqUCuvZwfkHOJfKXvbvGYhVZfbiYUC8XXFRy9EYLZqtIiJeGMG8tGqWchhYTG0+U0dhqwlWtIHLvT4zduROZ1YoACAY9HgdSWP/af9kaO8pWtsOXO7gfAjCgLNvuhZ/9G73wvxZmq0hDq/TSWi3S3w2tJmRWC2Ep21FU5HKoTVjuT93F9NZWu/MVBj1iRgaM6S7yMquFuad24dJFaDQmA/0r89iuS2ZHVhU7sqrOWcflxScdbk8qOQXAgIpcNCaD3T5VFyGyu26byE/KSyexNMsmhK4mPYmlWUzKS7cdc6mRWS1MyktnQEUumYF98frDPBRKBdtOVZIX3g+zTIHKarY7R2k1c+/+b5FbLb1uiBzdk3a6CqrGZGBIWfY56y0AGqsJR6ZlD6ZPr6JDLtTpZdfQtwl81zJk+1MJA8KA4albyIobyqapCwg7m0N+eD/WBA1iTUZpFx3QYxYEBARkiAhdepy2a12Em+5iuexFfvmBQp5fm3lB57zxw3KEtm5kOzLRyqzTe+1E3tHLPbxNMNqtlF/jhf+90pOF+2nSXPRKtb3FrFRzMiC6xzKGF5/s9hL2dE5PpIX1Z8zZY922p4cmAJAZ2LdbvYwOLPnO1w1wVzO+uaibEHZtCC4lDu9r2lqbYG8KG8r9vuEkVOV3u2ftwu9q0jOuJpe8CTKYMxtRFNGbrFhEyWo+W9tCdkUTNeShP7C6m+vHkUTplWqOBvcjoqHivJ+hJ4nruv3XkMLe1EUuWok7fZi4nKPILBZQq/lTSAi1L7/OmfJGEktPoTIZpYPF9j4J0NyMYe8+alZ8hzB3Ds98f4K8qmZWXhtDwKhhUFsLPj5w8iT4O8wU/LO57EMoZw4IYuaAwAs8q3curJ4sv67WX/sLfzlww/BwbhoV0atjOzeCckRbgycCR0Li0Ck1WBDQKTUcCYljZ3RSj2VoLCbbSyciuVF6Oqcnlo2eT73azeYWEIF6tRvLRs8HYGd0Urd6pYcmkBaagE6hxoLk8ijwCmZ31FAAKpsM7HELR69U212rawMUG+DGlPgAloyJ4prE0F7X2RFT8g8zsjKn232dlJcOgFUm5+qb/oZOobH7rF2xNuvY9OVmvtxfyImSRk4U1fLq3a9TNvc6FItvJuSfrxJXdJqmsEgMGhesgFnoLhkioFcoyY8ZyO6nX6fF1d3uur8fh/DPQwbILZa23r0B8vNxXbyIwU/eh7Jd4B2gNLSy/7V3GfvKT2w/XUltfjH+MRGI1dXSmER1NQQEQNW5e6IXy2VvyQd6aHj/5uE0tJo4W9OCTAYGs5XiulY2HCtlZ1YVJqOJeaVHmSdW4jlmBDU3LsD63B7kYoc1bxEENsSPs/0vs1pQWC1YZHLknbr
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd1hUV9rAf/dOn6EX6UhHsAAWFHtJs6QaY0zvvW/KJhs3dTe9mt67MZrERBNLjDUWQAQLiAqI0pvUgen3+2OGkYFRMTHZ/dz5PQ8PcHs773nP244gSRIePHjw4OH0RPxPX4AHDx48ePjz8Ah5Dx48eDiN8Qh5Dx48eDiN8Qh5Dx48eDiN8Qh5Dx48eDiN8Qh5Dx48eDiN8Qh5Dx48eDiN8Qh5D/8zCILQ0ePHJghCV4//L++17a+CIEiCIMh7LBMEQbhLEIQ9giDoBUGoFARhsSAIQx3rH3CsaxcE4aAgCA/81ffowUNv5CfexIOH0wNJkry6/xYEoRy4QZKkNb23cwh8hZtDvAbMBG4ENgMy4ELHst2AAFwF7ALigdWCIFRIkvT1qb0TDx76j+DJePXwv8ixhLwgCL5ALnZhvRVQSJJkEQQhESgGsiRJyunnOV7H3sbuPKUX78HDSeAx13jw4Mq/gbeB2l7LpwGVJyHgBWACUHhqL8+Dh5PDI+Q9eHAgCMJIYBywwM3qQKDmJA73OPb29fEfvzIPHn4/Hpu8Bw+AIAgi8BZwt8M803uTJiCsn8e6A7u5Z4IkScZTeqEePJwkHk3egwc7PsBIYJEgCLXY7fIAlYIgTAB+BSId2v4xEQThOuDvwDRJkir/zAv24KE/eDR5Dx7stALhPf6PAnKAEUCDJEkmQRDeAhYKgnAjsAW7knQBECNJ0rOOqJx/A1MkSSr7ay/fgwf3eDR5Dx4AyU5t9w/Q4FhVJ0mSyfH3XcAbwJtAC1CKPYRymWP909ht97k94u/f+evuwoOHvnhCKD148ODhNMajyXvw4MHDaYxHyHvw4MHDaYxHyHvw4MHDaYxHyHvw4MHDacx/VQhlUFCQFBMT85++DA8ePHj4f0VeXl6jJEnB7tb9Vwn5mJgYtm/f/p++DA8ePHj4f4UgCIeOtc5jrvHgwYOH0xiPkPfgwYOH0xiPkPfgwYOH0xiPkPfgwYOH0xiPkPfgwYOH0xiPkPfgwYOH0xiPkPfgwYOH0xiPkPfgwYOH05j/qmQoDx5OOVYrrFgB+fmQkQHTp4NM9p++qj+H3vd61lmwejXk5dnXyWQwYsTp+wz+l971SeAR8h6OcqJG8v+tEVmtcPbZkJ0Nej3odDB6NKxa9edet7vnZLXC00/D5s0wbhw8+igolSd1WEmSsEkgE/vMP+u8V8lxrza1BotcgWgxI+/qBEAA+zMYM6bvM7BaYflyWLzY/v+cOTBr1u9/Tn/0eCf7rZ3oXf9/+3ZPIf9Vk4aMHDlSOm3LGvy3f2Q9Gomk1yNptXSkjaDgg0V02UCBjZE3zUO3Mw+xsxO0OhidibB69X/XffRk+XKYNw86OpyLjHIlj182n8JRk/HTKgnQKvDXKfHXKhkYqGV4tD+R/hoEmw1WrMCYm4dxyDDEc85GtWYV5oWL6DJbKZ86k+px07CJIharhNUmYbRY0XeamH7/1YQW70Jh7MKkVFMak0r04X14dR29jja1FzMf/gZRrUYhE1DIRFRyEUEQqG8zUN1qOKlbnVqSw4Ifn0dnPrqfhEOw90LS6WDhQqSZs6hu7aKkppUh184hcPtW+32D/Z1OmoR15SoONndhsUnIRQFvtYJAnRK5zL2lV5Ik+zHOOgs2bACrFQlAFGkfM449Hy/GR6fGT6tggLcapdzNcaxW+/5btoDBAGo15tFjmDHzUarazZyVGkJquA8Gsw0ftZxdVa34rFnFA5887nL/BpWGQ298gN/s8xlw8XkIOTl/bWf/FyIIQp4kSW7nH/7f1OS7BW5uLuzbB3V1EBICyckwapRTAJssNpr0Rhraj/4YzFbkMhGFTEAuisgdDVQuOn47lqsUIn4aBb4aBb5KEfmM6SfUKCVJorXLTH27kfo2IxISoiAgCCAKguMHBEFAJtr/7r1eJgqo5Pbzq2QyVAoRpUxEdKP9SZJEY4eJNXvrKPloIfdt3IzObEAABL0eWW4On8x/m7UJmUwtyWF0bg6y7kak70C/aQsPXvEvCtLG46tRYLbasNok/LQK/LVK/ByCc4BWTlT2BgL2F9KQkEr5qIl02cBgtmG0WDFabDTrTfxW0kh9u7Hfr3F4tB9pUX4MCfclOdSbMF81WqUctcIuLMnPR9LrXQSd0mLilh/eZErIMGyi+wYu2qx89s0/yajeh9psxCJXYpbJ0Zi60Eg2NIDvD9/SFTWUq+Y+5XKcqSU5XFlYgMrxnNTGLhIP7ERhs7hch4+hgzmrPuOVCVf0+36Px+C6UtTm/j07m76TV174hjc2i85rfisv56iAB7BaMWzYxG3znmJtQuZJXcvUkhze2rAJtdUKODoamw1l9jY+eOSt4x4vJlDLdU27uWLdekTJcT0GA7KNG4kOXM+BxDEsLajmxx0VTC7LY3BdKS0h8UTWHuhz/wqjgWWfraBo5T4W/LYFncnx7XZ02NviihX20UUvrDaJbWVN/Ly7hro2AyarhFImopQLSBKE+qoZHu3P8IH+hPuqAZAk3Lax/wZOayFvsdrIr2ghUC1Dt/YXxIICDIOHEPbJe8hzsqGzk97jGJtaQ9fITJY88S5b3llIck0JhSHxrI8bcUyhcCKmluSwwCFAAejoQL9xM3df8jg2QWRwXekfPscf5c7C3X0aidpsJLW+jLUJmW6FiNpsJLbqAD8NHEFVS5fb43YLzHSHwIxWqBDDk7nqkiexiTLkFhN3bv2GqZVFDIxMZUHWJVjk/TNj7Djcwo7DLcdcP7XEwluiHLXV7FwmACEdTUwuyzumsJlclkdG9T7n+9JZjEgWo4uQlksSI6v29jmOu+eksFncnmdEZZHL/6LN6hRc/f0ewn3VjI4LpLAkHoNC5aLJHguDQkXRgDiXa1b2eEbdKK1m5/s/Gf7I8cqbOvFe9h2CZHNZLkg27tv0BZIgsjEmnU+WPOHshA0KFYf8wjDIlegsR599930OritFbXJ9J5JejyVvB/KZM9l0oJGCihbmZUYT7K3iwSW7+HZHJV4qOdEBWhRyEbPFhslqv6Z1++r5eHO52+uP8NPw6qXpjIoJcF3xHxzJn9ZC/oXV+3h/3QEXrcwskyPYrOD4iHr3vTJDF8K2bWReOp05LTXOjyi/h2A6WY4lIOf/+gFBnS2n5By9OVmBURjSV0j0FAYnWn+sc8ttVleBaTaQUb2PyWV5bIxJJ/eNq/Az2s0YYw/v4qq85Yy647N+C/rjsT5uBHVeAUS31rlq81YzQxoOHlPYDKktQdMPYelOaLl7TmZR7lbQR7XUIdqs2ERZn9FDf7+H6lYD3+dXIcaNID88+ej+ciUmmQKlzYK2571oteiysvho0T+RRJEusxXj90bIXQJG12/UJFO4fb8noiN1KOZtSlQWk8tyi0JJ6vSJrLtxMnVtBi59b9sxjuBeI05pKGfBj89zyC+Mgc3VToGuMxsY2FJDuX84Mb3a7Pq4EQB93kmnXMWdu62sffhn57KXf9mPt0pOu9H+rjqMFlo6TUxKHsC0QQOIC9YRHaBFslqp+uo7mjdns1wWysc+g5zvqKqliznvbAXg7slxXN9ejGJnPk2ff01gXSVqkwHhLzYXndZCfsaQMOq++tZFyMis5j7ae2/UFiMJTRUoHQ2zp2ByJxhOJFDdNnyZnJCOJqeWeaJznAy/R2Cs7y0kejWS9XEjaB02HHF3PkqToc96l3N//SiZVUXIbVZsgoCsl9+ne4SQVrMfP2OHS5P2M3Zw59ZvTokZwybKeHLajbz5w3Mu2nyXQs2e4Fi3+4g2K2fv33oMMeNKbyEo2qyIko1GrR+yjiYUVjNmlYbGlDRCinciN3Q5jysAQfpm5/vuM3o4xvfw/OxhyGUC9e1G3lxb4hRINlH
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hUVfr4P/dOn/SEJCSQEJIASaihhl4VUYqKAjas66KuZW271lWsX9eCinWtqKugKE1QSqTXhE4SSkggPSSQNpPp9/fHnUymBYLi7v5wPs8zTya3nnvm3Pe8533f8x5BkiQCBAgQIMDFifjfLkCAAAECBPj9CAj5AAECBLiICQj5AAECBLiICQj5AAECBLiICQj5AAECBLiICQj5AAECBLiICQj5AAECBLiICQj5AH8YBEFocvs4BEFodvv/BkEQbhYEIVcQhAZBEEoFQXhFEASl1zWuFwQhx3lOhSAIqwRBGOHc10sQhJ8FQagRBCEwASXA/wQBIR/gD4MkScEtH+AkMMVt21eAHngA6AAMAcYDD7ecLwjCg8A84EUgFkgE3gWmOQ+xAouA2/9DjxQgwDkRAjNeA/wREQShGLhDkqS1ZznmQWCsJElTBEEIA8qAWyVJ+vYc104FjkqSJFzIMgcI8GsIaPIBArTNKOCQ8/tQQAv88N8rToAA54/y3IcECPDHQxCE24CBwB3OTVFAjSRJtv9eqQIEOH8CmnyAAF4IgnAl8BIwSZKkGufmWqCDtyM2QID/dQJCPkAANwRBuAz4F7JT9oDbrm2AGbjyv1KwAAF+JQEhHyCAE0EQxgFfAdMlSdrpvk+SpHrgaeAdQRCuFARBLwiCShCESYIgvOI8XxAEQQuonf9rBUHQ/IcfI0AADwJDzwABWnkKCANWCoIrMGaTJEmTACRJek0QhErgSeTOoBHIBV5wHtsFKHK7XjNwAkj63UseIEAbBEIoAwQIEOAiJmCuCRAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmP+p6JoOHTpISUlJ/+1iBAgQIMD/V+Tm5tZIkhTtb9//lJBPSkoiJyfnv12MAAECBPj/CkEQTrS1L2CuCRAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmICQDxAgQICLmP+pyVABAvhgt8OqVbBnD2RmwqRJoFD8t0v1v4+z3uy5u1EM6O9Zb951eumlsHp1oI4vUgJC/o/CbxCWDoe85oAoCuc48jdit8OKFbBwIZSXQ8eOkJ8Px4+DwQBBQTB4MNx3H+zff17PkVN8mq2FtRgtdsw2OyEaJRFBaiL0audfFXFhOqJDzn8hJ0mSaLbaaTLbMJjtNJlszu82DBYbjSbnd7ONJrOdJrMVSQK1UkSlEFErRdQK+btKKbi+N5ltVNQ3U1lvprrRRGW9iTqjlYggFR1DtUQFa1CKAidqjZw8baTZagdAdNhZsOhpMssPo7WaMag07Ivvwb6PF3LXmFSYOBF27JDrVK8HtRrJYkEwGuU6HjIEfv5ZfrhztZmW3+zbb+X/r70WJk/+XTqJ/aV1ZBdUM7VvPLGhWvRqBYIg4HBIlJ5p5mh1I8dPGaioN1HZ0Ez/xAhuH9EVtwVg/pD8Ty0aMnDgQOl3S2vg3RivvhrOnIE77vA8TqEApRJiYmDYMLjuOpg8GYcgUtdsxWixoVKIKEUBpUJEpRBQivLf/9nGZLd7vtjuL7LzZaxuNPHzoSqy95cw+ut3GLdrDZGmRgRBoCQsltdHXM/G9OFEhuoI0aoI1ioJ0SoJ1aoI06kI18t/Q3UqwnXObRoFHTatQ7N/LwrJgahSYu3VB7PNjrB3H3abDRsiDRm9qcoaTa87ZhG+cwuCV5t0r1WHIOJQqVBYLRAUhOD1HP46sy93lfLkkoP/ocpuH0qbhXu3LWJgaR45nTN4e+gMbEr1Bbn2uGM7eXvZKwRZTa5tBpWWe6c+CuCzT8Kzjg0qLa/OfporNy2mW3EeGosJk0rD/k5p3HfLS3QIDyIqSE2kTsED/3cPSYdyEBwOQP59tif0YvUbC8irNlB4ysBpg6XNskboVWTEh5LWMZTuscF0jtCjVYkoxJZ3TH6/wnQq7vwihz0n6867PuLDtGTEh9KjYwj9EiKYkB7Trnf1+Kkm9pfW02iyIopy59tSnl6dwkiK0nte57846hQEIVeSpIF+9/0hhLzdLg9JN2yQv7cDye2vRVRSFBpLXUgEHw+cSnbqEBzi+f14ggAdgjVEB2voEKKhodnK3pLzb7AKUSBMp0IUQBAEFILg+i6KoBAENEoFGpWsHaoFiSlL/sX0NV+hdthc12l56bNTB7u2KW0Wdr11I+FWo997l4ZGM/72d7Godecsp0ubLM1Hb299yR0ISEi4154NkfKQDsQaTqN1K6M//AmkludouWf/iiNoLSaaVRr2xnXnk4FTyagu4lBsCuuTB5z3b9fW8405nkvPqsLzvq7SZmHX/NmEm5tan0Op4cEr/kqP2pLfXM57t3zNA5v/jYLWd9uOwBsjb0CQJJ993nVqR2BpxmgmHt3ut6NoaTPjju3k3SUvobVbPe5vUqi4+8rHPNrW/xIhGiWb/jaWcH1rp3q0qpGYEC1hehWnDRb+tCCH3BNnznqdyCA1mQnh7C2p40xjM98s/gf9Ko6gMjUj+FGkfk/OJuQvanPN0r1l3P/NXrkxbtyMtp0CHlobvQBoHDbS6sqgrowhJQfZmtCb2bOeP6+XUJLgVKOZU41mqDi/53DH7pDOqhm50yL0BpccQuUlPLVWMxnVxz1exHu3LSLcaqQtHadzwyl2v3UD/R74BoeoOKuQG3M8l8yyAoLsnmV1Fy4tqHCQ2Fjdrmfyxv05xhzPJbP8MHqnYAqymsgqOciAsnxUdhsmlYY98T2YPWPubxL03uaQ873uvdsWEW5u8qjnIJuZ95e+jITwm8t5KDYFk0rjIaDNag2xI7NIiNTD7iVgNLR5vkmlQUCuW3e820zPqkLUXgIeQG23+rQtb/omhLPvLErOb+lEz0Wj2Ua/uWvQqxWoFCL1za3PcHX/TkgSLgF/95gUMhMjSIrSo1UpsDkkjBYb+0rq2X3yDLtPnqHWYGHc8Vx6lhSgbqnzpiasW7dhX7YC7VXTqDNamLf2KLeP6Cr/Bv9BLmoh32yRhXrPqkLUtvYJRn8IXt8Hl+Ux5njur9JUfs9hujctQk/tRzs2qTTkxSR7bBtYmnfW6wlAkN3CXzd9xdii3aTWlqB02GhWaX2EUs+qQrQ281mv533ttsaUkvNjUahQOmwo3Uaf7s/Rs6rQRzCJkgOtXTYlBFlNZJYf9vvbxYVpkSSoa7ZgsjrOWtaWeg1y60y8r9sjNoS+CWH06BhKSnQQKoWIxe5gx/HTDPzat54FWtqZ1GY5FaLAvJn9kACj2UZZXTPzfzmG92B8ffIA9sT38OyE4nrwgtSFnppQElN7kXj0AApTM+j1CE6bPEYj6IPQDB7E2DvvRrhjl2ze81PXIHcmFoXKR5O3KFSciE/hjrpD9K0twpLRi9SYYDJOFSP160duxlCeWN52W/utnWh7MVrsgKfi9/3uMo//311f6PfcmBANswYl8MScYby3oRDNlq/RunWqAGKzkTff+I61JyI4XNWIJMGK/eW8e8MABneNvGDPcS4uaiE/a3AiEzJi2f5WGfYd3yFa2y90zobSYadnVSGi5OCKgs2AxIq0kfySMuisjdB7mD7s5H5m565g0F8WXHBBLzrsTMnfiM6r4UmAVVSyJ74H65MHeOzL6ZzBsJP7z3ntO3KWoHbYXZ1fkNXE4JJD3L/la94cfh0OUcGh2BTsohLFOcwv7rgLeg87PALLMkbxY9pIbslZRmbFEY+Xv+U5/Gmw3uisJibnb/TRDCvq2z7HG3+dibeWe7iqkcNVjf7L0I56dr+eu1b7w5Ed59RqHaKC2TPmMuZ4LhnVx8mLSZbPcQjkljYyYeLjjOnWum9jUj9GFe/1OJbdsKBDKpkWT0G
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hUVdrAf3d6Mum9d0ijQ0KHUKQjoKKIva6KrsoiKnbs7lpR7F0BRRGU3nuHkACpJIH03pPpc78/JhkyZEJRXPfD+T3PPMncdm6b97znbUcQRREHDhw4cHBlIvmrT8CBAwcOHPx5OIS8AwcOHFzBOIS8AwcOHFzBOIS8AwcOHFzBOIS8AwcOHFzBOIS8AwcOHFzBOIS8AwcOHFzBOIS8g78NgiA0d/iYBUHQdPh+kyAIPQRB2CAIQrUgCHYTSARBmC0IwuG2fcoEQVgnCMKwtnUX3N+Bg/82DiHv4G+DKIou7R+gEJjaYdn3gAH4EbjL3v6CIMwF3gFeAfyBMGAxMK1tk/Pu78DBX4HgyHh18HdEEITTwN2iKG62sy4GyBVFUeiwzB0oAe4QRXH5BY7daX8HDv4qHJq8AwcXx2BABfzyV5+IAweXgkPIO3BwcXgD1aIoGv/qE3Hg4FJwCHkHDi6OGsBHEATZX30iDhxcCg4h78DBxbEP0AHT/+oTceDgUnAIeQcO2hAsqABF23eVIAhKAFEUG4BngQ8EQZguCIKzIAhyQRAmCoLwxoX2d+Dgr8Ix9HTg4CzhQEGH7xrgDBABIIrim4IglANPA98DTcAR4OWL2d+Bg78CRwilAwcOHFzBOMw1Dhw4cHAF4xDyDhw4cHAF4xDyDhw4cHAF4xDyDhw4cHAF8z8VXePj4yNGRET81afhwIEDB/+vOHLkSLUoir721v1PCfmIiAgOHz78V5+GAwcOHPy/QhCEM12tc5hrHDhw4OAKxiHkHThw4OAKxiHkHThw4OAKxiHkHThw4OAKxiHkHThw4OAKxiHkHThw4OAKxiHkHThw4OAKxiHkHThw4OAK5n8qGcqBg8uOyQTr1kFqKvTtCxMnglT6V5/Vf4e/67X/Xa+7CxxC3sGVi8kE48fDgQPQ0gJqNQwcCBs2XNE/eq3BRFV9K67TJuNy7AhSTSuoVAhBQfDWWzBlyuW9/naheuSI5X+pFPr37yxcO25nMEB+PkgkMHPmxZ+TXg8vvQR79sDQofD006BQ2LbxN3zm5+N/atKQAQMGiI6yBn89ZrNIZZOOulY9GoMJjd6EIICbSo6LUoarSoarSo5C9j9u7Vu9Gm68EZqbrYs0SideveVZMvuPwNNZYfmoFXg6ywnycKJPqAchnk4IggBY7oXWaEIpkyKVCJjMIpVNWsobtOiNZkxmEaNZRGc006Iz0qwzWv+2/9+iM9GkM9KkNaDRm2jWGRFFUMokyKUSFDIJcqlg/V9nMFNSr6GkXvO7L330qYMs+vUN1AatdZkIoFLBkCEIGzeCVIpGbyK/upmyei1xga6EeDp3OlZ9q55GjRGpVEAlk+DprEAisdwfq1Ddv98iVNtxcYGkJBgyBHHvXrRJAxH37EWZegRJ69ntBABBAD8/SEmBmBiQy6FPH8sG6elWbfzzbblcO2UA7rqz+5sUSjY99TZrowaQU6UhfN9W3l75us11m9VqhKVLEaZO/d33838dQRCOiKI4wN46hyZ/HkRRpFFrpKpJR3WzjqomHRq9CZlUQCaVIJdY/sqkAnJJ21+pgEwiQSmX4O4kx91JjpNcahUaF8JsFmnQGKhp0dGkNQIgEQQkgoAgtP0vaV8GQts6Sds6qURAIZOglElQyqTIpcJ529YaTKQW1vNLajE/Hi6+pPujkJ29RncnOQqpBCeFFGeFFLVChpuTjChfF5wVUuRSCTKJgN5kRmswozWY0BpM6Ixmy8dgIrO8iUMFtWgMpgu2rZRJGB3nR68QD+ICXQn3csbfTYWTXHpWAKWmIra00PHqFTotbtknOeTd85Ku9f8biRV5qAw6m2UCgFZLy669PHTDQrbGJP/hdkafOsiinXtshCoAzc2Yt21D2LYNAFXbX7tvoigiVlTADz+cXSSRIJjNAJglUrQBQfi6huOus32eUr2OMS/8E9fQRG69fiHjyk51um6xpZUP3v6JHcXeJEV68Y8RUXg4KzgfVU06fksrJaOskWatEZlUQCGVIJUIKOUSegV70C/cgyhPJyQb1v9Pm4auaCFvNovkVzfjopQjkYBWb0ZEJNDdyaqFms0iOZVNGE0iaqUMmURAJhU4dLqOp1Ycp0ln/Iuv4n8XvdFMVZOl8/tvozOaWXeinHUnyrvcZvQpI4tkShsBpJUryfCLOu+xJWYTKflH6FGei1QUMQkSTgTEsD2qP2bJ/9YPGGBCYgDrT9reh5P+0Wjlys7CF1AZdCRU5l8WIW+vM2lHwFaon89mcK7wbxfwAFKzCefSIqZSZHc/hdlI39JsRuUdQmY2YZJIkZrP/m61ciWpXhEcPlPH4TN1fLg9j8U39WN0nB8quZSqJh2rjpXQP9yTvmGeHDlTy6xP9mMwiQS4qXBVyTCJonXk1qwz8t3+QiRmE0t+eo6exVmoDDq0ciVpwbGc/OJH7k7pdtGK3Z/NFS3k39mcw3tbT/3Vp+HgL2J7VH9Sg2LpW5qNyqBDp1ByKjKR1MRBoDPb3UdiNvHNj8/StyQLZ+NZ4dUqV5EaFMut1y/8nxP0G4+XMDr/CIkVeZz0j2Z7VH/rtQ8ozkBpMtgIUYmLmnnzb+CRiRNp0hpp0lpMS01aAydKGymua2VnThV5VS1dttnO+TqTy0n7+XfVUTgZtLy/8jUUZiOSDtu1P7ftUf1ttn/g+6MX1a5UIpAc6cWgKG/iAlwJ83ZGLpFYTFzf/UTfshyUbdeuNmjpXZLN5+9/y74PJEwTK+hz9SjCbrqWTdnVvLEhm1dm9CQ50uuSr/+PcEXb5A+fruXWLw7Sqr/w8N/BhVkwKY6DBbVszqz8q0/lomnXyhMq88nwi7qgNm7Plt1Oi1zFQ1fPt6sBt7fTUdC6OCsZHO3NiO6+hHg6o5JJeGtTDgcKan/39Tw1KR5fVyXNOiMrU0s4WlBt6ZTaOjKtXGntjJzkUoblHGDB5k/xb6pBbjKilSvJiUhg+cuf0TvCmz6hnnTzc2FHbhVpRfX0CnHvYE6z/C2t15Be3MDB053PW2I28cuqhcSePolSf/aeGQUJMtF8Xk3+cuq57cfueEy9IOHTYTfw5qBZf0rH/NCepTy6ZwmSDjLUBBS7B+DTWt/peZglUhQyCU9OjOO2wRFnzYqXgfPZ5K9oIQ8Wu3pacQMldRqkEgGNwUh+VQsrjpZQUq9BIkB3f1fGxvvj7aIgp6KJ7PImjhbWX9bz+DsT5asmKdwLmVTgaGE9mWWNf/UpdclDe5byyO4lSO3ojCYE3h5+E+8PmWWz3Kr92xG0f7bWb69TOrczutSOritUbbboARGeBLqrOHS6jl/TSjt0cKeQm0xMP7mN0Eb7isDFiDV7Arurbc53XBF4a9hsFg2dfRGtXhoSs4mH9yzl/v0/oehoGpLKkUgEFAa9zbJ1sUPZ3X8sB+OSKWq0rNv35GgC3Z0uy/n8rR2vgiDQJ9SDPqEeNsvnjo6m6PsV7Fq+iS25QSwutbz41/UP4dr+IVe8kJdLBbr5uZLRJnDvGhZJpI+aBo2BRo2BVr2JFp2RFaklf7gto0kkMdiNqb2C8FRbHF6VTVpOVTbz+ros0oobrNvOnxDLvrwaduVW/+F27WFP4+4o8M5nfujKnp+Sf4S+pdnWfdQGLX1Ls0nJP8L2qP7nbe+PYs8mrjJomZK506bNrTHJf9gGrzWYOXi6tpNGb5ZIrccffeogD+xf3knoXooqaRCkmCUSlCaDdZk9Ia5xdqUsKp7IEwe77BBkHWz7F3r2HZnSKxB/NxWf7y7otK6jSU9uNtqYhqqdPQhpqLDZXmkyMC1jO+Ny99t0/gcLapnWJ/h8t+KycMULebuYTAgTJhB24AA3trQwXXZW8/rpSDE/HblwlEl
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"data": {
|
|||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEICAYAAAC6fYRZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydd3hUVf6H3zt90nvvlRQgISGE3kGKyipgBbuurl3Xsuuurl3X3sW1V0QsiKB0pEMCoSSQSkjvPdNn7u+PmUwyKTTL7g/nfZ55AnfuvXPr55zzbUcQRREnTpw4cXJuIvlvH4ATJ06cOPntcIq8EydOnJzDOEXeiRMnTs5hnCLvxIkTJ+cwTpF34sSJk3MYp8g7ceLEyTmMU+SdOHHi5BzGKfJO/jAIgtDV52MRBEHb5/9X9Ft3oyAIoiAIsj7LBEEQbhcE4YggCN2CIFQJgrBCEIThtu+nCoKwWRCEdkEQyn/n03PiZFCcIu/kD4Moim49H6ACOL/Psk971rMJvnyQXbwM3AHcDvgACcC3wDzb993Ae8Bff8PTcOLkjJCdehUnTv44CILgCTwMLAV29VkeD/wFGCuK4t4+m9gbB9vyvYIgzPidDteJk1Pi7Mk7ceLIk8CbQF2/5dOBqn4C78TJ/zxOkXfixIYgCJnAeODVQb72BWp/3yNy4uSX4xR5J04AQRAkwBvAHaIomgZZpRkI/n2PyomTX45T5J04seIBZALLBUGoA/bZllcJgjAR2AiE2Xr7Tpz8v8HpeHXixEo7ENLn/+HAXiADaBRF0SAIwhvA54Ig3ADsxNpJWgBEiaL4tG00oMAamSMIgqACLKIoGn7PE3HipC9OkXfiBBCtEyvYna02gQao72O+ud32eR2IBlqB7cCjtu8nAZv77FYLbAWm/GYH7sTJKRCck4Y4ceLEybmL0ybvxIkTJ+cwTpF34sSJk3MYp8g7ceLEyTmMU+SdOHHi5Bzmfyq6xs/PT4yKivpvH4YTJ06c/L8iNze3SRRF/8G++58S+aioKHJycv7bh+HEiRMn/68QBOHEUN85zTVOnDhxcg7jFHknTpw4OYdxirwTJ06cnMM4Rd6JEydOzmGcIu/EiRMn5zBOkXfixImTcxinyDtx4sTJOYxT5J04ceLkHOZ/KhnKiRMnZ4nZDKtXw4oV1v8vWgTz54NU+t89Lif/dZwi78TJ/yPMFhGpRACDAR5/HHbsgLFjrX+3bUM0m60rfv45TRlj+ebp95iWFEDo7q2oDh9EGDUK5sxxiv8fCKfIOxmU1m4D1W1auvQmtEYzep2RsN1b8CnKh/R0ZPPm4u6qQiWXIAjCb3swZjOsXQsHDkB6+hmLlCiKrMipYnlOJZ06I1qjGZVMirerAh8XBd6ucrxdFPi5KRkR5klqqCcquXX/ZotIm8aATCpBKZNgMFuobtVS265FZ7RgsoiYLRZMZhGjWURjMNGlN9GlM9FtMNGlN9OlM9KtN9OqMaAxmNEZzZhFEblUgkIqQSGTIJcKKGQSZBIJ3RodsbnbiK8uIT8wli0xGVgkvecrMxnY99pSvPRd1gWbNgEg2D4AWCy479/Lntc/Jjn3e4JrCrEY9eDigiUrC+m6dUjkMlq7DZQ0dlHTpiU52IO4ALcB97O6TUuH1ohMIqCSS/F3V9qvz2DXuv/2RrOFunYdrRoDXToTbioZXmoF/u5K1IrTu48Wi8i1H+5jS2Ejc1KDGB7miSiCu0rG0doOjlW1ErRzC3HVxfZrFh/sxb8uTCHUS02Yt/q3f07/R/lji3zPEHf5cqipgZAQuOQSh2Gu3mSmqctAY6fe/jFZLMgkEmRSAblUQCaR2P9al0mQSQRcFDI81XI8XeS4K2VIJKd+yERRpFVjpKFTR3OXdWpQQQCJINg+INj+SiXWZf2/7xEkpUyCUi5FIbUe32APucls4USLhu8OVPPF7nJSD+0kpb7UQVwkFjMffflPomsKURt1mCQySnzDGbf0eVAo7efoqpCh1emZUJxDSn0xrhIBlUqBOGoULZOmI1NYHzed0YLOaLb/1Zusf0sbu9hZ2ozZ0jtbWc9vp9cUojLq0cmVHAhJZOniRwn1c2NkmBcjw7yIC3QjMdAdN5UMF7kUmdTqbrJYRD7Zc4J/fpf/ix+X34OTnW+P0N+260u89F30vZuDze+mMBuZd2w76TWFuBp11oWabnQ7dnHbZY+xKS7rNz+fMyXIQ8WEeD+CPFTW90gqoJRJcFHI2FLYCMDaI3WsPWKfqfGk1+zSZbsB8FTLifV3Ze7wYJaOjUIhG9odKYoiuSdaWX2oluo2LQaTBbntnRIRCXBXMSrSm4xIb0I8VbZtOK33+7/BOS/yepOZ/JoO/N2UiCLoTGYkgkCYhwLVvDmwdWvvEBdg+XK6sydQ99V3fH2wjje2lP73Dv535GQvypSyXAehUFhMJDUe57uP7uH8q1+iuVukudvQu4/qY7iY9PZ9a1aqBgjV6dL/t12NOtJrCplSlssmSRaVLVpWH6r99S7E74zEYmZKWS4p9aUUBEQzsqaQMRWHkYsWoN/52kQ5s6rgtPZtkMoRAJVR77BcZdST3FD2PynydR06vsqtOqNtTvqM2M6xXWtkf0UbeeXN7Hz5Q56ONuE1PgvF+fPJrWpne3EzizLDCPFS88QPR/nP9uOo5VIifV1sIzgRg8mqE5uPNfLBzvJBjyXMW80Li9PIivY5+4vwK3POi/xfVxxi1cGaAcunlezlja3bUJnNDj0iLBZke3bx5G0v/k++BH2RWMxMK9nDdTmr8NR1sSFuDK+MuwSTTHHG+zrZi5JSX4q6pydoQwDimysdXiT7PkyOojLYS3e6pNSXnpVI9RXPvqOS8XG+tGmMVDRr6NSbhtz+90BiMfPx8n+QUX0UhdmIyODhbiqjzuF8c8KSGVdxaMB6FnrNNRZBQm5oEj8Mm8Cs4t29PXlA4ubKvfddwr3z56EzmumwmZNEUaSiRcOJZg3bihvZcLThrM4rxs+VAA8llS1au4lKJhFIj/AiK9qH81KCUcgk3PRxDuXNmrP6jb6c7jMyoCPzspK9fTofL24owl0psz8XWqPVtDY6yocpif7E+rsR5q0G4FhdJ7knWll1sIbcE63236hq1bL47V0A3DkjnqVjo5BJBR5ZlU9du453rxp92iaqX4tzXuQXZ4YPKvIp9aUozMZBt1GYjfYHZCix+G8jsZj56IuHGFd52P5iJzWWs3T/ajJv/fiMhf5kL0p+YCwmiQyFxVEUpRaTw4s02D767+tMRT4/MBadXOkgUjq5koKAmCG3OdmoZEdJ8xn9/i9ldJQ3k2N8GF24h8CSo2x1C+MxfSiTyvM4v2Ar2RWHkQ5qbOnFKJXbz/e2aXHELnwKw7S1KDrb7eu0KVy5f85tnFdsNU+sGTaBHYljEEXIC0kkrc+1KAxK4LOuEEbuPkF6uBeJQe7sO15Pfk07I8O98HKRMys5iInx/uiMZpq69Bysamfv8ZZBjy8r2gd/dyVms0heZRtlTd2UNXWjlEkYEebJqAhvEgLd8XKRc7ypm6ve30tLt+FXusJDPyONscPwcpHTprG+56fT4+/f8Jc3ayhvPsHHuwev5Ds7JdD+7/5a8YrFzEsbih3WP/+17bxyaTrJIR6//MRPk3Ne5CfE+1H25FwOVbfT3KVHIgh06U3o1RUYdq9AZRr4sJnlCvQpw0/LPvrfYkpZLlnVBQN6fp76bm7b9SUvTrxy0O2GarROJqZbYjIo8Q0nqfG4w6hHJ3MU28H20YNJpWLYeRO5KiWSD3cNWfp6AFtiMjgQkjjgHmyJyRhym9N5mX8vcsuauOOpmxlp661fIpFxgVyJwmJGbdRxKiuuCNS7+SCIFm7b8Tn5JbG8HpOB5Mb3uW3Xl2RUHyU3NIlXxy7GJFOwbtiE3o2tFh+W2ExuyQ1l9vtpyavlq7zTN3PJJAJuShlp4V5kRnkT4qlma3EjPxyqHVL89SYL+8pb2VfeOuj3J8NdKSPYS0VRfdcp1x3qGVnhPxyLprcjd7ajwpPxU3490NuxyKgrQqH
|
|||
|
|
"text/plain": [
|
|||
|
|
"<Figure size 432x288 with 1 Axes>"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
"metadata": {
|
|||
|
|
"needs_background": "light"
|
|||
|
|
},
|
|||
|
|
"output_type": "display_data"
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"source": [
|
|||
|
|
"for key in action.keys():\n",
|
|||
|
|
" if key[0] != 'T':\n",
|
|||
|
|
" continue\n",
|
|||
|
|
" spikes = action[key]['cellTS']\n",
|
|||
|
|
" spikes = spikes[(spikes > t[0]) & (spikes <= t[-1])]\n",
|
|||
|
|
" plt.figure()\n",
|
|||
|
|
" plt.plot(x,y)\n",
|
|||
|
|
" plt.scatter(sx(spikes), sy(spikes), marker='.', s=100, c='r', zorder=100)\n",
|
|||
|
|
" plt.title(key)"
|
|||
|
|
]
|
|||
|
|
},
|
|||
|
|
{
|
|||
|
|
"cell_type": "code",
|
|||
|
|
"execution_count": null,
|
|||
|
|
"metadata": {},
|
|||
|
|
"outputs": [],
|
|||
|
|
"source": []
|
|||
|
|
}
|
|||
|
|
],
|
|||
|
|
"metadata": {
|
|||
|
|
"kernelspec": {
|
|||
|
|
"display_name": "Python 3",
|
|||
|
|
"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.6.8"
|
|||
|
|
}
|
|||
|
|
},
|
|||
|
|
"nbformat": 4,
|
|||
|
|
"nbformat_minor": 4
|
|||
|
|
}
|