Nebraska_past_years_phenology.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>2002_025_NDVI</th>\n",
       "      <th>2002_033_NDVI</th>\n",
       "      <th>2002_041_NDVI</th>\n",
       "      <th>2002_049_NDVI</th>\n",
       "      <th>2002_057_NDVI</th>\n",
       "      <th>2002_065_NDVI</th>\n",
       "      <th>2002_073_NDVI</th>\n",
       "      <th>2002_081_NDVI</th>\n",
       "      <th>2002_089_NDVI</th>\n",
       "      <th>2002_097_NDVI</th>\n",
       "      <th>...</th>\n",
       "      <th>2002_297_NDVI</th>\n",
       "      <th>2002_305_NDVI</th>\n",
       "      <th>2002_313_NDVI</th>\n",
       "      <th>2002_321_NDVI</th>\n",
       "      <th>2002_329_NDVI</th>\n",
       "      <th>2002_337_NDVI</th>\n",
       "      <th>2002_345_NDVI</th>\n",
       "      <th>2002_353_NDVI</th>\n",
       "      <th>2002_361_NDVI</th>\n",
       "      <th>GEOID</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>199.979369</td>\n",
       "      <td>195.341170</td>\n",
       "      <td>200.847186</td>\n",
       "      <td>202.875689</td>\n",
       "      <td>NaN</td>\n",
       "      <td>194.099093</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>208.873180</td>\n",
       "      <td>221.068704</td>\n",
       "      <td>...</td>\n",
       "      <td>281.463133</td>\n",
       "      <td>NaN</td>\n",
       "      <td>249.743801</td>\n",
       "      <td>237.984343</td>\n",
       "      <td>229.278986</td>\n",
       "      <td>229.250037</td>\n",
       "      <td>219.821214</td>\n",
       "      <td>214.984558</td>\n",
       "      <td>212.731740</td>\n",
       "      <td>31007</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>203.351533</td>\n",
       "      <td>201.146636</td>\n",
       "      <td>203.924538</td>\n",
       "      <td>216.894404</td>\n",
       "      <td>NaN</td>\n",
       "      <td>203.419017</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>191.155802</td>\n",
       "      <td>204.691938</td>\n",
       "      <td>...</td>\n",
       "      <td>301.913183</td>\n",
       "      <td>266.251160</td>\n",
       "      <td>239.458615</td>\n",
       "      <td>250.944560</td>\n",
       "      <td>231.718090</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>214.261857</td>\n",
       "      <td>219.209543</td>\n",
       "      <td>31045</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\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",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>...</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",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>31075</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>212.819195</td>\n",
       "      <td>213.835711</td>\n",
       "      <td>224.030314</td>\n",
       "      <td>216.884286</td>\n",
       "      <td>208.877590</td>\n",
       "      <td>210.909783</td>\n",
       "      <td>213.257699</td>\n",
       "      <td>215.908769</td>\n",
       "      <td>215.358000</td>\n",
       "      <td>221.722039</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>253.638928</td>\n",
       "      <td>238.831443</td>\n",
       "      <td>234.096765</td>\n",
       "      <td>226.862497</td>\n",
       "      <td>223.787780</td>\n",
       "      <td>222.051742</td>\n",
       "      <td>227.336760</td>\n",
       "      <td>221.165678</td>\n",
       "      <td>31113</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>221.803229</td>\n",
       "      <td>NaN</td>\n",
       "      <td>208.278847</td>\n",
       "      <td>234.653056</td>\n",
       "      <td>220.686890</td>\n",
       "      <td>232.823441</td>\n",
       "      <td>223.750504</td>\n",
       "      <td>231.863896</td>\n",
       "      <td>232.966169</td>\n",
       "      <td>244.505615</td>\n",
       "      <td>...</td>\n",
       "      <td>259.154760</td>\n",
       "      <td>220.426673</td>\n",
       "      <td>226.173553</td>\n",
       "      <td>235.380445</td>\n",
       "      <td>228.814447</td>\n",
       "      <td>223.299799</td>\n",
       "      <td>231.816201</td>\n",
       "      <td>232.165852</td>\n",
       "      <td>228.895232</td>\n",
       "      <td>31185</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>229.916099</td>\n",
       "      <td>226.190997</td>\n",
       "      <td>227.608255</td>\n",
       "      <td>222.942409</td>\n",
       "      <td>212.259014</td>\n",
       "      <td>207.636691</td>\n",
       "      <td>NaN</td>\n",
       "      <td>206.848260</td>\n",
       "      <td>214.297155</td>\n",
       "      <td>227.841844</td>\n",
       "      <td>...</td>\n",
       "      <td>296.791499</td>\n",
       "      <td>342.314950</td>\n",
       "      <td>304.395678</td>\n",
       "      <td>285.335320</td>\n",
       "      <td>253.459906</td>\n",
       "      <td>242.168471</td>\n",
       "      <td>236.214690</td>\n",
       "      <td>230.999540</td>\n",
       "      <td>234.907088</td>\n",
       "      <td>31005</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>235.391085</td>\n",
       "      <td>214.114682</td>\n",
       "      <td>226.964055</td>\n",
       "      <td>216.516103</td>\n",
       "      <td>219.143350</td>\n",
       "      <td>NaN</td>\n",
       "      <td>215.304442</td>\n",
       "      <td>213.190887</td>\n",
       "      <td>222.066034</td>\n",
       "      <td>239.968478</td>\n",
       "      <td>...</td>\n",
       "      <td>298.770726</td>\n",
       "      <td>274.771954</td>\n",
       "      <td>255.840324</td>\n",
       "      <td>246.847215</td>\n",
       "      <td>249.119410</td>\n",
       "      <td>235.243378</td>\n",
       "      <td>237.865975</td>\n",
       "      <td>233.151520</td>\n",
       "      <td>226.999504</td>\n",
       "      <td>31071</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>212.339754</td>\n",
       "      <td>212.275480</td>\n",
       "      <td>213.647377</td>\n",
       "      <td>210.994667</td>\n",
       "      <td>NaN</td>\n",
       "      <td>207.719078</td>\n",
       "      <td>217.716226</td>\n",
       "      <td>220.356139</td>\n",
       "      <td>212.470068</td>\n",
       "      <td>230.302013</td>\n",
       "      <td>...</td>\n",
       "      <td>302.585162</td>\n",
       "      <td>261.858709</td>\n",
       "      <td>251.620920</td>\n",
       "      <td>240.537196</td>\n",
       "      <td>226.203985</td>\n",
       "      <td>227.249151</td>\n",
       "      <td>225.747537</td>\n",
       "      <td>217.432238</td>\n",
       "      <td>217.896430</td>\n",
       "      <td>31101</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>226.551550</td>\n",
       "      <td>NaN</td>\n",
       "      <td>225.780873</td>\n",
       "      <td>237.574696</td>\n",
       "      <td>227.863166</td>\n",
       "      <td>221.063138</td>\n",
       "      <td>236.405947</td>\n",
       "      <td>239.055175</td>\n",
       "      <td>236.669350</td>\n",
       "      <td>257.192697</td>\n",
       "      <td>...</td>\n",
       "      <td>269.260749</td>\n",
       "      <td>248.306612</td>\n",
       "      <td>239.132158</td>\n",
       "      <td>243.130186</td>\n",
       "      <td>240.258149</td>\n",
       "      <td>235.511602</td>\n",
       "      <td>238.513087</td>\n",
       "      <td>239.084113</td>\n",
       "      <td>236.272009</td>\n",
       "      <td>31151</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>222.940246</td>\n",
       "      <td>NaN</td>\n",
       "      <td>216.315210</td>\n",
       "      <td>238.430340</td>\n",
       "      <td>224.973504</td>\n",
       "      <td>236.699749</td>\n",
       "      <td>229.727361</td>\n",
       "      <td>238.107216</td>\n",
       "      <td>238.114405</td>\n",
       "      <td>250.240022</td>\n",
       "      <td>...</td>\n",
       "      <td>243.716207</td>\n",
       "      <td>226.287687</td>\n",
       "      <td>217.557427</td>\n",
       "      <td>227.190027</td>\n",
       "      <td>225.817783</td>\n",
       "      <td>223.190388</td>\n",
       "      <td>228.394813</td>\n",
       "      <td>229.302264</td>\n",
       "      <td>226.193111</td>\n",
       "      <td>31059</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>213.087390</td>\n",
       "      <td>200.290834</td>\n",
       "      <td>197.167042</td>\n",
       "      <td>211.494392</td>\n",
       "      <td>203.471166</td>\n",
       "      <td>NaN</td>\n",
       "      <td>200.169625</td>\n",
       "      <td>NaN</td>\n",
       "      <td>204.162609</td>\n",
       "      <td>218.450447</td>\n",
       "      <td>...</td>\n",
       "      <td>276.328730</td>\n",
       "      <td>241.475141</td>\n",
       "      <td>241.116998</td>\n",
       "      <td>234.251233</td>\n",
       "      <td>232.884580</td>\n",
       "      <td>228.291128</td>\n",
       "      <td>224.764699</td>\n",
       "      <td>220.751375</td>\n",
       "      <td>220.066632</td>\n",
       "      <td>31119</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>215.664492</td>\n",
       "      <td>208.997795</td>\n",
       "      <td>209.394583</td>\n",
       "      <td>217.713067</td>\n",
       "      <td>195.850396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>216.781553</td>\n",
       "      <td>NaN</td>\n",
       "      <td>212.394788</td>\n",
       "      <td>217.624835</td>\n",
       "      <td>...</td>\n",
       "      <td>296.039235</td>\n",
       "      <td>267.350001</td>\n",
       "      <td>251.110718</td>\n",
       "      <td>246.170961</td>\n",
       "      <td>238.863715</td>\n",
       "      <td>228.058503</td>\n",
       "      <td>233.839241</td>\n",
       "      <td>231.545832</td>\n",
       "      <td>229.893855</td>\n",
       "      <td>31183</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>229.262530</td>\n",
       "      <td>230.118864</td>\n",
       "      <td>238.688525</td>\n",
       "      <td>227.655732</td>\n",
       "      <td>245.322250</td>\n",
       "      <td>NaN</td>\n",
       "      <td>217.796205</td>\n",
       "      <td>221.179550</td>\n",
       "      <td>227.242232</td>\n",
       "      <td>247.042382</td>\n",
       "      <td>...</td>\n",
       "      <td>311.530339</td>\n",
       "      <td>283.933175</td>\n",
       "      <td>275.990437</td>\n",
       "      <td>268.327353</td>\n",
       "      <td>262.677335</td>\n",
       "      <td>451.712677</td>\n",
       "      <td>260.192161</td>\n",
       "      <td>250.888916</td>\n",
       "      <td>232.393577</td>\n",
       "      <td>31009</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>228.198784</td>\n",
       "      <td>217.080894</td>\n",
       "      <td>227.626282</td>\n",
       "      <td>231.060463</td>\n",
       "      <td>NaN</td>\n",
       "      <td>168.647194</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>220.692757</td>\n",
       "      <td>251.722046</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>351.530609</td>\n",
       "      <td>479.382004</td>\n",
       "      <td>367.240463</td>\n",
       "      <td>391.930695</td>\n",
       "      <td>NaN</td>\n",
       "      <td>353.646042</td>\n",
       "      <td>425.709457</td>\n",
       "      <td>398.537155</td>\n",
       "      <td>31091</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>NaN</td>\n",
       "      <td>220.111292</td>\n",
       "      <td>228.955993</td>\n",
       "      <td>218.872942</td>\n",
       "      <td>221.335526</td>\n",
       "      <td>NaN</td>\n",
       "      <td>216.722434</td>\n",
       "      <td>216.859332</td>\n",
       "      <td>220.246166</td>\n",
       "      <td>236.874035</td>\n",
       "      <td>...</td>\n",
       "      <td>305.912349</td>\n",
       "      <td>279.553929</td>\n",
       "      <td>268.892267</td>\n",
       "      <td>253.810679</td>\n",
       "      <td>250.104062</td>\n",
       "      <td>242.322301</td>\n",
       "      <td>238.880964</td>\n",
       "      <td>236.060772</td>\n",
       "      <td>228.101122</td>\n",
       "      <td>31115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>223.113451</td>\n",
       "      <td>216.537375</td>\n",
       "      <td>205.879802</td>\n",
       "      <td>220.499701</td>\n",
       "      <td>211.254640</td>\n",
       "      <td>NaN</td>\n",
       "      <td>227.043816</td>\n",
       "      <td>197.070824</td>\n",
       "      <td>214.770376</td>\n",
       "      <td>225.288235</td>\n",
       "      <td>...</td>\n",
       "      <td>293.473458</td>\n",
       "      <td>260.204741</td>\n",
       "      <td>253.929321</td>\n",
       "      <td>244.339834</td>\n",
       "      <td>240.340041</td>\n",
       "      <td>234.601654</td>\n",
       "      <td>232.839708</td>\n",
       "      <td>227.979829</td>\n",
       "      <td>225.623739</td>\n",
       "      <td>31003</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>204.493015</td>\n",
       "      <td>205.783988</td>\n",
       "      <td>214.464836</td>\n",
       "      <td>208.363281</td>\n",
       "      <td>204.308503</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>208.424568</td>\n",
       "      <td>201.091952</td>\n",
       "      <td>225.252987</td>\n",
       "      <td>...</td>\n",
       "      <td>288.203969</td>\n",
       "      <td>262.108939</td>\n",
       "      <td>255.520184</td>\n",
       "      <td>245.465956</td>\n",
       "      <td>241.372719</td>\n",
       "      <td>236.413167</td>\n",
       "      <td>233.739821</td>\n",
       "      <td>232.305500</td>\n",
       "      <td>224.970702</td>\n",
       "      <td>31027</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>223.015067</td>\n",
       "      <td>216.410609</td>\n",
       "      <td>222.997272</td>\n",
       "      <td>215.303024</td>\n",
       "      <td>199.178759</td>\n",
       "      <td>222.789958</td>\n",
       "      <td>225.640728</td>\n",
       "      <td>224.922155</td>\n",
       "      <td>224.525629</td>\n",
       "      <td>237.756216</td>\n",
       "      <td>...</td>\n",
       "      <td>293.017282</td>\n",
       "      <td>275.189193</td>\n",
       "      <td>259.419199</td>\n",
       "      <td>245.260745</td>\n",
       "      <td>239.860045</td>\n",
       "      <td>238.887649</td>\n",
       "      <td>234.719141</td>\n",
       "      <td>230.851799</td>\n",
       "      <td>227.575202</td>\n",
       "      <td>31029</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>211.463358</td>\n",
       "      <td>NaN</td>\n",
       "      <td>219.454761</td>\n",
       "      <td>231.509915</td>\n",
       "      <td>218.950347</td>\n",
       "      <td>230.677036</td>\n",
       "      <td>223.316829</td>\n",
       "      <td>227.962304</td>\n",
       "      <td>233.859957</td>\n",
       "      <td>248.770612</td>\n",
       "      <td>...</td>\n",
       "      <td>255.214874</td>\n",
       "      <td>236.707776</td>\n",
       "      <td>224.221904</td>\n",
       "      <td>228.630716</td>\n",
       "      <td>225.810371</td>\n",
       "      <td>222.145856</td>\n",
       "      <td>228.175999</td>\n",
       "      <td>227.118750</td>\n",
       "      <td>225.065908</td>\n",
       "      <td>31035</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>222.050867</td>\n",
       "      <td>209.426620</td>\n",
       "      <td>212.194563</td>\n",
       "      <td>221.793128</td>\n",
       "      <td>221.055657</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>212.024915</td>\n",
       "      <td>215.886264</td>\n",
       "      <td>242.860799</td>\n",
       "      <td>...</td>\n",
       "      <td>272.912082</td>\n",
       "      <td>246.144415</td>\n",
       "      <td>248.588913</td>\n",
       "      <td>243.150218</td>\n",
       "      <td>243.498574</td>\n",
       "      <td>238.942913</td>\n",
       "      <td>234.542641</td>\n",
       "      <td>234.752078</td>\n",
       "      <td>227.161276</td>\n",
       "      <td>31039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>197.230012</td>\n",
       "      <td>196.013566</td>\n",
       "      <td>198.147255</td>\n",
       "      <td>196.737671</td>\n",
       "      <td>NaN</td>\n",
       "      <td>202.529293</td>\n",
       "      <td>200.720407</td>\n",
       "      <td>NaN</td>\n",
       "      <td>209.780270</td>\n",
       "      <td>218.362723</td>\n",
       "      <td>...</td>\n",
       "      <td>297.675962</td>\n",
       "      <td>242.170773</td>\n",
       "      <td>232.205259</td>\n",
       "      <td>224.833441</td>\n",
       "      <td>213.230526</td>\n",
       "      <td>213.055090</td>\n",
       "      <td>211.939743</td>\n",
       "      <td>205.739942</td>\n",
       "      <td>208.560988</td>\n",
       "      <td>31049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>207.651102</td>\n",
       "      <td>201.204848</td>\n",
       "      <td>185.729315</td>\n",
       "      <td>207.247734</td>\n",
       "      <td>213.281543</td>\n",
       "      <td>NaN</td>\n",
       "      <td>176.103800</td>\n",
       "      <td>209.581853</td>\n",
       "      <td>200.597896</td>\n",
       "      <td>227.037783</td>\n",
       "      <td>...</td>\n",
       "      <td>290.014531</td>\n",
       "      <td>263.704238</td>\n",
       "      <td>257.658566</td>\n",
       "      <td>243.814411</td>\n",
       "      <td>239.592426</td>\n",
       "      <td>239.531701</td>\n",
       "      <td>231.712770</td>\n",
       "      <td>229.048020</td>\n",
       "      <td>222.688131</td>\n",
       "      <td>31051</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>242.666308</td>\n",
       "      <td>231.599703</td>\n",
       "      <td>235.794515</td>\n",
       "      <td>225.961006</td>\n",
       "      <td>234.321865</td>\n",
       "      <td>230.243851</td>\n",
       "      <td>232.960386</td>\n",
       "      <td>228.276336</td>\n",
       "      <td>236.611572</td>\n",
       "      <td>249.283510</td>\n",
       "      <td>...</td>\n",
       "      <td>388.105254</td>\n",
       "      <td>295.129787</td>\n",
       "      <td>284.616600</td>\n",
       "      <td>264.080364</td>\n",
       "      <td>259.291434</td>\n",
       "      <td>255.999083</td>\n",
       "      <td>244.524268</td>\n",
       "      <td>245.147904</td>\n",
       "      <td>237.348237</td>\n",
       "      <td>31057</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>219.450249</td>\n",
       "      <td>212.416231</td>\n",
       "      <td>226.055764</td>\n",
       "      <td>219.538792</td>\n",
       "      <td>216.250584</td>\n",
       "      <td>223.426848</td>\n",
       "      <td>226.041834</td>\n",
       "      <td>223.943916</td>\n",
       "      <td>225.558735</td>\n",
       "      <td>245.639428</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>238.529693</td>\n",
       "      <td>229.845925</td>\n",
       "      <td>222.964434</td>\n",
       "      <td>220.352473</td>\n",
       "      <td>219.302174</td>\n",
       "      <td>221.513899</td>\n",
       "      <td>218.431285</td>\n",
       "      <td>214.112700</td>\n",
       "      <td>31063</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>232.153714</td>\n",
       "      <td>NaN</td>\n",
       "      <td>231.175874</td>\n",
       "      <td>237.683450</td>\n",
       "      <td>223.718070</td>\n",
       "      <td>224.490121</td>\n",
       "      <td>235.426831</td>\n",
       "      <td>240.613085</td>\n",
       "      <td>236.556450</td>\n",
       "      <td>268.807481</td>\n",
       "      <td>...</td>\n",
       "      <td>284.272495</td>\n",
       "      <td>260.207541</td>\n",
       "      <td>254.054105</td>\n",
       "      <td>254.585254</td>\n",
       "      <td>251.166482</td>\n",
       "      <td>243.519961</td>\n",
       "      <td>248.215218</td>\n",
       "      <td>248.104461</td>\n",
       "      <td>243.690696</td>\n",
       "      <td>31067</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>211.420660</td>\n",
       "      <td>204.102347</td>\n",
       "      <td>219.626292</td>\n",
       "      <td>214.845444</td>\n",
       "      <td>205.364904</td>\n",
       "      <td>220.773557</td>\n",
       "      <td>221.179931</td>\n",
       "      <td>189.675670</td>\n",
       "      <td>211.559741</td>\n",
       "      <td>219.927054</td>\n",
       "      <td>...</td>\n",
       "      <td>293.056082</td>\n",
       "      <td>261.224720</td>\n",
       "      <td>252.441535</td>\n",
       "      <td>245.945802</td>\n",
       "      <td>240.673315</td>\n",
       "      <td>228.889383</td>\n",
       "      <td>232.114628</td>\n",
       "      <td>232.490905</td>\n",
       "      <td>227.547505</td>\n",
       "      <td>31077</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>230.244046</td>\n",
       "      <td>NaN</td>\n",
       "      <td>222.098803</td>\n",
       "      <td>224.550970</td>\n",
       "      <td>215.997692</td>\n",
       "      <td>224.721818</td>\n",
       "      <td>224.347877</td>\n",
       "      <td>218.873843</td>\n",
       "      <td>223.009762</td>\n",
       "      <td>239.527747</td>\n",
       "      <td>...</td>\n",
       "      <td>298.416446</td>\n",
       "      <td>263.235406</td>\n",
       "      <td>260.235061</td>\n",
       "      <td>251.679767</td>\n",
       "      <td>239.487084</td>\n",
       "      <td>234.138533</td>\n",
       "      <td>240.196273</td>\n",
       "      <td>236.796861</td>\n",
       "      <td>231.038189</td>\n",
       "      <td>31093</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>228.113689</td>\n",
       "      <td>224.909362</td>\n",
       "      <td>225.451196</td>\n",
       "      <td>228.623511</td>\n",
       "      <td>230.275897</td>\n",
       "      <td>217.735140</td>\n",
       "      <td>188.953720</td>\n",
       "      <td>NaN</td>\n",
       "      <td>230.458329</td>\n",
       "      <td>246.522200</td>\n",
       "      <td>...</td>\n",
       "      <td>281.771326</td>\n",
       "      <td>264.266931</td>\n",
       "      <td>248.872840</td>\n",
       "      <td>242.603717</td>\n",
       "      <td>229.880533</td>\n",
       "      <td>233.534608</td>\n",
       "      <td>226.502091</td>\n",
       "      <td>223.341168</td>\n",
       "      <td>220.525165</td>\n",
       "      <td>31123</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>229.464527</td>\n",
       "      <td>NaN</td>\n",
       "      <td>231.358043</td>\n",
       "      <td>240.091513</td>\n",
       "      <td>226.052552</td>\n",
       "      <td>213.447005</td>\n",
       "      <td>235.330268</td>\n",
       "      <td>255.875778</td>\n",
       "      <td>241.406905</td>\n",
       "      <td>270.577516</td>\n",
       "      <td>...</td>\n",
       "      <td>304.303369</td>\n",
       "      <td>288.503341</td>\n",
       "      <td>277.611749</td>\n",
       "      <td>274.740986</td>\n",
       "      <td>269.581529</td>\n",
       "      <td>266.464531</td>\n",
       "      <td>261.839946</td>\n",
       "      <td>265.305399</td>\n",
       "      <td>263.990195</td>\n",
       "      <td>31133</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>212.331221</td>\n",
       "      <td>208.906524</td>\n",
       "      <td>204.555088</td>\n",
       "      <td>214.181301</td>\n",
       "      <td>209.629666</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>207.101442</td>\n",
       "      <td>208.294738</td>\n",
       "      <td>226.756304</td>\n",
       "      <td>...</td>\n",
       "      <td>297.018123</td>\n",
       "      <td>262.781072</td>\n",
       "      <td>261.065115</td>\n",
       "      <td>246.609795</td>\n",
       "      <td>241.701473</td>\n",
       "      <td>236.014225</td>\n",
       "      <td>232.692441</td>\n",
       "      <td>228.584283</td>\n",
       "      <td>225.481419</td>\n",
       "      <td>31139</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",
       "      <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",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>193.814649</td>\n",
       "      <td>206.311882</td>\n",
       "      <td>215.171245</td>\n",
       "      <td>203.209433</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>219.551935</td>\n",
       "      <td>233.150780</td>\n",
       "      <td>...</td>\n",
       "      <td>291.935394</td>\n",
       "      <td>291.021334</td>\n",
       "      <td>280.745745</td>\n",
       "      <td>274.320888</td>\n",
       "      <td>259.884211</td>\n",
       "      <td>252.452640</td>\n",
       "      <td>253.940006</td>\n",
       "      <td>261.101928</td>\n",
       "      <td>238.774878</td>\n",
       "      <td>31103</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>64</th>\n",
       "      <td>214.917312</td>\n",
       "      <td>205.954823</td>\n",
       "      <td>222.972146</td>\n",
       "      <td>214.983491</td>\n",
       "      <td>207.958192</td>\n",
       "      <td>210.273061</td>\n",
       "      <td>217.328152</td>\n",
       "      <td>213.658577</td>\n",
       "      <td>218.705533</td>\n",
       "      <td>231.241189</td>\n",
       "      <td>...</td>\n",
       "      <td>289.013780</td>\n",
       "      <td>264.161719</td>\n",
       "      <td>246.881307</td>\n",
       "      <td>237.388786</td>\n",
       "      <td>230.002543</td>\n",
       "      <td>228.818851</td>\n",
       "      <td>227.221174</td>\n",
       "      <td>221.606303</td>\n",
       "      <td>224.714297</td>\n",
       "      <td>31111</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>220.803198</td>\n",
       "      <td>NaN</td>\n",
       "      <td>219.848205</td>\n",
       "      <td>232.277201</td>\n",
       "      <td>222.299859</td>\n",
       "      <td>230.929156</td>\n",
       "      <td>NaN</td>\n",
       "      <td>227.791331</td>\n",
       "      <td>233.589550</td>\n",
       "      <td>245.372966</td>\n",
       "      <td>...</td>\n",
       "      <td>282.279564</td>\n",
       "      <td>247.758407</td>\n",
       "      <td>248.450048</td>\n",
       "      <td>244.555951</td>\n",
       "      <td>234.962705</td>\n",
       "      <td>226.276053</td>\n",
       "      <td>237.354817</td>\n",
       "      <td>234.900983</td>\n",
       "      <td>229.521826</td>\n",
       "      <td>31121</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>66</th>\n",
       "      <td>242.433150</td>\n",
       "      <td>NaN</td>\n",
       "      <td>227.342188</td>\n",
       "      <td>231.498883</td>\n",
       "      <td>222.243880</td>\n",
       "      <td>221.989210</td>\n",
       "      <td>225.755707</td>\n",
       "      <td>231.777910</td>\n",
       "      <td>232.180752</td>\n",
       "      <td>250.745227</td>\n",
       "      <td>...</td>\n",
       "      <td>269.413014</td>\n",
       "      <td>252.192130</td>\n",
       "      <td>242.588541</td>\n",
       "      <td>238.737482</td>\n",
       "      <td>236.149806</td>\n",
       "      <td>237.692255</td>\n",
       "      <td>236.785699</td>\n",
       "      <td>235.181311</td>\n",
       "      <td>233.594761</td>\n",
       "      <td>31129</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>67</th>\n",
       "      <td>221.300050</td>\n",
       "      <td>NaN</td>\n",
       "      <td>203.774410</td>\n",
       "      <td>233.625819</td>\n",
       "      <td>223.150836</td>\n",
       "      <td>231.907039</td>\n",
       "      <td>NaN</td>\n",
       "      <td>228.982839</td>\n",
       "      <td>232.746413</td>\n",
       "      <td>246.538325</td>\n",
       "      <td>...</td>\n",
       "      <td>264.688497</td>\n",
       "      <td>232.829392</td>\n",
       "      <td>235.869030</td>\n",
       "      <td>238.202210</td>\n",
       "      <td>233.001977</td>\n",
       "      <td>224.561481</td>\n",
       "      <td>236.080358</td>\n",
       "      <td>234.148641</td>\n",
       "      <td>231.898808</td>\n",
       "      <td>31143</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>68</th>\n",
       "      <td>219.407627</td>\n",
       "      <td>203.312203</td>\n",
       "      <td>220.341442</td>\n",
       "      <td>220.498925</td>\n",
       "      <td>221.442048</td>\n",
       "      <td>212.509602</td>\n",
       "      <td>221.038805</td>\n",
       "      <td>216.600120</td>\n",
       "      <td>221.534680</td>\n",
       "      <td>244.483284</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>243.781136</td>\n",
       "      <td>234.991093</td>\n",
       "      <td>228.454418</td>\n",
       "      <td>224.015630</td>\n",
       "      <td>223.315872</td>\n",
       "      <td>223.337294</td>\n",
       "      <td>220.213319</td>\n",
       "      <td>217.213538</td>\n",
       "      <td>31145</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>212.861784</td>\n",
       "      <td>NaN</td>\n",
       "      <td>217.958204</td>\n",
       "      <td>233.059434</td>\n",
       "      <td>217.592601</td>\n",
       "      <td>222.831417</td>\n",
       "      <td>NaN</td>\n",
       "      <td>227.692886</td>\n",
       "      <td>223.013312</td>\n",
       "      <td>247.533964</td>\n",
       "      <td>...</td>\n",
       "      <td>275.360540</td>\n",
       "      <td>243.966690</td>\n",
       "      <td>251.070115</td>\n",
       "      <td>247.055581</td>\n",
       "      <td>242.484065</td>\n",
       "      <td>230.180067</td>\n",
       "      <td>238.788935</td>\n",
       "      <td>238.341592</td>\n",
       "      <td>231.487470</td>\n",
       "      <td>31155</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70</th>\n",
       "      <td>223.332273</td>\n",
       "      <td>NaN</td>\n",
       "      <td>203.036219</td>\n",
       "      <td>235.870952</td>\n",
       "      <td>223.276035</td>\n",
       "      <td>222.781461</td>\n",
       "      <td>231.215785</td>\n",
       "      <td>234.651766</td>\n",
       "      <td>234.085998</td>\n",
       "      <td>254.890991</td>\n",
       "      <td>...</td>\n",
       "      <td>268.962974</td>\n",
       "      <td>237.391334</td>\n",
       "      <td>242.550593</td>\n",
       "      <td>243.988495</td>\n",
       "      <td>237.599357</td>\n",
       "      <td>229.323270</td>\n",
       "      <td>236.509636</td>\n",
       "      <td>237.002648</td>\n",
       "      <td>233.973439</td>\n",
       "      <td>31159</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>197.877592</td>\n",
       "      <td>203.440447</td>\n",
       "      <td>200.271968</td>\n",
       "      <td>203.080493</td>\n",
       "      <td>218.074661</td>\n",
       "      <td>204.468616</td>\n",
       "      <td>NaN</td>\n",
       "      <td>181.880683</td>\n",
       "      <td>197.219207</td>\n",
       "      <td>211.149236</td>\n",
       "      <td>...</td>\n",
       "      <td>279.002013</td>\n",
       "      <td>232.825632</td>\n",
       "      <td>243.402558</td>\n",
       "      <td>224.848668</td>\n",
       "      <td>214.510684</td>\n",
       "      <td>208.763332</td>\n",
       "      <td>206.483264</td>\n",
       "      <td>204.298354</td>\n",
       "      <td>210.421736</td>\n",
       "      <td>31165</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72</th>\n",
       "      <td>206.045084</td>\n",
       "      <td>198.668973</td>\n",
       "      <td>188.038388</td>\n",
       "      <td>210.902233</td>\n",
       "      <td>209.045800</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>203.862209</td>\n",
       "      <td>204.054247</td>\n",
       "      <td>228.091318</td>\n",
       "      <td>...</td>\n",
       "      <td>272.902348</td>\n",
       "      <td>241.825806</td>\n",
       "      <td>241.467507</td>\n",
       "      <td>234.108402</td>\n",
       "      <td>234.985950</td>\n",
       "      <td>228.883492</td>\n",
       "      <td>227.054333</td>\n",
       "      <td>223.771350</td>\n",
       "      <td>221.138301</td>\n",
       "      <td>31173</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73</th>\n",
       "      <td>217.701407</td>\n",
       "      <td>222.907841</td>\n",
       "      <td>234.379332</td>\n",
       "      <td>226.419808</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>233.794670</td>\n",
       "      <td>248.161668</td>\n",
       "      <td>...</td>\n",
       "      <td>339.556515</td>\n",
       "      <td>296.289542</td>\n",
       "      <td>284.574216</td>\n",
       "      <td>278.489692</td>\n",
       "      <td>262.741244</td>\n",
       "      <td>229.943484</td>\n",
       "      <td>249.324856</td>\n",
       "      <td>257.722801</td>\n",
       "      <td>228.635842</td>\n",
       "      <td>31017</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74</th>\n",
       "      <td>200.953448</td>\n",
       "      <td>206.060110</td>\n",
       "      <td>216.937576</td>\n",
       "      <td>210.223460</td>\n",
       "      <td>206.205303</td>\n",
       "      <td>215.623962</td>\n",
       "      <td>210.083161</td>\n",
       "      <td>206.938707</td>\n",
       "      <td>211.231238</td>\n",
       "      <td>226.236504</td>\n",
       "      <td>...</td>\n",
       "      <td>293.819170</td>\n",
       "      <td>254.291571</td>\n",
       "      <td>242.991075</td>\n",
       "      <td>234.911629</td>\n",
       "      <td>231.250376</td>\n",
       "      <td>222.969831</td>\n",
       "      <td>225.899550</td>\n",
       "      <td>224.839268</td>\n",
       "      <td>217.664702</td>\n",
       "      <td>31041</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>213.306554</td>\n",
       "      <td>214.694366</td>\n",
       "      <td>209.325423</td>\n",
       "      <td>218.592480</td>\n",
       "      <td>211.894921</td>\n",
       "      <td>NaN</td>\n",
       "      <td>212.138548</td>\n",
       "      <td>218.947000</td>\n",
       "      <td>213.892217</td>\n",
       "      <td>227.047417</td>\n",
       "      <td>...</td>\n",
       "      <td>320.738291</td>\n",
       "      <td>287.761072</td>\n",
       "      <td>281.759466</td>\n",
       "      <td>263.291155</td>\n",
       "      <td>255.083175</td>\n",
       "      <td>246.266212</td>\n",
       "      <td>242.163655</td>\n",
       "      <td>240.368819</td>\n",
       "      <td>230.440891</td>\n",
       "      <td>31107</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>208.200017</td>\n",
       "      <td>NaN</td>\n",
       "      <td>211.188890</td>\n",
       "      <td>221.151065</td>\n",
       "      <td>210.418216</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>212.462624</td>\n",
       "      <td>205.740099</td>\n",
       "      <td>219.780829</td>\n",
       "      <td>...</td>\n",
       "      <td>264.865892</td>\n",
       "      <td>239.964399</td>\n",
       "      <td>242.297777</td>\n",
       "      <td>236.567732</td>\n",
       "      <td>235.464354</td>\n",
       "      <td>227.669941</td>\n",
       "      <td>225.798336</td>\n",
       "      <td>226.513042</td>\n",
       "      <td>221.618504</td>\n",
       "      <td>31021</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>220.128334</td>\n",
       "      <td>NaN</td>\n",
       "      <td>216.613270</td>\n",
       "      <td>232.370776</td>\n",
       "      <td>223.800229</td>\n",
       "      <td>235.779874</td>\n",
       "      <td>218.016180</td>\n",
       "      <td>228.569673</td>\n",
       "      <td>232.864794</td>\n",
       "      <td>244.416168</td>\n",
       "      <td>...</td>\n",
       "      <td>258.754519</td>\n",
       "      <td>222.692285</td>\n",
       "      <td>223.154147</td>\n",
       "      <td>230.851573</td>\n",
       "      <td>225.182302</td>\n",
       "      <td>219.624896</td>\n",
       "      <td>229.315692</td>\n",
       "      <td>229.157956</td>\n",
       "      <td>224.787351</td>\n",
       "      <td>31081</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>222.951654</td>\n",
       "      <td>NaN</td>\n",
       "      <td>225.046722</td>\n",
       "      <td>231.514478</td>\n",
       "      <td>223.368601</td>\n",
       "      <td>222.556728</td>\n",
       "      <td>229.929421</td>\n",
       "      <td>234.766896</td>\n",
       "      <td>227.221551</td>\n",
       "      <td>260.528451</td>\n",
       "      <td>...</td>\n",
       "      <td>291.217604</td>\n",
       "      <td>264.728758</td>\n",
       "      <td>262.457621</td>\n",
       "      <td>258.267617</td>\n",
       "      <td>254.196341</td>\n",
       "      <td>243.819473</td>\n",
       "      <td>248.228371</td>\n",
       "      <td>245.972857</td>\n",
       "      <td>239.716396</td>\n",
       "      <td>31109</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>217.982175</td>\n",
       "      <td>184.018642</td>\n",
       "      <td>198.955474</td>\n",
       "      <td>217.596965</td>\n",
       "      <td>211.986353</td>\n",
       "      <td>NaN</td>\n",
       "      <td>225.969618</td>\n",
       "      <td>203.002012</td>\n",
       "      <td>211.518187</td>\n",
       "      <td>223.541789</td>\n",
       "      <td>...</td>\n",
       "      <td>251.286455</td>\n",
       "      <td>228.920519</td>\n",
       "      <td>230.119585</td>\n",
       "      <td>230.566416</td>\n",
       "      <td>229.161100</td>\n",
       "      <td>226.782593</td>\n",
       "      <td>222.754595</td>\n",
       "      <td>222.980816</td>\n",
       "      <td>219.718988</td>\n",
       "      <td>31141</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>219.152881</td>\n",
       "      <td>NaN</td>\n",
       "      <td>215.921434</td>\n",
       "      <td>220.432367</td>\n",
       "      <td>207.625292</td>\n",
       "      <td>224.495043</td>\n",
       "      <td>216.487530</td>\n",
       "      <td>215.904011</td>\n",
       "      <td>219.848048</td>\n",
       "      <td>238.306663</td>\n",
       "      <td>...</td>\n",
       "      <td>301.366033</td>\n",
       "      <td>270.745413</td>\n",
       "      <td>263.500023</td>\n",
       "      <td>255.723882</td>\n",
       "      <td>256.227198</td>\n",
       "      <td>249.694215</td>\n",
       "      <td>253.619100</td>\n",
       "      <td>249.240000</td>\n",
       "      <td>241.869105</td>\n",
       "      <td>31181</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>215.424432</td>\n",
       "      <td>220.790961</td>\n",
       "      <td>213.146858</td>\n",
       "      <td>215.894533</td>\n",
       "      <td>217.583988</td>\n",
       "      <td>208.400401</td>\n",
       "      <td>210.195921</td>\n",
       "      <td>214.509242</td>\n",
       "      <td>219.416691</td>\n",
       "      <td>230.945432</td>\n",
       "      <td>...</td>\n",
       "      <td>267.739736</td>\n",
       "      <td>213.056601</td>\n",
       "      <td>236.075802</td>\n",
       "      <td>230.081746</td>\n",
       "      <td>216.143480</td>\n",
       "      <td>213.371444</td>\n",
       "      <td>211.261474</td>\n",
       "      <td>206.096197</td>\n",
       "      <td>207.776389</td>\n",
       "      <td>31033</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>220.615949</td>\n",
       "      <td>NaN</td>\n",
       "      <td>218.590481</td>\n",
       "      <td>230.384907</td>\n",
       "      <td>215.910038</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>228.329349</td>\n",
       "      <td>221.630601</td>\n",
       "      <td>247.241145</td>\n",
       "      <td>...</td>\n",
       "      <td>279.716475</td>\n",
       "      <td>250.523619</td>\n",
       "      <td>251.056605</td>\n",
       "      <td>244.642549</td>\n",
       "      <td>242.791943</td>\n",
       "      <td>225.357400</td>\n",
       "      <td>234.996204</td>\n",
       "      <td>234.952977</td>\n",
       "      <td>229.293571</td>\n",
       "      <td>31055</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>228.968772</td>\n",
       "      <td>NaN</td>\n",
       "      <td>226.616642</td>\n",
       "      <td>234.626317</td>\n",
       "      <td>220.825939</td>\n",
       "      <td>218.560240</td>\n",
       "      <td>236.098840</td>\n",
       "      <td>242.943334</td>\n",
       "      <td>241.149231</td>\n",
       "      <td>255.129856</td>\n",
       "      <td>...</td>\n",
       "      <td>322.789305</td>\n",
       "      <td>300.595578</td>\n",
       "      <td>282.576496</td>\n",
       "      <td>275.210792</td>\n",
       "      <td>270.903471</td>\n",
       "      <td>259.716230</td>\n",
       "      <td>259.992881</td>\n",
       "      <td>262.652732</td>\n",
       "      <td>257.772853</td>\n",
       "      <td>31097</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84</th>\n",
       "      <td>217.122714</td>\n",
       "      <td>NaN</td>\n",
       "      <td>220.588667</td>\n",
       "      <td>224.793355</td>\n",
       "      <td>216.715185</td>\n",
       "      <td>201.254937</td>\n",
       "      <td>NaN</td>\n",
       "      <td>223.499036</td>\n",
       "      <td>222.254405</td>\n",
       "      <td>246.718066</td>\n",
       "      <td>...</td>\n",
       "      <td>297.772530</td>\n",
       "      <td>269.055244</td>\n",
       "      <td>263.614309</td>\n",
       "      <td>253.572545</td>\n",
       "      <td>250.900770</td>\n",
       "      <td>234.525439</td>\n",
       "      <td>240.677298</td>\n",
       "      <td>237.778291</td>\n",
       "      <td>233.653822</td>\n",
       "      <td>31153</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>85</th>\n",
       "      <td>190.786923</td>\n",
       "      <td>194.258377</td>\n",
       "      <td>205.254833</td>\n",
       "      <td>206.168567</td>\n",
       "      <td>204.464054</td>\n",
       "      <td>NaN</td>\n",
       "      <td>193.423589</td>\n",
       "      <td>200.672699</td>\n",
       "      <td>196.276665</td>\n",
       "      <td>222.467141</td>\n",
       "      <td>...</td>\n",
       "      <td>283.313799</td>\n",
       "      <td>248.855616</td>\n",
       "      <td>248.210045</td>\n",
       "      <td>237.006558</td>\n",
       "      <td>235.965975</td>\n",
       "      <td>231.553569</td>\n",
       "      <td>227.879886</td>\n",
       "      <td>220.657325</td>\n",
       "      <td>219.124405</td>\n",
       "      <td>31043</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86</th>\n",
       "      <td>226.903136</td>\n",
       "      <td>NaN</td>\n",
       "      <td>217.298721</td>\n",
       "      <td>232.690471</td>\n",
       "      <td>226.693745</td>\n",
       "      <td>223.294759</td>\n",
       "      <td>232.200741</td>\n",
       "      <td>237.514793</td>\n",
       "      <td>234.626024</td>\n",
       "      <td>250.220679</td>\n",
       "      <td>...</td>\n",
       "      <td>308.759891</td>\n",
       "      <td>275.880833</td>\n",
       "      <td>271.737378</td>\n",
       "      <td>264.195565</td>\n",
       "      <td>259.092307</td>\n",
       "      <td>249.159207</td>\n",
       "      <td>251.684872</td>\n",
       "      <td>250.699740</td>\n",
       "      <td>247.217644</td>\n",
       "      <td>31131</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>87</th>\n",
       "      <td>222.816864</td>\n",
       "      <td>NaN</td>\n",
       "      <td>219.793983</td>\n",
       "      <td>233.542474</td>\n",
       "      <td>218.842753</td>\n",
       "      <td>234.712354</td>\n",
       "      <td>233.139654</td>\n",
       "      <td>251.989267</td>\n",
       "      <td>235.842948</td>\n",
       "      <td>247.520496</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>260.048591</td>\n",
       "      <td>254.777771</td>\n",
       "      <td>249.707581</td>\n",
       "      <td>243.314528</td>\n",
       "      <td>237.824675</td>\n",
       "      <td>243.525288</td>\n",
       "      <td>241.774603</td>\n",
       "      <td>240.562567</td>\n",
       "      <td>31147</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>88</th>\n",
       "      <td>223.846551</td>\n",
       "      <td>NaN</td>\n",
       "      <td>210.670390</td>\n",
       "      <td>232.879060</td>\n",
       "      <td>225.152223</td>\n",
       "      <td>215.899163</td>\n",
       "      <td>NaN</td>\n",
       "      <td>232.789696</td>\n",
       "      <td>231.620330</td>\n",
       "      <td>248.684776</td>\n",
       "      <td>...</td>\n",
       "      <td>264.927246</td>\n",
       "      <td>233.460208</td>\n",
       "      <td>239.091626</td>\n",
       "      <td>238.116686</td>\n",
       "      <td>233.232793</td>\n",
       "      <td>224.158004</td>\n",
       "      <td>232.405534</td>\n",
       "      <td>231.976091</td>\n",
       "      <td>230.493777</td>\n",
       "      <td>31023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>216.186418</td>\n",
       "      <td>NaN</td>\n",
       "      <td>216.345459</td>\n",
       "      <td>216.085532</td>\n",
       "      <td>216.390580</td>\n",
       "      <td>214.583651</td>\n",
       "      <td>214.641276</td>\n",
       "      <td>217.125215</td>\n",
       "      <td>222.123433</td>\n",
       "      <td>242.547982</td>\n",
       "      <td>...</td>\n",
       "      <td>260.608103</td>\n",
       "      <td>227.217177</td>\n",
       "      <td>218.290770</td>\n",
       "      <td>211.707467</td>\n",
       "      <td>210.252033</td>\n",
       "      <td>211.601584</td>\n",
       "      <td>212.450116</td>\n",
       "      <td>211.870831</td>\n",
       "      <td>208.910666</td>\n",
       "      <td>31083</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>90</th>\n",
       "      <td>206.472624</td>\n",
       "      <td>NaN</td>\n",
       "      <td>212.636108</td>\n",
       "      <td>215.304622</td>\n",
       "      <td>195.865098</td>\n",
       "      <td>206.871345</td>\n",
       "      <td>217.701428</td>\n",
       "      <td>214.915451</td>\n",
       "      <td>218.273763</td>\n",
       "      <td>236.500500</td>\n",
       "      <td>...</td>\n",
       "      <td>233.420079</td>\n",
       "      <td>220.608166</td>\n",
       "      <td>212.326673</td>\n",
       "      <td>207.429167</td>\n",
       "      <td>207.235747</td>\n",
       "      <td>205.763075</td>\n",
       "      <td>210.048064</td>\n",
       "      <td>209.705287</td>\n",
       "      <td>206.342136</td>\n",
       "      <td>31073</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>224.170237</td>\n",
       "      <td>215.308356</td>\n",
       "      <td>225.367747</td>\n",
       "      <td>213.631549</td>\n",
       "      <td>217.226432</td>\n",
       "      <td>201.458804</td>\n",
       "      <td>209.102443</td>\n",
       "      <td>201.641332</td>\n",
       "      <td>220.279505</td>\n",
       "      <td>230.839008</td>\n",
       "      <td>...</td>\n",
       "      <td>303.016449</td>\n",
       "      <td>268.867140</td>\n",
       "      <td>258.528966</td>\n",
       "      <td>245.403212</td>\n",
       "      <td>239.077052</td>\n",
       "      <td>228.004627</td>\n",
       "      <td>233.510711</td>\n",
       "      <td>232.537016</td>\n",
       "      <td>226.737241</td>\n",
       "      <td>31175</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92</th>\n",
       "      <td>205.750940</td>\n",
       "      <td>NaN</td>\n",
       "      <td>214.851819</td>\n",
       "      <td>207.865020</td>\n",
       "      <td>210.824867</td>\n",
       "      <td>214.473123</td>\n",
       "      <td>210.327268</td>\n",
       "      <td>213.663282</td>\n",
       "      <td>216.601491</td>\n",
       "      <td>237.016307</td>\n",
       "      <td>...</td>\n",
       "      <td>275.939034</td>\n",
       "      <td>249.405079</td>\n",
       "      <td>239.725840</td>\n",
       "      <td>233.203012</td>\n",
       "      <td>225.368129</td>\n",
       "      <td>222.835122</td>\n",
       "      <td>227.445873</td>\n",
       "      <td>226.993683</td>\n",
       "      <td>221.399025</td>\n",
       "      <td>31019</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>93 rows × 44 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    2002_025_NDVI  2002_033_NDVI  2002_041_NDVI  2002_049_NDVI  2002_057_NDVI  \\\n",
       "0      199.979369     195.341170     200.847186     202.875689            NaN   \n",
       "1      203.351533     201.146636     203.924538     216.894404            NaN   \n",
       "2             NaN            NaN            NaN            NaN            NaN   \n",
       "3      212.819195     213.835711     224.030314     216.884286     208.877590   \n",
       "4      221.803229            NaN     208.278847     234.653056     220.686890   \n",
       "5      229.916099     226.190997     227.608255     222.942409     212.259014   \n",
       "6      235.391085     214.114682     226.964055     216.516103     219.143350   \n",
       "7      212.339754     212.275480     213.647377     210.994667            NaN   \n",
       "8      226.551550            NaN     225.780873     237.574696     227.863166   \n",
       "9      222.940246            NaN     216.315210     238.430340     224.973504   \n",
       "10     213.087390     200.290834     197.167042     211.494392     203.471166   \n",
       "11     215.664492     208.997795     209.394583     217.713067     195.850396   \n",
       "12     229.262530     230.118864     238.688525     227.655732     245.322250   \n",
       "13     228.198784     217.080894     227.626282     231.060463            NaN   \n",
       "14            NaN     220.111292     228.955993     218.872942     221.335526   \n",
       "15     223.113451     216.537375     205.879802     220.499701     211.254640   \n",
       "16     204.493015     205.783988     214.464836     208.363281     204.308503   \n",
       "17     223.015067     216.410609     222.997272     215.303024     199.178759   \n",
       "18     211.463358            NaN     219.454761     231.509915     218.950347   \n",
       "19     222.050867     209.426620     212.194563     221.793128     221.055657   \n",
       "20     197.230012     196.013566     198.147255     196.737671            NaN   \n",
       "21     207.651102     201.204848     185.729315     207.247734     213.281543   \n",
       "22     242.666308     231.599703     235.794515     225.961006     234.321865   \n",
       "23     219.450249     212.416231     226.055764     219.538792     216.250584   \n",
       "24     232.153714            NaN     231.175874     237.683450     223.718070   \n",
       "25     211.420660     204.102347     219.626292     214.845444     205.364904   \n",
       "26     230.244046            NaN     222.098803     224.550970     215.997692   \n",
       "27     228.113689     224.909362     225.451196     228.623511     230.275897   \n",
       "28     229.464527            NaN     231.358043     240.091513     226.052552   \n",
       "29     212.331221     208.906524     204.555088     214.181301     209.629666   \n",
       "..            ...            ...            ...            ...            ...   \n",
       "63     193.814649     206.311882     215.171245     203.209433            NaN   \n",
       "64     214.917312     205.954823     222.972146     214.983491     207.958192   \n",
       "65     220.803198            NaN     219.848205     232.277201     222.299859   \n",
       "66     242.433150            NaN     227.342188     231.498883     222.243880   \n",
       "67     221.300050            NaN     203.774410     233.625819     223.150836   \n",
       "68     219.407627     203.312203     220.341442     220.498925     221.442048   \n",
       "69     212.861784            NaN     217.958204     233.059434     217.592601   \n",
       "70     223.332273            NaN     203.036219     235.870952     223.276035   \n",
       "71     197.877592     203.440447     200.271968     203.080493     218.074661   \n",
       "72     206.045084     198.668973     188.038388     210.902233     209.045800   \n",
       "73     217.701407     222.907841     234.379332     226.419808            NaN   \n",
       "74     200.953448     206.060110     216.937576     210.223460     206.205303   \n",
       "75     213.306554     214.694366     209.325423     218.592480     211.894921   \n",
       "76     208.200017            NaN     211.188890     221.151065     210.418216   \n",
       "77     220.128334            NaN     216.613270     232.370776     223.800229   \n",
       "78     222.951654            NaN     225.046722     231.514478     223.368601   \n",
       "79     217.982175     184.018642     198.955474     217.596965     211.986353   \n",
       "80     219.152881            NaN     215.921434     220.432367     207.625292   \n",
       "81     215.424432     220.790961     213.146858     215.894533     217.583988   \n",
       "82     220.615949            NaN     218.590481     230.384907     215.910038   \n",
       "83     228.968772            NaN     226.616642     234.626317     220.825939   \n",
       "84     217.122714            NaN     220.588667     224.793355     216.715185   \n",
       "85     190.786923     194.258377     205.254833     206.168567     204.464054   \n",
       "86     226.903136            NaN     217.298721     232.690471     226.693745   \n",
       "87     222.816864            NaN     219.793983     233.542474     218.842753   \n",
       "88     223.846551            NaN     210.670390     232.879060     225.152223   \n",
       "89     216.186418            NaN     216.345459     216.085532     216.390580   \n",
       "90     206.472624            NaN     212.636108     215.304622     195.865098   \n",
       "91     224.170237     215.308356     225.367747     213.631549     217.226432   \n",
       "92     205.750940            NaN     214.851819     207.865020     210.824867   \n",
       "\n",
       "    2002_065_NDVI  2002_073_NDVI  2002_081_NDVI  2002_089_NDVI  2002_097_NDVI  \\\n",
       "0      194.099093            NaN            NaN     208.873180     221.068704   \n",
       "1      203.419017            NaN            NaN     191.155802     204.691938   \n",
       "2             NaN            NaN            NaN            NaN            NaN   \n",
       "3      210.909783     213.257699     215.908769     215.358000     221.722039   \n",
       "4      232.823441     223.750504     231.863896     232.966169     244.505615   \n",
       "5      207.636691            NaN     206.848260     214.297155     227.841844   \n",
       "6             NaN     215.304442     213.190887     222.066034     239.968478   \n",
       "7      207.719078     217.716226     220.356139     212.470068     230.302013   \n",
       "8      221.063138     236.405947     239.055175     236.669350     257.192697   \n",
       "9      236.699749     229.727361     238.107216     238.114405     250.240022   \n",
       "10            NaN     200.169625            NaN     204.162609     218.450447   \n",
       "11            NaN     216.781553            NaN     212.394788     217.624835   \n",
       "12            NaN     217.796205     221.179550     227.242232     247.042382   \n",
       "13     168.647194            NaN            NaN     220.692757     251.722046   \n",
       "14            NaN     216.722434     216.859332     220.246166     236.874035   \n",
       "15            NaN     227.043816     197.070824     214.770376     225.288235   \n",
       "16            NaN            NaN     208.424568     201.091952     225.252987   \n",
       "17     222.789958     225.640728     224.922155     224.525629     237.756216   \n",
       "18     230.677036     223.316829     227.962304     233.859957     248.770612   \n",
       "19            NaN            NaN     212.024915     215.886264     242.860799   \n",
       "20     202.529293     200.720407            NaN     209.780270     218.362723   \n",
       "21            NaN     176.103800     209.581853     200.597896     227.037783   \n",
       "22     230.243851     232.960386     228.276336     236.611572     249.283510   \n",
       "23     223.426848     226.041834     223.943916     225.558735     245.639428   \n",
       "24     224.490121     235.426831     240.613085     236.556450     268.807481   \n",
       "25     220.773557     221.179931     189.675670     211.559741     219.927054   \n",
       "26     224.721818     224.347877     218.873843     223.009762     239.527747   \n",
       "27     217.735140     188.953720            NaN     230.458329     246.522200   \n",
       "28     213.447005     235.330268     255.875778     241.406905     270.577516   \n",
       "29            NaN            NaN     207.101442     208.294738     226.756304   \n",
       "..            ...            ...            ...            ...            ...   \n",
       "63            NaN            NaN            NaN     219.551935     233.150780   \n",
       "64     210.273061     217.328152     213.658577     218.705533     231.241189   \n",
       "65     230.929156            NaN     227.791331     233.589550     245.372966   \n",
       "66     221.989210     225.755707     231.777910     232.180752     250.745227   \n",
       "67     231.907039            NaN     228.982839     232.746413     246.538325   \n",
       "68     212.509602     221.038805     216.600120     221.534680     244.483284   \n",
       "69     222.831417            NaN     227.692886     223.013312     247.533964   \n",
       "70     222.781461     231.215785     234.651766     234.085998     254.890991   \n",
       "71     204.468616            NaN     181.880683     197.219207     211.149236   \n",
       "72            NaN            NaN     203.862209     204.054247     228.091318   \n",
       "73            NaN            NaN            NaN     233.794670     248.161668   \n",
       "74     215.623962     210.083161     206.938707     211.231238     226.236504   \n",
       "75            NaN     212.138548     218.947000     213.892217     227.047417   \n",
       "76            NaN            NaN     212.462624     205.740099     219.780829   \n",
       "77     235.779874     218.016180     228.569673     232.864794     244.416168   \n",
       "78     222.556728     229.929421     234.766896     227.221551     260.528451   \n",
       "79            NaN     225.969618     203.002012     211.518187     223.541789   \n",
       "80     224.495043     216.487530     215.904011     219.848048     238.306663   \n",
       "81     208.400401     210.195921     214.509242     219.416691     230.945432   \n",
       "82            NaN            NaN     228.329349     221.630601     247.241145   \n",
       "83     218.560240     236.098840     242.943334     241.149231     255.129856   \n",
       "84     201.254937            NaN     223.499036     222.254405     246.718066   \n",
       "85            NaN     193.423589     200.672699     196.276665     222.467141   \n",
       "86     223.294759     232.200741     237.514793     234.626024     250.220679   \n",
       "87     234.712354     233.139654     251.989267     235.842948     247.520496   \n",
       "88     215.899163            NaN     232.789696     231.620330     248.684776   \n",
       "89     214.583651     214.641276     217.125215     222.123433     242.547982   \n",
       "90     206.871345     217.701428     214.915451     218.273763     236.500500   \n",
       "91     201.458804     209.102443     201.641332     220.279505     230.839008   \n",
       "92     214.473123     210.327268     213.663282     216.601491     237.016307   \n",
       "\n",
       "    ...    2002_297_NDVI  2002_305_NDVI  2002_313_NDVI  2002_321_NDVI  \\\n",
       "0   ...       281.463133            NaN     249.743801     237.984343   \n",
       "1   ...       301.913183     266.251160     239.458615     250.944560   \n",
       "2   ...              NaN            NaN            NaN            NaN   \n",
       "3   ...              NaN     253.638928     238.831443     234.096765   \n",
       "4   ...       259.154760     220.426673     226.173553     235.380445   \n",
       "5   ...       296.791499     342.314950     304.395678     285.335320   \n",
       "6   ...       298.770726     274.771954     255.840324     246.847215   \n",
       "7   ...       302.585162     261.858709     251.620920     240.537196   \n",
       "8   ...       269.260749     248.306612     239.132158     243.130186   \n",
       "9   ...       243.716207     226.287687     217.557427     227.190027   \n",
       "10  ...       276.328730     241.475141     241.116998     234.251233   \n",
       "11  ...       296.039235     267.350001     251.110718     246.170961   \n",
       "12  ...       311.530339     283.933175     275.990437     268.327353   \n",
       "13  ...              NaN     351.530609     479.382004     367.240463   \n",
       "14  ...       305.912349     279.553929     268.892267     253.810679   \n",
       "15  ...       293.473458     260.204741     253.929321     244.339834   \n",
       "16  ...       288.203969     262.108939     255.520184     245.465956   \n",
       "17  ...       293.017282     275.189193     259.419199     245.260745   \n",
       "18  ...       255.214874     236.707776     224.221904     228.630716   \n",
       "19  ...       272.912082     246.144415     248.588913     243.150218   \n",
       "20  ...       297.675962     242.170773     232.205259     224.833441   \n",
       "21  ...       290.014531     263.704238     257.658566     243.814411   \n",
       "22  ...       388.105254     295.129787     284.616600     264.080364   \n",
       "23  ...              NaN     238.529693     229.845925     222.964434   \n",
       "24  ...       284.272495     260.207541     254.054105     254.585254   \n",
       "25  ...       293.056082     261.224720     252.441535     245.945802   \n",
       "26  ...       298.416446     263.235406     260.235061     251.679767   \n",
       "27  ...       281.771326     264.266931     248.872840     242.603717   \n",
       "28  ...       304.303369     288.503341     277.611749     274.740986   \n",
       "29  ...       297.018123     262.781072     261.065115     246.609795   \n",
       "..  ...              ...            ...            ...            ...   \n",
       "63  ...       291.935394     291.021334     280.745745     274.320888   \n",
       "64  ...       289.013780     264.161719     246.881307     237.388786   \n",
       "65  ...       282.279564     247.758407     248.450048     244.555951   \n",
       "66  ...       269.413014     252.192130     242.588541     238.737482   \n",
       "67  ...       264.688497     232.829392     235.869030     238.202210   \n",
       "68  ...              NaN     243.781136     234.991093     228.454418   \n",
       "69  ...       275.360540     243.966690     251.070115     247.055581   \n",
       "70  ...       268.962974     237.391334     242.550593     243.988495   \n",
       "71  ...       279.002013     232.825632     243.402558     224.848668   \n",
       "72  ...       272.902348     241.825806     241.467507     234.108402   \n",
       "73  ...       339.556515     296.289542     284.574216     278.489692   \n",
       "74  ...       293.819170     254.291571     242.991075     234.911629   \n",
       "75  ...       320.738291     287.761072     281.759466     263.291155   \n",
       "76  ...       264.865892     239.964399     242.297777     236.567732   \n",
       "77  ...       258.754519     222.692285     223.154147     230.851573   \n",
       "78  ...       291.217604     264.728758     262.457621     258.267617   \n",
       "79  ...       251.286455     228.920519     230.119585     230.566416   \n",
       "80  ...       301.366033     270.745413     263.500023     255.723882   \n",
       "81  ...       267.739736     213.056601     236.075802     230.081746   \n",
       "82  ...       279.716475     250.523619     251.056605     244.642549   \n",
       "83  ...       322.789305     300.595578     282.576496     275.210792   \n",
       "84  ...       297.772530     269.055244     263.614309     253.572545   \n",
       "85  ...       283.313799     248.855616     248.210045     237.006558   \n",
       "86  ...       308.759891     275.880833     271.737378     264.195565   \n",
       "87  ...              NaN     260.048591     254.777771     249.707581   \n",
       "88  ...       264.927246     233.460208     239.091626     238.116686   \n",
       "89  ...       260.608103     227.217177     218.290770     211.707467   \n",
       "90  ...       233.420079     220.608166     212.326673     207.429167   \n",
       "91  ...       303.016449     268.867140     258.528966     245.403212   \n",
       "92  ...       275.939034     249.405079     239.725840     233.203012   \n",
       "\n",
       "    2002_329_NDVI  2002_337_NDVI  2002_345_NDVI  2002_353_NDVI  2002_361_NDVI  \\\n",
       "0      229.278986     229.250037     219.821214     214.984558     212.731740   \n",
       "1      231.718090            NaN            NaN     214.261857     219.209543   \n",
       "2             NaN            NaN            NaN            NaN            NaN   \n",
       "3      226.862497     223.787780     222.051742     227.336760     221.165678   \n",
       "4      228.814447     223.299799     231.816201     232.165852     228.895232   \n",
       "5      253.459906     242.168471     236.214690     230.999540     234.907088   \n",
       "6      249.119410     235.243378     237.865975     233.151520     226.999504   \n",
       "7      226.203985     227.249151     225.747537     217.432238     217.896430   \n",
       "8      240.258149     235.511602     238.513087     239.084113     236.272009   \n",
       "9      225.817783     223.190388     228.394813     229.302264     226.193111   \n",
       "10     232.884580     228.291128     224.764699     220.751375     220.066632   \n",
       "11     238.863715     228.058503     233.839241     231.545832     229.893855   \n",
       "12     262.677335     451.712677     260.192161     250.888916     232.393577   \n",
       "13     391.930695            NaN     353.646042     425.709457     398.537155   \n",
       "14     250.104062     242.322301     238.880964     236.060772     228.101122   \n",
       "15     240.340041     234.601654     232.839708     227.979829     225.623739   \n",
       "16     241.372719     236.413167     233.739821     232.305500     224.970702   \n",
       "17     239.860045     238.887649     234.719141     230.851799     227.575202   \n",
       "18     225.810371     222.145856     228.175999     227.118750     225.065908   \n",
       "19     243.498574     238.942913     234.542641     234.752078     227.161276   \n",
       "20     213.230526     213.055090     211.939743     205.739942     208.560988   \n",
       "21     239.592426     239.531701     231.712770     229.048020     222.688131   \n",
       "22     259.291434     255.999083     244.524268     245.147904     237.348237   \n",
       "23     220.352473     219.302174     221.513899     218.431285     214.112700   \n",
       "24     251.166482     243.519961     248.215218     248.104461     243.690696   \n",
       "25     240.673315     228.889383     232.114628     232.490905     227.547505   \n",
       "26     239.487084     234.138533     240.196273     236.796861     231.038189   \n",
       "27     229.880533     233.534608     226.502091     223.341168     220.525165   \n",
       "28     269.581529     266.464531     261.839946     265.305399     263.990195   \n",
       "29     241.701473     236.014225     232.692441     228.584283     225.481419   \n",
       "..            ...            ...            ...            ...            ...   \n",
       "63     259.884211     252.452640     253.940006     261.101928     238.774878   \n",
       "64     230.002543     228.818851     227.221174     221.606303     224.714297   \n",
       "65     234.962705     226.276053     237.354817     234.900983     229.521826   \n",
       "66     236.149806     237.692255     236.785699     235.181311     233.594761   \n",
       "67     233.001977     224.561481     236.080358     234.148641     231.898808   \n",
       "68     224.015630     223.315872     223.337294     220.213319     217.213538   \n",
       "69     242.484065     230.180067     238.788935     238.341592     231.487470   \n",
       "70     237.599357     229.323270     236.509636     237.002648     233.973439   \n",
       "71     214.510684     208.763332     206.483264     204.298354     210.421736   \n",
       "72     234.985950     228.883492     227.054333     223.771350     221.138301   \n",
       "73     262.741244     229.943484     249.324856     257.722801     228.635842   \n",
       "74     231.250376     222.969831     225.899550     224.839268     217.664702   \n",
       "75     255.083175     246.266212     242.163655     240.368819     230.440891   \n",
       "76     235.464354     227.669941     225.798336     226.513042     221.618504   \n",
       "77     225.182302     219.624896     229.315692     229.157956     224.787351   \n",
       "78     254.196341     243.819473     248.228371     245.972857     239.716396   \n",
       "79     229.161100     226.782593     222.754595     222.980816     219.718988   \n",
       "80     256.227198     249.694215     253.619100     249.240000     241.869105   \n",
       "81     216.143480     213.371444     211.261474     206.096197     207.776389   \n",
       "82     242.791943     225.357400     234.996204     234.952977     229.293571   \n",
       "83     270.903471     259.716230     259.992881     262.652732     257.772853   \n",
       "84     250.900770     234.525439     240.677298     237.778291     233.653822   \n",
       "85     235.965975     231.553569     227.879886     220.657325     219.124405   \n",
       "86     259.092307     249.159207     251.684872     250.699740     247.217644   \n",
       "87     243.314528     237.824675     243.525288     241.774603     240.562567   \n",
       "88     233.232793     224.158004     232.405534     231.976091     230.493777   \n",
       "89     210.252033     211.601584     212.450116     211.870831     208.910666   \n",
       "90     207.235747     205.763075     210.048064     209.705287     206.342136   \n",
       "91     239.077052     228.004627     233.510711     232.537016     226.737241   \n",
       "92     225.368129     222.835122     227.445873     226.993683     221.399025   \n",
       "\n",
       "    GEOID  \n",
       "0   31007  \n",
       "1   31045  \n",
       "2   31075  \n",
       "3   31113  \n",
       "4   31185  \n",
       "5   31005  \n",
       "6   31071  \n",
       "7   31101  \n",
       "8   31151  \n",
       "9   31059  \n",
       "10  31119  \n",
       "11  31183  \n",
       "12  31009  \n",
       "13  31091  \n",
       "14  31115  \n",
       "15  31003  \n",
       "16  31027  \n",
       "17  31029  \n",
       "18  31035  \n",
       "19  31039  \n",
       "20  31049  \n",
       "21  31051  \n",
       "22  31057  \n",
       "23  31063  \n",
       "24  31067  \n",
       "25  31077  \n",
       "26  31093  \n",
       "27  31123  \n",
       "28  31133  \n",
       "29  31139  \n",
       "..    ...  \n",
       "63  31103  \n",
       "64  31111  \n",
       "65  31121  \n",
       "66  31129  \n",
       "67  31143  \n",
       "68  31145  \n",
       "69  31155  \n",
       "70  31159  \n",
       "71  31165  \n",
       "72  31173  \n",
       "73  31017  \n",
       "74  31041  \n",
       "75  31107  \n",
       "76  31021  \n",
       "77  31081  \n",
       "78  31109  \n",
       "79  31141  \n",
       "80  31181  \n",
       "81  31033  \n",
       "82  31055  \n",
       "83  31097  \n",
       "84  31153  \n",
       "85  31043  \n",
       "86  31131  \n",
       "87  31147  \n",
       "88  31023  \n",
       "89  31083  \n",
       "90  31073  \n",
       "91  31175  \n",
       "92  31019  \n",
       "\n",
       "[93 rows x 44 columns]"
      ]
     },
     "execution_count": 57,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.read_csv(\"county_NDVI_2002_corn.csv\").drop(\"Unnamed: 0\",axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Index([  31007, '31045', '31075', '31113', '31185', '31005', '31071', '31101',\n",
      "       '31151', '31059', '31119', '31183', '31009', '31091', '31115', '31003',\n",
      "       '31027', '31029', '31035', '31039', '31049', '31051', '31057', '31063',\n",
      "       '31067', '31077', '31093', '31123', '31133', '31139', '31149', '31163',\n",
      "       '31169', '31179', '31061', '31065', '31069', '31079', '31085', '31001',\n",
      "       '31011', '31015', '31089', '31095', '31099', '31105', '31117', '31125',\n",
      "       '31127', '31135', '31137', '31157', '31161', '31167', '31171', '31177',\n",
      "       '31013', '31025', '31031', '31037', '31047', '31053', '31087', '31103',\n",
      "       '31111', '31121', '31129', '31143', '31145', '31155', '31159', '31165',\n",
      "       '31173', '31017', '31041', '31107', '31021', '31081', '31109', '31141',\n",
      "       '31181', '31033', '31055', '31097', '31153', '31043', '31131', '31147',\n",
      "       '31023', '31083', '31073', '31175', '31019'],\n",
      "      dtype='object')\n",
      "Index([  31007, '31045', '31075', '31113', '31185', '31005', '31071', '31101',\n",
      "       '31151', '31059', '31119', '31183', '31009', '31091', '31115', '31003',\n",
      "       '31027', '31029', '31035', '31039', '31049', '31051', '31057', '31063',\n",
      "       '31067', '31077', '31093', '31123', '31133', '31139', '31149', '31163',\n",
      "       '31169', '31179', '31061', '31065', '31069', '31079', '31085', '31001',\n",
      "       '31011', '31015', '31089', '31095', '31099', '31105', '31117', '31125',\n",
      "       '31127', '31135', '31137', '31157', '31161', '31167', '31171', '31177',\n",
      "       '31013', '31025', '31031', '31037', '31047', '31053', '31087', '31103',\n",
      "       '31111', '31121', '31129', '31143', '31145', '31155', '31159', '31165',\n",
      "       '31173', '31017', '31041', '31107', '31021', '31081', '31109', '31141',\n",
      "       '31181', '31033', '31055', '31097', '31153', '31043', '31131', '31147',\n",
      "       '31023', '31083', '31073', '31175', '31019', '31007'],\n",
      "      dtype='object')\n"
     ]
    }
   ],
   "source": [
    "template=pd.read_csv('template.csv').drop('Unnamed: 0',axis=1)\n",
    "def crop_phenologies(file_path,crop):\n",
    "    p=pd.read_csv(file_path).drop(\"Unnamed: 0\",axis=1)\n",
    "    num_countys=p.shape[0]\n",
    "    num_days=p.shape[1]-1\n",
    "    i=0\n",
    "    num_days_missing=46-num_days\n",
    "    q=pd.melt(pd.Series.to_frame(p.iloc[i])).drop('variable',axis=1).drop([43],axis=0)\n",
    "    q.columns=[p.GEOID[i]]\n",
    "    #print(p.iloc[92])\n",
    "    #print(pd.melt(pd.Series.to_frame(p.iloc[0])).drop('variable',axis=1).drop([43],axis=0))\n",
    "    for i in range(1,num_countys):\n",
    "        q=pd.concat([q,pd.melt(pd.Series.to_frame(p.iloc[i])).drop('variable',axis=1).drop([43],axis=0)],axis=1)\n",
    "        q=q.rename(columns = {'value':str(p.GEOID[i])})\n",
    "    for i in range(num_days_missing):\n",
    "        q.loc[-1] = [np.nan]*93\n",
    "        q.index+=1\n",
    "    q.sort_index(inplace=True)\n",
    "    \n",
    "    print(q.columns)\n",
    "    df=pd.DataFrame()\n",
    "    for i in range(len(q)):\n",
    "        df=pd.concat([df,q.iloc[[i]],template])\n",
    "    df.index=range(1,1+len(df.index))\n",
    "    df=df[:365]\n",
    "    \n",
    "    if crop=='corn':\n",
    "            gl=140\n",
    "            gu=176\n",
    "            sl=213\n",
    "            su=305\n",
    "    elif crop=='soy':\n",
    "            gl=156\n",
    "            gu=196\n",
    "            sl=227\n",
    "            su=290\n",
    "    elif crop=='wheat':\n",
    "            gl=10\n",
    "            gu=60\n",
    "            sl=135\n",
    "            su=196\n",
    "    pairs=[]\n",
    "    print(df.columns)\n",
    "    #get_PTD(pd.Series.to_frame(df.iloc[:,93]),gl,gu,sl,su)\n",
    "    #for i in range(1,1+len(df.columns)):\n",
    "        #print(i)\n",
    "        #pairs.append(get_PTD(pd.Series.to_frame(df.iloc[:,i]),gl,gu,sl,su))\n",
    "    #return pairs\n",
    "crop_phenologies(\"county_NDVI_2002_corn.csv\",'corn')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import pdb\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from scipy.optimize import leastsq\n",
    "\n",
    "\n",
    "def model_fourier(params, agdd, n_harm):\n",
    "    \"\"\"\n",
    "    Fourier model\n",
    "    :param params:\n",
    "    :param agdd:\n",
    "    :param n_harm:\n",
    "    :return:\n",
    "    \"\"\"\n",
    "    integration_time = len(agdd)\n",
    "    t = np.arange(1, integration_time + 1)\n",
    "    result = t*.0 + params[0]\n",
    "    w = 1\n",
    "\n",
    "    for i in range(1, n_harm * 4, 4):\n",
    "        result = result + params[i] * np.cos(2.0 * np.pi * w * t / integration_time + params[i+1]) \\\n",
    "                 + params[i+2]*np.sin(2.0 * np.pi * w * t / integration_time + params[i+3])\n",
    "        w += 1\n",
    "\n",
    "    return result\n",
    "\n",
    "\n",
    "def mismatch_function(params, func_phenology, ndvi, agdd):\n",
    "    \"\"\"\n",
    "    The NDVI/Phenology model mismatch function\n",
    "    :param params:\n",
    "    :param func_phenology:\n",
    "    :param ndvi:\n",
    "    :param agdd:\n",
    "    :param years:\n",
    "    :return:\n",
    "    \"\"\"\n",
    "    # output stores the predictions\n",
    "    output = []\n",
    "\n",
    "    oot = ndvi - func_phenology(params, agdd, n_harm=8)\n",
    "\n",
    "    [output.append(x) for x in oot]\n",
    "\n",
    "    return np.array(output).squeeze()\n",
    "\n",
    "\n",
    "def do_fourier(ndvi, gdd, n_harm=8, init_params=None):\n",
    "    \"\"\"\n",
    "    :param ndvi:\n",
    "    :param gdd:\n",
    "    :param n_harm:\n",
    "    :param init_params:\n",
    "    :return:\n",
    "    \"\"\"\n",
    "    n_params = 1 + n_harm * 4\n",
    "\n",
    "    if init_params is None:\n",
    "        init_params = [.25, ] * n_params\n",
    "        (xsol, mesg) = leastsq(mismatch_function, init_params, args=(model_fourier, ndvi, gdd), maxfev=1000000)\n",
    "        model_fitted = model_fourier(xsol, gdd, n_harm)\n",
    "\n",
    "    return model_fitted\n",
    "\n",
    "\n",
    "def get_PTD(df,gl,gu,sl,su):\n",
    "    \"\"\"\n",
    "    Get phenological transition dates (greenup, senescence)\n",
    "    :param df:\n",
    "    :return:\n",
    "    \"\"\"\n",
    "    # Input dataframe has an index comprised of day of year and remaining columns signify NDVI\n",
    "    # Linearly interpolate dataframe columns to fill in missing values\n",
    "    \n",
    "    if df.isnull().values.all():\n",
    "        return 0,0\n",
    "    df = df.apply(pd.Series.interpolate)\n",
    "    first_valid_index=df.first_valid_index()\n",
    "    df=df.loc[first_valid_index:,:]\n",
    "    #print(df)\n",
    "    plt.plot(list(range(len(df))),df)\n",
    "    plt.xlabel(\"DOY\")\n",
    "    plt.ylabel(\"NDVI\")\n",
    "    plt.title(\"NDVI plot\")\n",
    "    #plt.show()\n",
    "    \n",
    "    # Now compute mean of all columns and get the smoothened NDVI\n",
    "    arr_smooth = do_fourier(df.iloc[:,0], [8.0] * len(df))\n",
    "    #print(arr_smooth)\n",
    "    plt.plot(list(range(len(arr_smooth))),arr_smooth)\n",
    "    plt.xlabel(\"DOY\")\n",
    "    plt.ylabel(\"NDVI\")\n",
    "    plt.title(\"NDVI smoothened plot\")\n",
    "    #plt.show()\n",
    "    \n",
    "    # For all other crops and regions, take differential\n",
    "    # To get doy_green, find the last occurrence of the max differential\n",
    "    diff_green = np.diff(arr_smooth[:365 + 1])\n",
    "    \n",
    "    plt.plot(list(range(len(diff_green))),diff_green)\n",
    "    plt.xlabel(\"DOY\")\n",
    "    plt.ylabel(\"NDVI differential\")\n",
    "    plt.title(\"NDVI diff_green plot\")\n",
    "    #plt.show()\n",
    "    #doy_green = np.where(diff_green == diff_green[gl:gu+1].max())[0][-1]\n",
    "    #doy_senesc = np.diff(arr_smooth[:365 + 1]).argmin()\n",
    "    green_indices=np.where(diff_green == diff_green[gl:gu+1].max())[0]\n",
    "    for i in green_indices:\n",
    "        if i<=gu and i>=gl:\n",
    "            doy_green = i\n",
    "            break\n",
    "    sen_indices=np.where(diff_green == diff_green[sl:su+1].min())[0]\n",
    "    for i in sen_indices:\n",
    "        if i<=su and i>=sl:\n",
    "            doy_senesc = i\n",
    "            break\n",
    "    return doy_green+first_valid_index, doy_senesc+first_valid_index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'df' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-15-db4b79ce6093>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0msl\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m213\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      4\u001b[0m \u001b[0msu\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m305\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 5\u001b[0;31m \u001b[0mget_PTD\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpd\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSeries\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mto_frame\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0miloc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m12\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mgl\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mgu\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0msl\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0msu\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m: name 'df' is not defined"
     ]
    }
   ],
   "source": [
    "gl=140\n",
    "gu=176\n",
    "sl=213\n",
    "su=305\n",
    "get_PTD(pd.Series.to_frame(df.iloc[:,12]),gl,gu,sl,su)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "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.5.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}