import mikeio
import matplotlib.pyplot as plt
ds = mikeio.read("../tests/testdata/waves.dfs2")
ds
<mikeio.Dataset> Dimensions: (time:3, y:31, x:31) Time: 2004-01-01 00:00:00 - 2004-01-03 00:00:00 (3 records) Items: 0: Sign. Wave Height <Significant wave height> (meter) 1: Peak Wave Period <Wave period> (second) 2: Mean Wave Direction <Mean Wave Direction> (degree)
ds.geometry.find_index(x=1400, y=0)
(array([28]), array([0]))
ds[0].plot()
plt.axvline(1400,c='k', linestyle='dashed');
ds_col = ds.isel(28, axis=2)
ds_col
<mikeio.Dataset> Dimensions: (time:3, y:31) Time: 2004-01-01 00:00:00 - 2004-01-03 00:00:00 (3 records) Items: 0: Sign. Wave Height <Significant wave height> (meter) 1: Peak Wave Period <Wave period> (second) 2: Mean Wave Direction <Mean Wave Direction> (degree)
ds_col[0].plot();
ds_col.to_dfs("waves_x1400.dfs1")
dsnew = mikeio.read("waves_x1400.dfs1")
dsnew
<mikeio.Dataset> Dimensions: (time:3, x:31) Time: 2004-01-01 00:00:00 - 2004-01-03 00:00:00 (3 records) Items: 0: Sign. Wave Height <Significant wave height> (meter) 1: Peak Wave Period <Wave period> (second) 2: Mean Wave Direction <Mean Wave Direction> (degree)
import os
os.remove("waves_x1400.dfs1")