Large_Melt_Map-020822-Alaska.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Large_Melt_Map\n",
    "\n",
    "Calculate SIR and GRD melt onset dates for selected subset areas.\n",
    "\n",
    "Also calculates average across years.\n",
    "\n",
    "Saves SIR and GRD MOD data in pickle files\n",
    "\n",
    "Makes maps of various annual and/or average results.\n",
    "\n",
    "FIXME: \n",
    "The pixel overlays are simply a scatter plot with rectangular markers and are not using EASE-grid projection\n",
    "\n",
    "FIXME:\n",
    "Add a location for where to find the RGI overlays\n",
    "and read/display the RGI shapefiles"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Load in all the modules needed"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "%pylab notebook\n",
    "# check if a windows machine, it needs special attention\n",
    "# this extra step will bypass an error from mpl_toolkits.basemap\n",
    "import os\n",
    "if os.name == 'nt':\n",
    "    os.environ[\"PROJ_LIB\"] = os.path.join(os.environ[\"CONDA_PREFIX\"], \"Library\", \"share\")\n",
    "    os.environ[\"GDAL_DATA\"] = os.path.join(os.environ[\"CONDA_PREFIX\"], \"Library\", \"share\", \"gdal\")\n",
    "import matplotlib.pyplot as plt\n",
    "from mpl_toolkits.axes_grid1 import make_axes_locatable\n",
    "from netCDF4 import Dataset, num2date\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from pathlib import Path\n",
    "import re\n",
    "from cetbtools.ease2conv import Ease2Transform\n",
    "from mpl_toolkits.basemap import Basemap\n",
    "import glob\n",
    "from pathlib import Path"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "pd.set_option('display.max_rows', 999)\n",
    "pd.set_option('display.min_rows', 200)\n",
    "pd.set_option('display.max_columns', 10)\n",
    "pd.set_option('display.width', 120)\n",
    "pd.set_option('display.max_colwidth', 10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/mnt/data3/cetb/ipynb_melt_onset/scripts\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "('/mnt/data3/cetb/nsidc0630_v1/', 'Joan')"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Define the local machine location of CETB data cubes\n",
    "# This directory is expected to contain subdirectories in the following hierarchy\n",
    "# that duplicates the hierarchy on the Google Shared Drive NSIDC-SD-CETB/v1/, \n",
    "# for example:\n",
    "# dataDir/F13_SSMI/N/nc_cubes/cubes_<regionName>\n",
    "user = 'Joan' #Mariah #MJWindows\n",
    "if ('Joan' == user):\n",
    "    dataDir = '/mnt/data3/cetb/nsidc0630_v1/' #jmr machine fringe \n",
    "    scriptDir = '/mnt/data3/cetb/ipynb_melt_onset/scripts'\n",
    "elif ('Mariah' == user):\n",
    "    dataDir = Path(Path.home(), 'nsidc0630_v1') # Mariah's PC or Mary Jo's Mac\n",
    "    scriptDir = Path(Path.home(), 'ipynb_melt_onset', 'scripts')\n",
    "elif ('MJWindows' == user):\n",
    "    dataDir = Path('Z:/mj On My Mac/nsidc0630_v1') # Mary Jo's Windows machine\n",
    "    scriptDir = Path(Path.home(), 'ipynb_melt_onset', 'scripts')\n",
    "else:\n",
    "    raise ValueError(\"unknown user= %s\\n\" % (user) )\n",
    "    \n",
    "%cd $scriptDir\n",
    "dataDir, user"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# load the custom functions\n",
    "from CETB_IO import read_Tb_whole\n",
    "from CETB_IO import coords\n",
    "from CETB_algorithms import calc_DAV\n",
    "from CETB_IO import find_cube_offset\n",
    "from CETB_IO import grid_locations_of_subset\n",
    "from CETB_IO import years_for\n",
    "from CETB_algorithms import DAV_MOD\n",
    "from CETB_analysis import MOD_array"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Specify inputs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Specify region, satellite, sensor, channel, and image reconstruction algorithm of interest in file name\n",
    "# this notebook will read in 2 CETB datasets so that channels/algorithms/sensors can be compared\n",
    "region='AKYukon'  #'GLaIL'  #make this the same syntax as cubefilenames and sub-directory\n",
    "sat_GRD='F18'   #'AQUA' for AMSRE, 'F13','F14','F15'... for SSMI\n",
    "sat_SIR= 'F18'\n",
    "sensor_GRD='SSMIS'  #'AMSRE', 'SSMI', etc.\n",
    "sensor_SIR='SSMIS'\n",
    "channel_GRD='37V'  #'36V','36H', '18V','18H', etc. '19V','19H' and '37V','37H' for SSMI)\n",
    "channel_SIR='37V'\n",
    "alg_GRD='GRD'   #SIR or GRD\n",
    "alg_SIR='SIR'\n",
    "\n",
    "# set the sir to grd factor, depends on the channel\n",
    "if (re.match('^[389]', channel_GRD)):\n",
    "    sir_2_grd_factor = 8 # assume 3.125 km to 25 km\n",
    "elif (re.match('^[12]', channel_GRD)):\n",
    "    sir_2_grd_factor = 4 # assume 6.25 km to 25 km\n",
    "else:\n",
    "    raise ValueError(\"Cannot determine sir_2_grd_factor from channel %s\\n\" % (channel_GRD) )\n",
    "\n",
    "cubeType_GRD = channel_GRD + '-' + alg_GRD\n",
    "cubeType_SIR = channel_SIR + '-' + alg_SIR\n",
    "  \n",
    "if ('SSMI' == sensor_GRD) or ('SSMIS' == sensor_GRD):\n",
    "    provider='CSU' \n",
    "    version='v1.*'\n",
    "elif 'AMSRE' == sensor_GRD:\n",
    "    provider='RSS'\n",
    "    version='v1.3'\n",
    "\n",
    "hemName = 'N'    \n",
    "\n",
    "\n",
    "# on MJ's machine\n",
    "datadir_GRD = \"%s/%s_%s/%s/nc_cubes/cubes_%s/\" % (\n",
    "    dataDir, sat_GRD, sensor_GRD, hemName, region )\n",
    "datadir_SIR = \"%s/%s_%s/%s/nc_cubes/cubes_%s/\" % (\n",
    "    dataDir, sat_SIR, sensor_SIR, hemName, region )\n",
    "\n",
    "# prefix filepath\n",
    "prefix_GRD = 'CETB.cubefile.'+region+'.'+sat_GRD+'_'+sensor_GRD+'-'+channel_GRD+'-'+alg_GRD+'-'+provider+'-'+version\n",
    "prefix_SIR = 'CETB.cubefile.'+region+'.'+sat_SIR+'_'+sensor_SIR+'-'+channel_SIR+'-'+alg_SIR+'-'+provider+'-'+version\n",
    "\n",
    "Years=years_for(sat_GRD)\n",
    "#might want to truncate Years to subset if very slow during testing\n",
    "#if we give it more years than available what do we want it to do? \n",
    "#warn me but return what it finds\n",
    "\n",
    "#Enter a site name for titles of plots\n",
    "Site='Barrow, AK' #'Great Slave Lake, Canada'  #'Vatnajokull, Iceland'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# SPECIFY latitude and longitude in decimal degrees, need to choose lat/lon corners so that we will load\n",
    "# in a rectangle of pixels within the corners of these coordinates\n",
    "areaname= 'kuparuk' #'bathurst_range' #'gsl' #'hunza' #'vatna' #'barrow' #'kuparak'\n",
    "\n",
    "if ('vatna' == areaname):\n",
    "    lat_start=63.75  \n",
    "    lat_end=64.88    \n",
    "    lon_start=-20 \n",
    "    lon_end=-15  \n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Vatnajokull, Iceland'\n",
    "elif 'hunza' == areaname:\n",
    "    lat_start=35.9  \n",
    "    lat_end=37.1   \n",
    "    lon_start=74 \n",
    "    lon_end=76 \n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Hunza Basin'\n",
    "elif 'gsl' == areaname:\n",
    "    lat_start=59.00  \n",
    "    lat_end=67.00   \n",
    "    lon_start=-119.00 \n",
    "    lon_end=-107.00\n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Great Slave Lake, Canada'\n",
    "elif 'bathurst_range' == areaname:\n",
    "    lat_start=58.00  \n",
    "    lat_end=69.00   \n",
    "    lon_start=-125.00 \n",
    "    lon_end=-106.00\n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Bathurst Caribou Range, NWT'\n",
    "elif ('barrow' == areaname):\n",
    "    lat_start=69.50  \n",
    "    lat_end=71.50    \n",
    "    lon_start=-158 \n",
    "    lon_end=-152  \n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Barrow/Utkiagvik, AK'\n",
    "elif ('kuparuk' == areaname):\n",
    "    lat_start=68.50  \n",
    "    lat_end=70.50    \n",
    "    lon_start=-151 \n",
    "    lon_end=-148  \n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Kuparuk Basin, AK'\n",
    "elif ('sza' == areaname):\n",
    "    lat_start=81  \n",
    "    lat_end=80    \n",
    "    lon_start=92 \n",
    "    lon_end=98  \n",
    "    #Enter a site name for titles of plots\n",
    "    Site='Komsomolets Island AWS, SZ'\n",
    "else: \n",
    "    raise ValueError(\"Unknown area name=%s\" % (areaname) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading offset information from cubeFile=/mnt/data3/cetb/nsidc0630_v1//F18_SSMIS/N/nc_cubes/cubes_AKYukon/CETB.cubefile.AKYukon.F18_SSMIS-37V-SIR-CSU-v1.3.2010.TB.nc...\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(1791.9999999243946, 2063.999999943305)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# get cube offset for finding row/col\n",
    "# function is region specific\n",
    "find_cube_offset(region, cubeDir=datadir_SIR, cubeType=cubeType_SIR, verbose=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(53, 62, 52, 60)\n",
      "(424, 496, 416, 480)\n"
     ]
    }
   ],
   "source": [
    "# get the GRD pixel IDs for the lat/lon rectangle chosen\n",
    "# and then calculate the corrsponding SIR pixel row/col numbers\n",
    "rows_cols_GRD=coords(datadir_GRD, prefix_GRD, lat_start, lat_end, lon_start, lon_end)\n",
    "rows_cols_env = tuple(np.array(rows_cols_GRD) * sir_2_grd_factor)\n",
    "print(rows_cols_GRD)\n",
    "print(rows_cols_env)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Next filename=/mnt/data3/cetb/nsidc0630_v1//F18_SSMIS/N/nc_cubes/cubes_AKYukon/CETB.cubefile.AKYukon.F18_SSMIS-37V-SIR-CSU-v1.4.2018.TB.nc...\n",
      "Next filename=/mnt/data3/cetb/nsidc0630_v1//F18_SSMIS/N/nc_cubes/cubes_AKYukon/CETB.cubefile.AKYukon.F18_SSMIS-37V-SIR-CSU-v1.4.2019.TB.nc...\n",
      "Next filename=/mnt/data3/cetb/nsidc0630_v1//F18_SSMIS/N/nc_cubes/cubes_AKYukon/CETB.cubefile.AKYukon.F18_SSMIS-37V-GRD-CSU-v1.4.2018.TB.nc...\n",
      "Next filename=/mnt/data3/cetb/nsidc0630_v1//F18_SSMIS/N/nc_cubes/cubes_AKYukon/CETB.cubefile.AKYukon.F18_SSMIS-37V-GRD-CSU-v1.4.2019.TB.nc...\n"
     ]
    }
   ],
   "source": [
    "# load in SIR TB data\n",
    "# Truncate Years here for speed \n",
    "subYears = Years[8:10]\n",
    "\n",
    "data_SIR = read_Tb_whole(datadir_SIR, prefix_SIR, subYears,\n",
    "                         rows_cols_env[0], rows_cols_env[1], rows_cols_env[2], rows_cols_env[3])\n",
    "#CETB_SIR = data_SIR['TB']   # 3-D Tb time-series array of TB\n",
    "#cal_date = data_SIR['cal_date']    # 1-D array of dates, these will get passed to later functions\n",
    "#cal_year = data_SIR['cal_year']    # 1-D array of years\n",
    "#cal_month = data_SIR['cal_month']   # 1-D array of months\n",
    "\n",
    "# load GRD Tb data\n",
    "data_GRD = read_Tb_whole(datadir_GRD, prefix_GRD, subYears,\n",
    "                         rows_cols_GRD[0], rows_cols_GRD[1], rows_cols_GRD[2], rows_cols_GRD[3])\n",
    "#CETB_GRD=data_GRD['TB']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# calculate DAV for the Tb data that was imported\n",
    "DAV_SIR = calc_DAV(data_SIR['TB'])\n",
    "DAV_GRD = calc_DAV(data_GRD['TB'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## create arrays of MOD to be used for plotting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "newdata.shape (1460, 4608)\n",
      "moving flag array to newdata...\n",
      "number of days = 1460\n",
      "Next d = 0\n",
      "Next d = 100\n",
      "Next d = 200\n",
      "Next d = 300\n",
      "Next d = 400\n",
      "Next d = 500\n",
      "Next d = 600\n",
      "Next d = 700\n",
      "Next d = 800\n",
      "Next d = 900\n",
      "Next d = 1000\n",
      "Next d = 1100\n",
      "Next d = 1200\n",
      "Next d = 1300\n",
      "Next d = 1400\n",
      "dataFrame is ready with flag data\n",
      "doing rolling sums...\n",
      "Next year = 2018...\n",
      "Next year = 2019...\n",
      "Getting average MOD at each pixel...\n",
      "Setting geolocation information...\n"
     ]
    },
    {
     "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>pixel</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>row</th>\n",
       "      <th>column</th>\n",
       "      <th>2018</th>\n",
       "      <th>2019</th>\n",
       "      <th>Avg</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>424,416</td>\n",
       "      <td>-1248437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.189531</td>\n",
       "      <td>-148.94...</td>\n",
       "      <td>424</td>\n",
       "      <td>416</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>424,417</td>\n",
       "      <td>-1245312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.204235</td>\n",
       "      <td>-149.01...</td>\n",
       "      <td>424</td>\n",
       "      <td>417</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>424,418</td>\n",
       "      <td>-1242187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.218912</td>\n",
       "      <td>-149.07...</td>\n",
       "      <td>424</td>\n",
       "      <td>418</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>424,419</td>\n",
       "      <td>-1239062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.233561</td>\n",
       "      <td>-149.13...</td>\n",
       "      <td>424</td>\n",
       "      <td>419</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>424,420</td>\n",
       "      <td>-1235937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.248182</td>\n",
       "      <td>-149.20...</td>\n",
       "      <td>424</td>\n",
       "      <td>420</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>424,421</td>\n",
       "      <td>-1232812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.262776</td>\n",
       "      <td>-149.26...</td>\n",
       "      <td>424</td>\n",
       "      <td>421</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>424,422</td>\n",
       "      <td>-1229687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.277342</td>\n",
       "      <td>-149.32...</td>\n",
       "      <td>424</td>\n",
       "      <td>422</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>424,423</td>\n",
       "      <td>-1226562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.291880</td>\n",
       "      <td>-149.39...</td>\n",
       "      <td>424</td>\n",
       "      <td>423</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>424,424</td>\n",
       "      <td>-1223437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.306391</td>\n",
       "      <td>-149.45...</td>\n",
       "      <td>424</td>\n",
       "      <td>424</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>424,425</td>\n",
       "      <td>-1220312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.320873</td>\n",
       "      <td>-149.52...</td>\n",
       "      <td>424</td>\n",
       "      <td>425</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>424,426</td>\n",
       "      <td>-1217187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.335328</td>\n",
       "      <td>-149.58...</td>\n",
       "      <td>424</td>\n",
       "      <td>426</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>424,427</td>\n",
       "      <td>-1214062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.349755</td>\n",
       "      <td>-149.64...</td>\n",
       "      <td>424</td>\n",
       "      <td>427</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>424,428</td>\n",
       "      <td>-1210937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.364154</td>\n",
       "      <td>-149.71...</td>\n",
       "      <td>424</td>\n",
       "      <td>428</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>424,429</td>\n",
       "      <td>-1207812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.378525</td>\n",
       "      <td>-149.77...</td>\n",
       "      <td>424</td>\n",
       "      <td>429</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>424,430</td>\n",
       "      <td>-1204687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.392867</td>\n",
       "      <td>-149.84...</td>\n",
       "      <td>424</td>\n",
       "      <td>430</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>424,431</td>\n",
       "      <td>-1201562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.407182</td>\n",
       "      <td>-149.90...</td>\n",
       "      <td>424</td>\n",
       "      <td>431</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>424,432</td>\n",
       "      <td>-1198437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.421468</td>\n",
       "      <td>-149.97...</td>\n",
       "      <td>424</td>\n",
       "      <td>432</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>424,433</td>\n",
       "      <td>-1195312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.435726</td>\n",
       "      <td>-150.03...</td>\n",
       "      <td>424</td>\n",
       "      <td>433</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>424,434</td>\n",
       "      <td>-1192187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.449956</td>\n",
       "      <td>-150.10...</td>\n",
       "      <td>424</td>\n",
       "      <td>434</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>424,435</td>\n",
       "      <td>-1189062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.464157</td>\n",
       "      <td>-150.16...</td>\n",
       "      <td>424</td>\n",
       "      <td>435</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>424,436</td>\n",
       "      <td>-1185937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.478330</td>\n",
       "      <td>-150.23...</td>\n",
       "      <td>424</td>\n",
       "      <td>436</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>424,437</td>\n",
       "      <td>-1182812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.492475</td>\n",
       "      <td>-150.29...</td>\n",
       "      <td>424</td>\n",
       "      <td>437</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>424,438</td>\n",
       "      <td>-1179687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.506591</td>\n",
       "      <td>-150.36...</td>\n",
       "      <td>424</td>\n",
       "      <td>438</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>424,439</td>\n",
       "      <td>-1176562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.520678</td>\n",
       "      <td>-150.42...</td>\n",
       "      <td>424</td>\n",
       "      <td>439</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>424,440</td>\n",
       "      <td>-1173437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.534737</td>\n",
       "      <td>-150.49...</td>\n",
       "      <td>424</td>\n",
       "      <td>440</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>424,441</td>\n",
       "      <td>-1170312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.548767</td>\n",
       "      <td>-150.55...</td>\n",
       "      <td>424</td>\n",
       "      <td>441</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>424,442</td>\n",
       "      <td>-1167187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.562769</td>\n",
       "      <td>-150.62...</td>\n",
       "      <td>424</td>\n",
       "      <td>442</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>424,443</td>\n",
       "      <td>-1164062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.576741</td>\n",
       "      <td>-150.68...</td>\n",
       "      <td>424</td>\n",
       "      <td>443</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>424,444</td>\n",
       "      <td>-1160937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.590685</td>\n",
       "      <td>-150.75...</td>\n",
       "      <td>424</td>\n",
       "      <td>444</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>424,445</td>\n",
       "      <td>-1157812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.604600</td>\n",
       "      <td>-150.82...</td>\n",
       "      <td>424</td>\n",
       "      <td>445</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>424,446</td>\n",
       "      <td>-1154687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.618486</td>\n",
       "      <td>-150.88...</td>\n",
       "      <td>424</td>\n",
       "      <td>446</td>\n",
       "      <td>117</td>\n",
       "      <td>90</td>\n",
       "      <td>103.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>424,447</td>\n",
       "      <td>-1151562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.632343</td>\n",
       "      <td>-150.95...</td>\n",
       "      <td>424</td>\n",
       "      <td>447</td>\n",
       "      <td>117</td>\n",
       "      <td>90</td>\n",
       "      <td>103.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>424,448</td>\n",
       "      <td>-1148437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.646170</td>\n",
       "      <td>-151.01...</td>\n",
       "      <td>424</td>\n",
       "      <td>448</td>\n",
       "      <td>117</td>\n",
       "      <td>90</td>\n",
       "      <td>103.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>424,449</td>\n",
       "      <td>-1145312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.659969</td>\n",
       "      <td>-151.08...</td>\n",
       "      <td>424</td>\n",
       "      <td>449</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>424,450</td>\n",
       "      <td>-1142187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.673739</td>\n",
       "      <td>-151.15...</td>\n",
       "      <td>424</td>\n",
       "      <td>450</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>424,451</td>\n",
       "      <td>-1139062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.687479</td>\n",
       "      <td>-151.21...</td>\n",
       "      <td>424</td>\n",
       "      <td>451</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>424,452</td>\n",
       "      <td>-1135937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.701191</td>\n",
       "      <td>-151.28...</td>\n",
       "      <td>424</td>\n",
       "      <td>452</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>424,453</td>\n",
       "      <td>-1132812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.714872</td>\n",
       "      <td>-151.35...</td>\n",
       "      <td>424</td>\n",
       "      <td>453</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>424,454</td>\n",
       "      <td>-1129687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.728525</td>\n",
       "      <td>-151.41...</td>\n",
       "      <td>424</td>\n",
       "      <td>454</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>424,455</td>\n",
       "      <td>-1126562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.742148</td>\n",
       "      <td>-151.48...</td>\n",
       "      <td>424</td>\n",
       "      <td>455</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40</th>\n",
       "      <td>424,456</td>\n",
       "      <td>-1123437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.755742</td>\n",
       "      <td>-151.55...</td>\n",
       "      <td>424</td>\n",
       "      <td>456</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>424,457</td>\n",
       "      <td>-1120312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.769306</td>\n",
       "      <td>-151.61...</td>\n",
       "      <td>424</td>\n",
       "      <td>457</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>424,458</td>\n",
       "      <td>-1117187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.782840</td>\n",
       "      <td>-151.68...</td>\n",
       "      <td>424</td>\n",
       "      <td>458</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>424,459</td>\n",
       "      <td>-1114062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.796345</td>\n",
       "      <td>-151.75...</td>\n",
       "      <td>424</td>\n",
       "      <td>459</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>44</th>\n",
       "      <td>424,460</td>\n",
       "      <td>-1110937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.809820</td>\n",
       "      <td>-151.81...</td>\n",
       "      <td>424</td>\n",
       "      <td>460</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>424,461</td>\n",
       "      <td>-1107812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.823266</td>\n",
       "      <td>-151.88...</td>\n",
       "      <td>424</td>\n",
       "      <td>461</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>46</th>\n",
       "      <td>424,462</td>\n",
       "      <td>-1104687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.836681</td>\n",
       "      <td>-151.95...</td>\n",
       "      <td>424</td>\n",
       "      <td>462</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>47</th>\n",
       "      <td>424,463</td>\n",
       "      <td>-1101562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.850067</td>\n",
       "      <td>-152.01...</td>\n",
       "      <td>424</td>\n",
       "      <td>463</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48</th>\n",
       "      <td>424,464</td>\n",
       "      <td>-1098437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.863423</td>\n",
       "      <td>-152.08...</td>\n",
       "      <td>424</td>\n",
       "      <td>464</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>49</th>\n",
       "      <td>424,465</td>\n",
       "      <td>-1095312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.876749</td>\n",
       "      <td>-152.15...</td>\n",
       "      <td>424</td>\n",
       "      <td>465</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>424,466</td>\n",
       "      <td>-1092187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.890045</td>\n",
       "      <td>-152.22...</td>\n",
       "      <td>424</td>\n",
       "      <td>466</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>51</th>\n",
       "      <td>424,467</td>\n",
       "      <td>-1089062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.903311</td>\n",
       "      <td>-152.28...</td>\n",
       "      <td>424</td>\n",
       "      <td>467</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>52</th>\n",
       "      <td>424,468</td>\n",
       "      <td>-1085937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.916547</td>\n",
       "      <td>-152.35...</td>\n",
       "      <td>424</td>\n",
       "      <td>468</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>424,469</td>\n",
       "      <td>-1082812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.929753</td>\n",
       "      <td>-152.42...</td>\n",
       "      <td>424</td>\n",
       "      <td>469</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>54</th>\n",
       "      <td>424,470</td>\n",
       "      <td>-1079687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.942928</td>\n",
       "      <td>-152.49...</td>\n",
       "      <td>424</td>\n",
       "      <td>470</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55</th>\n",
       "      <td>424,471</td>\n",
       "      <td>-1076562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.956073</td>\n",
       "      <td>-152.56...</td>\n",
       "      <td>424</td>\n",
       "      <td>471</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>56</th>\n",
       "      <td>424,472</td>\n",
       "      <td>-1073437.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.969188</td>\n",
       "      <td>-152.62...</td>\n",
       "      <td>424</td>\n",
       "      <td>472</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>57</th>\n",
       "      <td>424,473</td>\n",
       "      <td>-1070312.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.982273</td>\n",
       "      <td>-152.69...</td>\n",
       "      <td>424</td>\n",
       "      <td>473</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td>424,474</td>\n",
       "      <td>-1067187.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>68.995327</td>\n",
       "      <td>-152.76...</td>\n",
       "      <td>424</td>\n",
       "      <td>474</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>424,475</td>\n",
       "      <td>-1064062.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>69.008350</td>\n",
       "      <td>-152.83...</td>\n",
       "      <td>424</td>\n",
       "      <td>475</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>424,476</td>\n",
       "      <td>-1060937.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>69.021343</td>\n",
       "      <td>-152.90...</td>\n",
       "      <td>424</td>\n",
       "      <td>476</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>424,477</td>\n",
       "      <td>-1057812.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>69.034306</td>\n",
       "      <td>-152.97...</td>\n",
       "      <td>424</td>\n",
       "      <td>477</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>424,478</td>\n",
       "      <td>-1054687.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>69.047238</td>\n",
       "      <td>-153.03...</td>\n",
       "      <td>424</td>\n",
       "      <td>478</td>\n",
       "      <td>115</td>\n",
       "      <td>119</td>\n",
       "      <td>117.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>424,479</td>\n",
       "      <td>-1051562.5</td>\n",
       "      <td>2073437.5</td>\n",
       "      <td>69.060139</td>\n",
       "      <td>-153.10...</td>\n",
       "      <td>424</td>\n",
       "      <td>479</td>\n",
       "      <td>116</td>\n",
       "      <td>119</td>\n",
       "      <td>117.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>64</th>\n",
       "      <td>425,416</td>\n",
       "      <td>-1248437.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.213969</td>\n",
       "      <td>-148.90...</td>\n",
       "      <td>425</td>\n",
       "      <td>416</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>425,417</td>\n",
       "      <td>-1245312.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.228689</td>\n",
       "      <td>-148.97...</td>\n",
       "      <td>425</td>\n",
       "      <td>417</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>66</th>\n",
       "      <td>425,418</td>\n",
       "      <td>-1242187.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.243381</td>\n",
       "      <td>-149.03...</td>\n",
       "      <td>425</td>\n",
       "      <td>418</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>67</th>\n",
       "      <td>425,419</td>\n",
       "      <td>-1239062.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.258046</td>\n",
       "      <td>-149.09...</td>\n",
       "      <td>425</td>\n",
       "      <td>419</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>68</th>\n",
       "      <td>425,420</td>\n",
       "      <td>-1235937.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.272683</td>\n",
       "      <td>-149.16...</td>\n",
       "      <td>425</td>\n",
       "      <td>420</td>\n",
       "      <td>129</td>\n",
       "      <td>121</td>\n",
       "      <td>125.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>425,421</td>\n",
       "      <td>-1232812.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.287292</td>\n",
       "      <td>-149.22...</td>\n",
       "      <td>425</td>\n",
       "      <td>421</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>70</th>\n",
       "      <td>425,422</td>\n",
       "      <td>-1229687.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.301874</td>\n",
       "      <td>-149.29...</td>\n",
       "      <td>425</td>\n",
       "      <td>422</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>71</th>\n",
       "      <td>425,423</td>\n",
       "      <td>-1226562.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.316428</td>\n",
       "      <td>-149.35...</td>\n",
       "      <td>425</td>\n",
       "      <td>423</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>72</th>\n",
       "      <td>425,424</td>\n",
       "      <td>-1223437.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.330954</td>\n",
       "      <td>-149.41...</td>\n",
       "      <td>425</td>\n",
       "      <td>424</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>73</th>\n",
       "      <td>425,425</td>\n",
       "      <td>-1220312.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.345452</td>\n",
       "      <td>-149.48...</td>\n",
       "      <td>425</td>\n",
       "      <td>425</td>\n",
       "      <td>129</td>\n",
       "      <td>120</td>\n",
       "      <td>124.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74</th>\n",
       "      <td>425,426</td>\n",
       "      <td>-1217187.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.359922</td>\n",
       "      <td>-149.54...</td>\n",
       "      <td>425</td>\n",
       "      <td>426</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>425,427</td>\n",
       "      <td>-1214062.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.374365</td>\n",
       "      <td>-149.61...</td>\n",
       "      <td>425</td>\n",
       "      <td>427</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>425,428</td>\n",
       "      <td>-1210937.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.388779</td>\n",
       "      <td>-149.67...</td>\n",
       "      <td>425</td>\n",
       "      <td>428</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>425,429</td>\n",
       "      <td>-1207812.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.403165</td>\n",
       "      <td>-149.74...</td>\n",
       "      <td>425</td>\n",
       "      <td>429</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>425,430</td>\n",
       "      <td>-1204687.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.417524</td>\n",
       "      <td>-149.80...</td>\n",
       "      <td>425</td>\n",
       "      <td>430</td>\n",
       "      <td>129</td>\n",
       "      <td>119</td>\n",
       "      <td>124.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>425,431</td>\n",
       "      <td>-1201562.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.431854</td>\n",
       "      <td>-149.87...</td>\n",
       "      <td>425</td>\n",
       "      <td>431</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>80</th>\n",
       "      <td>425,432</td>\n",
       "      <td>-1198437.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.446156</td>\n",
       "      <td>-149.93...</td>\n",
       "      <td>425</td>\n",
       "      <td>432</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>425,433</td>\n",
       "      <td>-1195312.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.460429</td>\n",
       "      <td>-149.99...</td>\n",
       "      <td>425</td>\n",
       "      <td>433</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>82</th>\n",
       "      <td>425,434</td>\n",
       "      <td>-1192187.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.474674</td>\n",
       "      <td>-150.06...</td>\n",
       "      <td>425</td>\n",
       "      <td>434</td>\n",
       "      <td>129</td>\n",
       "      <td>90</td>\n",
       "      <td>109.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>83</th>\n",
       "      <td>425,435</td>\n",
       "      <td>-1189062.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.488891</td>\n",
       "      <td>-150.12...</td>\n",
       "      <td>425</td>\n",
       "      <td>435</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>84</th>\n",
       "      <td>425,436</td>\n",
       "      <td>-1185937.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.503080</td>\n",
       "      <td>-150.19...</td>\n",
       "      <td>425</td>\n",
       "      <td>436</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>85</th>\n",
       "      <td>425,437</td>\n",
       "      <td>-1182812.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.517240</td>\n",
       "      <td>-150.25...</td>\n",
       "      <td>425</td>\n",
       "      <td>437</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>86</th>\n",
       "      <td>425,438</td>\n",
       "      <td>-1179687.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.531371</td>\n",
       "      <td>-150.32...</td>\n",
       "      <td>425</td>\n",
       "      <td>438</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>87</th>\n",
       "      <td>425,439</td>\n",
       "      <td>-1176562.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.545474</td>\n",
       "      <td>-150.39...</td>\n",
       "      <td>425</td>\n",
       "      <td>439</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>88</th>\n",
       "      <td>425,440</td>\n",
       "      <td>-1173437.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.559548</td>\n",
       "      <td>-150.45...</td>\n",
       "      <td>425</td>\n",
       "      <td>440</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>425,441</td>\n",
       "      <td>-1170312.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.573594</td>\n",
       "      <td>-150.52...</td>\n",
       "      <td>425</td>\n",
       "      <td>441</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>90</th>\n",
       "      <td>425,442</td>\n",
       "      <td>-1167187.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.587611</td>\n",
       "      <td>-150.58...</td>\n",
       "      <td>425</td>\n",
       "      <td>442</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>425,443</td>\n",
       "      <td>-1164062.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.601599</td>\n",
       "      <td>-150.65...</td>\n",
       "      <td>425</td>\n",
       "      <td>443</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92</th>\n",
       "      <td>425,444</td>\n",
       "      <td>-1160937.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.615558</td>\n",
       "      <td>-150.71...</td>\n",
       "      <td>425</td>\n",
       "      <td>444</td>\n",
       "      <td>128</td>\n",
       "      <td>90</td>\n",
       "      <td>109.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>93</th>\n",
       "      <td>425,445</td>\n",
       "      <td>-1157812.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.629488</td>\n",
       "      <td>-150.78...</td>\n",
       "      <td>425</td>\n",
       "      <td>445</td>\n",
       "      <td>117</td>\n",
       "      <td>90</td>\n",
       "      <td>103.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>94</th>\n",
       "      <td>425,446</td>\n",
       "      <td>-1154687.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.643389</td>\n",
       "      <td>-150.85...</td>\n",
       "      <td>425</td>\n",
       "      <td>446</td>\n",
       "      <td>117</td>\n",
       "      <td>90</td>\n",
       "      <td>103.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95</th>\n",
       "      <td>425,447</td>\n",
       "      <td>-1151562.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.657261</td>\n",
       "      <td>-150.91...</td>\n",
       "      <td>425</td>\n",
       "      <td>447</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>96</th>\n",
       "      <td>425,448</td>\n",
       "      <td>-1148437.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.671105</td>\n",
       "      <td>-150.98...</td>\n",
       "      <td>425</td>\n",
       "      <td>448</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>97</th>\n",
       "      <td>425,449</td>\n",
       "      <td>-1145312.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.684919</td>\n",
       "      <td>-151.04...</td>\n",
       "      <td>425</td>\n",
       "      <td>449</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>98</th>\n",
       "      <td>425,450</td>\n",
       "      <td>-1142187.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.698704</td>\n",
       "      <td>-151.11...</td>\n",
       "      <td>425</td>\n",
       "      <td>450</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>99</th>\n",
       "      <td>425,451</td>\n",
       "      <td>-1139062.5</td>\n",
       "      <td>2070312.5</td>\n",
       "      <td>68.712460</td>\n",
       "      <td>-151.18...</td>\n",
       "      <td>425</td>\n",
       "      <td>451</td>\n",
       "      <td>115</td>\n",
       "      <td>90</td>\n",
       "      <td>102.5</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>4508</th>\n",
       "      <td>494,444</td>\n",
       "      <td>-1160937.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.305772</td>\n",
       "      <td>-147.95...</td>\n",
       "      <td>494</td>\n",
       "      <td>444</td>\n",
       "      <td>116</td>\n",
       "      <td>89</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4509</th>\n",
       "      <td>494,445</td>\n",
       "      <td>-1157812.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.320841</td>\n",
       "      <td>-148.02...</td>\n",
       "      <td>494</td>\n",
       "      <td>445</td>\n",
       "      <td>116</td>\n",
       "      <td>89</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4510</th>\n",
       "      <td>494,446</td>\n",
       "      <td>-1154687.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.335881</td>\n",
       "      <td>-148.09...</td>\n",
       "      <td>494</td>\n",
       "      <td>446</td>\n",
       "      <td>116</td>\n",
       "      <td>89</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4511</th>\n",
       "      <td>494,447</td>\n",
       "      <td>-1151562.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.350890</td>\n",
       "      <td>-148.16...</td>\n",
       "      <td>494</td>\n",
       "      <td>447</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4512</th>\n",
       "      <td>494,448</td>\n",
       "      <td>-1148437.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.365870</td>\n",
       "      <td>-148.23...</td>\n",
       "      <td>494</td>\n",
       "      <td>448</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4513</th>\n",
       "      <td>494,449</td>\n",
       "      <td>-1145312.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.380820</td>\n",
       "      <td>-148.30...</td>\n",
       "      <td>494</td>\n",
       "      <td>449</td>\n",
       "      <td>134</td>\n",
       "      <td>89</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4514</th>\n",
       "      <td>494,450</td>\n",
       "      <td>-1142187.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.395741</td>\n",
       "      <td>-148.37...</td>\n",
       "      <td>494</td>\n",
       "      <td>450</td>\n",
       "      <td>134</td>\n",
       "      <td>89</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4515</th>\n",
       "      <td>494,451</td>\n",
       "      <td>-1139062.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.410631</td>\n",
       "      <td>-148.44...</td>\n",
       "      <td>494</td>\n",
       "      <td>451</td>\n",
       "      <td>133</td>\n",
       "      <td>89</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4516</th>\n",
       "      <td>494,452</td>\n",
       "      <td>-1135937.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.425492</td>\n",
       "      <td>-148.51...</td>\n",
       "      <td>494</td>\n",
       "      <td>452</td>\n",
       "      <td>133</td>\n",
       "      <td>89</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4517</th>\n",
       "      <td>494,453</td>\n",
       "      <td>-1132812.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.440322</td>\n",
       "      <td>-148.58...</td>\n",
       "      <td>494</td>\n",
       "      <td>453</td>\n",
       "      <td>133</td>\n",
       "      <td>89</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4518</th>\n",
       "      <td>494,454</td>\n",
       "      <td>-1129687.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.455122</td>\n",
       "      <td>-148.65...</td>\n",
       "      <td>494</td>\n",
       "      <td>454</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4519</th>\n",
       "      <td>494,455</td>\n",
       "      <td>-1126562.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.469892</td>\n",
       "      <td>-148.72...</td>\n",
       "      <td>494</td>\n",
       "      <td>455</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4520</th>\n",
       "      <td>494,456</td>\n",
       "      <td>-1123437.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.484632</td>\n",
       "      <td>-148.79...</td>\n",
       "      <td>494</td>\n",
       "      <td>456</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4521</th>\n",
       "      <td>494,457</td>\n",
       "      <td>-1120312.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.499342</td>\n",
       "      <td>-148.86...</td>\n",
       "      <td>494</td>\n",
       "      <td>457</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4522</th>\n",
       "      <td>494,458</td>\n",
       "      <td>-1117187.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.514021</td>\n",
       "      <td>-148.93...</td>\n",
       "      <td>494</td>\n",
       "      <td>458</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4523</th>\n",
       "      <td>494,459</td>\n",
       "      <td>-1114062.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.528670</td>\n",
       "      <td>-149.00...</td>\n",
       "      <td>494</td>\n",
       "      <td>459</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4524</th>\n",
       "      <td>494,460</td>\n",
       "      <td>-1110937.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.543288</td>\n",
       "      <td>-149.07...</td>\n",
       "      <td>494</td>\n",
       "      <td>460</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4525</th>\n",
       "      <td>494,461</td>\n",
       "      <td>-1107812.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.557875</td>\n",
       "      <td>-149.14...</td>\n",
       "      <td>494</td>\n",
       "      <td>461</td>\n",
       "      <td>132</td>\n",
       "      <td>89</td>\n",
       "      <td>110.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4526</th>\n",
       "      <td>494,462</td>\n",
       "      <td>-1104687.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.572432</td>\n",
       "      <td>-149.22...</td>\n",
       "      <td>494</td>\n",
       "      <td>462</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4527</th>\n",
       "      <td>494,463</td>\n",
       "      <td>-1101562.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.586959</td>\n",
       "      <td>-149.29...</td>\n",
       "      <td>494</td>\n",
       "      <td>463</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4528</th>\n",
       "      <td>494,464</td>\n",
       "      <td>-1098437.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.601454</td>\n",
       "      <td>-149.36...</td>\n",
       "      <td>494</td>\n",
       "      <td>464</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4529</th>\n",
       "      <td>494,465</td>\n",
       "      <td>-1095312.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.615919</td>\n",
       "      <td>-149.43...</td>\n",
       "      <td>494</td>\n",
       "      <td>465</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4530</th>\n",
       "      <td>494,466</td>\n",
       "      <td>-1092187.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.630352</td>\n",
       "      <td>-149.50...</td>\n",
       "      <td>494</td>\n",
       "      <td>466</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4531</th>\n",
       "      <td>494,467</td>\n",
       "      <td>-1089062.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.644755</td>\n",
       "      <td>-149.57...</td>\n",
       "      <td>494</td>\n",
       "      <td>467</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4532</th>\n",
       "      <td>494,468</td>\n",
       "      <td>-1085937.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.659127</td>\n",
       "      <td>-149.65...</td>\n",
       "      <td>494</td>\n",
       "      <td>468</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4533</th>\n",
       "      <td>494,469</td>\n",
       "      <td>-1082812.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.673468</td>\n",
       "      <td>-149.72...</td>\n",
       "      <td>494</td>\n",
       "      <td>469</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4534</th>\n",
       "      <td>494,470</td>\n",
       "      <td>-1079687.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.687777</td>\n",
       "      <td>-149.79...</td>\n",
       "      <td>494</td>\n",
       "      <td>470</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4535</th>\n",
       "      <td>494,471</td>\n",
       "      <td>-1076562.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.702055</td>\n",
       "      <td>-149.86...</td>\n",
       "      <td>494</td>\n",
       "      <td>471</td>\n",
       "      <td>132</td>\n",
       "      <td>90</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4536</th>\n",
       "      <td>494,472</td>\n",
       "      <td>-1073437.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.716302</td>\n",
       "      <td>-149.93...</td>\n",
       "      <td>494</td>\n",
       "      <td>472</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4537</th>\n",
       "      <td>494,473</td>\n",
       "      <td>-1070312.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.730518</td>\n",
       "      <td>-150.01...</td>\n",
       "      <td>494</td>\n",
       "      <td>473</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4538</th>\n",
       "      <td>494,474</td>\n",
       "      <td>-1067187.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.744702</td>\n",
       "      <td>-150.08...</td>\n",
       "      <td>494</td>\n",
       "      <td>474</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4539</th>\n",
       "      <td>494,475</td>\n",
       "      <td>-1064062.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.758855</td>\n",
       "      <td>-150.15...</td>\n",
       "      <td>494</td>\n",
       "      <td>475</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4540</th>\n",
       "      <td>494,476</td>\n",
       "      <td>-1060937.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.772976</td>\n",
       "      <td>-150.22...</td>\n",
       "      <td>494</td>\n",
       "      <td>476</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4541</th>\n",
       "      <td>494,477</td>\n",
       "      <td>-1057812.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.787066</td>\n",
       "      <td>-150.30...</td>\n",
       "      <td>494</td>\n",
       "      <td>477</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4542</th>\n",
       "      <td>494,478</td>\n",
       "      <td>-1054687.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.801124</td>\n",
       "      <td>-150.37...</td>\n",
       "      <td>494</td>\n",
       "      <td>478</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4543</th>\n",
       "      <td>494,479</td>\n",
       "      <td>-1051562.5</td>\n",
       "      <td>1854687.5</td>\n",
       "      <td>70.815150</td>\n",
       "      <td>-150.44...</td>\n",
       "      <td>494</td>\n",
       "      <td>479</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4544</th>\n",
       "      <td>495,416</td>\n",
       "      <td>-1248437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.895671</td>\n",
       "      <td>-146.00...</td>\n",
       "      <td>495</td>\n",
       "      <td>416</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4545</th>\n",
       "      <td>495,417</td>\n",
       "      <td>-1245312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.911561</td>\n",
       "      <td>-146.07...</td>\n",
       "      <td>495</td>\n",
       "      <td>417</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4546</th>\n",
       "      <td>495,418</td>\n",
       "      <td>-1242187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.927423</td>\n",
       "      <td>-146.14...</td>\n",
       "      <td>495</td>\n",
       "      <td>418</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4547</th>\n",
       "      <td>495,419</td>\n",
       "      <td>-1239062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.943258</td>\n",
       "      <td>-146.20...</td>\n",
       "      <td>495</td>\n",
       "      <td>419</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4548</th>\n",
       "      <td>495,420</td>\n",
       "      <td>-1235937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.959064</td>\n",
       "      <td>-146.27...</td>\n",
       "      <td>495</td>\n",
       "      <td>420</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4549</th>\n",
       "      <td>495,421</td>\n",
       "      <td>-1232812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.974842</td>\n",
       "      <td>-146.34...</td>\n",
       "      <td>495</td>\n",
       "      <td>421</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4550</th>\n",
       "      <td>495,422</td>\n",
       "      <td>-1229687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>69.990593</td>\n",
       "      <td>-146.41...</td>\n",
       "      <td>495</td>\n",
       "      <td>422</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4551</th>\n",
       "      <td>495,423</td>\n",
       "      <td>-1226562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.006315</td>\n",
       "      <td>-146.47...</td>\n",
       "      <td>495</td>\n",
       "      <td>423</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4552</th>\n",
       "      <td>495,424</td>\n",
       "      <td>-1223437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.022009</td>\n",
       "      <td>-146.54...</td>\n",
       "      <td>495</td>\n",
       "      <td>424</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4553</th>\n",
       "      <td>495,425</td>\n",
       "      <td>-1220312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.037674</td>\n",
       "      <td>-146.61...</td>\n",
       "      <td>495</td>\n",
       "      <td>425</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4554</th>\n",
       "      <td>495,426</td>\n",
       "      <td>-1217187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.053312</td>\n",
       "      <td>-146.67...</td>\n",
       "      <td>495</td>\n",
       "      <td>426</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4555</th>\n",
       "      <td>495,427</td>\n",
       "      <td>-1214062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.068920</td>\n",
       "      <td>-146.74...</td>\n",
       "      <td>495</td>\n",
       "      <td>427</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4556</th>\n",
       "      <td>495,428</td>\n",
       "      <td>-1210937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.084501</td>\n",
       "      <td>-146.81...</td>\n",
       "      <td>495</td>\n",
       "      <td>428</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4557</th>\n",
       "      <td>495,429</td>\n",
       "      <td>-1207812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.100053</td>\n",
       "      <td>-146.88...</td>\n",
       "      <td>495</td>\n",
       "      <td>429</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4558</th>\n",
       "      <td>495,430</td>\n",
       "      <td>-1204687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.115576</td>\n",
       "      <td>-146.95...</td>\n",
       "      <td>495</td>\n",
       "      <td>430</td>\n",
       "      <td>120</td>\n",
       "      <td>89</td>\n",
       "      <td>104.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4559</th>\n",
       "      <td>495,431</td>\n",
       "      <td>-1201562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.131071</td>\n",
       "      <td>-147.01...</td>\n",
       "      <td>495</td>\n",
       "      <td>431</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4560</th>\n",
       "      <td>495,432</td>\n",
       "      <td>-1198437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.146536</td>\n",
       "      <td>-147.08...</td>\n",
       "      <td>495</td>\n",
       "      <td>432</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4561</th>\n",
       "      <td>495,433</td>\n",
       "      <td>-1195312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.161973</td>\n",
       "      <td>-147.15...</td>\n",
       "      <td>495</td>\n",
       "      <td>433</td>\n",
       "      <td>115</td>\n",
       "      <td>65</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4562</th>\n",
       "      <td>495,434</td>\n",
       "      <td>-1192187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.177382</td>\n",
       "      <td>-147.22...</td>\n",
       "      <td>495</td>\n",
       "      <td>434</td>\n",
       "      <td>50</td>\n",
       "      <td>65</td>\n",
       "      <td>57.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4563</th>\n",
       "      <td>495,435</td>\n",
       "      <td>-1189062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.192761</td>\n",
       "      <td>-147.29...</td>\n",
       "      <td>495</td>\n",
       "      <td>435</td>\n",
       "      <td>50</td>\n",
       "      <td>65</td>\n",
       "      <td>57.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4564</th>\n",
       "      <td>495,436</td>\n",
       "      <td>-1185937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.208111</td>\n",
       "      <td>-147.36...</td>\n",
       "      <td>495</td>\n",
       "      <td>436</td>\n",
       "      <td>50</td>\n",
       "      <td>83</td>\n",
       "      <td>66.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4565</th>\n",
       "      <td>495,437</td>\n",
       "      <td>-1182812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.223433</td>\n",
       "      <td>-147.42...</td>\n",
       "      <td>495</td>\n",
       "      <td>437</td>\n",
       "      <td>65</td>\n",
       "      <td>89</td>\n",
       "      <td>77.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4566</th>\n",
       "      <td>495,438</td>\n",
       "      <td>-1179687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.238725</td>\n",
       "      <td>-147.49...</td>\n",
       "      <td>495</td>\n",
       "      <td>438</td>\n",
       "      <td>50</td>\n",
       "      <td>83</td>\n",
       "      <td>66.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4567</th>\n",
       "      <td>495,439</td>\n",
       "      <td>-1176562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.253988</td>\n",
       "      <td>-147.56...</td>\n",
       "      <td>495</td>\n",
       "      <td>439</td>\n",
       "      <td>115</td>\n",
       "      <td>89</td>\n",
       "      <td>102.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4568</th>\n",
       "      <td>495,440</td>\n",
       "      <td>-1173437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.269222</td>\n",
       "      <td>-147.63...</td>\n",
       "      <td>495</td>\n",
       "      <td>440</td>\n",
       "      <td>116</td>\n",
       "      <td>89</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4569</th>\n",
       "      <td>495,441</td>\n",
       "      <td>-1170312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.284426</td>\n",
       "      <td>-147.70...</td>\n",
       "      <td>495</td>\n",
       "      <td>441</td>\n",
       "      <td>116</td>\n",
       "      <td>89</td>\n",
       "      <td>102.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4570</th>\n",
       "      <td>495,442</td>\n",
       "      <td>-1167187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.299601</td>\n",
       "      <td>-147.77...</td>\n",
       "      <td>495</td>\n",
       "      <td>442</td>\n",
       "      <td>116</td>\n",
       "      <td>120</td>\n",
       "      <td>118.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4571</th>\n",
       "      <td>495,443</td>\n",
       "      <td>-1164062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.314747</td>\n",
       "      <td>-147.84...</td>\n",
       "      <td>495</td>\n",
       "      <td>443</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4572</th>\n",
       "      <td>495,444</td>\n",
       "      <td>-1160937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.329863</td>\n",
       "      <td>-147.91...</td>\n",
       "      <td>495</td>\n",
       "      <td>444</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4573</th>\n",
       "      <td>495,445</td>\n",
       "      <td>-1157812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.344950</td>\n",
       "      <td>-147.98...</td>\n",
       "      <td>495</td>\n",
       "      <td>445</td>\n",
       "      <td>120</td>\n",
       "      <td>120</td>\n",
       "      <td>120.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4574</th>\n",
       "      <td>495,446</td>\n",
       "      <td>-1154687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.360007</td>\n",
       "      <td>-148.05...</td>\n",
       "      <td>495</td>\n",
       "      <td>446</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4575</th>\n",
       "      <td>495,447</td>\n",
       "      <td>-1151562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.375035</td>\n",
       "      <td>-148.12...</td>\n",
       "      <td>495</td>\n",
       "      <td>447</td>\n",
       "      <td>134</td>\n",
       "      <td>89</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4576</th>\n",
       "      <td>495,448</td>\n",
       "      <td>-1148437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.390032</td>\n",
       "      <td>-148.19...</td>\n",
       "      <td>495</td>\n",
       "      <td>448</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4577</th>\n",
       "      <td>495,449</td>\n",
       "      <td>-1145312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.405000</td>\n",
       "      <td>-148.26...</td>\n",
       "      <td>495</td>\n",
       "      <td>449</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4578</th>\n",
       "      <td>495,450</td>\n",
       "      <td>-1142187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.419938</td>\n",
       "      <td>-148.33...</td>\n",
       "      <td>495</td>\n",
       "      <td>450</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4579</th>\n",
       "      <td>495,451</td>\n",
       "      <td>-1139062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.434846</td>\n",
       "      <td>-148.40...</td>\n",
       "      <td>495</td>\n",
       "      <td>451</td>\n",
       "      <td>134</td>\n",
       "      <td>89</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4580</th>\n",
       "      <td>495,452</td>\n",
       "      <td>-1135937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.449724</td>\n",
       "      <td>-148.47...</td>\n",
       "      <td>495</td>\n",
       "      <td>452</td>\n",
       "      <td>133</td>\n",
       "      <td>89</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4581</th>\n",
       "      <td>495,453</td>\n",
       "      <td>-1132812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.464572</td>\n",
       "      <td>-148.54...</td>\n",
       "      <td>495</td>\n",
       "      <td>453</td>\n",
       "      <td>133</td>\n",
       "      <td>89</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4582</th>\n",
       "      <td>495,454</td>\n",
       "      <td>-1129687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.479390</td>\n",
       "      <td>-148.61...</td>\n",
       "      <td>495</td>\n",
       "      <td>454</td>\n",
       "      <td>134</td>\n",
       "      <td>120</td>\n",
       "      <td>127.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4583</th>\n",
       "      <td>495,455</td>\n",
       "      <td>-1126562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.494178</td>\n",
       "      <td>-148.68...</td>\n",
       "      <td>495</td>\n",
       "      <td>455</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4584</th>\n",
       "      <td>495,456</td>\n",
       "      <td>-1123437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.508935</td>\n",
       "      <td>-148.75...</td>\n",
       "      <td>495</td>\n",
       "      <td>456</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4585</th>\n",
       "      <td>495,457</td>\n",
       "      <td>-1120312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.523662</td>\n",
       "      <td>-148.82...</td>\n",
       "      <td>495</td>\n",
       "      <td>457</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4586</th>\n",
       "      <td>495,458</td>\n",
       "      <td>-1117187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.538359</td>\n",
       "      <td>-148.89...</td>\n",
       "      <td>495</td>\n",
       "      <td>458</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4587</th>\n",
       "      <td>495,459</td>\n",
       "      <td>-1114062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.553025</td>\n",
       "      <td>-148.96...</td>\n",
       "      <td>495</td>\n",
       "      <td>459</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4588</th>\n",
       "      <td>495,460</td>\n",
       "      <td>-1110937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.567661</td>\n",
       "      <td>-149.03...</td>\n",
       "      <td>495</td>\n",
       "      <td>460</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4589</th>\n",
       "      <td>495,461</td>\n",
       "      <td>-1107812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.582266</td>\n",
       "      <td>-149.10...</td>\n",
       "      <td>495</td>\n",
       "      <td>461</td>\n",
       "      <td>132</td>\n",
       "      <td>120</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4590</th>\n",
       "      <td>495,462</td>\n",
       "      <td>-1104687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.596841</td>\n",
       "      <td>-149.17...</td>\n",
       "      <td>495</td>\n",
       "      <td>462</td>\n",
       "      <td>132</td>\n",
       "      <td>90</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4591</th>\n",
       "      <td>495,463</td>\n",
       "      <td>-1101562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.611384</td>\n",
       "      <td>-149.25...</td>\n",
       "      <td>495</td>\n",
       "      <td>463</td>\n",
       "      <td>132</td>\n",
       "      <td>90</td>\n",
       "      <td>111.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4592</th>\n",
       "      <td>495,464</td>\n",
       "      <td>-1098437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.625897</td>\n",
       "      <td>-149.32...</td>\n",
       "      <td>495</td>\n",
       "      <td>464</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4593</th>\n",
       "      <td>495,465</td>\n",
       "      <td>-1095312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.640380</td>\n",
       "      <td>-149.39...</td>\n",
       "      <td>495</td>\n",
       "      <td>465</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4594</th>\n",
       "      <td>495,466</td>\n",
       "      <td>-1092187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.654831</td>\n",
       "      <td>-149.46...</td>\n",
       "      <td>495</td>\n",
       "      <td>466</td>\n",
       "      <td>132</td>\n",
       "      <td>83</td>\n",
       "      <td>107.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4595</th>\n",
       "      <td>495,467</td>\n",
       "      <td>-1089062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.669251</td>\n",
       "      <td>-149.53...</td>\n",
       "      <td>495</td>\n",
       "      <td>467</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4596</th>\n",
       "      <td>495,468</td>\n",
       "      <td>-1085937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.683640</td>\n",
       "      <td>-149.60...</td>\n",
       "      <td>495</td>\n",
       "      <td>468</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4597</th>\n",
       "      <td>495,469</td>\n",
       "      <td>-1082812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.697999</td>\n",
       "      <td>-149.68...</td>\n",
       "      <td>495</td>\n",
       "      <td>469</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4598</th>\n",
       "      <td>495,470</td>\n",
       "      <td>-1079687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.712326</td>\n",
       "      <td>-149.75...</td>\n",
       "      <td>495</td>\n",
       "      <td>470</td>\n",
       "      <td>132</td>\n",
       "      <td>119</td>\n",
       "      <td>125.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4599</th>\n",
       "      <td>495,471</td>\n",
       "      <td>-1076562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.726621</td>\n",
       "      <td>-149.82...</td>\n",
       "      <td>495</td>\n",
       "      <td>471</td>\n",
       "      <td>133</td>\n",
       "      <td>119</td>\n",
       "      <td>126.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4600</th>\n",
       "      <td>495,472</td>\n",
       "      <td>-1073437.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.740886</td>\n",
       "      <td>-149.89...</td>\n",
       "      <td>495</td>\n",
       "      <td>472</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4601</th>\n",
       "      <td>495,473</td>\n",
       "      <td>-1070312.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.755119</td>\n",
       "      <td>-149.96...</td>\n",
       "      <td>495</td>\n",
       "      <td>473</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4602</th>\n",
       "      <td>495,474</td>\n",
       "      <td>-1067187.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.769320</td>\n",
       "      <td>-150.04...</td>\n",
       "      <td>495</td>\n",
       "      <td>474</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4603</th>\n",
       "      <td>495,475</td>\n",
       "      <td>-1064062.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.783491</td>\n",
       "      <td>-150.11...</td>\n",
       "      <td>495</td>\n",
       "      <td>475</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4604</th>\n",
       "      <td>495,476</td>\n",
       "      <td>-1060937.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.797629</td>\n",
       "      <td>-150.18...</td>\n",
       "      <td>495</td>\n",
       "      <td>476</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4605</th>\n",
       "      <td>495,477</td>\n",
       "      <td>-1057812.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.811736</td>\n",
       "      <td>-150.26...</td>\n",
       "      <td>495</td>\n",
       "      <td>477</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4606</th>\n",
       "      <td>495,478</td>\n",
       "      <td>-1054687.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.825811</td>\n",
       "      <td>-150.33...</td>\n",
       "      <td>495</td>\n",
       "      <td>478</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4607</th>\n",
       "      <td>495,479</td>\n",
       "      <td>-1051562.5</td>\n",
       "      <td>1851562.5</td>\n",
       "      <td>70.839855</td>\n",
       "      <td>-150.40...</td>\n",
       "      <td>495</td>\n",
       "      <td>479</td>\n",
       "      <td>133</td>\n",
       "      <td>90</td>\n",
       "      <td>111.5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>4608 rows × 10 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        pixel          x          y   latitude  longitude  row  column  2018  2019    Avg\n",
       "0     424,416 -1248437.5  2073437.5  68.189531 -148.94...  424     416   129   121  125.0\n",
       "1     424,417 -1245312.5  2073437.5  68.204235 -149.01...  424     417   129   121  125.0\n",
       "2     424,418 -1242187.5  2073437.5  68.218912 -149.07...  424     418   129   121  125.0\n",
       "3     424,419 -1239062.5  2073437.5  68.233561 -149.13...  424     419   129   121  125.0\n",
       "4     424,420 -1235937.5  2073437.5  68.248182 -149.20...  424     420   129   121  125.0\n",
       "5     424,421 -1232812.5  2073437.5  68.262776 -149.26...  424     421   129   121  125.0\n",
       "6     424,422 -1229687.5  2073437.5  68.277342 -149.32...  424     422   129   121  125.0\n",
       "7     424,423 -1226562.5  2073437.5  68.291880 -149.39...  424     423   129   120  124.5\n",
       "8     424,424 -1223437.5  2073437.5  68.306391 -149.45...  424     424   129   120  124.5\n",
       "9     424,425 -1220312.5  2073437.5  68.320873 -149.52...  424     425   129   120  124.5\n",
       "10    424,426 -1217187.5  2073437.5  68.335328 -149.58...  424     426   129   120  124.5\n",
       "11    424,427 -1214062.5  2073437.5  68.349755 -149.64...  424     427   129   120  124.5\n",
       "12    424,428 -1210937.5  2073437.5  68.364154 -149.71...  424     428   129   119  124.0\n",
       "13    424,429 -1207812.5  2073437.5  68.378525 -149.77...  424     429   129   119  124.0\n",
       "14    424,430 -1204687.5  2073437.5  68.392867 -149.84...  424     430   129   119  124.0\n",
       "15    424,431 -1201562.5  2073437.5  68.407182 -149.90...  424     431   129   119  124.0\n",
       "16    424,432 -1198437.5  2073437.5  68.421468 -149.97...  424     432   129    90  109.5\n",
       "17    424,433 -1195312.5  2073437.5  68.435726 -150.03...  424     433   129    90  109.5\n",
       "18    424,434 -1192187.5  2073437.5  68.449956 -150.10...  424     434   129    90  109.5\n",
       "19    424,435 -1189062.5  2073437.5  68.464157 -150.16...  424     435   129    90  109.5\n",
       "20    424,436 -1185937.5  2073437.5  68.478330 -150.23...  424     436   129    90  109.5\n",
       "21    424,437 -1182812.5  2073437.5  68.492475 -150.29...  424     437   128    90  109.0\n",
       "22    424,438 -1179687.5  2073437.5  68.506591 -150.36...  424     438   128    90  109.0\n",
       "23    424,439 -1176562.5  2073437.5  68.520678 -150.42...  424     439   128    90  109.0\n",
       "24    424,440 -1173437.5  2073437.5  68.534737 -150.49...  424     440   128    90  109.0\n",
       "25    424,441 -1170312.5  2073437.5  68.548767 -150.55...  424     441   128    90  109.0\n",
       "26    424,442 -1167187.5  2073437.5  68.562769 -150.62...  424     442   128    90  109.0\n",
       "27    424,443 -1164062.5  2073437.5  68.576741 -150.68...  424     443   128    90  109.0\n",
       "28    424,444 -1160937.5  2073437.5  68.590685 -150.75...  424     444   128    90  109.0\n",
       "29    424,445 -1157812.5  2073437.5  68.604600 -150.82...  424     445   128    90  109.0\n",
       "30    424,446 -1154687.5  2073437.5  68.618486 -150.88...  424     446   117    90  103.5\n",
       "31    424,447 -1151562.5  2073437.5  68.632343 -150.95...  424     447   117    90  103.5\n",
       "32    424,448 -1148437.5  2073437.5  68.646170 -151.01...  424     448   117    90  103.5\n",
       "33    424,449 -1145312.5  2073437.5  68.659969 -151.08...  424     449   115    90  102.5\n",
       "34    424,450 -1142187.5  2073437.5  68.673739 -151.15...  424     450   115    90  102.5\n",
       "35    424,451 -1139062.5  2073437.5  68.687479 -151.21...  424     451   115    90  102.5\n",
       "36    424,452 -1135937.5  2073437.5  68.701191 -151.28...  424     452   115    90  102.5\n",
       "37    424,453 -1132812.5  2073437.5  68.714872 -151.35...  424     453   115    90  102.5\n",
       "38    424,454 -1129687.5  2073437.5  68.728525 -151.41...  424     454   115    90  102.5\n",
       "39    424,455 -1126562.5  2073437.5  68.742148 -151.48...  424     455   115    90  102.5\n",
       "40    424,456 -1123437.5  2073437.5  68.755742 -151.55...  424     456   115    90  102.5\n",
       "41    424,457 -1120312.5  2073437.5  68.769306 -151.61...  424     457   115    90  102.5\n",
       "42    424,458 -1117187.5  2073437.5  68.782840 -151.68...  424     458   115    90  102.5\n",
       "43    424,459 -1114062.5  2073437.5  68.796345 -151.75...  424     459   115    90  102.5\n",
       "44    424,460 -1110937.5  2073437.5  68.809820 -151.81...  424     460   115    90  102.5\n",
       "45    424,461 -1107812.5  2073437.5  68.823266 -151.88...  424     461   115    90  102.5\n",
       "46    424,462 -1104687.5  2073437.5  68.836681 -151.95...  424     462   115    90  102.5\n",
       "47    424,463 -1101562.5  2073437.5  68.850067 -152.01...  424     463   115    90  102.5\n",
       "48    424,464 -1098437.5  2073437.5  68.863423 -152.08...  424     464   115    90  102.5\n",
       "49    424,465 -1095312.5  2073437.5  68.876749 -152.15...  424     465   115    90  102.5\n",
       "50    424,466 -1092187.5  2073437.5  68.890045 -152.22...  424     466   115    90  102.5\n",
       "51    424,467 -1089062.5  2073437.5  68.903311 -152.28...  424     467   115    90  102.5\n",
       "52    424,468 -1085937.5  2073437.5  68.916547 -152.35...  424     468   115    90  102.5\n",
       "53    424,469 -1082812.5  2073437.5  68.929753 -152.42...  424     469   115    90  102.5\n",
       "54    424,470 -1079687.5  2073437.5  68.942928 -152.49...  424     470   115    90  102.5\n",
       "55    424,471 -1076562.5  2073437.5  68.956073 -152.56...  424     471   115    90  102.5\n",
       "56    424,472 -1073437.5  2073437.5  68.969188 -152.62...  424     472   115    90  102.5\n",
       "57    424,473 -1070312.5  2073437.5  68.982273 -152.69...  424     473   115    90  102.5\n",
       "58    424,474 -1067187.5  2073437.5  68.995327 -152.76...  424     474   115    90  102.5\n",
       "59    424,475 -1064062.5  2073437.5  69.008350 -152.83...  424     475   115    90  102.5\n",
       "60    424,476 -1060937.5  2073437.5  69.021343 -152.90...  424     476   115    90  102.5\n",
       "61    424,477 -1057812.5  2073437.5  69.034306 -152.97...  424     477   115    90  102.5\n",
       "62    424,478 -1054687.5  2073437.5  69.047238 -153.03...  424     478   115   119  117.0\n",
       "63    424,479 -1051562.5  2073437.5  69.060139 -153.10...  424     479   116   119  117.5\n",
       "64    425,416 -1248437.5  2070312.5  68.213969 -148.90...  425     416   129   121  125.0\n",
       "65    425,417 -1245312.5  2070312.5  68.228689 -148.97...  425     417   129   121  125.0\n",
       "66    425,418 -1242187.5  2070312.5  68.243381 -149.03...  425     418   129   121  125.0\n",
       "67    425,419 -1239062.5  2070312.5  68.258046 -149.09...  425     419   129   121  125.0\n",
       "68    425,420 -1235937.5  2070312.5  68.272683 -149.16...  425     420   129   121  125.0\n",
       "69    425,421 -1232812.5  2070312.5  68.287292 -149.22...  425     421   129   120  124.5\n",
       "70    425,422 -1229687.5  2070312.5  68.301874 -149.29...  425     422   129   120  124.5\n",
       "71    425,423 -1226562.5  2070312.5  68.316428 -149.35...  425     423   129   120  124.5\n",
       "72    425,424 -1223437.5  2070312.5  68.330954 -149.41...  425     424   129   120  124.5\n",
       "73    425,425 -1220312.5  2070312.5  68.345452 -149.48...  425     425   129   120  124.5\n",
       "74    425,426 -1217187.5  2070312.5  68.359922 -149.54...  425     426   129   119  124.0\n",
       "75    425,427 -1214062.5  2070312.5  68.374365 -149.61...  425     427   129   119  124.0\n",
       "76    425,428 -1210937.5  2070312.5  68.388779 -149.67...  425     428   129   119  124.0\n",
       "77    425,429 -1207812.5  2070312.5  68.403165 -149.74...  425     429   129   119  124.0\n",
       "78    425,430 -1204687.5  2070312.5  68.417524 -149.80...  425     430   129   119  124.0\n",
       "79    425,431 -1201562.5  2070312.5  68.431854 -149.87...  425     431   129    90  109.5\n",
       "80    425,432 -1198437.5  2070312.5  68.446156 -149.93...  425     432   129    90  109.5\n",
       "81    425,433 -1195312.5  2070312.5  68.460429 -149.99...  425     433   129    90  109.5\n",
       "82    425,434 -1192187.5  2070312.5  68.474674 -150.06...  425     434   129    90  109.5\n",
       "83    425,435 -1189062.5  2070312.5  68.488891 -150.12...  425     435   128    90  109.0\n",
       "84    425,436 -1185937.5  2070312.5  68.503080 -150.19...  425     436   128    90  109.0\n",
       "85    425,437 -1182812.5  2070312.5  68.517240 -150.25...  425     437   128    90  109.0\n",
       "86    425,438 -1179687.5  2070312.5  68.531371 -150.32...  425     438   128    90  109.0\n",
       "87    425,439 -1176562.5  2070312.5  68.545474 -150.39...  425     439   128    90  109.0\n",
       "88    425,440 -1173437.5  2070312.5  68.559548 -150.45...  425     440   128    90  109.0\n",
       "89    425,441 -1170312.5  2070312.5  68.573594 -150.52...  425     441   128    90  109.0\n",
       "90    425,442 -1167187.5  2070312.5  68.587611 -150.58...  425     442   128    90  109.0\n",
       "91    425,443 -1164062.5  2070312.5  68.601599 -150.65...  425     443   128    90  109.0\n",
       "92    425,444 -1160937.5  2070312.5  68.615558 -150.71...  425     444   128    90  109.0\n",
       "93    425,445 -1157812.5  2070312.5  68.629488 -150.78...  425     445   117    90  103.5\n",
       "94    425,446 -1154687.5  2070312.5  68.643389 -150.85...  425     446   117    90  103.5\n",
       "95    425,447 -1151562.5  2070312.5  68.657261 -150.91...  425     447   115    90  102.5\n",
       "96    425,448 -1148437.5  2070312.5  68.671105 -150.98...  425     448   115    90  102.5\n",
       "97    425,449 -1145312.5  2070312.5  68.684919 -151.04...  425     449   115    90  102.5\n",
       "98    425,450 -1142187.5  2070312.5  68.698704 -151.11...  425     450   115    90  102.5\n",
       "99    425,451 -1139062.5  2070312.5  68.712460 -151.18...  425     451   115    90  102.5\n",
       "...       ...        ...        ...        ...        ...  ...     ...   ...   ...    ...\n",
       "4508  494,444 -1160937.5  1854687.5  70.305772 -147.95...  494     444   116    89  102.5\n",
       "4509  494,445 -1157812.5  1854687.5  70.320841 -148.02...  494     445   116    89  102.5\n",
       "4510  494,446 -1154687.5  1854687.5  70.335881 -148.09...  494     446   116    89  102.5\n",
       "4511  494,447 -1151562.5  1854687.5  70.350890 -148.16...  494     447   120    89  104.5\n",
       "4512  494,448 -1148437.5  1854687.5  70.365870 -148.23...  494     448   120    89  104.5\n",
       "4513  494,449 -1145312.5  1854687.5  70.380820 -148.30...  494     449   134    89  111.5\n",
       "4514  494,450 -1142187.5  1854687.5  70.395741 -148.37...  494     450   134    89  111.5\n",
       "4515  494,451 -1139062.5  1854687.5  70.410631 -148.44...  494     451   133    89  111.0\n",
       "4516  494,452 -1135937.5  1854687.5  70.425492 -148.51...  494     452   133    89  111.0\n",
       "4517  494,453 -1132812.5  1854687.5  70.440322 -148.58...  494     453   133    89  111.0\n",
       "4518  494,454 -1129687.5  1854687.5  70.455122 -148.65...  494     454   132    89  110.5\n",
       "4519  494,455 -1126562.5  1854687.5  70.469892 -148.72...  494     455   132    89  110.5\n",
       "4520  494,456 -1123437.5  1854687.5  70.484632 -148.79...  494     456   132    89  110.5\n",
       "4521  494,457 -1120312.5  1854687.5  70.499342 -148.86...  494     457   132    89  110.5\n",
       "4522  494,458 -1117187.5  1854687.5  70.514021 -148.93...  494     458   132   120  126.0\n",
       "4523  494,459 -1114062.5  1854687.5  70.528670 -149.00...  494     459   132   120  126.0\n",
       "4524  494,460 -1110937.5  1854687.5  70.543288 -149.07...  494     460   132    89  110.5\n",
       "4525  494,461 -1107812.5  1854687.5  70.557875 -149.14...  494     461   132    89  110.5\n",
       "4526  494,462 -1104687.5  1854687.5  70.572432 -149.22...  494     462   132    83  107.5\n",
       "4527  494,463 -1101562.5  1854687.5  70.586959 -149.29...  494     463   132    83  107.5\n",
       "4528  494,464 -1098437.5  1854687.5  70.601454 -149.36...  494     464   132    83  107.5\n",
       "4529  494,465 -1095312.5  1854687.5  70.615919 -149.43...  494     465   132    83  107.5\n",
       "4530  494,466 -1092187.5  1854687.5  70.630352 -149.50...  494     466   132    83  107.5\n",
       "4531  494,467 -1089062.5  1854687.5  70.644755 -149.57...  494     467   132    83  107.5\n",
       "4532  494,468 -1085937.5  1854687.5  70.659127 -149.65...  494     468   132    83  107.5\n",
       "4533  494,469 -1082812.5  1854687.5  70.673468 -149.72...  494     469   132   119  125.5\n",
       "4534  494,470 -1079687.5  1854687.5  70.687777 -149.79...  494     470   132   119  125.5\n",
       "4535  494,471 -1076562.5  1854687.5  70.702055 -149.86...  494     471   132    90  111.0\n",
       "4536  494,472 -1073437.5  1854687.5  70.716302 -149.93...  494     472   132   119  125.5\n",
       "4537  494,473 -1070312.5  1854687.5  70.730518 -150.01...  494     473   133    90  111.5\n",
       "4538  494,474 -1067187.5  1854687.5  70.744702 -150.08...  494     474   133    90  111.5\n",
       "4539  494,475 -1064062.5  1854687.5  70.758855 -150.15...  494     475   133    90  111.5\n",
       "4540  494,476 -1060937.5  1854687.5  70.772976 -150.22...  494     476   133    90  111.5\n",
       "4541  494,477 -1057812.5  1854687.5  70.787066 -150.30...  494     477   133    90  111.5\n",
       "4542  494,478 -1054687.5  1854687.5  70.801124 -150.37...  494     478   133    90  111.5\n",
       "4543  494,479 -1051562.5  1854687.5  70.815150 -150.44...  494     479   133    90  111.5\n",
       "4544  495,416 -1248437.5  1851562.5  69.895671 -146.00...  495     416   115    89  102.0\n",
       "4545  495,417 -1245312.5  1851562.5  69.911561 -146.07...  495     417   115    89  102.0\n",
       "4546  495,418 -1242187.5  1851562.5  69.927423 -146.14...  495     418   115    89  102.0\n",
       "4547  495,419 -1239062.5  1851562.5  69.943258 -146.20...  495     419   115    89  102.0\n",
       "4548  495,420 -1235937.5  1851562.5  69.959064 -146.27...  495     420   120    89  104.5\n",
       "4549  495,421 -1232812.5  1851562.5  69.974842 -146.34...  495     421   120    89  104.5\n",
       "4550  495,422 -1229687.5  1851562.5  69.990593 -146.41...  495     422   120    89  104.5\n",
       "4551  495,423 -1226562.5  1851562.5  70.006315 -146.47...  495     423   120    89  104.5\n",
       "4552  495,424 -1223437.5  1851562.5  70.022009 -146.54...  495     424   120    89  104.5\n",
       "4553  495,425 -1220312.5  1851562.5  70.037674 -146.61...  495     425   120    89  104.5\n",
       "4554  495,426 -1217187.5  1851562.5  70.053312 -146.67...  495     426   120    89  104.5\n",
       "4555  495,427 -1214062.5  1851562.5  70.068920 -146.74...  495     427   120    89  104.5\n",
       "4556  495,428 -1210937.5  1851562.5  70.084501 -146.81...  495     428   120    89  104.5\n",
       "4557  495,429 -1207812.5  1851562.5  70.100053 -146.88...  495     429   120    89  104.5\n",
       "4558  495,430 -1204687.5  1851562.5  70.115576 -146.95...  495     430   120    89  104.5\n",
       "4559  495,431 -1201562.5  1851562.5  70.131071 -147.01...  495     431   115    89  102.0\n",
       "4560  495,432 -1198437.5  1851562.5  70.146536 -147.08...  495     432   115    89  102.0\n",
       "4561  495,433 -1195312.5  1851562.5  70.161973 -147.15...  495     433   115    65   90.0\n",
       "4562  495,434 -1192187.5  1851562.5  70.177382 -147.22...  495     434    50    65   57.5\n",
       "4563  495,435 -1189062.5  1851562.5  70.192761 -147.29...  495     435    50    65   57.5\n",
       "4564  495,436 -1185937.5  1851562.5  70.208111 -147.36...  495     436    50    83   66.5\n",
       "4565  495,437 -1182812.5  1851562.5  70.223433 -147.42...  495     437    65    89   77.0\n",
       "4566  495,438 -1179687.5  1851562.5  70.238725 -147.49...  495     438    50    83   66.5\n",
       "4567  495,439 -1176562.5  1851562.5  70.253988 -147.56...  495     439   115    89  102.0\n",
       "4568  495,440 -1173437.5  1851562.5  70.269222 -147.63...  495     440   116    89  102.5\n",
       "4569  495,441 -1170312.5  1851562.5  70.284426 -147.70...  495     441   116    89  102.5\n",
       "4570  495,442 -1167187.5  1851562.5  70.299601 -147.77...  495     442   116   120  118.0\n",
       "4571  495,443 -1164062.5  1851562.5  70.314747 -147.84...  495     443   134   120  127.0\n",
       "4572  495,444 -1160937.5  1851562.5  70.329863 -147.91...  495     444   134   120  127.0\n",
       "4573  495,445 -1157812.5  1851562.5  70.344950 -147.98...  495     445   120   120  120.0\n",
       "4574  495,446 -1154687.5  1851562.5  70.360007 -148.05...  495     446   134   120  127.0\n",
       "4575  495,447 -1151562.5  1851562.5  70.375035 -148.12...  495     447   134    89  111.5\n",
       "4576  495,448 -1148437.5  1851562.5  70.390032 -148.19...  495     448   134   120  127.0\n",
       "4577  495,449 -1145312.5  1851562.5  70.405000 -148.26...  495     449   134   120  127.0\n",
       "4578  495,450 -1142187.5  1851562.5  70.419938 -148.33...  495     450   134   120  127.0\n",
       "4579  495,451 -1139062.5  1851562.5  70.434846 -148.40...  495     451   134    89  111.5\n",
       "4580  495,452 -1135937.5  1851562.5  70.449724 -148.47...  495     452   133    89  111.0\n",
       "4581  495,453 -1132812.5  1851562.5  70.464572 -148.54...  495     453   133    89  111.0\n",
       "4582  495,454 -1129687.5  1851562.5  70.479390 -148.61...  495     454   134   120  127.0\n",
       "4583  495,455 -1126562.5  1851562.5  70.494178 -148.68...  495     455   132   120  126.0\n",
       "4584  495,456 -1123437.5  1851562.5  70.508935 -148.75...  495     456   132   120  126.0\n",
       "4585  495,457 -1120312.5  1851562.5  70.523662 -148.82...  495     457   132   120  126.0\n",
       "4586  495,458 -1117187.5  1851562.5  70.538359 -148.89...  495     458   132   120  126.0\n",
       "4587  495,459 -1114062.5  1851562.5  70.553025 -148.96...  495     459   132   120  126.0\n",
       "4588  495,460 -1110937.5  1851562.5  70.567661 -149.03...  495     460   132   120  126.0\n",
       "4589  495,461 -1107812.5  1851562.5  70.582266 -149.10...  495     461   132   120  126.0\n",
       "4590  495,462 -1104687.5  1851562.5  70.596841 -149.17...  495     462   132    90  111.0\n",
       "4591  495,463 -1101562.5  1851562.5  70.611384 -149.25...  495     463   132    90  111.0\n",
       "4592  495,464 -1098437.5  1851562.5  70.625897 -149.32...  495     464   132   119  125.5\n",
       "4593  495,465 -1095312.5  1851562.5  70.640380 -149.39...  495     465   132    83  107.5\n",
       "4594  495,466 -1092187.5  1851562.5  70.654831 -149.46...  495     466   132    83  107.5\n",
       "4595  495,467 -1089062.5  1851562.5  70.669251 -149.53...  495     467   132   119  125.5\n",
       "4596  495,468 -1085937.5  1851562.5  70.683640 -149.60...  495     468   132   119  125.5\n",
       "4597  495,469 -1082812.5  1851562.5  70.697999 -149.68...  495     469   132   119  125.5\n",
       "4598  495,470 -1079687.5  1851562.5  70.712326 -149.75...  495     470   132   119  125.5\n",
       "4599  495,471 -1076562.5  1851562.5  70.726621 -149.82...  495     471   133   119  126.0\n",
       "4600  495,472 -1073437.5  1851562.5  70.740886 -149.89...  495     472   133    90  111.5\n",
       "4601  495,473 -1070312.5  1851562.5  70.755119 -149.96...  495     473   133    90  111.5\n",
       "4602  495,474 -1067187.5  1851562.5  70.769320 -150.04...  495     474   133    90  111.5\n",
       "4603  495,475 -1064062.5  1851562.5  70.783491 -150.11...  495     475   133    90  111.5\n",
       "4604  495,476 -1060937.5  1851562.5  70.797629 -150.18...  495     476   133    90  111.5\n",
       "4605  495,477 -1057812.5  1851562.5  70.811736 -150.26...  495     477   133    90  111.5\n",
       "4606  495,478 -1054687.5  1851562.5  70.825811 -150.33...  495     478   133    90  111.5\n",
       "4607  495,479 -1051562.5  1851562.5  70.839855 -150.40...  495     479   133    90  111.5\n",
       "\n",
       "[4608 rows x 10 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# get array of average MODs for SIR pixels for all the years loaded\n",
    "#Colorado (Johnson et al 2020) used 5 times in 7 day window\n",
    "#Patagonia (Monahan and Ramage 2010) and Yukon (Semmens et al 2013?) used 3 times in 5 day window\n",
    "window = 1   # window for MOD algorithm, '10' would be 5 days (2 measurements per day)\n",
    "count = 1    # number of Tb/DAV exceedances to trigger MOD\n",
    "#From Johnson et al 2020 AMSRE rSIR Tb >= 249 DAV>=13 and AMSRE GRD Tb>=243 DAV>=14\n",
    "#From Johnson et al 2020 SSMI rSIR and GRD Tb>=247 DAV>=10\n",
    "DAV_threshold = 10\n",
    "Tb_threshold = 247\n",
    "\n",
    "# sir MOD array - MOD will be in day of year (DOY) #changeback to Years for all years\n",
    "MOD_DOY_array, MOD_DOY_df, meltflag_df = MOD_array(datadir_SIR, prefix_SIR, data_SIR, DAV_SIR, \n",
    "                                      rows_cols_env, subYears, window, count, \n",
    "                                      DAV_threshold, Tb_threshold)\n",
    "MOD_DOY_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'num_lats_SIR' 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/tmp/ipykernel_865226/1555237545.py\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0;31m#Make a matching data frame with lat long for each pixel\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0marray\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mzeros\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mnum_lats_SIR\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      3\u001b[0m \u001b[0mgeolocation_SIR_df\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpd\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mDataFrame\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0marray\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolumns\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmeltflag_df\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcolumns\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mindex\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'Latitude'\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m'Longitude'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      4\u001b[0m \u001b[0mgeolocation_SIR_df\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mloc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'Latitude'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mpixel_lats\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      5\u001b[0m \u001b[0mgeolocation_SIR_df\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mloc\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'Longitude'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mpixel_lons\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mNameError\u001b[0m: name 'num_lats_SIR' is not defined"
     ]
    }
   ],
   "source": [
    "#Make a matching data frame with lat long for each pixel\n",
    "array=np.zeros((2,num_lats_SIR))\n",
    "geolocation_SIR_df = pd.DataFrame(array, columns = meltflag_df.columns, index=['Latitude','Longitude'])\n",
    "geolocation_SIR_df.loc['Latitude']=pixel_lats.data\n",
    "geolocation_SIR_df.loc['Longitude']=pixel_lons.data\n",
    "geolocation_SIR_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# MOD of the GRD pixel - avg all years\n",
    "MOD_DOY_array_GRD, MOD_DOY_GRD_df, meltflag_GRD_df = MOD_array(datadir_GRD, prefix_GRD, CETB_GRD, DAV_GRD, \n",
    "                                              rows_cols_GRD, cal_date, subYears, window, count, \n",
    "                                              DAV_threshold, Tb_threshold)\n",
    "MOD_DOY_GRD_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "array=np.zeros((2,num_lats_GRD))\n",
    "geolocation_GRD_df = pd.DataFrame(array, columns = meltflag_GRD_df.columns, index=['Latitude','Longitude'])\n",
    "geolocation_GRD_df.loc['Latitude'] = pixel_lats_GRD.data\n",
    "geolocation_GRD_df.loc['Longitude'] = pixel_lons_GRD.data\n",
    "geolocation_GRD_df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Notes about changes from original notebooks:\n",
    "\n",
    "The old notebooks used to call MOD_array to get the average MOD for a set of years \n",
    "and then call MOD_array_year for a given year of interest.\n",
    "\n",
    "Now we just call MOD_array for SIR data and GRD data, and get back a data frame with \n",
    "MOD columns for each individual year, and one column for the avg MOD for all the years.\n",
    "\n",
    "This will run much faster, and can be saved and just re-read from a saved file on disk.\n",
    "\n",
    "N.B. Eventually we should just be able to read saved pickle files with the data we have already saved with:\n",
    "\n",
    "new = pd.read_pickle(MOD_DOY_filename)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Save the MOD by year data frames for SIR and GRD to pickle files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "dataDir"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "modsDir = \"%s/MODs\" % dataDir\n",
    "if not os.path.isdir(modsDir):\n",
    "    os.makedirs(modsDir)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.MOD.pkl\" % (\n",
    "    modsDir, region, areaname, sat_SIR, channel_SIR, alg_SIR)\n",
    "MOD_DOY_df.to_pickle(filename)\n",
    "print(\"MOD_DOY dataframe saved to %s\\n\" % filename)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.MOD.pkl\" % (\n",
    "    modsDir, region, areaname, sat_GRD, channel_GRD, alg_GRD)\n",
    "MOD_DOY_GRD_df.to_pickle(filename)\n",
    "print(\"MOD_DOY_GRD dataframe saved to %s\\n\" % filename)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.meltflag.pkl\" % (\n",
    "    modsDir, region, areaname, sat_SIR, channel_SIR, alg_SIR)\n",
    "meltflag_df.to_pickle(filename)\n",
    "print(\"meltflag_df dataframe saved to %s\\n\" % filename)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.meltflag.pkl\" % (\n",
    "    modsDir, region, areaname, sat_GRD, channel_GRD, alg_GRD)\n",
    "meltflag_GRD_df.to_pickle(filename)\n",
    "print(\"meltflag_GRD_df dataframe saved to %s\\n\" % filename)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.geolocation.pkl\" % (\n",
    "    modsDir, region, areaname, sat_SIR, channel_SIR, alg_SIR)\n",
    "geolocation_SIR_df.to_pickle(filename)\n",
    "print(\"geolocation_SIR_df dataframe saved to %s\\n\" % filename)\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.geolocation.pkl\" % (\n",
    "    modsDir, region, areaname, sat_GRD, channel_GRD, alg_GRD)\n",
    "geolocation_GRD_df.to_pickle(filename)\n",
    "print(\"geolocation_GRD_df dataframe saved to %s\\n\" % filename)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## create the maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#how to create subplots with basemap (w maps) https://basemaptutorial.readthedocs.io/en/latest/subplots.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "# PLOT the MOD MAP - avg of all years or one year of interest - SIR\n",
    "# Choose plot type - 'year' for one year of interest, 'all' for average of all years\n",
    "#plot_type='year'  \n",
    "#year=2010\n",
    "plot_type='all'\n",
    "\n",
    "\n",
    "if plot_type == 'all':\n",
    "    array = MOD_DOY_df['Avg']\n",
    "    label = 'AvgDOY'+str(subYears[0])+'-'+str(subYears[-1])\n",
    "    title = sensor_SIR + '-' + Site + ' - MOD - (Avg DOY '+str(subYears[0])+'-'+str(subYears[-1])+')'\n",
    "elif plot_type == 'year':\n",
    "    array = MOD_DOY_df[year]\n",
    "    label = str(year)+'DOY' \n",
    "    title = sensor_SIR + '-' + Site + ' - MOD - (' + str(year) + ' DOY)'\n",
    "\n",
    "# Set a few common things\n",
    "graticule_fontsize = 4\n",
    "\n",
    "# create figure and axes instances\n",
    "fig = plt.figure(figsize=(8,8))\n",
    "ax = fig.add_axes([0.1,0.1,0.8,0.8])\n",
    "#ax = fig.add_subplot(121)\n",
    "# create polar stereographic Basemap instance.\n",
    "m = Basemap(projection='stere',\n",
    "            lon_0=(lon_start+lon_end)/2,\n",
    "            lat_0=(lat_start+lat_end)/2,\n",
    "            lat_ts=90.,\n",
    "            llcrnrlat=(lat_start-.3),\n",
    "            urcrnrlat=(lat_end+.5),\n",
    "            llcrnrlon=(lon_start-.3),\n",
    "            urcrnrlon=(lon_end+.3),\n",
    "            rsphere=6371200.,\n",
    "            resolution='l',\n",
    "            area_thresh=10000, \n",
    "            epsg=3857)\n",
    "m.drawcoastlines()\n",
    "m.drawcountries()\n",
    "m.drawstates()\n",
    "#m.etopo(scale=5, alpha=0.7)\n",
    "#m.bluemarble(scale=1)\n",
    "parallels = np.arange(0.,90,0.5)\n",
    "m.drawparallels(parallels, labels=[1,0,0,0], fontsize=graticule_fontsize)\n",
    "# draw meridians\n",
    "meridians = np.arange(0,360.,0.5)\n",
    "m.drawmeridians(meridians, labels=[0,0,0,1], fontsize=graticule_fontsize)\n",
    "m.arcgisimage(service='World_Physical_Map', xpixels = 1500, verbose= True)\n",
    "\n",
    "#m.readshapefile('/home/mij216/ExploringCETB/IN_Hunza_at_DainyorBridge', 'Hunza', color='red', linewidth=1)\n",
    "#m.readshapefile('/mnt/data3/rgi6.0/06_rgi60_Iceland', 'Iceland', color='red', linewidth=1)\n",
    "#m.readshapefile('/mnt/data3/Vicki-Kuparuk/aga_kuparuk_basin_hydrology', 'aga_kuparuk_basin_hydrology', color='red', linewidth=1)\n",
    "\n",
    "x,y = m(pixel_lons, pixel_lats)\n",
    "#m.scatter(x,y,c=array, s=30, marker='s',lw=0,cmap='RdBu', alpha=.6) #this line is what was here\n",
    "im = m.scatter(x, y, c=array, s=30, marker='s',\n",
    "              lw=0, cmap='BuPu_r', alpha=.6) # this line is for testing PuOr, rainbow,binary good for Tb\n",
    "plt.title(title)\n",
    "\n",
    "# Get the current axes and configure placement so colorbar will\n",
    "# be 5% of ax width and the padding will be 0.05 in.\n",
    "ax = plt.gca()\n",
    "divider = make_axes_locatable(ax)\n",
    "cax = divider.append_axes(\"right\", size=\"5%\", pad=0.05)\n",
    "plt.colorbar(im, cax=cax, label='DOY')\n",
    "plt.clim(75,175)  #color bar limits\n",
    "\n",
    "plt.show()\n",
    "\n",
    "filename = \"%s/%s.%s.%s.%s.%s.%s.png\" % (\n",
    "    modsDir, region, areaname, sat_SIR, channel_SIR, alg_SIR, label)\n",
    "plt.savefig(filename)\n",
    "print(\"avg saved to %s\\n\" % filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#color bars here https://matplotlib.org/stable/tutorials/colors/colormaps.html\n",
    "#documentation on basemap https://matplotlib.org/basemap/users/laea.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Trial-and-error decide on size for GRD vs. SIR scatter markers\n",
    "#np.sqrt((1200)/8)^2\n",
    "sGRD=700\n",
    "sSIR=12"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 2plots left (GRD) and right (SIR)\n",
    "#redo basemap to be in easegrid projection EASE2.0\n",
    "\n",
    "# Choose plot type - 'year' for one year of interest, 'all' for average of all years\n",
    "plot_type='year'  \n",
    "year=2018\n",
    "#plot_type='all'\n",
    "\n",
    "\n",
    "if plot_type == 'all':\n",
    "    sir_array = MOD_DOY_df['Avg']\n",
    "    sir_title = sensor_SIR + '-' + Site + ' rSIR MOD (Avg DOY '+str(subYears[0])+'-'+str(subYears[-1])+')'\n",
    "    grd_array = MOD_DOY_GRD_df['Avg']\n",
    "    grd_title = sensor_GRD + '-' + Site + ' GRD MOD (Avg DOY '+str(subYears[0])+'-'+str(subYears[-1])+')'\n",
    "elif plot_type == 'year':\n",
    "    sir_array = MOD_DOY_df[year]\n",
    "    sir_title = sensor_SIR + '-' + Site + ' rSIR MOD (' + str(year) + ' DOY)'\n",
    "    grd_array = MOD_DOY_GRD_df[year]\n",
    "    grd_title = sensor_GRD + '-' + Site + ' GRD MOD (' + str(year) + ' DOY)'\n",
    "\n",
    "fig, axes = plt.subplots(1, 2, figsize=(14,5))\n",
    "\n",
    "# Set a few common things\n",
    "graticule_fontsize = 8\n",
    "\n",
    "#min day of year to plot\n",
    "minday=75\n",
    "#max day of year to plot\n",
    "maxday=175\n",
    "\n",
    "# create polar stereographic Basemap instance.\n",
    "axes[0].set_title(grd_title)\n",
    "m = Basemap(ax=axes[0],projection='stere',\n",
    "            lon_0=(lon_start+lon_end)/2,\n",
    "            lat_0=(lat_start+lat_end)/2,lat_ts=90.,\n",
    "            llcrnrlat=(lat_start-.3),\n",
    "            urcrnrlat=(lat_end+.5),\n",
    "            llcrnrlon=(lon_start-.3),\n",
    "            urcrnrlon=(lon_end+.3),\n",
    "            rsphere=6371200.,\n",
    "            resolution='l',\n",
    "            area_thresh=10000, \n",
    "            epsg=3857)\n",
    "m.drawcoastlines()\n",
    "m.drawcountries()\n",
    "#m.drawstates()\n",
    "#m.etopo(scale=3, alpha=0.7)\n",
    "#m.bluemarble(scale=1)\n",
    "parallels = np.arange(0.,90,0.5)\n",
    "m.drawparallels(parallels, labels=[1,0,0,0], fontsize=graticule_fontsize)\n",
    "# draw meridians\n",
    "meridians = np.arange(0,360.,0.5)\n",
    "m.drawmeridians(meridians, labels=[0,0,0,1], fontsize=graticule_fontsize)\n",
    "m.arcgisimage(service='World_Physical_Map', xpixels = 1500, verbose= True)\n",
    "\n",
    "#m.readshapefile('/home/mij216/ExploringCETB/IN_Hunza_at_DainyorBridge', 'Hunza', color='red', linewidth=1)\n",
    "#m.readshapefile('/mnt/data3/rgi6.0/06_rgi60_Iceland', 'Vatnajokull', color='red', linewidth=1)\n",
    "\n",
    "x,y = m(pixel_lons_GRD, pixel_lats_GRD)\n",
    "im0 = m.scatter(x, y, c=grd_array, \n",
    "                s=sGRD, marker='s',lw=0,cmap='BuPu_r', alpha=.6, vmin=minday, vmax=maxday) #RdBu, \n",
    "\n",
    "# Get the current axes and configure placement so colorbar will\n",
    "# be 5% of ax width and the padding will be 0.08 in.\n",
    "divider = make_axes_locatable(axes[0])\n",
    "cax = divider.append_axes(\"right\", size=\"5%\", pad=0.08)\n",
    "cbar = plt.colorbar(im0, cax=cax, label='DOY')\n",
    "cbar.set_clim(minday, maxday)  #color bar limits\n",
    "\n",
    "axes[1].set_title(sir_title)\n",
    "m = Basemap(ax=axes[1],projection='stere',\n",
    "            lon_0=(lon_start+lon_end)/2,\n",
    "            lat_0=(lat_start+lat_end)/2,lat_ts=90.,\n",
    "            llcrnrlat=(lat_start-.3),\n",
    "            urcrnrlat=(lat_end+.5),\n",
    "            llcrnrlon=(lon_start-.3),\n",
    "            urcrnrlon=(lon_end+.3),\n",
    "            rsphere=6371200.,\n",
    "            resolution='l',\n",
    "            area_thresh=10000, \n",
    "            epsg=3857)\n",
    "m.drawcoastlines()\n",
    "m.drawcountries()\n",
    "#m.drawstates()\n",
    "#m.etopo(scale=3, alpha=0.7)\n",
    "#m.bluemarble(scale=1)\n",
    "parallels = np.arange(0.,90,0.5)\n",
    "m.drawparallels(parallels, labels=[1,0,0,0], fontsize=10)\n",
    "# draw meridians\n",
    "meridians = np.arange(0,360.,0.5)\n",
    "m.drawmeridians(meridians, labels=[0,0,0,1], fontsize=10)\n",
    "m.arcgisimage(service='World_Physical_Map', xpixels = 1500, verbose= True)\n",
    "\n",
    "#m.readshapefile('/home/mij216/ExploringCETB/IN_Hunza_at_DainyorBridge', 'Hunza', color='red', linewidth=1)\n",
    "#m.readshapefile('/mnt/data3/rgi6.0/06_rgi60_Iceland', 'Vatnajokull', color='red', linewidth=1)\n",
    "\n",
    "x,y = m(pixel_lons, pixel_lats)\n",
    "im1 = m.scatter(x, y, c=sir_array, \n",
    "                s=sSIR, marker='s', lw=0,cmap='BuPu_r', alpha=.6,vmin=minday,vmax=maxday) #RdBu, \n",
    "#fig.colorbar(im1, ax=axes[1],label='DOY')\n",
    "#plt.clim(75,175)  #color bar limits\n",
    "\n",
    "# Get the current axes and configure placement so colorbar will\n",
    "# be 5% of ax width and the padding will be 0.05 in.\n",
    "divider = make_axes_locatable(axes[1])\n",
    "cax = divider.append_axes(\"right\", size=\"5%\", pad=0.05)\n",
    "cbar = plt.colorbar(im1, cax=cax, label='DOY')\n",
    "cbar.set_clim(minday, maxday)  #color bar limits\n",
    "\n",
    "plt.tight_layout(pad=3.0)\n",
    "plt.show()\n",
    "filename = \"%s/%s.%s.%s.%s.%s.%s2.png\" % (\n",
    "    modsDir, region, areaname, sat_SIR, channel_SIR, alg_SIR, label)\n",
    "plt.savefig(filename)\n",
    "print(\"avg saved to %s\\n\" % filename)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "# PLOT the MOD MAP - avg of all years or one year of interest - GRD\n",
    "# Choose plot type - 'year' for one year of interest, 'all' for average of all years\n",
    "plot_type='year'  \n",
    "#plot_type='all'\n",
    "year=2018\n",
    "\n",
    "if plot_type=='all':\n",
    "    array=MOD_DOY_array_GRD\n",
    "    title=sensor_SIR+'-'+Site+' - MOD - (Avg DOY '+str(Years[0])+'-'+str(Years[-1])+')'\n",
    "elif plot_type=='year':\n",
    "    array=MOD_DOY_array_GRD_year\n",
    "    title=sensor_SIR+'-'+Site+' - MOD - ('+str(year)+' DOY)'\n",
    "\n",
    "fig = plt.figure(figsize=(8,8))\n",
    "# create figure and axes instances\n",
    "#fig = plt.figure()\n",
    "ax = fig.add_axes([0.1,0.1,0.8,0.8])\n",
    "#ax = fig.add_subplot(121)\n",
    "# create polar stereographic Basemap instance.\n",
    "m = Basemap(projection='stere',lon_0=(lon_start+lon_end)/2,lat_0=(lat_start+lat_end)/2,lat_ts=90.,\\\n",
    "            llcrnrlat=(lat_start-.3),urcrnrlat=(lat_end+.5),\\\n",
    "            llcrnrlon=(lon_start-.3),urcrnrlon=(lon_end+.3),\\\n",
    "            rsphere=6371200.,resolution='l',area_thresh=10000, epsg=3857)\n",
    "m.drawcoastlines()\n",
    "m.drawcountries()\n",
    "#m.drawstates()\n",
    "#m.etopo(scale=3, alpha=0.7)\n",
    "#m.bluemarble(scale=1)\n",
    "parallels = np.arange(0.,90,0.5)\n",
    "m.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)\n",
    "# draw meridians\n",
    "meridians = np.arange(0,360.,0.5)\n",
    "m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)\n",
    "m.arcgisimage(service='World_Physical_Map', xpixels = 1500, verbose= True)\n",
    "\n",
    "#m.readshapefile('/home/mij216/ExploringCETB/IN_Hunza_at_DainyorBridge', 'Hunza', color='red', linewidth=1)\n",
    "m.readshapefile('/mnt/data3/rgi6.0/06_rgi60_Iceland', 'Vatnajokull', color='red', linewidth=1)\n",
    "\n",
    "x,y=m(pixel_lons_GRD, pixel_lats_GRD)\n",
    "m.scatter(x,y,c=array, s=1200, marker='s',lw=0,cmap='BuPu_r', alpha=.6) #RdBu, \n",
    "\n",
    "plt.title(title)\n",
    "\n",
    "plt.colorbar(label='DOY')\n",
    "plt.clim(75,175)  #color bar limits\n",
    "\n",
    "plt.show()\n",
    "\n",
    "\n",
    "# plot another blank map with no pixel overlay\n",
    "fig = plt.figure(figsize=(8,8))\n",
    "ax = fig.add_axes([0.1,0.1,0.8,0.8])\n",
    "\n",
    "# create polar stereographic Basemap instance.\n",
    "m = Basemap(projection='stere',lon_0=(lon_start+lon_end)/2,lat_0=(lat_start+lat_end)/2,lat_ts=90.,\\\n",
    "            llcrnrlat=(lat_start-.3),urcrnrlat=(lat_end+.5),\\\n",
    "            llcrnrlon=(lon_start-.3),urcrnrlon=(lon_end+.3),\\\n",
    "            rsphere=6371200.,resolution='l',area_thresh=10000, epsg=3857)\n",
    "m.drawcoastlines()\n",
    "m.drawcountries()\n",
    "m.drawstates()\n",
    "#m.etopo(scale=5, alpha=0.7)\n",
    "#m.bluemarble(scale=1)\n",
    "parallels = np.arange(0.,90,0.5)\n",
    "m.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)\n",
    "# draw meridians\n",
    "meridians = np.arange(180.,360.,0.5)\n",
    "m.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)\n",
    "m.arcgisimage(service='World_Physical_Map', xpixels = 1500, verbose= True)\n",
    "#m.readshapefile('/home/mij216/ExploringCETB/IN_Hunza_at_DainyorBridge', 'Hunza', color='red', linewidth=1)\n",
    "m.readshapefile('/mnt/data3/rgi6.0/06_rgi60_Iceland', 'Vatnajokull', color='red', linewidth=1)\n",
    "\n",
    "plt.title(title)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#setup a dataframe so that all the data from ssmi and ssmis are read into one array for long term average\n",
    "#need to decide which sensor has priority when there is overlap\n",
    "#can create a reader that loops through enough sensors to get a full time series \n",
    "#(skip F10- orbit very elliptical and F19 - short)\n",
    "#create that time series for the analysis, save for analysis\n",
    "#\n",
    "#save dataframes with MOD so that we can do other analysys\n",
    "#pickle the data frame that is inside the MOD_array_DOY function \n",
    "#\n",
    "#check in CETB_analysis.py that when it does the count/window it is in fact \n",
    "#saving the first occurrence of the count, even though it doesn't know until it gets to 3x"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Add to this\n",
    "#Plot of which pixels have melted as of a specific day"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#figure out refreeze? then from MOD to Freeze, count # of dates with Tb>=threshold (with or without DAV? might depend on site)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Add to this\n",
    "#Plot of how many days (or occurrences) experienced melt \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Add to this\n",
    "#Plot of how many short (1-2 occurrences) melt days EMEs prior to the first MOD\n",
    "#From Jan 1 - MOD, Number of days with Tb>= threshold and DAV>= threshold (same thresholds)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}