diff --git a/yt_common/yt_common/__init__.py b/yt_common/yt_common/__init__.py
index c024e69..3007da5 100644
--- a/yt_common/yt_common/__init__.py
+++ b/yt_common/yt_common/__init__.py
@@ -1 +1 @@
-from . import antialiasing, config, automation, logging, source  # noqa: F401
+from . import antialiasing, automation, config, deband, logging, source  # noqa: F401
diff --git a/yt_common/yt_common/deband.py b/yt_common/yt_common/deband.py
new file mode 100644
index 0000000..7effd72
--- /dev/null
+++ b/yt_common/yt_common/deband.py
@@ -0,0 +1,95 @@
+import vapoursynth as vs
+
+from enum import Enum
+from functools import partial
+from typing import Any, Callable, List, Optional, Tuple, Union
+
+import vsutil
+
+
+core = vs.core
+
+
+class XxpandModes(Enum):
+    SQUARE = 1
+    ELLIPSE = 2
+    LOSANGE = 3
+
+
+def mt_xxpand_multi(clip: vs.VideoNode, sw: int = 1, sh: Optional[int] = None,
+                    mode: XxpandModes = XxpandModes.SQUARE,
+                    planes: Union[int, List[int], None] = None, start: int = 0,
+                    m__imum: Callable[..., vs.VideoNode] = core.std.Maximum, **params: Any) -> List[vs.VideoNode]:
+    """
+    blame zastin for this
+    """
+    sh = vsutil.fallback(sh, sw)
+    assert clip.format is not None
+    planes = list(range(clip.format.num_planes)) if planes is None else [planes] if isinstance(planes, int) else planes
+
+    if mode == XxpandModes.ELLIPSE:
+        coordinates = [[1]*8, [0, 1, 0, 1, 1, 0, 1, 0], [0, 1, 0, 1, 1, 0, 1, 0]]
+    elif mode == XxpandModes.LOSANGE:
+        coordinates = [[0, 1, 0, 1, 1, 0, 1, 0]] * 3
+    else:
+        coordinates = [[1]*8] * 3
+
+    clips = [clip]
+
+    end = min(sw, sh) + start
+
+    for x in range(start, end):
+        clips += [m__imum(clips[-1], coordinates=coordinates[x % 3], planes=planes, **params)]
+
+    for x in range(end, end + sw - sh):
+        clips += [m__imum(clips[-1], coordinates=[0, 0, 0, 1, 1, 0, 0, 0], planes=planes, **params)]
+
+    for x in range(end, end + sh - sw):
+        clips += [m__imum(clips[-1], coordinates=[0, 1, 0, 0, 0, 0, 1, 0], planes=planes, **params)]
+
+    return clips
+
+
+def morpho_mask(clip: vs.VideoNode) -> Tuple[vs.VideoNode, vs.VideoNode]:
+    """
+    blame zastin for this
+    """
+    y = vsutil.get_y(clip)
+    assert y.format is not None
+
+    maxm = partial(mt_xxpand_multi, m__imum=core.std.Maximum)
+    minm = partial(mt_xxpand_multi, m__imum=core.std.Minimum)
+
+    ymax = maxm(y, sw=40, mode='ellipse')
+    ymin = minm(y, sw=40, mode='ellipse')
+
+    thr = round(0.0125 * ((1 << y.format.bits_per_sample) - 1))
+    ypw0 = y.std.Prewitt()
+    ypw = ypw0.std.Binarize(thr).rgvs.RemoveGrain(11)
+
+    rad = 3
+    thr = round(0.0098 * ((1 << y.format.bits_per_sample) - 1))
+    yrangesml = core.std.Expr([ymax[3], ymin[3]], 'x y - abs')
+    yrangesml = yrangesml.std.Binarize(thr).std.BoxBlur(0, 2, 1, 2, 1)
+
+    rad = 16
+    thr = round(0.0156 * ((1 << y.format.bits_per_sample) - 1))
+    yrangebig0 = core.std.Expr([ymax[rad], ymin[rad]], 'x y - abs')
+    yrangebig = yrangebig0.std.Binarize(thr)
+    yrangebig = minm(yrangebig,
+                     sw=rad * 3 // 4,
+                     threshold=(1 << y.format.bits_per_sample)//((rad * 3 // 4) + 1),
+                     mode='ellipse')[-1]
+    yrangebig = yrangebig.std.BoxBlur(0, rad//4, 1, rad//4, 1)
+
+    rad = 30
+    thr = round(0.0039 * ((1 << y.format.bits_per_sample) - 1))
+    ymph = core.std.Expr([y, maxm(ymin[rad], sw=rad, mode='ellipse')[rad],
+                          minm(ymax[rad], sw=rad, mode='ellipse')[rad]], 'x y - z x - max')
+    ymph = ymph.std.Binarize(round(0.00586 * ((1 << y.format.bits_per_sample) - 1)))
+    ymph = ymph.std.Minimum().std.Maximum()
+    ymph = ymph.std.BoxBlur(0, 4, 1, 4, 1)
+
+    grad_mask = core.std.Expr([ymph, yrangesml, ypw], expr="x y z max max")
+
+    return grad_mask, yrangebig