{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Soil water-holding capacity - gridded for HiResIreland\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", " | rlon | \n", "rlat | \n", "sr | \n", "ni | \n", "geometry | \n", "
|---|---|---|---|---|---|
| 0 | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((417558.169 590305.235, 417549.771 59... | \n", "
| 1 | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "
| 2 | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0.35 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "
| 3 | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "
| 4 | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "ni | \n", "mean | \n", "count | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.234, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "10.436507 | \n", "15 | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "15 | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "15 | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "15 | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "5.218254 | \n", "15 | \n", "
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "ni | \n", "mean | \n", "count | \n", "
|---|
| \n", " | geometry | \n", "rlon | \n", "rlat | \n", "sr | \n", "ni | \n", "whc | \n", "
|---|---|---|---|---|---|---|
| 0 | \n", "POLYGON ((417558.169 590305.234, 417549.771 59... | \n", "-1.680 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "10.436507 | \n", "
| 1 | \n", "POLYGON ((417549.771 594200.519, 417541.440 59... | \n", "-1.680 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "
| 2 | \n", "POLYGON ((421531.348 559152.004, 421522.478 56... | \n", "-1.645 | \n", "-1.595 | \n", "0.763591 | \n", "0.35 | \n", "0.000000 | \n", "
| 3 | \n", "POLYGON ((421462.259 590312.894, 421453.924 59... | \n", "-1.645 | \n", "-1.315 | \n", "0.917024 | \n", "0.35 | \n", "0.000000 | \n", "
| 4 | \n", "POLYGON ((421453.924 594208.111, 421445.655 59... | \n", "-1.645 | \n", "-1.280 | \n", "0.917024 | \n", "0.35 | \n", "5.218254 | \n", "