{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Choropleth Plot of Sentinel 5P Data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Loading Data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import glob\n", "\n", "import geopandas\n", "import pandas\n", "import s5a\n", "\n", "# Read world map\n", "world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))\n", "\n", "# Read in all files\n", "data = []\n", "for filename in glob.glob('data/S5P_OFFL_L2__NO2____/*.nc'):\n", " data.append(\n", " s5a.filter_by_quality(\n", " s5a.load_ncfile(filename)))\n", "\n", "# Combine points\n", "data = pandas.concat(data, ignore_index=True)\n", "\n", "# Reduce points\n", "data = s5a.point_to_h3(data, resolution=5)\n", "data = s5a.aggregate_h3(data)\n", "data = s5a.h3_to_point(data)\n", "\n", "# Create geopandas dataframe\n", "geometry = geopandas.points_from_xy(data.longitude, data.latitude)\n", "data = geopandas.GeoDataFrame(data, geometry=geometry, crs={'init' :'epsg:4326'})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculate mean value per Country" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "scrolled": false }, "outputs": [ { "data": { "text/html": [ "
\n", " | pop_est | \n", "continent | \n", "name | \n", "iso_a3 | \n", "gdp_md_est | \n", "geometry | \n", "value | \n", "
---|---|---|---|---|---|---|---|
0 | \n", "920938 | \n", "Oceania | \n", "Fiji | \n", "FJI | \n", "8374.0 | \n", "MULTIPOLYGON (((180.00000 -16.06713, 180.00000... | \n", "0.000005 | \n", "
1 | \n", "53950935 | \n", "Africa | \n", "Tanzania | \n", "TZA | \n", "150600.0 | \n", "POLYGON ((33.90371 -0.95000, 34.07262 -1.05982... | \n", "0.000007 | \n", "
2 | \n", "603253 | \n", "Africa | \n", "W. Sahara | \n", "ESH | \n", "906.5 | \n", "POLYGON ((-8.66559 27.65643, -8.66512 27.58948... | \n", "0.000009 | \n", "
3 | \n", "35623680 | \n", "North America | \n", "Canada | \n", "CAN | \n", "1674000.0 | \n", "MULTIPOLYGON (((-122.84000 49.00000, -122.9742... | \n", "0.000002 | \n", "
4 | \n", "326625791 | \n", "North America | \n", "United States of America | \n", "USA | \n", "18560000.0 | \n", "MULTIPOLYGON (((-122.84000 49.00000, -120.0000... | \n", "0.000016 | \n", "
... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
172 | \n", "7111024 | \n", "Europe | \n", "Serbia | \n", "SRB | \n", "101800.0 | \n", "POLYGON ((18.82982 45.90887, 18.82984 45.90888... | \n", "0.000028 | \n", "
173 | \n", "642550 | \n", "Europe | \n", "Montenegro | \n", "MNE | \n", "10610.0 | \n", "POLYGON ((20.07070 42.58863, 19.80161 42.50009... | \n", "0.000010 | \n", "
174 | \n", "1895250 | \n", "Europe | \n", "Kosovo | \n", "-99 | \n", "18490.0 | \n", "POLYGON ((20.59025 41.85541, 20.52295 42.21787... | \n", "0.000023 | \n", "
175 | \n", "1218208 | \n", "North America | \n", "Trinidad and Tobago | \n", "TTO | \n", "43570.0 | \n", "POLYGON ((-61.68000 10.76000, -61.10500 10.890... | \n", "0.000014 | \n", "
176 | \n", "13026129 | \n", "Africa | \n", "S. Sudan | \n", "SSD | \n", "20880.0 | \n", "POLYGON ((30.83385 3.50917, 29.95350 4.17370, ... | \n", "0.000032 | \n", "
177 rows × 7 columns
\n", "