{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Track comparison for global wave model "
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from matplotlib_inline.backend_inline import set_matplotlib_formats\n",
"set_matplotlib_formats('png')\n",
"%matplotlib inline\n",
"\n",
"from mikeio.eum import EUMType, ItemInfo\n",
"from fmskill import ModelResult, TrackObservation, Connector"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Big data\n",
"Run the download.ipynb first"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"fn = '../data/SW_gwm_3a_extracted_2018.dfs0'\n",
"mr = ModelResult(fn, name='GWM')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"fn = '../data/altimetry_3a_2018_filter1.dfs0'\n",
"o1 = TrackObservation(fn, item=2, name='3a')"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"con = Connector(o1, mr[2])"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"cc = con.extract()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" n | \n",
" bias | \n",
" rmse | \n",
" urmse | \n",
" mae | \n",
" cc | \n",
" si | \n",
" r2 | \n",
"
\n",
" \n",
" | observation | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 3a | \n",
" 372399 | \n",
" -0.475266 | \n",
" 0.633188 | \n",
" 0.418388 | \n",
" 0.510794 | \n",
" 0.940363 | \n",
" 0.116908 | \n",
" 0.719939 | \n",
"
\n",
" \n",
"
\n",
"
\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" n | \n",
" bias | \n",
" rmse | \n",
" urmse | \n",
" mae | \n",
" cc | \n",
" si | \n",
" r2 | \n",
"
\n",
" \n",
" | observation | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 3a | \n",
" 9111762 | \n",
" -0.490248 | \n",
" 0.647898 | \n",
" 0.42359 | \n",
" 0.521119 | \n",
" 0.942401 | \n",
" 0.116337 | \n",
" 0.719179 | \n",
"
\n",
" \n",
"
\n",
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"ss.ds['bias'].attrs = dict(long_name=\"Bias of significant wave height, Hm0\",units=\"m\")\n",
"ss.ds['n'].attrs = dict(long_name=\"N of significant wave height\",units=\"-\")\n",
"fig, axes = plt.subplots(ncols=1, nrows=2, figsize = (8, 10))\n",
"ss.plot('n', ax=axes[0])\n",
"ss.plot('bias', ax=axes[1]);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Multiple bins - spatial skill for wave height \n",
"\n",
"Use all_df to obtain and df argument to pass customized data back to comparer."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"
\n",
" \n",
" \n",
" | \n",
" model | \n",
" observation | \n",
" x | \n",
" y | \n",
" mod_val | \n",
" obs_val | \n",
" xBin | \n",
" yBin | \n",
" val_cat | \n",
"
\n",
" \n",
" | datetime | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" | 2018-01-01 00:00:00 | \n",
" GWM | \n",
" 3a | \n",
" -33.706020 | \n",
" 23.181158 | \n",
" 2.292599 | \n",
" 2.611 | \n",
" -33.5 | \n",
" 23.5 | \n",
" Hm0[m]=[2, 5) | \n",
"
\n",
" \n",
" | 2018-01-01 00:00:01 | \n",
" GWM | \n",
" 3a | \n",
" -33.720741 | \n",
" 23.240074 | \n",
" 2.292612 | \n",
" 2.608 | \n",
" -33.5 | \n",
" 23.5 | \n",
" Hm0[m]=[2, 5) | \n",
"
\n",
" \n",
" | 2018-01-01 00:00:02 | \n",
" GWM | \n",
" 3a | \n",
" -33.735474 | \n",
" 23.298990 | \n",
" 2.292624 | \n",
" 2.518 | \n",
" -33.5 | \n",
" 23.5 | \n",
" Hm0[m]=[2, 5) | \n",
"
\n",
" \n",
" | 2018-01-01 00:00:03 | \n",
" GWM | \n",
" 3a | \n",
" -33.750214 | \n",
" 23.357904 | \n",
" 2.292637 | \n",
" 2.729 | \n",
" -33.5 | \n",
" 23.5 | \n",
" Hm0[m]=[2, 5) | \n",
"
\n",
" \n",
" | 2018-01-01 00:00:04 | \n",
" GWM | \n",
" 3a | \n",
" -33.764965 | \n",
" 23.416819 | \n",
" 2.292650 | \n",
" 2.593 | \n",
" -33.5 | \n",
" 23.5 | \n",
" Hm0[m]=[2, 5) | \n",
"
\n",
" \n",
"
\n",
"