Skip to content

M
mooc-rr
Commit 0326bd24 authored by cb1e5ba91280d02583d0604166e71c1c's avatar cb1e5ba91280d02583d0604166e71c1c
Browse files
Options

no commit message

parent 28899bba
Show whitespace changes
Inline Side-by-side
Showing with 218 additions and 1261 deletions
module3/exo3/exercice.ipynb
View file @ 0326bd24
......@@ -9,7 +9,7 @@
},
{
"cell_type": "code",
"execution_count": 27,
"execution_count": 60,
"metadata": {},
"outputs": [],
"source": [
......@@ -34,7 +34,7 @@
},
{
"cell_type": "code",
"execution_count": 29,
"execution_count": 61,
"metadata": {},
"outputs": [
{
......@@ -156,7 +156,7 @@
"4 315.06 317.51 314.71 "
]
},
"execution_count": 29,
"execution_count": 61,
"metadata": {},
"output_type": "execute_result"
}
......@@ -175,8 +175,7 @@
"else:\n",
" raw_data.to_csv('monthly_in_situ_co2_mlo',index=True)\n",
"\n",
"raw_data.head(5)\n",
"#raw_data.tail(5)"
"raw_data.head(5)"
]
},
{
......@@ -188,18 +187,28 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": 62,
"metadata": {},
"outputs": [],
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Index([' ', ' ', ' Excel', ' ', ' [ppm]', ' [ppm] ',\n",
" ' [ppm]', ' [ppm]', ' [ppm]', ' [ppm].1'],\n",
" dtype='object')\n"
]
}
],
"source": [
"raw_data.columns\n",
"raw_data.columns = ['Yr','Mn','Date 1','Date 2','s1','s2','s3','s4','s5','s6']\n",
"data = raw_data"
"raw_data_1 = raw_data.copy()\n",
"print(raw_data.columns)\n",
"raw_data_1.columns = ['Yr','Mn','Date 1','Date 2','s1','s2','s3','s4','s5','s6']"
]
},
{
"cell_type": "code",
"execution_count": 17,
"execution_count": 63,
"metadata": {},
"outputs": [
{
......@@ -237,102 +246,11 @@
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>738</th>\n",
" <td>2019</td>\n",
" <td>7</td>\n",
" <td>43661</td>\n",
" <td>2019.5370</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" <td>412.36</td>\n",
" <td>411.57</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>739</th>\n",
" <td>2019</td>\n",
" <td>8</td>\n",
" <td>43692</td>\n",
" <td>2019.6219</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" <td>410.23</td>\n",
" <td>411.81</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" </tr>\n",
" <tr>\n",
" <th>740</th>\n",
" <td>2019</td>\n",
" <td>9</td>\n",
" <td>43723</td>\n",
" <td>2019.7068</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" <td>408.52</td>\n",
" <td>412.05</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>741</th>\n",
" <td>2019</td>\n",
" <td>10</td>\n",
" <td>43753</td>\n",
" <td>2019.7890</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" <td>408.67</td>\n",
" <td>412.29</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" </tr>\n",
" <tr>\n",
" <th>742</th>\n",
" <td>2019</td>\n",
" <td>11</td>\n",
" <td>43784</td>\n",
" <td>2019.8740</td>\n",
" <td>410.28</td>\n",
" <td>412.56</td>\n",
" <td>410.29</td>\n",
" <td>412.54</td>\n",
" <td>410.28</td>\n",
" <td>412.56</td>\n",
" </tr>\n",
" <tr>\n",
" <th>743</th>\n",
" <td>2019</td>\n",
" <td>12</td>\n",
" <td>43814</td>\n",
" <td>2019.9562</td>\n",
" <td>411.85</td>\n",
" <td>412.78</td>\n",
" <td>411.88</td>\n",
" <td>412.79</td>\n",
" <td>411.85</td>\n",
" <td>412.78</td>\n",
" </tr>\n",
" <tr>\n",
" <th>744</th>\n",
" <td>2020</td>\n",
" <th>0</th>\n",
" <td>1958</td>\n",
" <td>1</td>\n",
" <td>43845</td>\n",
" <td>2020.0410</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" <td>413.11</td>\n",
" <td>413.05</td>\n",
" <td>413.37</td>\n",
" <td>413.33</td>\n",
" </tr>\n",
" <tr>\n",
" <th>745</th>\n",
" <td>2020</td>\n",
" <td>2</td>\n",
" <td>43876</td>\n",
" <td>2020.1257</td>\n",
" <td>21200</td>\n",
" <td>1958.0411</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
......@@ -341,11 +259,11 @@
" <td>-99.99</td>\n",
" </tr>\n",
" <tr>\n",
" <th>746</th>\n",
" <td>2020</td>\n",
" <td>3</td>\n",
" <td>43905</td>\n",
" <td>2020.2049</td>\n",
" <th>1</th>\n",
" <td>1958</td>\n",
" <td>2</td>\n",
" <td>21231</td>\n",
" <td>1958.1260</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
......@@ -354,105 +272,76 @@
" <td>-99.99</td>\n",
" </tr>\n",
" <tr>\n",
" <th>747</th>\n",
" <td>2020</td>\n",
" <th>2</th>\n",
" <td>1958</td>\n",
" <td>3</td>\n",
" <td>21259</td>\n",
" <td>1958.2027</td>\n",
" <td>315.70</td>\n",
" <td>314.44</td>\n",
" <td>316.19</td>\n",
" <td>314.91</td>\n",
" <td>315.70</td>\n",
" <td>314.44</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>1958</td>\n",
" <td>4</td>\n",
" <td>43936</td>\n",
" <td>2020.2896</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>-99.99</td>\n",
" <td>21290</td>\n",
" <td>1958.2877</td>\n",
" <td>317.45</td>\n",
" <td>315.16</td>\n",
" <td>317.30</td>\n",
" <td>314.99</td>\n",
" <td>317.45</td>\n",
" <td>315.16</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>1958</td>\n",
" <td>5</td>\n",
" <td>21320</td>\n",
" <td>1958.3699</td>\n",
" <td>317.51</td>\n",
" <td>314.71</td>\n",
" <td>317.86</td>\n",
" <td>315.06</td>\n",
" <td>317.51</td>\n",
" <td>314.71</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 \\\n",
"738 2019 7 43661 2019.5370 411.85 411.03 412.36 411.57 411.85 \n",
"739 2019 8 43692 2019.6219 410.08 411.62 410.23 411.81 410.08 \n",
"740 2019 9 43723 2019.7068 408.55 412.06 408.52 412.05 408.55 \n",
"741 2019 10 43753 2019.7890 408.43 412.07 408.67 412.29 408.43 \n",
"742 2019 11 43784 2019.8740 410.28 412.56 410.29 412.54 410.28 \n",
"743 2019 12 43814 2019.9562 411.85 412.78 411.88 412.79 411.85 \n",
"744 2020 1 43845 2020.0410 413.37 413.33 413.11 413.05 413.37 \n",
"745 2020 2 43876 2020.1257 -99.99 -99.99 -99.99 -99.99 -99.99 \n",
"746 2020 3 43905 2020.2049 -99.99 -99.99 -99.99 -99.99 -99.99 \n",
"747 2020 4 43936 2020.2896 -99.99 -99.99 -99.99 -99.99 -99.99 \n",
"\n",
" s6 \n",
"738 411.03 \n",
"739 411.62 \n",
"740 412.06 \n",
"741 412.07 \n",
"742 412.56 \n",
"743 412.78 \n",
"744 413.33 \n",
"745 -99.99 \n",
"746 -99.99 \n",
"747 -99.99 "
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 s6\n",
"0 1958 1 21200 1958.0411 -99.99 -99.99 -99.99 -99.99 -99.99 -99.99\n",
"1 1958 2 21231 1958.1260 -99.99 -99.99 -99.99 -99.99 -99.99 -99.99\n",
"2 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 315.70 314.44\n",
"3 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 317.45 315.16\n",
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 314.71"
]
},
"execution_count": 17,
"execution_count": 63,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data.tail(10)"
"raw_data_1.head()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Les données vide sont représentés par la valeur -99.99, nous remplaçons cette valeur par une valeur plus adéquate NaN dans une autre DataFrame\n",
"while k>=i:\n",
" data.loc[k]\n",
" if data.loc[k].isnull().any():\n",
" print(k)\n",
" print(data.loc[k])\n",
" k = k -1"
]
},
{
"cell_type": "code",
"execution_count": 32,
"metadata": {},
"outputs": [],
"source": [
"date = datetime.date.today()\n",
"data = data.drop(data[(data.Yr == date.year) & (data.Mn > date.month)].index)\n",
"data = data.replace(-99.99,np.NaN);\n",
"d = data[(data.Yr == date.year)].index\n",
"i,k = d[0], d[-1]\n",
"\n",
"while k>=i:\n",
" if data.loc[k].isnull().any():\n",
" data = data.drop(k)\n",
" else:\n",
" break\n",
" k = k - 1\n",
" \n",
"annee, mois = data.Yr[0], data.Mn[0]\n",
"d = data[(data.Yr == annee)].index\n",
"i,k = d[0], d[-1]\n",
"\n",
"while i<=k:\n",
" if data.loc[i].isnull().any():\n",
" data = data.drop(i)\n",
" else:\n",
" break\n",
" i = i + 1\n",
"\n"
"Les données vide sont représentés par la valeur -99.99, nous remplaçons cette valeur par une valeur plus adéquate NaN dans une autre DataFrame"
]
},
{
"cell_type": "code",
"execution_count": 33,
"execution_count": 71,
"metadata": {},
"outputs": [
{
......@@ -528,721 +417,55 @@
" <td>317.51</td>\n",
" <td>314.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
" <td>21351</td>\n",
" <td>1958.4548</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>1958</td>\n",
" <td>7</td>\n",
" <td>21381</td>\n",
" <td>1958.5370</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" <td>315.86</td>\n",
" <td>315.22</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>7</th>\n",
" <td>1958</td>\n",
" <td>8</td>\n",
" <td>21412</td>\n",
" <td>1958.6219</td>\n",
" <td>314.93</td>\n",
" <td>316.19</td>\n",
" <td>314.00</td>\n",
" <td>315.29</td>\n",
" <td>314.93</td>\n",
" <td>316.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>8</th>\n",
" <td>1958</td>\n",
" <td>9</td>\n",
" <td>21443</td>\n",
" <td>1958.7068</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" <td>312.46</td>\n",
" <td>315.35</td>\n",
" <td>313.21</td>\n",
" <td>316.08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>9</th>\n",
" <td>1958</td>\n",
" <td>10</td>\n",
" <td>21473</td>\n",
" <td>1958.7890</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>312.44</td>\n",
" <td>315.40</td>\n",
" <td>312.44</td>\n",
" <td>315.40</td>\n",
" </tr>\n",
" <tr>\n",
" <th>10</th>\n",
" <td>1958</td>\n",
" <td>11</td>\n",
" <td>21504</td>\n",
" <td>1958.8740</td>\n",
" <td>313.33</td>\n",
" <td>315.20</td>\n",
" <td>313.62</td>\n",
" <td>315.46</td>\n",
" <td>313.33</td>\n",
" <td>315.20</td>\n",
" </tr>\n",
" <tr>\n",
" <th>11</th>\n",
" <td>1958</td>\n",
" <td>12</td>\n",
" <td>21534</td>\n",
" <td>1958.9562</td>\n",
" <td>314.67</td>\n",
" <td>315.43</td>\n",
" <td>314.77</td>\n",
" <td>315.51</td>\n",
" <td>314.67</td>\n",
" <td>315.43</td>\n",
" </tr>\n",
" <tr>\n",
" <th>12</th>\n",
" <td>1959</td>\n",
" <td>1</td>\n",
" <td>21565</td>\n",
" <td>1959.0411</td>\n",
" <td>315.58</td>\n",
" <td>315.54</td>\n",
" <td>315.62</td>\n",
" <td>315.57</td>\n",
" <td>315.58</td>\n",
" <td>315.54</td>\n",
" </tr>\n",
" <tr>\n",
" <th>13</th>\n",
" <td>1959</td>\n",
" <td>2</td>\n",
" <td>21596</td>\n",
" <td>1959.1260</td>\n",
" <td>316.49</td>\n",
" <td>315.86</td>\n",
" <td>316.27</td>\n",
" <td>315.63</td>\n",
" <td>316.49</td>\n",
" <td>315.86</td>\n",
" </tr>\n",
" <tr>\n",
" <th>14</th>\n",
" <td>1959</td>\n",
" <td>3</td>\n",
" <td>21624</td>\n",
" <td>1959.2027</td>\n",
" <td>316.65</td>\n",
" <td>315.38</td>\n",
" <td>316.98</td>\n",
" <td>315.69</td>\n",
" <td>316.65</td>\n",
" <td>315.38</td>\n",
" </tr>\n",
" <tr>\n",
" <th>15</th>\n",
" <td>1959</td>\n",
" <td>4</td>\n",
" <td>21655</td>\n",
" <td>1959.2877</td>\n",
" <td>317.72</td>\n",
" <td>315.42</td>\n",
" <td>318.09</td>\n",
" <td>315.77</td>\n",
" <td>317.72</td>\n",
" <td>315.42</td>\n",
" </tr>\n",
" <tr>\n",
" <th>16</th>\n",
" <td>1959</td>\n",
" <td>5</td>\n",
" <td>21685</td>\n",
" <td>1959.3699</td>\n",
" <td>318.29</td>\n",
" <td>315.49</td>\n",
" <td>318.65</td>\n",
" <td>315.85</td>\n",
" <td>318.29</td>\n",
" <td>315.49</td>\n",
" </tr>\n",
" <tr>\n",
" <th>17</th>\n",
" <td>1959</td>\n",
" <td>6</td>\n",
" <td>21716</td>\n",
" <td>1959.4548</td>\n",
" <td>318.15</td>\n",
" <td>316.03</td>\n",
" <td>318.04</td>\n",
" <td>315.94</td>\n",
" <td>318.15</td>\n",
" <td>316.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>18</th>\n",
" <td>1959</td>\n",
" <td>7</td>\n",
" <td>21746</td>\n",
" <td>1959.5370</td>\n",
" <td>316.54</td>\n",
" <td>315.86</td>\n",
" <td>316.67</td>\n",
" <td>316.03</td>\n",
" <td>316.54</td>\n",
" <td>315.86</td>\n",
" </tr>\n",
" <tr>\n",
" <th>19</th>\n",
" <td>1959</td>\n",
" <td>8</td>\n",
" <td>21777</td>\n",
" <td>1959.6219</td>\n",
" <td>314.80</td>\n",
" <td>316.06</td>\n",
" <td>314.83</td>\n",
" <td>316.12</td>\n",
" <td>314.80</td>\n",
" <td>316.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>20</th>\n",
" <td>1959</td>\n",
" <td>9</td>\n",
" <td>21808</td>\n",
" <td>1959.7068</td>\n",
" <td>313.84</td>\n",
" <td>316.72</td>\n",
" <td>313.32</td>\n",
" <td>316.22</td>\n",
" <td>313.84</td>\n",
" <td>316.72</td>\n",
" </tr>\n",
" <tr>\n",
" <th>21</th>\n",
" <td>1959</td>\n",
" <td>10</td>\n",
" <td>21838</td>\n",
" <td>1959.7890</td>\n",
" <td>313.33</td>\n",
" <td>316.32</td>\n",
" <td>313.33</td>\n",
" <td>316.30</td>\n",
" <td>313.33</td>\n",
" <td>316.32</td>\n",
" </tr>\n",
" <tr>\n",
" <th>22</th>\n",
" <td>1959</td>\n",
" <td>11</td>\n",
" <td>21869</td>\n",
" <td>1959.8740</td>\n",
" <td>314.81</td>\n",
" <td>316.68</td>\n",
" <td>314.54</td>\n",
" <td>316.39</td>\n",
" <td>314.81</td>\n",
" <td>316.68</td>\n",
" </tr>\n",
" <tr>\n",
" <th>23</th>\n",
" <td>1959</td>\n",
" <td>12</td>\n",
" <td>21899</td>\n",
" <td>1959.9562</td>\n",
" <td>315.58</td>\n",
" <td>316.35</td>\n",
" <td>315.72</td>\n",
" <td>316.47</td>\n",
" <td>315.58</td>\n",
" <td>316.35</td>\n",
" </tr>\n",
" <tr>\n",
" <th>24</th>\n",
" <td>1960</td>\n",
" <td>1</td>\n",
" <td>21930</td>\n",
" <td>1960.0410</td>\n",
" <td>316.43</td>\n",
" <td>316.39</td>\n",
" <td>316.61</td>\n",
" <td>316.56</td>\n",
" <td>316.43</td>\n",
" <td>316.39</td>\n",
" </tr>\n",
" <tr>\n",
" <th>25</th>\n",
" <td>1960</td>\n",
" <td>2</td>\n",
" <td>21961</td>\n",
" <td>1960.1257</td>\n",
" <td>316.98</td>\n",
" <td>316.35</td>\n",
" <td>317.27</td>\n",
" <td>316.64</td>\n",
" <td>316.98</td>\n",
" <td>316.35</td>\n",
" </tr>\n",
" <tr>\n",
" <th>26</th>\n",
" <td>1960</td>\n",
" <td>3</td>\n",
" <td>21990</td>\n",
" <td>1960.2049</td>\n",
" <td>317.58</td>\n",
" <td>316.28</td>\n",
" <td>318.03</td>\n",
" <td>316.71</td>\n",
" <td>317.58</td>\n",
" <td>316.28</td>\n",
" </tr>\n",
" <tr>\n",
" <th>27</th>\n",
" <td>1960</td>\n",
" <td>4</td>\n",
" <td>22021</td>\n",
" <td>1960.2896</td>\n",
" <td>319.03</td>\n",
" <td>316.70</td>\n",
" <td>319.14</td>\n",
" <td>316.79</td>\n",
" <td>319.03</td>\n",
" <td>316.70</td>\n",
" </tr>\n",
" <tr>\n",
" <th>28</th>\n",
" <td>1960</td>\n",
" <td>5</td>\n",
" <td>22051</td>\n",
" <td>1960.3716</td>\n",
" <td>320.04</td>\n",
" <td>317.22</td>\n",
" <td>319.67</td>\n",
" <td>316.86</td>\n",
" <td>320.04</td>\n",
" <td>317.22</td>\n",
" </tr>\n",
" <tr>\n",
" <th>29</th>\n",
" <td>1960</td>\n",
" <td>6</td>\n",
" <td>22082</td>\n",
" <td>1960.4563</td>\n",
" <td>319.59</td>\n",
" <td>317.48</td>\n",
" <td>319.01</td>\n",
" <td>316.93</td>\n",
" <td>319.59</td>\n",
" <td>317.48</td>\n",
" </tr>\n",
" <tr>\n",
" <th>30</th>\n",
" <td>1960</td>\n",
" <td>7</td>\n",
" <td>22112</td>\n",
" <td>1960.5383</td>\n",
" <td>318.18</td>\n",
" <td>317.53</td>\n",
" <td>317.60</td>\n",
" <td>316.98</td>\n",
" <td>318.18</td>\n",
" <td>317.53</td>\n",
" </tr>\n",
" <tr>\n",
" <th>31</th>\n",
" <td>1960</td>\n",
" <td>8</td>\n",
" <td>22143</td>\n",
" <td>1960.6230</td>\n",
" <td>315.90</td>\n",
" <td>317.20</td>\n",
" <td>315.69</td>\n",
" <td>317.01</td>\n",
" <td>315.90</td>\n",
" <td>317.20</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>715</th>\n",
" <td>2017</td>\n",
" <td>8</td>\n",
" <td>42962</td>\n",
" <td>2017.6219</td>\n",
" <td>405.24</td>\n",
" <td>406.77</td>\n",
" <td>405.20</td>\n",
" <td>406.76</td>\n",
" <td>405.24</td>\n",
" <td>406.77</td>\n",
" </tr>\n",
" <tr>\n",
" <th>716</th>\n",
" <td>2017</td>\n",
" <td>9</td>\n",
" <td>42993</td>\n",
" <td>2017.7068</td>\n",
" <td>403.27</td>\n",
" <td>406.76</td>\n",
" <td>403.41</td>\n",
" <td>406.92</td>\n",
" <td>403.27</td>\n",
" <td>406.76</td>\n",
" </tr>\n",
" <tr>\n",
" <th>717</th>\n",
" <td>2017</td>\n",
" <td>10</td>\n",
" <td>43023</td>\n",
" <td>2017.7890</td>\n",
" <td>403.64</td>\n",
" <td>407.25</td>\n",
" <td>403.47</td>\n",
" <td>407.07</td>\n",
" <td>403.64</td>\n",
" <td>407.25</td>\n",
" </tr>\n",
" <tr>\n",
" <th>718</th>\n",
" <td>2017</td>\n",
" <td>11</td>\n",
" <td>43054</td>\n",
" <td>2017.8740</td>\n",
" <td>405.17</td>\n",
" <td>407.43</td>\n",
" <td>404.98</td>\n",
" <td>407.22</td>\n",
" <td>405.17</td>\n",
" <td>407.43</td>\n",
" </tr>\n",
" <tr>\n",
" <th>719</th>\n",
" <td>2017</td>\n",
" <td>12</td>\n",
" <td>43084</td>\n",
" <td>2017.9562</td>\n",
" <td>406.75</td>\n",
" <td>407.68</td>\n",
" <td>406.46</td>\n",
" <td>407.37</td>\n",
" <td>406.75</td>\n",
" <td>407.68</td>\n",
" </tr>\n",
" <tr>\n",
" <th>720</th>\n",
" <td>2018</td>\n",
" <td>1</td>\n",
" <td>43115</td>\n",
" <td>2018.0411</td>\n",
" <td>408.05</td>\n",
" <td>408.00</td>\n",
" <td>407.58</td>\n",
" <td>407.52</td>\n",
" <td>408.05</td>\n",
" <td>408.00</td>\n",
" </tr>\n",
" <tr>\n",
" <th>721</th>\n",
" <td>2018</td>\n",
" <td>2</td>\n",
" <td>43146</td>\n",
" <td>2018.1260</td>\n",
" <td>408.34</td>\n",
" <td>407.59</td>\n",
" <td>408.44</td>\n",
" <td>407.68</td>\n",
" <td>408.34</td>\n",
" <td>407.59</td>\n",
" </tr>\n",
" <tr>\n",
" <th>722</th>\n",
" <td>2018</td>\n",
" <td>3</td>\n",
" <td>43174</td>\n",
" <td>2018.2027</td>\n",
" <td>409.25</td>\n",
" <td>407.71</td>\n",
" <td>409.39</td>\n",
" <td>407.82</td>\n",
" <td>409.25</td>\n",
" <td>407.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>723</th>\n",
" <td>2018</td>\n",
" <td>4</td>\n",
" <td>43205</td>\n",
" <td>2018.2877</td>\n",
" <td>410.30</td>\n",
" <td>407.51</td>\n",
" <td>410.82</td>\n",
" <td>408.00</td>\n",
" <td>410.30</td>\n",
" <td>407.51</td>\n",
" </tr>\n",
" <tr>\n",
" <th>724</th>\n",
" <td>2018</td>\n",
" <td>5</td>\n",
" <td>43235</td>\n",
" <td>2018.3699</td>\n",
" <td>411.30</td>\n",
" <td>407.90</td>\n",
" <td>411.60</td>\n",
" <td>408.20</td>\n",
" <td>411.30</td>\n",
" <td>407.90</td>\n",
" </tr>\n",
" <tr>\n",
" <th>725</th>\n",
" <td>2018</td>\n",
" <td>6</td>\n",
" <td>43266</td>\n",
" <td>2018.4548</td>\n",
" <td>410.88</td>\n",
" <td>408.30</td>\n",
" <td>410.97</td>\n",
" <td>408.42</td>\n",
" <td>410.88</td>\n",
" <td>408.30</td>\n",
" </tr>\n",
" <tr>\n",
" <th>726</th>\n",
" <td>2018</td>\n",
" <td>7</td>\n",
" <td>43296</td>\n",
" <td>2018.5370</td>\n",
" <td>408.90</td>\n",
" <td>408.08</td>\n",
" <td>409.43</td>\n",
" <td>408.65</td>\n",
" <td>408.90</td>\n",
" <td>408.08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>727</th>\n",
" <td>2018</td>\n",
" <td>8</td>\n",
" <td>43327</td>\n",
" <td>2018.6219</td>\n",
" <td>407.10</td>\n",
" <td>408.63</td>\n",
" <td>407.33</td>\n",
" <td>408.91</td>\n",
" <td>407.10</td>\n",
" <td>408.63</td>\n",
" </tr>\n",
" <tr>\n",
" <th>728</th>\n",
" <td>2018</td>\n",
" <td>9</td>\n",
" <td>43358</td>\n",
" <td>2018.7068</td>\n",
" <td>405.59</td>\n",
" <td>409.09</td>\n",
" <td>405.66</td>\n",
" <td>409.18</td>\n",
" <td>405.59</td>\n",
" <td>409.09</td>\n",
" </tr>\n",
" <tr>\n",
" <th>729</th>\n",
" <td>2018</td>\n",
" <td>10</td>\n",
" <td>43388</td>\n",
" <td>2018.7890</td>\n",
" <td>405.99</td>\n",
" <td>409.62</td>\n",
" <td>405.83</td>\n",
" <td>409.44</td>\n",
" <td>405.99</td>\n",
" <td>409.62</td>\n",
" </tr>\n",
" <tr>\n",
" <th>730</th>\n",
" <td>2018</td>\n",
" <td>11</td>\n",
" <td>43419</td>\n",
" <td>2018.8740</td>\n",
" <td>408.12</td>\n",
" <td>410.39</td>\n",
" <td>407.47</td>\n",
" <td>409.72</td>\n",
" <td>408.12</td>\n",
" <td>410.39</td>\n",
" </tr>\n",
" <tr>\n",
" <th>731</th>\n",
" <td>2018</td>\n",
" <td>12</td>\n",
" <td>43449</td>\n",
" <td>2018.9562</td>\n",
" <td>409.23</td>\n",
" <td>410.16</td>\n",
" <td>409.07</td>\n",
" <td>409.97</td>\n",
" <td>409.23</td>\n",
" <td>410.16</td>\n",
" </tr>\n",
" <tr>\n",
" <th>732</th>\n",
" <td>2019</td>\n",
" <td>1</td>\n",
" <td>43480</td>\n",
" <td>2019.0411</td>\n",
" <td>410.92</td>\n",
" <td>410.87</td>\n",
" <td>410.29</td>\n",
" <td>410.23</td>\n",
" <td>410.92</td>\n",
" <td>410.87</td>\n",
" </tr>\n",
" <tr>\n",
" <th>733</th>\n",
" <td>2019</td>\n",
" <td>2</td>\n",
" <td>43511</td>\n",
" <td>2019.1260</td>\n",
" <td>411.66</td>\n",
" <td>410.90</td>\n",
" <td>411.24</td>\n",
" <td>410.47</td>\n",
" <td>411.66</td>\n",
" <td>410.90</td>\n",
" </tr>\n",
" <tr>\n",
" <th>734</th>\n",
" <td>2019</td>\n",
" <td>3</td>\n",
" <td>43539</td>\n",
" <td>2019.2027</td>\n",
" <td>412.00</td>\n",
" <td>410.45</td>\n",
" <td>412.25</td>\n",
" <td>410.68</td>\n",
" <td>412.00</td>\n",
" <td>410.45</td>\n",
" </tr>\n",
" <tr>\n",
" <th>735</th>\n",
" <td>2019</td>\n",
" <td>4</td>\n",
" <td>43570</td>\n",
" <td>2019.2877</td>\n",
" <td>413.52</td>\n",
" <td>410.72</td>\n",
" <td>413.73</td>\n",
" <td>410.91</td>\n",
" <td>413.52</td>\n",
" <td>410.72</td>\n",
" </tr>\n",
" <tr>\n",
" <th>736</th>\n",
" <td>2019</td>\n",
" <td>5</td>\n",
" <td>43600</td>\n",
" <td>2019.3699</td>\n",
" <td>414.83</td>\n",
" <td>411.42</td>\n",
" <td>414.54</td>\n",
" <td>411.13</td>\n",
" <td>414.83</td>\n",
" <td>411.42</td>\n",
" </tr>\n",
" <tr>\n",
" <th>737</th>\n",
" <td>2019</td>\n",
" <td>6</td>\n",
" <td>43631</td>\n",
" <td>2019.4548</td>\n",
" <td>413.96</td>\n",
" <td>411.38</td>\n",
" <td>413.91</td>\n",
" <td>411.35</td>\n",
" <td>413.96</td>\n",
" <td>411.38</td>\n",
" </tr>\n",
" <tr>\n",
" <th>738</th>\n",
" <td>2019</td>\n",
" <td>7</td>\n",
" <td>43661</td>\n",
" <td>2019.5370</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" <td>412.36</td>\n",
" <td>411.57</td>\n",
" <td>411.85</td>\n",
" <td>411.03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>739</th>\n",
" <td>2019</td>\n",
" <td>8</td>\n",
" <td>43692</td>\n",
" <td>2019.6219</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" <td>410.23</td>\n",
" <td>411.81</td>\n",
" <td>410.08</td>\n",
" <td>411.62</td>\n",
" </tr>\n",
" <tr>\n",
" <th>740</th>\n",
" <td>2019</td>\n",
" <td>9</td>\n",
" <td>43723</td>\n",
" <td>2019.7068</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" <td>408.52</td>\n",
" <td>412.05</td>\n",
" <td>408.55</td>\n",
" <td>412.06</td>\n",
" </tr>\n",
" <tr>\n",
" <th>741</th>\n",
" <td>2019</td>\n",
" <td>10</td>\n",
" <td>43753</td>\n",
" <td>2019.7890</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" <td>408.67</td>\n",
" <td>412.29</td>\n",
" <td>408.43</td>\n",
" <td>412.07</td>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 s6\n",
"2 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 315.70 314.44\n",
"3 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 317.45 315.16\n",
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 314.71"
]
},
"execution_count": 71,
"metadata": {},
"output_type": "execute_result"
},
{
"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>Yr</th>\n",
" <th>Mn</th>\n",
" <th>Date 1</th>\n",
" <th>Date 2</th>\n",
" <th>s1</th>\n",
" <th>s2</th>\n",
" <th>s3</th>\n",
" <th>s4</th>\n",
" <th>s5</th>\n",
" <th>s6</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>742</th>\n",
" <td>2019</td>\n",
......@@ -1284,147 +507,54 @@
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>743 rows × 10 columns</p>\n",
"</div>"
],
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 \\\n",
"2 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 315.70 \n",
"3 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 317.45 \n",
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 \n",
"5 1958 6 21351 1958.4548 NaN NaN 317.24 315.14 317.24 \n",
"6 1958 7 21381 1958.5370 315.86 315.19 315.86 315.22 315.86 \n",
"7 1958 8 21412 1958.6219 314.93 316.19 314.00 315.29 314.93 \n",
"8 1958 9 21443 1958.7068 313.21 316.08 312.46 315.35 313.21 \n",
"9 1958 10 21473 1958.7890 NaN NaN 312.44 315.40 312.44 \n",
"10 1958 11 21504 1958.8740 313.33 315.20 313.62 315.46 313.33 \n",
"11 1958 12 21534 1958.9562 314.67 315.43 314.77 315.51 314.67 \n",
"12 1959 1 21565 1959.0411 315.58 315.54 315.62 315.57 315.58 \n",
"13 1959 2 21596 1959.1260 316.49 315.86 316.27 315.63 316.49 \n",
"14 1959 3 21624 1959.2027 316.65 315.38 316.98 315.69 316.65 \n",
"15 1959 4 21655 1959.2877 317.72 315.42 318.09 315.77 317.72 \n",
"16 1959 5 21685 1959.3699 318.29 315.49 318.65 315.85 318.29 \n",
"17 1959 6 21716 1959.4548 318.15 316.03 318.04 315.94 318.15 \n",
"18 1959 7 21746 1959.5370 316.54 315.86 316.67 316.03 316.54 \n",
"19 1959 8 21777 1959.6219 314.80 316.06 314.83 316.12 314.80 \n",
"20 1959 9 21808 1959.7068 313.84 316.72 313.32 316.22 313.84 \n",
"21 1959 10 21838 1959.7890 313.33 316.32 313.33 316.30 313.33 \n",
"22 1959 11 21869 1959.8740 314.81 316.68 314.54 316.39 314.81 \n",
"23 1959 12 21899 1959.9562 315.58 316.35 315.72 316.47 315.58 \n",
"24 1960 1 21930 1960.0410 316.43 316.39 316.61 316.56 316.43 \n",
"25 1960 2 21961 1960.1257 316.98 316.35 317.27 316.64 316.98 \n",
"26 1960 3 21990 1960.2049 317.58 316.28 318.03 316.71 317.58 \n",
"27 1960 4 22021 1960.2896 319.03 316.70 319.14 316.79 319.03 \n",
"28 1960 5 22051 1960.3716 320.04 317.22 319.67 316.86 320.04 \n",
"29 1960 6 22082 1960.4563 319.59 317.48 319.01 316.93 319.59 \n",
"30 1960 7 22112 1960.5383 318.18 317.53 317.60 316.98 318.18 \n",
"31 1960 8 22143 1960.6230 315.90 317.20 315.69 317.01 315.90 \n",
".. ... .. ... ... ... ... ... ... ... \n",
"715 2017 8 42962 2017.6219 405.24 406.77 405.20 406.76 405.24 \n",
"716 2017 9 42993 2017.7068 403.27 406.76 403.41 406.92 403.27 \n",
"717 2017 10 43023 2017.7890 403.64 407.25 403.47 407.07 403.64 \n",
"718 2017 11 43054 2017.8740 405.17 407.43 404.98 407.22 405.17 \n",
"719 2017 12 43084 2017.9562 406.75 407.68 406.46 407.37 406.75 \n",
"720 2018 1 43115 2018.0411 408.05 408.00 407.58 407.52 408.05 \n",
"721 2018 2 43146 2018.1260 408.34 407.59 408.44 407.68 408.34 \n",
"722 2018 3 43174 2018.2027 409.25 407.71 409.39 407.82 409.25 \n",
"723 2018 4 43205 2018.2877 410.30 407.51 410.82 408.00 410.30 \n",
"724 2018 5 43235 2018.3699 411.30 407.90 411.60 408.20 411.30 \n",
"725 2018 6 43266 2018.4548 410.88 408.30 410.97 408.42 410.88 \n",
"726 2018 7 43296 2018.5370 408.90 408.08 409.43 408.65 408.90 \n",
"727 2018 8 43327 2018.6219 407.10 408.63 407.33 408.91 407.10 \n",
"728 2018 9 43358 2018.7068 405.59 409.09 405.66 409.18 405.59 \n",
"729 2018 10 43388 2018.7890 405.99 409.62 405.83 409.44 405.99 \n",
"730 2018 11 43419 2018.8740 408.12 410.39 407.47 409.72 408.12 \n",
"731 2018 12 43449 2018.9562 409.23 410.16 409.07 409.97 409.23 \n",
"732 2019 1 43480 2019.0411 410.92 410.87 410.29 410.23 410.92 \n",
"733 2019 2 43511 2019.1260 411.66 410.90 411.24 410.47 411.66 \n",
"734 2019 3 43539 2019.2027 412.00 410.45 412.25 410.68 412.00 \n",
"735 2019 4 43570 2019.2877 413.52 410.72 413.73 410.91 413.52 \n",
"736 2019 5 43600 2019.3699 414.83 411.42 414.54 411.13 414.83 \n",
"737 2019 6 43631 2019.4548 413.96 411.38 413.91 411.35 413.96 \n",
"738 2019 7 43661 2019.5370 411.85 411.03 412.36 411.57 411.85 \n",
"739 2019 8 43692 2019.6219 410.08 411.62 410.23 411.81 410.08 \n",
"740 2019 9 43723 2019.7068 408.55 412.06 408.52 412.05 408.55 \n",
"741 2019 10 43753 2019.7890 408.43 412.07 408.67 412.29 408.43 \n",
"742 2019 11 43784 2019.8740 410.28 412.56 410.29 412.54 410.28 \n",
"743 2019 12 43814 2019.9562 411.85 412.78 411.88 412.79 411.85 \n",
"744 2020 1 43845 2020.0410 413.37 413.33 413.11 413.05 413.37 \n",
"\n",
" s6 \n",
"2 314.44 \n",
"3 315.16 \n",
"4 314.71 \n",
"5 315.14 \n",
"6 315.19 \n",
"7 316.19 \n",
"8 316.08 \n",
"9 315.40 \n",
"10 315.20 \n",
"11 315.43 \n",
"12 315.54 \n",
"13 315.86 \n",
"14 315.38 \n",
"15 315.42 \n",
"16 315.49 \n",
"17 316.03 \n",
"18 315.86 \n",
"19 316.06 \n",
"20 316.72 \n",
"21 316.32 \n",
"22 316.68 \n",
"23 316.35 \n",
"24 316.39 \n",
"25 316.35 \n",
"26 316.28 \n",
"27 316.70 \n",
"28 317.22 \n",
"29 317.48 \n",
"30 317.53 \n",
"31 317.20 \n",
".. ... \n",
"715 406.77 \n",
"716 406.76 \n",
"717 407.25 \n",
"718 407.43 \n",
"719 407.68 \n",
"720 408.00 \n",
"721 407.59 \n",
"722 407.71 \n",
"723 407.51 \n",
"724 407.90 \n",
"725 408.30 \n",
"726 408.08 \n",
"727 408.63 \n",
"728 409.09 \n",
"729 409.62 \n",
"730 410.39 \n",
"731 410.16 \n",
"732 410.87 \n",
"733 410.90 \n",
"734 410.45 \n",
"735 410.72 \n",
"736 411.42 \n",
"737 411.38 \n",
"738 411.03 \n",
"739 411.62 \n",
"740 412.06 \n",
"741 412.07 \n",
"742 412.56 \n",
"743 412.78 \n",
"744 413.33 \n",
"\n",
"[743 rows x 10 columns]"
"744 413.33 "
]
},
"execution_count": 33,
"execution_count": 71,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"#data = data.replace(-99.99,np.NaN);\n",
"data"
"date = datetime.date.today()\n",
"data = raw_data_1.copy()\n",
"data = data.drop(data[(data.Yr == date.year) & (data.Mn > date.month)].index)\n",
"data = data.replace(-99.99,np.NaN);\n",
"d = data[(data.Yr == date.year)].index\n",
"i,k = d[0], d[-1]\n",
"\n",
"while k>=i:\n",
" if data.loc[k].isnull().any():\n",
" data = data.drop(k)\n",
" else:\n",
" break\n",
" k = k - 1\n",
" \n",
"annee, mois = data.Yr[0], data.Mn[0]\n",
"d = data[(data.Yr == annee)].index\n",
"i,k = d[0], d[-1]\n",
"\n",
"while i<=k:\n",
" if data.loc[i].isnull().any():\n",
" data = data.drop(i)\n",
" else:\n",
" break\n",
" i = i + 1\n",
"InteractiveShell.ast_node_interactivity = \"all\" \n",
"data.head(3)\n",
"data.tail(3)\n",
"InteractiveShell.ast_node_interactivity = \"last_expr\" "
]
},
{
......@@ -1436,7 +566,7 @@
},
{
"cell_type": "code",
"execution_count": 128,
"execution_count": 72,
"metadata": {},
"outputs": [
{
......@@ -1474,32 +604,6 @@
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1958</td>\n",
" <td>1</td>\n",
" <td>21200</td>\n",
" <td>1958.0411</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1958</td>\n",
" <td>2</td>\n",
" <td>21231</td>\n",
" <td>1958.1260</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
......@@ -1564,176 +668,20 @@
" <td>321.83</td>\n",
" <td>319.45</td>\n",
" </tr>\n",
" <tr>\n",
" <th>745</th>\n",
" <td>2020</td>\n",
" <td>2</td>\n",
" <td>43876</td>\n",
" <td>2020.1257</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>746</th>\n",
" <td>2020</td>\n",
" <td>3</td>\n",
" <td>43905</td>\n",
" <td>2020.2049</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>747</th>\n",
" <td>2020</td>\n",
" <td>4</td>\n",
" <td>43936</td>\n",
" <td>2020.2896</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>748</th>\n",
" <td>2020</td>\n",
" <td>5</td>\n",
" <td>43966</td>\n",
" <td>2020.3716</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>749</th>\n",
" <td>2020</td>\n",
" <td>6</td>\n",
" <td>43997</td>\n",
" <td>2020.4563</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>750</th>\n",
" <td>2020</td>\n",
" <td>7</td>\n",
" <td>44027</td>\n",
" <td>2020.5383</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>751</th>\n",
" <td>2020</td>\n",
" <td>8</td>\n",
" <td>44058</td>\n",
" <td>2020.6230</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>752</th>\n",
" <td>2020</td>\n",
" <td>9</td>\n",
" <td>44089</td>\n",
" <td>2020.7077</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>753</th>\n",
" <td>2020</td>\n",
" <td>10</td>\n",
" <td>44119</td>\n",
" <td>2020.7896</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>754</th>\n",
" <td>2020</td>\n",
" <td>11</td>\n",
" <td>44150</td>\n",
" <td>2020.8743</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>755</th>\n",
" <td>2020</td>\n",
" <td>12</td>\n",
" <td>44180</td>\n",
" <td>2020.9563</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 s6\n",
"0 1958 1 21200 1958.0411 NaN NaN NaN NaN NaN NaN\n",
"1 1958 2 21231 1958.1260 NaN NaN NaN NaN NaN NaN\n",
"5 1958 6 21351 1958.4548 NaN NaN 317.24 315.14 317.24 315.14\n",
"9 1958 10 21473 1958.7890 NaN NaN 312.44 315.40 312.44 315.40\n",
"73 1964 2 23422 1964.1257 NaN NaN 320.01 319.36 320.01 319.36\n",
"74 1964 3 23451 1964.2049 NaN NaN 320.74 319.41 320.74 319.41\n",
"75 1964 4 23482 1964.2896 NaN NaN 321.83 319.45 321.83 319.45\n",
"745 2020 2 43876 2020.1257 NaN NaN NaN NaN NaN NaN\n",
"746 2020 3 43905 2020.2049 NaN NaN NaN NaN NaN NaN\n",
"747 2020 4 43936 2020.2896 NaN NaN NaN NaN NaN NaN\n",
"748 2020 5 43966 2020.3716 NaN NaN NaN NaN NaN NaN\n",
"749 2020 6 43997 2020.4563 NaN NaN NaN NaN NaN NaN\n",
"750 2020 7 44027 2020.5383 NaN NaN NaN NaN NaN NaN\n",
"751 2020 8 44058 2020.6230 NaN NaN NaN NaN NaN NaN\n",
"752 2020 9 44089 2020.7077 NaN NaN NaN NaN NaN NaN\n",
"753 2020 10 44119 2020.7896 NaN NaN NaN NaN NaN NaN\n",
"754 2020 11 44150 2020.8743 NaN NaN NaN NaN NaN NaN\n",
"755 2020 12 44180 2020.9563 NaN NaN NaN NaN NaN NaN"
"75 1964 4 23482 1964.2896 NaN NaN 321.83 319.45 321.83 319.45"
]
},
"execution_count": 128,
"execution_count": 72,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1744,7 +692,7 @@
},
{
"cell_type": "code",
"execution_count": 129,
"execution_count": 73,
"metadata": {},
"outputs": [
{
......@@ -1782,32 +730,6 @@
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1958</td>\n",
" <td>1</td>\n",
" <td>21200</td>\n",
" <td>1958.0411</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1958</td>\n",
" <td>2</td>\n",
" <td>21231</td>\n",
" <td>1958.1260</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>1958</td>\n",
" <td>3</td>\n",
......@@ -1846,20 +768,46 @@
" <td>317.51</td>\n",
" <td>314.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
" <td>1958</td>\n",
" <td>6</td>\n",
" <td>21351</td>\n",
" <td>1958.4548</td>\n",
" <td>NaN</td>\n",
" <td>NaN</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" <td>317.24</td>\n",
" <td>315.14</td>\n",
" </tr>\n",
" <tr>\n",
" <th>6</th>\n",
" <td>1958</td>\n",
" <td>7</td>\n",
" <td>21381</td>\n",
" <td>1958.5370</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" <td>315.86</td>\n",
" <td>315.22</td>\n",
" <td>315.86</td>\n",
" <td>315.19</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 s5 s6\n",
"0 1958 1 21200 1958.0411 NaN NaN NaN NaN NaN NaN\n",
"1 1958 2 21231 1958.1260 NaN NaN NaN NaN NaN NaN\n",
"2 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 315.70 314.44\n",
"3 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 317.45 315.16\n",
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 314.71"
"4 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 317.51 314.71\n",
"5 1958 6 21351 1958.4548 NaN NaN 317.24 315.14 317.24 315.14\n",
"6 1958 7 21381 1958.5370 315.86 315.19 315.86 315.22 315.86 315.19"
]
},
"execution_count": 129,
"execution_count": 73,
"metadata": {},
"output_type": "execute_result"
}
......@@ -1878,7 +826,7 @@
},
{
"cell_type": "code",
"execution_count": 130,
"execution_count": 74,
"metadata": {
"scrolled": true
},
......@@ -1931,7 +879,7 @@
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1958-03-01</th>\n",
" <th>1958-05-01</th>\n",
" <td>1958</td>\n",
" <td>3</td>\n",
" <td>21259</td>\n",
......@@ -1944,7 +892,7 @@
" <td>314.44</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-04-01</th>\n",
" <th>1958-06-01</th>\n",
" <td>1958</td>\n",
" <td>4</td>\n",
" <td>21290</td>\n",
......@@ -1957,7 +905,7 @@
" <td>315.16</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-05-01</th>\n",
" <th>1958-07-01</th>\n",
" <td>1958</td>\n",
" <td>5</td>\n",
" <td>21320</td>\n",
......@@ -1970,7 +918,7 @@
" <td>314.71</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-07-01</th>\n",
" <th>1958-09-01</th>\n",
" <td>1958</td>\n",
" <td>7</td>\n",
" <td>21381</td>\n",
......@@ -1983,7 +931,7 @@
" <td>315.19</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1958-08-01</th>\n",
" <th>1958-10-01</th>\n",
" <td>1958</td>\n",
" <td>8</td>\n",
" <td>21412</td>\n",
......@@ -2002,22 +950,22 @@
"text/plain": [
" Yr Mn Date 1 Date 2 s1 s2 s3 s4 \\\n",
"period \n",
"1958-03-01 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 \n",
"1958-04-01 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 \n",
"1958-05-01 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 \n",
"1958-07-01 1958 7 21381 1958.5370 315.86 315.19 315.86 315.22 \n",
"1958-08-01 1958 8 21412 1958.6219 314.93 316.19 314.00 315.29 \n",
"1958-05-01 1958 3 21259 1958.2027 315.70 314.44 316.19 314.91 \n",
"1958-06-01 1958 4 21290 1958.2877 317.45 315.16 317.30 314.99 \n",
"1958-07-01 1958 5 21320 1958.3699 317.51 314.71 317.86 315.06 \n",
"1958-09-01 1958 7 21381 1958.5370 315.86 315.19 315.86 315.22 \n",
"1958-10-01 1958 8 21412 1958.6219 314.93 316.19 314.00 315.29 \n",
"\n",
" s5 s6 \n",
"period \n",
"1958-03-01 315.70 314.44 \n",
"1958-04-01 317.45 315.16 \n",
"1958-05-01 317.51 314.71 \n",
"1958-07-01 315.86 315.19 \n",
"1958-08-01 314.93 316.19 "
"1958-05-01 315.70 314.44 \n",
"1958-06-01 317.45 315.16 \n",
"1958-07-01 317.51 314.71 \n",
"1958-09-01 315.86 315.19 \n",
"1958-10-01 314.93 316.19 "
]
},
"execution_count": 130,
"execution_count": 74,
"metadata": {},
"output_type": "execute_result"
}
......@@ -2168,6 +1116,15 @@
"ax2.plot(data_cube['co2']); "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous recherchons l'oscillation lente, pour cela nous allons opéré en 2 etapes :\n",
"- Recherche de la fréquence lente, par FFT. Pour cela nous devons faire une interpolation pour les quelque points manquant. Ce qui nous donnera une gamme de fréquence possible (échantillonne de la FFT)\n",
"- Pour la gamme de fréquence spécifiées, nous estimons tous les paramètres pour une approximation sinusoïdale (moyenne, amplitude et phase) par les moindres carrés ordinaires \n"
]
},
{
"cell_type": "code",
"execution_count": null,
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment