{ "cells": [ { "cell_type": "markdown", "id": "80d20575-9c17-4a25-9d58-97e80601eadd", "metadata": {}, "source": [ "In this notebook we do spatial interpolation of CML rainfall data and compare IDW and Kriging, each with different parameters. Before we can do the spatial interpolation we first have to process the example CML data." ] }, { "cell_type": "code", "execution_count": 1, "id": "cbbdcae7-6419-4dec-9de4-5f64f5bcc32e", "metadata": {}, "outputs": [], "source": [ "import pycomlink as pycml\n", "import matplotlib.pyplot as plt\n", "import matplotlib as mpl\n", "import numpy as np\n", "import xarray as xr" ] }, { "cell_type": "markdown", "id": "04459b9b-6479-4d19-82ae-52212f933bef", "metadata": {}, "source": [ "# Load example CML data" ] }, { "cell_type": "code", "execution_count": 2, "id": "2f7d34b5-2559-499c-9f93-d518bdd62e61", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (time: 15840, cml_id: 500, channel_id: 2)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2018-05-10 ... 2018-05-20T23:59:00\n", " * cml_id (cml_id) object '0' '1' '2' '3' ... '497' '498' '499'\n", " length (cml_id) float64 6.179 5.673 7.52 ... 8.472 14.57 4.994\n", " site_a_latitude (cml_id) float64 58.26 58.09 58.19 ... 57.51 57.77 57.07\n", " site_a_longitude (cml_id) float64 1.388 1.637 1.359 ... 2.022 1.471 2.09\n", " site_b_latitude (cml_id) float64 58.25 58.13 58.21 ... 57.47 57.83 57.07\n", " site_b_longitude (cml_id) float64 1.304 1.59 1.461 ... 2.118 1.298 2.023\n", " * channel_id (channel_id) object 'channel_1' 'channel_2'\n", " frequency (cml_id, channel_id) float64 2.491e+10 ... 2.598e+10\n", " polarization (cml_id, channel_id) object 'V' 'V' 'H' ... 'V' 'V' 'V'\n", "Data variables:\n", " rsl (channel_id, cml_id, time) float32 ...\n", " tsl (channel_id, cml_id, time) float32 ...
<xarray.Dataset>\n", "Dimensions: (cml_id: 500, channel_id: 2, time: 264)\n", "Coordinates:\n", " * cml_id (cml_id) object '0' '1' '2' '3' ... '497' '498' '499'\n", " length (cml_id) float64 6.179 5.673 7.52 ... 8.472 14.57 4.994\n", " site_a_latitude (cml_id) float64 58.26 58.09 58.19 ... 57.51 57.77 57.07\n", " site_a_longitude (cml_id) float64 1.388 1.637 1.359 ... 2.022 1.471 2.09\n", " site_b_latitude (cml_id) float64 58.25 58.13 58.21 ... 57.47 57.83 57.07\n", " site_b_longitude (cml_id) float64 1.304 1.59 1.461 ... 2.118 1.298 2.023\n", " * channel_id (channel_id) object 'channel_1' 'channel_2'\n", " frequency (cml_id, channel_id) float64 2.491e+10 ... 2.598e+10\n", " polarization (cml_id, channel_id) object 'V' 'V' 'H' ... 'V' 'V' 'V'\n", " * time (time) datetime64[ns] 2018-05-10T01:00:00 ... 2018-05-21\n", "Data variables:\n", " R (channel_id, cml_id, time) float64 0.0 0.0 0.0 ... 0.0 0.0