4-Land-Use-Calcs.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "import _config"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "input_taz_shapefile_folder = _config.input_taz_shapefile_folder\n",
    "taz_shapefile              = _config.taz_shapefile\n",
    "se_years                   = _config.se_years\n",
    "input_model_se_folder      = _config.input_model_se_folder\n",
    "se_base_year               = _config.se_base_year\n",
    "df_centers                 = _config.df_centers\n",
    "df_projects                = _config.df_projects\n",
    "hh_factor                  = _config.hh_factor"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Read in data from previous notebooks"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>nearest_AreaType</th>\n",
       "      <th>TAZID</th>\n",
       "      <th>TAZACRES</th>\n",
       "      <th>sqm</th>\n",
       "      <th>acres</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1472</td>\n",
       "      <td>42.036264</td>\n",
       "      <td>49740.404222</td>\n",
       "      <td>12.291122</td>\n",
       "      <td>POLYGON ((424043.145 4505847.017, 424036.495 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1473</td>\n",
       "      <td>44.089344</td>\n",
       "      <td>52165.559080</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>POLYGON ((424421.505 4505865.604, 424414.855 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1481</td>\n",
       "      <td>68.022708</td>\n",
       "      <td>13056.790069</td>\n",
       "      <td>3.226403</td>\n",
       "      <td>POLYGON ((424229.000 4505722.832, 424209.282 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1482</td>\n",
       "      <td>39.922698</td>\n",
       "      <td>11968.229688</td>\n",
       "      <td>2.957414</td>\n",
       "      <td>POLYGON ((424396.265 4505812.237, 424384.505 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>15183</td>\n",
       "      <td>1</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>176</td>\n",
       "      <td>288.616036</td>\n",
       "      <td>35706.046305</td>\n",
       "      <td>8.823156</td>\n",
       "      <td>POLYGON ((412496.126 4566424.387, 412502.776 4...</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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>173</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1383</td>\n",
       "      <td>82.068433</td>\n",
       "      <td>163541.758588</td>\n",
       "      <td>40.412049</td>\n",
       "      <td>POLYGON ((418599.665 4505156.000, 418601.602 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>174</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1383</td>\n",
       "      <td>82.068433</td>\n",
       "      <td>10615.218005</td>\n",
       "      <td>2.623077</td>\n",
       "      <td>POLYGON ((418812.341 4505510.828, 418778.475 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>175</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>176</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>177</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1396</td>\n",
       "      <td>98.491949</td>\n",
       "      <td>129814.582427</td>\n",
       "      <td>32.077882</td>\n",
       "      <td>POLYGON ((419041.436 4504755.602, 419002.000 4...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>178 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         N  project_id     nearest_AreaType  TAZID    TAZACRES            sqm  \\\n",
       "0    15086           2  Neighborhood Center   1472   42.036264   49740.404222   \n",
       "1    15086           2  Neighborhood Center   1473   44.089344   52165.559080   \n",
       "2    15086           2  Neighborhood Center   1481   68.022708   13056.790069   \n",
       "3    15086           2  Neighborhood Center   1482   39.922698   11968.229688   \n",
       "4    15183           1  Neighborhood Center    176  288.616036   35706.046305   \n",
       "..     ...         ...                  ...    ...         ...            ...   \n",
       "173  15069           2         Urban Center   1383   82.068433  163541.758588   \n",
       "174  23506           2         Urban Center   1383   82.068433   10615.218005   \n",
       "175  15069           2         Urban Center   1384   96.889876  208196.248028   \n",
       "176  23506           2         Urban Center   1384   96.889876   34533.423293   \n",
       "177  15069           2         Urban Center   1396   98.491949  129814.582427   \n",
       "\n",
       "         acres                                           geometry  \n",
       "0    12.291122  POLYGON ((424043.145 4505847.017, 424036.495 4...  \n",
       "1    12.890390  POLYGON ((424421.505 4505865.604, 424414.855 4...  \n",
       "2     3.226403  POLYGON ((424229.000 4505722.832, 424209.282 4...  \n",
       "3     2.957414  POLYGON ((424396.265 4505812.237, 424384.505 4...  \n",
       "4     8.823156  POLYGON ((412496.126 4566424.387, 412502.776 4...  \n",
       "..         ...                                                ...  \n",
       "173  40.412049  POLYGON ((418599.665 4505156.000, 418601.602 4...  \n",
       "174   2.623077  POLYGON ((418812.341 4505510.828, 418778.475 4...  \n",
       "175  51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  \n",
       "176   8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  \n",
       "177  32.077882  POLYGON ((419041.436 4504755.602, 419002.000 4...  \n",
       "\n",
       "[178 rows x 8 columns]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP0AAAGsCAYAAAAWkKEaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAtgUlEQVR4nO3dfVxUZd4/8M+BgQFhGB4EGYRETcVAlLDCNVF36OEn6Va7bRbttm7d271Lre7+3J/LvQ/FbhttWTe5rbs3rLvdSEW1WLdld6EluvnMCIkaSD4xDig+MQMiZ4C5fn+YkxOgM4gNw/V5v17n9eJc5zpnvpfOZ86Zc2bOKEIIASKShp+3CyCirxdDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTySZYRf6zZs3Y/78+YiLi4OiKHjnnXc83oYQAsuXL8fEiROh1WqRkJCAZ555ZvCLJfICjbcLGGznzp3D1KlTsWjRInz7298e0DYWL16MiooKLF++HFOmTIHVasWpU6cGuVIiLxHDGADx9ttvu7Spqip+8YtfiLi4ODFixAhx8803i40bNzqX79+/X2g0GlFXV/f1Fkv0NRl2h/dXsmjRImzZsgVlZWXYs2cP7rvvPtx5551oaGgAALz77rsYN24c3nvvPYwdOxaJiYl49NFHcebMGS9XTjQ4pAr9wYMH8frrr+Ott97CrFmzMH78eCxduhS33nor/vGPfwAADh06hKNHj+Ktt95CSUkJXnnlFZhMJnznO9/xcvVEg2PYvae/nN27d0MIgYkTJ7q0q6qKqKgoAIDD4YCqqigpKXH2W7VqFdLT01FfX49JkyZ97XUTDSapQu9wOODv7w+TyQR/f3+XZaGhoQAAg8EAjUbj8sIwefJkAEBjYyNDTz5PqtCnpaWhp6cHLS0tmDVrVp99Zs6cie7ubhw8eBDjx48HABw4cAAAMGbMmK+tVqJrRRFieN0uq729HZ9//jmACyF/8cUXMXfuXERGRuK6667DQw89hC1btuCFF15AWloaTp06hY8//hhTpkzBvHnz4HA4cNNNNyE0NBSFhYVwOBzIzc1FWFgYKioqvDw6okHg5asHg27jxo0CQK/p4YcfFkIIYbfbxW9/+1uRmJgoAgICRGxsrLjnnnvEnj17nNuwWCzi3nvvFaGhoWLUqFHiBz/4gTh9+rSXRkQ0uIbdnp6ILk+qS3ZEdJWhLygogKIoWLJkSb99KisroShKr6murs7ZZ86cOX32yc7OvpryiKgPAz57v2vXLhQVFSE1NdWt/vX19QgLC3POR0dHO/9es2YN7Ha7c/706dOYOnUq7rvvPrfrcTgcaGpqgk6ng6Iobq9H5MuEEGhra0NcXBz8/Nzbhw8o9O3t7cjJyUFxcTGefvppt9aJiYlBeHh4n8siIyNd5svKyjBixAiPQt/U1ISEhAS3+xMNJ2azGfHx8W71HVDoc3NzkZ2djaysLLdDn5aWhs7OTtxwww349a9/jblz5/bbd9WqVVi4cCFCQkL67aOqKlRVdc5fPB9pNptdjiiIhjObzYaEhATodDq31/E49GVlZTCZTKiqqnKrv8FgQFFREdLT06GqKlavXg2j0YjKykpkZmb26r9z507s3bsXq1atuux2CwoKkJ+f36s9LCyMoSfpePKW1qNLdmazGdOnT0dFRQWmTp0K4MJJuGnTpqGwsNDtB50/fz4URcHatWt7LXvsscewdetW1NbWXnYbX93TX3zFs1qtDD1Jw2azQa/Xe/S89+jsvclkQktLC9LT06HRaKDRaLBp0yasWLECGo0GPT09bm0nIyPD+VXWS3V0dKCsrAyPPvroFbeh1Wqde3Xu3Ync59HhvdFo7LUHXrRoEZKSkrBs2bJeX2LpT3V1NQwGQ6/2N998E6qq4qGHHvKkLCLygEeh1+l0SElJcWkLCQlBVFSUsz0vLw8WiwUlJSUAgMLCQiQmJiI5ORl2ux2lpaUoLy9HeXl5r+2vWrUKd999t/NrrkQ0+Ab9W3bNzc1obGx0ztvtdixduhQWiwXBwcFITk7GunXrMG/ePJf1Dhw4gE8++YRfaiG6xobNZ+8HckKDyNdd8xN5ROT7GHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGalCbz7TgZc/bsDZc/YrdyYapqS67/3Kys/x+k4zAvz98Njs8d4uh8grpAr992ckItDfDwumxXm7FCKvkerwfrIhDPnfSoFBH9xrWbvajd2NZzGQTyULIXDk1Dk4HL77ieYTtk68VWWG+UyHt0uha0yq0F9O3ppa3LtyKz7cd8LjdVdvP4o5yyvxl00Hr9jX4RD42Rs1+I+3+79JyGs7jsL4QiU+b2nrc/n/1FjwWEkVjls7+1y+fv9x/O7d/ejscr2/wdaDp/DT16thPnvOpb3HIbD0rU+xevsRPFi8HR1q9xXHQb6Lof9CxrhI3GAIw/UxoR6vG6cPRmRIIBIiR1yxb7u9G+9+2oS3d1vQ3ePos8+BE+04ePIcGvvZ675dbcGH+0+g1mLtc/l/bTqEtZ9aUGtpdV1vtwVrP21CZf0pl/auHge2HzqNPcdsMJ89j3Y7Qz+c8Vt2XrDXYkWAvx8mxfZ9M8MOtRuHT59Dcpy+z+UnbJ341NyKrMmj4OfX+95o1Y1nsbfJhgdvvg7+lyw/YetEZX0LvjVtNIICXG94Um46htd3NuKuVAN+MHPsVYyOvk4Ded4z9HTVOrt68NwHdai1WDFj/EgsNk5webGha4dfrSWv+EvlQfx9yxF09wis+KgBb1aZvV0SXQZDT1ftVLuK1Hg9AjQKwkcE4NDJdm+XRJfB0NNVmzUhGnuOWVF7zAqHQ+DOlN43PaWhQ6oP59C1cWdKLMp/PAO7j7bi1glRmGzo+wQkDQ0MPQ2K9DGRSB8TeeWO5HU8vCeSDENPJBke3tOg6+5xoGjzIby6oxHTEsIxY1wkcjLGePQji3TtcE9Pg+5PH3+O5z6sR7QuECfbO/Gm6Rhe3dF45RXpa8HQ06D7nxrLhT+EAiEAf0X5so28jof3NOiiQrU4Z+9BZ3cPznb04LrIERip03q7LPoC9/Q06O5Lj8fJNhXW8104eqYD57t6MG8KP7AzVHBPT4Nu4c3XYYRWg48/O44eoeDetNGYmxTj7bLoCww9XRMLpsZhwVTelmwo4uE9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJ5qpCX1BQAEVRsGTJkn77VFZWQlGUXlNdXZ1Lv9bWVuTm5sJgMCAoKAiTJ0/G+++/fzXlSa/HIfCHdfvx0N92wHT0jLfLoSFiwF+t3bVrF4qKipCamupW//r6epcf2IuOjnb+bbfbcdtttyEmJgb//Oc/ER8fD7PZDJ2u7191JfdsO3gKH33WgviIYBRuOIA/PXAjllfUw9phx89vn4SxIz3/WW7yfQPa07e3tyMnJwfFxcWIiIhwa52YmBjExsY6J3//L38q+e9//zvOnDmDd955BzNnzsSYMWNw6623YurUqQMpj74QExaEiBGB2NxwCkEBGryxy4zS7Y040WbHSxsavF0eecmAQp+bm4vs7GxkZWW5vU5aWhoMBgOMRiM2btzosmzt2rWYMWMGcnNzMWrUKKSkpOCZZ55BT09Pv9tTVRU2m81lIlcTR+nwm7sm44/fnoL/vH8axkSNQNIoHXYePoProkK8XR55iceH92VlZTCZTKiqqnKrv8FgQFFREdLT06GqKlavXg2j0YjKykpkZmYCAA4dOoSPP/4YOTk5eP/999HQ0IDc3Fx0d3fjt7/9bZ/bLSgoQH5+vqflS2fadRGYdt2Fo7E7UwzQBWlgPd+NO5JjvVwZeYsihBDudjabzZg+fToqKiqch95z5szBtGnTUFhY6PaDzp8/H4qiYO3atQCAiRMnorOzE4cPH3Ye9r/44ot4/vnn0dzc3Oc2VFWFqqrOeZvNhoSEBFitVpdzB0TDmc1mg16v9+h579Ge3mQyoaWlBenp6c62np4ebN68GS+//DJUVXV5r96fjIwMlJaWOucNBgMCAgJc1p08eTKOHz8Ou92OwMDAXtvQarXQanlbZSJPeRR6o9GI2tpal7ZFixYhKSkJy5YtcyvwAFBdXQ2D4ctbIs+cOROvvfYaHA4H/PwunGY4cOAADAZDn4EnooHzKPQ6nQ4pKSkubSEhIYiKinK25+XlwWKxoKSkBABQWFiIxMREJCcnw263o7S0FOXl5SgvL3du48c//jH+9Kc/YfHixXjiiSfQ0NCAZ555Bj/96U+vdnzDnhAC1eaz0GkDMGEUL3HSlQ36LbCbm5vR2Pjl75bZ7XYsXboUFosFwcHBSE5Oxrp16zBv3jxnn4SEBFRUVOBnP/sZUlNTMXr0aCxevBjLli0b7PKGnVWfHMbKyoMI1Wrw7L1T8I3rR3q7JBriPDqRN5QN5ITGcPD0e/tQa7Gis8uBe28cjYe/MdbbJdHXaCDPe3723sctmDoaAgoCNX74Pyn86Si6Mv7CjY9LTQjHm4/N8HYZ5EO4pyeSDENPJBmGnkgyDD2RZBj6YaDx9DlYO7q8XQb5CIbexx09fQ5L3qhB7mu7vV0K+QhesvNx+uAARIYEIiYsyNulkI9g6H1c+IhA/O3hm7xdBvkQHt4TSYahJ5IMQ08kGYaeSDIMPZFkGHof19phx3Mf1OFfDSe9XQr5CIbex+04dBp1x23439rj3i6FfASv0/u4OUkxaD3fjfQx4d4uhXwEQ+/jtBp/3H9TgrfLIB/Cw3siyTD0RJJh6Ikkw9ATSYahHwY67N3eLoF8CEPv4z7Y24yHV+3AXzcd9HYp5CMYeh+nKAoUPz8o3i6EfAav0/u4O5JjMWvCSIwI5H8luYd7+mGAgSdPMPREkmHoiSTD0BNJhqEfBnp6HN4ugXwIzwD5MCEEHn+9GqfbVSTFhuGpBcneLol8APf0PuzY2fM42aYiVKvBviYrOvnJPHID9/Q+zKAPQoe9G0IIjAjUIIiX7sgNfJb4MI2/H97+yUwcPX0OY0eGersc8hEMvY8L8PfD9TE6b5dBPoTv6Ykkw9ATSYahJ5IMQ+/jdhw6jdXbjqDF1untUshHMPQ+7Mipc3jwbztQvtuCx0pN3i6HfMRVhb6goACKomDJkiX99qmsrLxwo4evTHV1dc4+r7zySp99Oju597ocW2cX9MEBCNT4wU/hbTTIPQO+ZLdr1y4UFRUhNTXVrf719fUICwtzzkdHR7ssDwsLQ319vUtbUFDQQMuTwpTReuTOvR5HTrXj7rR4b5dDPmJAoW9vb0dOTg6Ki4vx9NNPu7VOTEwMwsPD+12uKApiY2MHUo60FEXBI7eO9XYZ5GMGdHifm5uL7OxsZGVlub1OWloaDAYDjEYjNm7c2Gt5e3s7xowZg/j4eNx1112orq6+7PZUVYXNZnOZaPCcaVfxTrUFWz8/BSGEt8uhQeTxnr6srAwmkwlVVVVu9TcYDCgqKkJ6ejpUVcXq1athNBpRWVmJzMxMAEBSUhJeeeUVTJkyBTabDS+99BJmzpyJTz/9FBMmTOhzuwUFBcjPz/e0fHJDa4cdua+acLqjCwdOtOP/3jYRTxj7/n8g36MID17GzWYzpk+fjoqKCkydOhUAMGfOHEybNg2FhYVuP+j8+fOhKArWrl3b53KHw4Ebb7wRmZmZWLFiRZ99VFWFqqrOeZvNhoSEBFitVpdzB+S59/Y04fHXqjEmcgTMZzsQGxaErXlGb5dFfbDZbNDr9R497z3a05tMJrS0tCA9Pd3Z1tPTg82bN+Pll1+Gqqrw9/e/4nYyMjJQWlra73I/Pz/cdNNNaGho6LePVquFVqv1pHxyU2JUCADg6JkO3JQYAY2fH063q1AUBZEhgV6ujq6WR6E3Go2ora11aVu0aBGSkpKwbNkytwIPANXV1TAYDP0uF0KgpqYGU6ZM8aQ86Vg7uvCj1VXwUxT8aPY4zJ0U02e/vZZWxIQFIUbX99UQe7cDjWfOYXx0KBRFQcpoPf5wTwpWbzuKUWFBGBkSiPSnN8BPAZbdmYTHZo+/lsOia8yj0Ot0OqSkpLi0hYSEICoqytmel5cHi8WCkpISAEBhYSESExORnJwMu92O0tJSlJeXo7y83LmN/Px8ZGRkYMKECbDZbFixYgVqamrw5z//+WrHN6xVHmhBV48DwYF+eLfG0mfoD5yw4ffvfYawoAAUPzy9z+288GE93vnUgvwFKbgz5cIVlJxbxiDnljFoaevET1/bDX8F6BHAcx/U4cFbroMuKOCajo2unUH/am1zczMaGxud83a7HUuXLoXFYkFwcDCSk5Oxbt06zJs3z9mntbUVP/rRj3D8+HHo9XqkpaVh8+bNuPnmmwe7vGFlUqwOe5tsiNMHIS0hss8+o3TBiNZpMSZqRL/biY8Mxg2GMMTqex8JKFBQY7YiMlSLk20q/P35IU5f59GJvKFsICc0hoOm1vM41aZiSrweyjX6VN5LGxrwyecnsbuxFb+8Mwn/ljnumjwOeW4gz3uGntxy5pwKfz8/6IN5WD+UXPOz9ySvyBBeKRku+AaNSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDJXFfqCggIoioIlS5b026eyshKKovSa6urq+uxfVlYGRVFw9913X01pRNQPzUBX3LVrF4qKipCamupW//r6eoSFhTnno6Oje/U5evQoli5dilmzZg20LCK6ggHt6dvb25GTk4Pi4mJERES4tU5MTAxiY2Odk7+/v8vynp4e5OTkID8/H+PGjRtIWUTkhgGFPjc3F9nZ2cjKynJ7nbS0NBgMBhiNRmzcuLHX8t/97neIjo7GI4884tb2VFWFzWZzmYjoyjw+vC8rK4PJZEJVVZVb/Q0GA4qKipCeng5VVbF69WoYjUZUVlYiMzMTALBlyxasWrUKNTU1btdRUFCA/Px8T8snkp5HoTebzVi8eDEqKioQFBTk1jqTJk3CpEmTnPMzZsyA2WzG8uXLkZmZiba2Njz00EMoLi7GyJEj3a4lLy8PP//5z53zNpsNCQkJ7g+GSFIehd5kMqGlpQXp6enOtp6eHmzevBkvv/wyVFXt9V69LxkZGSgtLQUAHDx4EEeOHMH8+fOdyx0Ox4XiNBrU19dj/Pjxvbah1Wqh1Wo9KZ+I4GHojUYjamtrXdoWLVqEpKQkLFu2zK3AA0B1dTUMBgMAICkpqdc2f/3rX6OtrQ0vvfQS995Eg8yj0Ot0OqSkpLi0hYSEICoqytmel5cHi8WCkpISAEBhYSESExORnJwMu92O0tJSlJeXo7y8HAAQFBTUa5vh4eEA0KudiK7egK/T96e5uRmNjY3OebvdjqVLl8JisSA4OBjJyclYt24d5s2bN9gPTURuUIQQwttFDAabzQa9Xg+r1eryISCi4Wwgz3t+9p5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWQYeiLJMPREkmHoiSTD0BNJhqEnkgxDTyQZhp5IMgw9kWSuKvQFBQVQFAVLlizpt09lZSUURek11dXVOfusWbMG06dPR3h4OEJCQjBt2jSsXr36akojon5oBrrirl27UFRUhNTUVLf619fXIywszDkfHR3t/DsyMhK/+tWvkJSUhMDAQLz33ntYtGgRYmJicMcddwy0RCLqw4D29O3t7cjJyUFxcTEiIiLcWicmJgaxsbHOyd/f37lszpw5uOeeezB58mSMHz8eixcvRmpqKj755JOBlEdElzGg0Ofm5iI7OxtZWVlur5OWlgaDwQCj0YiNGzf2208IgY8++gj19fXIzMzst5+qqrDZbC4TEV2Zx4f3ZWVlMJlMqKqqcqu/wWBAUVER0tPToaoqVq9eDaPRiMrKSpdQW61WjB49Gqqqwt/fHytXrsRtt93W73YLCgqQn5/vaflE0lOEEMLdzmazGdOnT0dFRQWmTp0K4MKh+bRp01BYWOj2g86fPx+KomDt2rXONofDgUOHDqG9vR0fffQRfv/73+Odd97BnDlz+tyGqqpQVdU5b7PZkJCQAKvV6nLugGg4s9ls0Ov1Hj3vPQr9O++8g3vuucfl/XhPTw8URYGfn59zL30lf/jDH1BaWorPPvus3z6PPvoozGYzPvzwQ7dqG8jgiXzdQJ73Hh3eG41G1NbWurQtWrQISUlJWLZsmVuBB4Dq6moYDIbL9hFCuOzJiWhweBR6nU6HlJQUl7aQkBBERUU52/Py8mCxWFBSUgIAKCwsRGJiIpKTk2G321FaWory8nKUl5c7t1FQUIDp06dj/PjxsNvteP/991FSUoK//OUvVzs+IvqKAV+n709zczMaGxud83a7HUuXLoXFYkFwcDCSk5Oxbt06zJs3z9nn3Llz+MlPfoJjx44hODgYSUlJKC0txf333z/Y5RFJz6P39EMZ39OTjAbyvOdn74kkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSeaqQl9QUABFUbBkyZJ++1RWVkJRlF5TXV2ds09xcTFmzZqFiIgIREREICsrCzt37rya0oioH5qBrrhr1y4UFRUhNTXVrf719fUICwtzzkdHRzv/rqysxAMPPIBvfOMbCAoKwnPPPYfbb78d+/btw+jRowdaIhH1YUB7+vb2duTk5KC4uBgRERFurRMTE4PY2Fjn5O/v71z26quv4ic/+QmmTZuGpKQkFBcXw+Fw4KOPPhpIeUR0GQMKfW5uLrKzs5GVleX2OmlpaTAYDDAajdi4ceNl+3Z0dKCrqwuRkZH99lFVFTabzWUioivz+PC+rKwMJpMJVVVVbvU3GAwoKipCeno6VFXF6tWrYTQaUVlZiczMzD7X+eUvf4nRo0df9kWloKAA+fn5npZPRMIDjY2NIiYmRtTU1DjbZs+eLRYvXuzJZsRdd90l5s+f3+eyP/7xjyIiIkJ8+umnl91GZ2ensFqtzslsNgsAwmq1elQLkS+zWq0eP+892tObTCa0tLQgPT3d2dbT04PNmzfj5ZdfhqqqLu/V+5ORkYHS0tJe7cuXL8czzzyDDRs2XPEEoVarhVar9aR8IoKHh/dGoxG1tbUubYsWLUJSUhKWLVvmVuABoLq6GgaDwaXt+eefx9NPP40PP/wQ06dP96QsIvKAR6HX6XRISUlxaQsJCUFUVJSzPS8vDxaLBSUlJQCAwsJCJCYmIjk5GXa7HaWlpSgvL0d5eblzG8899xx+85vf4LXXXkNiYiKOHz8OAAgNDUVoaOhVDZCIXA34On1/mpub0djY6Jy32+1YunQpLBYLgoODkZycjHXr1mHevHnOPitXroTdbsd3vvMdl209+eSTeOqppwa7RCKpKUII4e0iBoPNZoNer4fVanX5EBDRcDaQ5z0/e08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDJXFfqCggIoioIlS5b026eyshKKovSa6urqnH327duHb3/720hMTISiKCgsLLyasojoMjQDXXHXrl0oKipCamqqW/3r6+sRFhbmnI+Ojnb+3dHRgXHjxuG+++7Dz372s4GWRERuGNCevr29HTk5OSguLkZERIRb68TExCA2NtY5+fv7O5fddNNNeP7557Fw4UJotdqBlEREbhpQ6HNzc5GdnY2srCy310lLS4PBYIDRaMTGjRsH8rAuVFWFzWZzmYjoyjw+vC8rK4PJZEJVVZVb/Q0GA4qKipCeng5VVbF69WoYjUZUVlYiMzPT44IvKigoQH5+/oDXJ5KVR6E3m81YvHgxKioqEBQU5NY6kyZNwqRJk5zzM2bMgNlsxvLly68q9Hl5efj5z3/unLfZbEhISBjw9ohk4dHhvclkQktLC9LT06HRaKDRaLBp0yasWLECGo0GPT09bm0nIyMDDQ0NAyr4Iq1Wi7CwMJeJiK7Moz290WhEbW2tS9uiRYuQlJSEZcuWuZycu5zq6moYDAZPHpqIBolHodfpdEhJSXFpCwkJQVRUlLM9Ly8PFosFJSUlAIDCwkIkJiYiOTkZdrsdpaWlKC8vR3l5uXMbdrsd+/fvd/5tsVhQU1OD0NBQXH/99Vc1QCJyNeDr9P1pbm5GY2Ojc95ut2Pp0qWwWCwIDg5GcnIy1q1bh3nz5jn7NDU1IS0tzTm/fPlyLF++HLNnz0ZlZeVgl0gkNUUIIbxdxGCw2WzQ6/WwWq18f0/SGMjznp+9J5IMQ08kGYaeSDIMPZFkGHoiyTD0RJJh6L8GVUfO4NDJdm+XQQSAob/mGk934KUNB/DU2n3eLoUIwDX4RB65ignTYly0Dvpg/lPT0MBn4jUWFOCP/G8le7sMIice3g8BH9edwCtbD+PsObu3SyEJMPReZjp6Fj98pQqrtx3FL/65x9vlkAQYei+zne+CTqtBcIA/BIbFd59oiGPovWz2xGj8KHMcpl0XjmV3Jnm7nGHp4Ml2tKvdAAAhBD4/0Ybzdte7PDW1nkdnl2ubtcPeqw0A1O4eOByOa1fwNcbQe5mfn4InjBPw9N1TEB8RjIMn29HXt50dDoFT7aoXKhyYHofAp+ZWmM909FrW1tmFQyfb4XBc3ZFNV48DpduPYnfjWWebwyHw1Np9+MeWwwCA2mNW5L66G4+tvnAj1yfX7sP/K9+DrBc3wdbZBQB4oaIej/z3LnxzeSVOf/Fv/H5tM+546V/IKPgIB060Obf/dvUxpDz5IX74ShVOtvX9/7Gl4SQ2HzjpnG9qPY91e5rQ1eNw1r2pvgVtXzz+pVo77Kg6cqbP58BgYeiHiNYOO7Je2ATjC5vw2/9xvaZv73bgnpVbMP3pDchbU9vn+g0n2vDnjZ/j3Bd7tEuZz3TgzSoz7N2ue6djZzvwcd2Jyz7B9hw7C0vr+X6Xd/U4UH/c5rKN7h4Hcv62Hd/68xZkPr8RpduPOpd90nAKGc98hCVlNbjvv7b12uMCgK2zC/+1+SDerDJftjbTkbP4YO9xlGw94mw7c86OPcdasWH/CTgcDozUBeIGQxhmXT8SAHCqTUXEiEDEhQfh8xMXPjDVcKINhrBgxEeMQO0xKwBgx6HTGD8yBJNG6bD90Gnn9jcfOIUbr4uAxl/BtoNftl+kdvWgcMMBrPiowRnqv1R+jreqzNiw/wQA4O3dx/DKliP466aDvdb/w7rPkLemFls+P9XvuK8WL9kNEXXH2xAXHoy48GBYWl33jgdPtkMb4I9bxkbi6Olzfa7/6o5G7G+yIj4iGN+aNtpl2d/+dQh7m6wI1fpj3pQ4Z/t/rm/AybbzCArwxzfGj+y1zcMn2/H7dz/DCK0//vuHt/T5uC9//DlMR8/i3htH494b4wEA2w+dxoET7YgcEQgHBP7+yWE8lDEGAPDnjZ9DF6RBoL+C+uM2fLCvGfekxbts8/fv7se/Gk7iuE1FSKA/slPjej0uAEy7Lhw3xIVhxrgoZ9tInRaLsyYiVOsPPz8/GPTBePH+ac7laWMisH7fCWg1/kgZrQcA3Dw2Cu/UWGDvduDGxAs/3jJjfBQ21regrbMbfxj/5S3ikuPC8O6nTUiN1+PGMeG9atIG+GPmhGio3Q6Eai/EK3NiDP7VcBJp113Y9pT4cOw6chbpY3r/UExyXBgaWtoxNjq0zzEPBoZ+iJiWEA4/PwXVjWd7vbe/PiYUPQ4B09Gz+P6MMX2uv2DahWDMmhDda9mcSTE409GF6WMiXdpvHhuB6kYFk0bp+txmTFgQRkdceCHqT0LkCOxtsmL0JX1CtRqcueTyY+SIQOffuiANjttURIUG4oY4Pfz9eh9s9giBEYEa6IMdON/V/3vnoAB//Me8yb3aZ0/s/W9w0b/NGodHZibC75LH/eGtY/HgzQnQBvhDURQAwJ0pBtx6/UgEaPyg1Xx5w9dHZ43DXalxCAvWYERg3/FZkjXRZf62G0bhthtGOecnG8Lw/H1T+1z3BzPH4gczx/Zb/2Dg7bKGGIdDwM9P6dUuhIDtfDf0IwK8UJVnhBB47sN6rPrXYUSGBOI/75+GGeOj0Gw9jxkFHwO4sEfb32RD1a+zEBXq+lNmx62dWPXJIYRqA5A7dzw0/nwX2p+BPO+5px9i+go8ACiK4hOBBy7UuuzOJCy9fRL8FDj3ngF+fkiICIYDgBBAgMYPwYG9b5seqw/Cr7Jv+JqrlgdfQuma8fdTnIEHLrzffvgbiTjZpiJAo+D330ru9xCZrh0e3tPXTu3ugSKAwAD3fhyF+se74ZJP0Gr8GfgreGXLYfzpo4Zrcr2eoScaYtSuHvzv3uN4Yf0BnOvjcwxXi2+oiIYYbYA/fjx7PB65dazzWv9gYuiJhqA5STHXbNs8vCeSDENPJBmGnkgyDD2RZBh6Iskw9ESSYeiJJMPQE0mGoSeSDENPJBmGnkgyDD2RZBh6Iskw9ESSGTZfrb14hxGbzeblSoi+Phef757cYWfYhL6t7cJPDyUkJHi5EqKvX1tbG/R6vVt9h82NMR0OB5qamqDT6VzuwGqz2ZCQkACz2TxsbpjJMfmGr2NMQgi0tbUhLi7O5Qc8LmfY7On9/PwQHx/f7/KwsLBh82S6iGPyDdd6TO7u4S/iiTwiyTD0RJIZ9qHXarV48sknodVqr9zZR3BMvmGojmnYnMgjIvcM+z09Ebli6Ikkw9ATSYahJ5LMkAp9QUEBFEXBkiVLnG1CCDz11FOIi4tDcHAw5syZg3379rmsp6oqnnjiCYwcORIhISFYsGABjh075tLn7Nmz+N73vge9Xg+9Xo/vfe97aG1tdenT2NiI+fPnIyQkBCNHjsRPf/pT2O12lz61tbWYPXs2goODMXr0aPzud7+77Oee+xrTmjVrcMcdd2DkyJFQFAU1NTW91vOlMXV1dWHZsmWYMmUKQkJCEBcXh+9///toamrymTH1NS4AeOqpp5CUlISQkBBEREQgKysLO3bs8Klx9SKGiJ07d4rExESRmpoqFi9e7Gx/9tlnhU6nE+Xl5aK2tlbcf//9wmAwCJvN5uzz7//+72L06NFi/fr1Yvfu3WLu3Lli6tSporu729nnzjvvFCkpKWLr1q1i69atIiUlRdx1113O5d3d3SIlJUXMnTtX7N69W6xfv17ExcWJxx9/3NnHarWKUaNGiYULF4ra2lpRXl4udDqdWL58uUdjKikpEfn5+aK4uFgAENXV1b3W9aUxtba2iqysLPHGG2+Iuro6sW3bNnHLLbeI9PR0nxjT5f6vXn31VbF+/Xpx8OBBsXfvXvHII4+IsLAw0dLS4hPj6suQCH1bW5uYMGGCWL9+vZg9e7bzH93hcIjY2Fjx7LPPOvt2dnYKvV4v/vrXvwohLjzhAgICRFlZmbOPxWIRfn5+4oMPPhBCCLF//34BQGzfvt3ZZ9u2bQKAqKurE0II8f777ws/Pz9hsVicfV5//XWh1WqF1WoVQgixcuVKodfrRWdnp7NPQUGBiIuLEw6Hw60xXerw4cN9ht6Xx3TRzp07BQBx9OjRIT0mT8dltVoFALFhw4YhP67+DInD+9zcXGRnZyMrK8ul/fDhwzh+/Dhuv/12Z5tWq8Xs2bOxdetWAIDJZEJXV5dLn7i4OKSkpDj7bNu2DXq9HrfccouzT0ZGBvR6vUuflJQUxMXFOfvccccdUFUVJpPJ2Wf27NkuH7a444470NTUhCNHjrg1JncMhzFZrVYoioLw8PAhPSZPxmW321FUVAS9Xo+pU6cO+XH1x+tfuCkrK4PJZEJVVVWvZcePHwcAjBo1yqV91KhROHr0qLNPYGAgIiIievW5uP7x48cRE9P7p39jYmJc+nz1cSIiIhAYGOjSJzExsdfjXFw2duzYK47JHb4+ps7OTvzyl7/Egw8+6PyiyVAck7vjeu+997Bw4UJ0dHTAYDBg/fr1GDly5JAe1+V4NfRmsxmLFy9GRUUFgoKC+u136VdlgQsn977a9lVf7dNX/8HoI744iXKx3d0xDYQvjKmrqwsLFy6Ew+HAypUrh+yYAPfHNXfuXNTU1ODUqVMoLi7Gd7/7XezYsaPPIA+FcV2JVw/vTSYTWlpakJ6eDo1GA41Gg02bNmHFihXQaDQur2KXamlpcS6LjY2F3W7H2bNnL9vnxIkTvR7/5MmTLn2++jhnz55FV1fXZfu0tLQA+PIV90pj6unpueK/i6+OqaurC9/97ndx+PBhrF+/3uXrpENtTJ6MKyQkBNdffz0yMjKwatUqaDQarFq1asiO64rcfvd/DdhsNlFbW+syTZ8+XTz00EOitrbWeSLvj3/8o3MdVVX7PJH3xhtvOPs0NTX1eSJlx44dzj7bt2/v80RKU1OTs09ZWVmvEynh4eFCVVVnn2effdblRMqVxnSpK53I86Ux2e12cffdd4vk5GSXM9tDdUye/l9davz48eLJJ58csuO6kiFx9v5SXz17+uyzzwq9Xi/WrFkjamtrxQMPPNDnJbv4+HixYcMGsXv3bvHNb36zz0smqampYtu2bWLbtm1iypQpfV4yMRqNYvfu3WLDhg0iPj7e5ZJJa2urGDVqlHjggQdEbW2tWLNmjQgLC7viJZOvjun06dOiurparFu3TgAQZWVlorq6WjQ3N/vkmLq6usSCBQtEfHy8qKmpEc3Nzc7p0ifoUB/TV8fV3t4u8vLyxLZt28SRI0eEyWQSjzzyiNBqtWLv3r0+Na5LDfnQOxwO8eSTT4rY2Fih1WpFZmZmr1fh8+fPi8cff1xERkaK4OBgcdddd4nGxkaXPqdPnxY5OTlCp9MJnU4ncnJyxNmzZ136HD16VGRnZ4vg4GARGRkpHn/8cZfLI0IIsWfPHjFr1iyh1WpFbGyseOqpp674KvvVMf3jH/8QAHpNF/cevjami0csfU0bN270mTF9dVznz58X99xzj4iLixOBgYHCYDCIBQsWiJ07d7qs4wvjuhS/WkskmSFxnZ6Ivj4MPZFkGHoiyTD0RJJh6Ikkw9ATSYahJ5IMQ08kGYaeSDIMPZFkGHoiyTD0RJL5/wcVXJ3hGAhIAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gdf_buffered_stops_with_centers_taz = gpd.read_file(\"intermediate/buffered-stops-with-centers-taz.geojson\")\n",
    "display(gdf_buffered_stops_with_centers_taz)\n",
    "gdf_buffered_stops_with_centers_taz.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>TAZID</th>\n",
       "      <th>pctu_ch</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>10036</td>\n",
       "      <td>3.0</td>\n",
       "      <td>657.0</td>\n",
       "      <td>0.918003</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>10036</td>\n",
       "      <td>3.0</td>\n",
       "      <td>658.0</td>\n",
       "      <td>0.733657</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>10036</td>\n",
       "      <td>3.0</td>\n",
       "      <td>663.0</td>\n",
       "      <td>0.991408</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>10036</td>\n",
       "      <td>3.0</td>\n",
       "      <td>664.0</td>\n",
       "      <td>0.194214</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>10042</td>\n",
       "      <td>3.0</td>\n",
       "      <td>497.0</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>166</th>\n",
       "      <td>27931</td>\n",
       "      <td>3.0</td>\n",
       "      <td>508.0</td>\n",
       "      <td>0.054408</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>167</th>\n",
       "      <td>28117</td>\n",
       "      <td>3.0</td>\n",
       "      <td>468.0</td>\n",
       "      <td>0.572808</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>168</th>\n",
       "      <td>28117</td>\n",
       "      <td>3.0</td>\n",
       "      <td>473.0</td>\n",
       "      <td>0.996042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>169</th>\n",
       "      <td>28117</td>\n",
       "      <td>3.0</td>\n",
       "      <td>497.0</td>\n",
       "      <td>0.791458</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>170</th>\n",
       "      <td>28117</td>\n",
       "      <td>3.0</td>\n",
       "      <td>500.0</td>\n",
       "      <td>0.702747</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>171 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         N  project_id  TAZID   pctu_ch\n",
       "0    10036         3.0  657.0  0.918003\n",
       "1    10036         3.0  658.0  0.733657\n",
       "2    10036         3.0  663.0  0.991408\n",
       "3    10036         3.0  664.0  0.194214\n",
       "4    10042         3.0  497.0  0.000000\n",
       "..     ...         ...    ...       ...\n",
       "166  27931         3.0  508.0  0.054408\n",
       "167  28117         3.0  468.0  0.572808\n",
       "168  28117         3.0  473.0  0.996042\n",
       "169  28117         3.0  497.0  0.791458\n",
       "170  28117         3.0  500.0  0.702747\n",
       "\n",
       "[171 rows x 4 columns]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>TAZID</th>\n",
       "      <th>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>TOTEMP</th>\n",
       "      <th>TOTHH</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>165.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.991571</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>165.0</td>\n",
       "      <td>15191.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.008429</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>172.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.971864</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>172.0</td>\n",
       "      <td>15185.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.028136</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>173.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.996895</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>318</th>\n",
       "      <td>1530.0</td>\n",
       "      <td>25381.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000236</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>319</th>\n",
       "      <td>1531.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.601742</td>\n",
       "      <td>0.974230</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>320</th>\n",
       "      <td>1531.0</td>\n",
       "      <td>25381.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0.398258</td>\n",
       "      <td>0.025770</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>321</th>\n",
       "      <td>1679.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.999530</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>322</th>\n",
       "      <td>1679.0</td>\n",
       "      <td>25381.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000470</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>323 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      TAZID        N  project_id    TOTEMP     TOTHH\n",
       "0     165.0     -1.0        -1.0  0.991571  1.000000\n",
       "1     165.0  15191.0         1.0  0.008429  0.000000\n",
       "2     172.0     -1.0        -1.0  0.000000  0.971864\n",
       "3     172.0  15185.0         1.0  0.000000  0.028136\n",
       "4     173.0     -1.0        -1.0  0.000000  0.996895\n",
       "..      ...      ...         ...       ...       ...\n",
       "318  1530.0  25381.0         2.0  0.000000  0.000236\n",
       "319  1531.0     -1.0        -1.0  0.601742  0.974230\n",
       "320  1531.0  25381.0         2.0  0.398258  0.025770\n",
       "321  1679.0     -1.0        -1.0  0.000000  0.999530\n",
       "322  1679.0  25381.0         2.0  0.000000  0.000470\n",
       "\n",
       "[323 rows x 5 columns]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# read proportions from REMM data\n",
    "df_percent_changeable = pd.read_csv('intermediate/percent-changeable.csv')\n",
    "display(df_percent_changeable)\n",
    "df_area_ratio_of_taz = pd.read_csv('intermediate/parameter-area-ratio-of-taz.csv')\n",
    "display(df_area_ratio_of_taz)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "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>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>nearest_AreaType</th>\n",
       "      <th>TAZID</th>\n",
       "      <th>TAZACRES</th>\n",
       "      <th>sqm</th>\n",
       "      <th>acres</th>\n",
       "      <th>geometry</th>\n",
       "      <th>pctu_ch</th>\n",
       "      <th>parameter</th>\n",
       "      <th>taz_ratio</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1472</td>\n",
       "      <td>42.036264</td>\n",
       "      <td>49740.404222</td>\n",
       "      <td>12.291122</td>\n",
       "      <td>POLYGON ((424043.145 4505847.017, 424036.495 4...</td>\n",
       "      <td>0.776678</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.238638</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1473</td>\n",
       "      <td>44.089344</td>\n",
       "      <td>52165.559080</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>POLYGON ((424421.505 4505865.604, 424414.855 4...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.327186</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1481</td>\n",
       "      <td>68.022708</td>\n",
       "      <td>13056.790069</td>\n",
       "      <td>3.226403</td>\n",
       "      <td>POLYGON ((424229.000 4505722.832, 424209.282 4...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.041407</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1482</td>\n",
       "      <td>39.922698</td>\n",
       "      <td>11968.229688</td>\n",
       "      <td>2.957414</td>\n",
       "      <td>POLYGON ((424396.265 4505812.237, 424384.505 4...</td>\n",
       "      <td>0.567510</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.039493</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>15183</td>\n",
       "      <td>1</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>176</td>\n",
       "      <td>288.616036</td>\n",
       "      <td>35706.046305</td>\n",
       "      <td>8.823156</td>\n",
       "      <td>POLYGON ((412496.126 4566424.387, 412502.776 4...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.104515</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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>351</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1383</td>\n",
       "      <td>82.068433</td>\n",
       "      <td>163541.758588</td>\n",
       "      <td>40.412049</td>\n",
       "      <td>POLYGON ((418599.665 4505156.000, 418601.602 4...</td>\n",
       "      <td>0.289566</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.355068</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>352</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1383</td>\n",
       "      <td>82.068433</td>\n",
       "      <td>10615.218005</td>\n",
       "      <td>2.623077</td>\n",
       "      <td>POLYGON ((418812.341 4505510.828, 418778.475 4...</td>\n",
       "      <td>0.603012</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.039094</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>353</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.628618</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>354</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.060703</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>355</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1396</td>\n",
       "      <td>98.491949</td>\n",
       "      <td>129814.582427</td>\n",
       "      <td>32.077882</td>\n",
       "      <td>POLYGON ((419041.436 4504755.602, 419002.000 4...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.333237</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>356 rows × 11 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         N  project_id     nearest_AreaType  TAZID    TAZACRES            sqm  \\\n",
       "0    15086           2  Neighborhood Center   1472   42.036264   49740.404222   \n",
       "1    15086           2  Neighborhood Center   1473   44.089344   52165.559080   \n",
       "2    15086           2  Neighborhood Center   1481   68.022708   13056.790069   \n",
       "3    15086           2  Neighborhood Center   1482   39.922698   11968.229688   \n",
       "4    15183           1  Neighborhood Center    176  288.616036   35706.046305   \n",
       "..     ...         ...                  ...    ...         ...            ...   \n",
       "351  15069           2         Urban Center   1383   82.068433  163541.758588   \n",
       "352  23506           2         Urban Center   1383   82.068433   10615.218005   \n",
       "353  15069           2         Urban Center   1384   96.889876  208196.248028   \n",
       "354  23506           2         Urban Center   1384   96.889876   34533.423293   \n",
       "355  15069           2         Urban Center   1396   98.491949  129814.582427   \n",
       "\n",
       "         acres                                           geometry   pctu_ch  \\\n",
       "0    12.291122  POLYGON ((424043.145 4505847.017, 424036.495 4...  0.776678   \n",
       "1    12.890390  POLYGON ((424421.505 4505865.604, 424414.855 4...  1.000000   \n",
       "2     3.226403  POLYGON ((424229.000 4505722.832, 424209.282 4...  1.000000   \n",
       "3     2.957414  POLYGON ((424396.265 4505812.237, 424384.505 4...  0.567510   \n",
       "4     8.823156  POLYGON ((412496.126 4566424.387, 412502.776 4...  0.000000   \n",
       "..         ...                                                ...       ...   \n",
       "351  40.412049  POLYGON ((418599.665 4505156.000, 418601.602 4...  0.289566   \n",
       "352   2.623077  POLYGON ((418812.341 4505510.828, 418778.475 4...  0.603012   \n",
       "353  51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "354   8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "355  32.077882  POLYGON ((419041.436 4504755.602, 419002.000 4...  0.000000   \n",
       "\n",
       "    parameter  taz_ratio  \n",
       "0      TOTEMP   0.238638  \n",
       "1      TOTEMP   0.327186  \n",
       "2      TOTEMP   0.041407  \n",
       "3      TOTEMP   0.039493  \n",
       "4      TOTEMP   0.104515  \n",
       "..        ...        ...  \n",
       "351     TOTHH   0.355068  \n",
       "352     TOTHH   0.039094  \n",
       "353     TOTHH   0.628618  \n",
       "354     TOTHH   0.060703  \n",
       "355     TOTHH   0.333237  \n",
       "\n",
       "[356 rows x 11 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# join tables\n",
    "gdf_buffered_stops_with_centers_taz_changeable = pd.merge(gdf_buffered_stops_with_centers_taz, df_percent_changeable, on=['N','project_id','TAZID'], how='left')\n",
    "gdf_buffered_stops_with_centers_taz_changeable_ratios = pd.merge(gdf_buffered_stops_with_centers_taz_changeable, df_area_ratio_of_taz, on=['N','project_id','TAZID'], how='left')\n",
    "gdf_buffered_stops_with_centers_taz_changeable_ratios.fillna(0, inplace=True)\n",
    "gdf_buffered_stops_with_centers_taz_changeable_ratios = gdf_buffered_stops_with_centers_taz_changeable_ratios.melt(id_vars=['N','project_id','nearest_AreaType','TAZID','TAZACRES','sqm','acres','geometry','pctu_ch'],var_name='parameter',value_name='taz_ratio')\n",
    "gdf_buffered_stops_with_centers_taz_changeable_ratios"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Prep socioeconomic data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>TAZID</th>\n",
       "      <th>parameter</th>\n",
       "      <th>value_B</th>\n",
       "      <th>se_year_B</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30001.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30002.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30003.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30004.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30005.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92191</th>\n",
       "      <td>3542</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92192</th>\n",
       "      <td>3543</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92193</th>\n",
       "      <td>3544</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92194</th>\n",
       "      <td>3545</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92195</th>\n",
       "      <td>3546</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>92196 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       TAZID parameter  value_B  se_year_B\n",
       "0          1  CO_TAZID  30001.0       2023\n",
       "1          2  CO_TAZID  30002.0       2023\n",
       "2          3  CO_TAZID  30003.0       2023\n",
       "3          4  CO_TAZID  30004.0       2023\n",
       "4          5  CO_TAZID  30005.0       2023\n",
       "...      ...       ...      ...        ...\n",
       "92191   3542   CO_FIPS     49.0       2023\n",
       "92192   3543   CO_FIPS     49.0       2023\n",
       "92193   3544   CO_FIPS     49.0       2023\n",
       "92194   3545   CO_FIPS     49.0       2023\n",
       "92195   3546   CO_FIPS     49.0       2023\n",
       "\n",
       "[92196 rows x 4 columns]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "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>TAZID</th>\n",
       "      <th>parameter</th>\n",
       "      <th>value_F</th>\n",
       "      <th>se_year_F</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30001.0</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30002.0</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30003.0</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30004.0</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>CO_TAZID</td>\n",
       "      <td>30005.0</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92191</th>\n",
       "      <td>3542</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92192</th>\n",
       "      <td>3543</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92193</th>\n",
       "      <td>3544</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92194</th>\n",
       "      <td>3545</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>92195</th>\n",
       "      <td>3546</td>\n",
       "      <td>CO_FIPS</td>\n",
       "      <td>49.0</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>276588 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       TAZID parameter  value_F  se_year_F\n",
       "0          1  CO_TAZID  30001.0       2032\n",
       "1          2  CO_TAZID  30002.0       2032\n",
       "2          3  CO_TAZID  30003.0       2032\n",
       "3          4  CO_TAZID  30004.0       2032\n",
       "4          5  CO_TAZID  30005.0       2032\n",
       "...      ...       ...      ...        ...\n",
       "92191   3542   CO_FIPS     49.0       2050\n",
       "92192   3543   CO_FIPS     49.0       2050\n",
       "92193   3544   CO_FIPS     49.0       2050\n",
       "92194   3545   CO_FIPS     49.0       2050\n",
       "92195   3546   CO_FIPS     49.0       2050\n",
       "\n",
       "[276588 rows x 4 columns]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# get SE data\n",
    "df_se_data = pd.DataFrame()\n",
    "\n",
    "for se_year in se_years:\n",
    "    file = \"SE_\" + str(se_year) + '.csv'\n",
    "    se_file = os.path.join(input_model_se_folder, file)\n",
    "\n",
    "    _dbf = pd.read_csv(se_file)\n",
    "\n",
    "    # Identify non-numeric columns\n",
    "    non_numeric_columns = _dbf.select_dtypes(include=[int, float]).columns.tolist()\n",
    "    non_numeric_columns.remove(';TAZID')  # Keep 'TAZID' as identifier variable\n",
    "\n",
    "    _df_pivot = _dbf.melt(id_vars=[';TAZID'], value_vars=non_numeric_columns)\n",
    "\n",
    "    _df_pivot['se_year'] = se_year\n",
    "\n",
    "    if se_year==se_base_year:\n",
    "        df_se_data_base = _df_pivot\n",
    "    else:\n",
    "        df_se_data = pd.concat([df_se_data,_df_pivot])\n",
    "\n",
    "df_se_data_base.columns = df_se_data_base.columns + '_B'\n",
    "df_se_data.columns = df_se_data.columns + '_F'\n",
    "\n",
    "df_se_data_base.rename(columns={';TAZID_B':'TAZID','variable_B':'parameter'}, inplace=True)\n",
    "df_se_data.rename(columns={';TAZID_F':'TAZID','variable_F':'parameter'}, inplace=True)\n",
    "\n",
    "display(df_se_data_base)\n",
    "display(df_se_data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>TAZID</th>\n",
       "      <th>parameter</th>\n",
       "      <th>value_F</th>\n",
       "      <th>se_year_F</th>\n",
       "      <th>value_B</th>\n",
       "      <th>se_year_B</th>\n",
       "      <th>growth</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>10638</th>\n",
       "      <td>1</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2032</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10639</th>\n",
       "      <td>1</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2042</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10640</th>\n",
       "      <td>1</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2050</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10641</th>\n",
       "      <td>2</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>17.4</td>\n",
       "      <td>2032</td>\n",
       "      <td>13.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>4.4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10642</th>\n",
       "      <td>2</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>29.8</td>\n",
       "      <td>2042</td>\n",
       "      <td>13.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>16.8</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",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53185</th>\n",
       "      <td>3545</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2042</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53186</th>\n",
       "      <td>3545</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2050</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53187</th>\n",
       "      <td>3546</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2032</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53188</th>\n",
       "      <td>3546</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2042</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53189</th>\n",
       "      <td>3546</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2050</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2023</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>21276 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       TAZID parameter  value_F  se_year_F  value_B  se_year_B  growth\n",
       "10638      1     TOTHH      0.0       2032      0.0       2023     0.0\n",
       "10639      1     TOTHH      0.0       2042      0.0       2023     0.0\n",
       "10640      1     TOTHH      0.0       2050      0.0       2023     0.0\n",
       "10641      2     TOTHH     17.4       2032     13.0       2023     4.4\n",
       "10642      2     TOTHH     29.8       2042     13.0       2023    16.8\n",
       "...      ...       ...      ...        ...      ...        ...     ...\n",
       "53185   3545    TOTEMP      0.0       2042      0.0       2023     0.0\n",
       "53186   3545    TOTEMP      0.0       2050      0.0       2023     0.0\n",
       "53187   3546    TOTEMP      0.0       2032      0.0       2023     0.0\n",
       "53188   3546    TOTEMP      0.0       2042      0.0       2023     0.0\n",
       "53189   3546    TOTEMP      0.0       2050      0.0       2023     0.0\n",
       "\n",
       "[21276 rows x 7 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_se_growth = pd.merge(df_se_data, df_se_data_base, on=['TAZID','parameter'])\n",
    "df_se_growth['growth'] = df_se_growth['value_F'] - df_se_growth['value_B']\n",
    "df_se_growth_hhemp = df_se_growth[df_se_growth['parameter'].isin(['TOTHH','TOTEMP'])]\n",
    "df_se_growth_hhemp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_se_growth.to_csv('intermediate/se_growth.csv', index=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get max densities by center type"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>center_type</th>\n",
       "      <th>parameter</th>\n",
       "      <th>perac_ch</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>18</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>City Center</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>39</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Urban Center</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>59</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>91</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>City Center</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Urban Center</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>20</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>30</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           center_type parameter  perac_ch\n",
       "0  Neighborhood Center     TOTHH        18\n",
       "1          City Center     TOTHH        39\n",
       "2         Urban Center     TOTHH        59\n",
       "3  Metropolitan Center     TOTHH        91\n",
       "4  Neighborhood Center    TOTEMP         6\n",
       "5          City Center    TOTEMP        13\n",
       "6         Urban Center    TOTEMP        20\n",
       "7  Metropolitan Center    TOTEMP        30"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# units are TOTHH and emps are TOTEMP\n",
    "_df = df_centers.copy()\n",
    "_df.rename(columns={'final_units_per_acre':'TOTHH','final_emps_per_acre':'TOTEMP'}, inplace=True)\n",
    "df_centers_max_per_acre = _df.melt(id_vars='center_type',var_name='parameter', value_vars=['TOTHH','TOTEMP'], value_name='perac_ch')\n",
    "df_centers_max_per_acre['perac_ch'] = df_centers_max_per_acre['perac_ch'].astype(int)\n",
    "df_centers_max_per_acre"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculate additional hh/emp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_stop_taz_calcs = pd.merge(gdf_buffered_stops_with_centers_taz_changeable_ratios, df_centers_max_per_acre, left_on=('nearest_AreaType', 'parameter'), right_on=('center_type','parameter'))\n",
    "gdf_stop_taz_calcs_se = pd.merge(gdf_stop_taz_calcs, df_se_growth, on=['TAZID','parameter'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>nearest_AreaType</th>\n",
       "      <th>TAZID</th>\n",
       "      <th>TAZACRES</th>\n",
       "      <th>sqm</th>\n",
       "      <th>acres</th>\n",
       "      <th>geometry</th>\n",
       "      <th>pctu_ch</th>\n",
       "      <th>parameter</th>\n",
       "      <th>...</th>\n",
       "      <th>acres_ch</th>\n",
       "      <th>s_ch</th>\n",
       "      <th>s_tot</th>\n",
       "      <th>s_diff</th>\n",
       "      <th>s_grow</th>\n",
       "      <th>b_tot_den</th>\n",
       "      <th>f_tot_den</th>\n",
       "      <th>s_tot_den</th>\n",
       "      <th>s_grow_den</th>\n",
       "      <th>s_diff_den</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1472</td>\n",
       "      <td>42.036264</td>\n",
       "      <td>49740.404222</td>\n",
       "      <td>12.291122</td>\n",
       "      <td>POLYGON ((424043.145 4505847.017, 424036.495 4...</td>\n",
       "      <td>0.776678</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>9.546246</td>\n",
       "      <td>57.277477</td>\n",
       "      <td>97.353816</td>\n",
       "      <td>-82.101744</td>\n",
       "      <td>-75.738072</td>\n",
       "      <td>14.082676</td>\n",
       "      <td>14.600422</td>\n",
       "      <td>7.920662</td>\n",
       "      <td>-6.162015</td>\n",
       "      <td>-6.679760</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1472</td>\n",
       "      <td>42.036264</td>\n",
       "      <td>49740.404222</td>\n",
       "      <td>12.291122</td>\n",
       "      <td>POLYGON ((424043.145 4505847.017, 424036.495 4...</td>\n",
       "      <td>0.776678</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>9.546246</td>\n",
       "      <td>57.277477</td>\n",
       "      <td>98.828256</td>\n",
       "      <td>-87.229614</td>\n",
       "      <td>-74.263632</td>\n",
       "      <td>14.082676</td>\n",
       "      <td>15.137583</td>\n",
       "      <td>8.040621</td>\n",
       "      <td>-6.042055</td>\n",
       "      <td>-7.096961</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1472</td>\n",
       "      <td>42.036264</td>\n",
       "      <td>49740.404222</td>\n",
       "      <td>12.291122</td>\n",
       "      <td>POLYGON ((424043.145 4505847.017, 424036.495 4...</td>\n",
       "      <td>0.776678</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>9.546246</td>\n",
       "      <td>57.277477</td>\n",
       "      <td>100.204991</td>\n",
       "      <td>-92.017686</td>\n",
       "      <td>-72.886896</td>\n",
       "      <td>14.082676</td>\n",
       "      <td>15.639149</td>\n",
       "      <td>8.152632</td>\n",
       "      <td>-5.930044</td>\n",
       "      <td>-7.486517</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1473</td>\n",
       "      <td>44.089344</td>\n",
       "      <td>52165.559080</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>POLYGON ((424421.505 4505865.604, 424414.855 4...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>77.342342</td>\n",
       "      <td>77.342342</td>\n",
       "      <td>-168.865158</td>\n",
       "      <td>-139.145758</td>\n",
       "      <td>16.794534</td>\n",
       "      <td>19.100081</td>\n",
       "      <td>6.000000</td>\n",
       "      <td>-10.794534</td>\n",
       "      <td>-13.100081</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>15086</td>\n",
       "      <td>2</td>\n",
       "      <td>Neighborhood Center</td>\n",
       "      <td>1473</td>\n",
       "      <td>44.089344</td>\n",
       "      <td>52165.559080</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>POLYGON ((424421.505 4505865.604, 424414.855 4...</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>12.890390</td>\n",
       "      <td>77.342342</td>\n",
       "      <td>77.342342</td>\n",
       "      <td>-206.764208</td>\n",
       "      <td>-139.145758</td>\n",
       "      <td>16.794534</td>\n",
       "      <td>22.040182</td>\n",
       "      <td>6.000000</td>\n",
       "      <td>-10.794534</td>\n",
       "      <td>-16.040182</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>1063</th>\n",
       "      <td>10046</td>\n",
       "      <td>1</td>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>390</td>\n",
       "      <td>37.583728</td>\n",
       "      <td>146896.322244</td>\n",
       "      <td>36.298872</td>\n",
       "      <td>POLYGON ((418169.004 4563560.006, 418123.854 4...</td>\n",
       "      <td>0.588706</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>21.369354</td>\n",
       "      <td>1944.611210</td>\n",
       "      <td>1996.818707</td>\n",
       "      <td>1869.884050</td>\n",
       "      <td>1908.744848</td>\n",
       "      <td>2.426353</td>\n",
       "      <td>3.496931</td>\n",
       "      <td>55.010490</td>\n",
       "      <td>52.584137</td>\n",
       "      <td>51.513558</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1064</th>\n",
       "      <td>10046</td>\n",
       "      <td>1</td>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>390</td>\n",
       "      <td>37.583728</td>\n",
       "      <td>146896.322244</td>\n",
       "      <td>36.298872</td>\n",
       "      <td>POLYGON ((418169.004 4563560.006, 418123.854 4...</td>\n",
       "      <td>0.588706</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>21.369354</td>\n",
       "      <td>1944.611210</td>\n",
       "      <td>2005.989409</td>\n",
       "      <td>1856.757573</td>\n",
       "      <td>1917.915550</td>\n",
       "      <td>2.426353</td>\n",
       "      <td>4.111198</td>\n",
       "      <td>55.263134</td>\n",
       "      <td>52.836781</td>\n",
       "      <td>51.151936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1065</th>\n",
       "      <td>10046</td>\n",
       "      <td>1</td>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>392</td>\n",
       "      <td>26.797133</td>\n",
       "      <td>55531.069672</td>\n",
       "      <td>13.722026</td>\n",
       "      <td>POLYGON ((418384.296 4563850.293, 418366.563 4...</td>\n",
       "      <td>0.842694</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>11.563467</td>\n",
       "      <td>1052.275512</td>\n",
       "      <td>1052.950453</td>\n",
       "      <td>1048.659832</td>\n",
       "      <td>1052.950453</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.312681</td>\n",
       "      <td>76.734328</td>\n",
       "      <td>76.734328</td>\n",
       "      <td>76.421646</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1066</th>\n",
       "      <td>10046</td>\n",
       "      <td>1</td>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>392</td>\n",
       "      <td>26.797133</td>\n",
       "      <td>55531.069672</td>\n",
       "      <td>13.722026</td>\n",
       "      <td>POLYGON ((418384.296 4563850.293, 418366.563 4...</td>\n",
       "      <td>0.842694</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>11.563467</td>\n",
       "      <td>1052.275512</td>\n",
       "      <td>1054.300336</td>\n",
       "      <td>1041.428471</td>\n",
       "      <td>1054.300336</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.938044</td>\n",
       "      <td>76.832701</td>\n",
       "      <td>76.832701</td>\n",
       "      <td>75.894657</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1067</th>\n",
       "      <td>10046</td>\n",
       "      <td>1</td>\n",
       "      <td>Metropolitan Center</td>\n",
       "      <td>392</td>\n",
       "      <td>26.797133</td>\n",
       "      <td>55531.069672</td>\n",
       "      <td>13.722026</td>\n",
       "      <td>POLYGON ((418384.296 4563850.293, 418366.563 4...</td>\n",
       "      <td>0.842694</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>11.563467</td>\n",
       "      <td>1052.275512</td>\n",
       "      <td>1056.601271</td>\n",
       "      <td>1029.102288</td>\n",
       "      <td>1056.601271</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.004003</td>\n",
       "      <td>77.000383</td>\n",
       "      <td>77.000383</td>\n",
       "      <td>74.996380</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1068 rows × 32 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          N  project_id     nearest_AreaType  TAZID   TAZACRES            sqm  \\\n",
       "0     15086           2  Neighborhood Center   1472  42.036264   49740.404222   \n",
       "1     15086           2  Neighborhood Center   1472  42.036264   49740.404222   \n",
       "2     15086           2  Neighborhood Center   1472  42.036264   49740.404222   \n",
       "3     15086           2  Neighborhood Center   1473  44.089344   52165.559080   \n",
       "4     15086           2  Neighborhood Center   1473  44.089344   52165.559080   \n",
       "...     ...         ...                  ...    ...        ...            ...   \n",
       "1063  10046           1  Metropolitan Center    390  37.583728  146896.322244   \n",
       "1064  10046           1  Metropolitan Center    390  37.583728  146896.322244   \n",
       "1065  10046           1  Metropolitan Center    392  26.797133   55531.069672   \n",
       "1066  10046           1  Metropolitan Center    392  26.797133   55531.069672   \n",
       "1067  10046           1  Metropolitan Center    392  26.797133   55531.069672   \n",
       "\n",
       "          acres                                           geometry   pctu_ch  \\\n",
       "0     12.291122  POLYGON ((424043.145 4505847.017, 424036.495 4...  0.776678   \n",
       "1     12.291122  POLYGON ((424043.145 4505847.017, 424036.495 4...  0.776678   \n",
       "2     12.291122  POLYGON ((424043.145 4505847.017, 424036.495 4...  0.776678   \n",
       "3     12.890390  POLYGON ((424421.505 4505865.604, 424414.855 4...  1.000000   \n",
       "4     12.890390  POLYGON ((424421.505 4505865.604, 424414.855 4...  1.000000   \n",
       "...         ...                                                ...       ...   \n",
       "1063  36.298872  POLYGON ((418169.004 4563560.006, 418123.854 4...  0.588706   \n",
       "1064  36.298872  POLYGON ((418169.004 4563560.006, 418123.854 4...  0.588706   \n",
       "1065  13.722026  POLYGON ((418384.296 4563850.293, 418366.563 4...  0.842694   \n",
       "1066  13.722026  POLYGON ((418384.296 4563850.293, 418366.563 4...  0.842694   \n",
       "1067  13.722026  POLYGON ((418384.296 4563850.293, 418366.563 4...  0.842694   \n",
       "\n",
       "     parameter  ...   acres_ch         s_ch        s_tot       s_diff  \\\n",
       "0       TOTEMP  ...   9.546246    57.277477    97.353816   -82.101744   \n",
       "1       TOTEMP  ...   9.546246    57.277477    98.828256   -87.229614   \n",
       "2       TOTEMP  ...   9.546246    57.277477   100.204991   -92.017686   \n",
       "3       TOTEMP  ...  12.890390    77.342342    77.342342  -168.865158   \n",
       "4       TOTEMP  ...  12.890390    77.342342    77.342342  -206.764208   \n",
       "...        ...  ...        ...          ...          ...          ...   \n",
       "1063     TOTHH  ...  21.369354  1944.611210  1996.818707  1869.884050   \n",
       "1064     TOTHH  ...  21.369354  1944.611210  2005.989409  1856.757573   \n",
       "1065     TOTHH  ...  11.563467  1052.275512  1052.950453  1048.659832   \n",
       "1066     TOTHH  ...  11.563467  1052.275512  1054.300336  1041.428471   \n",
       "1067     TOTHH  ...  11.563467  1052.275512  1056.601271  1029.102288   \n",
       "\n",
       "           s_grow  b_tot_den  f_tot_den  s_tot_den  s_grow_den  s_diff_den  \n",
       "0      -75.738072  14.082676  14.600422   7.920662   -6.162015   -6.679760  \n",
       "1      -74.263632  14.082676  15.137583   8.040621   -6.042055   -7.096961  \n",
       "2      -72.886896  14.082676  15.639149   8.152632   -5.930044   -7.486517  \n",
       "3     -139.145758  16.794534  19.100081   6.000000  -10.794534  -13.100081  \n",
       "4     -139.145758  16.794534  22.040182   6.000000  -10.794534  -16.040182  \n",
       "...           ...        ...        ...        ...         ...         ...  \n",
       "1063  1908.744848   2.426353   3.496931  55.010490   52.584137   51.513558  \n",
       "1064  1917.915550   2.426353   4.111198  55.263134   52.836781   51.151936  \n",
       "1065  1052.950453   0.000000   0.312681  76.734328   76.734328   76.421646  \n",
       "1066  1054.300336   0.000000   0.938044  76.832701   76.832701   75.894657  \n",
       "1067  1056.601271   0.000000   2.004003  77.000383   77.000383   74.996380  \n",
       "\n",
       "[1068 rows x 32 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "_gdf = gdf_stop_taz_calcs_se.copy()\n",
    "\n",
    "# Calculate expected growth for given polygon area\n",
    "# taz_ratio is ratio of parcels allowing hh/emp growth compared to rest of TAZ\n",
    "# value_B, value_F are base/future TAZ values\n",
    "\n",
    "# total base and future\n",
    "_gdf['b_tot'] = _gdf['taz_ratio'] * _gdf['value_B']\n",
    "_gdf['f_tot'] = _gdf['taz_ratio'] * _gdf['value_F']\n",
    "\n",
    "# changeable base and future\n",
    "_gdf['b_ch'] = _gdf['b_tot'] * _gdf['pctu_ch']\n",
    "_gdf['f_ch'] = _gdf['f_tot'] * _gdf['pctu_ch']\n",
    "\n",
    "# changeable acres\n",
    "_gdf['acres_ch'] = _gdf['acres'] * _gdf['pctu_ch']\n",
    "\n",
    "# changeable is the acres changeable multiplied by max per acre density\n",
    "_gdf['s_ch'] = _gdf['acres_ch'] * _gdf['perac_ch']\n",
    "\n",
    "# scenario total is forecast total - forecast changeable + scenario changeable\n",
    "_gdf['s_tot'] = _gdf['f_tot'] - _gdf['f_ch'] + _gdf['s_ch']\n",
    "\n",
    "# additional for scenario is the difference\n",
    "_gdf['s_diff'] = _gdf['s_tot'] - _gdf['f_tot']\n",
    "\n",
    "# growth from base\n",
    "_gdf['s_grow'] = _gdf['s_tot'] - _gdf['b_tot']\n",
    "\n",
    "# calculate per acre densities for polygon area\n",
    "_gdf['b_tot_den'] = _gdf['b_tot'] / _gdf['acres']\n",
    "_gdf['f_tot_den'] = _gdf['f_tot'] / _gdf['acres']\n",
    "_gdf['s_tot_den'] = _gdf['s_tot'] / _gdf['acres']\n",
    "\n",
    "_gdf['s_grow_den'] = (_gdf['s_tot'] - _gdf['b_tot']) / _gdf['acres']\n",
    "_gdf['s_diff_den'] = (_gdf['s_tot'] - _gdf['f_tot']) / _gdf['acres']\n",
    "\n",
    "_gdf['b_tot_den'] = (_gdf['taz_ratio'] * _gdf['value_B' ]) / _gdf['acres']\n",
    "_gdf['f_tot_den'] = (_gdf['taz_ratio'] * _gdf['value_F' ]) / _gdf['acres']\n",
    "\n",
    "gdf_se_calcs = _gdf\n",
    "\n",
    "gdf_se_calcs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>nearest_AreaType</th>\n",
       "      <th>TAZID</th>\n",
       "      <th>TAZACRES</th>\n",
       "      <th>sqm</th>\n",
       "      <th>acres</th>\n",
       "      <th>geometry</th>\n",
       "      <th>pctu_ch</th>\n",
       "      <th>parameter</th>\n",
       "      <th>...</th>\n",
       "      <th>acres_ch</th>\n",
       "      <th>s_ch</th>\n",
       "      <th>s_tot</th>\n",
       "      <th>s_diff</th>\n",
       "      <th>s_grow</th>\n",
       "      <th>b_tot_den</th>\n",
       "      <th>f_tot_den</th>\n",
       "      <th>s_tot_den</th>\n",
       "      <th>s_grow_den</th>\n",
       "      <th>s_diff_den</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>495</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>985.600584</td>\n",
       "      <td>1030.960347</td>\n",
       "      <td>-46.224958</td>\n",
       "      <td>-29.828973</td>\n",
       "      <td>20.619306</td>\n",
       "      <td>20.938006</td>\n",
       "      <td>20.039499</td>\n",
       "      <td>-0.579807</td>\n",
       "      <td>-0.898507</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>496</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>985.600584</td>\n",
       "      <td>1061.026497</td>\n",
       "      <td>-730.157570</td>\n",
       "      <td>0.237177</td>\n",
       "      <td>20.619306</td>\n",
       "      <td>34.816500</td>\n",
       "      <td>20.623916</td>\n",
       "      <td>0.004610</td>\n",
       "      <td>-14.192585</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>497</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>985.600584</td>\n",
       "      <td>1063.979268</td>\n",
       "      <td>-797.325995</td>\n",
       "      <td>3.189947</td>\n",
       "      <td>20.619306</td>\n",
       "      <td>36.179495</td>\n",
       "      <td>20.681311</td>\n",
       "      <td>0.062005</td>\n",
       "      <td>-15.498184</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>498</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>149.543844</td>\n",
       "      <td>162.418692</td>\n",
       "      <td>58.398719</td>\n",
       "      <td>59.982022</td>\n",
       "      <td>12.004211</td>\n",
       "      <td>12.189753</td>\n",
       "      <td>19.033304</td>\n",
       "      <td>7.029093</td>\n",
       "      <td>6.843551</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>499</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>149.543844</td>\n",
       "      <td>170.952623</td>\n",
       "      <td>-2.015683</td>\n",
       "      <td>68.515953</td>\n",
       "      <td>12.004211</td>\n",
       "      <td>20.269578</td>\n",
       "      <td>20.033366</td>\n",
       "      <td>8.029155</td>\n",
       "      <td>-0.236211</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>500</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>149.543844</td>\n",
       "      <td>171.790733</td>\n",
       "      <td>-7.948929</td>\n",
       "      <td>69.354063</td>\n",
       "      <td>12.004211</td>\n",
       "      <td>21.063090</td>\n",
       "      <td>20.131582</td>\n",
       "      <td>8.127371</td>\n",
       "      <td>-0.931508</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1029</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>2907.521723</td>\n",
       "      <td>2914.465893</td>\n",
       "      <td>2749.558542</td>\n",
       "      <td>2900.845845</td>\n",
       "      <td>0.264742</td>\n",
       "      <td>3.205420</td>\n",
       "      <td>56.650517</td>\n",
       "      <td>56.385774</td>\n",
       "      <td>53.445097</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1030</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>2907.521723</td>\n",
       "      <td>2915.952669</td>\n",
       "      <td>2715.737963</td>\n",
       "      <td>2902.332621</td>\n",
       "      <td>0.264742</td>\n",
       "      <td>3.891714</td>\n",
       "      <td>56.679416</td>\n",
       "      <td>56.414674</td>\n",
       "      <td>52.787703</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1031</th>\n",
       "      <td>15069</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>208196.248028</td>\n",
       "      <td>51.446413</td>\n",
       "      <td>POLYGON ((419387.011 4505039.208, 419373.709 4...</td>\n",
       "      <td>0.957890</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>49.280029</td>\n",
       "      <td>2907.521723</td>\n",
       "      <td>2916.526203</td>\n",
       "      <td>2702.691449</td>\n",
       "      <td>2902.906155</td>\n",
       "      <td>0.264742</td>\n",
       "      <td>4.156456</td>\n",
       "      <td>56.690564</td>\n",
       "      <td>56.425822</td>\n",
       "      <td>52.534108</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1032</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>441.154340</td>\n",
       "      <td>443.125363</td>\n",
       "      <td>427.200844</td>\n",
       "      <td>441.810124</td>\n",
       "      <td>0.154129</td>\n",
       "      <td>1.866141</td>\n",
       "      <td>51.928380</td>\n",
       "      <td>51.774251</td>\n",
       "      <td>50.062239</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1033</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>441.154340</td>\n",
       "      <td>443.547368</td>\n",
       "      <td>424.213342</td>\n",
       "      <td>442.232128</td>\n",
       "      <td>0.154129</td>\n",
       "      <td>2.265690</td>\n",
       "      <td>51.977833</td>\n",
       "      <td>51.823705</td>\n",
       "      <td>49.712143</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1034</th>\n",
       "      <td>23506</td>\n",
       "      <td>2</td>\n",
       "      <td>Urban Center</td>\n",
       "      <td>1384</td>\n",
       "      <td>96.889876</td>\n",
       "      <td>34533.423293</td>\n",
       "      <td>8.533395</td>\n",
       "      <td>POLYGON ((419457.973 4505447.709, 419448.379 4...</td>\n",
       "      <td>0.876227</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>...</td>\n",
       "      <td>7.477192</td>\n",
       "      <td>441.154340</td>\n",
       "      <td>443.710159</td>\n",
       "      <td>423.060893</td>\n",
       "      <td>442.394919</td>\n",
       "      <td>0.154129</td>\n",
       "      <td>2.419818</td>\n",
       "      <td>51.996910</td>\n",
       "      <td>51.842781</td>\n",
       "      <td>49.577092</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>12 rows × 32 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          N  project_id nearest_AreaType  TAZID   TAZACRES            sqm  \\\n",
       "495   15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "496   15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "497   15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "498   23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "499   23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "500   23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "1029  15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "1030  15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "1031  15069           2     Urban Center   1384  96.889876  208196.248028   \n",
       "1032  23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "1033  23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "1034  23506           2     Urban Center   1384  96.889876   34533.423293   \n",
       "\n",
       "          acres                                           geometry   pctu_ch  \\\n",
       "495   51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "496   51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "497   51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "498    8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "499    8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "500    8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "1029  51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "1030  51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "1031  51.446413  POLYGON ((419387.011 4505039.208, 419373.709 4...  0.957890   \n",
       "1032   8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "1033   8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "1034   8.533395  POLYGON ((419457.973 4505447.709, 419448.379 4...  0.876227   \n",
       "\n",
       "     parameter  ...   acres_ch         s_ch        s_tot       s_diff  \\\n",
       "495     TOTEMP  ...  49.280029   985.600584  1030.960347   -46.224958   \n",
       "496     TOTEMP  ...  49.280029   985.600584  1061.026497  -730.157570   \n",
       "497     TOTEMP  ...  49.280029   985.600584  1063.979268  -797.325995   \n",
       "498     TOTEMP  ...   7.477192   149.543844   162.418692    58.398719   \n",
       "499     TOTEMP  ...   7.477192   149.543844   170.952623    -2.015683   \n",
       "500     TOTEMP  ...   7.477192   149.543844   171.790733    -7.948929   \n",
       "1029     TOTHH  ...  49.280029  2907.521723  2914.465893  2749.558542   \n",
       "1030     TOTHH  ...  49.280029  2907.521723  2915.952669  2715.737963   \n",
       "1031     TOTHH  ...  49.280029  2907.521723  2916.526203  2702.691449   \n",
       "1032     TOTHH  ...   7.477192   441.154340   443.125363   427.200844   \n",
       "1033     TOTHH  ...   7.477192   441.154340   443.547368   424.213342   \n",
       "1034     TOTHH  ...   7.477192   441.154340   443.710159   423.060893   \n",
       "\n",
       "           s_grow  b_tot_den  f_tot_den  s_tot_den  s_grow_den  s_diff_den  \n",
       "495    -29.828973  20.619306  20.938006  20.039499   -0.579807   -0.898507  \n",
       "496      0.237177  20.619306  34.816500  20.623916    0.004610  -14.192585  \n",
       "497      3.189947  20.619306  36.179495  20.681311    0.062005  -15.498184  \n",
       "498     59.982022  12.004211  12.189753  19.033304    7.029093    6.843551  \n",
       "499     68.515953  12.004211  20.269578  20.033366    8.029155   -0.236211  \n",
       "500     69.354063  12.004211  21.063090  20.131582    8.127371   -0.931508  \n",
       "1029  2900.845845   0.264742   3.205420  56.650517   56.385774   53.445097  \n",
       "1030  2902.332621   0.264742   3.891714  56.679416   56.414674   52.787703  \n",
       "1031  2902.906155   0.264742   4.156456  56.690564   56.425822   52.534108  \n",
       "1032   441.810124   0.154129   1.866141  51.928380   51.774251   50.062239  \n",
       "1033   442.232128   0.154129   2.265690  51.977833   51.823705   49.712143  \n",
       "1034   442.394919   0.154129   2.419818  51.996910   51.842781   49.577092  \n",
       "\n",
       "[12 rows x 32 columns]"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_1384 = gdf_se_calcs[(gdf_se_calcs['TAZID']==1384)]\n",
    "df_1384 = pd.DataFrame(gdf_1384)\n",
    "df_1384"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "\n",
    "# Assuming gdf_se_calcs is your existing GeoDataFrame\n",
    "\n",
    "# Step 1: Filter rows where parameter is TOTHH\n",
    "gdf_tothh = gdf_se_calcs[gdf_se_calcs['parameter'] == 'TOTHH'].copy()\n",
    "\n",
    "# Step 2: Identify numeric columns that don't contain 'pct' or 'acres'\n",
    "columns_to_multiply = [col for col in gdf_tothh.columns if col not in ['parameter','project_id','nearest_AreaType','TAZID','TAZACRES','taz_ratio','sqm','geometry','N'] and 'year' not in col and 'pct' not in col and 'acres' not in col and gdf_tothh[col].dtype != 'object']\n",
    "\n",
    "# Step 3: Multiply selected columns by 1.8\n",
    "gdf_tothh[columns_to_multiply] *= hh_factor\n",
    "\n",
    "# Step 4: Filter rows where parameter is TOTEMP\n",
    "gdf_totemp = gdf_se_calcs[gdf_se_calcs['parameter'] == 'TOTEMP'].copy()\n",
    "\n",
    "# Step 5: Ensure the indices match for correct addition\n",
    "gdf_tothh.index = gdf_totemp.index\n",
    "\n",
    "# Step 6: Create a new DataFrame for the sum of TOTHH * 1.8 and TOTEMP\n",
    "gdf_tothhemp = gdf_totemp.copy()\n",
    "gdf_tothhemp[columns_to_multiply] += gdf_tothh[columns_to_multiply]\n",
    "gdf_tothhemp['parameter'] = 'TOTHHEMP'\n",
    "\n",
    "# Step 7: Append the new rows to the original GeoDataFrame\n",
    "gdf_se_calcs2 = pd.concat([gdf_se_calcs, gdf_tothhemp], ignore_index=True)\n",
    "\n",
    "# Now gdf_se_calcs should have the new TOTHHEMP rows added\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>se_year_B</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1597</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1598</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1599</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1600</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1601</th>\n",
       "      <td>2023</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1602 rows × 1 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      se_year_B\n",
       "0          2023\n",
       "1          2023\n",
       "2          2023\n",
       "3          2023\n",
       "4          2023\n",
       "...         ...\n",
       "1597       2023\n",
       "1598       2023\n",
       "1599       2023\n",
       "1600       2023\n",
       "1601       2023\n",
       "\n",
       "[1602 rows x 1 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_se_calcs2[['se_year_B']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "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>TAZID</th>\n",
       "      <th>parameter</th>\n",
       "      <th>N</th>\n",
       "      <th>project_id</th>\n",
       "      <th>TAZACRES</th>\n",
       "      <th>sqm</th>\n",
       "      <th>acres</th>\n",
       "      <th>pctu_ch</th>\n",
       "      <th>taz_ratio</th>\n",
       "      <th>perac_ch</th>\n",
       "      <th>...</th>\n",
       "      <th>acres_ch</th>\n",
       "      <th>s_ch</th>\n",
       "      <th>s_tot</th>\n",
       "      <th>s_diff</th>\n",
       "      <th>s_grow</th>\n",
       "      <th>b_tot_den</th>\n",
       "      <th>f_tot_den</th>\n",
       "      <th>s_tot_den</th>\n",
       "      <th>s_grow_den</th>\n",
       "      <th>s_diff_den</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>165</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>45573</td>\n",
       "      <td>3</td>\n",
       "      <td>14313.894679</td>\n",
       "      <td>380792.949177</td>\n",
       "      <td>94.095987</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.025288</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>19.401624</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>11.086040</td>\n",
       "      <td>0.265120</td>\n",
       "      <td>0.618569</td>\n",
       "      <td>0.618569</td>\n",
       "      <td>0.353449</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>165</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>45573</td>\n",
       "      <td>3</td>\n",
       "      <td>14313.894679</td>\n",
       "      <td>380792.949177</td>\n",
       "      <td>94.095987</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>54</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>172</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>45555</td>\n",
       "      <td>3</td>\n",
       "      <td>1165.075927</td>\n",
       "      <td>112874.323214</td>\n",
       "      <td>27.891853</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>172</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>45555</td>\n",
       "      <td>3</td>\n",
       "      <td>1165.075927</td>\n",
       "      <td>112874.323214</td>\n",
       "      <td>27.891853</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.084408</td>\n",
       "      <td>54</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>15.423307</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>8.712832</td>\n",
       "      <td>0.721767</td>\n",
       "      <td>1.658905</td>\n",
       "      <td>1.658905</td>\n",
       "      <td>0.937137</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>173</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>45552</td>\n",
       "      <td>3</td>\n",
       "      <td>1515.829447</td>\n",
       "      <td>32155.649119</td>\n",
       "      <td>7.945834</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</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>305</th>\n",
       "      <td>1530</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>76143</td>\n",
       "      <td>6</td>\n",
       "      <td>944.282336</td>\n",
       "      <td>27988.527121</td>\n",
       "      <td>6.916116</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000709</td>\n",
       "      <td>54</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.532118</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.049821</td>\n",
       "      <td>0.209206</td>\n",
       "      <td>0.230816</td>\n",
       "      <td>0.230816</td>\n",
       "      <td>0.021611</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>306</th>\n",
       "      <td>1531</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>76143</td>\n",
       "      <td>6</td>\n",
       "      <td>1477.929288</td>\n",
       "      <td>323888.238224</td>\n",
       "      <td>80.034527</td>\n",
       "      <td>1.877119</td>\n",
       "      <td>1.194773</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>50.078116</td>\n",
       "      <td>300.468696</td>\n",
       "      <td>679.914376</td>\n",
       "      <td>-333.850448</td>\n",
       "      <td>-352.568557</td>\n",
       "      <td>38.701407</td>\n",
       "      <td>37.999781</td>\n",
       "      <td>25.485790</td>\n",
       "      <td>-13.215617</td>\n",
       "      <td>-12.513991</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>307</th>\n",
       "      <td>1531</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>76143</td>\n",
       "      <td>6</td>\n",
       "      <td>1477.929288</td>\n",
       "      <td>323888.238224</td>\n",
       "      <td>80.034527</td>\n",
       "      <td>1.877119</td>\n",
       "      <td>0.077309</td>\n",
       "      <td>54</td>\n",
       "      <td>...</td>\n",
       "      <td>50.078116</td>\n",
       "      <td>901.406088</td>\n",
       "      <td>929.263505</td>\n",
       "      <td>854.836860</td>\n",
       "      <td>866.269933</td>\n",
       "      <td>2.361240</td>\n",
       "      <td>2.789795</td>\n",
       "      <td>34.832348</td>\n",
       "      <td>32.471108</td>\n",
       "      <td>32.042553</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>308</th>\n",
       "      <td>1679</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>76143</td>\n",
       "      <td>6</td>\n",
       "      <td>232.749424</td>\n",
       "      <td>28916.183832</td>\n",
       "      <td>7.145345</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>7.145345</td>\n",
       "      <td>42.872068</td>\n",
       "      <td>42.872068</td>\n",
       "      <td>42.872068</td>\n",
       "      <td>42.872068</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>18.000000</td>\n",
       "      <td>18.000000</td>\n",
       "      <td>18.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>309</th>\n",
       "      <td>1679</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>76143</td>\n",
       "      <td>6</td>\n",
       "      <td>232.749424</td>\n",
       "      <td>28916.183832</td>\n",
       "      <td>7.145345</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>0.001409</td>\n",
       "      <td>54</td>\n",
       "      <td>...</td>\n",
       "      <td>7.145345</td>\n",
       "      <td>128.616204</td>\n",
       "      <td>128.616204</td>\n",
       "      <td>128.304706</td>\n",
       "      <td>128.328198</td>\n",
       "      <td>0.120920</td>\n",
       "      <td>0.130783</td>\n",
       "      <td>54.000000</td>\n",
       "      <td>53.879080</td>\n",
       "      <td>53.869217</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>310 rows × 29 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     TAZID parameter      N  project_id      TAZACRES            sqm  \\\n",
       "0      165    TOTEMP  45573           3  14313.894679  380792.949177   \n",
       "1      165     TOTHH  45573           3  14313.894679  380792.949177   \n",
       "2      172    TOTEMP  45555           3   1165.075927  112874.323214   \n",
       "3      172     TOTHH  45555           3   1165.075927  112874.323214   \n",
       "4      173    TOTEMP  45552           3   1515.829447   32155.649119   \n",
       "..     ...       ...    ...         ...           ...            ...   \n",
       "305   1530     TOTHH  76143           6    944.282336   27988.527121   \n",
       "306   1531    TOTEMP  76143           6   1477.929288  323888.238224   \n",
       "307   1531     TOTHH  76143           6   1477.929288  323888.238224   \n",
       "308   1679    TOTEMP  76143           6    232.749424   28916.183832   \n",
       "309   1679     TOTHH  76143           6    232.749424   28916.183832   \n",
       "\n",
       "         acres   pctu_ch  taz_ratio  perac_ch  ...   acres_ch        s_ch  \\\n",
       "0    94.095987  0.000000   0.025288        18  ...   0.000000    0.000000   \n",
       "1    94.095987  0.000000   0.000000        54  ...   0.000000    0.000000   \n",
       "2    27.891853  0.000000   0.000000        18  ...   0.000000    0.000000   \n",
       "3    27.891853  0.000000   0.084408        54  ...   0.000000    0.000000   \n",
       "4     7.945834  0.000000   0.000000        18  ...   0.000000    0.000000   \n",
       "..         ...       ...        ...       ...  ...        ...         ...   \n",
       "305   6.916116  0.000000   0.000709        54  ...   0.000000    0.000000   \n",
       "306  80.034527  1.877119   1.194773        18  ...  50.078116  300.468696   \n",
       "307  80.034527  1.877119   0.077309        54  ...  50.078116  901.406088   \n",
       "308   7.145345  3.000000   0.000000        18  ...   7.145345   42.872068   \n",
       "309   7.145345  3.000000   0.001409        54  ...   7.145345  128.616204   \n",
       "\n",
       "          s_tot      s_diff      s_grow  b_tot_den  f_tot_den  s_tot_den  \\\n",
       "0     19.401624    0.000000   11.086040   0.265120   0.618569   0.618569   \n",
       "1      0.000000    0.000000    0.000000   0.000000   0.000000   0.000000   \n",
       "2      0.000000    0.000000    0.000000   0.000000   0.000000   0.000000   \n",
       "3     15.423307    0.000000    8.712832   0.721767   1.658905   1.658905   \n",
       "4      0.000000    0.000000    0.000000   0.000000   0.000000   0.000000   \n",
       "..          ...         ...         ...        ...        ...        ...   \n",
       "305    0.532118    0.000000    0.049821   0.209206   0.230816   0.230816   \n",
       "306  679.914376 -333.850448 -352.568557  38.701407  37.999781  25.485790   \n",
       "307  929.263505  854.836860  866.269933   2.361240   2.789795  34.832348   \n",
       "308   42.872068   42.872068   42.872068   0.000000   0.000000  18.000000   \n",
       "309  128.616204  128.304706  128.328198   0.120920   0.130783  54.000000   \n",
       "\n",
       "     s_grow_den  s_diff_den  \n",
       "0      0.353449    0.000000  \n",
       "1      0.000000    0.000000  \n",
       "2      0.000000    0.000000  \n",
       "3      0.937137    0.000000  \n",
       "4      0.000000    0.000000  \n",
       "..          ...         ...  \n",
       "305    0.021611    0.000000  \n",
       "306  -13.215617  -12.513991  \n",
       "307   32.471108   32.042553  \n",
       "308   18.000000   18.000000  \n",
       "309   53.879080   53.869217  \n",
       "\n",
       "[310 rows x 29 columns]"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "_gdf.groupby(['TAZID','parameter'], as_index=False).sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "# filter for map year\n",
    "gdf_se_calcs_to_map = gdf_se_calcs2[gdf_se_calcs2['se_year_F']==_config.map_year]\n",
    "\n",
    "# Project to WGS 84 (EPSG:4326)\n",
    "gdf_se_calcs_to_map = gdf_se_calcs_to_map.to_crs(epsg=4326)\n",
    "\n",
    "# Export to GeoJSON\n",
    "gdf_se_calcs_to_map.to_file('_site/data/se_calcs.geojson', driver='GeoJSON')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "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>taz_ratio</th>\n",
       "      <th>parameter</th>\n",
       "      <th>N</th>\n",
       "      <th>se_year_F</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>534</th>\n",
       "      <td>0.419705</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>15086</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>535</th>\n",
       "      <td>0.419705</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>15086</td>\n",
       "      <td>2042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>536</th>\n",
       "      <td>0.419705</td>\n",
       "      <td>TOTHH</td>\n",
       "      <td>15086</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1068</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTHHEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2032</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1069</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTHHEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1070</th>\n",
       "      <td>0.238638</td>\n",
       "      <td>TOTHHEMP</td>\n",
       "      <td>15086</td>\n",
       "      <td>2050</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      taz_ratio parameter      N  se_year_F\n",
       "0      0.238638    TOTEMP  15086       2032\n",
       "1      0.238638    TOTEMP  15086       2042\n",
       "2      0.238638    TOTEMP  15086       2050\n",
       "534    0.419705     TOTHH  15086       2032\n",
       "535    0.419705     TOTHH  15086       2042\n",
       "536    0.419705     TOTHH  15086       2050\n",
       "1068   0.238638  TOTHHEMP  15086       2032\n",
       "1069   0.238638  TOTHHEMP  15086       2042\n",
       "1070   0.238638  TOTHHEMP  15086       2050"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_se_calcs2[(gdf_se_calcs2['N']==15086) & (gdf_se_calcs2['TAZID']==1472)][['taz_ratio','parameter','N','se_year_F']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create HTML table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<h3>Project Name: West Weber Rail</h3>\n",
      "<table border=\"1\" class=\"dataframe\">\n",
      "  <thead>\n",
      "    <tr style=\"text-align: right;\">\n",
      "      <th>parameter</th>\n",
      "      <th>b_tot</th>\n",
      "      <th>f_tot</th>\n",
      "      <th>s_tot</th>\n",
      "    </tr>\n",
      "  </thead>\n",
      "  <tbody>\n",
      "    <tr>\n",
      "      <td>Employees</td>\n",
      "      <td>4,597</td>\n",
      "      <td>5,082</td>\n",
      "      <td>7,143</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>Households</td>\n",
      "      <td>144</td>\n",
      "      <td>662</td>\n",
      "      <td>18,230</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>HH+Emp Intensity</td>\n",
      "      <td>4,857</td>\n",
      "      <td>6,274</td>\n",
      "      <td>39,957</td>\n",
      "    </tr>\n",
      "  </tbody>\n",
      "</table>\n",
      "<h3>Project Name: 3300/3500 South</h3>\n",
      "<table border=\"1\" class=\"dataframe\">\n",
      "  <thead>\n",
      "    <tr style=\"text-align: right;\">\n",
      "      <th>parameter</th>\n",
      "      <th>b_tot</th>\n",
      "      <th>f_tot</th>\n",
      "      <th>s_tot</th>\n",
      "    </tr>\n",
      "  </thead>\n",
      "  <tbody>\n",
      "    <tr>\n",
      "      <td>Employees</td>\n",
      "      <td>15,590</td>\n",
      "      <td>16,786</td>\n",
      "      <td>11,987</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>Households</td>\n",
      "      <td>4,172</td>\n",
      "      <td>6,301</td>\n",
      "      <td>24,646</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>HH+Emp Intensity</td>\n",
      "      <td>23,099</td>\n",
      "      <td>28,127</td>\n",
      "      <td>56,350</td>\n",
      "    </tr>\n",
      "  </tbody>\n",
      "</table>\n",
      "<h3>Project Name: Roy to Clearfield via 3500 West</h3>\n",
      "<table border=\"1\" class=\"dataframe\">\n",
      "  <thead>\n",
      "    <tr style=\"text-align: right;\">\n",
      "      <th>parameter</th>\n",
      "      <th>b_tot</th>\n",
      "      <th>f_tot</th>\n",
      "      <th>s_tot</th>\n",
      "    </tr>\n",
      "  </thead>\n",
      "  <tbody>\n",
      "    <tr>\n",
      "      <td>Employees</td>\n",
      "      <td>5,257</td>\n",
      "      <td>7,459</td>\n",
      "      <td>8,077</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>Households</td>\n",
      "      <td>1,804</td>\n",
      "      <td>2,901</td>\n",
      "      <td>21,290</td>\n",
      "    </tr>\n",
      "    <tr>\n",
      "      <td>HH+Emp Intensity</td>\n",
      "      <td>8,505</td>\n",
      "      <td>12,681</td>\n",
      "      <td>46,400</td>\n",
      "    </tr>\n",
      "  </tbody>\n",
      "</table>\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\bhereth\\AppData\\Local\\Temp\\ipykernel_24104\\579806420.py:29: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_project[['b_tot', 'f_tot', 's_tot']] = df_project[['b_tot', 'f_tot', 's_tot']].applymap(format_numbers)\n",
      "C:\\Users\\bhereth\\AppData\\Local\\Temp\\ipykernel_24104\\579806420.py:29: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_project[['b_tot', 'f_tot', 's_tot']] = df_project[['b_tot', 'f_tot', 's_tot']].applymap(format_numbers)\n",
      "C:\\Users\\bhereth\\AppData\\Local\\Temp\\ipykernel_24104\\579806420.py:29: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame.\n",
      "Try using .loc[row_indexer,col_indexer] = value instead\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  df_project[['b_tot', 'f_tot', 's_tot']] = df_project[['b_tot', 'f_tot', 's_tot']].applymap(format_numbers)\n"
     ]
    }
   ],
   "source": [
    "# summary table\n",
    "_df = gdf_se_calcs2[gdf_se_calcs2['se_year_F']==2050]\n",
    "\n",
    "_df = _df[['project_id','center_type','b_tot','f_tot','s_tot','parameter',]]\n",
    "\n",
    "df_grouped = _df.groupby(['project_id','parameter'], as_index=False).sum()\n",
    "\n",
    "# Replacing values in the 'parameter' column\n",
    "df_grouped['parameter'] = df_grouped['parameter'].replace({\n",
    "    'TOTHH': 'Households',\n",
    "    'TOTEMP': 'Employees',\n",
    "    'TOTHHEMP': 'HH+Emp Intensity'\n",
    "})\n",
    "# Merging df_grouped with df_projects to get project_name\n",
    "df_grouped = df_grouped.merge(df_projects[['project_id', 'project_name']], on='project_id')\n",
    "\n",
    "# Getting unique project_names\n",
    "project_names = df_grouped['project_name'].unique()\n",
    "\n",
    "# Function to format numbers\n",
    "def format_numbers(value):\n",
    "    return \"{:,.0f}\".format(value)\n",
    "\n",
    "# Iterate over each project_name and create a separate HTML table\n",
    "for project_name in project_names:\n",
    "    df_project = df_grouped[df_grouped['project_name'] == project_name]\n",
    "    \n",
    "    # Apply number formatting\n",
    "    df_project[['b_tot', 'f_tot', 's_tot']] = df_project[['b_tot', 'f_tot', 's_tot']].applymap(format_numbers)\n",
    "    \n",
    "    # Drop project_id and project_name column for the table\n",
    "    df_project = df_project.drop(columns=['project_id', 'project_name'])\n",
    "    \n",
    "    # Convert to HTML\n",
    "    html_table = df_project.to_html(index=False, escape=False)\n",
    "    \n",
    "    # Add h3 header for project_name\n",
    "    html_output = f\"<h3>Project Name: {project_name}</h3>\\n{html_table}\"\n",
    "    \n",
    "    print(html_output)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_se_calcs2.to_file('intermediate/se_calcs2.geojson', format='GeoJSON')\n",
    "df_se_data.to_csv('intermediate/se_data.csv', index=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "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></th>\n",
       "      <th>s_tot</th>\n",
       "      <th>s_ch</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>parameter</th>\n",
       "      <th>se_year_F</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"3\" valign=\"top\">TOTEMP</th>\n",
       "      <th>2032</th>\n",
       "      <td>26691.841753</td>\n",
       "      <td>20005.133482</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2042</th>\n",
       "      <td>27065.919043</td>\n",
       "      <td>20005.133482</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2050</th>\n",
       "      <td>27207.466663</td>\n",
       "      <td>20005.133482</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"3\" valign=\"top\">TOTHH</th>\n",
       "      <th>2032</th>\n",
       "      <td>63239.683618</td>\n",
       "      <td>59766.854637</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2042</th>\n",
       "      <td>63715.339044</td>\n",
       "      <td>59766.854637</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2050</th>\n",
       "      <td>64166.329236</td>\n",
       "      <td>59766.854637</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"3\" valign=\"top\">TOTHHEMP</th>\n",
       "      <th>2032</th>\n",
       "      <td>140523.272266</td>\n",
       "      <td>127585.471829</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2042</th>\n",
       "      <td>141753.529322</td>\n",
       "      <td>127585.471829</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2050</th>\n",
       "      <td>142706.859288</td>\n",
       "      <td>127585.471829</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                             s_tot           s_ch\n",
       "parameter se_year_F                              \n",
       "TOTEMP    2032        26691.841753   20005.133482\n",
       "          2042        27065.919043   20005.133482\n",
       "          2050        27207.466663   20005.133482\n",
       "TOTHH     2032        63239.683618   59766.854637\n",
       "          2042        63715.339044   59766.854637\n",
       "          2050        64166.329236   59766.854637\n",
       "TOTHHEMP  2032       140523.272266  127585.471829\n",
       "          2042       141753.529322  127585.471829\n",
       "          2050       142706.859288  127585.471829"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_se_calcs2[['parameter','se_year_F','s_tot','s_ch']].groupby(['parameter','se_year_F']).sum()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "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.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}