Lillian-Violet/stembureau-meting
1
0
Fork 0
mirror of https://github.com/openstate/stembureau-meting.git synced 2024-11-14 02:34:02 +01:00
Code Issues Packages Projects Releases Wiki Activity
stembureau-meting/jupyter/.ipynb_checkpoints/stembureau_data-checkpoint.ipynb

2090 lines
323 KiB
Plaintext
Raw Normal View History

Filtered out municipalities without elections
2022-04-13 15:08:29 +02:00
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Data gemeenteraadsverkiezingen 2022 Nederland\n",
"\n",
"Dit notebook is voor het verwerken van de data van de gemeenteraadsverkiezingen van 2022. Er zal hier stap voor stap door de data gelopen worden om het proces reproduceerbaar te maken voor latere verkiezingen. De eerste stap was de data ophalen van de bronnen, zowel de overheid als waar is mijn stemlokaal (voor geografische data van de stemlokalen). De bronnen gebruikt voor de data zijn voor de verkiezingen van 2022 is als volgt:\n",
"\n",
"- [Verkiezingsuitslagen Gemeenteraad 2022](https://data.overheid.nl/dataset/08b04bec-3332-4c76-bb0c-68bfaeb5df43)\n",
" - [Directe link naar uitslagen per gemeente CSV](https://data.overheid.nl/sites/default/files/dataset/08b04bec-3332-4c76-bb0c-68bfaeb5df43/resources/GR2022_2022-03-29T15.14.zip)\n",
" - [Directe link naar kandidatenlijst met uitslagen CSV](https://data.overheid.nl/sites/default/files/dataset/08b04bec-3332-4c76-bb0c-68bfaeb5df43/resources/GR2022_alle-kandidaten_2022-02-22T08.34.csv)\n",
"- [Waar is mijn stemlokaal stembureau data](https://waarismijnstemlokaal.nl/data)\n",
" - [Directe link naar waar is mijn stemlokaal gemeenteraad 2022 CSV (CKAN)](https://ckan.dataplatform.nl/datastore/dump/d6a1b4c4-73c8-457b-9b75-a38428bded68)\n",
" - [Verkiezingsuitslagen gemeenteraadsverkiezingen 2022 geodata (Volkskrant)](https://data.openstate.eu/dataset/verkiezingsuitslagen-gemeenteraadsverkiezingen-2022)\n",
" - [Directe link naar GEOJSON bestand](https://data.openstate.eu/dataset/a1767f1b-bf0c-409b-b3b1-3af9954b57f4/resource/413be255-5070-48f4-b631-895097976abb/download/2022gr.geo.json)\n",
"- [CBS Wijk- en buurtkaart 2021](https://www.cbs.nl/nl-nl/dossier/nederland-regionaal/geografische-data/wijk-en-buurtkaart-2021)\n",
" - [Directe link naar zip bestand met SHP bestand er in](https://www.cbs.nl/-/media/cbs/dossiers/nederland-regionaal/wijk-en-buurtstatistieken/wijkbuurtkaart_2021_v1.zip)\n",
"- [CBS bevolkingsdichtheid kaart 100 bij 100 meter](https://www.cbs.nl/nl-nl/dossier/nederland-regionaal/geografische-data/kaart-van-100-meter-bij-100-meter-met-statistieken)\n",
" - [Directe link naar 7z bestand met SHP bestand er in](https://www.cbs.nl/-/media/cbs/dossiers/nederland-regionaal/vierkanten/100/nl_vierkant_100meter_bij_100meter.7z)\n",
"- [CBS bevolkingsdichtheid kaart 500 bij 500 meter](https://www.cbs.nl/nl-nl/dossier/nederland-regionaal/geografische-data/kaart-van-500-meter-bij-500-meter-met-statistieken)\n",
" - [Directe link naar 7z bestand met SHP bestand er in](https://www.cbs.nl/-/media/cbs/dossiers/nederland-regionaal/vierkanten/500/2021-cbs_vk500_2020_v1.zip)\n",
"\n",
"De eerste stap die we moeten maken is de data importeren voor de analyse, daarna kunnen we kijken hoe goed de data is, hoe we het aan kunnen vullen, en wat er mee te doen. De makkelijkste structuur die we vonden was het geojson bestand van open state en de Volkskrant, daar staan alle stembureaus al in een lijst, en we hebben een makkelijk framework om het te importeren; geopandas. We laden deze dan ook als eerste in."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/lillian/.local/lib/python3.9/site-packages/fiona/collection.py:208: FeatureWarning: Empty field name at index 61\n",
" self._schema = self.session.get_schema()\n",
"/home/lillian/anaconda3/lib/python3.9/site-packages/geopandas/geodataframe.py:600: UserWarning: Empty field name at index 61\n",
" for feature in features_lst:\n"
]
}
],
"source": [
"import pandas as pd\n",
"import geopandas as gpd\n",
"\n",
"crs = {'init':'EPSG:4326'}\n",
"df_geojson = gpd.read_file(r'../data/2022gr.geo.json', crs=crs)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We importeren nu de gemeente shapefile kaart van het CBS, om de overlay van onze stemlokalen eroverheen te kunnen doen. "
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAOQAAAD4CAYAAAD8Q0ptAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB9nUlEQVR4nO29d5hk6V3f+3nPOXUqx845TM65Z1eriBAIAcKA4QpzsSUbC/yQrn0v2BgHHPDFxnCJtpEJNkmAAZkggYSEtCvtandyjj3TOXdXzqfOOfePU3Wmqququ2dnZqdnpz7Po0fT1V3Vp3rrd973/YXvV5imSYsWLbYH0tO+gBYtWjygFZAtWmwjWgHZosU2ohWQLVpsI1oB2aLFNkJ52hfQiPb2dnN4ePhpX0aLFk+M8+fPr5qm2bH+8W0ZkMPDw5w7d+5pX0aLFk8MIcRUo8dbW9YWLbYRrYBs0WIb0QrIFi22Ea2AbNFiG9EKyBYtthGtgGzRYhvRCsgWLbYRrYBs0WIb0QrIFm8LDMOkWDIwjMbzvYZh8itfHOfqbOItvrKHY1t26rRo8TC8Nr7Kv/v0TW4uJJEEDETczMXyyJJAkQSyJABI5kv80t/c5d98+ADfeXIAIcRTvvJ6thSQQohJIAXoQMk0zZNCiH8HfAtgAMvAR03TnN/Kcx/Ppbd43plYTfNTn77F528u2Y8dGwyTymv0hdxMRbN4vSrDbR6yxRLJxTR5zeCf/vFV/ubWMv/hWw/R5nM+xXdQz8OskO8zTXO16uufMU3zXwIIIX4Y+FfA92/xuS1aNKUiK9NoBSvpBjOxHC/fXuZzN5Y4Oxmt+b4soM3rZGI1hUMWFDWdmbUMmUKp5uc+e32J81MxfvrbDvO1+7ue3Jt5SN70ltU0zWTVl16gJc7T4rHwhZvL/Iv/fY0XRiMc7AuynCpwfyXN/dUM02sZIl4n3UEXbofEiaEwyVyJdr/KK3dWuTGfpFAy0MpnSU3XMYG8Xv/xXE0X+d7fOsfHXhrmx79hH6ry9FMqWw1IE/icEMIEftU0zU8ACCF+Cvi7QAJ438M8t0WL9ZimyVqmyB+emyGWKTAfzzOfyFMs6aTyJTKFEkIICiUDl0Mmni0SzxaJZjWcjiAA6aJe97qZBo9V85uvTnJuMsYvfOQoox2+J/LetorYiuqcEKLXNM15IUQn8NfAD5mm+UrV938ccJmm+a8f9rlVP/dx4OMAg4ODJ6amGk6ntNgm3FlK8a/+9BoCwXe/MMjX7e/e0gpTKOksJvLMxXK8dm+N+ytpLs7EGW7zcH4qTlE36A+7afOqXC5nRIUAWQi8ToXRdi/X5hJohsmp4TBnJ2OP7T35XQqf/IcvcLAv+NhesxlCiPON8ilbCsh1L/STQNo0zf9c9dgQ8GnTNA8+7HMbcfLkSbM1D7k9mVrL8OeX5/n5z9+lVN4WHu4PYhgm79rdQZtXJVPQyRZLpAslskVrdVtO5kgXSxRLJook0A2TmViOw/1BYpki7T4nc7Esy+kie7p8yJLA6ZBxyBIl3WAtXWQ1XahZ7U6PRHhjItrsUreEJKC6UtLuc/I73zvG3u7AI73uZjQLyE23rEIILyCZppkq//vrgH8rhNhlmubd8o99GLi11ec+yhtp8dazlMzz2euL/O+Lc1yYjgMwNhLhzESUE0Nhzk9Zq9R8Ikeb14nXqSAJKJQMMgUrMONZDSHgcH8IWRLoukEir+FSZGRJcHEmzuG+AMvpIpNrGVRFJpUvbXBVD5+0UGXBvt4AqiyRyGkspwrEsxqjHV7avCqGAflSiW//L68hhOB7Xhzkx75+71taHtnKGbIL+FT5ohTg90zT/CshxB8LIfZglT2mKGdYhRC9wK+ZpvmhZs99/G+jxeNmJprls9cX+czVBaKZIp1+J4ok2dvEi9MxBiNuLs/E7edEMxoRrxMBxLJFVtNFgm4Hox0+8kWdibWMHbw+VWZPT4AzVVlS1SEDUCiZFErNg1GR4EBvkHi2uKX3EnYr7OoOML6c4vJMfWPA/ZUM91cydAdc9IXchL0qhmny+v0on7+5zAfewizsQ29Z3wpaW9a3FtM0WUjkuTIb58psgq/eX0MrGbhVmZlYjrV0AVWW2NfjRzfg0kwcEzg5FEYzDDIFnUyhxEDEzc35JKmCjmDjFexgb4Br88max9q8KgNhD9OxLNFMbbAJ4EBvAJdD5vZiilShPmB9qkx30EXEa9UWS4bBUrLAfCLHw37M93T6WMkUeWlnO7/0Xcce7slb4E1vWVu8/VhJFezguzqX4NZCkoDbQdDtIJHTuL+aoVgyAFAkwWiHl7ymc24qXvM6s7Eci8m8/fVC4sG/m33+hyIeuoOump+tsJYpslYOxO6Ai91dPjSthJBldMPEME0uz8QplksYOzq8hD0OllNFllN50kWd8ZUMrGQe6u/hkAS6abKvJ0ChZKDKEpmCRjRT5JU7K+iGaXf7PGlaAfmcsZYu8NJP/w0hj4OhNg+Zgs5KusB8VYCMdngJuhxougEC0jmN6Wiu/sVE/ZceVUbTDTtoAHZ3+Qh5VOZiWaai1v8A+sNueoIuDBNyxRLRjEbJMPA5FVwOmdV0gXafSk4zuLOUJpHT6PA76Qm6cCoS56ZiD73yrSfscTAY8bCcKnC9asU+MRhmKpojkdO4Pp/gcH/o0X7RFmkF5HNErqjz85+/y5GBIGcnYyynCjXfVxWJAz0Bbi0mua9tvsqsXzP29fi5sZDCqUgc6PVzYz7JscEQF6bj9IfcdAVcdAVcCASz8SyzsRyzsfpAX01bq+S7d7WzksyRyOt0BZwkchorqQIr6677zeJSJIJuh11eafbm/sdrk/zcdx59LL9zM1oB+RxQ0g1+49UJPvHKffwuB7Mxa4VyOST2dQdQFYliSefafJKLVUmah0Uqb+sKJYO7y2l6w26WkwVODoW5OBNnNt5gld0A0wSfW6U9IDMXy6FIwi61vFkckkBVJFRFos3nZHw53fDnqhNGf3Jhjr8zNsjJ4cgj/e6t8PR7hVo8cW4tpvh/P3MLw4SQx4GmmxzpD6KVDC7OxHljIsrFmQRag/ayjaheITv9TgolncpRq92n0uV3MhvPcW4qhv4mAqmgG6iKhGSauByCTv+jNYKfHomAsDp3skW9aTACzMfzNe/vp//yFm9FArQVkM8B1+YSdAdduBwSl6bj9ARdXJ5N8JDxV0/5E9sdcCIJwd2lDC/uaGNfj5+FRN6uWb5Z5mI5UvkS99ay9AQ8dAZcvGNHBJfj4T+2xwZDnJ+K2TedQjlp1YycpnNsMISr3H10birGTKNz9GOmtWV9Drg+n6zJajbKcL4ZBIK+kAtNN2j3OekKOCnpJjcXUg/1OmGPA5dDJqfpxLOa/fh8IkdO04lmiqylixRKOof6goy2+7ixUFsyOdwXJJYrohsm8/H696fK0kNvdy9Mxzk5FObOUopkvsSfX5nnw0d6kSXBQjzLpekYi6kCxZJBKl/ia/d38aFDvQ/1O9bTCsi3KRenY/zxhVlev7fG+EqG3V0+XIrMlbnHMzHf7lPpD7ko6FatrxIgQxE3RwdCXGpyFpUlwfHBEMvJAmGvynw8x2Cbh3OTMY6XE0AVfKpCKm8FaE6zWuaEENxYSHBqOIwJLCXyqIrElbkEfSE38azG4f4gV9Ylai5Oxzg1HObSTBxNN5EEKLJEh9/JXIPEUoXz0zHavU5ODvn5y6sL/Mxnb9vf29flIZrVyWs6iXwJSZJaAdmiFsMw+a8v3+NTF+fwOGRShRJjIxHOTUY5NhDi9GiEN+4/Wv/n7i4fDkliKVVgci1b872paA6/29H0uQd6AnZD+Ewsi2FiZ3vX1/oaFf+d5S1ko6byuXLSaCVVwCEJewQLoKibnJ2McXIwVN6qm0ysZQm5HWTKrX19ITdz8Rx9ITemaTKfyGOasJIusJIu4JAFxwZCqIqECcQyRbqCDjv4J1cfrv7ZiFZAvo1YSRX4sT+6QrqgMbma4Uh/iFRO48xElJ2dPoq6gf4YEhMdfifnJmMNz2GH+gJcnUsS8ars6PBybc5aOQ/1B7kxn0Q
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"nl_map = gpd.read_file(r'../data/shape/Netherlands_shapefile/gemeente_2021_v1.shp')\n",
"nl_map.to_crs(epsg=4326).plot()\n"
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>GM_CODE</th>\n",
" <th>GM_NAAM</th>\n",
" <th>H2O</th>\n",
" <th>OAD</th>\n",
" <th>STED</th>\n",
" <th>BEV_DICHTH</th>\n",
" <th>AANT_INW</th>\n",
" <th>AANT_MAN</th>\n",
" <th>AANT_VROUW</th>\n",
" <th>P_00_14_JR</th>\n",
" <th>...</th>\n",
" <th>P_TURKIJE</th>\n",
" <th>P_OVER_NW</th>\n",
" <th>OPP_TOT</th>\n",
" <th>OPP_LAND</th>\n",
" <th>OPP_WATER</th>\n",
" <th>JRSTATCODE</th>\n",
" <th>JAAR</th>\n",
" <th>Shape_Leng</th>\n",
" <th>Shape_Area</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>GM0034</td>\n",
" <td>Almere</td>\n",
" <td>JA</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>...</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>2021GM0034</td>\n",
" <td>2021</td>\n",
" <td>122665.358635</td>\n",
" <td>1.095623e+08</td>\n",
" <td>MULTIPOLYGON (((150213.998 479503.726, 150087....</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>GM0050</td>\n",
" <td>Zeewolde</td>\n",
" <td>JA</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>...</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>2021GM0050</td>\n",
" <td>2021</td>\n",
" <td>83290.773617</td>\n",
" <td>1.633050e+07</td>\n",
" <td>MULTIPOLYGON (((170588.413 486792.192, 170570....</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
" <td>GM0060</td>\n",
" <td>Ameland</td>\n",
" <td>JA</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>...</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>2021GM0060</td>\n",
" <td>2021</td>\n",
" <td>156406.894023</td>\n",
" <td>2.085914e+08</td>\n",
" <td>POLYGON ((196000.000 610000.000, 196000.000 60...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>3</th>\n",
" <td>GM0072</td>\n",
" <td>Harlingen</td>\n",
" <td>JA</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>...</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>2021GM0072</td>\n",
" <td>2021</td>\n",
" <td>104274.683035</td>\n",
" <td>3.618323e+08</td>\n",
" <td>MULTIPOLYGON (((158392.775 580357.500, 158387....</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>GM0088</td>\n",
" <td>Schiermonnikoog</td>\n",
" <td>JA</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>...</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>-99999999</td>\n",
" <td>2021GM0088</td>\n",
" <td>2021</td>\n",
" <td>152568.858991</td>\n",
" <td>1.624101e+08</td>\n",
" <td>POLYGON ((219000.000 616567.418, 219000.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>430</th>\n",
" <td>GM1966</td>\n",
" <td>Het Hogeland</td>\n",
" <td>NEE</td>\n",
" <td>414</td>\n",
" <td>5</td>\n",
" <td>99</td>\n",
" <td>47834</td>\n",
" <td>24052</td>\n",
" <td>23782</td>\n",
" <td>15</td>\n",
" <td>...</td>\n",
" <td>0</td>\n",
" <td>3</td>\n",
" <td>90308</td>\n",
" <td>48249</td>\n",
" <td>42059</td>\n",
" <td>2021GM1966</td>\n",
" <td>2021</td>\n",
" <td>321301.166576</td>\n",
" <td>4.875998e+08</td>\n",
" <td>MULTIPOLYGON (((217037.735 601967.991, 217043....</td>\n",
" </tr>\n",
" <tr>\n",
" <th>431</th>\n",
" <td>GM1969</td>\n",
" <td>Westerkwartier</td>\n",
" <td>NEE</td>\n",
" <td>476</td>\n",
" <td>5</td>\n",
" <td>176</td>\n",
" <td>63678</td>\n",
" <td>32034</td>\n",
" <td>31644</td>\n",
" <td>17</td>\n",
" <td>...</td>\n",
" <td>0</td>\n",
" <td>3</td>\n",
" <td>36887</td>\n",
" <td>36269</td>\n",
" <td>618</td>\n",
" <td>2021GM1969</td>\n",
" <td>2021</td>\n",
" <td>99030.762281</td>\n",
" <td>3.688457e+08</td>\n",
" <td>POLYGON ((215186.661 595044.495, 215190.788 59...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>432</th>\n",
" <td>GM1970</td>\n",
" <td>Noardeast-Fryslân</td>\n",
" <td>NEE</td>\n",
" <td>463</td>\n",
" <td>5</td>\n",
" <td>120</td>\n",
" <td>45481</td>\n",
" <td>22879</td>\n",
" <td>22602</td>\n",
" <td>17</td>\n",
" <td>...</td>\n",
" <td>0</td>\n",
" <td>2</td>\n",
" <td>51645</td>\n",
" <td>37783</td>\n",
" <td>13862</td>\n",
" <td>2021GM1970</td>\n",
" <td>2021</td>\n",
" <td>238426.959926</td>\n",
" <td>3.837930e+08</td>\n",
" <td>MULTIPOLYGON (((207769.575 603284.528, 207773....</td>\n",
" </tr>\n",
" <tr>\n",
" <th>433</th>\n",
" <td>GM1978</td>\n",
" <td>Molenlanden</td>\n",
" <td>NEE</td>\n",
" <td>393</td>\n",
" <td>5</td>\n",
" <td>243</td>\n",
" <td>44130</td>\n",
" <td>22317</td>\n",
" <td>21813</td>\n",
" <td>18</td>\n",
" <td>...</td>\n",
" <td>0</td>\n",
" <td>2</td>\n",
" <td>19158</td>\n",
" <td>18173</td>\n",
" <td>986</td>\n",
" <td>2021GM1978</td>\n",
" <td>2021</td>\n",
" <td>88570.376691</td>\n",
" <td>1.915841e+08</td>\n",
" <td>POLYGON ((123569.533 440132.167, 123576.153 44...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>434</th>\n",
" <td>GM1979</td>\n",
" <td>Eemsdelta</td>\n",
" <td>NEE</td>\n",
" <td>692</td>\n",
" <td>4</td>\n",
" <td>170</td>\n",
" <td>45587</td>\n",
" <td>22841</td>\n",
" <td>22746</td>\n",
" <td>14</td>\n",
" <td>...</td>\n",
" <td>2</td>\n",
" <td>5</td>\n",
" <td>36407</td>\n",
" <td>26789</td>\n",
" <td>9618</td>\n",
" <td>2021GM1979</td>\n",
" <td>2021</td>\n",
" <td>173966.493452</td>\n",
" <td>2.727771e+08</td>\n",
" <td>MULTIPOLYGON (((269190.000 594253.046, 269190....</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>435 rows × 38 columns</p>\n",
"</div>"
],
"text/plain": [
" GM_CODE GM_NAAM H2O OAD STED BEV_DICHTH AANT_INW \\\n",
"0 GM0034 Almere JA -99999999 -99999999 -99999999 -99999999 \n",
"1 GM0050 Zeewolde JA -99999999 -99999999 -99999999 -99999999 \n",
"2 GM0060 Ameland JA -99999999 -99999999 -99999999 -99999999 \n",
"3 GM0072 Harlingen JA -99999999 -99999999 -99999999 -99999999 \n",
"4 GM0088 Schiermonnikoog JA -99999999 -99999999 -99999999 -99999999 \n",
".. ... ... ... ... ... ... ... \n",
"430 GM1966 Het Hogeland NEE 414 5 99 47834 \n",
"431 GM1969 Westerkwartier NEE 476 5 176 63678 \n",
"432 GM1970 Noardeast-Fryslân NEE 463 5 120 45481 \n",
"433 GM1978 Molenlanden NEE 393 5 243 44130 \n",
"434 GM1979 Eemsdelta NEE 692 4 170 45587 \n",
"\n",
" AANT_MAN AANT_VROUW P_00_14_JR ... P_TURKIJE P_OVER_NW OPP_TOT \\\n",
"0 -99999999 -99999999 -99999999 ... -99999999 -99999999 -99999999 \n",
"1 -99999999 -99999999 -99999999 ... -99999999 -99999999 -99999999 \n",
"2 -99999999 -99999999 -99999999 ... -99999999 -99999999 -99999999 \n",
"3 -99999999 -99999999 -99999999 ... -99999999 -99999999 -99999999 \n",
"4 -99999999 -99999999 -99999999 ... -99999999 -99999999 -99999999 \n",
".. ... ... ... ... ... ... ... \n",
"430 24052 23782 15 ... 0 3 90308 \n",
"431 32034 31644 17 ... 0 3 36887 \n",
"432 22879 22602 17 ... 0 2 51645 \n",
"433 22317 21813 18 ... 0 2 19158 \n",
"434 22841 22746 14 ... 2 5 36407 \n",
"\n",
" OPP_LAND OPP_WATER JRSTATCODE JAAR Shape_Leng Shape_Area \\\n",
"0 -99999999 -99999999 2021GM0034 2021 122665.358635 1.095623e+08 \n",
"1 -99999999 -99999999 2021GM0050 2021 83290.773617 1.633050e+07 \n",
"2 -99999999 -99999999 2021GM0060 2021 156406.894023 2.085914e+08 \n",
"3 -99999999 -99999999 2021GM0072 2021 104274.683035 3.618323e+08 \n",
"4 -99999999 -99999999 2021GM0088 2021 152568.858991 1.624101e+08 \n",
".. ... ... ... ... ... ... \n",
"430 48249 42059 2021GM1966 2021 321301.166576 4.875998e+08 \n",
"431 36269 618 2021GM1969 2021 99030.762281 3.688457e+08 \n",
"432 37783 13862 2021GM1970 2021 238426.959926 3.837930e+08 \n",
"433 18173 986 2021GM1978 2021 88570.376691 1.915841e+08 \n",
"434 26789 9618 2021GM1979 2021 173966.493452 2.727771e+08 \n",
"\n",
" geometry \n",
"0 MULTIPOLYGON (((150213.998 479503.726, 150087.... \n",
"1 MULTIPOLYGON (((170588.413 486792.192, 170570.... \n",
"2 POLYGON ((196000.000 610000.000, 196000.000 60... \n",
"3 MULTIPOLYGON (((158392.775 580357.500, 158387.... \n",
"4 POLYGON ((219000.000 616567.418, 219000.000 61... \n",
".. ... \n",
"430 MULTIPOLYGON (((217037.735 601967.991, 217043.... \n",
"431 POLYGON ((215186.661 595044.495, 215190.788 59... \n",
"432 MULTIPOLYGON (((207769.575 603284.528, 207773.... \n",
"433 POLYGON ((123569.533 440132.167, 123576.153 44... \n",
"434 MULTIPOLYGON (((269190.000 594253.046, 269190.... \n",
"\n",
"[435 rows x 38 columns]"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"nl_map"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAGMAAAI/CAYAAABj6GcwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAANdklEQVR4nO3cb2xd913H8fenTsu60nWJaluhXXE1RWiI/alklaI8gYVMEUEkmlRpE0N+UC1PmLRqSMwwJDSJCfNk4glPAkxYasdaaeSP1qnMMlQTqCt1SGhaWvA2vNAkit2WUlindm2/PPCJ8BI799rX9/oTn89Lqs7/nJ/y7jnnHic3qirCww1bPYD4f4lhJDGMJIaRxDCSGEZ2DPJkt99+e42NjQ3ylHZOnTr1UlUNr7ZtoDHGxsaYm5sb5CntSPrBWttymzKSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMIwP92yHXs+Onz/O5R87wzhXrd777Rl59/cf8zHtv5vyrP7rquE/ddxd/dPiDXZ0jMdbwm3/+JP/4vVc67vdfr/8YYNUQAA995xxAV0Fym1pFtyG69dB3znH35GMcP33+mvvlyljh+OnzfP7rz/DGW1fejHpXwIOPnOGGm9+za619EqNx/PR5HnzkTN/Ps+O2kbvX2pbbVGMQITpJDCOJAR0frIOSGMAXjp3d6iEAiQHAD998e6uHACTGwL3z5o9eW2tb62OMTT420PO99cr5+bW2tTrGoEN00toYbiGgpTEcQ0ALY/zil2a2eghratXPplyviMtaEWP/l59gfvGHWz0MFqYOoj9Ze/u2j/GhP3yc197Y+pe6hamDHffZ1s+MPzh+1iLE3vev+UcYP2Fbx3i4+SPPrbT3/bt4+NO/1NW+2/o2tVX/4N+7hsQLX/q1dR+3ra+MrbKREJAYVhLDyLaO8Z6fGhr4OUdvvWnDx27rB/gzXzzQ9XvGnpFbmPncL1+1fr1v7U99Yf+69l9pW8eA5SCXrfYmfvONN/DHH/8Qh++5Y9XjF6YODuzHKNs+xkqr/Z/vZFs/MzbLoJ49idGFlbe6a+nl4Q1dxpC0IOmspDOS5pp1uyTNSJpvpjt7Gom5hamD6BrbR2+9qaeHN6zvmfErVfXSiuVJYLaqpiRNNsuf72k05v6ji5+89qKX29QhYLqZnwYO9zyalus2RgHfknRK0pFm3WhVXQRopiP9GGCbdHub2ltVFySNADOSXuj2BE28IwB33XXXBobYHl1dGVV1oZkuAseAe4FLknYDNNPFNY49WlXjVTU+PLzqvzodjY4xJN0i6dbL88DHgGeBk8BEs9sEcKJfg2yLbm5To8AxSZf3/2pVPS7paeBRSQ8A54D7+zfMdugYo6q+D3x4lfUvA/v6Mai2yhu4kcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGuo4haUjSaUnfaJZ3SZqRNN9Md/ZvmO2wnivjs8DzK5Yngdmq2gPMNsvRg65iSLoTOAj8xYrVh4DpZn4aOLypI2uhbq+MPwV+F3hnxbrRqroI0ExHNndo7dMxhqRfBxar6tRGTiDpiKQ5SXNLS0sb+SVao5srYy/wG5IWgK8BH5X0EHBJ0m6AZrq42sFVdbSqxqtqfHh4eJOGvT11jFFVv1dVd1bVGPAJ4O+q6lPASWCi2W0CONG3UbZEL+8ZU8B+SfPA/mY5erBjPTtX1RPAE838y8C+zR9Se+UN3EhiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMI4lhJDGMJIaRxDCSGEYSw0hiGEkMIx1jSHqXpH+S9C+SnpP0xWb9Lkkzkuab6c7+D3d76+bKeAP4aFV9GPgIcEDSfcAkMFtVe4DZZjl60DFGLfvfZvHG5r8CDgHTzfpp4HA/BtgmXT0zJA1JOgMsAjNV9RQwWlUXAZrpSN9G2RJdxaiqt6vqI8CdwL2SfqHbE0g6ImlO0tzS0tIGh9kO6/o0VVWvAk8AB4BLknYDNNPFNY45WlXjVTU+PDzc22i3uW4+TQ1Lem8zfzPwq8ALwElgotltAjjRpzG2xo4u9tkNTEsaYjneo1X1DUlPAo9KegA4B9zfx3G2QscYVfUMcM8q618G9vVjUG2VN3AjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpHEMJIYRhLDSGIYSQwjiWEkMYwkhpGOMSS9T9LfS3pe0nOSPtus3yVpRtJ8M93Z/+Fub91cGW8Bv1NVHwDuA35b0s8Dk8BsVe0BZpvl6EHHGFV1sar+uZn/H+B54A7gEDDd7DYNHO7TGFtjXc8MSWPAPcBTwGhVXYTlYMDIpo+uZbqOIemnga8DD1bVa+s47oikOUlzS0tLGxlja3QVQ9KNLId4uKr+pll9SdLuZvtuYHG1Y6vqaFWNV9X48PDwZox52+rm05SAvwSer6ovr9h0Epho5ieAE5s/vHbZ0cU+e4HfAs5KOtOs+31gCnhU0gPAOeD+voywRTrGqKp/ALTG5n2bO5x2yxu4kcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGOsaQ9BVJi5KeXbFul6QZSfPNdGd/h9kO3VwZfwUcuGLdJDBbVXuA2WY5etQxRlV9G3jlitWHgOlmfho4vLnDaqeNPjNGq+oiQDMd2bwhtVffH+CSjkiakzS3tLTU79Nd1zYa45Kk3QDNdHGtHavqaFWNV9X48PDwBk/XDhuNcRKYaOYngBObM5x26+aj7V8DTwI/J+lFSQ8AU8B+SfPA/mY5erSj0w5V9ck1Nu3b5LG0Xt7AjSSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQwkhhGEsNIYhhJDCOJYSQxjCSGkcQw0lM
"text/plain": [
"<Figure size 720x720 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import matplotlib.pyplot as plt\n",
"\n",
"fig, ax = plt.subplots(figsize = (10,10))\n",
"nl_map.to_crs(epsg=4326).plot(ax=ax, color='lightgrey')\n",
"df_geojson.plot(ax=ax)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We hebben hier een probleem dat sommige punten de coördinaten 0,0 hebben, dit is niet iets wat we willen weergeven. Laten we eerst kijken waarom dit het geval is in de data voordat we het er helemaal uithalen."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Stembureau 32\n",
"Adres \n",
"Locatie SB32\n",
"description Stembureau Mobiel Stembureau 1\n",
"Geldige stemmen 99\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 2520, Length: 973, dtype: object\n",
"Stembureau 33\n",
"Adres \n",
"Locatie SB33\n",
"description Stembureau Mobiel Stembureau 2\n",
"Geldige stemmen 115\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 2521, Length: 973, dtype: object\n",
"Stembureau 157\n",
"Adres \n",
"Locatie SB157\n",
"description Stembureau Stembureau Stembus\n",
"Geldige stemmen 72\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 2918, Length: 973, dtype: object\n",
"Stembureau 171\n",
"Adres \n",
"Locatie SB171\n",
"description Stembureau Stembureau Stembureau op locatie al...\n",
"Geldige stemmen 174\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 3579, Length: 973, dtype: object\n",
"Stembureau 172\n",
"Adres \n",
"Locatie SB172\n",
"description Stembureau Stembureau Stembureau op locatie al...\n",
"Geldige stemmen 66\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 3580, Length: 973, dtype: object\n",
"Stembureau 173\n",
"Adres \n",
"Locatie SB173\n",
"description Stembureau Stembureau Stembureau op locatie al...\n",
"Geldige stemmen 162\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 3581, Length: 973, dtype: object\n",
"Stembureau 174\n",
"Adres \n",
"Locatie SB174\n",
"description Stembureau Stembureau Stembureau op locatie al...\n",
"Geldige stemmen 65\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 3582, Length: 973, dtype: object\n",
"Stembureau 175\n",
"Adres \n",
"Locatie SB175\n",
"description Stembureau Stembureau Stembureau op locatie al...\n",
"Geldige stemmen 121\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 3583, Length: 973, dtype: object\n",
"Stembureau 39\n",
"Adres \n",
"Locatie SB39\n",
"description Stembureau Mobiel stembureau met beperkte toeg...\n",
"Geldige stemmen 56\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 5494, Length: 973, dtype: object\n",
"Stembureau 34\n",
"Adres \n",
"Locatie SB34\n",
"description Stembureau Stembureau Mobiel stembureau verzor...\n",
"Geldige stemmen 119\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 5784, Length: 973, dtype: object\n",
"Stembureau 35\n",
"Adres \n",
"Locatie SB35\n",
"description Stembureau Stembureau Mobiel stembureau verzor...\n",
"Geldige stemmen 151\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 5785, Length: 973, dtype: object\n",
"Stembureau 11\n",
"Adres \n",
"Locatie SB11\n",
"description Stembureau Mobiel stembureau Ommedijk\n",
"Geldige stemmen 14\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 6260, Length: 973, dtype: object\n",
"Stembureau 12\n",
"Adres \n",
"Locatie SB12\n",
"description Stembureau Mobiel stembureau Van Alphenstaete ...\n",
"Geldige stemmen 68\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 6261, Length: 973, dtype: object\n",
"Stembureau 70\n",
"Adres \n",
"Locatie SB70\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 34\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 6738, Length: 973, dtype: object\n",
"Stembureau 71\n",
"Adres \n",
"Locatie SB71\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 43\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 6739, Length: 973, dtype: object\n",
"Stembureau 699\n",
"Adres \n",
"Locatie SB699\n",
"description Stembureau Stembureau Mobiel stembureau\n",
"Geldige stemmen 59\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 6863, Length: 973, dtype: object\n",
"Stembureau 32\n",
"Adres \n",
"Locatie SB32\n",
"description Stembureau Mobielstembureau\n",
"Geldige stemmen 63\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7348, Length: 973, dtype: object\n",
"Stembureau 90\n",
"Adres \n",
"Locatie SB90\n",
"description Stembureau Roll and Vote\n",
"Geldige stemmen 339\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7364, Length: 973, dtype: object\n",
"Stembureau 19\n",
"Adres \n",
"Locatie SB19\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 77\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7857, Length: 973, dtype: object\n",
"Stembureau 1428\n",
"Adres \n",
"Locatie SB1428\n",
"description Stembureau Drive-thru parkeerplaats Rusheuvel 1\n",
"Geldige stemmen 163\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7912, Length: 973, dtype: object\n",
"Stembureau 1437\n",
"Adres \n",
"Locatie SB1437\n",
"description Stembureau Drive-thru parkeerplaats Rusheuvel 2\n",
"Geldige stemmen 239\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7913, Length: 973, dtype: object\n",
"Stembureau 1528\n",
"Adres \n",
"Locatie SB1528\n",
"description Stembureau Drive-thru parkeerplaats Rusheuvel 1\n",
"Geldige stemmen 323\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7923, Length: 973, dtype: object\n",
"Stembureau 1537\n",
"Adres \n",
"Locatie SB1537\n",
"description Stembureau Drive-thru parkeerplaats Rusheuvel 2\n",
"Geldige stemmen 208\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7924, Length: 973, dtype: object\n",
"Stembureau 1541\n",
"Adres \n",
"Locatie SB1541\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 125\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7926, Length: 973, dtype: object\n",
"Stembureau 200\n",
"Adres \n",
"Locatie SB200\n",
"description Stembureau Mobiel stembureau Ma\n",
"Geldige stemmen 30\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7958, Length: 973, dtype: object\n",
"Stembureau 201\n",
"Adres \n",
"Locatie SB201\n",
"description Stembureau Mobiel stembureau Di\n",
"Geldige stemmen 30\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7959, Length: 973, dtype: object\n",
"Stembureau 202\n",
"Adres \n",
"Locatie SB202\n",
"description Stembureau Mobiel stembureau Woe\n",
"Geldige stemmen 44\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 7960, Length: 973, dtype: object\n",
"Stembureau 30\n",
"Adres \n",
"Locatie SB30\n",
"description Stembureau Stembureau voor de telling stembure...\n",
"Geldige stemmen 2004\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8108, Length: 973, dtype: object\n",
"Stembureau 40\n",
"Adres \n",
"Locatie SB40\n",
"description Stembureau Stembureau voor de telling stembure...\n",
"Geldige stemmen 2953\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8109, Length: 973, dtype: object\n",
"Stembureau 990\n",
"Adres \n",
"Locatie SB990\n",
"description Stembureau Mobiel stembureau (meerdere plaatsen)\n",
"Geldige stemmen 47\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8161, Length: 973, dtype: object\n",
"Stembureau 991\n",
"Adres \n",
"Locatie SB991\n",
"description Stembureau Mobiel stembureau (meerdere plaatsen)\n",
"Geldige stemmen 70\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8162, Length: 973, dtype: object\n",
"Stembureau 992\n",
"Adres \n",
"Locatie SB992\n",
"description Stembureau Mobiel stembureau (meerdere plaatsen)\n",
"Geldige stemmen 39\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8163, Length: 973, dtype: object\n",
"Stembureau 17\n",
"Adres \n",
"Locatie SB17\n",
"description Stembureau Mobiel Stembureau dinsdag\n",
"Geldige stemmen 115\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 8318, Length: 973, dtype: object\n",
"Stembureau 228\n",
"Adres \n",
"Locatie SB228\n",
"description Stembureau Stembureau Drive-in\n",
"Geldige stemmen 80\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9093, Length: 973, dtype: object\n",
"Stembureau 30\n",
"Adres \n",
"Locatie SB30\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 25\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9151, Length: 973, dtype: object\n",
"Stembureau 31\n",
"Adres \n",
"Locatie SB31\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 8\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9152, Length: 973, dtype: object\n",
"Stembureau 32\n",
"Adres \n",
"Locatie SB32\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 24\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9153, Length: 973, dtype: object\n",
"Stembureau 76\n",
"Adres \n",
"Locatie SB76\n",
"description Stembureau Rondvaartboot Toerist VI\n",
"Geldige stemmen 148\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9661, Length: 973, dtype: object\n",
"Stembureau 68\n",
"Adres \n",
"Locatie SB68\n",
"description Stembureau Drive Inn (evenemententerrein Het L...\n",
"Geldige stemmen 285\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 9853, Length: 973, dtype: object\n",
"Stembureau 14\n",
"Adres \n",
"Locatie SB14\n",
"description Stembureau Mobiel Stembureau\n",
"Geldige stemmen 24\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10427, Length: 973, dtype: object\n",
"Stembureau 18\n",
"Adres \n",
"Locatie SB18\n",
"description Stembureau Mobiel stembureau\n",
"Geldige stemmen 50\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10431, Length: 973, dtype: object\n",
"Stembureau 23\n",
"Adres \n",
"Locatie SB23\n",
"description Stembureau Stembureau Mobiel: Wkp Zout Veere\n",
"Geldige stemmen 55\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10497, Length: 973, dtype: object\n",
"Stembureau 24\n",
"Adres \n",
"Locatie SB24\n",
"description Stembureau Stembureau Mobiel: Okp Dom\n",
"Geldige stemmen 35\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10498, Length: 973, dtype: object\n",
"Stembureau 115\n",
"Adres \n",
"Locatie SB115\n",
"description Stembureau Tent Nijverheidslaan\n",
"Geldige stemmen 258\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10525, Length: 973, dtype: object\n",
"Stembureau 116\n",
"Adres \n",
"Locatie SB116\n",
"description Stembureau Tent Nijverheidslaan\n",
"Geldige stemmen 377\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10526, Length: 973, dtype: object\n",
"Stembureau 117\n",
"Adres \n",
"Locatie SB117\n",
"description Stembureau Tent Peter Benenson park\n",
"Geldige stemmen 285\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10527, Length: 973, dtype: object\n",
"Stembureau 118\n",
"Adres \n",
"Locatie SB118\n",
"description Stembureau Tent Peter Benenson park\n",
"Geldige stemmen 329\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10528, Length: 973, dtype: object\n",
"Stembureau 310\n",
"Adres \n",
"Locatie SB310\n",
"description Stembureau Mobiel Bijzonder stembureau 1\n",
"Geldige stemmen 26\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10623, Length: 973, dtype: object\n",
"Stembureau 311\n",
"Adres \n",
"Locatie SB311\n",
"description Stembureau Mobiel Bijzonder stembureau 2\n",
"Geldige stemmen 55\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10624, Length: 973, dtype: object\n",
"Stembureau 88\n",
"Adres \n",
"Locatie SB88\n",
"description Stembureau Mobiel stembureau 1\n",
"Geldige stemmen 78\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10894, Length: 973, dtype: object\n",
"Stembureau 89\n",
"Adres \n",
"Locatie SB89\n",
"description Stembureau Mobiel stembureau 2\n",
"Geldige stemmen 35\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 10895, Length: 973, dtype: object\n",
"Stembureau 31\n",
"Adres \n",
"Locatie SB31\n",
"description Stembureau Stembureau Mobiel Stembureau\n",
"Geldige stemmen 32\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 11032, Length: 973, dtype: object\n",
"Stembureau 0\n",
"Adres \n",
"Locatie SB0\n",
"description Stembureau -\n",
"Geldige stemmen 0\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 11101, Length: 973, dtype: object\n",
"Stembureau 132\n",
"Adres \n",
"Locatie SB132\n",
"description Stembureau Stembureau IJsselhallen drive thoug...\n",
"Geldige stemmen 52\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12038, Length: 973, dtype: object\n",
"Stembureau 133\n",
"Adres \n",
"Locatie SB133\n",
"description Stembureau Stembureau IJsselhallen drive throu...\n",
"Geldige stemmen 48\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12039, Length: 973, dtype: object\n",
"Stembureau 232\n",
"Adres \n",
"Locatie SB232\n",
"description Stembureau Stembureau IJsselhallen drive thoug...\n",
"Geldige stemmen 113\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12053, Length: 973, dtype: object\n",
"Stembureau 233\n",
"Adres \n",
"Locatie SB233\n",
"description Stembureau Stembureau IJsselhallen drive throu...\n",
"Geldige stemmen 73\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12054, Length: 973, dtype: object\n",
"Stembureau 358\n",
"Adres \n",
"Locatie SB358\n",
"description Stembureau Stembureau Mobiel stembureau 1\n",
"Geldige stemmen 70\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12118, Length: 973, dtype: object\n",
"Stembureau 359\n",
"Adres \n",
"Locatie SB359\n",
"description Stembureau Stembureau Mobiel stembureau 2\n",
"Geldige stemmen 169\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12119, Length: 973, dtype: object\n",
"Stembureau 360\n",
"Adres \n",
"Locatie SB360\n",
"description Stembureau Stembureau Mobiel stembureau 3\n",
"Geldige stemmen 89\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12120, Length: 973, dtype: object\n",
"Stembureau 361\n",
"Adres \n",
"Locatie SB361\n",
"description Stembureau Stembureau Mobiel stembureau 4\n",
"Geldige stemmen 36\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie NaN\n",
"RAADSGROEPERING ''BOSCH-BELANG'' NaN\n",
"gewoon ge-DREVEN NaN\n",
"VOOR Den Bosch Joep Gersjes NaN\n",
"geometry POINT (0 0)\n",
"Name: 12121, Length: 973, dtype: object\n",
"Stembureau 114\n",
"Adres \n",
"Locatie SB114\n",
"description Stembureau Stembureau Mobiel stembureau (beper...\n",
"Geldige stemmen 109\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie 11.0\n",
"RAADSGROEPERING ''BOSCH-BELANG'' 8.0\n",
"gewoon ge-DREVEN 0.0\n",
"VOOR Den Bosch Joep Gersjes 1.0\n",
"geometry POINT (0 0)\n",
"Name: 12504, Length: 973, dtype: object\n",
"Stembureau 115\n",
"Adres \n",
"Locatie SB115\n",
"description Stembureau Stembureau Mobiel Stembureau (beper...\n",
"Geldige stemmen 80\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie 9.0\n",
"RAADSGROEPERING ''BOSCH-BELANG'' 5.0\n",
"gewoon ge-DREVEN 0.0\n",
"VOOR Den Bosch Joep Gersjes 3.0\n",
"geometry POINT (0 0)\n",
"Name: 12505, Length: 973, dtype: object\n",
"Stembureau 116\n",
"Adres \n",
"Locatie SB116\n",
"description Stembureau Stembureau Mobiel Stembureau (beper...\n",
"Geldige stemmen 119\n",
" ... \n",
"\"Leefbaar 's-Hertogenbosch\" Paul Kagie 3.0\n",
"RAADSGROEPERING ''BOSCH-BELANG'' 3.0\n",
"gewoon ge-DREVEN 20.0\n",
"VOOR Den Bosch Joep Gersjes 3.0\n",
"geometry POINT (0 0)\n",
"Name: 12506, Length: 973, dtype: object\n"
]
}
],
"source": [
"for index, row in df_geojson.iterrows():\n",
" if (row.geometry.centroid.x == 0.00000):\n",
" print(row)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Het zijn dus duidelijk allemaal mobiele stembureaus, die kunnen we in een aparte dataset houden voor later, voor nu zijn ze onbelangrijk voor de kaartweergave, laten we de dataset dus opsplitsen in drie dataframes, de originele, één frame zonder de mobiele stembureaus, en de mobiele stembureaus alleen."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/tmp/ipykernel_58008/291552514.py:3: UserWarning: Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
"\n",
" mobiel_mask = df_geojson['geometry'].centroid.x == 0\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"2520\n",
"2521\n",
"2918\n",
"3579\n",
"3580\n",
"3581\n",
"3582\n",
"3583\n",
"5494\n",
"5784\n",
"5785\n",
"6260\n",
"6261\n",
"6738\n",
"6739\n",
"6863\n",
"7348\n",
"7364\n",
"7857\n",
"7912\n",
"7913\n",
"7923\n",
"7924\n",
"7926\n",
"7958\n",
"7959\n",
"7960\n",
"8108\n",
"8109\n",
"8161\n",
"8162\n",
"8163\n",
"8318\n",
"9093\n",
"9151\n",
"9152\n",
"9153\n",
"9661\n",
"9853\n",
"10427\n",
"10431\n",
"10497\n",
"10498\n",
"10525\n",
"10526\n",
"10527\n",
"10528\n",
"10623\n",
"10624\n",
"10894\n",
"10895\n",
"11032\n",
"11101\n",
"12038\n",
"12039\n",
"12053\n",
"12054\n",
"12118\n",
"12119\n",
"12120\n",
"12121\n",
"12504\n",
"12505\n",
"12506\n"
]
}
],
"source": [
"#de waarschuwing over de projectie die niet klopt kan genegeerd worden, we zoeken naar data die x = 0 is, de projectie gaat geen effect hebben op die data in Nederland\n",
"#filter de data en print welke waardes 0 hebben op x\n",
"mobiel_mask = df_geojson['geometry'].centroid.x == 0\n",
"i = 0\n",
"for item in mobiel_mask:\n",
" if item:\n",
" print(i)\n",
" i = i + 1"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" Stembureau Adres Locatie \\\n",
"0 1 9461BH SB1 \n",
"1 2 9461DA SB2 \n",
"2 3 9461JA SB3 \n",
"3 4 9451KD SB4 \n",
"4 6 9454PL SB6 \n",
"... ... ... ... \n",
"12591 703 5391AR SB703 \n",
"12592 705 5391AR SB705 \n",
"12593 750 5382KE SB750 \n",
"12594 751 5382KE SB751 \n",
"12595 752 5283KE SB752 \n",
"\n",
" description Geldige stemmen \\\n",
"0 Stembureau Gemeentehuis Gieten (postcode: 9461... 784 \n",
"1 Stembureau OBS Gieten (postcode: 9461 DA) 562 \n",
"2 Stembureau Zorgcentrum Dekelhem (postcode: 946... 566 \n",
"3 Stembureau Ontmoetingscentrum Boerhorn Rolde (... 1495 \n",
"4 Stembureau Dropshuis de Eekhof (postcode: 9454... 347 \n",
"... ... ... \n",
"12591 Stembureau Stembureau Gemeenschapshuis de Meen... 268 \n",
"12592 Stembureau Stembureau Gemeenschapshuis De Meen... 398 \n",
"12593 Stembureau Stembureau Gemeenschapshuis 't Zijl... 663 \n",
"12594 Stembureau Stembureau Gemeenschapshuis 't Zijl... 170 \n",
"12595 Stembureau Stembureau Gemeenschapshuis 't Zijl... 222 \n",
"\n",
" Opgeroepen Ongeldig Blanco Geldige stempassen \\\n",
"0 2780 3 3 700 \n",
"1 1396 0 0 518 \n",
"2 1409 2 2 516 \n",
"3 2209 2 4 1335 \n",
"4 477 0 2 298 \n",
"... ... ... ... ... \n",
"12591 0 0 1 237 \n",
"12592 0 1 0 359 \n",
"12593 2321 2 0 552 \n",
"12594 0 0 0 151 \n",
"12595 0 0 0 193 \n",
"\n",
" Geldige volmachtbewijzen ... \\\n",
"0 90 ... \n",
"1 44 ... \n",
"2 54 ... \n",
"3 166 ... \n",
"4 51 ... \n",
"... ... ... \n",
"12591 32 ... \n",
"12592 40 ... \n",
"12593 113 ... \n",
"12594 19 ... \n",
"12595 29 ... \n",
"\n",
" Nationale Bond tegen Overheidszaken - DH Haags Belang INL Den Haag \\\n",
"0 NaN NaN NaN \n",
"1 NaN NaN NaN \n",
"2 NaN NaN NaN \n",
"3 NaN NaN NaN \n",
"4 NaN NaN NaN \n",
"... ... ... ... \n",
"12591 NaN NaN NaN \n",
"12592 NaN NaN NaN \n",
"12593 NaN NaN NaN \n",
"12594 NaN NaN NaN \n",
"12595 NaN NaN NaN \n",
"\n",
" Rosmalens Belang De Bossche Groenen \\\n",
"0 NaN NaN \n",
"1 NaN NaN \n",
"2 NaN NaN \n",
"3 NaN NaN \n",
"4 NaN NaN \n",
"... ... ... \n",
"12591 105.0 2.0 \n",
"12592 174.0 6.0 \n",
"12593 47.0 19.0 \n",
"12594 8.0 11.0 \n",
"12595 10.0 2.0 \n",
"\n",
" \"Leefbaar 's-Hertogenbosch\" Paul Kagie \\\n",
"0 NaN \n",
"1 NaN \n",
"2 NaN \n",
"3 NaN \n",
"4 NaN \n",
"... ... \n",
"12591 6.0 \n",
"12592 2.0 \n",
"12593 53.0 \n",
"12594 13.0 \n",
"12595 22.0 \n",
"\n",
" RAADSGROEPERING ''BOSCH-BELANG'' gewoon ge-DREVEN \\\n",
"0 NaN NaN \n",
"1 NaN NaN \n",
"2 NaN NaN \n",
"3 NaN NaN \n",
"4 NaN NaN \n",
"... ... ... \n",
"12591 0.0 33.0 \n",
"12592 3.0 31.0 \n",
"12593 2.0 62.0 \n",
"12594 1.0 22.0 \n",
"12595 2.0 28.0 \n",
"\n",
" VOOR Den Bosch Joep Gersjes geometry \n",
"0 NaN POINT (6.75899 53.00524) \n",
"1 NaN POINT (6.75990 52.99975) \n",
"2 NaN POINT (6.76600 53.00494) \n",
"3 NaN POINT (6.64736 52.98281) \n",
"4 NaN POINT (6.60459 52.95269) \n",
"... ... ... \n",
"12591 0.0 POINT (5.43290 51.72810) \n",
"12592 2.0 POINT (5.43290 51.72810) \n",
"12593 5.0 POINT (5.45919 51.70595) \n",
"12594 1.0 POINT (5.45919 51.70595) \n",
"12595 0.0 POINT (5.45919 51.70595) \n",
"\n",
"[12532 rows x 973 columns]\n"
]
}
],
"source": [
"df_geojson_clean = df_geojson[~mobiel_mask]\n",
"df_geojson_mobiel = df_geojson[mobiel_mask]\n",
"print(df_geojson_clean)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfoAAAI/CAYAAACf2BvTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAEAAElEQVR4nOz9eZwd1Xnnj39OVd2qu3Wr1epuoR0JCZBJItmWECCDLQthEEtM4niJAVvYaGzCJDYkM/Ew85v5zQzjeEFgh+CJwJYN8hLHiRKMWSxwY7AEQgJL8TiAJdS01Fp6Vav7rnVv1fn+UX2q69atU8u9txdJ5/16+WXUfZe6davrOed5Ps/nIZRSCAQCgUAgODuRpvoABAKBQCAQTBwi0AsEAoFAcBYjAr1AIBAIBGcxItALBAKBQHAWIwK9QCAQCARnMSLQCwQCgUBwFqNM9QF40dbWRs8///ypPgyBQCAQCM4IXnvttQFKabvX76ZloD///POxb9++qT4MgUAgEAjOCAgh3bzfidS9QCAQCARnMSLQCwQCgUBwFiMCvUAgEAgEZzEi0AsEAoFAcBYjAr1AIBAIBGcxItALBAKBQHAWIwK9QCAQCARnMSLQCwQCgUBwFiMCvUAgEAgEZzEi0AsEAoFAcBYjAr1AIBAIBGcxItALBAKBQHAWIwK9QCAQCARnMSLQCwQCgUBwFiMCvUAgEAgEZzEi0AsEAoFAcBYjAr1AIBAIBGcxItALBAKBQHAWIwK9QCAQCARnMSLQCwQCgUBwFiMCvUAgEAgEZzEi0AsEAoFAcBYjAr1AIBAIJhVKKSilU30Y5wzKVB+AQCAQCM4dTNPEsWPHQCnF/PnzIUlivznRiEAvEAgEAk8KhQLK5TIAgBACQkhd/63rOo4cOYJCoQAAeOedd7Bw4UIoighFE4k4uwKBQCCooFQqobe3F8PDwxP6PrlcDm+//TYWLFiAZDI5oe91LiMCvUAgEAgAWGn1gYEB9Pf3+9bQE4kEyuUySqVS3e9ZKpVw+PBhdHR0oL293c4CCBqHCPQCgUBwjkMpxfDwMHp7e+1UvReEEMTjcZimCcMwEIvFGhLsAaCvrw+ZTAbz58+HqqoNeU2BhQj0AoFAcA6TzWZx4sQJu27uRywWQz6fRyKRQDweh2EY3McqigJZlgFYC4RSqeT7eMBK5R86dAhz5sxBS0uL2N03CCF3FAgEgnOQcrmMEydOoKurK1SQZ88BgHw+j3K5DFmWkUwmEY/HKx6nKApUVYUkSTAMA5IkIRaLharDM1X+0aNHfbMLgvCIHb1AIBCcY2QyGfT09EROkWuahnw+j1QqBUqpncJngZyl8WVZRi6Xs5/HAjbb4YdhZGQE2WwW8+bNQ3Nzc6TjFFQiAr1AIBCc5VBKUS6XUSwWMTo6isHBQQCIHOjz+TxisZj9Wjycv2Otc+VyOTB178YwDBw5cgQtLS2YM2dOpIWCYJxQgZ4Q8g6AUQAGgDKldBUh5H8B+EMAJoA+AJ+mlB4P89zGHLpAIBCcvVBKceTIESQSCaTTaSQSicCatWEYKBaL0HW96v9N06x6fLlctoO9u5edLQ7cYrtYLBbJ1U5RFEiSVFcafnh4GNlsFgsXLkQikaj5dc5VSJgvbCxYr6KUDjh+1kwpHRn77z8H8C5K6efCPDeIVatW0X379oV9uEAgEJx1jI6Ooru72/63JElIJpNIpVJIJpMwDMMO5CyYhw2msVgMsVgMAOwUPHsPtuuWZdkz0BNCpsy+lhCC8847D62trUKo54IQ8hpvI11z6p4F+TFSAIRxsUAgEDQIt1mNaZrIZDLIZDKhX0OWZft/hBCYpolSqWT/D7CCezweB6UUuq7bgT6RSHi2zk2lRz2lFCdOnLBr9yKVH46wgZ4C+DkhhAL4e0rpVgAghNwH4DYApwGsi/JcgUAgEHhjGAZGRka4vyeE2O1rzCue7cyZQI79v1ddnBneSJKEYrFYIZxzvsd0ZWRkBPl8XjjqhSRse91aSul7AFwH4M8IIVcBAKX0XkrpAgDfB3BXlOe6IYRsJoTsI4Ts6+/vj/YpBAKB4CyBUorTp09X7JxZ2p7tYCmlKJVKKBQKyOVyyOVyyOfzKBaLKJVKnvV4AEgmk7ZynvW1szY4N7zXmC4wR72+vj4xCS+AUDX6iicQ8j8AZCilX3f8bBGAn1FKfy/qc70QNXqBQBAGtrObOXPmtHVTo5Qin88jn8/bP2O7ZXcqvVQqVdTZJUmCpmn4+e9O4Tv7BjCQNdCWknHbyhasW5wO9f6apkGWZRSLxardfSKRsI+L1exlWYZhGA1zvJsMkskkFi1adE6n8uuq0RNCUgAkSuno2H9fA+B/EkKWUUoPjj3sJgBvhn1urR9EIBAIAEDXdZw4cQKjo6MAgP7+fjQ1NaG1tRXpdHrK086maSKbzWJ0dBQjIyO2ul3X9YrHxeNxrlmNqqqQZRk/f2sID+7qR9GwNmX9WQMPvTIEABXB3i2SYyp6vzY4Jyywn0kBnpHL5XDkyBEsWrRIjL31IEyNfjaAHWN/OAqAH1BKnyGE/BMh5CJY7XXdAD4HAISQuQAepZRu5D238R9DIBCcC/gNXRkdHcXo6ChisRhaW1sxc+bMCRt/ylLnvP95tbN5ZU/9ghI79u+8NmgHeUbRoHhs/7Ad6CVJAqXU/n9N01AsFs+plHY2m8XRo0excOHCKV/oTTcC/woopYcBrPD4+R9zHn8cwEa/5woEAkFURkdHceLEiapdsRs2YrWvrw/Nzc1obW1FMpn0vPlTSitEbE4hm/Pffun1WmH96JIkcevhhBAMZL1NZgayht0m52yRMwwjtKXt2cbo6CiOHz+OuXPnimDvQDjjCQSCaYtpmhgZGbFr8UFp5Xg8DkIICCEoFAo4ffo0Tp8+DU3ToCiKZzCvB0IIJEmq2Jmz1+XtpuPxOHRdt48nHo97qt4Ba5fanlLQl61eWLSnlIo2OYHFqVOnIMsyOjo6RBp/DBHoBQLBtIJSimw2i+HhYYyMjNjBmBCCZDLJDYpev4vH4/bNngX3UqnkGYQ7uzJ4bP8wBrIG2tMKPrO6HVcvnWHvDL3a1yilni1sTETn3lnHYjFb6c5EcLz2MHaMt66cgYf2DKFYHj9mTSa4deUM7xN4hiPLsu2mx3QHpmnaiyq3MRAhxDb3YQwMDODUqVNobm62n9fc3BzKXfBsRAR6gUAw5VBK7R348PCwZ2qcUopcLlehFGfsOqbj0Vd7qlTpvBS2JEn45TtZbHttCP3ZMpo0CVndBCuF92XK2PLiSei67qtuj8ViFToAQohtQ+sO/qxMEIvFKo6Ll1Vgn3Hd4jQkScJ3Xx+qSXXvx3So4bOuAObCx+v9dxKPx+0FlyzL0HUdyWQShULBPs+GYeDUqVP2cwYGBtDa2orzzjvvnNvpR26vmwxEe51AcPbDAjdLzUdJQTuD/a5jOra8eLJCsEZQadWpycBdl82yg2NnVwYPvTJUJXJz056Sse3m+VU/Z2NYedkFJ7FYzDam8WIqLWX9VP8TCesocDrx1QI7d8lksiLbwrwC3Oc1Ho+jubkZHR0d9X6EaceEWOAKBAJBVFhangV3r507c3tj6XEvnMHh26/2VwVsd9gsGsCWXdbEtnWL03hs/3BgkAdgC+HYzp0NemH/CyKZTHoGnIpjpXRKAq6fCHAiIITYVruN+qzsvDKXP2B8RK6iKCCEVCwgi8Ui+vv70draOmEdGdORc+eTCgSCKYH1lJ8+fRqjo6NVOzhWR2dDWpypW6Yqd++cdV23d/X9HkI1LyiAx/efxnXL27hKdjftactmNqroTdM0O2MRhslKJbuDbVAHQ6NgaXV3yaURVGgrUgo+c2k7PrikCcViEbIs262G7P8BYGho6Kzc1fMQgV4gEDQc0zRts5jR0dGKnSMTqgHWDsu9u2O94KxXPRaL2W5tTtiOmqdK96I/W0ahUEBbSkZ/QLDXZIJbV8yoObUcJYg2ol2Phzu4T0SwdSJJkm2rywblhF3w8HAGc6dGwV2C6ctWaitYGh9AhWve8PCwCPQCgUAQhXK5bNu85vN
"text/plain": [
"<Figure size 720x720 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize = (10,10))\n",
"nl_map.to_crs(epsg=4326).plot(ax=ax, color='lightgrey')\n",
"df_geojson_clean.plot(ax=ax)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nu we een mooie kaart hebben van nederland met de stemlokalen erop, kunnen we gaan kijken naar hoe dit verhoud met de bevolkingsdichtheid als eerste voorbeeld. De kaart is opgehaald van het CBS, eerst gaan we kijken hoe de kaart er op zichzelf met dezelfde projectie uitziet."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAOQAAAD4CAYAAAD8Q0ptAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB/IklEQVR4nO39eZxe53HfiX6rnvM2dgIEQYIkSHBfRFIUJVGStViyZNmWLFl2bMdLnM3Zc+PMJLmfJHaSmTs3mcz1jOfexDc3i2U769jjceIocWxZkiVbm20tpESJ+75gIQhi37vPearuH1XPeZsSSAEEQDSItxJaQKPffk93nzpVz69+9fuJuzOLWcxiaYSe7QuYxSxmMY1ZQs5iFksoZgk5i1ksoZgl5CxmsYRilpCzmMUSiu5sX8DxYsOGDX711Vef7cuYxSzOSNx999273P3i4/3bkkzIq6++mrvuuutsX8YsZnFGQkSefrF/m7Wss5jFEopZQs5iFksoZgk5i1ksoZgl5CxmsYRilpCzmMUSillCzmIWSyhmCTmLWSyhmCXkLGaxhGJJEgNmMYuTjad3H+b/+KOnmAj8hXddz/rVy+ircfTYAsuXzTHXnRu154QSUkSeAg4CFRjc/U4R+UfA9wMG7AT+rLtvP5HXnp5Ln8UspvHwzoNs2XsYd+eXPvc4f+f9t/Af/vBJvvzELtatXs7f/O6buGTN8rN9md8yTqZCvtvddy36+8+5+/8AICL/HfA/An/lBF87i1mclnhg+36++sw+fuLbruL373+O116xlnfcEDTRH3zjlSz0lTdft+GcSEY4hZbV3Q8s+usqYKYFMovTElv2HMYdHnvuIHsOL+AYqsK+Iz3P7jnKc4eOUgfnLddfwiVrl/HIs3Er/swHbmHVso6iAsC6lXP82LddPf79XIgTTUgHPiEiDvyCu38YQET+MfCngf3Au0/mtd8YIvKXgL8EsHnz5hP/Dmbxqorf+PIz/O4Dz+IC4k6thoiAKkWVVcs63nb9JVx3yRrefM1FABw+OgBwwYrJN329dSvnXtHrP9WQExG5EpHL3X27iFwC/C7w1939s4v+/WeA5e7+/zjZ1x4v7rzzTp9te5wf8eD2Azz6/AE+9LorAPhb/+dXObzQgwAmmPeoFpZ3E370267kbdcdd2vpnAoRufvFsJQTqpANrHH3nSLyEeDNwOKk+lXgt4FvSsgTeO0szsE4cHSBf/n7j3Fovkfcef/tm7hqwyqEmKV96ck9LFQDcR5+9hB9ragKuKAK775pIyuWdew+NM/OPUf41S8+yZ94yzWgBgKXX7CcK9av5Mr1K7n5srWsnOvYsGbZKV3zE88f4kuP72HX4WMc6ys79h9j1bLCbZvWcmihMinC99x6GRefxfPmt0xIEVkFqLsfzD9/N/APReQGd380P+1DwEMn+trTd/mzONNxdKHy5K5D8Rd3brpsLUWFX/7M4zy+8yAIiMC/+fwTiJBIgqBimMc/usOyIty8aS1uTleUWzatZWEwXr/5Qv7hY7t588VrAHjPjRv51IM7+VvffRNrVp5aAgI8u/con3n0OR7afpDtew/h7rgL3aTgwMFjC2zbe5hSFHfhS0/s5jtu3sgPvuHKU37vlxMnUiE3Ah8Rkfb5v+ruHxOR3xCRm4ixx9MkwioilwO/5O7f+2KvPf3fxixOd9y3dR+ff/R5Htqxn6O94R6J94NvuILvue1y1q6aICqRf+6ICmuWd1y8eo7nD/ZMOmH5RNm8fiVv3HwRmzesetHz3I+/ZTM3bIyE/ODrr+CDr7/iZV2zu3P/tv08f3CeJ54/yBPPH2bfsQWGwXDAi8a1ujP0FRGotTKZzOEY4nCsH/ijx3Yt3YR09yeA1x3n4z/0Ip+/Hfjel3rtLJZebN93hAefPcCf/3d3A/D+Wy6mAoIz1xVed8VaLlq9nOsvWQ3Aj7zpKqo7jz57iN6cN1+7nqsuWsXNl13AiknHirlywu/dkvHp3Ye54sKVFBV2H5xnvlZUoFMFhHUrJ3RF6Qdj9+FjfP2ZvWzbN8/+Yz1P7T6E18rR6giAC078WUTAQc0wEaCgagw+IBPFbcAqaCmIKFlAzkrMmDrnYQyD8dyho3zukZ08teswO/Ye4siC4+aLPmdABESUv/VdN3PNxatf8DVWL5/wF995Awt9xYFlkxNPwGN95etb97LQG/dvP0BR4fDCwENb93P9pWuYdMpD2/dRLT7fFcThsgtW8PabLuaT9+9gz5F5xJ12ydYPaNehpYtPZpqIbo6IYCIIjkulutNpARGsQHEwDK/OwgIcWRhYOffKp8csIc+jGKrxb//wSe7fto9D8z04iAu1GloKJXPqs3/nO3h69xHA2bx+FVddtOpFv+bct0hEM+fDv/8QB44NaL7Bk7uO0Juj7jgGCJK06keeO0itFS2KIHkudRzYfuAY/+nLW3AcPEYhWoVqlTKJ5BpqT4m0Q1VBC4rhonG+zSSFHqTDXRAxUB0z/+iCsevgMTZftPr439QZjFlCnkfx8fu38+XHd+PiiEvO96IKLivKTZev5V/92bcAsHn9iyfhicbDOw7w7/7gCXbuP4xqwXHco1qJO9UFdwMtlOKA4yJkzwkCJoI6+THi383R+GxEDFFHtKAQ75MzTACGnh4QhaKKu2NulFKw8TEQZ2Ex4iMK9Sx53swS8jyJzz+2k49+7VlEAwrtivD6zespIlxz8Uqu37iWKy5ceVre695te/nEvc/x2M4DVHNECqIJqBAFyrINLQiuAa64A2aAU00ppcTniyDmoIaY4CqogeHx6RoPGCN7XAQRjcQUpfOKiVNrpS4MlLkOrGQaRuIrYCoYjgh84fHdXLNhzWn5eZxMzBLyPIl9hxZwAXfjopXL+MDrLuftN1xy2r7+sb7yqQef476te3nq+cNxHkNwHxAtYGQjCZqgiSPRSpqjIgjGgFNKAYTqNSq55etcUI3xhEmMLxCnaMHy393zCClR8XDDVRFXwOjmSiZyjzto1yGtA8YjiQ3ufmoPb7/+Yja/RLt+JuKEmDqvdMyYOmcmnnr+EM8dOMqd12w4bfzOu57aw+ce3snTuw5zbOiJtCLucAGzGmfH1naaZ/sZ50NrFGiL4111QyVQVTwSxHHUBRHHKogY7sJQK5NJh6JYa1NzGCoqU5Aq8pbB4qzcElzwaJkRrBrdZDJ+HTeY6+AvvOsG7ti8/rT8rFqcMlNnFq+OuPri1Vx98akDFWbGFx57no/eu52dB+fjTIhgxLnUEVQymQzIh7470TJ7/CWa06hsSPBWUUGz/RQYT3lIVDhVcOlQd8Qr1QzzaFWlKOqauZwnxEx6q45oVk3iHEsmvhnoRBhsQEyo6hRRFgx+8bOPc+3FO3nz1eu58dIL2Lj2zLJ4Zgn5Koyn9xziN+96hgMLlRs2ruX779hEUaErp76ke+DIAj//yYfZkigsRbCh5iEM2h+qASrUvse7Ls5ynU5XgpRMvMgLywqIKFVAXRO0cSRPnt4SycGqMSmT/HrxOWYVE3CzoOhNshgbuCTqClG9yQdES/nWDovh1agYKkIvzqM7DvDQswdYMREuW7eSGzeuYfu+Y7znNRu5ddPaU/6ZLo5ZQr6K4uCxno/c/Qx3PbWXI8cW0FLYdXCBe7fsoajy/W+4ktdvvvBlf/1Hdhzgn33yIeaz2rgkU6eU6BRVEDegJFJpTObmoiAKDFZxM4oo5tE0GlGNVAWXqI4xPIxq6llNA3FVtAExaoBGG5tnUpE4SxYKjmFDsHH6fmBu2TIkE7OhvHGmjtpuCe6Kw6SU8QHhtWLmoMpCFZ7edZindh3CBd5104ZT+XUdN2YJ+SqJh7Yf4Bc/9yiHjvV0Wrhp0zre+5pLWb96Gf/hD5/ih954Ba+5/OU/zQ/ND/zbP3ySY31sX7gHaqqQZ0MgmTEuyZEx8C7qm2b76l3ccupxhpNMyqEfmIhiMkQyFkVMEM3Kicbnu1DNECnZGgtqjmuMQzoC8AFFVTGLc2a1IccnAiVaVcm21iSuD3NcDbxEoqugJtA5Yk4dDBehK4V1Kybcevm6U/iNHT9
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"nl_map_cbs = gpd.read_file(r'../data/shape/Netherlands_shapefile/CBS_vk500_2020_v1.shp')\n",
"nl_map_cbs.to_crs(epsg=4326).plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Het verwerken van de 500m bij 500m vierkanten kaar duurde heel lang (langer dan een half uur) op de laptop, en is dus niet ideaal tenzij we dit niveau van detail écht nodig hebben. Er is namelijk nog een andere kaart met dezelfde gegevens beschikbaar, maar met een grid van 500 bij 500 meter in plaats van 100 bij 100. Laten we dus kijken hoe vaak een stemlokaal dichterbij dan 100 meter van de dichtstbijzijnde andere is. De makkelijkste manier om dat te doen zonder alle punten met alle andere te vergelijken (heel veel moeite), is een extra dataframe maken als een kopie, alle indexen 1 opschuiven (want alle stemlokalen zijn al in een volgorde van clustering), en dan de laagste afstand bekijken."
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"11165 0.002608\n",
"1021 0.004853\n",
"3196 0.004875\n",
"1923 2.693202\n",
"3757 3.609362\n",
" ... \n",
"7444 253385.654588\n",
"9711 263408.187856\n",
"4392 266925.200994\n",
"9297 271954.076837\n",
"0 NaN\n",
"Name: distance, Length: 11670, dtype: float64\n"
]
}
],
"source": [
"df_shifted = df_geojson_clean.to_crs('EPSG:28992')\n",
"\n",
"df_shifted['geometry (shifted)'] = df_shifted['geometry'].shift(periods=1)\n",
"df_shifted['distance'] = df_shifted['geometry'].distance(df_shifted['geometry (shifted)'])\n",
"df_shifted.sort_values(['distance'], inplace=True, ascending=True)\n",
"zero_mask_booth = df_shifted['distance'] == 0.000000\n",
"df_shifted = df_shifted[~zero_mask_booth]\n",
"print(df_shifted['distance'])"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We kunnen dus vaststellen dat afstanden van onder de 500 meter waarschijnlijk zeldzaam zijn, en we verder kunnen gaan met de 500 meter bij 500 meter kaart. Laten we nu dus een projectie proberen te maken met de bevolkingsdichtheid erop om het te vergelijken met de stemlokalen en hun posities en clustering."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"['c28992r500', 'INWONER', 'MAN', 'VROUW', 'INW_014', 'INW_1524', 'INW_2544', 'INW_4564', 'INW_65PL', 'P_NL_ACHTG', 'P_WE_MIG_A', 'P_NW_MIG_A', 'AANTAL_HH', 'TOTHH_EENP', 'TOTHH_MPZK', 'HH_EENOUD', 'HH_TWEEOUD', 'GEM_HH_GR', 'WONING', 'WONVOOR45', 'WON_4564', 'WON_6574', 'WON_7584', 'WON_8594', 'WON_9504', 'WON_0514', 'WON_1524', 'WON_MRGEZ', 'P_KOOPWON', 'P_HUURWON', 'WON_HCORP', 'WON_NBEW', 'WOZWONING', 'UITKMINAOW', 'OAD', 'STED', 'geometry']\n"
]
},
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>c28992r500</th>\n",
" <th>INWONER</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>267</th>\n",
" <td>E2050N6110</td>\n",
" <td>5</td>\n",
" <td>POLYGON ((205000.000 611500.000, 205500.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>269</th>\n",
" <td>E2060N6110</td>\n",
" <td>15</td>\n",
" <td>POLYGON ((206000.000 611500.000, 206500.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>292</th>\n",
" <td>E2055N6105</td>\n",
" <td>20</td>\n",
" <td>POLYGON ((205500.000 611000.000, 206000.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>293</th>\n",
" <td>E2060N6105</td>\n",
" <td>185</td>\n",
" <td>POLYGON ((206000.000 611000.000, 206500.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>294</th>\n",
" <td>E2065N6105</td>\n",
" <td>340</td>\n",
" <td>POLYGON ((206500.000 611000.000, 207000.000 61...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>...</th>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" <td>...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>151071</th>\n",
" <td>E1970N3075</td>\n",
" <td>10</td>\n",
" <td>POLYGON ((197000.000 308000.000, 197500.000 30...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>151073</th>\n",
" <td>E1980N3075</td>\n",
" <td>90</td>\n",
" <td>POLYGON ((198000.000 308000.000, 198500.000 30...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>151074</th>\n",
" <td>E1985N3075</td>\n",
" <td>15</td>\n",
" <td>POLYGON ((198500.000 308000.000, 199000.000 30...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>151100</th>\n",
" <td>E1980N3070</td>\n",
" <td>5</td>\n",
" <td>POLYGON ((198000.000 307500.000, 198500.000 30...</td>\n",
" </tr>\n",
" <tr>\n",
" <th>151106</th>\n",
" <td>E1920N3065</td>\n",
" <td>20</td>\n",
" <td>POLYGON ((192000.000 307000.000, 192500.000 30...</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"<p>70656 rows × 3 columns</p>\n",
"</div>"
],
"text/plain": [
" c28992r500 INWONER geometry\n",
"267 E2050N6110 5 POLYGON ((205000.000 611500.000, 205500.000 61...\n",
"269 E2060N6110 15 POLYGON ((206000.000 611500.000, 206500.000 61...\n",
"292 E2055N6105 20 POLYGON ((205500.000 611000.000, 206000.000 61...\n",
"293 E2060N6105 185 POLYGON ((206000.000 611000.000, 206500.000 61...\n",
"294 E2065N6105 340 POLYGON ((206500.000 611000.000, 207000.000 61...\n",
"... ... ... ...\n",
"151071 E1970N3075 10 POLYGON ((197000.000 308000.000, 197500.000 30...\n",
"151073 E1980N3075 90 POLYGON ((198000.000 308000.000, 198500.000 30...\n",
"151074 E1985N3075 15 POLYGON ((198500.000 308000.000, 199000.000 30...\n",
"151100 E1980N3070 5 POLYGON ((198000.000 307500.000, 198500.000 30...\n",
"151106 E1920N3065 20 POLYGON ((192000.000 307000.000, 192500.000 30...\n",
"\n",
"[70656 rows x 3 columns]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"print(list(nl_map_cbs))\n",
"nl_map_dst = nl_map_cbs.drop(columns=['MAN', 'VROUW', 'INW_014', 'INW_1524', 'INW_2544', 'INW_4564', 'INW_65PL', 'P_NL_ACHTG', 'P_WE_MIG_A', 'P_NW_MIG_A', 'AANTAL_HH', 'TOTHH_EENP', 'TOTHH_MPZK', 'HH_EENOUD', 'HH_TWEEOUD', 'GEM_HH_GR', 'WONING', 'WONVOOR45', 'WON_4564', 'WON_6574', 'WON_7584', 'WON_8594', 'WON_9504', 'WON_0514', 'WON_1524', 'WON_MRGEZ', 'P_KOOPWON', 'P_HUURWON', 'WON_HCORP', 'WON_NBEW', 'WOZWONING', 'UITKMINAOW', 'OAD', 'STED'])\n",
"zero_mask_pop = nl_map_dst['INWONER'] == -99997\n",
"nl_map_dst = nl_map_dst[~zero_mask_pop]\n",
"nl_map_dst"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {},
"outputs": [],
"source": [
"import folium\n",
"from folium import plugins\n",
"\n",
"map = folium.Map(location = [52.155, 5.3875], zoom_start = 9, tiles=\"cartodbdark_matter\", prefer_canvas=True)\n",
"\n",
"nl_map_dst.to_crs(epsg=4326)\n",
"#Make sure the index is a string so folium can read it correctly as a key.\n",
"nl_map_dst['c28992r500'] = nl_map_dst['c28992r500'].apply(lambda x: str(x))\n",
"\n",
"folium.Choropleth(\n",
" geo_data = nl_map_dst,\n",
" name=\"Bevolkingsdichtheid\",\n",
" data = nl_map_dst,\n",
" columns = [\"c28992r500\", \"INWONER\"],\n",
" key_on = 'feature.properties.c28992r500',\n",
" fill_color = 'RdPu',\n",
" nan_fill_color= 'white',\n",
" fill_opacity = 0.7,\n",
" nan_fill_opacity = 0.7,\n",
" line_opacity = 0,\n",
" legend_name = 'Bevolkingsdichtheid',\n",
" smooth_factor = 1.0,\n",
" show=False\n",
").add_to(map)\n",
"\n",
"# Renders the map to an HTML file and displays it in an embed.\n",
"def embed_map(m):\n",
" #from IPython.display import IFrame\n",
" m.save('index.html')\n",
" #return IFrame('index.html', width='100%', height='750px')"
]
},
{
"cell_type": "code",
"execution_count": 48,
"metadata": {},
"outputs": [],
"source": [
"df_geojson_clean_hmp = [[point.xy[1][0], point.xy[0][0]] for point in df_geojson_clean.geometry]\n",
"\n",
"plugins.HeatMap(df_geojson_clean_hmp, name=\"Stemlokalen Heatmap\").add_to(map)\n",
"\n",
"folium.LayerControl().add_to(map)\n",
"\n",
"\n",
"\n",
"embed_map(map)\n",
"\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>c28992r500</th>\n",
" <th>INWONER</th>\n",
" <th>geometry</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>267</th>\n",
" <td>E2050N6110</td>\n",
" <td>5</td>\n",
" <td>POINT (205250.000 611250.000)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>269</th>\n",
" <td>E2060N6110</td>\n",
" <td>15</td>\n",
" <td>POINT (206250.000 611250.000)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>292</th>\n",
" <td>E2055N6105</td>\n",
" <td>20</td>\n",
" <td>POINT (205750.000 610750.000)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>293</th>\n",
" <td>E2060N6105</td>\n",
" <td>185</td>\n",
" <td>POINT (206250.000 610750.000)</td>\n",
" </tr>\n",
" <tr>\n",
" <th>294</th>\n",
" <td>E2065N6105</td>\n",
" <td>340</td>\n",
" <td>POINT (206750.000 610750.000)</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" c28992r500 INWONER geometry\n",
"267 E2050N6110 5 POINT (205250.000 611250.000)\n",
"269 E2060N6110 15 POINT (206250.000 611250.000)\n",
"292 E2055N6105 20 POINT (205750.000 610750.000)\n",
"293 E2060N6105 185 POINT (206250.000 610750.000)\n",
"294 E2065N6105 340 POINT (206750.000 610750.000)"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"nl_map_dst_points = nl_map_dst\n",
"nl_map_dst_points['geometry'] = nl_map_dst_points['geometry'].centroid\n",
"nl_map_dst_points.head()"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:ylabel='Frequency'>"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA24AAANZCAYAAAB3Lm1yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAicElEQVR4nO3dfbCmd13f8c+XLEjCQ6lmsTQBF2wMUkZIWCia1kqQNhiF2ooNI9YylPQBKbS2Ghim2j86E2daHpwiNTwozwjhoZQUEBREOxDYABaSwEBDgAU0y5MBZAiBb/84d+AkObs52ZzrnO+efb1mzux9X+c+1/XdyW92897ruq+7ujsAAADMdbudHgAAAIAjE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADDcuHCrqhdV1TVV9eFNvv7nq+qKqrq8ql6x9HwAAADbraZ9jltV/XiSryZ5SXff/xZee1qSVyc5u7u/VFV37+5rtmNOAACA7TLujFt3vyvJF9dvq6ofrKq3VNVlVfUnVXXf1beemOS53f2l1c+KNgAAYNcZF26HcVGSJ3f3g5L8hyS/vdr+Q0l+qKr+T1W9p6rO2bEJAQAAFrJnpwe4JVV15yQ/luQ1VXXD5u9Z/bonyWlJfiLJqUn+pKru391f3uYxAQAAFjM+3LJ2VvDL3f3ADb53MMl7uvubST5RVR/NWsi9bxvnAwAAWNT4SyW7+9qsRdljkqTWPGD17Tckedhq+8lZu3Tyqp2YEwAAYCnjwq2qXpnk3UlOr6qDVfWEJL+Q5AlV9WdJLk/y6NXL35rkC1V1RZJ3JPmP3f2FnZgbAABgKeM+DgAAAIAbG3fGDQAAgBsbdXOSk08+ufft27fTYwAAAOyIyy677PPdvfem20eF2759+3LgwIGdHgMAAGBHVNUnN9ruUkkAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMtFm5VdXpVfXDd17VV9dSljgcAALBb7Vlqx9390SQPTJKqOiHJZ5K8fqnjAQAA7Fbbdankw5P8v+7+5DYdDwAAYNfYrnA7L8krN/pGVZ1fVQeq6sChQ4e2aRwAAIBjx+LhVlV3SPKoJK/Z6PvdfVF37+/u/Xv37l16HAAAgGPOdpxxe2SS93f3X2zDsQAAAHad7Qi3x+Ywl0kCAABwyxYNt6o6KckjkrxuyeMAAADsZot9HECSdPdfJfm+JY8BAACw223XXSUBAAA4SsINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGG7PTg9wLNh3wSVbtq+rLzx3y/YFAAAcH5xxAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIbbs9MDHG/2XXDJlu3r6gvP3bJ9AQAAcznjBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADDcouFWVXerqour6iNVdWVV/eiSxwMAANiN9iy8/+ckeUt3/1xV3SHJSQsfDwAAYNdZLNyq6q5JfjzJP0+S7r4uyXVLHQ8AAGC3WvJSyfskOZTkd6vqA1X1gqq6001fVFXnV9WBqjpw6NChBccBAAA4Ni0ZbnuSnJnked19RpKvJbngpi/q7ou6e39379+7d++C4wAAAByblgy3g0kOdvelq+cXZy3kAAAAuBUWC7fu/vMkn66q01ebHp7kiqWOBwAAsFstfVfJJyd5+eqOklclefzCxwMAANh1Fg237v5gkv1LHgMAAGC3W/QDuAEAALjthBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOH2LLnzqro6yVeSfCvJ9d29f8njAQAA7EaLhtvKw7r789twHAAAgF3JpZIAAADDLR1uneQPquqyqjp/4WMBAADsSktfKnlWd3+2qu6e5G1V9ZHuftf6F6yC7vwkude97rXwOAAAAMeeRc+4dfdnV79ek+T1SR6ywWsu6u793b1/7969S44DAABwTFos3KrqTlV1lxseJ/kHST681PEAAAB2qyUvlfz+JK+vqhuO84rufsuCxwMAANiVFgu37r4qyQOW2j8AAMDxwscBAAAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQAAYDjhBgAAMJxwAwAAGE64AQAADCfcAAAAhhNuAAAAwwk3AACA4YQbAADAcMINAABgOOEGAAAwnHADAAAYTrgBAAAMJ9wAAACGE24AAADDCTcAAIDhhBsAAMBwwg0AAGA44QYAADCccAMAABhOuAEAAAwn3AAAAIYTbgAAAMMJNwAAgOGEGwAAwHDCDQA
"text/plain": [
"<Figure size 1080x1080 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"import numpy as np\n",
"\n",
"from scipy.spatial import cKDTree\n",
"from shapely.geometry import Point\n",
"\n",
"def ckdnearest(gdA, gdB):\n",
"\n",
" nA = np.array(list(gdA.geometry.apply(lambda x: (x.x, x.y))))\n",
" nB = np.array(list(gdB.geometry.apply(lambda x: (x.x, x.y))))\n",
" btree = cKDTree(nB)\n",
" dist, idx = btree.query(nA, k=1)\n",
" gdB_nearest = gdB.iloc[idx].drop(columns=\"geometry\").reset_index(drop=True)\n",
" gdf = pd.concat(\n",
" [\n",
" gdA.reset_index(drop=True),\n",
" gdB_nearest,\n",
" pd.Series(dist, name='dist')\n",
" ], \n",
" axis=1)\n",
"\n",
" return gdf\n",
"\n",
"#error bar is +- sqrt(250²*2)m ≈ +-353.553390593m\n",
"df_nearest = ckdnearest(nl_map_dst_points.to_crs(epsg=28992), df_geojson_clean.to_crs(epsg=28992))\n",
"\n",
"#ax = df_nearest.plot.bar(x='INWONER', y='dist')\n",
"df_nearest['dist'].plot(kind=\"hist\", bins=40, weights=df_nearest['INWONER'], figsize=(15,15))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"ename": "SyntaxError",
"evalue": "invalid syntax (745234620.py, line 3)",
"output_type": "error",
"traceback": [
"\u001b[0;36m Input \u001b[0;32mIn [13]\u001b[0;36m\u001b[0m\n\u001b[0;31m dfwims.\u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
]
}
],
"source": [
"df_wims = pd.read_csv(r'../data/wims.csv')\n",
"print(df_wims)\n",
"dfwims."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"#de eerste vijf rijen zijn meta-data en kunnen we gerust weghalen voor nu \n",
"df_gr_gr = pd.read_csv(r'../data/stemmen/01_Groningen/osv4-3_telling_gr2022_groningen.csv', skiprows=5, header=None, delimiter=';')\n",
"print(df_gr_gr)\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"De data is heel gecompliceerd, dus we moeten het bruikbaar maken in pandas. Met hoe de data eruit ziet in een CSV in libre office willen we een constructie maken van de verschillende onderdelen. We hebben de data van de stemlokalen apart al, die kunnen we dus negeren. We willen de data van de lijsten apart hebben. "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"TODO:mobiele stemlokalen eruit filteren die NIET 0,0 zijn\n",
"TODO:gooi dichtheid data van niet verkiezingsgemeente weg\n",
"TODO:datastandaard duitse verkiezingen stemlokaalafstand opzoeken"
]
}
],
"metadata": {
"interpreter": {
"hash": "3d1ccb9efb47efc90847d5703ae4e5663eea6663e461c6dfd943b50a653adc65"
},
"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.10.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Reference in a new issue Copy permalink