Add GUFunc support for image to image filters

This is a very elegant method, suggested by @mrocklin, to enable application of a scalar ITK image filter on a NumPy on Dask array that has additional dimensions that corresponding to channels and time, for example.

More information:

- https://blog.dask.org/2019/08/09/image-itk
- https://docs.scipy.org/doc/numpy/reference/c-api.generalized-ufuncs.html
- https://blog.dask.org/2019/04/09/numba-stencil
- https://github.com/numpy/numpy/issues/14020