def test_hillshade_simple_transfer_function():
"""
Assert Simple Hillshade transfer function
"""
da_gaussian = xr.DataArray(data_gaussian)
da_gaussian_shade = geo.hillshade(da_gaussian)
assert da_gaussian_shade.mean() > 0
assert da_gaussian_shade[60, 60] > 0
def test_hillshade_simple_transfer_function():
"""
Assert Simple Hillshade transfer function
"""
da_gaussian = xr.DataArray(data_gaussian)
da_gaussian_shade = geo.hillshade(da_gaussian)
assert da_gaussian_shade.mean() > 0
assert da_gaussian_shade[60,60] > 0
def hillshade(ds, groupby=None, elev_var='topography__elevation', **kwargs):
elev = ds[elev_var]
if groupby is not None:
hshade = elev.groupby(groupby).apply(dgeo.hillshade,
shortcut=True,
**kwargs)
else:
hshade = dgeo.hillshade(elev, shortcut=True, **kwargs)
return (hshade.rename(dim_0='y', dim_1='x').assign_coords(x=ds.x, y=ds.y))
def hillshade(ds, groupby=None, elev_var='topography__elevation', **kwargs):
elev = ds[elev_var]
if groupby is not None:
# TODO: use shortcut=True
# https://github.com/holoviz/datashader/issues/871
hshade = elev.groupby(groupby).apply(dgeo.hillshade, shortcut=False,
**kwargs)
else:
hshade = dgeo.hillshade(elev, **kwargs)
return hshade
def _hillshade_one_run(ds, time_dim, elev_var, **kwargs):
elev = ds[elev_var]
if time_dim is not None:
# TODO: use shortcut=True
# https://github.com/holoviz/datashader/issues/871
hshade = elev.groupby(time_dim).apply(dgeo.hillshade,
shortcut=False,
**kwargs)
else:
hshade = dgeo.hillshade(elev, **kwargs)
return hshade