DataArray - Dfs123 plotting¶
A DataArray with gridded data, can be plotted in many different ways.
In [1]:
import matplotlib.pyplot as plt
import mikeio
Dfs1¶
In [2]:
ds = mikeio.read("../tests/testdata/vu_tide_hourly.dfs1")
ds = ds.rename({"Tidal current component (geographic East)":"Tidal current u-comp"})
da = ds["Tidal current u-comp"]
da
Out[2]:
<mikeio.DataArray> name: Tidal current u-comp dims: (time:721, x:11) time: 2021-08-01 00:00:00 - 2021-08-31 00:00:00 (721 records) geometry: Grid1D (n=11, dx=0.09818)
In [3]:
da.geometry
Out[3]:
<mikeio.Grid1D> x: [0, 0.09818, ..., 0.9818] (nx=11, dx=0.09818)
In [4]:
da[0:10:2].plot();
In [5]:
da[0:10:2].plot.line()
plt.legend(da.time[0:10:2]);
In [6]:
# plot all points on line as time series
da.plot.timeseries();
In [7]:
# first 48 hours...
da[:49].plot.pcolormesh();
In [8]:
# single point on line as timeseries
da.sel(x=0.5).sel(time=slice("2021-08-01","2021-08-03")).plot();
In [9]:
# all data as histogram
da.plot.hist(bins=40);
Dfs2¶
In [10]:
da = mikeio.read("../tests/testdata/gebco_sound.dfs2")[0]
da
Out[10]:
<mikeio.DataArray> name: Elevation dims: (time:1, y:264, x:216) time: 2020-05-15 11:04:52 (time-invariant) geometry: Grid2D (ny=264, nx=216)
In [11]:
da.geometry
Out[11]:
<mikeio.Grid2D> x: [12.2, 12.21, ..., 13.1] (nx=216, dx=0.004167) y: [55.2, 55.21, ..., 56.3] (ny=264, dy=0.004167) projection: LONG/LAT
In [12]:
da.plot(figsize=(10,6));
It is also possible to customize the labels of the axes as well as the color bar, e.g. for localized adaption.
In [13]:
da.plot.contourf(figsize=(10,6), levels=4, label="Højde (m)")
plt.xlabel("Længdekreds (°)")
plt.ylabel("Breddekreds (°)")
Out[13]:
Text(0, 0.5, 'Breddekreds (°)')
In [14]:
ax = da.plot.contour(figsize=(8,8), cmap="plasma")
ax.set_xlim([12.5, 12.9])
ax.set_ylim([55.8, 56]);
In [15]:
da.plot.hist(bins=20);
Dfs3¶
In [16]:
fn = "../tests/testdata/test_dfs3.dfs3"
dfs = mikeio.open(fn)
dfs
Out[16]:
<mikeio.Dfs3> geometry: Grid3D(nz=34, ny=17, nx=21) items: 0: Temperature <Temperature> (degree Kelvin) time: 2 steps start time: 2020-12-30 00:00:00
In [17]:
dfs.geometry
Out[17]:
<mikeio.Grid3D> x: [0, 0.25, ..., 5] (nx=21, dx=0.25) y: [0, 0.25, ..., 4] (ny=17, dy=0.25) z: [0, 1, ..., 33] (nz=34, dz=1) origin: (10, 54), orientation: 0.000 projection: LONG/LAT
In [18]:
ds=dfs.read()
ds
Out[18]:
<mikeio.Dataset> dims: (time:2, z:34, y:17, x:21) time: 2020-12-30 00:00:00 - 2020-12-30 01:00:00 (2 records) geometry: Grid3D(nz=34, ny=17, nx=21) items: 0: Temperature <Temperature> (degree Kelvin)
In [19]:
ds.Temperature.plot();
In [20]:
ax = ds.Temperature[:,0,:,:].plot.contourf()
ax.grid()
In [21]:
ds.Temperature[:,0,:,0].plot();
In [22]:
ds=dfs.read(layers=0)
ds
Out[22]:
<mikeio.Dataset> dims: (time:2, y:17, x:21) time: 2020-12-30 00:00:00 - 2020-12-30 01:00:00 (2 records) geometry: Grid2D (ny=17, nx=21) items: 0: Temperature <Temperature> (degree Kelvin)
In [23]:
ds.Temperature.plot();
In [ ]: