{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "!curl --location --remote-name --remote-header-name \\\n",
    "    http://ftp.eea.europa.eu/www/eprtr/v17/E-PRTR_database_v17_xls.zip"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "!unzip E-PRTR_database_v17_xls.zip 'Pollutant releases.xlsx'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas\n",
    "\n",
    "pollutant_releases = pandas.read_excel('Pollutant releases.xlsx')\n",
    "pollutant_releases.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas\n",
    "\n",
    "# http://geopandas.org/gallery/create_geopandas_from_pandas.html\n",
    "pollutant_releases_geo = geopandas.GeoDataFrame(\n",
    "    pollutant_releases,\n",
    "    geometry=geopandas.points_from_xy(\n",
    "        pollutant_releases['Long'],\n",
    "        pollutant_releases['Lat'],\n",
    "    ),\n",
    ")\n",
    "pollutant_releases_geo.geometry.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "ax = world[(world.continent == 'Europe')].plot(\n",
    "    edgecolor='black',\n",
    "    color='white',\n",
    "    figsize=(20, 12),\n",
    ")\n",
    "ax.set_xlim(-28, 32)\n",
    "ax.set_ylim(32, 75)\n",
    "ax.set_title('E-PRTR_database_v17_xls.zip / Pollutant releases.xlsx')\n",
    "pollutant_releases_geo.plot(ax=ax, markersize=0.1);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "european_facilities = geopandas.GeoDataFrame(pollutant_releases_geo[\n",
    "    (pollutant_releases_geo.Long > -28)\n",
    "    & (pollutant_releases_geo.Long < 32)\n",
    "    & (pollutant_releases_geo.Lat > 32)\n",
    "    & (pollutant_releases_geo.Lat < 75)\n",
    "].groupby('FacilityID').first())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geoplot\n",
    "\n",
    "ax = geoplot.kdeplot(\n",
    "    european_facilities,\n",
    "    clip=world.geometry,\n",
    "    shade=True,\n",
    "    cmap='Reds',\n",
    "    projection=geoplot.crs.AlbersEqualArea(),\n",
    "    figsize=(24, 16),\n",
    ")\n",
    "ax.set_title('E-PRTR_database_v17_xls.zip / Pollutant releases.xlsx / kdeplot')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}