def test_calculate_delta_height():
"""
Check that calculating change in elevation works.
"""
atl11_dataset: xr.Dataset = catalog.test_data.atl11_test_case.to_dask()
delta_height = calculate_delta(dataset=atl11_dataset,
oldcyclenum=3,
newcyclenum=4,
variable="h_corr")
assert isinstance(delta_height, xr.DataArray)
assert delta_height.shape == (1404, )
npt.assert_allclose(actual=delta_height.min().data, desired=-3.10612352)
npt.assert_allclose(actual=delta_height.mean().data, desired=-0.90124122)
npt.assert_allclose(actual=delta_height.max().data, desired=9.49908442)
def test_calculate_delta_time():
"""
Check that calculating change in time works.
"""
atl11_dataset: xr.Dataset = catalog.test_data.atl11_test_case.to_dask()
delta_time = calculate_delta(dataset=atl11_dataset,
oldcyclenum=3,
newcyclenum=4,
variable="delta_time")
assert isinstance(delta_time, xr.DataArray)
assert delta_time.shape == (1404, )
npt.assert_equal(actual=np.asarray(delta_time.min()),
desired=np.timedelta64(7846786703322903))
npt.assert_equal(actual=np.asarray(delta_time.mean()),
desired=np.timedelta64(7846786865357197)),
npt.assert_equal(actual=np.asarray(delta_time.max()),
desired=np.timedelta64(7846787022726588))
rasterize=True,
hover=True,
)
# %%
# %% [markdown]
# # Calculate Elevation Change (dh) over ICESAT-2 cycles!!
#
# Let's take a look at the change in elevation over one recent ICESat-2 cycle.
# From our loaded dataset (ds), we'll select Cycles 6 and 5,
# and subtract the height (h_corr) between them to get a height difference (dh).
# %%
dh: xr.DataArray = deepicedrain.calculate_delta(dataset=ds,
oldcyclenum=5,
newcyclenum=6,
variable="h_corr")
# %%
# Persist data in memory
dh = dh.persist()
# %%
delta_h: xr.Dataset = dh.dropna(dim="ref_pt").to_dataset(name="delta_height")
delta_h
# %%
df_dh: pd.DataFrame = delta_h.to_dataframe()
df_dh.head()
# %%