{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Nitrogen nutritional index - gridded for HiResIreland\n", "\n", "Soil chemical properties based on LUCAS topsoil data (Ballabio et al., 2019;\n", "European Commission, n.d.; Panagos et al., 2022; Panagos et al., 2012):\n", "https://esdac.jrc.ec.europa.eu/content/chemical-properties-european-scale-based-lucas-topsoil-data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "import geopandas as gpd\n", "import matplotlib.pyplot as plt\n", "import rioxarray as rxr\n", "from rasterstats import zonal_stats" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "DATA_DIR = os.path.join(\n", " \"data\",\n", " \"soil\",\n", " \"chemical-properties-european-scale-based-lucas-topsoil-data\",\n", ")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "DATA_FILE = os.path.join(DATA_DIR, \"IE_N.tif\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "data = rxr.open_rasterio(DATA_FILE, chunks=\"auto\", masked=True)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.DataArray (band: 1, y: 818, x: 835)>\n",
"dask.array<open_rasterio-fc87d89a5ce8db95ff48491a55ae131c<this-array>, shape=(1, 818, 835), dtype=float32, chunksize=(1, 818, 835), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * band (band) int64 1\n",
" * x (x) float64 2.918e+06 2.919e+06 ... 3.335e+06 3.335e+06\n",
" * y (y) float64 3.728e+06 3.728e+06 3.727e+06 ... 3.32e+06 3.32e+06\n",
" spatial_ref int64 0\n",
"Attributes:\n",
" AREA_OR_POINT: Area\n",
" DataType: Generic\n",
" RepresentationType: ATHEMATIC\n",
" STATISTICS_COVARIANCES: 0.9927304438511346\n",
" STATISTICS_MAXIMUM: 15.40726184845\n",
" STATISTICS_MEAN: 2.1118209123088\n",
" STATISTICS_MINIMUM: 0\n",
" STATISTICS_SKIPFACTORX: 1\n",
" STATISTICS_SKIPFACTORY: 1\n",
" STATISTICS_STDDEV: 0.99635859199945\n",
" scale_factor: 1.0\n",
" add_offset: 0.0| \n", " | rlon | \n", "rlat | \n", "sr | \n", "geometry | \n", "
|---|---|---|---|---|
| 0 | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "
| 1 | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "
| 2 | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "
| 3 | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "
| 4 | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "POLYGON ((421453.924 594208.112, 421445.655 59... | \n", "
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "count | \n", "mean | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "count | \n", "mean | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 6103 | \n", "POLYGON ((756735.711 882973.333, 756728.348 88... | \n", "1.365 | \n", "1.310 | \n", "2.021384 | \n", "0 | \n", "NaN | \n", "
| 6104 | \n", "POLYGON ((760717.627 840119.590, 760710.879 84... | \n", "1.400 | \n", "0.925 | \n", "2.126344 | \n", "0 | \n", "NaN | \n", "
| 6113 | \n", "POLYGON ((764608.869 847919.894, 764601.966 85... | \n", "1.435 | \n", "0.995 | \n", "2.685055 | \n", "0 | \n", "NaN | \n", "
| 6114 | \n", "POLYGON ((764601.966 851816.404, 764595.002 85... | \n", "1.435 | \n", "1.030 | \n", "2.685055 | \n", "0 | \n", "NaN | \n", "
| 6117 | \n", "POLYGON ((764580.889 863506.071, 764573.741 86... | \n", "1.435 | \n", "1.135 | \n", "2.497536 | \n", "0 | \n", "NaN | \n", "
290 rows × 6 columns
\n", "| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "count | \n", "mean | \n", "ni | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "count | \n", "mean | \n", "ni | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "