{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Nitrogen nutritional index - gridded for MÉRA\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-9a5b296a8d3bf80d52b3a819e76bdfdf<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", " | x | \n", "y | \n", "sr | \n", "geometry | \n", "
|---|---|---|---|---|
| 0 | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "
| 1 | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "POLYGON ((424513.435 557931.478, 424058.108 56... | \n", "
| 2 | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "
| 3 | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "
| 4 | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "count | \n", "mean | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "count | \n", "mean | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 14483 | \n", "POLYGON ((759655.600 876544.928, 759200.942 87... | \n", "805000.0 | \n", "712500.0 | \n", "1.581323 | \n", "0 | \n", "NaN | \n", "
| 14484 | \n", "POLYGON ((759200.942 879003.081, 758746.267 88... | \n", "805000.0 | \n", "715000.0 | \n", "0.900000 | \n", "0 | \n", "NaN | \n", "
| 14485 | \n", "POLYGON ((758746.267 881461.207, 758291.577 88... | \n", "805000.0 | \n", "717500.0 | \n", "2.021384 | \n", "0 | \n", "NaN | \n", "
| 14486 | \n", "POLYGON ((758291.577 883919.305, 757836.870 88... | \n", "805000.0 | \n", "720000.0 | \n", "2.021384 | \n", "0 | \n", "NaN | \n", "
| 14487 | \n", "POLYGON ((765750.631 857333.286, 765296.086 85... | \n", "807500.0 | \n", "692500.0 | \n", "2.497536 | \n", "0 | \n", "NaN | \n", "
729 rows × 6 columns
\n", "| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "count | \n", "mean | \n", "ni | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "NaN | \n", "
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "count | \n", "mean | \n", "ni | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0 | \n", "NaN | \n", "0.35 | \n", "