ta_analysis.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import arcpy\n",
    "import pandas as pd # pandas upgrade required to use .query, PROPY env does not like that version \n",
    "import geopandas as gpd\n",
    "import os\n",
    "import csv\n",
    "\n",
    "import sys\n",
    "sys.path.insert(0,'M:/code/USFS')\n",
    "from helpers import set_years, MakeIterDir"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'M:\\\\projects\\\\data_processing'"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "os.getcwd()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "ename": "ModuleNotFoundError",
     "evalue": "No module named 'ee'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mModuleNotFoundError\u001b[0m                       Traceback (most recent call last)",
      "In  \u001b[0;34m[13]\u001b[0m:\nLine \u001b[0;34m1\u001b[0m:     \u001b[34mimport\u001b[39;49;00m \u001b[04m\u001b[36mee\u001b[39;49;00m\n",
      "\u001b[0;31mModuleNotFoundError\u001b[0m: No module named 'ee'\n\u001b[0;31m---------------------------------------------------------------------------\u001b[0m"
     ]
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1  counties in codf\n"
     ]
    }
   ],
   "source": [
    "lcdates_csv = r\"M:\\projects\\tc_fact_sheet\\data\\county_list_legendloc.csv\"\n",
    "codf = pd.read_csv(lcdates_csv)\n",
    "\n",
    "omit = ['wash_11001']\n",
    "maintain = ['bedf_51515', 'lanc_51103']\n",
    "codf = codf[~codf['cofips'].isin(omit) & codf['cofips'].isin(maintain)]\n",
    "codf = codf[codf.Size.isin(['s','xs'])]\n",
    "print(len(codf), \" counties in codf\")\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "#WITH SAFETIES\n",
    "\n",
    "pngDIR = \"M:/projects/tc_fact_sheet/map_exports/auto\"\n",
    "IterationDir = MakeIterDir(pngDIR, \"test\")\n",
    "\n",
    "fail_list = []\n",
    "counter = 0\n",
    "\n",
    "for index, row in codf.iterrows():\n",
    "    cf = row['cofips']\n",
    "    T1 = row['T1']\n",
    "    T2 = row['T2']\n",
    "    QA = row['QA']\n",
    "    size = row['Size']\n",
    "    \n",
    "    if size == \n",
    "    \n",
    "    FIPS = cf.split('_')[1]\n",
    "    counter += 1\n",
    "#     if counter > 5: # for testing\n",
    "#         break\n",
    "    print(f\"{cf} {counter}/{len(codf)}\", end='\\r')\n",
    "\n",
    "    FIPS = cf.split('_')[1]\n",
    "    tempDir = f\"M:/projects/tc_fact_sheet/data/temp\"\n",
    "    coDir = f\"M:/projects/tc_fact_sheet/data/prod/{cf}_{hex_ac}\"\n",
    "    CoBoundsPath = r\"J:/GIS/CBP_Obj_1/data/county_sets/tl_2017_bay_county.shp\"\n",
    "\n",
    "    if not os.path.isdir(coDir):\n",
    "        os.mkdir(coDir)\n",
    "\n",
    "\n",
    "    CoSel = os.path.join(tempDir,f\"{cf}_bounds.shp\")\n",
    "\n",
    "    where_clause = f\"\"\"GEOID = '{FIPS}'\"\"\"\n",
    "    print(where_clause)\n",
    "\n",
    "    BayHexPath = f\"M:/projects/data_processing/data_processing.gdb/baywide_{hex_ac}\"\n",
    "    CoHexPath = os.path.join(coDir,f\"{cf}_hex.shp\")\n",
    "    calcCSVPath = os.path.join(coDir,f\"{cf}_ta_calc.csv\")\n",
    "    ta_dbf_path = os.path.join(coDir,f\"{cf}_hex_ta.dbf\")\n",
    "    ta_xlsx_path = ta_dbf_path.replace(\".dbf\",\".xlsx\")\n",
    "    lulc_chgPath = f\"B:/landuse/version2/{cf}/output/{cf}_landusechange_{T1}{T2}.tif\"\n",
    "\n",
    "    t2luPath = f\"B:/landuse/version2/{cf}/output/{cf}_lu_{T2}.tif\"\n",
    "    t2_dbf_path = os.path.join(coDir,f\"{cf}_hex_t2_ta.dbf\")\n",
    "    t2_xlsx_path = t2_dbf_path.replace(\".dbf\",\".xlsx\")\n",
    "    t2_csv_path = t2_xlsx_path.replace('xlsx','csv')\n",
    "\n",
    "    tajoined3Path = os.path.join(coDir,f\"{cf}_hex_tajoined3.shp\")\n",
    "\n",
    "    if not arcpy.Exists(CoSel):\n",
    "        sel_temp = arcpy.management.SelectLayerByAttribute(CoBoundsPath, \"NEW_SELECTION\", where_clause)\n",
    "        arcpy.FeatureClassToFeatureClass_conversion(sel_temp, tempDir, os.path.basename(CoSel))\n",
    "    else:\n",
    "        print(f\"{CoSel} exists.\")\n",
    "\n",
    "    if not arcpy.Exists(CoHexPath):\n",
    "        tempHex = arcpy.SelectLayerByLocation_management(BayHexPath, 'INTERSECT', CoSel)\n",
    "        arcpy.FeatureClassToFeatureClass_conversion(tempHex, coDir, os.path.basename(CoHexPath))\n",
    "        del tempHex\n",
    "\n",
    "    \n",
    "    if not os.path.exists(ta_dbf_path):\n",
    "        print(f\"TA - {ta_dbf_path}\")\n",
    "        arcpy.sa.TabulateArea(CoHexPath, \n",
    "                            \"GRID_ID\", \n",
    "                            lulc_chgPath, \n",
    "                            \"Value\", \n",
    "                            ta_dbf_path, \n",
    "                            lulc_chgPath, \n",
    "                            \"CLASSES_AS_FIELDS\")\n",
    "\n",
    "\n",
    "    if not os.path.exists(t2_dbf_path):\n",
    "        print(f\"TA - {ta_dbf_path}\")\n",
    "        arcpy.sa.TabulateArea(CoHexPath, \n",
    "                            \"GRID_ID\", \n",
    "                            t2luPath, \n",
    "                            \"Value\", \n",
    "                            t2_dbf_path, \n",
    "                            t2luPath, \n",
    "                            \"CLASSES_AS_FIELDS\")\n",
    "\n",
    "    # Convert to excel format to visually inspect easier and read into pandas\n",
    "    if not os.path.exists(ta_xlsx_path):\n",
    "        arcpy.conversion.TableToExcel(ta_dbf_path, ta_xlsx_path, \"NAME\", \"CODE\")\n",
    "    if not os.path.exists(t2_xlsx_path):   \n",
    "        arcpy.conversion.TableToExcel(t2_dbf_path, t2_xlsx_path, \"NAME\", \"CODE\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "alex_51510\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\jczawlytko\\.conda\\envs\\open-gis\\lib\\site-packages\\ipykernel_launcher.py:89: UserWarning: Column names longer than 10 characters will be truncated when saved to ESRI Shapefile.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "complete\n",
      "[]\n",
      "arli_51013\n",
      "complete\n",
      "[]\n",
      "balt_24510\n",
      "complete\n",
      "[]\n",
      "buen_51530\n",
      "complete\n",
      "[]\n",
      "char_51540\n",
      "complete\n",
      "[]\n",
      "colo_51570\n",
      "complete\n",
      "[]\n",
      "covi_51580\n",
      "complete\n",
      "[]\n",
      "fair_51600\n",
      "complete\n",
      "[]\n",
      "fall_51610\n",
      "complete\n",
      "[]\n",
      "fred_51630\n",
      "complete\n",
      "[]\n",
      "hamp_51650\n",
      "complete\n",
      "[]\n",
      "harr_51660\n",
      "complete\n",
      "[]\n",
      "hope_51670\n",
      "complete\n",
      "[]\n",
      "lexi_51678\n",
      "complete\n",
      "[]\n",
      "lync_51680\n",
      "complete\n",
      "[]\n",
      "mana_51683\n",
      "complete\n",
      "[]\n",
      "mana_51685\n",
      "complete\n",
      "[]\n",
      "mont_42093\n",
      "complete\n",
      "[]\n",
      "newp_51700\n",
      "complete\n",
      "[]\n",
      "norf_51710\n",
      "complete\n",
      "[]\n",
      "pete_51730\n",
      "complete\n",
      "[]\n",
      "poqu_51735\n",
      "complete\n",
      "[]\n",
      "port_51740\n",
      "complete\n",
      "[]\n",
      "rich_51760\n",
      "complete\n",
      "[]\n",
      "roan_51770\n",
      "complete\n",
      "[]\n",
      "sale_51775\n",
      "complete\n",
      "[]\n",
      "stau_51790\n",
      "complete\n",
      "[]\n",
      "wayn_51820\n",
      "complete\n",
      "[]\n",
      "will_51830\n",
      "complete\n",
      "[]\n"
     ]
    }
   ],
   "source": [
    "elist = [] #clear elist\n",
    "counter = 0\n",
    "\n",
    "tc = [24,25,26,27,41,42]\n",
    "dev = [21,22,23,28,29,33,34,35,36,37,38,51,52,53]\n",
    "wat = [11]\n",
    "\n",
    "for index, row in codf.iterrows():\n",
    "    cf = row['cofips']\n",
    "    print(cf)\n",
    "    T1 = row['T1']\n",
    "    T2 = row['T2']\n",
    "    QA = row['QA']\n",
    "    size = row['Size']\n",
    "    FIPS = cf.split('_')[1]\n",
    "    T1,T2 = set_years(cf)\n",
    "\n",
    "    counter += 1    \n",
    "\n",
    "    # BUILD CLASS so you don't have to define these for each loop\n",
    "    coDir = f\"M:/projects/tc_fact_sheet/data/prod/{cf}_{hex_ac}\"\n",
    "    CoHexPath = os.path.join(coDir,f\"{cf}_hex.shp\")\n",
    "    BayHexPath = f\"M:/projects/data_processing/data_processing.gdb/baywide_{hex_ac}\"\n",
    "    calcCSVPath = os.path.join(coDir,f\"{cf}_ta_calc.csv\")\n",
    "    ta_dbf_path = os.path.join(coDir,f\"{cf}_hex_ta.dbf\")\n",
    "    ta_xlsx_path = ta_dbf_path.replace(\".dbf\",\".xlsx\")\n",
    "    lulc_chgPath = f\"B:/landuse/version2/{cf}/output/{cf}_landusechange_{T1}{T2}.tif\"\n",
    "\n",
    "    t2luPath = f\"B:/landuse/version2/{cf}/output/{cf}_lu_{T2}.tif\"\n",
    "    t2_dbf_path = os.path.join(coDir,f\"{cf}_hex_t2_ta.dbf\")\n",
    "    t2_xlsx_path = t2_dbf_path.replace(\".dbf\",\".xlsx\")\n",
    "    t2_csv_path = t2_xlsx_path.replace('xlsx','csv')\n",
    "    tajoined3Path = os.path.join(coDir,f\"{cf}_hex_tajoined3.shp\")\n",
    "\n",
    "    # #read excel files to be edited\n",
    "    # tadf = pd.read_excel(ta_xlsx_path)\n",
    "    # t2df = pd.read_excel(t2_xlsx_path)\n",
    "\n",
    "    # #modify ta results for LULC Change\n",
    "    # tadf['tc_chg'] = 0\n",
    "    # tavals = [col for col in tadf.columns if col[0:5] == \"VALUE\"]\n",
    "    # loss_vals = []\n",
    "    # gain_vals = []\n",
    "    # # for tv in list(tadf.VALUE.unique()):\n",
    "    # for tv in tavals:\n",
    "    # ## If it contains the VALUE prefix use this section\n",
    "    #     T1val = int(tv[6:8])\n",
    "    #     T2val = int(tv[8:10])\n",
    "    # #     T1val = int(str(tv)[0:2])\n",
    "    # #     T2val = int(str(tv)[2:4])\n",
    "    #     if T1val in tc and T2val in dev:\n",
    "    #         loss_vals.append(tv)\n",
    "    #     elif T1val in dev and T2val in tc:\n",
    "    #         gain_vals.append(tv)\n",
    "    #     else:\n",
    "    #         pass\n",
    "\n",
    "    # tadf['tc_loss']=tadf[loss_vals].sum(axis=1)\n",
    "    # tadf['tc_gain']=tadf[gain_vals].sum(axis=1)\n",
    "    # tadf['tc_chg']=tadf['tc_gain']-tadf['tc_loss']\n",
    "\n",
    "    # tadf_clean = tadf[['GRID_ID','tc_chg','tc_gain','tc_loss']]\n",
    "    # tadf_clean.to_csv(calcCSVPath,index=True)\n",
    "\n",
    "    # # modify ta results for T2 LULC\n",
    "\n",
    "    # if not os.path.exists(t2_csv_path):  \n",
    "    #     t2df = pd.read_excel(t2_xlsx_path)\n",
    "    #     t2vals = [col for col in t2df.columns if col[0:5] == \"VALUE\"]\n",
    "        \n",
    "    #     possible_wat_vals = [\"VALUE_11\",\"VALUE_12\",\"VALUE_13\",\"VALUE_14\",\"VALUE_15\"]\n",
    "    #     wat_vals = [i for i in t2df.columns if i in possible_wat_vals]\n",
    "\n",
    "    #     t2df['wat_lu'] = t2df[wat_vals].sum(axis=1)\n",
    "    #     t2df['all_lu']= t2df[t2vals].sum(axis=1)\n",
    "    #     t2df.to_csv(t2_csv_path)\n",
    "    # else:\n",
    "    #     t2df = pd.read_excel(t2_xlsx_path)\n",
    "\n",
    "\n",
    "    # if not os.path.exists(tajoined3Path):  \n",
    "    t2_csv_df = pd.read_csv(t2_csv_path)\n",
    "    ta_csv_df = pd.read_csv(calcCSVPath)\n",
    "    hex_df = gpd.read_file(CoHexPath)\n",
    "\n",
    "    hex_df = hex_df.merge(t2_csv_df,how=\"left\",left_on=\"GRID_ID\",right_on=\"GRID_ID\")\n",
    "    hex_df = hex_df.merge(ta_csv_df,how=\"left\",left_on=\"GRID_ID\",right_on=\"GRID_ID\")\n",
    "\n",
    "    hex_df.to_file(tajoined3Path)\n",
    "\n",
    "\n",
    "print('complete')\n",
    "print(elist)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Shape_Leng</th>\n",
       "      <th>Shape_Area</th>\n",
       "      <th>GRID_ID</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Unnamed: 0_x</th>\n",
       "      <th>OID</th>\n",
       "      <th>VALUE_13</th>\n",
       "      <th>VALUE_14</th>\n",
       "      <th>VALUE_15</th>\n",
       "      <th>VALUE_21</th>\n",
       "      <th>...</th>\n",
       "      <th>VALUE_74</th>\n",
       "      <th>VALUE_83</th>\n",
       "      <th>VALUE_84</th>\n",
       "      <th>VALUE_85</th>\n",
       "      <th>wat_lu</th>\n",
       "      <th>all_lu</th>\n",
       "      <th>Unnamed: 0_y</th>\n",
       "      <th>tc_chg</th>\n",
       "      <th>tc_gain</th>\n",
       "      <th>tc_loss</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2368.014715</td>\n",
       "      <td>404685.999091</td>\n",
       "      <td>UF-897</td>\n",
       "      <td>POLYGON ((1517507.388 1932039.878, 1517310.054...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>573</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15323</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2368.014715</td>\n",
       "      <td>404685.999091</td>\n",
       "      <td>UF-896</td>\n",
       "      <td>POLYGON ((1517507.388 1932723.465, 1517310.054...</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>6882</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>22671</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>214977</td>\n",
       "      <td>0.0</td>\n",
       "      <td>114.0</td>\n",
       "      <td>279.0</td>\n",
       "      <td>165.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2368.014642</td>\n",
       "      <td>404685.974073</td>\n",
       "      <td>UD-895</td>\n",
       "      <td>POLYGON ((1516323.381 1933407.052, 1516126.046...</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>16</td>\n",
       "      <td>0</td>\n",
       "      <td>31160</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>149593</td>\n",
       "      <td>0</td>\n",
       "      <td>16</td>\n",
       "      <td>334349</td>\n",
       "      <td>1.0</td>\n",
       "      <td>907.0</td>\n",
       "      <td>930.0</td>\n",
       "      <td>23.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2368.014542</td>\n",
       "      <td>404685.939892</td>\n",
       "      <td>UE-895</td>\n",
       "      <td>POLYGON ((1516915.385 1933065.258, 1516718.050...</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>5220</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>187669</td>\n",
       "      <td>2.0</td>\n",
       "      <td>286.0</td>\n",
       "      <td>286.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2368.014542</td>\n",
       "      <td>404685.939892</td>\n",
       "      <td>UF-895</td>\n",
       "      <td>POLYGON ((1517507.388 1933407.052, 1517310.054...</td>\n",
       "      <td>4</td>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>18557</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>404849</td>\n",
       "      <td>3.0</td>\n",
       "      <td>1936.0</td>\n",
       "      <td>2214.0</td>\n",
       "      <td>278.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74</th>\n",
       "      <td>2368.014715</td>\n",
       "      <td>404685.999092</td>\n",
       "      <td>UF-885</td>\n",
       "      <td>POLYGON ((1517507.388 1940242.921, 1517310.054...</td>\n",
       "      <td>74</td>\n",
       "      <td>74</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8782</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>72</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>101799</td>\n",
       "      <td>72.0</td>\n",
       "      <td>-195.0</td>\n",
       "      <td>105.0</td>\n",
       "      <td>300.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75</th>\n",
       "      <td>2368.014815</td>\n",
       "      <td>404686.033272</td>\n",
       "      <td>UG-885</td>\n",
       "      <td>POLYGON ((1518099.392 1939901.128, 1517902.057...</td>\n",
       "      <td>75</td>\n",
       "      <td>75</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>62955</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>336695</td>\n",
       "      <td>73.0</td>\n",
       "      <td>69.0</td>\n",
       "      <td>1000.0</td>\n",
       "      <td>931.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>76</th>\n",
       "      <td>2368.014815</td>\n",
       "      <td>404686.033272</td>\n",
       "      <td>UH-885</td>\n",
       "      <td>POLYGON ((1518691.396 1940242.921, 1518494.061...</td>\n",
       "      <td>76</td>\n",
       "      <td>76</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>256</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77</th>\n",
       "      <td>2368.014815</td>\n",
       "      <td>404686.033271</td>\n",
       "      <td>UC-884</td>\n",
       "      <td>POLYGON ((1515731.377 1940584.715, 1515534.043...</td>\n",
       "      <td>77</td>\n",
       "      <td>77</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>2901</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3699</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>78</th>\n",
       "      <td>2368.014815</td>\n",
       "      <td>404686.033272</td>\n",
       "      <td>UG-884</td>\n",
       "      <td>POLYGON ((1518099.392 1940584.715, 1517902.057...</td>\n",
       "      <td>78</td>\n",
       "      <td>78</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>13023</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>97113</td>\n",
       "      <td>74.0</td>\n",
       "      <td>99.0</td>\n",
       "      <td>219.0</td>\n",
       "      <td>120.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>79 rows × 38 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     Shape_Leng     Shape_Area GRID_ID  \\\n",
       "0   2368.014715  404685.999091  UF-897   \n",
       "1   2368.014715  404685.999091  UF-896   \n",
       "2   2368.014642  404685.974073  UD-895   \n",
       "3   2368.014542  404685.939892  UE-895   \n",
       "4   2368.014542  404685.939892  UF-895   \n",
       "..          ...            ...     ...   \n",
       "74  2368.014715  404685.999092  UF-885   \n",
       "75  2368.014815  404686.033272  UG-885   \n",
       "76  2368.014815  404686.033272  UH-885   \n",
       "77  2368.014815  404686.033271  UC-884   \n",
       "78  2368.014815  404686.033272  UG-884   \n",
       "\n",
       "                                             geometry  Unnamed: 0_x  OID  \\\n",
       "0   POLYGON ((1517507.388 1932039.878, 1517310.054...             0    0   \n",
       "1   POLYGON ((1517507.388 1932723.465, 1517310.054...             1    1   \n",
       "2   POLYGON ((1516323.381 1933407.052, 1516126.046...             2    2   \n",
       "3   POLYGON ((1516915.385 1933065.258, 1516718.050...             3    3   \n",
       "4   POLYGON ((1517507.388 1933407.052, 1517310.054...             4    4   \n",
       "..                                                ...           ...  ...   \n",
       "74  POLYGON ((1517507.388 1940242.921, 1517310.054...            74   74   \n",
       "75  POLYGON ((1518099.392 1939901.128, 1517902.057...            75   75   \n",
       "76  POLYGON ((1518691.396 1940242.921, 1518494.061...            76   76   \n",
       "77  POLYGON ((1515731.377 1940584.715, 1515534.043...            77   77   \n",
       "78  POLYGON ((1518099.392 1940584.715, 1517902.057...            78   78   \n",
       "\n",
       "    VALUE_13  VALUE_14  VALUE_15  VALUE_21  ...  VALUE_74  VALUE_83  VALUE_84  \\\n",
       "0          0         0         0         0  ...         0         0       573   \n",
       "1          0         0         0      6882  ...         0         0     22671   \n",
       "2          0        16         0     31160  ...         0         0    149593   \n",
       "3          0         0         0      5220  ...         0         0         0   \n",
       "4          0         0         0     18557  ...         0         0         0   \n",
       "..       ...       ...       ...       ...  ...       ...       ...       ...   \n",
       "74         0         0         0      8782  ...         0         0        72   \n",
       "75         0         0         0     62955  ...         0         0         0   \n",
       "76         0         0         0        18  ...         0         0         0   \n",
       "77         0         0         0         0  ...         0         0      2901   \n",
       "78         0         0         0     13023  ...         0         0         0   \n",
       "\n",
       "    VALUE_85  wat_lu  all_lu  Unnamed: 0_y  tc_chg  tc_gain  tc_loss  \n",
       "0          0       0   15323           NaN     NaN      NaN      NaN  \n",
       "1          0       0  214977           0.0   114.0    279.0    165.0  \n",
       "2          0      16  334349           1.0   907.0    930.0     23.0  \n",
       "3          0       0  187669           2.0   286.0    286.0      0.0  \n",
       "4          0       0  404849           3.0  1936.0   2214.0    278.0  \n",
       "..       ...     ...     ...           ...     ...      ...      ...  \n",
       "74         5       0  101799          72.0  -195.0    105.0    300.0  \n",
       "75         0       0  336695          73.0    69.0   1000.0    931.0  \n",
       "76         0       0     256           NaN     NaN      NaN      NaN  \n",
       "77         0       0    3699           NaN     NaN      NaN      NaN  \n",
       "78         0       0   97113          74.0    99.0    219.0    120.0  \n",
       "\n",
       "[79 rows x 38 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "hex_df"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ArcGISPro",
   "language": "Python",
   "name": "python3"
  },
  "language_info": {
   "file_extension": ".py",
   "name": "python",
   "version": "3"
  },
  "vscode": {
   "interpreter": {
    "hash": "779ee9392788e65623d032f1ff917715589ea2f4728dac80fd4ab630ca72779b"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}