{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Metocean gridded spatial skill\n", "\n", "Demonstrate plotting functionality etc by the SpatialSkill class" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import modelskill as ms" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\jem\\Source\\modelskill\\modelskill\\timeseries\\_track.py:134: UserWarning: Removed 22 duplicate timestamps with keep=first\n", " warnings.warn(\n" ] } ], "source": [ "fn = '../tests/testdata/NorthSeaHD_and_windspeed.dfsu'\n", "mr = ms.model_result(fn, name='HD', item=0)\n", "fn = '../tests/testdata/altimetry_NorthSea_20171027.csv'\n", "df = pd.read_csv(fn, index_col=0, parse_dates=True)\n", "o1 = ms.TrackObservation(df, item=2, name='alti')\n", "cc = ms.match(o1, mr)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ss = cc.gridded_skill()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n",
"Dimensions: (y: 5, x: 5)\n",
"Coordinates:\n",
" * y (y) float64 50.6 51.66 52.7 53.75 54.8\n",
" * x (x) float64 -0.436 1.543 3.517 5.492 7.466\n",
" observation <U4 'alti'\n",
"Data variables:\n",
" n (y, x) int32 3 17 0 0 0 0 49 9 0 ... 70 33 15 27 37 79 0 20 75\n",
" bias (y, x) float64 -0.02626 -0.188 nan nan ... nan -0.2475 -0.1156\n",
" rmse (y, x) float64 0.03291 0.2105 nan nan ... nan 0.2635 0.1475\n",
" urmse (y, x) float64 0.01984 0.09468 nan nan ... nan 0.09066 0.09174\n",
" mae (y, x) float64 0.02626 0.1923 nan nan ... nan 0.2475 0.1217\n",
" cc (y, x) float64 0.9973 0.996 nan nan ... nan 0.6293 0.9842\n",
" si (y, x) float64 0.04669 0.08355 nan nan ... nan 0.0844 0.1475\n",
" r2 (y, x) float64 -1.061 0.9003 nan nan ... nan -18.78 0.9078