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