{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - Use a rasterize function when creating a raster" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "import geopandas\n", "\n", "from functools import partial\n", "from geocube.rasterize import rasterize_image\n", "\n", "from geocube.api.core import make_geocube\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "gdf = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 22, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEXCAYAAACpuuMDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA03ElEQVR4nO3deZxcVZn/8c+3k5BAQtjCEhNiwhIUEFQaEBdAthAcQUcRcAPEiSA66G9UQGYUf4iDoqM4ihgUBUU2F8hPIRGQxQWEhCUQMBIgQiCShH0NSffz++Oc27ldXcutqltd2/N+veqVu99T3ZU6fc9zznNkZjjnnHOJnmYXwDnnXGvxisE559wgXjE455wbxCsG55xzg3jF4JxzbhCvGJxzzg3iFYNrWZLOk/RfGY+9UdLHG1CGEyQ9IekFSZvlfX3nWpFXDK4lSDpG0p/S28zseDM7o4llGgX8D3CQmY0zsydzvPb+kv4m6SVJN0h6bV7Xdq5eXjG4riVpZIVDtgTGAItyvu8E4NfAfwGbAvOBy/K8h3P18IrBVU3SyZIek/S8pMWS9o/bT5f0S0mXxX13SNo1dd4pkh6M++6T9N64/fXAecBescnmmbj9p5K+Gpc3kfRbSSslPR2XJ1dZ7mMk/VnStyU9BZwuabSkb0p6JDYZnSdpfUnTgcXx1Gck/aHen1vKvwKLzOwKM3sFOB3YVdLrcryHczXzisFVRdIOwKeA3c1sQ2AGsDR1yGHAFYS/hH8BXBmbZAAeBN4BbAR8Bfi5pIlmdj9wPHBLbLLZuMite4CfAK8FpgAvA9+r4S3sCTwEbAGcCXwdmA68EdgOmAR8ycz+DuwUz9nYzPYrdjFJz5R5nVKiDDsBdycrZvYi4WezU4njnRtWXjG4avUBo4EdJY0ys6Vm9mBq/wIz+6WZrSG0z48B3gIQ/0J+3Mz6zewy4AFgjyw3NbMnzexXZvaSmT1P+FLfp4byP25m/2tma4FXgH8DPmtmT8Xrfg04MuvFzGzjMq+zSpw2Dni2YNuzwIY1vB/XQJIukLRC0r0Zjn2tpOslLYydIap6om0lXjG4qpjZEuAzhOaPFZIulfSa1CGPpo7tB5YBrwGQ9FFJdyV/UQM7AxOy3FfSBpJ+KOkfkp4DbgY2ljSiyrfwaGp5c2ADYEGqTHPj9kZ6ARhfsG088HyD7+uq91Pg4IzHfhO4yMx2Af4v8N+NKlSjecXgqmZmvzCztxOadYzQHJPYOlmQ1ANMBh6PvW7OJzRDbRabi+4FlFy2wm3/A9gB2NPMxgN7J7eptvip5VWEJqmdUn/lb2Rm47JeLMZESr2+WOK0RUA69jIW2Jacg9yufmZ2M/BUepukbSXNlbRA0h9TsaEdgevj8g2EZtW25BWDq4qkHSTtJ2k0oSnmZULzUmI3Sf8ae/x8BlgN3AqMJXwpr4zXOZbwxJB4Apgsab0St94w3usZSZsCX673vcQnmvOBb0vaIpZrkqQZVVxjXJnX10qc9htgZ0nvkzQG+BKw0Mz+Vu97csNiNvBpM9sN+Bxwbtx+N/C+uPxeYMN2HfviFYOr1mjgLMJf2/8kBHHTfxlfBRwBPA18BPhXM1tjZvcB3wJuIVQCbwD+nDrvD4S/mP8paVWR+34HWD/e91ZCk08eTgaWALfGJqrrCE8mDWNmKwlfIGcSfk57UkVcwzWPpHHAW4ErJN0F/BCYGHd/DthH0p2E+NdjwNpmlLNe8ol6XF4knQ5sZ2YfbnZZnMuLpKnAb81sZ0njgcVmNrHCOeOAv5lZWwag/YnBOecyMrPngIclHQ6gYNe4PCHG1QBOBS5oUjHr5hWD6yhxgFqxQPB5zS6baz+SLiE0f+4gaZmk44APAcdJupvQ/JkEmfcFFkv6O2HU/JlNKHIuvCnJOefcIP7E4JxzbpBKScTaxoQJE2zq1KnNLkZLeeDOh3O5zvZvmpb52L8veCiXe6ZN322b3K/putuCBQtWmVldAxlnvHOsPflUX+UDgQULV88zs6wD5ZquYyqGqVOnMn/+/GYXo6XMGPvRXK4zb/5FmY89sOfwXO6Zdu38K3K/putukv5R7zVWPdXHX+dl63Q0auKDmUb4t4qmNyVJ+qykRZLulXSJpDGSNpV0raQH4r+bNLuczjk3mNFn/Zle7aapFYOkScC/A71mtjMwgjDQ5xTgejPbnjDEvFSWSuecawoD+rFMr3bTCk1JI4H1Ja0hJDR7nNAHeN+4/0LgRsIIVVdEPU1G814s3UyUNAs1onmoGsn9r+1vTJNSsfeX3KvR926GmdOH/le65u9fL3Kkq6Sf9nsayKKpTwxm9hghI+EjwHLgWTP7PbClmS2PxywnpF0YQtIsSfMlzV+5cuVwFds55zCMPsv2ajfNbkrahDA4ZBohNfNYSZnTKZjZbDPrNbPezTdvdKZk55xbx4A19Gd6tZtmNyUdADwck4oh6deEBFVPxJm9lkuaCKxoZiGdc66YdowfZNHsiuER4C2SNiCkVN6fMDH6i8DRhCyeRxMydnaMwphAsXb+5JjCfaW2VytLbKFVNLp9v9z1WyG2kMQE8ooDlLtOsfhDlvO6kUFbNhNl0dSKwcz+KumXwB2E9LR3EnKdjwMuj3lJHgFa65vKOeegDRuJsmn2EwNm9mWGTrqymvD04JxzLckw+rwpyTnnXMIM1nRmveAVQ6tKYgCVYg1ZrpH1+ESWdvVWi0O0q3Sbfqk2/Gra9svFCOpV6trp8tUSD8k7hjJ8RF/VU463B68YnHOuBgb0+xODc865tE59YuiYiXp6e3utluyqex92dsVjbr7q87UUqWrVNg8VnlfY/FSrarrC5t2k1ArdQ4spfJ/DWc5yKSwa2XTUyuptdpK0wMx667nGTrusZ5f+rmhShiF2mfJY3fcbTv7E4JxzNeq3znxi8IrBOedq0I94lRHNLkZDeMXgnHM18ieGFvfAvcuYOf3kTG2PWeIKWY5f/7p7qrpOHsrFD6qJURTrBptsy5KyIy+tFFM4cMQRpXdqcL7JLG37xT6L3RoTyFsr/ByNzg0+d0zF4Jxzw0v0WdMnwWwIrxicc64GYQY3rxicc86leFNSi9t+58lcMz9b3+b0uIRq4w2dplxcIktsoXAKzHLH1KrRYwgGrq+hf/1d23dZOKZc/KGEZqSncOUlcR/pG3Vfy0ysMe+V5JxzLgrBZ29Kcs45N8CDz84551I8+NxAkjYGfgTsTPhZfwxYDFwGTAWWAh8ws6fzumenxhUaOd1nJVliDfVeOy9DyphzbME1TrEcUVnyRjUqpXdfhw5wa4Xq7hxgrpm9DtgVuB84BbjezLYHro/rzjnXMgyxxkZmelUi6QJJKyTdW2L/hyQtjK+/SNo19zeU0tSKQdJ4YG/gxwBm9qqZPQMcBlwYD7sQeE8zyuecc6Ukwecsrwx+ChxcZv/DwD5mtgtwBjC77jdQRrObkrYBVgI/iTXgAuAkYEszWw5gZsslFc1tK2kWMAtgypQpmW9aLo120sxUTartdNNUM9Jk5KVU+u5qmpianeLi4AmzAJi7qvL/m4HmrzLNRQNNSrH5aeR229RbRFcgr2ae4Z4BzlBuTUlmdrOkqWX2/yW1eiswOZcbl9DspqSRwJuBH5jZm4AXqaLZyMxmm1mvmfVuvvnmjSqjc84V1U9PphcwQdL81GtWHbc9Drgmn3dQXLOfGJYBy8zsr3H9l4SK4QlJE+PTwkRgRdNK6JxzRZhRTXfVVXlM1CPpnYSK4e31Xqucpj4xmNk/gUcl7RA37Q/cB8wBjo7bjgauakLxnHOuDNGf8ZXL3aRdCD04DzOzJ3O5aAnNfmIA+DRwsaT1gIeAYwkV1uWSjgMeAfLvA1lCLdN4FjsniTsMZ8whr3TZhddJrzcyBfdwKRUHSeIJaYVdW9c+uBSAkdtOzbtYXaOeWEArpQIx4NUMPY7yIGkK8GvgI2b290bfr+kVg5ndBRR7xNp/mIvinHOZGcptoh5JlwD7EmIRy4AvA6MAzOw84EvAZsC5kgDWNnIO6aZXDM45167yypVkZkdV2P9x4OO53CwDrxicc64GBvR7riRXjSTukGW6zUYpdu96YwS1jG1otHQcIIkfJOMZEv3PvTCwPGLTjQftG7n5ZkDxn1fP+uuH819ZPWh7EmuA/OMNfUsfrfqcEVO3zrUM1fjdzb8Zsu1de7+35PGl4gTVxB6Ge8xCcfL5GJxzzq3jTwzOOecG8Yl6nHPODeHzMbiaFGuLb7W4QymtFEcoJ0t+pp7x40ruW7ty6FihnjGjB62PGLtByfPtn4MH5muroqm9gNriB1mUu24Sf1i75KGarl2YHyq5zrzH7y55ThJ3KBdrqEVrxBaCMB+Dxxicc84N8BncnHPOpYTgsz8xuJy8fMAbhv2e7ZwOPG9J09GIjcYP+rfvmWcGjkm6p2pECC7amrXh37VrgMHNVwePPxaAuc/9BMiWtuH3r/4CgIPW+2BtbyKDuY/Mj0vh3xmvqW1ul1qboKB8V9YsM6+1sjBRjwefnXPOpficz8455waEtNvelOSccy7FYwyurSVxjSTWUC7OUUvq8WZK2vhhXTt/OUl67YM3/beKx1pf3+Bzi3SNzXLPRLnYgkYV/++YxDeyWBdXGKpY99Jq4g7luqdWku62OqTL6do+2lHIrupNSc4551K6MleSpDkZrvGUmR2TT3Gcc649GGJtf3f2Sno95XOAC/h+vYWQNILQp+4xM/sXSZsClwFTgaXAB8zs6Xrv45xzeerWkc+nmdlN5Q6Q9JUcynEScD8wPq6fAlxvZmdJOiWut2dn5xZTGFtot3hCMX0vvjSwnI43QPn2/7lPnT9oPR1zGBjToNCGXDjt58ytPllVGZOxEweOOGLQ9ZJ1GBrPSNJyFIs9FKYBn/fYnVWVZ+C8OuIGiZ6t1s002f/P6QDMnLp7uP4rFwNwTZnJKK956Jt1l6EZOrlXUtnIiZldXukCWY4pR9Jk4F2ESa4ThwEXxuULgffUcw/nnGuEfuvJ9Go3mYLPkqYDnwdemz7HzPbLoQzfAb4AbJjatqWZLY/3WC6paFYySbOAWQBTpkzJoSjOOZdNnnM+t5qsvZKuAM4Dzgdy61sm6V+AFWa2QNK+1Z5vZrOB2QC9vb2WV7mccy6Lbo0xJNaa2Q8acP+3AYdKOgQYA4yX9HPgCUkT49PCRGBF2as4V4cktlAYc4DBMQCoPraQSKYPzbI9iUckcYTCFOBptcYW8pTEFWBdvEGjP1bxvJnbfA5o4xgDdGyvpLKNX5I2jT2E/p+kT0qamGyL2+tiZqea2WQzmwocCfzBzD4MzAGOjocdDVxV772ccy5XFpqSsrzaTaUnhgWEijF5Z+kuLAZsM+SMfJwFXC7pOOAR4PAKxzvn3LDKc6IeSRcASdP6zkX2CzgHOAR4CTjGzO7I5ebFymNWuWle0hgze6XStmbq7e21+fNLpwNoRXsfdnbT7t0J3VSLGZihTpX/w8574cKKx1SjMM2F1hs1sNwzbmzm6/Q//0LFY2yn7QD45ZWhM9+4njGZr5+3dHfVUpJmo7RmNiFJWmBmvfVcY5PXbWH7/jjb36xXvv3csveTtDfwAnBRiYrhEODThIphT+AcM9uzpoJnkLUf1V8ybnPOua6QTNSTR1OSmd0MPFXmkMMIlYaZ2a3AxjH+2hCVUmJsBUwC1pf0JtY1KY0HSk+C65xzHS6kxMg8RmGCpHSTxuzYqzKrSUB6Yu9lcdvyKq6RWaUYwwzgGGAy8C3WVQzPAV9sRIGcc65dVBFjWFVn01WxGzWsi37ZisHMLpT0M+AoM7u4UYXoVkk7f7H21xffsCUAY+95YtB6XvfsJANxhSoVdkUtTHuRPmbkxPDzv2bZd0ter1TabID+F14Eqos1FEriCq2mWHfVRLt3SS3LhnU+hmXA1qn1ycDjjbpZxecgM+sHPtGoAjjnXDvKM8aQwRzgowreAjybZIdohKwD3K6V9DlCxtMXk41mVi5Y4pxzHS2vJwZJlwD7EmIRy4AvA6MAzOw84GpCj6QlhO6qxxa/Uj6yVgzJMMYTU9saOY7BOedaWp65kszsqAr7jcHfvw2VqWIws2mNLkg3K9b+Whh3SGINiXTMIUvcoJljJhqlmjELWRTGHIqZOfnfh2wrGXcoMkbomn+eG65TJrVGz4bjAOibslXJY5LxC4kX+sOQosO33WfdvR7+a9Fzk2Oh8eMfOjK2kNLXhplTs8iaXXUUcAKwd9x0I/BDM1vToHI551xLs+ENPg+rrE1JPyC0d50b1z8St5Wb3c055zqadXnFsLuZ7Zpa/4Ok+qd+cs65ttWeCfKyyFox9Ena1sweBJC0DTnOy+CGStpmi41xgKExh0o6afxCElvof/nlIft6Nhj+AfnF4g6lJGWf9+JFpa9XY2pvgCseXDcT7wv9lY8/ZIe3A3D14j/VfM8suZI6Vbc/MXweuEHSQ4QReK+lwd2lnHOulSXjGDpR1l5J10vaHtiBUDH8zcxWVzjNOec6l0FfN1cM0W7A1HjOrpIws9LPwy4XpZqUau0GeNCoIwH4/ZpL6ytYExVrQkrYq4M7yqVTXxfqKUhhYX3rWketv/o0NEnzUJYUHeWaw67tv6LoOc89Xt+85uNf8wgAB29Uena1d+02Y8i23y2YN2i9m5uO0owub0qK+ZK2Be5iXWzBAK8YnHNdyoPPvcCOlmVWnypI2ppQuWwF9BNS0Z4Tpw29jPCEshT4gJk9nee9nXOuXvl+I7aOrMP27iV8eedtLfAfZvZ64C3AiZJ2BE4Brjez7YHr47pzzrUUM2V6tZusTwwTgPsk3QYMBJ3N7NB6bh6zAy6Py89Lup8w+cRhhIRSABcSRlqfXM+9Ok065lAYbyiX6rhcbCFJCVEs/XS7SmIOg2INa9dWPE89g/8z9z/zbOZ7VhNrSCsVW0gkMQLIFm9IH1/JITvvC0Df0ytCWVKfgSSttscWBjODvuwT9bSVrBXD6Y0sBICkqcCbgL8CWyYpZc1suaQtSpwzC5gFMGVKfYE555yrVqc2JWXtrnpTuf2SbjGzvWothKRxwK+Az5jZc8qYFC1OjTcboLe3t0N/Rc65VtWOzURZ5PUcVHOKxpig71fAxWb267j5iWSi6/jvivqL6Jxz+TGyxRfasfKoZhxDOTX9ta7waPBj4H4z+5/UrjnA0cBZ8d+r6i5hmysVRyhcBli79B813aPVYwsH9hxe+8kZ4gr1GhKjyTAWolJcoVrVxBWKGbHJRsDgFOTX9pWPLczc+qSB5WsePaeu+7ebTm2myKtiqNXbCJla75F0V9z2RUKFcLmk44BHgDq+EZxzrgGsc5uS8qoYavrpmNmfypy7f+3Fcc65xrP+zqwYMsUYJI2V1BOXp0s6NMYGEh9pSOmcc66FmWV7tZusTww3A++QtAlhwNl84AjgQwBmdm9jiudKKZcrKWkfHtxO3Nrxg2poxIi4UPnvmoHxCDttt27joiXZ7zV6dOZjR06aCMDMKZ8pecy8l3+W+XrlJLGEWvInzX32giHbDp4wC4D+514AUj/jDLotrpDo5FxJWXslycxeAv4V+F8zey+wY+OK5ZxzLc4AU7ZXm8lcMUjai/CE8Lu4rdmBa+eca6pub0r6DHAq8BszWxRncLuhYaVydUmajbJ070ynYC5MWZ00e6Svk3f3ymok905Sh2NxirIiTUpJU4itjSkx0jvTzUpQtmmp/6WXgPIzw/Vsuknxc1e/MmRbq6Ydmbtq9qD1mZM+3aSStBPlFnyWdDBwDjAC+JGZnVWwfyPg58AUwvf2N83sJ7ncvIhqRj7fJGlsXH8IyD6foXPOdaIcngYkjQC+DxwILANulzTHzO5LHXYicJ+ZvVvS5sBiSReb2av1l2CorL2S9pJ0H3B/XN9V0rmNKJBzzrUFyy276h7AEjN7KH7RX0pIJFpwNzaMg4LHAU8RslMXJemkLNtKyRpj+A4wA3gSwMzuBvbOehPnnOtIlvEFEyTNT71mpa4yCXg0tb4sbkv7HvB64HHgHuAks6Qttaiji2w7JtubqiKAbGaPFiS36yt1rGsN5eIB5aZ3TMxYPwxP6RldcyqshkhShyexj571SnetTKbozNISnJ7ac8i+1dVPcT5ycuH/bVj72PKqr1NOvSkwSrnmsf+t6bwkPUb3dGHNHGNYZWa9VVyksJFqBmEGzf0Is2leK+mPZvbcoAtJRwEfBKZJmpPaNZ74h30WWSuGRyW9FTBJ6xHiC/dnvYlzznWkfHocLQO2Tq1PJjwZpB0LnBVn0Vwi6WHgdcBtBcf9hTDHzQTgW6ntzwMLsxYoa1PS8YTgxyTCm3hjXHfOue5kQL+yvcq7Hdhe0rT4h/eRhESiaY8Q0wRJ2hLYAXhoSJHM/mFmNwIHAH+MHYeWEyqbzI83WXslrSKOcnbOORfkMUbBzNZK+hQwj9Bd9YI4LOD4uP884Azgp5LuIXzBnxy/l0spm62ikkwVg6TpwA8IM6vtLGkX4FAz+2qW8117SKa/LBzPkFY4NmJE7MNf2A9+OCQxlBlj1n3W571yccXzZuz25bBQZvzCyC02B2DtipUAjNhsUwD6nnxqyLH2bGjm1UbjB20vFk9otfEL++/33wBc/4dTM5+TxJ4S6TQf3RNbiHIavGZmVwNXF2w7L7X8OHBQFZeUmb0UM1T/r5l9Q9KdWU/O2pR0PmGA25pYyIWExx3nnOterZsSo65sFVkP3MDMbivoldT4mU+cc66FqXXTXXyGOrJVZK0YVknalvjgJOn9hICGc851p3VjFFpOvdkqslYMJwKzgddJegx4mAYHoyvlDnENtHPMJXTvksHrAAtChvWe3XYGWuP/RZa4wqDjF3wFyJZLKok1JJJYg7308rqNMQKZxBqKabXYQiKJLSSxhpG3hCwM6bhB4b61e+046NzulanHUVPEZqQfE0ZJT5G0K/AJM/tklvMrVgwxj8cJZnZArH16zOz5egqd8Z6Vcoc451xztcJfRsV9hzAobg6EbBWSMmerqBh8NrM+YLe4/GKjK4UoS+4Q55xrruwpMYa/aGaPFmzKnK0ia1PSnXF49RXAi6kb/zrrjapULHfInoUHxXwjswCmTKl+JqtudNB6Hxy03rN+mXQXO283ZFPShFRoxu5fGbJt3u1frq5w7aTYDGf9g1PXtGrzUTGFTUpZ9xVuL9W8lOWYtpNM1NOa6spWkbVi2JSQZ2O/1DYDGlUxZMkdgpnNJsQ+6O3tbd2HOudcR2rhXknHE2K0SbaK31NFtoqsI5+PralotcuSO8Q555qrRSuGStkqJJ1qZsUfAck+8vm7RTY/C8w3s6uyXKNKA7lDgMcIg+k+WP4U55wbXi38xFDJ4UB9FQMwhpDJL8nj/D5gEXCcpHea2WfqKWGhUrlD8rxHtxuY+jKd/qJITKEehXGHVos5JCk1km6rhV1Tq9YT+nLMffpH9V2nRZWKDWSJGaSPKYxV9Nx0F9BeMZkBrRtjqKRswbNWDNsB+5nZWgBJPyC0WR1ImDQid8VyhzjnXMto4QFuGZQtedZcSZOAsan1scBrYlfW6mcwcc65TtDC3VUryOWJ4RvAXZJujBfcG/haHPB2XV3Fc865NtXGMYbS0zuSvVfSjyVdTRh4JuCLMQ0swOfrK58bTirsf59zXKGcJObQarGGZOrS/meHjt3s2WjD4S5OUxWLFxTGBOodh1BqXMSBI44YWG6beEO5WZebKCbNOwfYi1DKW4DPxpxJmNnXyp2fqSlJIa3q/sCuZnYlMFLSHnWU2znn2pos+6sJfgFcDmwFvIbwhHBJ1pOzxhjOJdQ8R8X15wm5jJxzrnu18HwMZvYzM1sbXz+nimhH1hjDnmb25mQGIDN7Og6zds657tW6MYYbJJ1CyDNnhGk9fydpUwAzGzoVYUrWimFNzHiazMewOS3buubKSdIpD0zPeG9qessYb7C7QkoVvfH1w1q2dtL/4ksDyz1jN2hiSRqvUbmNOiFnUgsHn5OAzScKtn+M8D2+TbmTs1YM3wV+A2wh6Uzg/cB/VlFI55zrPC1aMZjZtHrOz9or6WJJCwgBaAHvMbPMmfqcc67jGKhF200kjQJOIAwtALgR+KGZrSl5UkrZiiFpj4pWkIpqS9q0UjuVa11DmpRgcLMS65qU0uppXmq1bqpZ9D31TMVjkmalpCtw8jNNz4KWmDEm5DWrdta5dpW837TC997WP5MWfWIAfgCMInQcAvhI3PbxLCdXemJYQHjrAqYAT8fljYFHgLoeV5xzrp3lFWPIMpWxpH0JM7ONAlaZ2T5lLrm7me2aWv+DpLuzlqdsd1Uzm2Zm2xCS2b3bzCaY2WbAv9C4uRicc65rpKYyngnsCBwlaceCYzYm/PV/qJntRMiOWk6fpG1T529DA2Zw293Mjk9WzOwaSWdkvYlzznWkfJ4YBqYyBpCUTGWcnuP+g8CvzewRADNbUeGanyd0WX0ork8FMs+rk7ViWCXpP4FkkMSHCTO6uTaXbgcfFG8ooeu7strQaKNGjhq0Xi620OmyvM9SMYX0uW0Rb8gv+JxlKuPpwKiYr25D4Bwzu6jMNf8M/JDQYYi4fEvWAmWtGI4CvkzosmrAzawbBe2cc90p+xPDBEnzU+uz49TEkG0q45HAboQv+vWBWyTdamZ/L3G/i4DngKRl5yjgZ1Rughq4WUWx99FJWY51zrluIKoKPq8ys94S+7JMZbwsXuNF4EVJNwO7AqUqhh0Kgs835BZ8lnR6pQtkOabEeWdL+pukhZJ+E4Mryb5TJS2RtFjSjFqu75xzDZfPfAwDUxnHVENHAnMKjrkKeIekkZI2IDQ1lRtLdqektyQrkvYkNC9lUumJ4eOSniuzX4Q3cXrWG6ZcC5wap/H8OnAqcHKMxh8J7ETICnidpOlxUiDXQEnb+EGjjmxySYZX/+pXMh87JG15SrHYwpBj2qHtvEq1xk/aPu6SU+bUUlMZSzo+7j/PzO6XNBdYSEhH9CMzu7fMZfcEPirpkbg+Bbhf0j3hkrZLuTJVqhjOJwQ6Kh1TNTP7fWr1VkKaDQjR+EvNbDXwsKQlhKh95sCJc84Ni5zGMRSbytjMzitYPxs4O+MlD66nPGUrBjP7Srn9OfoYkMzMMYlQUSSWxW1DSJoFzAKYMmVKI8vnnHNDtGpKDDP7Rz3nZ+2VVBNJ1xEmiih0mpldFY85DVgLJM/YWSL0YWOI6s8G6O3tbd3B6c65ztSh3zoNrRjM7IBy+yUdTRhFvb+ZJT/iLBF6N4ysf+inv5qpR1o9R9K1/WH62wN7Kvfks74Q6kqmAy2n7dvQSxiO99UW+ZOyBZbbUtYZ3HIXc4OcTBji/VJq1xzgSEmjJU0Dtgdua0YZnXOunBae2rMuWed8ni7pekn3xvVd4kjoenyPENi+VtJdks4DMLNFhLlK7wPmAid6jyTnXEvKp7tqy8nalHQ+IffGDwHMbKGkXwBfrfXGZrZdmX1nAmfWem1Xm4PW+yAAPeuvX/HYvjtCGpeeNw/K9dXyzUb1GrHRRkO2zX36R0WPbelmkDbR6k1K7fg0kEXWimEDM7tNGtSyvLYB5XHOufZgdOwEx9Uk0duWdXM+vx9Y3rBSOedcixPVdcJoJ1krhhMJ3UJfJ+kx4GFChlXnnOte3dyUFPOEHyBpLNBjZs83tliu1RXGFjpBOn7Q9+yzQ7Y5V6grYwyS/k+J7QCY2f80oEzOOdceurFiYF2epB2A3VmX8e/dhDkZnHOuO+U3UU/LyZQrSdLvgTcnTUgx1fYVDS+dc861si59YkhMAV5Nrb9KmEPUdan+BSHjb89uOwOdMX4hPR7h4E0+nvnYbpOMKejUlB/V6MoYQ8rPgNskJVN7vpcwdZxzznWvbq4YzOxMSdcA74ibjjWzOxtXLOeca31d/cQgaQqwCvhNepuZPVL6LOec62Btmgcpi6xNSb9j3Y9gfWAasJgw/abrYp0QWyinm2MJWdWaxyhLjKJVcyRBHPncjb2SEmb2hvS6pDcDn2hIiZxzrl10+RPDIGZ2h6Td8y6Mc861E1ln1gxZYwzpEdA9wJuBlQ0pkRtWSaptgJ4xozOfl8x61qm8Camyept5hrPbazI7X66fW48xDIyAhpBu+3fAr/IvjnPOtY9O7ZWUdWrP+8zsK/F1ppldTEiLUTdJn5Nkkiaktp0qaYmkxZJm5HEf55zLm/qzvSpeRzo4ft8tkXRKmeN2l9QXpz5omKwVw6kZt1VF0tbAgcAjqW07AkcSejwdDJwraUS993LOudzlMLVn/H77PjAT2BE4Kn4PFjvu68C83MpfQqXsqjOBQ4BJkr6b2jWefGZw+zbwBeCq1LbDgEvNbDXwsKQlwB7ALTncz9Vg7nM/aXYRXAdr5S6pZVluTUl7AEvi9AZIupTwPXhfwXGfJjThN7zjT6UnhseB+cArwILUaw5QVxOPpEOBx8zs7oJdk4BHU+vL4rZi15glab6k+StXeizcOTfMcnhiIMN3nqRJhFRE5+VQ6ooqZVe9G7hb0sVmVvUTgqTrgK2K7DoN+CJwULHTihWlRPlmE2aWo7e3t0PDQM65ViSqemKYIGl+an12/P5KLlWo8MrfAU42s75kPpxGqtSUdLmZfQC4Uxr6IzCzXcqdb2YHlLjuGwijp++Ob3IycIekPQi15dapwycTnlycc661ZB/HsMrMekvsy/Kd1wtcGr8vJwCHSFprZldmL2x2lbqrnhT//Zc8b2pm9wBbJOuSlgK9ZrZK0hzgF5L+B3gNsD1wW573d+v8/tVfDCwfPP7YQfs8tuDaWTJ2ARo07ia/iXpuB7aXNA14jND55oPpA8xsWrIs6afAbxtVKUCFGIOZLY+LnzSzf6RfwCcbUSAzWwRcTgi8zAVONLO+RtzLOefqkUd31dhM/ylCb6P7gcvNbJGk4yUd3/h3MVTWAW4HAicXbJtZZFtNzGxqwfqZwJl5XNs55xomp8immV0NXF2wrWig2cyOyeeupVWKMZxAeDLYRtLC1K4NgT83smDOOdfqOnXkc6Unhl8A1wD/DaRH4z1vZk81rFRu2NSaK8m5VtWQvEjFGNUEn9tKpe6qzwLPAkcBSNoCGAOMkzTOJ+pxznWzTn1iyJQSQ9K7JT0APAzcBCwlPEk451xXSibqySNXUqvJGnz+KvAW4Doze5OkdxKfIlzn8W6qrp0NW0p4s45tSsqaRG+NmT0J9EjqMbMbgDc2rljOOdf6ZNle7SbrE8MzksYBNwMXS1pBPkn0nHOufbXhl34WWZ8YDgNeBj5LGHT2IDnNx+Ccc+2qq58YzOzF1OqFDSqLa7L+V1Y3uwjOtQ8D+trwWz+DSgPcnqf4w5IAM7PxDSmVc861gXZ8Gsii0jiGDcvtd865rtahvZKyBp+dc84V6MonBtf50mm3nXNVyDY7W1vyisE552oQZnDrzJrBKwbnnKuRurFXknPOuRI6uCkp6wC3hpD0aUmLJS2S9I3U9lMlLYn7ZjSzjM45V5yty5dU6dVmmvbEEBPxHQbsYmarY0pvJO1ImPN0J8Kcz9dJmu7TezrnWk2n9kpq5hPDCcBZZrYawMxWxO2HAZea2WozexhYAuzRpDI651xpHfrE0MyKYTrwDkl/lXSTpN3j9knAo6njlsVtQ0iaJWm+pPkrV65scHGdcy7FfD6Gmki6DtiqyK7T4r03IczzsDtwuaRtCL3AChWtcs1sNjAboLe3t/2qZedce+vvzK+dhlYMZnZAqX2STgB+bWYG3CapH5hAeELYOnXoZODxRpbTOedq0anjGJrZlHQlsB+ApOnAesAqYA5wpKTRkqYB2wO3NauQzjlXUofGGJo5juEC4AJJ9wKvAkfHp4dFki4H7iNMBnSi90hqrgN7Dh9YHrZpE51rdQbkFD+QdDBwDjAC+JGZnVWw/0PAyXH1BeAEM7s7n7sP1bSKwcxeBT5cYt+ZwJnDWyLnnMtOWC5NSZJGAN8HDiQ0pd8uaY6Z3Zc67GFgHzN7WtJMQmx1z7pvXoKPfHbOuVr15/LIsAewxMweApB0KaHb/kDFYGZ/SR1/KyH22jBNHfnsnHNtK2lKyvKCCUnX+vialbpS5i760XHANTm9i6L8icENkcQUkniCxxXyk47XlFLPzzvL9dtBu3zmqmhKWmVmvaUuU2Rb0QvHjBHHAW/PeuNaeMXgnHO1yqfHUaYu+pJ2AX4EzDSzJ/O4cSnelOScczXJLYne7cD2kqZJWo+QK25O+gBJU4BfAx8xs7835O2k+BODc87VwsjlicHM1kr6FDCP0F31AjNbJOn4uP884EvAZsC5kgDWlmmaqptXDG5AYft0YayhUxS+r3Lt8lmOyVupe6V/D50SSyilmvfXzM9nXhP1mNnVwNUF285LLX8c+HguN8vAKwbnnKtVG45qzsIrBuecq4XhSfRcfuppoil2bjVdIGu5d6ulxMirGSXLdVqpyaaVytJKCn8uw/cZbc88SFl4xeCcc7XyisE559wgXjE455wbYAZ9nZn42SuGKuTdxlvP9ao9t1RX1GKaGUfwdnRXq6Z8bv2JwTnn3ADvleScc26IDn1iaFquJElvlHSrpLtiGto9UvtOlbRE0mJJM5pVRuecK8un9szdN4CvmNk1kg6J6/tK2pGQRGon4DXAdZKmN2p6zyz9+jux3Tvv9th6Yxad+DN2w6NpMbEODj43M7uqAePj8kasSzN7GHCpma02s4eBJYQZjpxzrrX4E0PuPgPMk/RNQgX11rh9EmHqukTJ2YziLEizAKZMmdKwgjrnXFFt+KWfRUMrBknXAVsV2XUasD/wWTP7laQPAD8GDqCK2YzMbDZhUmx6e3s78zfknGtR5r2SamFmB5TaJ+ki4KS4egVhZiLIOJtRXkq1T3Z6m3et76/w59Vu+YZc52laengDs/7hvecwaWaM4XFgn7i8H/BAXJ4DHClptKRpwPbAbU0on3POlddv2V5tppkxhn8DzpE0EniFGCuIMxddDtwHrAVObFSPJOecq1kH90qSdUjwZLw2tT21f0Ou3YxZvJxzjXOd/XJBvVNjbjRigu019t2Zjp33/E/rvt9w8pHPzjlXI+vvzBiDVwzOOVeT9hyjkIVXDM45VwtPotfdPLbgnCtkgHVo8NkrBuecq4UZdOg4Bq8YnHOuRuZNSc455wbp0CeGjhnHIGkl8I8GXX4CsKpB166VlykbL1N2rViuRpXptWa2eT0XkDSXUL4sVpnZwfXcbzh1TMXQSJLmt9rgFC9TNl6m7FqxXK1Ypm7QzFxJzjnnWpBXDM455wbxiiGb2c0uQBFepmy8TNm1YrlasUwdz2MMzjnnBvEnBuecc4N4xeCcc24QrxhSJB0uaZGkfkm9qe1TJb0s6a74Oi+1bzdJ90haIum7korNWZ17meK+U+N9F0uaMVxlKijD6ZIeS/1sDqlUvuEi6eB47yWSThnu+6fKsTT+Pu6SND9u21TStZIeiP9u0uAyXCBphaR7U9tKlmE4fnclytSyn6euYmb+ii/g9cAOwI1Ab2r7VODeEufcBuwFCLgGmDlMZdoRuBsYDUwDHgRGDEeZCsp3OvC5IttLlm+Yfpcj4j23AdaLZdmxSZ+rpcCEgm3fAE6Jy6cAX29wGfYG3pz+HJcqw3D97kqUqSU/T9328ieGFDO738wWZz1e0kRgvJndYuHTexHwnmEq02HApWa22sweBpYAewxHmTIqWr5hvP8ewBIze8jMXgUujWVqFYcBF8blC2nw78jMbgaeyliGYfndlShTKc3+PHUVrxiymybpTkk3SXpH3DYJWJY6ZlncNhwmAY8WuXczyvQpSQtj00DSHFGqfMOl2fdPM+D3khZImhW3bWlmywHiv1s0oVylytDsn10rfp66Stcl0ZN0HbBVkV2nmdlVJU5bDkwxsycl7QZcKWknQlNNoar7/9ZYplL3zqVMg25UpnzAD4Az4j3OAL4FfKwR5ahSs++f9jYze1zSFsC1kv7WpHJk1cyfXat+nrpK11UMZnZADeesBlbH5QWSHgSmE/5qmZw6dDLw+HCUKd576yL3zqVMaVnLJ+l84LcVyjdcmn3/AWb2ePx3haTfEJpAnpA00cyWx+a/FU0oWqkyNO1nZ2ZPJMst9nnqKt6UlIGkzSWNiMvbANsDD8XH7+clvSX2/PkoUOov/LzNAY6UNFrStFim24a7TPELJfFeIOlhUrR8jSpHEbcD20uaJmk94MhYpmElaaykDZNl4CDCz2gOcHQ87GiG73OTVqoMTfvdtfDnqbs0O/rdSi/CB3EZ4engCWBe3P4+YBGhV8QdwLtT5/QSPrwPAt8jjiZvdJnivtPifReT6nnU6DIVlO9nwD3AQsJ/3omVyjeMv89DgL/HMpzWpM/UNvFzc3f8DJ0Wt28GXA88EP/dtMHluITQJLomfp6OK1eG4fjdlShTy36euunlKTGcc84N4k1JzjnnBvGKwTnn3CBeMTjnnBvEKwbnnHODeMXgnHNuEK8YnHPODeIVQxeR9EIDrnloktJa0nsk7VjDNW4sTCme4fjFkg4tsm9qOo1zp5P0xdTy+jFV9auSJjSzXK69ecXg6mJmc8zsrLj6HkJ65OHwITNr6EjmZLR7ixuoGMzsZTN7I54qwtXJK4YupOBsSffGCWSOiNv3jX+N/1LS3yRdHNNqIOmQuO1PCpP//DZuP0bS9yS9FTgUODv+1bpt+klA0gRJS+Py+pIujRk0LwPWT5XtIEm3SLpD0hWSxmV4P7tJulvSLcCJqe0j4vu8Pd7rE3F7j6RzFSZA+q2kqyW9P+5bKulLkv4EHF6qPPGeN8WMqfOSVA6S/l3SffF+l5Yp81iF7KG3K2TtPSxunyrpj/F+d8SfK5ImSro5/mzvlfQOSWcByVPCxZl++c5l0eyh1/4avhfwQvz3fcC1hMlstgQeASYC+wLPEhKU9QC3AG8HxhBSHk+L518C/DYuHwN8Ly7/FHh/6n43EicXAiYAS+Py/wEuiMu7AGsJaTwmADcDY+O+k4EvFXkfA9eN6wuBfeLy2cSJX4BZwH/G5dHAfMIkL+8Hro7vcSvg6aTchEl1vpAq85DyAKOAvwCbx+1HpN7P48DouLxxmd/F14APJ8cRUneMBTYAxsTt2wPz4/J/sC6dxghgw/TvtODaSymYGMhf/qrm1XXZVR0QvuwvMbM+QobNm4DdgecIifiWAUi6izB73QuEpIEPx/MvIXzp1mpv4LsAZrZQ0sK4/S2Epqg/xweV9QiVU0mSNiJ8Ad8UN/0MmBmXDwJ2SZ4GgI0IX7ZvB64ws37gn5JuKLjsZRXKswOwMyGFNoQv6uXxnIXAxZKuBK4sU/SDgEMlfS6ujwGmECqW70l6I9BHyOILISngBZJGAVea2V1lru1cXbxi6E7l5oBenVruI3xGap0zei3rmivHFOwrlqRLwLVmdlQV91CJayX7Pm1m8wZtlN5V4ZovliuPpDcAi8xsryLnvotQ8R0K/JekncxsbYmyvc8KZueTdDohWeKuhJ/dKxBmO5O0d7z+zySdbWYXVXgfztXEYwzd6WbgiNgGvznhi6xcCuO/AdtImhrXjyhx3PPAhqn1pcBucfn9qe03Ax8CkLQzoTkJ4FbgbZK2i/s2kDSdMszsGeBZSW+Pmz6U2j0POCH+lY2k6Qqpr/8EvC/GGrYkNKEVU6o8i4HNJe0Vt4+StJOkHmBrM7sB+AKhiahUjGQe8OlUDOdNcftGwPL4NPMRwtMIkl4LrDCz84EfE+ZKBliTvD/n8uIVQ3f6DaHJ427gD4Q29X+WOtjMXgY+CcyNQdknCLGIQpcCn4/B1G2BbxK+mP9CaK9P/AAYF5uQvkCslMxsJSFmcUncdyvwugzv51jg+zH4/HJq+4+A+4A7FLqw/pDwBPQrQprnZNtfi72fUuWxMIf0+4GvS7obuAt4K+FL/OeS7gHuBL4dK65iziDEKhbGsp0Rt58LHC3pVkIzUvL0si9wl6Q7CTGic+L22fEaHnx2ufG02y4TSePM7IX4F+73gQfM7NtNKsuNwOfMbH4d10jez2aEiult5SrHdhJ7f/Wa2apml8W1J39icFn9WwxGLyI0d/ywiWV5Cvipigxwq8Jv4/v5I3BGJ1QKigPcCE8i/U0ujmtj/sTgXANJOhY4qWDzn83sxGLHO9cKvGJwzjk3iDclOeecG8QrBuecc4N4xeCcc24Qrxicc84N8v8B5OE01LsLFq0AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "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": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEXCAYAAACpuuMDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1sUlEQVR4nO3deZxcVZn/8c+3k5CELEAIayAkQIKyKwHEBRh2cAQdUcANECeCuP5GBWRGcRgc3MbBUcSojIDIpggZBCIgiwsICRCWIBAgQEIEIjuEkHQ/vz/Oud23q2u5VXVr7ef9etUrt+56qrpSp+55znmOzAznnHMu0dPqAjjnnGsvXjE455wbxCsG55xzg3jF4JxzbhCvGJxzzg3iFYNzzrlBvGJwbUvSOZL+LeO+N0n6RAPKcIKkpyW9Imn9vM/vXDvyisG1BUnHSPpjep2ZHW9mp7ewTKOA/wIOMLPxZvb3HM+9r6S/SnpN0o2Stsjr3M7VyysGN2xJGllhl42AMcD9OV93MnA58G/AJGA+cEme13CuHl4xuKpJOknSMkkvS3pQ0r5x/WmSfiXpkrjtTkk7pY47WdIjcdsiSe+L698MnAPsEZtsXojrfy7pP+LyepKukvSspOfj8mZVlvsYSX+S9D1JzwGnSRot6TuSnohNRudIGitpJvBgPPQFSb+v931L+SfgfjO7zMxeB04DdpL0phyv4VzNvGJwVZG0DfBpYFczmwAcCCxJ7XIYcBnhl/AvgStikwzAI8C7gHWArwO/kLSJmT0AHA/cGpts1i1y6R7gf4EtgKnASuAHNbyE3YFHgQ2BM4BvAjOBnYGtgSnAV83sIWC7eMy6ZrZPsZNJeqHM4+QSZdgOWJg8MbNXCe/NdiX2d66pvGJw1eoFRgPbShplZkvM7JHU9gVm9iszW01onx8DvA0g/kJ+ysz6zOwS4GFgtywXNbO/m9mvzew1M3uZ8KW+Vw3lf8rM/sfM1gCvA/8MfMHMnovn/QZwZNaTmdm6ZR5nljhsPPBiwboXgQk1vB7XQJLOlfSMpPsy7LuFpBsk3RM7Q1R1R9tOvGJwVTGzxcDnCc0fz0i6WNKmqV2eTO3bBywFNgWQ9DFJdye/qIHtgclZritpbUk/lvS4pJeAW4B1JY2o8iU8mVreAFgbWJAq07VxfSO9AkwsWDcReLnB13XV+zlwUMZ9vwOcb2Y7Av8O/GejCtVoXjG4qpnZL83snYRmHSM0xyQ2TxYk9QCbAU/FXjc/ITRDrR+bi+4DlJy2wmX/BdgG2N3MJgJ7Jpeptvip5RWEJqntUr/y1zGz8VlPFmMipR5fKXHY/UA69jIO2Iqcg9yufmZ2C/Bcep2krSRdK2mBpD+kYkPbAjfE5RsJzaodySsGVxVJ20jaR9JoQlPMSkLzUmIXSf8Ue/x8HlgF3AaMI3wpPxvPcyzhjiHxNLCZpLVKXHpCvNYLkiYBX6v3tcQ7mp8A35O0YSzXFEkHVnGO8WUe3yhx2G+A7SW9X9IY4KvAPWb213pfk2uKOcBnzGwX4IvA2XH9QuD9cfl9wIROHfviFYOr1mjgTMKv7b8RgrjpX8ZXAkcAzwMfBf7JzFab2SLgu8CthEpgB+BPqeN+T/jF/DdJK4pc97+BsfG6txGafPJwErAYuC02UV1PuDNpGDN7lvAFcgbhfdqdKuIarnUkjQfeDlwm6W7gx8AmcfMXgb0k3UWIfy0D1rSinPWST9Tj8iLpNGBrM/tIq8viXF4kTQOuMrPtJU0EHjSzTSocMx74q5l1ZADa7xiccy4jM3sJeEzSBwAU7BSXJ8e4GsApwLktKmbdvGJwXSUOUCsWCD6n1WVznUfSRYTmz20kLZV0HPBh4DhJCwnNn0mQeW/gQUkPEUbNn9GCIufCm5Kcc84N4ncMzjnnBqmURKxjTJ482aZNm9bqYrTEw3c91pDzWl97303O3GXLVhehIz1839Ih62Zs35Ex0potWLBghZnVNZDxwH8YZ39/rrfyjsCCe1bNM7OsA+VarmsqhmnTpjF//vxWF6MlDhz3sYact2/lyoacNy/Xzb+s1UXoSAfPPGnIumvmf7PInt1L0uP1nmPFc738ZV62CnXUJo9kGuHfLlrelCTpC5Lul3SfpIskjZE0SdJ1kh6O/67X6nI659xgRq/1ZXp0mpZWDJKmAJ8FZpnZ9sAIwkCfk4EbzGwGYYh5qSyVzjnXEgb0YZkenaYdmpJGAmMlrSYkNHuK0Ad477j9POAmwgjVrtaoJqFy2r25qJH27/lAzcde11e5GatYk02hax5qXhNOufKU29bMMnaaPjrvbiCLlt4xmNkyQkbCJ4DlwItm9jtgIzNbHvdZTki7MISk2ZLmS5r/7LPPNqvYzjmHYfRatkenaXVT0nqEwSHTCamZx0nKnE7BzOaY2Swzm7XBBo3OlOyccwMMWE1fpkenaXVT0n7AYzGpGJIuJySoejrO7LVc0ibAM60spHPOFdOJ8YMsWl0xPAG8TdLahJTK+xImRn8VOJqQxfNoQsbOrtGKWEI3yNKuX8v5qok1VFOGLG3zWeIQ1Zyv1nPXcp5aYw+1lKsd4xwGHdlMlEVLKwYz+4ukXwF3EtLT3kXIdT4euDTmJXkCqD1K6JxzDdJ5jUTZtPqOATP7GkMnXVlFuHtwzrm2ZBi93pTknHMuYQaru7Ne8IqhWdo1rtAzdmz/8rxXz6+4fz19/9tV3rGLLKoZ45Dsm1fMIC+NjI90BtFb9ZTjncErBuecq4EBbZ5nsmZeMTjnXI269Y6haybqWWfMxvb2qUdnumXd87Bv53rtsdffm+v52kW5pqVmNik1q6ln0GtShbGfqcRoI7cenP47726qLn/XPvytBWY2q55zbLfjWnbxb4smZRhix6nL6r5eM/kdg3PO1ajPuvOOwSsG55yrQR/iDUa0uhgN4RWDc87VyO8YOkS5tttX39xRkyg1RBI3yNJ9ttg+yfG1pJOoValr5BV72H/EEWGhUlwhrcy+Hj8YHozuDT53XcXgnHPNIXqt5ZNgNoRXDM45V4Mwg5tXDM4551K8Kcl1hXZIzVEuNtDUmEU1MQXnCpiJ1ea9kpxzzkUh+NydPy68YnDOuZp48Nk551yKB58bSNK6wE+B7Qnv9ceBB4FLgGnAEuCDZvZ8ufPM2H4zrplfPEdN3rmR2l2xHEf1xBbK5UzKMp6hmvEGjcqL1I3pwl3r9XbpALd2qO7OAq41szcBOwEPACcDN5jZDOCG+Nw559qGIVbbyEyPSiSdK+kZSfeV2P5hSffEx58l7ZT7C0ppacUgaSKwJ/AzADN7w8xeAA4Dzou7nQe8txXlc865UpLgc5ZHBj8HDiqz/TFgLzPbETgdmFP3CyijpWm3Je1MeIGLCHcLC4DPAcvMbN3Ufs+b2XpFjp8NzAaYOnXqLo8//njZ69XbpHTLlV/KfJ52SMWdd5NSqXO22kGTZ5fcdu2K8P+nqpTa5cR024Wptl37KZf+XFLdabCn7zDe/v3y7TPt+7GZf6l4PUnTgKvMrOxJJa0H3GdmU7KWtVqtbkoaCbwV+JGZvQV4lSqajcxsjpnNMrNZG2ywQaPK6JxzRfXRk+kBTJY0P/Uo/WumsuOAa/J5BcW1Ovi8FFhqZn+Jz39FqBielrSJmS2XtAnwTMtK6JxzRZhRTXfVFXlM1CPpHwgVwzvrPVc5Lb1jMLO/AU9K2iau2pfQrDQXODquOxq4sgXFc865MkRfxkcuV5N2JPTgPMzM/p7LSUto9R0DwGeACyWtBTwKHEuosC6VdBzwBJBLX8NqYgTJvtVuS869cr8dgNbGGhqR/iI5ZzvEGpK4wYhJQ8JPQ6S7wfan2S7QM2b0kHV9K1cW3XfNI0v6l0duNa3i9V1l5WICWVKZZ5lSNU8GvJGhx1EeJE0FLgc+amYPNfp6La8YzOxuoNgt1r5NLopzzmVmKLeJeiRdBOxNiEUsBb4GjAIws3OArwLrA2dLAljTyDmkW14xOOdcp8orV5KZHVVh+yeAT+RysQy8YnDOuRoY0Oe5krpHuRhBvdph/EIztFOsoZxkjEPfS6/0r9OIkCp5xKR1AVjzbIjj9b2+qn+fYvEGgBHjxw9ZZ38Lnea08Yb1Fzild8mTJbeNmLZ5LtdYs/jRouuzjNP47S2/qbjPu/d8X9VlqlazYwsD5PMxOOecG+B3DM455wbxiXqcc84N0a3zMbQ0V1KeZs2aZfPnz291MYZI2uKTcQ3N1Mx4R7vEGkrlTUrHGBKFMYaiYm4kjRw1aHWpGERauZhDEj8oFysoF2MoVOw8hTGAAzfNJyHnvKcW5nKecvGHRscN8siVtOl269rsS/bKtO/Xd5hb9/Waye8YnHOuJj6Dm3POuZQQfPZeSa4GSRNLK2aRK9Z8Vdi8VK6Jq5u63vZMTHUz7QvNp/bG6vA8NhcNkqTmLpGiu/eVgaapYl1Ys+pbtnzIOlu9purzZOk6WqwJqJrmpYY1Ia3pBeCaR7+Ty/mbJUzU48Fn55xzKT7ns3POuX4h7bY3JTnnnEvxGIPrCvV2m22XbqkAB008trYDRw7+2I9Yd10Ael94YWBlQdxhxKT1Bz2/5m9n9y8fvPGnwsLECUAqzcTioV1kk1QT1ts76N9iCrvEplN25CWvuEElRbumrhn82g/e8ov9y50QbwjZVb0pyTnnXMqwzJUkaW6GczxnZsfkUxznnOsMhljTNzx7Jb2Z8jnABfyw3kJIGgHMB5aZ2T9KmgRcAkwDlgAfNLPn672Oc87lKa9pO9tNpYrhVDO7udwOkr6eQzk+BzwATIzPTwZuMLMzJZ0cn1ee28/lopFpydtZEmuAgnhDEf1xhbSXXq54jfSUoJWUiynMW3ZX5vO00sHTdgVg3utDU1ykYwrQGXGFtG7ulVQ2cmJml1Y6QZZ9ypG0GfBuwiTXicOA8+LyecB767mGc841Qp/1ZHp0mkzBZ0kzgS8BW6SPMbN9cijDfwNfBiak1m1kZsvjNZZLKpqNTNJsYDbA1KlTcyiKc85lk+ecz+0ma6+ky4BzgJ8ApfvXVUnSPwLPmNkCSXtXe7yZzQHmQMiumle5nHMui+EaY0isMbMfNeD67wAOlXQIMAaYKOkXwNOSNol3C5sAzzTg2q5Ap8UWel99DYAR49auuK/1DYxLSP4r970cYgI9EyYMPaBEjqQsRm64wZB1ham9R26w/pB9Kh3TCZKYgkaPjv8W2acgttCpDLq2V1LZT7+kSbGH0P9J+pSkTZJ1cX1dzOwUM9vMzKYBRwK/N7OPAHOBo+NuRwNX1nst55zLlYWmpCyPTlPpjmEBoWJMXln6J6UBlWcMr82ZwKWSjgOeAD7QoOs451xNjPyakiSdCyRN69sX2S7gLOAQ4DXgGDO7M5eLFytPlhncJI0xs9crrWuldp3BrZhWpOAupdOaj4pJZskrSqX/4/atDB/f63ovKbnP/iOOGPQ82Xf/nqG/VZImpN7nXgiXXmtg1ree8eNKl7GEpCmp2Gxxtt3WAFw79xdVn7eRejZ+CICD1vk4AFq/csNCK7qp5jGD23pv2tD2/lm236xXvPPssteTtCfwCnB+iYrhEOAzhIphd+AsM9u9poJnkLUh9c8Z1znn3LCQTNSTR1OSmd0CPFdml8MIlYaZ2W3AujH+2hCVUmJsDEwBxkp6CwNNShOByhE/55zrUiElRuZOCpMlpZs05sRelVlNAdKTgC+N64bO9JSDSjGGA4FjgM2A7zJQMbwEfKURBXLOuU5RRYxhRZ1NV8Uu1LAu+mUrBjM7T9IFwFFmdmGjCjHcJO36tXTbe3WHjXItQycrG1uoQmEcoZyDN/ssACM3HXoX3/f8CwBoVPxvlSF+V06x2EKhgw79CNB+sYZrXzwXKP4Z77TUFyVZU+djWApsnnq+GfBUoy5W8T7IzPqATzaqAM4514nyjDFkMBf4mIK3AS8m2SEaIesAt+skfZGQ8fTVZKWZlQuWOOdcV8vrjkHSRcDehFjEUuBrwCgAMzsHuJrQI2kxobtqjbNUZZO1Yvh4/PfE1LpGjmNwzrm2lmeuJDM7qsJ2Y/D3b0NlqhjMbHqjCzIcJW2t1cQaxt37NFB/rCHph39d32V1nacVMsUWyoxfaIW+lwdP89kzYXzFY5J9kmOTsQtpv7oiJCU+ePpeAFzz2F/qKmfeuiaeUEJvB2ZOzSJrdtVRwAnAnnHVTcCPzWx1g8rlnHNtzZobfG6qrE1JPyK0dyUzoH80ris3u5tzznU1G+YVw65mtlPq+e8lLWxEgZxzrjN0ZoK8LLJWDL2StjKzRwAkbUmO8zK46iWxBsjWjpvkZxr9f7c3rEyNVjK20GbxhCySuEGWWEO52ELiskfCDLyvDGQX5wNb7TX4mqtCbqhfLx2IQ4zvGZOtwGUkqbYB5rVN9rTmGO53DF8CbpT0KGEE3hY0uLuUc861s2QcQzfK2ivpBkkzgG0IFcNfzaz0TOXOOdftDHqHc8UQ7QJMi8fsJAkzO78hpWqiwhTYrUgVUawpqJZ0GYXH9D65tH+5v8FgRPfMONW3cuWQdT1rty63Y8966wIDqTGySHdjTV5PYRfiA3b7d2Bo81ElSfPSEfuFtBkDHSvbq0trpzKGeVNSzJe0FXA3A7EFAzq+YnDOudp48HkWsK1lmdWnCpI2J1QuGwN9hFS0Z8VpQy8h3KEsAT5oZs/neW3nnKtXvt+I7SPrsL37CF/eeVsD/IuZvRl4G3CipG2Bk4EbzGwGcEN87pxzbcVMmR6dJuvUnjcCOwO3A/1BZzM7NNfCSFcCP4iPvc1seZyl6CYz26bcsfVO7dkOsYZiaok1rFnyOAAjp23Rvy6JYxww6siSx1mfJQuD1rdL2oxi02mWopFhWs309JpDrFlTclPf6uLbesZm6N5Z48/IUjGGxEtPTa3pvIkk1mDL/lZyHxWZgvS3C+YNep5M35k4ePPP9S9f8+RZ9RSxafKY2nPs1pvalt+dnWnfRe/9et3Xa6asTUmnNbIQAJKmAW8hRMY2SlLKxsphwxLHzAZmA0ydWt9/Guecq1a3NiVl7a56c7ntkm41sz1qLYSk8cCvgc+b2UvKOGApTo03B8IdQ63Xd865WnRiM1EWeaUGrHn4ZEzQ92vgQjO7PK5+OpnoOv77TP1FdM65/BjZ4gudWHlUM46hnJp+rSvcGvwMeMDM/iu1aS5wNHBm/PfKuktYQbvEFApVk5q7d+ngmf7WPD4wd3gyfeV1vZeE50XSblczxWXHKBNHKEc9g/8za3TlaTb7VoZ8EFmm5CTGc5I0FVA5ljNx0yf6l7PEG9L7A1yzKPx70DofH7Jv7ytx/q3478iNB1pv373LgeH4Zf8z6Jh0bGG46tZmirwqhlq9g5Cp9V5Jd8d1XyFUCJdKOg54AsgedXTOuWaw7m1KyqtiqOndMbM/ljl239qL45xzjWd93VkxZIoxSBonqScuz5R0aIwNJD7akNI551wbM8v26DRZ7xhuAd4laT3CgLP5wBHAhwHM7L7GFM8lsqTWPmCtDwEDffitt3Rm9KQ9O93e3LPWWgDMW3lBzeVsJhXL+6TBv3WSsRnaYUbpE92/OOxb5v0qZKuG5pAcOWUTANYsXTZkW8/owf0z+t54A8h/jEhhXKGYa188t385+fuPSMYvjMreiJCMWUjHGpLlThnPUI9uzpWUtVeSzOw14J+A/zGz9wHbNq5YzjnX5gwwZXt0mMwVg6Q9CHcIv43rWh24ds65lhruTUmfB04BfmNm98cZ3G5sWKlcTX73xi8HPS/W/bQwrURhEwfAgWMHh4yq6VLZFLG5qL+ZqCefX2TppqmkqadfkaajaqTfw3Cx+oYQJU1GSbfVw6e9M275UP8+hZ+HLLJ0y02U6646PJqUlFvwWdJBwFnACOCnZnZmwfZ1gF8AUwnf298xs//N5eJFVDPy+WZJ4+LzR4HPNqpQzjnXEXK4G5A0AvghsD+wFLhD0lwzW5Ta7URgkZm9R9IGwIOSLjSzN4qcsm5ZeyXtIWkR8EB8vpOksxtRIOec6wiWW3bV3YDFZvZo/KK/GDhs6NWYEAcFjweeI2SnLkrSkNu5YutKyXo/+9/AgcDfAcxsIbBn1os451xXsowPmCxpfuqRTss6BXgy9XxpXJf2A+DNwFPAvcDnzArSIA92dJF1x2R7UVUEkM3syYLkdtn79rmWSNJfQHUpq8spPE8zYw7JtQpjJ+lupj1rFXRhLdeWH7upDoknpA+vYSrUkZuF/9PFuq0m0n+blordU+3VkPLbnnsBgJ5J6/bvUpgKozAG1TN5UuPK1/YyxxhWlEm7XewkhY1UBxJm0NyHMJvmdZL+YGYvDTqRdBQh0DRd0tzUponEH/ZZZK0YnpT0dsAkrUWILzyQ9SLOOdeV8ulxtBTYPPV8M8KdQdqxwJlxFs3Fkh4D3kSYIyftz8ByYDLw3dT6l4F7shYoa1PS8YTgxxTCi9g5PnfOueHJgD5le5R3BzBD0vT4w/tIQiLRtCeIaYIkbQRsAzw6pEhmj5vZTcB+wB9ix6HlhMom8+1N1l5JK4ijnJ1zzgV5jFEwszWSPg3MI3RXPTcOCzg+bj8HOB34uaR7CV/wJ8Xv5VLKZquoJFPFIGkm8CPCzGrbS9oRONTM/iPL8a71CmMBxVIvJ+yN1RXPN2LSeuE8k0MMzaaHdvV5d3yt1iJmVpg6vGoxtpBFEr+oJ9aQtmbZ8qrPU85ACowwfiHL2IV99/nPgSc7Tgdg5K2LBu1TGFeAgdhCMiYjGQOTHqswPMYvpOQ0eM3MrgauLlh3Tmr5KeCAKk4pM3stZqj+HzP7lqS7sh6ctSnpJ4QBbqtjIe8h3O4459zw1b4pMerKVpF1x7XN7PaCXkm1zYDinHNdQu2b7uLz1JGtImvFsELSVsQbJ0mHEwIazjk3PA2MUWg79WaryFoxnAjMAd4kaRnwGA0ORlfKHeJqk6TmTvSMrTxdd88u2/cv9y1onwzr9Y6hyBKjGLnhBoOe22srS+5rLw7qUj4wXWZKo8YvVJMX6Ybfn9K/PCjeUMKQcQsxtrBmj22HnmPGxoPWpa/VfTL1OGqJ2Iz0M8Io6amSdgI+aWafynJ8xYoh5vE4wcz2i7VPj5m9XE+hM16zUu4Q55xrrTa9Y2AgW8VcCNkqJGXOVlEx+GxmvcAucfnVRlcKUZbcIc4511rZU2I0v2hmTxasypytImtT0l1xePVlQP/9sZldnvVCVSqWO2T3wp1ivpHZAFOnTm1QUYah7bcuuSlpVir1WT9w16+XPLYZXVmzqKUJKdG3cmhTUs/YsRXP1zYpMAokTT2FTUqDnu9R+5xcWZqqCsuSlqQ/acv3L5mopz3Vla0ia8UwiZBnY5/UOgMaVTFkyR2Cmc0hxD6YNWtW+97UOee6Uhv3SjqeEKNNslX8jiqyVWQd+XxsTUWrXZbcIc4511ptWjFUylYh6RQzK3k7l3Xk8/eLrH4RmG9mV2Y5R5X6c4cAywiD6T5U/hDnnGuuNr5jqOQDQH0VAzCGkMkv6R/4fuB+4DhJ/2Bmn6+nhIVK5Q7J8xrDTWGXw0Q6/YXWGtXYMhSJP7Qi7pB0cy31nlQtSZfRNzg9flu2i7exYlPRFm5ru/e0fWMMlZQteNaKYWtgHzNbAyDpR4Q2q/0Jk0bkrljuEOecaxttPMAtg7Ilz5oraQowLvV8HLBp7Mpa3yzpzjnXqdq4u2oFudwxfAu4W9JN8YR7At+IA96ur6t4zjnXoTo4xlA2bUDWXkk/k3Q1YeCZgK/ENLAAX6qvfK6ZyqWPtplbNLEkrTdv5QVA8VhD34uDx3H2rDMhLJSbKrQn6w14+ykcQ1DN+IO0npvuHPQ8iefUer5CxcagNHN62SHKzbrcQjFp3lnAHoRS3gp8IeZMwsy+Ue74TJ9khbSq+wI7mdkVwEhJu9VRbuec62iy7I8W+CVwKbAxsCnhDuGirAdn/YlzNqHmOSo+f5mQy8g554avNp6PwcwuMLM18fELqoh2ZI0x7G5mb01mADKz5+Mwa+ecG77aN8Zwo6STCXnmjDCt528lTQIws+fKHZy1YlgdM54m8zFsQNu2rrliyrWn97tv8JSXtiaMcdDOb25YuTpFYcwhzV7vvo555dJlVxMvOHBMGHxb7ItmzdtD3q2Rfw6p3Mt9oRSOX0jHGpLlVsQa2jj4nAwK+WTB+o8Tvse3LHdw1orh+8BvgA0lnQEcDvxrFYV0zrnu06YVg5lNr+f4rL2SLpS0gBCAFvBeM8ucqc8557qOgdq03UTSKOAEwtACgJuAH5vZ6pIHpZStGJL2qOgZUlFtSZMqtVO59pOpSamA3T3wGyCPZqV2Sb9dj3LdfotJmlQS816/MM/iNFW5ZqbC15nlPMkxPaMGvo4qvT/Fmo1a0qTUpncMwI+AUYSOQwAfjes+keXgSncMCwgvXcBU4Pm4vC7wBFDX7YpzznWyvGIMWaYylrQ3YWa2UcAKM9urzCl3NbOdUs9/L2lh1vKU7a5qZtPNbEtCMrv3mNlkM1sf+EcaNxeDc84NG6mpjA8GtgWOkrRtwT7rEn79H2pm2xGyo5bTK2mr1PFb0oAZ3HY1s+OTJ2Z2jaTTs17EOee6Uj53DP1TGQNISqYyTs9x/yHgcjN7AsDMnqlwzi8Ruqw+Gp9PAzLPq5O1Ylgh6V+BZJDERwgzurWNhxY8yv49H2jt8PhhIIk3dEMX1qriLL2Zf2wNnL+K9vZOVuvrbMX7k2Va18zyCz5nmcp4JjAq5qubAJxlZueXOeefgB8TOgwRl2/NWqCsFcNRwNcIXVYNuIWBUdDOOTc8Zb9jmCxpfur5nDg1MWSbyngksAvhi34scKuk28zsoRLXOx94CUhado4CLqByE1T/xSqKvY8+l2Vf55wbDkRVwecVZjarxLYsUxkvjed4FXhV0i3ATkCpimGbguDzjbkFnyWdVukEWfYpcdy3Jf1V0j2SfhODK8m2UyQtlvSgpANrOb9zzjVcPvMx9E9lHFMNHQnMLdjnSuBdkkZKWpvQ1FRuLNldkt6WPJG0O6F5KZNKdwyfkPRSme0ivIjTsl4w5TrglDiN5zeBU4CTYjT+SGA7QlbA6yXNjJMClTRzly25br7HF7JKxjOkJW3uSSqMhPUN/WTXkhas3cYv9K16veZjy41j6B8r0qUxhma8ruQaWcZ75Bo3qEZOmVNLTWUs6fi4/Rwze0DStcA9hOwhPzWz+8qcdnfgY5KeiM+nAg9Iujec0nYsV6ZKFcNPCIGOSvtUzcx+l3p6GyHNBoRo/MVmtgp4TNJiQtQ+c+DEOeeaIqdxDMWmMjazcwqefxv4dsZTHlRPecpWDGY2dPb2xvg4kGTJmkKoKBJL47ohJM0GZgNMnTq1keVzzrkh2jUlhpk9Xs/xWXsl1UTS9YSJIgqdamZXxn1OBdYAyT1jlgh9WBmi+nMAZs2a1b6D051z3alLv3UaWjGY2X7ltks6mjCKel8zS97iLBF61wDV9NXvuzOMvel567YV9mxfyZiXvNqoi8VtukkrYibVxBqySP7mYVLKOmULLHeklk1SG3ODnEQY4v1aatNc4EhJoyVNB2YAt7eijM45V04bT+1Zl6xzPs+UdIOk++LzHeNI6Hr8gBDYvk7S3ZLOATCz+wlzlS4CrgVOrNQjyTnnWiKf7qptJ2tT0k8IuTd+DGBm90j6JfAftV7YzLYus+0M4Ixaz+3qpPB7oWfMaAB6X32t3N4ltVv31DyMWGcdAHpffLF/Xc/oMZmP7+R0261UrkmplWlwOvFuIIusFcPaZnZ7QbvcmgaUxznnOoPRtRMcV5NEbysG5nw+HFjesFI551ybE7UN9OwEWSuGEwndQt8kaRnwGCHDqnPODV/DuSkp5gnfT9I4oMfMXm5ssVw7GTFu7f7lSvGGTowrlOu2WhhTSMcWEqW6qXo8ofsNyxiDpP9XYj0AZvZfDSiTc851huFYMTCQJ2kbYFcGMv69hzAng3PODU/5TdTTdjLlSpL0O+CtSRNSTLXtqUydc8PbML1jSEwF3kg9f4Mwh6jrQsn4haLbOjgFRiVJPKHctiTGUG7fbpXETLo1nXgthmWMIeUC4HZJydSe7yNMHeecc8PXcK4YzOwMSdcA74qrjjWzuxpXLOeca3/D+o5B0lRgBfCb9Doze6L0Uc4518U6NA9SFlmbkn7LwFswFpgOPEiYftN1sAPW+lBV+/ctGDybYM8u2wOdOX6h0LXP/7R/+aD1PlF0nyS2kN7XDVU4hqPeuEQtY0KKjUvJM6+SGKa9khJmtkP6uaS3Ap9sSImcc65TDPM7hkHM7E5Ju+ZdGOec6ySy7qwZssYY0iOge4C3As82pESurY0YP37Q82u7oAmpGG8qKq2WZp1ix7Si22tes/UBHmNgYAQ0hHTbvwV+nX9xnHOuc3Rrr6SsU3suMrOvx8cZZnYhIS1G3SR9UZJJmpxad4qkxZIelHRgHtdxzrm8qS/bo+J5pIPi991iSSeX2W9XSb1x6oOGyVoxnJJxXVUkbQ7sDzyRWrctcCShx9NBwNmSRtR7Leecy10OU3vG77cfAgcD2wJHxe/BYvt9E5iXW/lLqJRd9WDgEGCKpO+nNk0knxncvgd8Gbgyte4w4GIzWwU8JmkxsBtwaw7Xc3W69qX/bXURXJfI0qW1rVOXW25NSbsBi+P0Bki6mPA9uKhgv88QmvAb3vGn0h3DU8B84HVgQeoxF6iriUfSocAyM1tYsGkK8GTq+dK4rtg5ZkuaL2n+s896LNw512Q53DGQ4TtP0hRCKqJzcih1RZWyqy4EFkq60MyqvkOQdD2wcZFNpwJfAQ4odlixopQo3xzCzHLMmjWrS8NAzrl2JKq6Y5gsaX7q+Zz4/ZWcqlDhmf8bOMnMepP5cBqpUlPSpWb2QeAuaehbYGY7ljvezPYrcd4dCKOnF8YXuRlwp6TdCLXl5qndNyPcuTjnXHvJPo5hhZnNKrEty3feLODi+H05GThE0hozuyJ7YbOr1F31c/Hff8zzomZ2L7Bh8lzSEmCWma2QNBf4paT/AjYFZgC353l9Vx2PK7hOkGWMQpISI5df3flN1HMHMEPSdGAZofPNoFw1ZjY9WZb0c+CqRlUKUCHGYGbL4+KnzOzx9AP4VCMKZGb3A5cSAi/XAieaWW8jruWcc/XIo7tqbKb/NKG30QPApWZ2v6TjJR3f+FcxVNYBbvsDJxWsO7jIupqY2bSC52cAZ+Rxbueca5icIptmdjVwdcG6ooFmMzsmn6uWVinGcALhzmBLSfekNk0A/tTIgjnnXLvr1pHPle4YfglcA/wnkB6N97KZPdewUrmGqzbdtnON1sgxC3mm2+5nVBN87iiVuqu+CLwIHAUgaUNgDDBe0nifqMc5N5x16x1DppQYkt4j6WHgMeBmYAnhTsI554alZKKePHIltZuswef/AN4GXG9mb5H0D8S7CNc9esaMbnURXBGNnomsWzT9PTHr2qakrEn0VpvZ34EeST1mdiOwc+OK5Zxz7U+W7dFpst4xvCBpPHALcKGkZ8gniZ5zznWuDvzSzyLrHcNhwErgC4RBZ4+Q03wMzjnXqYb1HYOZvZp6el6DyuJarO/1VYOe/+6NXwLext1q/l63KQN6O/BbP4NKA9xepvjNkgAzs4kNKZVzznWATrwbyKLSOIYJ5bY759yw1qW9krIGn51zzhUYlncMrnsl8YNyqTEOmngs4G3czhWVbXa2juQVg3PO1SDM4NadNYNXDM45VyMNx15JzjnnSvCmpMaQ9BnCzEVrgN+a2Zfj+lOA44Be4LNmNq91pexuSazBOVet7s2V1LKKISbiOwzY0cxWxZTeSNqWMOfpdoQ5n6+XNNOn93TOtZtu7ZWUNSVGI5wAnGlmqwDM7Jm4/jDgYjNbZWaPAYuB3VpURuecKy3JsFrp0WFaWTHMBN4l6S+Sbpa0a1w/BXgytd/SuG4ISbMlzZc0/9lnn21wcZ1zLsV8PoaaSLoe2LjIplPjtdcjzPOwK3CppC0JvcAKFa1yzWwOMAdg1qxZnVctO+c6W193fu00tGIws/1KbZN0AnC5mRlwu6Q+YDLhDmHz1K6bAU81spzOOVeLbh3H0MqmpCuAfQAkzQTWAlYAc4EjJY2WNB2YAdzeqkI651xJXRpjaGV31XOBcyXdB7wBHB3vHu6XdCmwiNCN9UTvkdQaxdJtJzxNRn3aNZV5Uq5iZSksczXlbdfXWxcDcoofSDoIOAsYAfzUzM4s2P5h4KT49BXgBDNbmM/Vh2pZxWBmbwAfKbHtDOCM5pbIOeeyE5ZLU5KkEcAPgf0JTel3SJprZotSuz0G7GVmz0s6mBBb3b3ui5fgI5+dc65WfbncMuwGLDazRwEkXUzott9fMZjZn1P730aIvTZMK2MMzjnXuZKmpCwPmJx0rY+P2akzZe6iHx0HXJPTqyjK7xhcv3IxhUTHtws3QGG7fJb3sdx5sij8O9R6zXrKUu81azm+3T5/VTQlrTCzWaVOU2Rd0RPHjBHHAe/MeuFaeMXgnHO1yqfHUaYu+pJ2BH4KHGxmf8/jwqV4U5JzztUkY1fVypXHHcAMSdMlrUXIFTc3vYOkqcDlwEfN7KGGvJwUv2NwzrlaGLncMZjZGkmfBuYRuquea2b3Szo+bj8H+CqwPnC2JIA1ZZqm6uYVQ4fIe0xB3m3SrdBur6GZ5Wm3194stb7uRsUm8pqox8yuBq4uWHdOavkTwCdyuVgGXjE451ytOnBUcxZeMTjnXC0MT6LnhmqX2/l2KUejDZfX6TpFZ+ZBysIrBuecq5VXDM455wbxisE551w/M+jtzsTPXjFE3n5dWq3plUsd5++1a6Smps3wOwbnnHP9vFeSc865Ibr0jqFluZIk7SzpNkl3xzS0u6W2nSJpsaQHJR3YqjI651xZPrVn7r4FfN3MrpF0SHy+t6RtCUmktgM2Ba6XNDPP6T29jbs69aaRbrdUyc7loouDz63MrmrAxLi8DgNpZg8DLjazVWb2GLCYMMORc861F79jyN3ngXmSvkOooN4e108hTF2XKDmbUZwFaTbA1KlTG1ZQ55wrqgO/9LNoaMUg6Xpg4yKbTgX2Bb5gZr+W9EHgZ8B+VDGbkZnNIUyKzaxZs7rzL+Sca1PmvZJqYWb7ldom6Xzgc/HpZYSZiSDjbEZZFLaNJ23dxdq8Pe7QOP7eumYq9nlrSJzLwKwv//O2gVbGGJ4C9orL+wAPx+W5wJGSRkuaDswAbm9B+Zxzrrw+y/boMK2MMfwzcJakkcDrxFhBnLnoUmARsAY4Mc8eSc45l4su7pUk65LgyURNst21b6uL4ZzrANfbrxbUOzXmOiMm2x7j3pNp33kv/7zu6zWTj3x2zrkaWV93xhi8YnDOuZp05hiFLLxicM65WngSPeecc2kGWJcGn71icM65WphBl45j8IrBOedqZN6U5JxzbpAuvWPomnEMkp4FHm/Q6ScDKxp07lp5mbLxMmXXjuVqVJm2MLMN6jmBpGsJ5ctihZkdVM/1mqlrKoZGkjS/3QaneJmy8TJl147lascyDQetzJXknHOuDXnF4JxzbhCvGLKZ0+oCFOFlysbLlF07lqsdy9T1PMbgnHNuEL9jcM45N4hXDM455wbxiiFF0gck3S+pT9Ks1PppklZKujs+zklt20XSvZIWS/q+pGJzVudeprjtlHjdByUd2KwyFZThNEnLUu/NIZXK1yySDorXXizp5GZfP1WOJfHvcbek+XHdJEnXSXo4/rteg8twrqRnJN2XWleyDM3425UoU9t+noYVM/NHfABvBrYBbgJmpdZPA+4rccztwB6AgGuAg5tUpm2BhcBoYDrwCDCiGWUqKN9pwBeLrC9Zvib9LUfEa24JrBXLsm2LPldLgMkF674FnByXTwa+2eAy7Am8Nf05LlWGZv3tSpSpLT9Pw+3hdwwpZvaAmT2YdX9JmwATzexWC5/e84H3NqlMhwEXm9kqM3sMWAzs1owyZVS0fE28/m7AYjN71MzeAC6OZWoXhwHnxeXzaPDfyMxuAZ7LWIam/O1KlKmUVn+ehhWvGLKbLukuSTdLeldcNwVYmtpnaVzXDFOAJ4tcuxVl+rSke2LTQNIcUap8zdLq66cZ8DtJCyTNjus2MrPlAPHfDVtQrlJlaPV7146fp2Fl2CXRk3Q9sHGRTaea2ZUlDlsOTDWzv0vaBbhC0naEpppCVff/rbFMpa6dS5kGXahM+YAfAafHa5wOfBf4eCPKUaVWXz/tHWb2lKQNgesk/bVF5ciqle9du36ehpVhVzGY2X41HLMKWBWXF0h6BJhJ+NWyWWrXzYCnmlGmeO3Ni1w7lzKlZS2fpJ8AV1UoX7O0+vr9zOyp+O8zkn5DaAJ5WtImZrY8Nv8904KilSpDy947M3s6WW6zz9Ow4k1JGUjaQNKIuLwlMAN4NN5+vyzpbbHnz8eAUr/w8zYXOFLSaEnTY5lub3aZ4hdK4n1A0sOkaPkaVY4i7gBmSJouaS3gyFimppI0TtKEZBk4gPAezQWOjrsdTfM+N2mlytCyv10bf56Gl1ZHv9vpQfggLiXcHTwNzIvr3w/cT+gVcSfwntQxswgf3keAHxBHkze6THHbqfG6D5LqedToMhWU7wLgXuAewn/eTSqVr4l/z0OAh2IZTm3RZ2rL+LlZGD9Dp8b16wM3AA/Hfyc1uBwXEZpEV8fP03HlytCMv12JMrXt52k4PTwlhnPOuUG8Kck559wgXjE455wbxCsG55xzg3jF4JxzbhCvGJxzzg3iFYNzzrlBvGIYRiS90oBzHpqktJb0Xknb1nCOmwpTimfY/0FJhxbZNi2dxrnbSfpKanlsTFX9hqTJrSyX62xeMbi6mNlcMzszPn0vIT1yM3zYzBo6kjkZ7d7m+isGM1tpZjvjqSJcnbxiGIYUfFvSfXECmSPi+r3jr/FfSfqrpAtjWg0kHRLX/VFh8p+r4vpjJP1A0tuBQ4Fvx1+tW6XvBCRNlrQkLo+VdHHMoHkJMDZVtgMk3SrpTkmXSRqf4fXsImmhpFuBE1PrR8TXeUe81ifj+h5JZytMgHSVpKslHR63LZH0VUl/BD5QqjzxmjfHjKnzklQOkj4raVG83sVlyjxOIXvoHQpZew+L66dJ+kO83p3xfUXSJpJuie/tfZLeJelMILlLuDDTH9+5LFo99NofzXsAr8R/3w9cR5jMZiPgCWATYG/gRUKCsh7gVuCdwBhCyuPp8fiLgKvi8jHAD+Lyz4HDU9e7iTi5EDAZWBKX/x9wblzeEVhDSOMxGbgFGBe3nQR8tcjr6D9vfH4PsFdc/jZx4hdgNvCvcXk0MJ8wycvhwNXxNW4MPJ+UmzCpzpdTZR5SHmAU8Gdgg7j+iNTreQoYHZfXLfO3+AbwkWQ/QuqOccDawJi4fgYwPy7/CwPpNEYAE9J/04JzL6FgYiB/+KOax7DLruqA8GV/kZn1EjJs3gzsCrxESMS3FEDS3YTZ614hJA18LB5/EeFLt1Z7At8HMLN7JN0T17+N0BT1p3ijshahcipJ0jqEL+Cb46oLgIPj8gHAjsndALAO4cv2ncBlZtYH/E3SjQWnvaRCebYBtiek0IbwRb08HnMPcKGkK4AryhT9AOBQSV+Mz8cAUwkVyw8k7Qz0ErL4QkgKeK6kUcAVZnZ3mXM7VxevGIancnNAr0ot9xI+I7XOGb2GgebKMQXbiiXpEnCdmR1VxTVU4lzJts+Y2bxBK6V3Vzjnq+XKI2kH4H4z26PIse8mVHyHAv8maTszW1OibO+3gtn5JJ1GSJa4E+G9ex3CbGeS9oznv0DSt83s/Aqvw7maeIxheLoFOCK2wW9A+CIrl8L4r8CWkqbF50eU2O9lYELq+RJgl7h8eGr9LcCHASRtT2hOArgNeIekreO2tSXNpAwzewF4UdI746oPpzbPA06Iv7KRNFMh9fUfgffHWMNGhCa0YkqV50FgA0l7xPWjJG0nqQfY3MxuBL5MaCIqFSOZB3wmFcN5S1y/DrA83s18lHA3gqQtgGfM7CfAzwhzJQOsTl6fc3nximF4+g2hyWMh8HtCm/rfSu1sZiuBTwHXxqDs04RYRKGLgS/FYOpWwHcIX8x/JrTXJ34EjI9NSF8mVkpm9iwhZnFR3HYb8KYMr+dY4Icx+Lwytf6nwCLgToUurD8m3AH9mpDmOVn3l2Kvp1R5LMwhfTjwTUkLgbuBtxO+xH8h6V7gLuB7seIq5nRCrOKeWLbT4/qzgaMl3UZoRkruXvYG7pZ0FyFGdFZcPyeew4PPLjeedttlImm8mb0Sf+H+EHjYzL7XorLcBHzRzObXcY7k9axPqJjeUa5y7CSx99csM1vR6rK4zuR3DC6rf47B6PsJzR0/bmFZngN+riID3KpwVXw9fwBO74ZKQXGAG+FOpK/FxXEdzO8YnGsgSccCnytY/SczO7HY/s61A68YnHPODeJNSc455wbxisE559wgXjE455wbxCsG55xzg/x/EJ0pfad/K+EAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "geo_grid = make_geocube(\n", " vector_data=gdf,\n", " measurements=['pop_est'],\n", " resolution=(-3.0, 3.0),\n", " rasterize_function=partial(rasterize_image, all_touched=True),\n", ")\n", "geo_grid[\"pop_est\"].plot()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.13 ('main')", "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" }, "orig_nbformat": 4, "vscode": { "interpreter": { "hash": "f7286e25aede83087e9621ce0e75f07eda2aa917ee0205c1a0b96e7ba4f20042" } } }, "nbformat": 4, "nbformat_minor": 2 }