Merge d15932f into 16c5b75

dask_ndfilters/_conv.py

@@ -6,6 +6,34 @@
 import dask_ndfilters._utils as _utils
 
 
+@_utils._update_wrapper(scipy.ndimage.filters.convolve1d)
+def convolve1d(input,
+               weights,
+               axis=-1,
+               mode='reflect',
+               cval=0.0,
+               origin=0):
+    _utils._validate_axis(input.ndim, axis)
+    origin = _utils._get_origin(weights.shape, origin)
+    depth = _utils._get_depth(weights.shape, origin)
+    print(depth)
+    depth, boundary = _utils._get_depth_boundary(input.ndim, depth, "none")
+
+    result = input.map_overlap(
+        scipy.ndimage.filters.convolve1d,
+        depth=depth[0],
+        boundary=boundary[0],
+        dtype=input.dtype,
+        name=scipy.ndimage.filters.convolve1d.__name__,
+        weights=weights,
+        mode=mode,
+        cval=cval,
+        origin=origin[0]
+    )
+
+    return result
+
+
 @_utils._update_wrapper(scipy.ndimage.filters.convolve)
 def convolve(input,
              weights,

dask_ndfilters/_utils.py

@@ -65,6 +65,18 @@ def _updater(wrapper):
     return _updater
 
 
+def _validate_axis(ndim, axis):
+    if not isinstance(ndim, numbers.Integral):
+        raise TypeError("Expected integer value for `ndim`.")
+    if ndim <= 0:
+        raise ValueError("Expected positive value for `ndim`.")
+
+    if not isinstance(axis, numbers.Integral):
+        raise TypeError("Expected integer value for `axis`.")
+    if -ndim > axis or axis >= ndim:
+        raise ValueError("Need `axis` to be in range of `ndim`.")
+
+
 def _get_depth_boundary(ndim, depth, boundary=None):
     if not isinstance(ndim, numbers.Integral):
         raise TypeError("Expected integer value for `ndim`.")