{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Soil water-holding capacity - gridded for MÉRA\n", "\n", "European soil database derived data - total available water content (TAWC)\n", "for the topsoil [mm] (European Commission, n.d.; Hiederer, 2013a;\n", "Hiederer, 2013b):\n", "https://esdac.jrc.ec.europa.eu/content/european-soil-database-derived-data" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "import matplotlib.pyplot as plt\n", "import geopandas as gpd\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(\"data\", \"soil\", \"european-soil-database-derived-data\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "DATA_FILE = os.path.join(DATA_DIR, \"IE_TAWC.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: 409, x: 418)>\n",
"dask.array<open_rasterio-f60863f961dba41a45774b674cfc20e4<this-array>, shape=(1, 409, 418), dtype=float32, chunksize=(1, 409, 418), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * band (band) int64 1\n",
" * x (x) float64 2.918e+06 2.92e+06 2.92e+06 ... 3.334e+06 3.336e+06\n",
" * y (y) float64 3.728e+06 3.728e+06 ... 3.322e+06 3.32e+06\n",
" spatial_ref int64 0\n",
"Attributes:\n",
" AREA_OR_POINT: Area\n",
" units: milimeter\n",
" scale_factor: 1.0\n",
" add_offset: 0.0<xarray.DataArray (band: 1, y: 409, x: 418)>\n",
"dask.array<open_rasterio-f60863f961dba41a45774b674cfc20e4<this-array>, shape=(1, 409, 418), dtype=float32, chunksize=(1, 409, 418), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * band (band) int64 1\n",
" * x (x) float64 2.918e+06 2.92e+06 2.92e+06 ... 3.334e+06 3.336e+06\n",
" * y (y) float64 3.728e+06 3.728e+06 ... 3.322e+06 3.32e+06\n",
" spatial_ref int64 0\n",
"Attributes:\n",
" AREA_OR_POINT: Area\n",
" units: milimeter\n",
" scale_factor: 1.0\n",
" add_offset: 0.0| \n", " | x | \n", "y | \n", "sr | \n", "ni | \n", "geometry | \n", "
|---|---|---|---|---|---|
| 0 | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "
| 1 | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0.35 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "
| 2 | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0.35 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "
| 3 | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "
| 4 | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "ni | \n", "mean | \n", "count | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "8 | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "6 | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "8 | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0.35 | \n", "8.697089 | \n", "9 | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "6 | \n", "
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "ni | \n", "mean | \n", "count | \n", "
|---|
| \n", " | geometry | \n", "x | \n", "y | \n", "sr | \n", "ni | \n", "whc | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((415227.037 594343.036, 414771.907 59... | \n", "415000.0 | \n", "497500.0 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "
| 1 | \n", "POLYGON ((424513.435 557931.477, 424058.108 56... | \n", "417500.0 | \n", "460000.0 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "
| 2 | \n", "POLYGON ((424058.108 560389.011, 423602.795 56... | \n", "417500.0 | \n", "462500.0 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "
| 3 | \n", "POLYGON ((418595.321 589882.162, 418140.190 59... | \n", "417500.0 | \n", "492500.0 | \n", "0.917024 | \n", "0.35 | \n", "8.697089 | \n", "
| 4 | \n", "POLYGON ((418140.190 592340.144, 417685.076 59... | \n", "417500.0 | \n", "495000.0 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "