#!/usr/bin/env python
# coding: utf-8

# # Comparing Directional Data (e.g. wind direction)
# 
# Comparing directional data requires a little extra care and needs different axis on the plots. 

# In[1]:


import modelskill as ms
import mikeio


# In[2]:


fn = "../tests/testdata/wave_dir.dfs0"
ds = mikeio.read(fn)
ds


# In[3]:


cmp = ms.from_matched(ds, obs_item=3, mod_items=[2])
cmp


# In[4]:


cmp.plot.timeseries();


# We need to tell the plotter that the data is directional and that it should be plotted on a circular axis.

# In[5]:


cmp.plot.is_directional = True


# In[6]:


cmp.plot.timeseries();


# In[7]:


cmp.plot.kde();   # note: the KDE estimate is not directional! (yet)


# In[8]:


cmp.plot.scatter();  # note: regression line and Q-Q are not shown 


# In[ ]: