{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - Raster DataArray to GeoDataFrame (vectorize)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import geopandas\n", "\n", "from geocube.api.core import make_geocube\n", "from geocube.vector import vectorize\n", "\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Generate Raster DataArray" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAHFCAYAAAD7ZFORAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABtUklEQVR4nO3deViUVfsH8O+AMoAsIiiLIqJpLrgglqGVWAKSoVavS5q7ppELuebPUjQRdynNNXMptywtc6fMLbUU8X3d0swFQpBUBNxYZs7vD/NxBhiclZlhvp/req5mnuXM/QDJzTnnfo5MCCFAREREVAHZmTsAIiIiIlNhokNEREQVFhMdIiIiqrCY6BAREVGFxUSHiIiIKiwmOkRERFRhMdEhIiKiCouJDhEREVVYTHSIiIiowmKiQ2RlFi9ejNWrV5fYf/XqVchkslKPPY0h1xrq559/RqtWrVClShXIZDJ8//335fbZGzduRIsWLeDo6Ag/Pz/Exsbi7t275fb5RGR6lcwdABHpZvHixfDy8kL//v3V9vv6+uLo0aOoV6+eeQLTgxAC3bt3R4MGDbBt2zZUqVIFzz77bLl89rp16/DOO+9g8ODBWLBgAS5evIgJEybg3Llz2Lt3b7nEQESmx0SHqIKQy+V44YUXzBqDQqFAUVER5HK5Vudfv34dt2/fxhtvvIFXX33VxNE9oVAoMG7cOERERGDFihUAgPbt28PV1RW9e/fGrl27EBUVVW7xEJHpcOiKSE///PMP3n33Xfj7+0Mul6N69epo27YtfvrpJ+mcsLAwBAUF4dChQ3jhhRfg5OSEmjVr4uOPP4ZCoVBrb+rUqWjdujWqVasGNzc3tGzZEitXroTqurt16tTB2bNnceDAAchkMshkMtSpUwdA6cNPly5dwoABA1C/fn04OzujZs2aiI6OxunTpw2+/8efN3v2bEyfPh2BgYGQy+X45ZdfAAAnTpxA586dUa1aNTg6OiI4OBjffPONdH1cXBxq1aoFAJgwYYLavZjasWPHkJGRgQEDBqjt79atG1xcXLB169ZyiYOITI89OkR66tOnD06ePIn4+Hg0aNAAd+7cwcmTJ3Hr1i218zIzM9GzZ098+OGHmDZtGnbs2IHp06cjOzsbixYtks67evUqhg4ditq1awN49Mt4xIgRSE9Px+TJkwEAW7duxX/+8x+4u7tj8eLFAFBm78n169fh6emJmTNnonr16rh9+zbWrFmD1q1bIyUlxSjDRJ999hkaNGiAuXPnws3NDfXr18cvv/yCjh07onXr1li6dCnc3d2xceNG9OjRA/fv30f//v0xePBgNG/eHG+++SZGjBiBXr16PbUnqKioSKuY7O3tIZPJNB4/c+YMAKBZs2Zq+ytXroyGDRtKx4moAhBEpBcXFxcRGxtb5jnt2rUTAMQPP/ygtn/IkCHCzs5OXLt2rdTrFAqFKCwsFNOmTROenp5CqVRKx5o0aSLatWtX4porV64IAGLVqlUa4ykqKhIFBQWifv364oMPPtDpWk2fV69ePVFQUKB2rGHDhiI4OFgUFhaq7X/99deFr6+vUCgUam3MmTNHq88EoNX2tPuIj48XAERGRkaJYxEREaJBgwZaxUNElo9DV0R6ev7557F69WpMnz4dx44dQ2FhYannubq6onPnzmr7evXqBaVSiYMHD0r79u3bhw4dOsDd3R329vaoXLkyJk+ejFu3biErK0uvGIuKijBjxgw0btwYDg4OqFSpEhwcHPDnn3/i/PnzerVZXOfOnVG5cmXp/aVLl/DHH3+gd+/eUgyPt9deew0ZGRm4cOGCXp91/Phxrbbo6Git2tPU61NWbxCRqRw8eBDR0dHw8/PTuwLxm2++QYsWLeDs7IyAgADMmTPH+IFaGQ5dEelp06ZNmD59Or744gt8/PHHcHFxwRtvvIHZs2fDx8dHOs/b27vEtY+PPx7m+v333xEREYGwsDCsWLECtWrVgoODA77//nvEx8fjwYMHesU4evRofP7555gwYQLatWsHDw8P2NnZYfDgwXq3WZyvr6/a+xs3bgAAxo4di7Fjx5Z6zc2bN/X6rBYtWmh1nr29fZnHPT09ATz6+hf//ty+fRvVqlXTKz4iQ9y7dw/NmzfHgAED8NZbb+l8/a5du9C7d28sXLgQEREROH/+PAYPHgwnJycMHz7cBBFbByY6RHry8vJCYmIiEhMTkZqaim3btuHDDz9EVlYWdu/eLZ33+Be/qszMTABPfuFu3LgRlStXxvbt2+Ho6CidZ+gzZb7++mv07dsXM2bMUNt/8+ZNVK1a1aC2Hyve++Hl5QUAmDhxIt58881Sr9F3bpBqz1FZVq1aVaL8XlXTpk0BAKdPn0bjxo2l/UVFRfjjjz/w9ttv6xUfkSGioqLKrPYrKCjARx99hHXr1uHOnTsICgrCrFmzEBYWBgD46quv0LVrVwwbNgwAULduXUyYMAGzZs3C+++/b7M9lUx0iIygdu3aGD58OH7++Wf8+uuvasfy8vKwbds2teGr9evXw87ODi+//DKAR8lCpUqV1HoiHjx4gK+++qrEZ8nlcq17Y2QyWYkJvjt27EB6ejqeeeYZre9PF88++yzq16+P//73vyUSLEMdP35cq/MCAwPLPN66dWv4+vpi9erV6NGjh7T/22+/xd27dzUmaETmNGDAAFy9ehUbN26En58ftm7dio4dO+L06dOoX78+8vPz4ezsrHaNk5MT/v77b1y7dq3cqhotDRMdIj3k5OSgffv26NWrFxo2bAhXV1ccP34cu3fvLvFL0tPTE++99x5SU1PRoEED7Ny5EytWrMB7770nVVh16tQJ8+fPR69evfDuu+/i1q1bmDt3bqlVSE2bNsXGjRuxadMm1K1bF46OjlIPRXGvv/46Vq9ejYYNG6JZs2ZITk7GnDlzpLJuU1m2bBmioqIQGRmJ/v37o2bNmrh9+zbOnz+PkydPYvPmzXq126pVK6PEZ29vj9mzZ6NPnz4YOnQo3n77bfz5558YP348wsPD0bFjR6N8DpGx/PXXX9iwYQP+/vtv+Pn5AXg0PLx7926sWrUKM2bMQGRkJD744AP0798f7du3x6VLl5CYmAgAyMjIYKJDRNpzdHRE69at8dVXX+Hq1asoLCxE7dq1MWHCBIwfP17tXB8fH3z++ecYO3YsTp8+jWrVquH//u//MHXqVOmcV155BV9++SVmzZqF6Oho1KxZE0OGDEGNGjUwaNAgtfamTp2KjIwMDBkyBHl5eQgICMDVq1dLjfPTTz9F5cqVkZCQgLt376Jly5bYsmULPvroI6N/TVS1b98ev//+O+Lj4xEbG4vs7Gx4enqicePG6N69u0k/W1vvvPMO7O3tMXPmTKxevRrVqlVD3759ER8fb+7QiEo4efIkhBBo0KCB2v78/HxpCHzIkCH466+/8Prrr6OwsBBubm4YNWoU4uLinjpvrSKTCaHyNDIiMqqwsDDcvHmTz2UhIp3IZDJs3boVXbt2BfCo+KF37944e/ZsiaTFxcVFrQBCoVAgMzMT1atXx88//4zXXnsNN27cQI0aNcrzFiwGe3SIiIgsXHBwMBQKBbKysvDSSy+Vea69vT1q1qwJANiwYQNCQ0NtNskBmOgQUSme9gRiOzs72NnxMVxExnT37l1cunRJen/lyhWcOnUK1apVQ4MGDdC7d2/07dsX8+bNQ3BwMG7evIl9+/ahadOmeO2113Dz5k18++23CAsLw8OHD7Fq1Sps3rwZBw4cMONdmR+HrohIzdWrV59atTRlyhTExcWVT0BENmL//v1o3759if39+vXD6tWrUVhYiOnTp2Pt2rVIT0+Hp6cnQkNDMXXqVDRt2hQ3b96U1rITQiA0NBTx8fFo3bq1Ge7GcjDRISI1BQUF+N///lfmOX5+flLlBxGRJWOiQ0RERBUWB9mJiIiowuJk5GKUSiWuX78OV1dXm31cNhERaUcIgby8PPj5+Zlsgv7Dhw9RUFBglLYcHBzUlpmxBUx0irl+/Tr8/f3NHQYREVmRtLQ0kzxx/OHDhwgMcEFmlsIo7fn4+ODKlSs2leww0SnG1dUVwKMfWjc3NzNHQ6re8Hm33D5ra+Zyg67v4t7XSJHo5oectWb5XCJblZubC39/f+l3h7EVFBQgM0uBK8kBcHM1rMcoN0+JwJBrKCgoYKJjyx4PV7m5uTHRsTCVZA7l9lmGfu8rybRbZdvY+DNLZB6mnurg5mpncKJjq5joEBERWTiFUEJhYI20QiiNE4yVsZr0sKioCB999BECAwPh5OSEunXrYtq0aVAqn3zjhBCIi4uDn58fnJycEBYWhrNnz5oxaiIiIsMpIYyy2SKrSXRmzZqFpUuXYtGiRTh//jxmz56NOXPmYOHChdI5s2fPxvz587Fo0SIcP34cPj4+CA8PR15enhkjJyIiInOxmqGro0ePokuXLujUqRMAoE6dOtiwYQNOnDgB4FFvTmJiIiZNmoQ333wTALBmzRp4e3tj/fr1GDp0qNlip7JFVjHtxN09954+QTfcrluZ762FatxJys1mjOSRsr6OqvEVP88SYrcmUQ0mlLp/18VZ5RwJmYoSShg68GR4C9bJanp0XnzxRfz888+4ePEiAOC///0vDh8+jNdeew3Ao8XPMjMzERERIV0jl8vRrl07HDlyRGO7+fn5yM3NVduIiIgsiUIIo2y2yGp6dCZMmICcnBw0bNgQ9vb2UCgUiI+Px9tvvw0AyMzMBAB4e3urXeft7Y1r165pbDchIQFTp041XeBERERkNlbTo7Np0yZ8/fXXWL9+PU6ePIk1a9Zg7ty5WLNmjdp5xUv8hBBllv1NnDgROTk50paWlmaS+ImIiPTFycj6s5oenXHjxuHDDz9Ez549AQBNmzbFtWvXkJCQgH79+sHHxwfAo54dX19f6bqsrKwSvTyq5HI55HK5aYMnIiIygBICCgMTFVtNdKymR+f+/fsl1hGxt7eXyssDAwPh4+ODpKQk6XhBQQEOHDiANm3alGusREREZBmspkcnOjoa8fHxqF27Npo0aYKUlBTMnz8fAwcOBPBoyCo2NhYzZsxA/fr1Ub9+fcyYMQPOzs7o1auXmaMnIiLSnzGGnmy1R8dqEp2FCxfi448/RkxMDLKysuDn54ehQ4di8uTJ0jnjx4/HgwcPEBMTg+zsbLRu3Rp79+412RokFZGmUu+ySrRVrzHGeYbStm1rLSEvi6WVZWsbj6XFrS/VMu/yLO3W9rM0laHr2x6VH2NUTbHqysK5uroiMTERiYmJGs+RyWSIi4tDXFxcucVFRERkasp/N0PbsEVWM0eHiIiISFdW06NDRERkqxRGqLoy9HprxUSHiIjIwikEjLB6uXFisTYcuiIiIqIKiz06ZDDVKid9KrD0+Rx921ClT6VPRazUotLpU0FlaLWStlVRpqZNHMXv1dgVZ+aqYLNUnIysPyY6REREFk4JGRTQvJyRtm3YIg5dERERUYXFHh0iIiILpxSPNkPbsEVMdIiIiCycwghDV4Zeb604dEVEREQVlkwIG138QoPc3Fy4u7sjJycHbm5uRmnz5S5zdL7m4A/jjPLZpmKMqilNbWmq4jI2Y6y3ZY4qrIqyLlRZNH1dLf3ey6pWUq0cspTqKlthyqotU/zOKK39I2d94eJqWN/E3Twl2jTJMFmslopDV0RERBZOKWRQCgOrrgy83lox0SEiIrJwnKOjP87RISIiogqLPTpEREQWTgE7KAzsm1AYKRZrw0SHiIjIwgkjzNERNjpHh0NXREREVGGxR0eDt4Ino5K9XK+yRH3KyfVtw+mn0wZ/VnnRtlTc0NJ1bcvVy/ocY5SeG8rSS6n1EW7fQ7sTZaX/DaZPWXZZ/w+zzNs2mPL7XKTIN1nbqjgZWX9MdIiIiCycQthBIQyco2OjT83j0BURERFVWOzRISIisnBKyKA0sG9CCdvs0mGiQ0REZOE4R0d/HLoiIiKiCouLehZjKYt6asuaqq5MydCFQMtzgU9TVlOVFYM5qrhKxKOhmipJsUn9Og3VWZXq1TFGWEQGUa3kK69FPbf+tz6quNob1Na9PAXeaP4nF/UkIiIiy/Jojo6Bi3ra6NAVEx0iIiILpzTCEhC2OhmZc3SIiIiowmKPDhERkYUzzgMDbbNHh4kOERGRhVPCjs/R0ROHroiIiKjCYo+OCZiynNwWmHJBTVO2rVq+rW2puSljKE8a71dDOTmgXlKu9WKfRFpSLQEvvqhnWce0ac8cFEIGhTDwgYEGXm+tmOgQERFZOIURqq4UHLqyfOnp6XjnnXfg6ekJZ2dntGjRAsnJydJxIQTi4uLg5+cHJycnhIWF4ezZs2aMmIiIiMzJahKd7OxstG3bFpUrV8auXbtw7tw5zJs3D1WrVpXOmT17NubPn49Fixbh+PHj8PHxQXh4OPLy8swXOBERkYGUws4omy4OHjyI6Oho+Pn5QSaT4fvvvy/z/C1btiA8PBzVq1eHm5sbQkNDsWfPHgPu2jisJtGZNWsW/P39sWrVKjz//POoU6cOXn31VdSrVw/Ao96cxMRETJo0CW+++SaCgoKwZs0a3L9/H+vXrzdz9ERERPp7PHRl6KaLe/fuoXnz5li0aJFW5x88eBDh4eHYuXMnkpOT0b59e0RHRyMlJUWfWzYaq5mjs23bNkRGRqJbt244cOAAatasiZiYGAwZMgQAcOXKFWRmZiIiIkK6Ri6Xo127djhy5AiGDh1aarv5+fnIz8+X3ufm5pr2RoiIiKxAVFQUoqKitD4/MTFR7f2MGTPwww8/4Mcff0RwcLCRo9Oe1SQ6ly9fxpIlSzB69Gj83//9H37//XeMHDkScrkcffv2RWZmJgDA29tb7Tpvb29cu3ZNY7sJCQmYOnWqUWM9+MM4rc4rXp2l7XXatMfFPkunuuBnWQuBGlqdZa7qJ2Po6PWu9Hr3zeU6X69WfaZlNVWJRT1VqrgqPVNX5xjI+pVXlZO5q6m0pYThVVPKf/9b/A96uVwOuVxuUNulfp5Siby8PFSrVs3obevCaoaulEolWrZsiRkzZiA4OBhDhw7FkCFDsGTJErXzZDL1HwQhRIl9qiZOnIicnBxpS0tLM0n8RERE+nr8wEBDNwDw9/eHu7u7tCUkJJgk5nnz5uHevXvo3r27SdrXltX06Pj6+qJx48Zq+xo1aoTvvvsOAODj4wMAyMzMhK+vr3ROVlZWiV4eVabKZImIiIzFOEtAPLo+LS0Nbm5u0n5T/A7csGED4uLi8MMPP6BGjRpGb18XVtOj07ZtW1y4cEFt38WLFxEQEAAACAwMhI+PD5KSkqTjBQUFOHDgANq0aVOusRIREVkqNzc3tc3Yic6mTZswaNAgfPPNN+jQoYNR29aH1fTofPDBB2jTpg1mzJiB7t274/fff8fy5cuxfPmjOQQymQyxsbGYMWMG6tevj/r162PGjBlwdnZGr169zBw9ERGR/pSQQQlD5+iY/snIGzZswMCBA7FhwwZ06tTJ5J+nDatJdJ577jls3boVEydOxLRp0xAYGIjExET07t1bOmf8+PF48OABYmJikJ2djdatW2Pv3r1wdXU1Y+RERESGMebQlbbu3r2LS5cuSe+vXLmCU6dOoVq1aqhduzYmTpyI9PR0rF37qHhjw4YN6Nu3Lz799FO88MILUpGQk5MT3N3dDYrdEDIhbHTddg1yc3Ph7u6OnJwctTFMa1cRKrK0rYQqXkFVXp9bURhadaWqeNVV8eoq6byy1gZTWS+rUr06BsVDlsWUFU9lrWFlzM819e+Mx+0vONEGTi6G9U08uFuED1od0TrW/fv3o3379iX29+vXD6tXr0b//v1x9epV7N+/HwAQFhaGAwcOaDzfXKymR4eIiMhWGWetK92uDwsLQ1l9IcWTl8cJj6VhokNERGThlEIGpaHP0bHR1cutpuqKiIiISFfs0SEiIrJwSiMMXSlttG+DiQ4REZGF02f18dLasEW2eddERERkE9ijYyNUFww1dvl1eSkrblOWgGtaCLSiKF7arboop2qpeXHK3LvSa/tqVUs9p1J1T7X3mr6Hdk5O6m0/zC/1vKK/rqq3b4Zyc8VV466HZ1/H36jtlZcdB7dqPNbp5Te0aqOsEnBVhpaDW8vCnWVRQAaFgQ/8M/R6a8VEh4iIyMJx6Ep/THSIiIgsnAKG98gojBOK1bHN9I6IiIhsAnt0iIiILByHrvTHRIeIiMjCmWNRz4qCiY4NKqtyqCJWZGmjIlZTaUu1ykoXdm4uTz2n6J9bmq93lGs8Zl/FWasYRGZWqftlPjW0ut7YFVT60DYG1eqsokuXjRpDpWfqajym+ll7rv9Xq/ZUK7K0rcAytopQaUXGwUSHiIjIwgnIoDRwMrJgeTkRERFZIg5d6c8275qIiIhsAnt0iIiILJxSyKAUhg09GXq9tWKiQ0REZOEURli93NDrrRUTHVLzoENTc4egNaefTps7BPqXanWVvbtbqa8BQHHnjvRadT0rmb292nmisOjJ66JC6XXxCrGObgOk17tzV0mvtV1DSdXegvVq7yMceunchrHtTj2h8u7J60i/5kb9HGNXcanSdk2s4lVS+nwPiUrDRIeIiMjCcehKf0x0iIiILJwSdlAaOPRk6PXWiokOERGRhVMIGRQG9sgYer21ss30joiIiGwCe3SIiIgsHOfo6I+JDhERkYUTRli9XNjok5GZ6JDVUi2FVy0117ZE/uAP44wekzVSLdEG1Mu0tZWk2PSkvWpDdL5eKBSa2y5j0VF9YlWlWlJeVjm5rLJ2/1SqlsXrQ72cXLOyFtc0Rum5tot36qOsknI1RZp/Joh0wUSHiIjIwikgg8LARTkNvd5aMdEhIiKycEph+BwbpTBSMFbGNgfsiIiIyCawR4eIiMjCKY0wGdnQ660VEx0iIiILp4QMSgPn2Bh6vbWy2kQnISEB//d//4dRo0YhMTERACCEwNSpU7F8+XJkZ2ejdevW+Pzzz9GkSRPzBksmV1alFauryqa4d1/tffEqrMe0rXDafXuFxmOqFVmqC3xCpv6XpmoVl6oonxitYiiL6gKk4fY9NH6m6rGyqsLsHOXS67Kqs1QXMVW1Jz1Fc7B6MGXFlCo7n4tq75WZDaTXUXWeexLPw3Vq5+1Sv0yjXZfn6h9cBcQnI+vPKvuxjh8/juXLl6NZs2Zq+2fPno358+dj0aJFOH78OHx8fBAeHo68vDwzRUpERETmZHWJzt27d9G7d2+sWLECHh4e0n4hBBITEzFp0iS8+eabCAoKwpo1a3D//n2sX7++jBaJiIgs2+M5OoZutsjq7vr9999Hp06d0KFDB7X9V65cQWZmJiIiIqR9crkc7dq1w5EjRzS2l5+fj9zcXLWNiIjIkighk5aB0HvjHB3Lt3HjRpw8eRLHjx8vcSwzMxMA4O3trbbf29sb165d09hmQkICpk6datxAiYiIyCJYTY9OWloaRo0aha+//hqOjo4az5PJ1DNWIUSJfaomTpyInJwcaUtLSzNazERERMYg/q26MmQT7NGxbMnJycjKykJISIi0T6FQ4ODBg1i0aBEuXLgA4FHPjq+vr3ROVlZWiV4eVXK5HHK5XONxIiIic+Pq5fqzmkTn1VdfxenTp9X2DRgwAA0bNsSECRNQt25d+Pj4ICkpCcHBwQCAgoICHDhwALNmlbFwHBGZTPEFPjWVnquWchdnjJJyVZWqexp0nmp5OqBeNq5aal4WY5eUm4NqOTmgXm4ukw/Uub2oumPV3rO8nIzFahIdV1dXBAUFqe2rUqUKPD09pf2xsbGYMWMG6tevj/r162PGjBlwdnZGr16aVyUmIiKydHwysv6sJtHRxvjx4/HgwQPExMRIDwzcu3cvXF1dzR0aERGR3jh0pT+rTnT279+v9l4mkyEuLg5xcXFmiYeIiIgsi1UnOkRERLaAa13pzzYH7IiIiKyIwQ8L1GPo6+DBg4iOjoafnx9kMhm+//77p15z4MABhISEwNHREXXr1sXSpUv1vGPjYY8OqVFdAPPlLnPMGIn+uIinboovZhlZpe+TNyrPoIp06adVe3vurpFel7XAZ1kxqNqVuVirNiIcSi86kDlUVntv51Kl1PO0re7StrKqRBwhTxYXvqs8Kr12sdP8XDBLVnxRT1W7c77UeKx4ddVjrLIqmznm6Ny7dw/NmzfHgAED8NZbbz31/CtXruC1117DkCFD8PXXX+PXX39FTEwMqlevrtX1pqJVotOyZUudGpXJZNi2bRtq1qypV1BERERkXlFRUYiKitL6/KVLl6J27dpITEwEADRq1AgnTpzA3LlzLT/ROXXqFMaMGQMXF5enniuEwMyZM5Gfn//Uc4mIiOjpjNmjU3xNR2M9OPfo0aNq600CQGRkJFauXInCwkJUrlxZw5WmpfXQ1bhx41CjRg2tzp03b57eAREREZE6YyY6/v7+avunTJlilGrlzMzMUtebLCoqws2bN9VWLShPWiU6V65cQfXq1bVu9Ny5c/Dz89M7KCIiIjKNtLQ0uLm5Se+NuQxSaetNlra/PGmV6AQEBOjUaPFskYiIiPQnYHh5uPj3v25ubmqJjrH4+PggMzNTbV9WVhYqVaoET0/tll4xBb2qru7cuYPff/8dWVlZUCqVasf69u2r4SqyNqrVS5oqJQDgXtMnXZVVTt8odb+psdJKf2pVVkZQ1rpVmqqril9TyffJz86uvz/T6nNllbX750x59570WlMFljGIJs+YrG1LUNZaV6q4hpVxWMOTkUNDQ/Hjjz+q7du7dy9atWpltvk5gB6Jzo8//ojevXvj3r17cHV1VeuOkslkTHSIiIgqgLt37+LSpUvS+ytXruDUqVOoVq0aateujYkTJyI9PR1r164FAAwbNgyLFi3C6NGjMWTIEBw9ehQrV67Ehg0bzHULAPR4YOCYMWMwcOBA5OXl4c6dO8jOzpa227dvmyJGIiIim2aOBwaeOHECwcHBCA4OBgCMHj0awcHBmDx5MgAgIyMDqamp0vmBgYHYuXMn9u/fjxYtWuCTTz7BZ599ZtbSckCPHp309HSMHDkSzs7OpoiHiIiIijHH0FVYWJg0mbg0q1evLrGvXbt2OHnypK6hmZTOPTqRkZE4ceKEKWIhIiIiMiqtenS2bdsmve7UqRPGjRuHc+fOoWnTpiUmGHXu3Nm4ERIREdk4a5iMbKm0SnS6du1aYt+0adNK7JPJZFAoFAYHRURERE8IIYMwMFEx9HprpVWiU7yEnGxPWSWhmkrPVUvNi1MtPde3NNxaFx21BJoW7jS1skrPNYmqNVLjMa1Kz8uYY6C6YKjWi3q6qi+Fo6jto9V1337/Ran77yofSq+71WunHt+V37RqW7UNS1sklOXkxqGEzODn6Bh6vbXSeY7O2rVrS13HqqCgQCoxIyIiIrIEOic6AwYMQE5OTon9eXl5GDBggFGCIiIioifMUV5eUehcXi6EKHXNir///hvu7u5GCYqIiIie4Bwd/Wmd6AQHB0Mmk0Emk+HVV19FpUpPLlUoFLhy5Qo6duxokiCJiIiI9KF1ovO48urUqVOIjIyEi8uTCXkODg6oU6eO2Z9+SEREVBGxvFx/Wic6U6ZMgUKhQEBAACIjI+Hr62vKuMiKqFZVlLX4p6qyKrK0xYU8daNaaaV88EDjeXZW9NTzsiqytKH6NdlzT7tiCm2rs/Sx+a8Dau/v6lHw+tqzL0qvd144bGhIGmlaxJNMg0NX+tNpMrK9vT2GDRuGhw8fPv1kIiIiIjPTueqqadOmuHz5siliISIiolIII1RcsUdHS/Hx8Rg7diy2b9+OjIwM5Obmqm1ERERkXAKPnn1p0GbumzATncvLH1dWde7cWa3M/HHZOZeAICIiIkuhc6Lzyy+/mCIOIiIi0kAJGWRcAkIvOic67dq1e/pJZLO0rcAy9vo3EZV7Sq/3Fm40atsVRVmVVqpEQWGp+2UOlbW63q6y5n9WhEqPr1AatyNdtWpKbS0vLRW/RtPXK0m5Wav2cq/X1jkGbbn5paq97+g+UKvrOoVElrp/R/IejdewusoysOpKfzonOgBw584drFy5EufPn4dMJkPjxo0xcOBAPhmZiIjIBJRCBhmfo6MXnScjnzhxAvXq1cOCBQtw+/Zt3Lx5E/Pnz0e9evVw8uRJU8RIREREpBede3Q++OADdO7cGStWrJCWgSgqKsLgwYMRGxuLgwcPGj1IIiIiW/a4csrQNmyRzonOiRMn1JIcAKhUqRLGjx+PVq1aGTU4IiIi4hwdQ+g8dOXm5obU1NQS+9PS0uDq6mqUoEqTkJCA5557Dq6urqhRowa6du2KCxcuqJ0jhEBcXBz8/Pzg5OSEsLAwnD171mQxERERkWXTOdHp0aMHBg0ahE2bNiEtLQ1///03Nm7ciMGDB+Ptt982RYwAgAMHDuD999/HsWPHkJSUhKKiIkRERODevXvSObNnz8b8+fOxaNEiHD9+HD4+PggPD0deXp7J4iIiIjK1xz06hm62SOehq7lz50Imk6Fv374oKioCAFSuXBnvvfceZs6cafQAH9u9e7fa+1WrVqFGjRpITk7Gyy+/DCEEEhMTMWnSJLz55psAgDVr1sDb2xvr16/H0KFDTRYb6U619LysUnNtz9O2pDzcvof0OkmxSatr6JHiZedq5eb//lugC5md5n90lXdydG5PlaGl5sVpW1KuqngJuD7l5sXb0NVrQWFq7xXZWdJr1Z9/ZWYDtfNYUm55WHWlP517dBwcHPDpp58iOzsbp06dQkpKCm7fvo0FCxZALpebIsZS5eQ8+oewWrVqAIArV64gMzMTERER0jlyuRzt2rXDkSNHNLaTn5/PZSyIiIgqKL2eowMAzs7OaNq0qTFj0ZoQAqNHj8aLL76IoKAgAEBmZiYAwNvbW+1cb29vXLt2TWNbCQkJmDp1qumCJSIiMhCrrvSnc6Jz7949zJw5Ez///DOysrKgVCrVjpfHyubDhw/H//73Pxw+fLjEMdX1t4Ana3BpMnHiRIwePVp6n5ubC39/f+MFS0REZKBHiY6hVVdGCsbK6JzoDB48GAcOHECfPn3g6+tbZhJhCiNGjMC2bdtw8OBB1KpVS9rv4+MD4FHPjq+vr7Q/KyurRC+PKrlcXq5DbkRERFR+dE50du3ahR07dqBt27amiEcjIQRGjBiBrVu3Yv/+/QgMDFQ7HhgYCB8fHyQlJSE4OBgAUFBQgAMHDmDWrFnlGisREZEx8Tk6+tM50fHw8JAmAJen999/H+vXr8cPP/wAV1dXaU6Ou7s7nJycIJPJEBsbixkzZqB+/fqoX78+ZsyYAWdnZ/Tq1avc4yXtq6nKWvyz6Krm+VX6sKVKq3C7bqb9AD0qrUxJY0WdnouH6lNpZShDq6yext7jyXqE6l8v3ausovxHqb3flfap/oHRU4l/N0PbsEU6JzqffPIJJk+ejDVr1sDZ2dkUMZVqyZIlAICwsDC1/atWrUL//v0BAOPHj8eDBw8QExOD7OxstG7dGnv37jXpgwyJiIhMjT06+tM50Zk3bx7++usveHt7o06dOqhcubLacVMt7Cm0mEUlk8kQFxeHuLg4k8RARERE1kXnRKdr164mCIOIiIg04tiV3nROdKZMmaLVeRs2bEDnzp1RpUoVnYMiIiIiFcZYwsFGh650fjKytoYOHYobN26YqnkiIiKip9L7ychPo82cGiIiIno6PhlZfyZLdIjKUlbpuSrVEljV14BtlYobg8zeXuWN7p25JRbhbPLMk9dnL+kZ1b9tG+GhnZVqPnlQaFTtWK2u2fPgK4M/V1uqpeP6LPBZlt05X5a6v6PXu2rvlbl3pddqPw96YDl5+WLVlf5MNnRFREREZG7s0SEiIrJ0Qmb4ZGIb7dFhokNERGThOEdHfyYbugoICCjxMEEiIiKi8qRzj05aWhpkMpm0cvjvv/+O9evXo3Hjxnj33ScT386cOWO8KImIiGwZHxioN50TnV69euHdd99Fnz59kJmZifDwcDRp0gRff/01MjMzMXnyZFPESTZKtbJK30UqO7oPlF6LgsJSzylefaP6WeZY3NEYiscdUbnnkzdC+eR1GRVYqpU5okj9a6c22q9agaVKy2os5f37au/ttFxHz66ax9Pbzn+o8ZjGhUAriN03l2s8FlVzRDlGQoZi1ZX+dB66OnPmDJ5//nkAwDfffIOgoCAcOXIE69evx+rVq40dHxEREQFPenX03fSwePFiBAYGwtHRESEhITh06FCZ569btw7NmzeHs7MzfH19MWDAANy6dUu/DzcSnROdwsJCyP995sVPP/2Ezp07AwAaNmyIjIwM40ZHREREZrFp0ybExsZi0qRJSElJwUsvvYSoqCikpqaWev7hw4fRt29fDBo0CGfPnsXmzZtx/PhxDB48uJwjV6dzotOkSRMsXboUhw4dQlJSEjp27AgAuH79Ojw9PY0eIBERka17PHRl6KaL+fPnY9CgQRg8eDAaNWqExMRE+Pv7Y8mSJaWef+zYMdSpUwcjR45EYGAgXnzxRQwdOhQnTpwwxpdAbzonOrNmzcKyZcsQFhaGt99+G82bNwcAbNu2TRrSIiIiIiMydNhKZfgqNzdXbcvPzy/xcQUFBUhOTkZERITa/oiICBw5cqTUENu0aYO///4bO3fuhBACN27cwLfffotOnToZevcG0XkyclhYGG7evInc3Fx4eDyZCPjuu+/CWcsJhERERGQe/v7+au+nTJmCuLg4tX03b96EQqGAt7e32n5vb29kZmaW2m6bNm2wbt069OjRAw8fPkRRURE6d+6MhQsXGjV+Xen1wEAhBJKTk/HXX3+hV69ecHV1hYODAxMdMiltq59Uq6y0FenUR+29ndxR5zYs3d7CjdJr1aoyOwft1jwSSvXZjPrUbwiFQrvzSvkLU1+VatXUeKwo3TzzClXXvTKXXenG/eUT5T/qSdtcB8sEZNDv/7ribTx6TIybm5u0V17GWnMymfpnCiFK7Hvs3LlzGDlyJCZPnozIyEhkZGRg3LhxGDZsGFauXGlg7PrTOdG5du0aOnbsiNTUVOTn5yM8PByurq6YPXs2Hj58iKVLl5oiTiIiIttlxOfouLm5qSU6pfHy8oK9vX2J3pusrKwSvTyPJSQkoG3bthg3bhwAoFmzZqhSpQpeeuklTJ8+Hb6+vqVeZ2o6z9EZNWoUWrVqhezsbDg5OUn733jjDfz8889GDY6IiIjKn4ODA0JCQpCUlKS2PykpCW3atCn1mvv378POTj2tsP/3WVxCh/Un7O3tkZWVVWL/rVu3pPZ0oXOPzuHDh/Hrr7/CwcFBbX9AQADS09N1DoCIiIiewgxPRh49ejT69OmDVq1aITQ0FMuXL0dqaiqGDRsGAJg4cSLS09Oxdu1aAEB0dDSGDBmCJUuWSENXsbGxeP755+Hn56d9mBqSovz8/BK5hzZ0TnSUSiUUpYyz//3333B1ddU5ACIiInoKM6xe3qNHD9y6dQvTpk1DRkYGgoKCsHPnTgQEBAAAMjIy1J6p079/f+Tl5WHRokUYM2YMqlatildeeQWzZs3S6vM+++wzAI/mBX3xxRdwcXGRjikUChw8eBANGzbU6R4APRKd8PBwJCYmYvny5VJAd+/exZQpU/Daa6/pHAARERFZppiYGMTExJR6rLTVEEaMGIERI/RbXmTBggUAHvXoLF26VG2YysHBAXXq1NFrHrDOic6CBQvQvn17NG7cGA8fPkSvXr3w559/wsvLCxs2bNA5ACIiIiqbEI82Q9uwZFeuXAEAtG/fHlu2bFF7hI0hdE50/Pz8cOrUKWzYsAEnT56EUqnEoEGD0Lt3b7XJyUSWSOZQWXqtaYHP4spaTNReZVHJshZQtDSqpfqRjr3Vju15uE7n9iJDpjx5o+VCnpVqVJdeF2X9o3bM3rOa9Fpx67bGNkROrvRa5l56FUlZJeQVcSHPV19JkF7/vG+iwe0Vf/TCY8UXwmVJuYnZ0Orlv/zyi9p7hUKB06dPIyAgQK/kR6/n6Dg5OWHgwIEYOFD355UQERGRjswwR8dcYmNj0bRpUwwaNAgKhQIvv/wyjh49CmdnZ2zfvh1hYWE6tadzeTkAfPXVV3jxxRfh5+eHa9euAXg0pPXDDz/o0xwRERERAGDz5s3S8lI//vgjrl69ij/++ENaYFRXOic6S5YswejRoxEVFYXs7GypAsvDwwOJiYk6B0BERERlkwnjbNbg1q1b8PHxAQDs3LkT3bp1Q4MGDTBo0CCcPn1a5/Z0TnQWLlyIFStWYNKkSahU6cnIV6tWrfQKgIiIiJ7CiIt6Wjpvb2+cO3cOCoUCu3fvRocOHQA8eiBhuTww8MqVKwgODi6xXy6X4969ezoHQERERPTYgAED0L17d/j6+kImkyE8PBwA8Ntvv5XPc3QCAwNx6tQp6YFBj+3atQuNGzfWOQAiswl65snrM5c0H0s+I720CwlSO81K/kAqkz5VViXaSJ4qvS6rSk0T1Qqs4lQrsMT9B+oHVeplVSuwylIRK61UqVZaqVZgAUClo+ek16pVU2WdVxT65N91Y1RxkZ5saDJyXFwcgoKCkJaWhm7dukmLjtrb2+PDDz/UuT2dE51x48bh/fffx8OHDyGEwO+//44NGzYgISEBX3zxhc4BEBER0VPYUHk5APznP/8BADx8+FDa169fP73a0nmOzoABAzBlyhSMHz8e9+/fR69evbB06VJ8+umn6Nmzp15BGNvixYsRGBgIR0dHhISE4NChQ+YOiYiIiLSgUCjwySefoGbNmnBxccHly5cBAB9//DFWrlypc3s6JTpFRUVYs2YNoqOjce3aNWRlZSEzMxNpaWkYNGiQzh9uCps2bZJK0FJSUvDSSy8hKipKbT0OIiIiq2JDk5Hj4+OxevVqzJ49W20Rz6ZNm+o1cqRTolOpUiW89957yM/PBwB4eXmhRo0aOn+oKc2fPx+DBg3C4MGD0ahRIyQmJsLf3x9Lliwxd2hERET6saFEZ+3atVi+fDl69+6tVmXVrFkz/PHHHzq3p/PQVevWrZGSkqLzB5WHgoICJCcnIyIiQm1/REQEjhw5Uuo1+fn5yM3NVduIiIjIPNLT0/HMM8+U2K9UKlFYqN3SPap0nowcExODMWPG4O+//0ZISAiqVKmidrxZs2Y6B2EsN2/ehEKhgLe3t9p+b29vZGZmlnpNQkICpk6dWuoxsh4RDr00HrNzcnx6A0El/6eSri9WaaVJ5HOaf472HJ+i8RjpoKxnaCiVpe6u6FVWZSleJVW8ukrX88o6R9uKLGOvxWUzbKjqqkmTJjh06FCJ6u7NmzeX+nibp9E50enRowcAYOTIkdI+mUwGIQRkMpn0pGRzksnUv5mPYyvNxIkTMXr0aOl9bm4u/P39TRofERGRLozxZGNreTLylClT0KdPH6Snp0OpVGLLli24cOEC1q5di+3bt+vcnl4PDLRUXl5esLe3L9F7k5WVVaKX5zG5XC7V6BMREVkkGyovj46OxqZNmzBjxgzIZDJMnjwZLVu2xI8//ig9PFAXOic6xbuSLImDgwNCQkKQlJSEN954Q9qflJSELl26mDEyIiIi0lZkZCQiIyPLPGfDhg3o3LlziSk0xemc6Gzbtq3U/TKZDI6OjnjmmWcQGBioa7NGM3r0aPTp0wetWrVCaGgoli9fjtTUVAwbNsxsMREREZFxDR06FK1bt0bdunXLPE/nRKdr167SnBxVqvN0XnzxRXz//ffw8PDQtXmD9ejRA7du3cK0adOQkZGBoKAg7Ny506J7ooiIiMoigxHm6BglEstRPA/RROfy8qSkJDz33HNISkpCTk4OcnJykJSUhOeffx7bt2/HwYMHcevWLYwdO1bnoI0lJiYGV69eRX5+PpKTk/Hyyy+bLRYiIiIyH517dEaNGoXly5ejTZs20r5XX30Vjo6OePfdd3H27FkkJiZi4MCBRg2USFuyYiXIokDluQtllJGbkqbS84pYdp6k3Cy9Lr7AZ1mLdxrM7snfbbuzue6eIbQp+9a3NLysRUdV2R04Jb225UcESGyovNzYdE50/vrrL7i5uZXY7+bmJq1HUb9+fdy8edPw6IiIiMimqq6MTeehq5CQEIwbNw7//POPtO+ff/7B+PHj8dxzzwEA/vzzT9SqVct4URIRERHpQecenZUrV6JLly6oVasW/P39IZPJkJqairp16+KHH34AANy9excff/yx0YMlIiKySezRKSEgIACVK1d+6nk6JzrPPvsszp8/jz179uDixYsQQqBhw4YIDw+H3b9j5F27dtU5YCIiIiqdLT0Z+bETJ07g/PnzkMlkaNiwIVq1aqV2/MyZM1q1o3OiAzwqJe/YsSPCwsIgl8s1Lq9AREREpIu///4bb7/9Nn799VdUrVoVAHDnzh20adMGGzZs0HmZJp0THaVSifj4eCxduhQ3btzAxYsXUbduXXz88ceoU6cOBg0apGuTRAYrXmmlxkyVVtooXo1V0aqw7OTqC6oqc/JKP8/dtTzCsXmaKqXKqn4y5cKb2lZghdv3kF7bbAWWDQ1dDRw4EIWFhTh//jyeffZZAMCFCxcwcOBADBo0CHv37tWpPZ0nI0+fPh2rV6/G7Nmz4eDgIO1v2rQpvviCJZ1ERERGJ4y0WYFDhw5hyZIlUpIDPJo2s3DhQhw6dEjn9nROdNauXYvly5ejd+/esFf5K7pZs2b4448/dA6AiIiIyvZ4jo6hmzWoXbs2CgsLS+wvKipCzZo1dW5P50QnPT0dzzxTcihAqVSWGhgRERGRtmbPno0RI0bgxIkT0jIPJ06cwKhRozB37lyd29N5jk6TJk1w6NChEmtHbd68GcHBwToHQERERE9hQ09G7t+/P+7fv4/WrVujUqVHaUpRUREqVaqEgQMHqq28cPv27ae2p3OiM2XKFPTp0wfp6elQKpXYsmULLly4gLVr12L79u26NkdERERPY0OTkRMTE43ans6JTnR0NDZt2oQZM2ZAJpNh8uTJaNmyJX788UeEh4cbNTgiIiKyLf369TNqe3o9RycyMhKRkZFGDYTIEHsefCW9jnTqo37wzKUnr1VKzcWp89JrWYtGJouNjEd5777ae7sqzmaKpGIxZQm5NcVgyWztgYEKhQLff/+99MDAxo0bo3PnzmpFUNrSK9EhIiKicmRDQ1eXLl3Ca6+9hvT0dDz77LMQQuDixYvw9/fHjh07UK9ePZ3a0yrR8fDw0Prpx9pMDCIiIiIqzciRI1GvXj0cO3YM1apVAwDcunUL77zzDkaOHIkdO3bo1J5WiY7qxKBbt25h+vTpiIyMRGhoKADg6NGj2LNnDxfyJCIiMgVjPAfHSnp0Dhw4oJbkAICnpydmzpyJtm3b6tyeVomO6sSgt956C9OmTcPw4cOlfSNHjsSiRYvw008/4YMPPtA5CCIiIiqDDQ1dyeVy5OWVXC7m7t27aisyaEvnBwbu2bMHHTt2LLE/MjISP/30k84BEBERET32+uuv491338Vvv/0GIQSEEDh27BiGDRuGzp0769yezpORPT09sXXrVowbN05t//fffw9PT0+dAyAyNtUKLKBYFZZqBZYK1Qqs4kxZkVXRFvHUl+L2Hb2uU63CUl3YVfV7XvznQVWkY+8n5z1cp1cMVDbVr3FxZX3N+b0pxoZ6dD777DP069cPoaGhqFy5MgCgsLAQXbp0waeffqpzezonOlOnTsWgQYOwf/9+aY7OsWPHsHv3bi7qSUREZAK2VF5etWpV/PDDD7h06RLOnTsHAGjcuHGpy09pQ+dEp3///mjUqBE+++wzbNmyBUIING7cGL/++itat26tVxBEREREj61cuRILFizAn3/+CQCoX78+YmNjMXjwYJ3b0nmODgC0bt0a69atw8mTJ5GSkoJ169YxySEiIqpgFi9ejMDAQDg6OiIkJASHDh0q8/z8/HxMmjQJAQEBkMvlqFevHr788kudPvPjjz/GqFGjEB0djc2bN2Pz5s2Ijo7GBx98gI8++kjne9CqRyc3Nxdubm5aN5qXlwdXV1edgyEiIqJSmGGOzqZNmxAbG4vFixejbdu2WLZsGaKionDu3DnUrl271Gu6d++OGzduYOXKlXjmmWeQlZWFoqIinT53yZIlWLFiBd5++21pX+fOndGsWTOMGDEC06dP16k9rR8YmJGRgRo1amjVaM2aNXHq1CnUrVtXp2CIiIioJHPM0Zk/fz4GDRokDRclJiZiz549WLJkCRISEkqcv3v3bhw4cACXL1+WnoFTp04dneNUKBRo1apVif0hISE6J02AlomOEAJffPEFXFxctGq0sLBQ50CITKXMdbC0wDWxzEwoNR6SVapc6n5tK63IePT5umpbWcUKLOPKzc1Vey+XyyGXy9X2FRQUIDk5GR9++KHa/oiICBw5cqTUdrdt24ZWrVph9uzZ+Oqrr1ClShV07twZn3zyCZycnLSO75133sGSJUswf/58tf3Lly9H7966/5xplejUrl0bK1as0LpRHx8fqSSMiIiIjMBIVVP+/v5q76dMmYK4uDi1fTdv3oRCoYC3t7fafm9vb2RmZpba7uXLl3H48GE4Ojpi69atuHnzJmJiYnD79m2d5+msXLkSe/fuxQsvvADgUXV3Wloa+vbti9GjR0vnFU+GSqNVonP16lWdAiQiIiIjMuIcnbS0NLV5t8V7c1QVX+dSCKFx7UulUgmZTIZ169bB3d0dwKNE5D//+Q8+//xzrXt1zpw5g5YtWwIA/vrrLwBA9erVUb16dZw5c0ZjbJpw9XIiIiIb4ubm9tQCIy8vL9jb25fovcnKyirRy/OYr68vatasKSU5ANCoUSMIIfD333+jfv36WsX3yy+/aHWetvQqLyciIqLy83gysqGbthwcHBASEoKkpCS1/UlJSWjTpk2p17Rt2xbXr1/H3bt3pX0XL16EnZ0datWqpdd9GwMTHSIiIksnjLTpYPTo0fjiiy/w5Zdf4vz58/jggw+QmpqKYcOGAQAmTpyIvn37Suf36tULnp6eGDBgAM6dO4eDBw9i3LhxGDhwoE6TkY2NQ1dERERUQo8ePXDr1i1MmzYNGRkZCAoKws6dOxEQEAAAyMjIQGpqqnS+i4sLkpKSMGLECLRq1Qqenp7o3r27zs+9MTarSHSuXr2KTz75BPv27UNmZib8/PzwzjvvYNKkSWpLtqempuL999/Hvn374OTkhF69emHu3Ll6LetOFZNq2XFE5Z5mjMQ2KPMfGtyG6mKdZSmrpFzjNSxV1puxy/RZ9l82c611FRMTg5iYmFKPrV69usS+hg0blhjuMje9Ep1Dhw5h2bJl+Ouvv/Dtt9+iZs2a+OqrrxAYGIgXX3zR2DHijz/+gFKpxLJly/DMM8/gzJkzGDJkCO7du4e5c+cCePSAoU6dOqF69eo4fPgwbt26hX79+kEIgYULFxo9JiIionJjQ6uXG5vOc3S+++47REZGwsnJCSkpKcjPzwfwaNmHGTNmGD1AAOjYsSNWrVqFiIgI1K1bF507d8bYsWOxZcsW6Zy9e/fi3Llz+PrrrxEcHIwOHTpg3rx5WLFiRYmHIxEREVkVM8zRqSh0TnSmT5+OpUuXYsWKFWoPBWzTpg1Onjxp1ODKkpOTIz1iGgCOHj2KoKAg+Pn5SfsiIyORn5+P5ORkje3k5+cjNzdXbSMiIqKKQedE58KFC3j55ZdL7Hdzc8OdO3eMEdNT/fXXX1i4cKE08xsAMjMzS9T2e3h4wMHBQeNTHAEgISEB7u7u0lb8iZFERETmVt7l5RWJzomOr68vLl26VGL/4cOHdV7EMy4uDjKZrMztxIkTatdcv34dHTt2RLdu3aSFxh4r7SmJZT3FEXhUHpeTkyNtaWlpOt0DERGRyXHoSm86T0YeOnQoRo0ahS+//BIymQzXr1/H0aNHMXbsWEyePFmntoYPH46ePcuufFFd+fT69eto3749QkNDsXz5crXzfHx88Ntvv6nty87ORmFhocanOAKlL2ZGJJSl/4ug3QPHy7bn+BQjtGIdkpSb1d6H23XTuQ2hUEiv7eSOOl/Pah7DWNrXr3g8rJyjp9E50Rk/fjxycnLQvn17PHz4EC+//DLkcjnGjh2L4cOH69SWl5cXvLy8tDo3PT0d7du3R0hICFatWgU7O/XOqNDQUMTHxyMjIwO+vr4AHk1QlsvlCAkJ0SkuIiIii8KqK73pVV4eHx+PSZMm4dy5c1AqlWjcuDFcXFyMHZvk+vXrCAsLQ+3atTF37lz8888/0jEfHx8Aj5aOb9y4Mfr06YM5c+bg9u3bGDt2LIYMGfLUNT2IiIgsmbmeo1MR6P3AQGdnZ7Rq1cqYsWi0d+9eXLp0CZcuXSqxXoYQj75z9vb22LFjB2JiYtC2bVu1BwYSERGRbdIq0XnzzTe1blD12TbG0r9/f/Tv3/+p59WuXRvbt283+ucTERGZFYeu9KZVoqO65LoQAlu3boW7u7vUo5OcnIw7d+7olBARERGRdjh0pT+tEp1Vq1ZJrydMmIDu3btj6dKlsP93DRqFQoGYmBjOhSGLF+HQS3ptp8dquoqT59Te27VsXOp5tlRZZWr2Kn9oFbc7+4unXs+qnIpNtQqL32sqjc7P0fnyyy8xduxYKckBHs2PGT16NL788kujBkdERETgc3QMoHOiU1RUhPPnz5fYf/78eSiVSqMERURERCqY6OhN56qrAQMGYODAgbh06RJeeOEFAMCxY8cwc+ZMDBgwwOgBEhER2ToZDH9gqTEeeGqNdE505s6dCx8fHyxYsAAZGRkAHi0LMX78eIwZM8boARIRERHpS+dEx87ODuPHj8f48eOllb45CZmIiMiEWF6uN70fGAgwwSEiIioPLC/Xn86JTmBgYJmrgV++fNmggIgsmaZyctKOaqm4Iien1P1ERMakc6ITGxur9r6wsBApKSnYvXs3xo0bZ6y4iIiI6DEOXelN50Rn1KhRpe7//PPPceLECYMDIiIiolLYaKJiKJ2fo6NJVFQUvvvuO2M1R0RERGQwgyYjq/r2229RrVo1YzVHRERE/+JkZP3pnOgEBwerTUYWQiAzMxP//PMPFi9ebNTgiIiICJyjYwCdE50uXbqoJTp2dnaoXr06wsLC0LBhQ6MGR2RplMln1N7bhQRJr7mQ59OpLsLZ0WOwQddT+VBdKFN1AU0ia6FzohMXF2eCMIiIiEgTDl3pT+fJyPb29sjKyiqx/9atW2ormhMREZGRcFFPvencoyNE6V+p/Px8ODg4GBwQERERqWOPjv60TnQ+++wzAIBMJsMXX3wBFxcX6ZhCocDBgwc5R4eIiIgsitaJzoIFCwA86tFZunSp2jCVg4MD6tSpg6VLlxo/QiIiIlvHqiu9aZ3oXLlyBQDQvn17bNmyBR4eHiYLioiIiFQw0dGbznN0fvnlF1PEQWSVWFJuHCwbtz6qZefGoG/purHjoIpHq0Rn9OjR+OSTT1ClShWMHj26zHPnz59vlMCIiIjoEU5G1p9WiU5KSgoKCwsBACdPnlR7YCARERGZGIeu9KZVoqM6XLV//35TxUJERERkVDo/MHDgwIHIy8srsf/evXsYOHCgUYIiIiKiJ2RCGGWzRTonOmvWrMGDBw9K7H/w4AHWrl1rlKCIiIhIBZ+MrDetq65yc3MhhIAQAnl5eXB0dJSOKRQK7Ny5EzVq1DBJkESGiHDoJb22c5Qb1FaScrOh4dC/WGllfUxZ4VS8bUtbQDTcrpv0mv8OWBetE52qVatCJpNBJpOhQYMGJY7LZDJMnTrVqMERERERq64MoXWi88svv0AIgVdeeQXfffcdqlWrJh1zcHBAQEAA/Pz8TBIkERGRTWPVld60TnTatWsH4NETkv39/WFnp/P0HiIiItIDe3T0p3O2EhAQADs7O9y/fx9//PEH/ve//6ltppafn48WLVpAJpPh1KlTasdSU1MRHR2NKlWqwMvLCyNHjkRBQYHJYyIiIiLLpHOi888//+D111+Hq6srmjRpguDgYLXN1MaPH1/qEJlCoUCnTp1w7949HD58GBs3bsR3332HMWPGmDwmIiIikzJT1dXixYsRGBgIR0dHhISE4NChQ1pd9+uvv6JSpUpo0aKF7h9qZDqvdRUbG4vs7GwcO3YM7du3x9atW3Hjxg1Mnz4d8+bNM0WMkl27dmHv3r347rvvsGvXLrVje/fuxblz55CWliYlQvPmzUP//v0RHx8PNzc3k8ZGFcvu3FXmDoHIpnENK3XmGLratGkTYmNjsXjxYrRt2xbLli1DVFQUzp07h9q1a2u8LicnB3379sWrr76KGzduGBa0Eejco7Nv3z4sWLAAzz33HOzs7BAQEIB33nkHs2fPRkJCgiliBADcuHEDQ4YMwVdffQVnZ+cSx48ePYqgoCC13p7IyEjk5+cjOTlZY7v5+fnIzc1V24iIiGzd/PnzMWjQIAwePBiNGjVCYmIi/P39sWTJkjKvGzp0KHr16oXQ0NByirRsOic69+7dk56XU61aNfzzzz8AgKZNm+LkyZPGje5fQgj0798fw4YNQ6tWrUo9JzMzE97e3mr7PDw84ODggMzMTI1tJyQkwN3dXdr8/f2NGjsREZHBjDh0VfyP+/z8/BIfV1BQgOTkZERERKjtj4iIwJEjRzSGuWrVKvz111+YMmWKIXdrVDonOs8++ywuXLgAAGjRogWWLVuG9PR0LF26FL6+vjq1FRcXJz2bR9N24sQJLFy4ELm5uZg4cWKZ7ZW22KgQosxFSCdOnIicnBxpS0tL0+keiIiIysPj4St9t8f8/f3V/sAvbTTm5s2bUCgUJToQvL29NXYe/Pnnn/jwww+xbt06VKqk88wYk9Frjk5GRgYAYMqUKYiMjMS6devg4OCA1atX69TW8OHD0bNnzzLPqVOnDqZPn45jx45BLld/qm2rVq3Qu3dvrFmzBj4+Pvjtt9/UjmdnZ6OwsLDEN0qVXC4v0S4REVFFlZaWpjZvtazfgcU7CjR1HigUCvTq1QtTp04t9aHC5qRzotO795PHcgcHB+Pq1av4448/ULt2bXh5eenUlpeXl1bXfPbZZ5g+fbr0/vr164iMjMSmTZvQunVrAEBoaCji4+ORkZEh9Szt3bsXcrkcISEhOsVFRERkUYR4tBnaBgA3N7enFuh4eXnB3t6+RO9NVlZWqZ0HeXl5OHHiBFJSUjB8+HAAgFKphBAClSpVwt69e/HKK68YFr+eDO5bcnZ2RsuWLY0Ri0bFZ3e7uLgAAOrVq4datWoBeDRu2LhxY/Tp0wdz5szB7du3MXbsWAwZMoQVV0REZNXKu+rKwcEBISEhSEpKwhtvvCHtT0pKQpcuXUqc7+bmhtOnT6vtW7x4Mfbt24dvv/0WgYGBesdtKK0SndGjR2vd4Pz58/UOxhD29vbYsWMHYmJi0LZtWzg5OaFXr16YO3euWeIhy7G3YL30uqPbgFLPYTk5ERVn6wt5jh49Gn369EGrVq0QGhqK5cuXIzU1FcOGDQPwaI5reno61q5dCzs7OwQFBaldX6NGDTg6OpbYX960SnRSUlK0aqysSb/GVKdOHYhSuvBq166N7du3l0sMRERE5cYMa1316NEDt27dwrRp05CRkYGgoCDs3LkTAQEBAICMjAykpqYaGJTpyURpGYMNy83Nhbu7O3JycjjkVQGxR4eItKVNj46pf2c8bv+5N6ajUmVHg9oqKnyI41s/srnfb5ZT/0VERESl4+rleuMS5ERERFRhsUeHiIjIwpljrauKgokOVXgRDr2k13aOfDgkEZVOdU4OYGGVVkZ8jo6t4dAVERERVVjs0SEiIrJwHLrSHxMdIiIiS8eqK71x6IqIiIgqLPboEBERWTgOXemPiQ7ZLD4NmYhUWVSVVXGsutIbh66IiIiowmKPDhERkYXj0JX+mOgQERFZOlZd6Y2JDhERkYVjj47+OEeHiIiIKiz26BAREVk6pXi0GdqGDWKiQzZF+TDf3CEQEemOc3T0xqErIiIiqrDYo0NERGThZDDCZGSjRGJ9mOgQERFZOj4ZWW8cuiIiIqIKiz06REREFo7P0dEfEx2q8PYWrDd3CEREhmHVld44dEVEREQVFnt0iIiILJxMCMgMnExs6PXWiokOERGRpVP+uxnahg1iokNERGTh2KOjP87RISIiogqLPTpERESWjlVXemOiQ0REZOn4ZGS9ceiKiIiIKiyrSnR27NiB1q1bw8nJCV5eXnjzzTfVjqempiI6OhpVqlSBl5cXRo4ciYKCAjNFS0REZByPn4xs6GaLrGbo6rvvvsOQIUMwY8YMvPLKKxBC4PTp09JxhUKBTp06oXr16jh8+DBu3bqFfv36QQiBhQsXmjFyIiIiA3HoSm9WkegUFRVh1KhRmDNnDgYNGiTtf/bZZ6XXe/fuxblz55CWlgY/Pz8AwLx589C/f3/Ex8fDzc2t3OMmIiIi87KKoauTJ08iPT0ddnZ2CA4Ohq+vL6KionD27FnpnKNHjyIoKEhKcgAgMjIS+fn5SE5O1th2fn4+cnNz1TYiIiJLIlMaZ7NFVpHoXL58GQAQFxeHjz76CNu3b4eHhwfatWuH27dvAwAyMzPh7e2tdp2HhwccHByQmZmpse2EhAS4u7tLm7+/v+luhIiISB+Ph64M3WyQWROduLg4yGSyMrcTJ05AqXyUhk6aNAlvvfUWQkJCsGrVKshkMmzevFlqTyaTlfgMIUSp+x+bOHEicnJypC0tLc34N0pERERmYdY5OsOHD0fPnj3LPKdOnTrIy8sDADRu3FjaL5fLUbduXaSmpgIAfHx88Ntvv6ldm52djcLCwhI9Parkcjnkcrm+t0BERGR6fGCg3sya6Hh5ecHLy+up54WEhEAul+PChQt48cUXAQCFhYW4evUqAgICAAChoaGIj49HRkYGfH19ATyaoCyXyxESEmK6myAiIjIxrnWlP6uounJzc8OwYcMwZcoU+Pv7IyAgAHPmzAEAdOvWDQAQERGBxo0bo0+fPpgzZw5u376NsWPHYsiQIay4IiIi68bycr1ZRaIDAHPmzEGlSpXQp08fPHjwAK1bt8a+ffvg4eEBALC3t8eOHTsQExODtm3bwsnJCb169cLcuXPNHDkRERGZi0wIG03xNMjNzYW7uztycnLYE2TDwu26Sa+TlJvLOJOIbJmpf2c8br99y4moZO9oUFtFiof45WSCzf1+s4ryciIiIlv2eI6OoZuuFi9ejMDAQDg6OiIkJASHDh3SeO6WLVsQHh6O6tWrw83NDaGhodizZ48ht20UTHSIiIiohE2bNiE2NhaTJk1CSkoKXnrpJURFRUnVzsUdPHgQ4eHh2LlzJ5KTk9G+fXtER0cjJSWlnCNXx6GrYjh0RQCHrohIO+U1dPVKiw9Ryd6wR6EUKfKx79RMrWNt3bo1WrZsiSVLlkj7GjVqhK5duyIhIUGrz2zSpAl69OiByZMn6x23odijQ0REZOmM+GTk4sse5efnl/i4goICJCcnIyIiQm1/REQEjhw5olXISqUSeXl5qFatmuH3bwAmOkRERDbE399fbemj0npnbt68CYVCUeKBu97e3mUuq6Rq3rx5uHfvHrp3726UuPVlNeXlRERENksJQPNqRtq3ASAtLU1t6Kqs1QGKL6H0tGWVHtuwYQPi4uLwww8/oEaNGvrFayRMdMhmlTUPh/NybIPqz4C2TPmzoU88tsZW/9805pOR3dzcnjpHx8vLC/b29iV6b7KysspcVgl4NIl50KBB2Lx5Mzp06GBQzMbAoSsiIiJLV86rlzs4OCAkJARJSUlq+5OSktCmTRuN123YsAH9+/fH+vXr0alTJ71v15jYo0NEREQljB49Gn369EGrVq0QGhqK5cuXIzU1FcOGDQMATJw4Eenp6Vi7di2AR0lO37598emnn+KFF16QeoOcnJzg7u5utvtgokNERGTpzLDWVY8ePXDr1i1MmzYNGRkZCAoKws6dO6XFtDMyMtSeqbNs2TIUFRXh/fffx/vvvy/t79evH1avXm1Y7AZgokNERGTpzLSoZ0xMDGJiYko9Vjx52b9/vx5BmR7n6BAREVGFxR4dsimaqlqK77fVyg59aapg07aKqPjX25KrjyriPVkTY3wdrfL/byOWl9saJjpEREQWzpjl5baGQ1dERERUYbFHh4iIyNKZaTJyRcBEh4iIyNIpBSAzMFFR2maiw6ErIiIiqrDYo0Nqylr/yVRt61tFoakNY8Rtyq+DuZRX1Y8+n1MRK5Iq4j1VFJq+Nxb9/zqHrvTGRIeIiMjiGSHRARMdIiIiskTs0dEb5+gQERFRhcUeHSIiIkunFDB46MlGq66Y6BAREVk6oXy0GdqGDeLQFREREVVY7NGxAuYqUzXl5xqjbW0X6NTEoktJy8CyZSLjsZp/BzgZWW9MdIiIiCwd5+jojUNXREREVGGxR4eIiMjScehKb0x0iIiILJ2AERIdo0RidTh0RURERBWW1fToXLx4EePGjcOvv/6KgoICNG3aFNOnT0f79u2lc1JTU/H+++9j3759cHJyQq9evTB37lw4ODiYMXL9FohkZY3+zFVFYcpqL/48EBmP1VRaqeLQld6spkenU6dOKCoqwr59+5CcnIwWLVrg9ddfR2ZmJgBAoVCgU6dOuHfvHg4fPoyNGzfiu+++w5gxY8wcORERkYGUSuNsNsgqEp2bN2/i0qVL+PDDD9GsWTPUr18fM2fOxP3793H27FkAwN69e3Hu3Dl8/fXXCA4ORocOHTBv3jysWLECubm5Zr4DIiIiAzzu0TF0s0FWkeh4enqiUaNGWLt2Le7du4eioiIsW7YM3t7eCAkJAQAcPXoUQUFB8PPzk66LjIxEfn4+kpOTNbadn5+P3NxctY2IiIgqBquYoyOTyZCUlIQuXbrA1dUVdnZ28Pb2xu7du1G1alUAQGZmJry9vdWu8/DwgIODgzS8VZqEhARMnTrVlOETEREZhnN09GbWHp24uDjIZLIytxMnTkAIgZiYGNSoUQOHDh3C77//ji5duuD1119HRkaG1J5MJivxGUKIUvc/NnHiROTk5EhbWlqaSe6ViIhIb0phnM0GmbVHZ/jw4ejZs2eZ59SpUwf79u3D9u3bkZ2dDTc3NwDA4sWLkZSUhDVr1uDDDz+Ej48PfvvtN7Vrs7OzUVhYWKKnR5VcLodcLjf8ZsrASqvyZeyvY1nfP30+i99nIvPSpxKWrJdZEx0vLy94eXk99bz79+8DAOzs1Dug7OzsoPx3FnloaCji4+ORkZEBX19fAI8mKMvlcmkeDxERkTUSQgkhDKuaMvR6a2UVk5FDQ0Ph4eGBfv364b///a/0TJ0rV66gU6dOAICIiAg0btwYffr0QUpKCn7++WeMHTsWQ4YMkXqBiIiIrJIwwrAV5+hYLi8vL+zevRt3797FK6+8glatWuHw4cP44Ycf0Lx5cwCAvb09duzYAUdHR7Rt2xbdu3dH165dMXfuXDNHT0REROZiFVVXANCqVSvs2bOnzHNq166N7du3l1NERERE5UQIGLxYlY326FhNokNERGSzlEpAZuAcG87RISIiIqpY2KOjQRf3vqgkq2y2z1cteWQ5snnx609UcRn6/3eRKDRSJE/BoSu9MdEhIiKycEKphDBw6MpWy8uZ6BAREVk69ujojXN0iIiIqMJijw4REZGlUwpAxh4dfTDRISIisnRCADC0vJyJDlkQVvoQEREZjokOERGRhRNKAWHg0JVgjw4RERFZJKGE4UNXtllezqorIiIiqrDYo0NERGThOHSlPyY6RERElo5DV3pjolPM44y3CIUGP4SSiIgqtiI8WuvK1L0lxvid9DhWW8NEp5i8vDwAwGHsNHMkRERkLfLy8uDu7m70dh0cHODj44PDmcb5neTj4wMHBwejtGUtZMJWB+00UCqVuH79OlxdXSGTycwSQ25uLvz9/ZGWlgY3NzezxGAqvDfrxHuzThX53gDLuD8hBPLy8uDn5wc7O9PU9zx8+BAFBQVGacvBwQGOjo5GactasEenGDs7O9SqVcvcYQAA3NzcKuQ/TgDvzVrx3qxTRb43wPz3Z4qeHFWOjo42l5wYE8vLiYiIqMJiokNEREQVFhMdCySXyzFlyhTI5XJzh2J0vDfrxHuzThX53oCKf39kHJyMTERERBUWe3SIiIiowmKiQ0RERBUWEx0iIiKqsJjoEBERUYXFRMeM4uPj0aZNGzg7O6Nq1aqlniOTyUpsS5cuVTvn9OnTaNeuHZycnFCzZk1MmzbN7KvUanNvqampiI6ORpUqVeDl5YWRI0eWePqnJd5baerUqVPi+/Thhx+qnaPN/VqqxYsXIzAwEI6OjggJCcGhQ4fMHZJO4uLiSnx/fHx8pONCCMTFxcHPzw9OTk4ICwvD2bNnzRhx2Q4ePIjo6Gj4+flBJpPh+++/Vzuuzf3k5+djxIgR8PLyQpUqVdC5c2f8/fff5XgXpXvavfXv37/E9/KFF15QO8dS743Mg4mOGRUUFKBbt2547733yjxv1apVyMjIkLZ+/fpJx3JzcxEeHg4/Pz8cP34cCxcuxNy5czF//nxTh1+mp92bQqFAp06dcO/ePRw+fBgbN27Ed999hzFjxkjnWOq9aTJt2jS179NHH30kHdPmfi3Vpk2bEBsbi0mTJiElJQUvvfQSoqKikJqaau7QdNKkSRO178/p06elY7Nnz8b8+fOxaNEiHD9+HD4+PggPD5fWvrM09+7dQ/PmzbFo0aJSj2tzP7Gxsdi6dSs2btyIw4cP4+7du3j99dehUCjK6zZK9bR7A4COHTuqfS937lRfB8pS743MRJDZrVq1Sri7u5d6DIDYunWrxmsXL14s3N3dxcOHD6V9CQkJws/PTyiVSiNHqjtN97Zz505hZ2cn0tPTpX0bNmwQcrlc5OTkCCEs/95UBQQEiAULFmg8rs39Wqrnn39eDBs2TG1fw4YNxYcffmimiHQ3ZcoU0bx581KPKZVK4ePjI2bOnCnte/jwoXB3dxdLly4tpwj1V/zfCG3u586dO6Jy5cpi48aN0jnp6enCzs5O7N69u9xif5rS/v3r16+f6NKli8ZrrOXeqPywR8cKDB8+HF5eXnjuueewdOlSKJVK6djRo0fRrl07tQdmRUZG4vr167h69aoZotXO0aNHERQUBD8/P2lfZGQk8vPzkZycLJ1jTfc2a9YseHp6okWLFoiPj1cbltLmfi1RQUEBkpOTERERobY/IiICR44cMVNU+vnzzz/h5+eHwMBA9OzZE5cvXwYAXLlyBZmZmWr3KJfL0a5dO6u7R0C7+0lOTkZhYaHaOX5+fggKCrKKe96/fz9q1KiBBg0aYMiQIcjKypKOWfu9kfFxUU8L98knn+DVV1+Fk5MTfv75Z4wZMwY3b96UhkUyMzNRp04dtWu8vb2lY4GBgeUdslYyMzOlOB/z8PCAg4MDMjMzpXOs5d5GjRqFli1bwsPDA7///jsmTpyIK1eu4IsvvgCg3f1aops3b0KhUJSI3dvb26LjLq5169ZYu3YtGjRogBs3bmD69Olo06YNzp49K91Hafd47do1c4RrEG3uJzMzEw4ODvDw8ChxjqV/X6OiotCtWzcEBATgypUr+Pjjj/HKK68gOTkZcrncqu+NTIM9OkZW2qTH4tuJEye0bu+jjz5CaGgoWrRogTFjxmDatGmYM2eO2jkymUztvfh3sm7x/YYy9r2VFp8QQm1/ed1baXS53w8++ADt2rVDs2bNMHjwYCxduhQrV67ErVu3NN7L4/spj3sxVGnfB2uI+7GoqCi89dZbaNq0KTp06IAdO3YAANasWSOdY+33WJw+92MN99yjRw906tQJQUFBiI6Oxq5du3Dx4kXpe6qJNdwbmQZ7dIxs+PDh6NmzZ5nnFO+l0MULL7yA3Nxc3LhxA97e3vDx8SnxV8rjbtzif9EZypj35uPjg99++01tX3Z2NgoLC6W4y/PeSmPI/T6uArl06RI8PT21ul9L5OXlBXt7+1K/D5Yc99NUqVIFTZs2xZ9//omuXbsCeNTL4evrK51jrff4uJqsrPvx8fFBQUEBsrOz1Xo+srKy0KZNm/IN2EC+vr4ICAjAn3/+CaBi3RsZB3t0jMzLywsNGzYsc3N0dNS7/ZSUFDg6Okol26GhoTh48KDafJC9e/fCz8/PoISqNMa8t9DQUJw5cwYZGRlqccvlcoSEhJT7vZXGkPtNSUkBAOkXjTb3a4kcHBwQEhKCpKQktf1JSUlW/UsjPz8f58+fh6+vLwIDA+Hj46N2jwUFBThw4IBV3qM29xMSEoLKlSurnZORkYEzZ85Y3T3funULaWlp0v9rFeneyEjMNg2axLVr10RKSoqYOnWqcHFxESkpKSIlJUXk5eUJIYTYtm2bWL58uTh9+rS4dOmSWLFihXBzcxMjR46U2rhz547w9vYWb7/9tjh9+rTYsmWLcHNzE3PnzjXXbQkhnn5vRUVFIigoSLz66qvi5MmT4qeffhK1atUSw4cPl9qw1Hsr7siRI2L+/PkiJSVFXL58WWzatEn4+fmJzp07S+doc7+WauPGjaJy5cpi5cqV4ty5cyI2NlZUqVJFXL161dyhaW3MmDFi//794vLly+LYsWPi9ddfF66urtI9zJw5U7i7u4stW7aI06dPi7ffflv4+vqK3NxcM0deury8POn/KQDSz9+1a9eEENrdz7Bhw0StWrXETz/9JE6ePCleeeUV0bx5c1FUVGSu2xJClH1veXl5YsyYMeLIkSPiypUr4pdffhGhoaGiZs2aVnFvZB5MdMyoX79+AkCJ7ZdffhFCCLFr1y7RokUL4eLiIpydnUVQUJBITEwUhYWFau3873//Ey+99JKQy+XCx8dHxMXFmb38+mn3JsSjZKhTp07CyclJVKtWTQwfPlytlFwIy7y34pKTk0Xr1q2Fu7u7cHR0FM8++6yYMmWKuHfvntp52tyvpfr8889FQECAcHBwEC1bthQHDhwwd0g66dGjh/D19RWVK1cWfn5+4s033xRnz56VjiuVSjFlyhTh4+Mj5HK5ePnll8Xp06fNGHHZfvnll1L//+rXr58QQrv7efDggRg+fLioVq2acHJyEq+//rpITU01w92oK+ve7t+/LyIiIkT16tVF5cqVRe3atUW/fv1KxG2p90bmIRPCAh8zS0RERGQEnKNDREREFRYTHSIiIqqwmOgQERFRhcVEh4iIiCosJjpERERUYTHRISIiogqLiQ4RERFVWEx0qMIJCwtDbGxshfrc/v37S2sy6atOnTrSYqR37tzReN7q1aulJUbI+Pr37y99H77//ntzh0NU4THRITKSLVu24JNPPpHe16lTB4mJieYLqBTTpk1DRkYG3N3dzR1Khbd///5Sk8pPP/1Ubc0zIjItrl5OZCTVqlUzdwhP5erqKq1ubW6FhYWoXLmyucMod+7u7kw0icoRe3SowsvOzkbfvn3h4eEBZ2dnREVF4c8//5SOPx6q2bNnDxo1agQXFxd07NhR7a/uoqIijBw5ElWrVoWnpycmTJiAfv36qQ0nqQ5dhYWF4dq1a/jggw+kYQoAiIuLQ4sWLdTiS0xMVFuNXaFQYPTo0dJnjR8/HsVXahFCYPbs2ahbty6cnJzQvHlzfPvtt3p9fVavXo3atWvD2dkZb7zxBm7dulXinB9//BEhISFwdHRE3bp1MXXqVBQVFUnH//jjD7z44otwdHRE48aN8dNPP6kNzVy9ehUymQzffPMNwsLC4OjoiK+//hoAsGrVKjRq1AiOjo5o2LAhFi9erPbZ6enp6NGjBzw8PODp6YkuXbrg6tWr0vH9+/fj+eefR5UqVVC1alW0bdsW165d0+ren3Zf8+fPR9OmTVGlShX4+/sjJiYGd+/elY5fu3YN0dHR8PDwQJUqVdCkSRPs3LkTV69eRfv27QEAHh4ekMlk6N+/v1YxEZFxMdGhCq9///44ceIEtm3bhqNHj0IIgddeew2FhYXSOffv38fcuXPx1Vdf4eDBg0hNTcXYsWOl47NmzcK6deuwatUq/Prrr8jNzS1zfsWWLVtQq1YtaahIl6GKefPm4csvv8TKlStx+PBh3L59G1u3blU756OPPsKqVauwZMkSnD17Fh988AHeeecdHDhwQPsvDIDffvsNAwcORExMDE6dOoX27dtj+vTpaufs2bMH77zzDkaOHIlz585h2bJlWL16NeLj4wEASqUSXbt2hbOzM3777TcsX74ckyZNKvXzJkyYgJEjR+L8+fOIjIzEihUrMGnSJMTHx+P8+fOYMWMGPv74Y6xZswbAo+9L+/bt4eLigoMHD+Lw4cNSIlpQUICioiJ07doV7dq1w//+9z8cPXoU7777rpRYluVp9wUAdnZ2+Oyzz3DmzBmsWbMG+/btw/jx46Xj77//PvLz83Hw4EGcPn0as2bNgouLC/z9/fHdd98BAC5cuICMjAx8+umnOn1viMhIzLqkKJEJtGvXTowaNUoIIcTFixcFAPHrr79Kx2/evCmcnJzEN998I4QQYtWqVQKAuHTpknTO559/Lry9vaX33t7eYs6cOdL7oqIiUbt2bdGlS5dSP1cIIQICAsSCBQvUYpsyZYpo3ry52r4FCxaIgIAA6b2vr6+YOXOm9L6wsFDUqlVL+qy7d+8KR0dHceTIEbV2Bg0aJN5++22NX5fS4nn77bdFx44d1fb16NFDuLu7S+9feuklMWPGDLVzvvrqK+Hr6yuEEGLXrl2iUqVKIiMjQzqelJQkAIitW7cKIYS4cuWKACASExPV2vH39xfr169X2/fJJ5+I0NBQIYQQK1euFM8++6zaivX5+fnCyclJ7NmzR9y6dUsAEPv379d435o87b5K88033whPT0/pfdOmTUVcXFyp5z5ehTs7O7vU46pfHyIyHc7RoQrt/PnzqFSpElq3bi3t8/T0xLPPPovz589L+5ydnVGvXj3pva+vL7KysgAAOTk5uHHjBp5//nnpuL29PUJCQqBUKo0ab05ODjIyMhAaGirtq1SpElq1aiUNX507dw4PHz5EeHi42rUFBQUIDg7W6fPOnz+PN954Q21faGgodu/eLb1PTk7G8ePH1Xo6FAoFHj58iPv37+PChQvw9/dXm/uj+rVS1apVK+n1P//8g7S0NAwaNAhDhgyR9hcVFUlzWJKTk3Hp0iW4urqqtfPw4UP89ddfiIiIQP/+/REZGYnw8HB06NAB3bt3h6+v71Pv/Wn35ezsjF9++QUzZszAuXPnkJubi6KiIjx8+BD37t1DlSpVMHLkSLz33nvYu3cvOnTogLfeegvNmjV76mcTUflhokMVmig2t0V1v+rwRvFJsTKZrMS1xYdDNLVdFjs7uxLXqQ6haeNxcrVjxw7UrFlT7ZhcLtepLW3uQalUYurUqXjzzTdLHHN0dCzxtSxLlSpV1NoFgBUrVqglosCjRPLxOSEhIVi3bl2JtqpXrw7g0RyfkSNHYvfu3di0aRM++ugjJCUl4YUXXjDovq5du4bXXnsNw4YNwyeffIJq1arh8OHDGDRokPQ9Gzx4MCIjI7Fjxw7s3bsXCQkJmDdvHkaMGKHV14OITI+JDlVojRs3RlFREX777Te0adMGAHDr1i1cvHgRjRo10qoNd3d3eHt74/fff8dLL70E4NFf/ikpKSUmFqtycHCAQqFQ21e9enVkZmaqJQenTp1S+yxfX18cO3YML7/8MoBHPRzJyclo2bKldE9yuRypqalo166dVvegSePGjXHs2DG1fcXft2zZEhcuXMAzzzxTahsNGzZEamoqbty4AW9vbwDA8ePHn/rZ3t7eqFmzJi5fvozevXuXek7Lli2xadMm1KhRA25ubhrbCg4ORnBwMCZOnIjQ0FCsX7/+qYnO0+7rxIkTKCoqwrx582Bn92g64zfffFPiPH9/fwwbNgzDhg3DxIkTsWLFCowYMQIODg4AUOJngIjKFxMdqtDq16+PLl26YMiQIVi2bBlcXV3x4YcfombNmujSpYvW7YwYMQIJCQl45pln0LBhQyxcuBDZ2dll9mTUqVMHBw8eRM+ePSGXy+Hl5YWwsDD8888/mD17Nv7zn/9g9+7d2LVrl9ov8VGjRmHmzJmoX78+GjVqhPnz56s9i8XV1RVjx47FBx98AKVSiRdffBG5ubk4cuQIXFxc0K9fP63va+TIkWjTpg1mz56Nrl27Yu/evWrDVgAwefJkvP766/D390e3bt1gZ2eH//3vfzh9+jSmT5+O8PBw1KtXD/369cPs2bORl5cnTUZ+Wk9PXFwcRo4cCTc3N0RFRSE/Px8nTpxAdnY2Ro8ejd69e2POnDno0qULpk2bhlq1aiE1NRVbtmzBuHHjUFhYiOXLl6Nz587w8/PDhQsXcPHiRfTt2/ep9/60+6pXrx6KioqwcOFCREdH49dff8XSpUvV2oiNjUVUVBQaNGiA7Oxs7Nu3T0qgAwICIJPJsH37drz22mtwcnKCi4uL1t8bIjISs80OIjKR4pOCb9++Lfr06SPc3d2Fk5OTiIyMFBcvXpSOr1q1Sm3yrRBCbN26Vaj+71FYWCiGDx8u3NzchIeHh5gwYYLo1q2b6Nmzp8bPPXr0qGjWrJmQy+VqbS1ZskT4+/uLKlWqiL59+4r4+Hi1yciFhYVi1KhRws3NTVStWlWMHj1a9O3bV23is1KpFJ9++ql49tlnReXKlUX16tVFZGSkOHDggMavS2mTkYV4NOG3Vq1awsnJSURHR4u5c+eW+Hrs3r1btGnTRjg5OQk3Nzfx/PPPi+XLl0vHz58/L9q2bSscHBxEw4YNxY8//igAiN27dwshnkxGTklJKfH569atEy1atBAODg7Cw8NDvPzyy2LLli3S8YyMDNG3b1/h5eUl5HK5qFu3rhgyZIjIyckRmZmZomvXrsLX11c4ODiIgIAAMXnyZKFQKDR+HXS5r/nz5wtfX1/p52bt2rVqE4yHDx8u6tWrJ+Ryuahevbro06ePuHnzpnT9tGnThI+Pj5DJZKJfv35qnw1ORiYqFzIh9JhoQGTjlEolGjVqhO7du6s9DdmS1alTB7GxseWyPMavv/6KF198EZcuXVKb5E1PyGQybN261eClPYiobHyODpEWrl27hhUrVuDixYs4ffo03nvvPVy5cgW9evUyd2g6mTBhAlxcXJCTk2PUdrdu3YqkpCRcvXoVP/30E9599120bduWSU4phg0bxiEsonLEHh0iLaSlpaFnz544c+YMhBAICgrCzJkzpQnD1uDatWtStVDdunWlCbbGsHbtWnzyySdIS0uDl5cXOnTogHnz5sHT09Non6GrJk2aaHxC8rJlyzROgDa1rKws5ObmAnj0GAPVSjQiMj4mOkRUIakmdsV5e3uXeDYPEVVMTHSIiIiowuIcHSIiIqqwmOgQERFRhcVEh4iIiCosJjpERERUYTHRISIiogqLiQ4RERFVWEx0iIiIqMJiokNEREQV1v8DfbDKusPqH5IAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "gdf = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))\n", "geo_grid = make_geocube(\n", " vector_data=gdf,\n", " measurements=['pop_est'],\n", " resolution=(-3.0, 3.0),\n", ")\n", "geo_grid[\"pop_est\"].plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Convert to GeoDataFrame (vectorize)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAisAAAEbCAYAAAAbPMQYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAApDklEQVR4nO3df3CUVZ7v8U8HkiYgCT8C6UQggDPuOMaREl0ItYqohGQRf01xQbcUdgTXhagI7lwpmCWgiA7Ieq8OMiiV0ju6UHNXXadE+aGAawGKmL0L6jgw8iNKIsJAGhA6kD73D5c2z5Okf6T76X46/X5VdVWe5zl9nnO6Q+fLOd9z2mOMMQIAAHCprFQ3AAAAIByCFQAA4GoEKwAAwNUIVgAAgKsRrAAAAFcjWAEAAK5GsAIAAFyta6obEK9gMKjDhw+rZ8+e8ng8qW4OAACIgjFGJ0+eVHFxsbKywo+dpH2wcvjwYQ0cODDVzQAAAB1QV1enAQMGhC2T9sFKz549JX3f2by8vBS3BgAARMPv92vgwIGhv+PhpH2wcmHqJy8vj2AFAIA0E00KBwm2AADA1QhWAACAqxGsAAAAVyNYAQAArkawAgAAXI1gBQAAuBrBCgAAcDWCFQAA4GppvykcEItxPe5JyX3Xn345YXWNzZqYsLrisTH4+1Q3AUCGYGQFAAC4GsEKAABwNYIVAADgauSsIO0lKw8llrwTe16JW/JMEsneJzfmsER63Vu2OR3601lUXvo/27329p+eSmJLkC4YWQEAAK5GsAIAAFyNYAUAALiaxxhjnKp88ODBOnjwYKvzM2bM0G9+8xtNnTpVL730kuXaiBEjtGPHjqjv4ff7lZ+fr8bGRuXl5cXdZkQvUq5IuBwP+3PjKetUzkqse6N0xryUcMjpSD/2XJF0yA8Jl99ilw79wQ9i+fvtaILtzp071dzcHDres2ePxo4dq4kTf/hQr6ioUE1NTeg4JyfHySYBAIA042iw0q9fP8vxk08+qUsuuUSjR48OnfN6vfL5fFHXGQgEFAgEQsd+vz/+hgIAANdK2tLlpqYm/e53v9Ps2bPl8XhC57ds2aL+/furV69eGj16tBYvXqz+/fu3W8+SJUu0cOHCZDQZDoo0lRPLtFAi7xuPeKZFMm0KCR3Xclok1mmPRE6TxDI9kyzxTBk5OUUWz3uG7yUtwfaNN97QiRMnNHXq1NC5yspKvfLKK3rvvff09NNPa+fOnbrhhhssIyd2c+fOVWNjY+hRV1eXhNYDAIBUSdrIyurVq1VZWani4uLQuUmTJoV+Li0t1dVXX62SkhK99dZbuuOOO9qsx+v1yuv1Ot5eAADgDkkJVg4ePKhNmzbptddeC1uuqKhIJSUl2rt3bzKaBQAA0oCjS5cvqK6u1m9/+1vV1dWpa9f246Njx47p4osv1qpVq3TPPdHlEqTr0uXrbl3a4ee+/+//lMCWJEci80wi1d2yfLK24rffNx5uyV/J9KXJsWzVny5i2ebejTkpsEr3/JdY/n47nrMSDAZVU1OjKVOmWAKVU6dO6ZFHHtH27dt14MABbdmyRRMmTFBBQYFuv/12p5sFAADShOPTQJs2bdKhQ4f0i1/8wnK+S5cu2r17t15++WWdOHFCRUVFGjNmjNauXauePXs63SwAAJAmHA9WysvL1dZMU25urtavX+/07QEAQJpLSs6Kk5zOWUnU+vh4clTikbtpd0rum67C5bvEs0dLLPksycxZSce8i3iM7TIpcqEodb1kcIefG+6zJFW5IuSsIBwn8mNclbMCAAAQD4IVAADgagQrAADA1chZSZFk5bCQsxKbRO3Rkqg9V6TYcliSmYMSrl1uyYVp1UZP+/8/29i81vrcCPkt8eSsAG6XjD1cyFkBAACdBsEKAABwtaR9kWGmS9XS5UyUyCkYN9zHPqXCdvxWYV+PMNM+knXqJ5HLmgEnRFpeHs/yc7dv3c/ICgAAcDWCFQAA4GoEKwAAwNXIWUmS9//9n2Iq3zLHJdbnRluvxNLmeITbjj/SVv3xcEuuSCJVFNwX+vmdo6s6XE+r/J4Y8lBaLV225cJ0/dHQDrcLmS1V+SBuz0OJBSMrAADA1QhWAACAqxGsAAAAVyNnxaUSmacSS72dIacl3Db5Tu6NEu6+TuawZLpYcnjseSktRdq/5vyfD1iO2W4fFyQrNySWfVM6G0ZWAACAqxGsAAAAVyNYAQAAruYxxphUNyIesXzFNGIXLg+jMwq3V0qy7psJ7PkhLfNOWu650pag/5TluEufXlHfN3jyVORCF8qeDURd1i35K80H6hyru8vggY7VnQpvvf962Ovjr7vdkfvGk98S6buA0k0sf78dHVmprq6Wx+OxPHw+X+i6MUbV1dUqLi5Wbm6urr/+en366adONgkAAKQZx6eBLr/8ctXX14ceu3f/sLrk17/+tZYvX67nnntOO3fulM/n09ixY3Xy5EmnmwUAANKE48FK165d5fP5Qo9+/fpJ+n5U5ZlnntG8efN0xx13qLS0VC+99JK+++47vfrqq043CwAApAnH91nZu3eviouL5fV6NWLECD3xxBMaOnSo9u/fr4aGBpWXl4fKer1ejR49Wtu2bdM//MM/tFlfIBBQIPDDXLLf73e6CxktXC5FZ8xnSVSfMi0HJVbxfL9RVt5FUZc9/+2x8HV187Z7rUuP7lHfxzQcCXvd4+sfdV1O5p3EI5Z22fNbzu/7MtHNkRT5+5rs911/+P9FXbc9p8WpHJZYpHuOSjwcHVkZMWKEXn75Za1fv14vvPCCGhoaNGrUKB07dkwNDQ2SpMLCQstzCgsLQ9fasmTJEuXn54ceAwd2rqQvAABg5WiwUllZqZ///Oe64oordNNNN+mtt96SJL300kuhMh6Px/IcY0yrcy3NnTtXjY2NoUddnTv/FwIAABIjqdvt9+jRQ1dccYX27t2r2267TZLU0NCgoqKiUJkjR460Gm1pyev1yuttf+gWyXPmpitS3YS4pePXCcCq5VRPl3zr8kf7cfOJE5bjlsuTPV26WK6Zc+etx+fPWY4ty63z/t5y7R1/jeU4nm3SNzRZc/jKc+7qcF1OeefQx7Yz1uNxxVc6cl+nppfaEm6ps32KyD5dk8nb5CdKUjeFCwQC+vzzz1VUVKQhQ4bI5/Np48aNoetNTU3aunWrRo0alcxmAQAAF3N0ZOWRRx7RhAkTNGjQIB05ckSPP/64/H6/pkyZIo/Ho1mzZumJJ57Qj3/8Y/34xz/WE088oe7du+uuu9z3PwcAAJAajgYrX331le68804dPXpU/fr108iRI7Vjxw6VlJRIkn75y1/qzJkzmjFjho4fP64RI0Zow4YN6tmzp5PNAgAAacTRYGXNmjVhr3s8HlVXV6u6utrJZgDtsufdtMxhiTUn5/1//6eEtCkTRMrxiMXG5rU/1NtneofrMc3N4e8TZrl1PO23iyVHxZMd20e4PQ8nHq3zVNoXbslwIvNZYlmaHI9IOSqtnA//u4XI+CJDAADgagQrAADA1QhWAACAq3mMMSbVjYhHLF8xjcS67talqW5CUpGTkjhju0yyHIfb2j6R+SD2nBbLvise6//dWubCtKXSNyMhbYr0lQD2dthfu3DCfZ2AXcs9Z9qy/uvaqOtKB1m+P1mOgw2XWo4rB19jOV5/9hXH25RpYvn7zcgKAABwNYIVAADgagQrAADA1ZL63UAAkEz2HJV3/vJCu2Uj5YIkKkfFrmu/vgkrb89/seehxJLD0tlyVOzsOSr2HBaP9xcdrrty6COW47e/XNbhuvA9RlYAAICrEawAAABXY+kyEqazLWVmqXLyjOtxzw8HHk+H61l/6qUEtCZ+4bbIlyRPTnbo56yLeiTsvsGTpzr8XHP5jyzH//eNFy3HF2V163DdbmSf9omFfZrHjmmf6LB0GQAAdBoEKwAAwNUIVgAAgKuxdBkJY8/xCDeve/qKQstxj93fhL2eLOSpJIclRyWBIi0/DreFvv25XYusv4Nvf/W/o26HJzv6j9bgqdOW40TmsERiz1PJJJGWLtu1/DwjJyX5GFkBAACuRrACAABcjWAFAAC4GvusICUi7VMQjj2fJZ48E/veMOSsOKNVjkoce6mEEzxzNmF12XNWwrHns8STk7P+9MuW43i2+W8e5Iu6rH1flUgmXjI69PPb+z+M6bmngj+8T27ZvyWefVfQMeyzAgAAOg1Hg5UlS5bommuuUc+ePdW/f3/ddttt+uKLLyxlpk6dKo/HY3mMHDnSyWYBAIA04miwsnXrVs2cOVM7duzQxo0bdf78eZWXl+v0aetSvYqKCtXX14ce69atc7JZAAAgjSQ1Z+Xbb79V//79tXXrVl133XWSvh9ZOXHihN54440O1UnOSucQTw4Lex64jz1nI3jmTNjyWd27J+S+qcpZadWO4ycS1g57Dks49vwWJ3NWEuV/XHaT5XjdFx8k5b7kqKSea3NWGhsbJUl9+vSxnN+yZYv69++vSy+9VNOnT9eRI0farSMQCMjv91seAACg80pasGKM0ezZs/U3f/M3Ki0tDZ2vrKzUK6+8ovfee09PP/20du7cqRtuuEGBQKDNepYsWaL8/PzQY+DAgcnqAgAASIGkTQPNnDlTb731lj744AMNGDCg3XL19fUqKSnRmjVrdMcdd7S6HggELIGM3+/XwIEDmQbqZMJNCzk57VOePTn084Zzaxy7T2c3NmtiTOU9XbPbv5bT/rVWzp8Pe9k0N1uPg+1//MUzDZTIpcx24abUNgZ/H1Nd/sOD4m1OVPKKD1mOK/J/0W5Z+zTQ+OHj2i371q71Ye/LVI+7xTINlJTvBnrggQf05ptv6v333w8bqEhSUVGRSkpKtHfv3jave71eeb1eJ5oJAABcyNFgxRijBx54QK+//rq2bNmiIUOGRHzOsWPHVFdXp6KiIiebBgAA0oSjOSszZ87U7373O7366qvq2bOnGhoa1NDQoDP/PYx56tQpPfLII9q+fbsOHDigLVu2aMKECSooKNDtt9/uZNMAAECacDRnxdPOlto1NTWaOnWqzpw5o9tuu021tbU6ceKEioqKNGbMGD322GNRJ86ydLlzimUpc7gcFns9icx3GdtlkuV4Y/PahNWd7hKZs9KqrD2HJUKeSkv2nJWw94ljunn9qZfCXo8nh6VlzkqsOSrhxJO/Ys9JiSRczoqyrQP+zccbLcct/50FGy61XCNHJb24JmclUhyUm5ur9evDJ0gBAIDMxncDAQAAVyNYAQAArpaUpctArGLJQwmX33L+wMGEtcmOHBWrWPNUOiyGHJVUiZjPFGZ/l0gSmacSj1jzVDqqS+98y3HL13Zjc8dzVCoHPmQ5frvuf3W4LjiPkRUAAOBqBCsAAMDVCFYAAICrJe27gZzCPisIx547YEfeSeK0zFnxdOlivejp+P+LPFm2/Zou/5H1+NN9UdcVbp+VePZVscvq0zvs9eC3x9q9tv7M/0lYO+IRad8Vp3JWKgrusxwH/afaLbuh6VVH2oDkiOXvNyMrAADA1QhWAACAq7F0GZ2afZonnuW19i3CTdO5dsvah/Lt93XL8tNEatmn8uzJ1osmaD0OMy1kn0Iy562vc6sv8bBPC7UUwxRR8LvvLMdZ3btH/dxI0z6t7hU42+61TP8ah3eOrgp7vfLiB5LUErgJIysAAMDVCFYAAICrEawAAABXI2cFGSWWXJGwX2Mfwbjcuy3HWd5uHa4rHW04t8ZybM/ZycqxLW0Ow9i2pm+VsxKDcEuXW5UNBOK4U3hdB1zc7rXzX9c7dt9YJGs7/Vi9/fWzjtTL9vvuxsgKAABwNYIVAADgagQrAADA1chZATrIk5NtOQ6374pduP1eutj27Ii070Q6sOcKjev2d5bj9Wdf6XDd44Yv+OEghn1VJKlr/36hn88f+dZyrUvfPpbj5mN/abce0+i3HHvyw28dHi4vJdP2VbnxhiWW43ffm5uwuu25Yy3Z90IiR8XdGFkBAACuRrACAABcjWAFAAC4mityVlasWKGlS5eqvr5el19+uZ555hlde+21qW4WEJvSFt9Rs2df+9ckadcey2HW8NLQz9ZdRTqneHJUWtW1a2Ho53i++6ll/kpb7Dks5rszLQ6s75o9hyWcTMtRsbPnqNhzWLpu/8xy3DLXJFLZ82U/DXsvpI+Uj6ysXbtWs2bN0rx581RbW6trr71WlZWVOnTInRsSAQCA5Ep5sLJ8+XLde++9mjZtmi677DI988wzGjhwoJ5//vk2ywcCAfn9fssDAAB0XimdBmpqatKuXbv06KOPWs6Xl5dr27ZtbT5nyZIlWrhwYZvXgHiV59zV7rWs3Bi2zLdP+9jrajHtE8m4a8L/vq/fuSDsdTikS5ivDAgGwz4106d+wok0LdTRsnbxLJl2crk12pbSkZWjR4+qublZhYWFlvOFhYVqaGho8zlz585VY2Nj6FFXV5eMpgIAgBRxRYKtx2P9ajJjTKtzF3i9Xnm93mQ0CwAAuEBKR1YKCgrUpUuXVqMoR44caTXaAgAAMlNKR1ZycnI0fPhwbdy4Ubfffnvo/MaNG3XrrbemsGWA5LHlJLTaTj9CXkqyhMtpybR8Fvu2/valzJGWJydMlvX/ge8cfzE5981wseSOxJNnEmuuTNbW/wz9TL5Sx6R8Gmj27Nm6++67dfXVV6usrEyrVq3SoUOHdP/996e6aQAAwAVSHqxMmjRJx44d06JFi1RfX6/S0lKtW7dOJSUlqW4aAABwgZQHK5I0Y8YMzZgxI9XNAAAALuSKYAVwC3ueioVLclRiYc9nybQcliyvdW+cYOPJ9svm93S6OeiAcLklkXJFkrX/SSw5LGO7TLIck8MSnZTvYAsAABAOwQoAAHA1ghUAAOBq5KwALbT8+vlxuXdbL+7ZZz225bCY//w89LNn2GUJbxvcJXj6u9DPWT26p7Almcut38nj1nalM0ZWAACAqxGsAAAAVyNYAQAAruYxxphUNyIefr9f+fn5amxsVF5eXqqbg06sVQ6LjTl/Luz1lpKV05Jp+6rYRXrPWjLNzQm7b9j9emTNjbIb1+3vrGXPvpKQNiE+9velpUjvUcvn8n7+IJa/34ysAAAAVyNYAQAArsbSZSBK9qH7WKYY7Fouc5ZY6pwWTLDdS56u2WGfGsu0D1Ijnvchlqk7pvk6hpEVAADgagQrAADA1QhWAACAq5GzAnSQPQ+hPHtyilqC9gQDZxNWV6TlyC2Fy1GJ+FxyGJLCyVwh8pASj5EVAADgagQrAADA1QhWAACAq5GzAiSBCYb/VgtPgu6T6dvr220M/t5yPDZrYofrsm/Hn+Xt1uG6yGlIjnR4ndl3JTqOjawcOHBA9957r4YMGaLc3FxdcsklWrBggZqamizlPB5Pq8fKlSudahYAAEgzjo2s/PGPf1QwGNRvf/tb/ehHP9KePXs0ffp0nT59WsuWLbOUrampUUVFReg4Pz/fqWYBAIA041iwUlFRYQlAhg4dqi+++ELPP/98q2ClV69e8vl8TjUFAACksaTmrDQ2NqpPnz6tzldVVWnatGkaMmSI7r33Xt13333Kymp7hioQCCgQCISO/X6/Y+0FwinPuctynJWb2+G6mj/57Id6rvpp2LLkpaRGlzAjvu8cfzGmushLQHvIYWlb0oKVP//5z3r22Wf19NNPW84/9thjuvHGG5Wbm6t3331Xc+bM0dGjRzV//vw261myZIkWLlyYjCYDAAAXiDnBtrq6us2k2JaPjz/+2PKcw4cPq6KiQhMnTtS0adMs1+bPn6+ysjINGzZMc+bM0aJFi7R06dJ27z937lw1NjaGHnV1dbF2AQAApJGYR1aqqqo0eXL4bcUHDx4c+vnw4cMaM2aMysrKtGrVqoj1jxw5Un6/X998840KCwtbXfd6vfJ6vbE2G3C1SFM/SAz7VE5zY2O71wC4R8zBSkFBgQoKCqIq+/XXX2vMmDEaPny4ampq2s1Daam2tlbdunVTr169Ym0aAADohBzLWTl8+LCuv/56DRo0SMuWLdO3334bunZh5c8f/vAHNTQ0qKysTLm5udq8ebPmzZun++67j9ETAAAgycFgZcOGDdq3b5/27dunAQMGWK4Z8/1untnZ2VqxYoVmz56tYDCooUOHatGiRZo5c6ZTzQIAAGnGYy5EDmnK7/crPz9fjY2NysvLS3VzkEFaLV3u1vHRwOZTp36oZ3ip5RpLlZ1T0Xta5ELtiHW5MlIvHbbft+vMS5dj+fvNFxkCAABXI1gBAACuRrACAABcLanb7QOIjByV1CAHJbM4mQsST25MZ85RiQcjKwAAwNUIVgAAgKsRrAAAAFcjZwWIUiL3VbHbGPx9wupC9MhTySzJygex3ycd9ncZmzXRcuy2zyRGVgAAgKsRrAAAAFdjGghIgnf8NaluAoAUYTly/BhZAQAArkawAgAAXI1gBQAAuBo5K0CUNjS9ajmuyPv7dsuSowLAzdy+VNmOkRUAAOBqBCsAAMDVCFYAAICrkbMCRMnJ7fYBwGkt81TcnqNix8gKAABwNUeDlcGDB8vj8Vgejz76qKXMoUOHNGHCBPXo0UMFBQV68MEH1dTU5GSzAABAGnF8GmjRokWaPn166Piiiy4K/dzc3Kzx48erX79++uCDD3Ts2DFNmTJFxhg9++yzTjcNAACkAceDlZ49e8rn87V5bcOGDfrss89UV1en4uJiSdLTTz+tqVOnavHixcrLy3O6eUDCsLcKADdLtzyVlhzPWXnqqafUt29fDRs2TIsXL7ZM8Wzfvl2lpaWhQEWSxo0bp0AgoF27drVZXyAQkN/vtzwAAEDn5ejIykMPPaSrrrpKvXv31kcffaS5c+dq//79evHFFyVJDQ0NKiwstDynd+/eysnJUUNDQ5t1LlmyRAsXLnSy2QAAwEViHlmprq5ulTRrf3z88ceSpIcfflijR4/Wz372M02bNk0rV67U6tWrdezYsVB9Ho+n1T2MMW2el6S5c+eqsbEx9Kirq4u1C0BCBM8GLA8AgDNiHlmpqqrS5MmTw5YZPHhwm+dHjhwpSdq3b5/69u0rn8+nDz/80FLm+PHjOnfuXKsRlwu8Xq+8Xva3AAAgU8QcrBQUFKigoKBDN6utrZUkFRUVSZLKysq0ePFi1dfXh85t2LBBXq9Xw4cP79A9AABA5+JYzsr27du1Y8cOjRkzRvn5+dq5c6cefvhh3XLLLRo0aJAkqby8XD/96U919913a+nSpfrLX/6iRx55RNOnT2clEAAAkCR5jDHGiYo/+eQTzZgxQ3/84x8VCARUUlKiyZMn65e//KW6d+8eKnfo0CHNmDFD7733nnJzc3XXXXdp2bJlUU/1+P1+5efnq7GxkQAHAIA0Ecvfb8eClWQhWAEAIP3E8veb7wYCAACuRrACAABcjWAFAAC4GsEKAABwNYIVAADgagQrAADA1QhWAACAqxGsAAAAVyNYAQAArkawAgAAXM2xLzIE8IOxWRMtxxuDv09RSwAg/TCyAgAAXI1gBQAAuBrTQECCtJzqsU/zMO2DaNmnDGORzN+zeNqJH/DZEB1GVgAAgKsRrAAAAFcjWAEAAK7mMcaYVDciHn6/X/n5+WpsbFReXl6qm4MMEsucPfPSzgmXKxTPe5QJORkt+5wJ/e3s0u1zJpa/34ysAAAAVyNYAQAArkawAgAAXM2xnJUtW7ZozJgxbV776KOPdM0113zfAI+n1fXnn39e999/f1T3IWcldZK1hXyk+yRyX4pw+Q+R2hXPfTsb8h8A93Hb504sf78d2xRu1KhRqq+vt5z71a9+pU2bNunqq6+2nK+pqVFFRUXoOD8/36lmAQCANONYsJKTkyOfzxc6PnfunN58801VVVW1Gk3p1auXpSwAAMAFSctZefPNN3X06FFNnTq11bWqqioVFBTommuu0cqVKxUMBtutJxAIyO/3Wx4AAKDzSto+K3/7t38rSVq3bp3l/OOPP64bb7xRubm5evfdd/XP//zPmjt3rubPn99mPdXV1Vq4cGGr85mcs0J+gDu4bT44XvxeAenN7Z9Jju6zUl1dLY/HE/bx8ccfW57z1Vdfaf369br33ntb1Td//nyVlZVp2LBhmjNnjhYtWqSlS5e2e/+5c+eqsbEx9Kirq4u1CwAAII3EnLNSVVWlyZMnhy0zePBgy3FNTY369u2rW265JWL9I0eOlN/v1zfffKPCwsJW171er7xeb0xtBgAA6SvmYKWgoEAFBQVRlzfGqKamRvfcc4+ys7Mjlq+trVW3bt3Uq1evWJuWVuJZ9svwfHK4ZQg1mdv687sFpC+3fGY5wbHVQBe899572r9/f5tTQH/4wx/U0NCgsrIy5ebmavPmzZo3b57uu+8+Rk8AAICkJAQrq1ev1qhRo3TZZZe1upadna0VK1Zo9uzZCgaDGjp0qBYtWqSZM2c63SwAAJAmHA9WXn311XavVVRUWDaDAwAAsEva0mWnZMJ2++QRdD7h5pZ5vwEkgttzWBxdugwAAJBMBCsAAMDVCFYAAICrOZ5gm+7cmD9gn4d0YxsRHu8ZAKcl8nMm1fkvjKwAAABXI1gBAACuRrACAABcjZyVNES+AwAgkzCyAgAAXI1gBQAAuBrTQBGkerkWAACZjpEVAADgagQrAADA1QhWAACAqxGsAAAAVyNYAQAArkawAgAAXI1gBQAAuBrBCgAAcDWCFQAA4Gppv4OtMUaS5Pf7U9wSAAAQrQt/ty/8HQ8n7YOVkydPSpIGDhyY4pYAAIBYnTx5Uvn5+WHLeEw0IY2LBYNBHT58WD179pTH40l1cxzn9/s1cOBA1dXVKS8vL9XNSapM7Xum9lui75nY90ztt5R5fTfG6OTJkyouLlZWVvislLQfWcnKytKAAQNS3Yyky8vLy4hf5rZkat8ztd8Sfc/Evmdqv6XM6nukEZULSLAFAACuRrACAABcjWAlzXi9Xi1YsEBerzfVTUm6TO17pvZbou+Z2PdM7beU2X2PJO0TbAEAQOfGyAoAAHA1ghUAAOBqBCsAAMDVCFYAAICrEawAAABXI1hxscWLF2vUqFHq3r27evXq1WYZj8fT6rFy5UpLmd27d2v06NHKzc3VxRdfrEWLFkX1xVGpEk2/Dx06pAkTJqhHjx4qKCjQgw8+qKamJkuZdOt3WwYPHtzq/X300UctZaJ5LdLVihUrNGTIEHXr1k3Dhw/Xf/zHf6S6SQlVXV3d6v31+Xyh68YYVVdXq7i4WLm5ubr++uv16aefprDFHff+++9rwoQJKi4ulsfj0RtvvGG5Hk1fA4GAHnjgARUUFKhHjx665ZZb9NVXXyWxF7GL1O+pU6e2+h0YOXKkpUw69jvRCFZcrKmpSRMnTtQ//uM/hi1XU1Oj+vr60GPKlCmha36/X2PHjlVxcbF27typZ599VsuWLdPy5cudbn6HRep3c3Ozxo8fr9OnT+uDDz7QmjVr9G//9m+aM2dOqEw69rs9ixYtsry/8+fPD12L5rVIV2vXrtWsWbM0b9481dbW6tprr1VlZaUOHTqU6qYl1OWXX255f3fv3h269utf/1rLly/Xc889p507d8rn82ns2LGhL3BNJ6dPn9aVV16p5557rs3r0fR11qxZev3117VmzRp98MEHOnXqlG6++WY1Nzcnqxsxi9RvSaqoqLD8Dqxbt85yPR37nXAGrldTU2Py8/PbvCbJvP766+0+d8WKFSY/P9+cPXs2dG7JkiWmuLjYBIPBBLc0sdrr97p160xWVpb5+uuvQ+f+9V//1Xi9XtPY2GiMSe9+t1RSUmL+5V/+pd3r0bwW6eqv//qvzf33328595Of/MQ8+uijKWpR4i1YsMBceeWVbV4LBoPG5/OZJ598MnTu7NmzJj8/36xcuTJJLXSG/XMrmr6eOHHCZGdnmzVr1oTKfP311yYrK8u88847SWt7PNr6vJ4yZYq59dZb231OZ+h3IjCy0glUVVWpoKBA11xzjVauXKlgMBi6tn37do0ePdqyI+K4ceN0+PBhHThwIAWtjd/27dtVWlqq4uLi0Llx48YpEAho165doTKdpd9PPfWU+vbtq2HDhmnx4sWWKZ5oXot01NTUpF27dqm8vNxyvry8XNu2bUtRq5yxd+9eFRcXa8iQIZo8ebK+/PJLSdL+/fvV0NBgeQ28Xq9Gjx7d6V6DaPq6a9cunTt3zlKmuLhYpaWlaf96bNmyRf3799ell16q6dOn68iRI6FrnbnfsUj7b13OdI899phuvPFG5ebm6t1339WcOXN09OjR0FRBQ0ODBg8ebHlOYWFh6NqQIUOS3eS4NTQ0hPpwQe/evZWTk6OGhoZQmc7Q74ceekhXXXWVevfurY8++khz587V/v379eKLL0qK7rVIR0ePHlVzc3OrvhUWFqZ1v+xGjBihl19+WZdeeqm++eYbPf744xo1apQ+/fTTUD/beg0OHjyYiuY6Jpq+NjQ0KCcnR717925VJp1/JyorKzVx4kSVlJRo//79+tWvfqUbbrhBu3btktfr7bT9jhUjK0nWVkKd/fHxxx9HXd/8+fNVVlamYcOGac6cOVq0aJGWLl1qKePxeCzH5r+TTO3nnZTofrfVdmOM5bwb+t2WWF6Lhx9+WKNHj9bPfvYzTZs2TStXrtTq1at17NixUH3RvBbpqq33sDP064LKykr9/Oc/1xVXXKGbbrpJb731liTppZdeCpXp7K9BSx3pa7q/HpMmTdL48eNVWlqqCRMm6O2339af/vSn0O9Ce9K937FiZCXJqqqqNHny5LBl7CMCsRg5cqT8fr+++eYbFRYWyufztYq+Lwwx2v8X46RE9tvn8+nDDz+0nDt+/LjOnTsX6pNb+t2WeF6LC6sE9u3bp759+0b1WqSjgoICdenSpc33MJ37FUmPHj10xRVXaO/evbrtttskfT+iUFRUFCrTGV+DCyugwvXV5/OpqalJx48ft4wyHDlyRKNGjUpugx1UVFSkkpIS7d27V1Lm9DsSRlaSrKCgQD/5yU/CPrp169bh+mtra9WtW7fQkt+ysjK9//77ljyHDRs2qLi4OK6gKFaJ7HdZWZn27Nmj+vr60LkNGzbI6/Vq+PDhoTJu6Hdb4nktamtrJSn0gR7Na5GOcnJyNHz4cG3cuNFyfuPGjZ36AzoQCOjzzz9XUVGRhgwZIp/PZ3kNmpqatHXr1k73GkTT1+HDhys7O9tSpr6+Xnv27OlUr8exY8dUV1cX+jeeKf2OKGWpvYjo4MGDpra21ixcuNBcdNFFpra21tTW1pqTJ08aY4x58803zapVq8zu3bvNvn37zAsvvGDy8vLMgw8+GKrjxIkTprCw0Nx5551m9+7d5rXXXjN5eXlm2bJlqepWRJH6ff78eVNaWmpuvPFG88knn5hNmzaZAQMGmKqqqlAd6dhvu23btpnly5eb2tpa8+WXX5q1a9ea4uJic8stt4TKRPNapKs1a9aY7Oxss3r1avPZZ5+ZWbNmmR49epgDBw6kumkJM2fOHLNlyxbz5Zdfmh07dpibb77Z9OzZM9THJ5980uTn55vXXnvN7N6929x5552mqKjI+P3+FLc8didPngz9W5YU+t0+ePCgMSa6vt5///1mwIABZtOmTeaTTz4xN9xwg7nyyivN+fPnU9WtiML1++TJk2bOnDlm27ZtZv/+/Wbz5s2mrKzMXHzxxWnf70QjWHGxKVOmGEmtHps3bzbGGPP222+bYcOGmYsuush0797dlJaWmmeeecacO3fOUs9//dd/mWuvvdZ4vV7j8/lMdXW1q5fvRuq3Md8HNOPHjze5ubmmT58+pqqqyrJM2Zj067fdrl27zIgRI0x+fr7p1q2b+au/+iuzYMECc/r0aUu5aF6LdPWb3/zGlJSUmJycHHPVVVeZrVu3prpJCTVp0iRTVFRksrOzTXFxsbnjjjvMp59+GroeDAbNggULjM/nM16v11x33XVm9+7dKWxxx23evLnNf9dTpkwxxkTX1zNnzpiqqirTp08fk5uba26++WZz6NChFPQmeuH6/d1335ny8nLTr18/k52dbQYNGmSmTJnSqk/p2O9E8xiTZlt6AgCAjELOCgAAcDWCFQAA4GoEKwAAwNUIVgAAgKsRrAAAAFcjWAEAAK5GsAIAAFyNYAUAALgawQoAAHA1ghUAAOBqBCsAAMDV/j+uelwyo3kt/AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "grid_gdf = vectorize(geo_grid.pop_est.astype(\"float32\"))\n", "grid_gdf.plot(column=\"pop_est\")" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.15" }, "vscode": { "interpreter": { "hash": "f7286e25aede83087e9621ce0e75f07eda2aa917ee0205c1a0b96e7ba4f20042" } } }, "nbformat": 4, "nbformat_minor": 4 }