{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import seaborn as sns\n", "from matplotlib import pyplot as plt\n", "from sklearn import preprocessing\n", "from scipy.signal import savgol_filter\n", "from scipy.signal import find_peaks" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array([['cyto', 'RFP', 'area'],\n", " ['cyto', 'RFP', 'cell_id'],\n", " ['cyto', 'RFP', 'convex_area'],\n", " ['cyto', 'RFP', 'cv_intensity'],\n", " ['cyto', 'RFP', 'eccentricity'],\n", " ['cyto', 'RFP', 'major_axis_length'],\n", " ['cyto', 'RFP', 'max_intensity'],\n", " ['cyto', 'RFP', 'mean_intensity'],\n", " ['cyto', 'RFP', 'median_intensity'],\n", " ['cyto', 'RFP', 'min_intensity'],\n", " ['cyto', 'RFP', 'minor_axis_length'],\n", " ['cyto', 'RFP', 'num_seg'],\n", " ['cyto', 'RFP', 'orientation'],\n", " ['cyto', 'RFP', 'parent'],\n", " ['cyto', 'RFP', 'perimeter'],\n", " ['cyto', 'RFP', 'solidity'],\n", " ['cyto', 'RFP', 'std_intensity'],\n", " ['cyto', 'RFP', 'total_intensity'],\n", " ['cyto', 'RFP', 'x'],\n", " ['cyto', 'RFP', 'y'],\n", " ['nuc', 'RFP', 'area'],\n", " ['nuc', 'RFP', 'cell_id'],\n", " ['nuc', 'RFP', 'convex_area'],\n", " ['nuc', 'RFP', 'cv_intensity'],\n", " ['nuc', 'RFP', 'eccentricity'],\n", " ['nuc', 'RFP', 'major_axis_length'],\n", " ['nuc', 'RFP', 'max_intensity'],\n", " ['nuc', 'RFP', 'mean_intensity'],\n", " ['nuc', 'RFP', 'median_intensity'],\n", " ['nuc', 'RFP', 'min_intensity'],\n", " ['nuc', 'RFP', 'minor_axis_length'],\n", " ['nuc', 'RFP', 'num_seg'],\n", " ['nuc', 'RFP', 'orientation'],\n", " ['nuc', 'RFP', 'parent'],\n", " ['nuc', 'RFP', 'perimeter'],\n", " ['nuc', 'RFP', 'solidity'],\n", " ['nuc', 'RFP', 'std_intensity'],\n", " ['nuc', 'RFP', 'total_intensity'],\n", " ['nuc', 'RFP', 'x'],\n", " ['nuc', 'RFP', 'y']], dtype=object)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "data = np.load(\"/Users/sakaikeiichirou/CellTK/output/20201217-PKAKTR-mating/df.npz\",allow_pickle=True)\n", "\n", "arr = data['arr']\n", "labels = data['labels']\n", "time = data['time']\n", "labels" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "cell = range(103)# 0-10 細胞\n", "label1 = np.all((labels == 'cyto') | (labels == 'RFP') | (labels == 'mean_intensity'), axis=1)\n", "label2 = np.all((labels == 'nuc') | (labels == 'RFP') | (labels == 'mean_intensity'), axis=1)\n", "label3 = np.all((labels == 'cyto') | (labels == 'RFP') | (labels == 'cell_id'), axis=1)\n", "label4 = np.all((labels == 'nuc') | (labels == 'RFP') | (labels == 'area'), axis=1)\n", "label5 = np.all((labels == 'nuc') | (labels == 'RFP') | (labels == 'total_intensity'), axis=1)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | nuc | \n", "cyto | \n", "cell_id_CellTK | \n", "nuc_area | \n", "nuc_intensity | \n", "c/n | \n", "cell_id | \n", "time | \n", "
---|---|---|---|---|---|---|---|---|
0 | \n", "631.379395 | \n", "652.070496 | \n", "1.0 | \n", "98.0 | \n", "61875.179688 | \n", "1.032771 | \n", "0 | \n", "0 | \n", "
1 | \n", "610.097473 | \n", "651.362244 | \n", "1.0 | \n", "98.0 | \n", "59789.554688 | \n", "1.067636 | \n", "0 | \n", "1 | \n", "
2 | \n", "670.208740 | \n", "611.098999 | \n", "1.0 | \n", "99.0 | \n", "66350.664062 | \n", "0.911804 | \n", "0 | \n", "2 | \n", "
3 | \n", "625.002075 | \n", "586.168396 | \n", "1.0 | \n", "99.0 | \n", "61875.203125 | \n", "0.937866 | \n", "0 | \n", "3 | \n", "
4 | \n", "614.701599 | \n", "595.280884 | \n", "1.0 | \n", "96.0 | \n", "59011.351562 | \n", "0.968406 | \n", "0 | \n", "4 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
44 | \n", "803.187744 | \n", "528.949463 | \n", "113.0 | \n", "44.0 | \n", "35340.261719 | \n", "0.658563 | \n", "102 | \n", "44 | \n", "
45 | \n", "794.455017 | \n", "650.435425 | \n", "113.0 | \n", "57.0 | \n", "45283.937500 | \n", "0.818719 | \n", "102 | \n", "45 | \n", "
46 | \n", "780.113647 | \n", "579.592834 | \n", "113.0 | \n", "54.0 | \n", "42126.136719 | \n", "0.742959 | \n", "102 | \n", "46 | \n", "
47 | \n", "760.558350 | \n", "640.652222 | \n", "113.0 | \n", "57.0 | \n", "43351.824219 | \n", "0.842345 | \n", "102 | \n", "47 | \n", "
48 | \n", "732.805237 | \n", "605.260620 | \n", "113.0 | \n", "60.0 | \n", "43968.312500 | \n", "0.825950 | \n", "102 | \n", "48 | \n", "
5047 rows × 8 columns
\n", "\n", " | nuc | \n", "cyto | \n", "cell_id_CellTK | \n", "nuc_area | \n", "nuc_intensity | \n", "c/n | \n", "cell_id | \n", "time | \n", "sav | \n", "
---|---|---|---|---|---|---|---|---|---|
0 | \n", "631.379395 | \n", "652.070496 | \n", "1.0 | \n", "98.0 | \n", "61875.179688 | \n", "1.032771 | \n", "0 | \n", "0 | \n", "97.600000 | \n", "
1 | \n", "610.097473 | \n", "651.362244 | \n", "1.0 | \n", "98.0 | \n", "59789.554688 | \n", "1.067636 | \n", "0 | \n", "1 | \n", "98.800000 | \n", "
2 | \n", "670.208740 | \n", "611.098999 | \n", "1.0 | \n", "99.0 | \n", "66350.664062 | \n", "0.911804 | \n", "0 | \n", "2 | \n", "99.000000 | \n", "
3 | \n", "625.002075 | \n", "586.168396 | \n", "1.0 | \n", "99.0 | \n", "61875.203125 | \n", "0.937866 | \n", "0 | \n", "3 | \n", "97.971429 | \n", "
4 | \n", "614.701599 | \n", "595.280884 | \n", "1.0 | \n", "96.0 | \n", "59011.351562 | \n", "0.968406 | \n", "0 | \n", "4 | \n", "97.628571 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
5042 | \n", "803.187744 | \n", "528.949463 | \n", "113.0 | \n", "44.0 | \n", "35340.261719 | \n", "0.658563 | \n", "102 | \n", "44 | \n", "NaN | \n", "
5043 | \n", "794.455017 | \n", "650.435425 | \n", "113.0 | \n", "57.0 | \n", "45283.937500 | \n", "0.818719 | \n", "102 | \n", "45 | \n", "52.885714 | \n", "
5044 | \n", "780.113647 | \n", "579.592834 | \n", "113.0 | \n", "54.0 | \n", "42126.136719 | \n", "0.742959 | \n", "102 | \n", "46 | \n", "56.400000 | \n", "
5045 | \n", "760.558350 | \n", "640.652222 | \n", "113.0 | \n", "57.0 | \n", "43351.824219 | \n", "0.842345 | \n", "102 | \n", "47 | \n", "58.600000 | \n", "
5046 | \n", "732.805237 | \n", "605.260620 | \n", "113.0 | \n", "60.0 | \n", "43968.312500 | \n", "0.825950 | \n", "102 | \n", "48 | \n", "58.800000 | \n", "
5047 rows × 9 columns
\n", "\n", " | nuc | \n", "cyto | \n", "cell_id_CellTK | \n", "nuc_area | \n", "nuc_intensity | \n", "c/n | \n", "cell_id | \n", "time | \n", "sav | \n", "
---|---|---|---|---|---|---|---|---|---|
245 | \n", "778.302612 | \n", "565.512756 | \n", "7.0 | \n", "68.0 | \n", "52924.578125 | \n", "0.726598 | \n", "5 | \n", "0 | \n", "NaN | \n", "
246 | \n", "890.961548 | \n", "627.026733 | \n", "7.0 | \n", "69.0 | \n", "61476.347656 | \n", "0.703764 | \n", "5 | \n", "1 | \n", "NaN | \n", "
247 | \n", "931.428467 | \n", "665.144043 | \n", "7.0 | \n", "72.0 | \n", "67062.851562 | \n", "0.714112 | \n", "5 | \n", "2 | \n", "68.742857 | \n", "
248 | \n", "940.747498 | \n", "604.788391 | \n", "7.0 | \n", "66.0 | \n", "62089.335938 | \n", "0.642881 | \n", "5 | \n", "3 | \n", "70.885714 | \n", "
249 | \n", "846.890015 | \n", "675.476746 | \n", "7.0 | \n", "78.0 | \n", "66057.421875 | \n", "0.797597 | \n", "5 | \n", "4 | \n", "74.057143 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
1220 | \n", "764.211487 | \n", "632.399658 | \n", "26.0 | \n", "75.0 | \n", "57315.863281 | \n", "0.827519 | \n", "24 | \n", "44 | \n", "74.400000 | \n", "
1221 | \n", "763.351135 | \n", "612.782715 | \n", "26.0 | \n", "76.0 | \n", "58014.687500 | \n", "0.802753 | \n", "24 | \n", "45 | \n", "77.457143 | \n", "
1222 | \n", "739.729065 | \n", "600.724548 | \n", "26.0 | \n", "79.0 | \n", "58438.597656 | \n", "0.812087 | \n", "24 | \n", "46 | \n", "75.657143 | \n", "
1223 | \n", "712.811890 | \n", "634.782104 | \n", "26.0 | \n", "72.0 | \n", "51322.457031 | \n", "0.890532 | \n", "24 | \n", "47 | \n", "NaN | \n", "
1224 | \n", "662.110840 | \n", "612.549438 | \n", "26.0 | \n", "82.0 | \n", "54293.089844 | \n", "0.925146 | \n", "24 | \n", "48 | \n", "NaN | \n", "
392 rows × 9 columns
\n", "\n", " | nuc | \n", "cyto | \n", "cell_id_CellTK | \n", "nuc_area | \n", "nuc_intensity | \n", "c/n | \n", "cell_id | \n", "time | \n", "sav | \n", "
---|---|---|---|---|---|---|---|---|---|
0 | \n", "778.302612 | \n", "565.512756 | \n", "7.0 | \n", "68.0 | \n", "52924.578125 | \n", "0.726598 | \n", "5 | \n", "0 | \n", "NaN | \n", "
1 | \n", "890.961548 | \n", "627.026733 | \n", "7.0 | \n", "69.0 | \n", "61476.347656 | \n", "0.703764 | \n", "5 | \n", "1 | \n", "NaN | \n", "
2 | \n", "931.428467 | \n", "665.144043 | \n", "7.0 | \n", "72.0 | \n", "67062.851562 | \n", "0.714112 | \n", "5 | \n", "2 | \n", "68.742857 | \n", "
3 | \n", "940.747498 | \n", "604.788391 | \n", "7.0 | \n", "66.0 | \n", "62089.335938 | \n", "0.642881 | \n", "5 | \n", "3 | \n", "70.885714 | \n", "
4 | \n", "846.890015 | \n", "675.476746 | \n", "7.0 | \n", "78.0 | \n", "66057.421875 | \n", "0.797597 | \n", "5 | \n", "4 | \n", "74.057143 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
387 | \n", "764.211487 | \n", "632.399658 | \n", "26.0 | \n", "75.0 | \n", "57315.863281 | \n", "0.827519 | \n", "24 | \n", "44 | \n", "74.400000 | \n", "
388 | \n", "763.351135 | \n", "612.782715 | \n", "26.0 | \n", "76.0 | \n", "58014.687500 | \n", "0.802753 | \n", "24 | \n", "45 | \n", "77.457143 | \n", "
389 | \n", "739.729065 | \n", "600.724548 | \n", "26.0 | \n", "79.0 | \n", "58438.597656 | \n", "0.812087 | \n", "24 | \n", "46 | \n", "75.657143 | \n", "
390 | \n", "712.811890 | \n", "634.782104 | \n", "26.0 | \n", "72.0 | \n", "51322.457031 | \n", "0.890532 | \n", "24 | \n", "47 | \n", "NaN | \n", "
391 | \n", "662.110840 | \n", "612.549438 | \n", "26.0 | \n", "82.0 | \n", "54293.089844 | \n", "0.925146 | \n", "24 | \n", "48 | \n", "NaN | \n", "
392 rows × 9 columns
\n", "\n", " | 0 | \n", "1 | \n", "2 | \n", "3 | \n", "4 | \n", "5 | \n", "6 | \n", "7 | \n", "8 | \n", "9 | \n", "... | \n", "40 | \n", "41 | \n", "42 | \n", "43 | \n", "44 | \n", "45 | \n", "46 | \n", "47 | \n", "48 | \n", "sum | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "0.949597 | \n", "0.892541 | \n", "0.969320 | \n", "0.948296 | \n", "0.954160 | \n", "0.934641 | \n", "1.001853 | \n", "1.059008 | \n", "1.069426 | \n", "1.017764 | \n", "... | \n", "0.959127 | \n", "1.038254 | \n", "1.218195 | \n", "1.187555 | \n", "1.181965 | \n", "1.251580 | \n", "1.245025 | \n", "1.357071 | \n", "1.330591 | \n", "50.437228 | \n", "
1 | \n", "1.018857 | \n", "0.998410 | \n", "1.011859 | \n", "1.044959 | \n", "1.009956 | \n", "0.992243 | \n", "1.024513 | \n", "1.058028 | \n", "1.019956 | \n", "0.968173 | \n", "... | \n", "1.263065 | \n", "1.226765 | \n", "1.155211 | \n", "1.211423 | \n", "1.195607 | \n", "1.145404 | \n", "1.205466 | \n", "1.297535 | \n", "1.214744 | \n", "50.393267 | \n", "
2 | \n", "0.726598 | \n", "0.703764 | \n", "0.714112 | \n", "0.642881 | \n", "0.797597 | \n", "0.845219 | \n", "0.906810 | \n", "0.903888 | \n", "1.029487 | \n", "1.003702 | \n", "... | \n", "1.102932 | \n", "1.165754 | \n", "1.206428 | \n", "1.141587 | \n", "1.206251 | \n", "1.193377 | \n", "1.194637 | \n", "1.095980 | \n", "1.097980 | \n", "46.666868 | \n", "
3 | \n", "0.922765 | \n", "0.781348 | \n", "0.817657 | \n", "0.820006 | \n", "0.872587 | \n", "0.884103 | \n", "0.917995 | \n", "0.939567 | \n", "0.942608 | \n", "0.922485 | \n", "... | \n", "0.817053 | \n", "0.994573 | \n", "1.043170 | \n", "1.117318 | \n", "1.122524 | \n", "1.059465 | \n", "0.974191 | \n", "1.006437 | \n", "1.043743 | \n", "45.110730 | \n", "
4 | \n", "1.032771 | \n", "1.067636 | \n", "0.911804 | \n", "0.937866 | \n", "0.968406 | \n", "0.952117 | \n", "1.055387 | \n", "1.024936 | \n", "1.016990 | \n", "0.999471 | \n", "... | \n", "0.957918 | \n", "0.833480 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "1.094948 | \n", "0.921935 | \n", "1.062213 | \n", "1.080750 | \n", "44.544500 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
98 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "... | \n", "NaN | \n", "NaN | \n", "NaN | \n", "1.193810 | \n", "0.794049 | \n", "0.834925 | \n", "0.798205 | \n", "0.801760 | \n", "0.748455 | \n", "5.171204 | \n", "
99 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "... | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "5.091612 | \n", "
100 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "... | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "4.797146 | \n", "
101 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "... | \n", "NaN | \n", "NaN | \n", "NaN | \n", "0.905515 | \n", "0.658563 | \n", "0.818719 | \n", "0.742959 | \n", "0.842345 | \n", "0.825950 | \n", "4.794050 | \n", "
102 | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "NaN | \n", "... | \n", "NaN | \n", "NaN | \n", "NaN | \n", "0.801815 | \n", "0.815597 | \n", "0.698161 | \n", "0.671207 | \n", "0.727087 | \n", "0.715472 | \n", "4.429339 | \n", "
103 rows × 50 columns
\n", "\n", " | 0 | \n", "1 | \n", "2 | \n", "3 | \n", "4 | \n", "5 | \n", "6 | \n", "7 | \n", "8 | \n", "9 | \n", "... | \n", "40 | \n", "41 | \n", "42 | \n", "43 | \n", "44 | \n", "45 | \n", "46 | \n", "47 | \n", "48 | \n", "sum | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "0.949597 | \n", "0.892541 | \n", "0.969320 | \n", "0.948296 | \n", "0.954160 | \n", "0.934641 | \n", "1.001853 | \n", "1.059008 | \n", "1.069426 | \n", "1.017764 | \n", "... | \n", "0.959127 | \n", "1.038254 | \n", "1.218195 | \n", "1.187555 | \n", "1.181965 | \n", "1.251580 | \n", "1.245025 | \n", "1.357071 | \n", "1.330591 | \n", "50.437228 | \n", "
1 | \n", "1.018857 | \n", "0.998410 | \n", "1.011859 | \n", "1.044959 | \n", "1.009956 | \n", "0.992243 | \n", "1.024513 | \n", "1.058028 | \n", "1.019956 | \n", "0.968173 | \n", "... | \n", "1.263065 | \n", "1.226765 | \n", "1.155211 | \n", "1.211423 | \n", "1.195607 | \n", "1.145404 | \n", "1.205466 | \n", "1.297535 | \n", "1.214744 | \n", "50.393267 | \n", "
2 | \n", "0.726598 | \n", "0.703764 | \n", "0.714112 | \n", "0.642881 | \n", "0.797597 | \n", "0.845219 | \n", "0.906810 | \n", "0.903888 | \n", "1.029487 | \n", "1.003702 | \n", "... | \n", "1.102932 | \n", "1.165754 | \n", "1.206428 | \n", "1.141587 | \n", "1.206251 | \n", "1.193377 | \n", "1.194637 | \n", "1.095980 | \n", "1.097980 | \n", "46.666868 | \n", "
3 | \n", "0.922765 | \n", "0.781348 | \n", "0.817657 | \n", "0.820006 | \n", "0.872587 | \n", "0.884103 | \n", "0.917995 | \n", "0.939567 | \n", "0.942608 | \n", "0.922485 | \n", "... | \n", "0.817053 | \n", "0.994573 | \n", "1.043170 | \n", "1.117318 | \n", "1.122524 | \n", "1.059465 | \n", "0.974191 | \n", "1.006437 | \n", "1.043743 | \n", "45.110730 | \n", "
4 | \n", "1.077662 | \n", "1.138995 | \n", "1.042852 | \n", "1.120010 | \n", "1.126836 | \n", "1.000584 | \n", "0.956947 | \n", "1.022942 | \n", "0.927019 | \n", "0.899393 | \n", "... | \n", "0.899262 | \n", "0.909205 | \n", "1.009297 | \n", "0.906015 | \n", "0.817220 | \n", "0.712576 | \n", "0.902584 | \n", "0.813840 | \n", "0.807195 | \n", "44.254944 | \n", "
5 | \n", "1.130260 | \n", "1.076677 | \n", "1.006540 | \n", "1.008871 | \n", "1.029734 | \n", "0.910807 | \n", "0.907097 | \n", "0.886559 | \n", "0.860312 | \n", "0.898637 | \n", "... | \n", "0.904958 | \n", "0.751175 | \n", "0.756242 | \n", "0.786730 | \n", "0.907457 | \n", "0.804729 | \n", "0.984117 | \n", "1.026825 | \n", "0.979290 | \n", "44.145890 | \n", "
6 | \n", "0.836502 | \n", "0.871300 | \n", "0.914461 | \n", "0.890506 | \n", "0.916582 | \n", "0.924223 | \n", "0.813281 | \n", "0.933831 | \n", "0.886805 | \n", "0.921183 | \n", "... | \n", "0.922378 | \n", "0.891434 | \n", "0.746491 | \n", "0.842014 | \n", "0.825187 | \n", "0.912558 | \n", "0.982578 | \n", "0.897389 | \n", "0.934093 | \n", "42.164062 | \n", "
7 | \n", "0.623727 | \n", "0.695840 | \n", "0.529479 | \n", "0.572524 | \n", "0.564623 | \n", "0.624924 | \n", "0.819714 | \n", "0.817189 | \n", "0.959135 | \n", "1.009863 | \n", "... | \n", "0.746295 | \n", "0.810380 | \n", "0.771778 | \n", "0.786081 | \n", "0.827519 | \n", "0.802753 | \n", "0.812087 | \n", "0.890532 | \n", "0.925146 | \n", "38.513013 | \n", "
8 rows × 50 columns
\n", "