Pixels.cs

using System;
using System.Collections.Generic;

/// <summary>
/// Holds methods that operate on arrays of pixels held by images.
/// </summary>
public static class Pixels
{
    #region Generic

    /// <summary>
    /// Creates a checker pattern in an array of pixels.
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="a">first color/param>
    /// <param name="b">second color/param>
    /// <param name="cols">columns</param>
    /// <param name="rows">rows</param>
    /// <returns>checker pattern</returns>
    public static T[] Checker<T>(
        in T[] target,
        in int w, in int h,
        in T a, in T b,
        in int cols = 8, in int rows = 8)
    {
        const int limit = 2;
        int vCols = cols < 2 ? 2 : cols > w / limit ? w / limit : cols;
        int vRows = rows < 2 ? 2 : rows > h / limit ? h / limit : rows;
        int wch = w / vCols;
        int hchw = w * h / vRows;

        int trgLen = target.Length;
        for (int i = 0; i < trgLen; ++i)
        {
            // % 2 can be replaced by & 1 for even or odd.
            target[i] = (i % w / wch + i / hchw & 1) == 0 ? a : b;
        }
        return target;
    }

    /// <summary>
    /// Expands the image to the nearest greater power of two. Centers the
    /// source image within the expanded dimensions. Returns a tuple with the
    /// expanded pixels and the new width and height.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="uniform">use uniform size</param>
    /// <returns>tuple</returns>
    public static (T[] pixels, int w, int h) ExpandToPow2<T>(
        in T[] source, in int w, in int h, in bool uniform = false)
    {
        int wTrg;
        int hTrg;
        if (uniform)
        {
            int u = Utils.NextPowerOf2(Utils.Max(w, h));
            wTrg = u;
            hTrg = u;
        }
        else
        {
            wTrg = Utils.NextPowerOf2(w);
            hTrg = Utils.NextPowerOf2(h);
        }

        int trgLen = wTrg * hTrg;
        T[] target = new T[trgLen];
        int srcLen = source.Length;

        if (w == wTrg && h == hTrg)
        {
            System.Array.Copy(source, 0, target, 0, srcLen);
            return (pixels: target, w: wTrg, h: hTrg);
        }

        int xtl = (int)(wTrg * 0.5f + 0.5f) - (int)(w * 0.5f + 0.5f);
        int ytl = (int)(hTrg * 0.5f + 0.5f) - (int)(h * 0.5f + 0.5f);
        for (int i = 0; i < trgLen; ++i)
        {
            int xSrc = (i % wTrg) - xtl;
            int ySrc = (i / hTrg) - ytl;
            if (ySrc > -1 && ySrc < h && xSrc > -1 && xSrc < w)
            {
                target[i] = source[xSrc + ySrc * w];
            }
            else
            {
                target[i] = default;
            }
        }
        return (pixels: target, w: wTrg, h: hTrg);
    }

    /// <summary>
    /// Fills the pixels target array with a color.
    /// </summary>
    /// <param name="source">fill color</param>
    /// <param name="target">target pixels</param>
    /// <returns>filled array</returns>
    public static T[] Fill<T>(in T source, in T[] target)
    {
        int len = target.Length;
        for (int i = 0; i < len; ++i) { target[i] = source; }
        return target;
    }

    /// <summary>
    /// Transposes and flips the source array. The flipped transposition is
    /// stored in the target pixel array. The transposed image's width and
    /// height are swapped.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>transposed pixels</returns>
    public static T[] FlipAll<T>(
        in T[] source, in T[] target,
        in int w, in int h)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            int wn1 = w - 1;
            int hn1 = h - 1;
            for (int i = 0; i < srcLen; ++i)
            {
                // TODO: Multiply wn1 * h before the loop then just x * h within?
                target[(wn1 - i % w) * h + hn1 - i / w] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Flips the pixels source array horizontally, on the x axis, and stores
    /// the result in the target array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>flipped pixels</returns>
    public static T[] FlipX<T>(in T[] source, in T[] target, in int w, in int h)
    {
        int srcLen = source.Length;
        if (source == target)
        {
            int wd2 = w / 2;
            int wn1 = w - 1;
            int lenHalf = wd2 * h;
            for (int i = 0; i < lenHalf; ++i)
            {
                int x = i % wd2;
                int yw = w * (i / wd2);
                int idxSrc = x + yw;
                int idxTrg = yw + wn1 - x;
                (source[idxTrg], source[idxSrc]) = (source[idxSrc], source[idxTrg]);
            }
        }
        else if (srcLen == target.Length)
        {
            int wn1 = w - 1;
            for (int i = 0; i < srcLen; ++i)
            {
                target[i / w * w + wn1 - i % w] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Flips the pixels source array vertically, on the y axis, and stores the
    /// result in the target array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>flipped pixels</returns>
    public static T[] FlipY<T>(in T[] source, in T[] target, in int w, in int h)
    {
        int srcLen = source.Length;
        if (source == target)
        {
            int hd2 = h / 2;
            int hn1 = h - 1;
            int lenHalf = w * hd2;
            for (int i = 0; i < lenHalf; ++i)
            {
                int j = i % w + w * (hn1 - i / w);
                (source[j], source[i]) = (source[i], source[j]);
            }
        }
        else if (srcLen == target.Length)
        {
            int hn1 = h - 1;
            for (int i = 0; i < srcLen; ++i)
            {
                // TODO: Multiply hn1 * w before the loop then just y * w within?
                target[(hn1 - i / w) * w + i % w] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Gets a region defined by a top left and bottom right corner. May return
    /// an empty array if coordinates are invalid. Same as getting a cropped
    /// copy of an image.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="wSrc">source image width</param>
    /// <param name="hSrc">source image height</param>
    /// <param name="xtl">top left corner x</param>
    /// <param name="ytl">top left corner y</param>
    /// <param name="xbr">bottom right corner x</param>
    /// <param name="ybr">bottom right corner y</param>
    /// <returns>cropped region</returns>
    public static (T[] pixels, int w, int h, int x, int y) GetRegion<T>(
        in T[] source,
        in int wSrc, in int hSrc,
        in int xtl, in int ytl,
        in int xbr, in int ybr)
    {
        int xtlVrf = xtl < 0 ? 0 : xtl > wSrc - 1 ? wSrc - 1 : xtl;
        int ytlVrf = ytl < 0 ? 0 : ytl > hSrc - 1 ? hSrc - 1 : ytl;
        int xbrVrf = xbr < 0 ? 0 : xbr > wSrc - 1 ? wSrc - 1 : xbr;
        int ybrVrf = ybr < 0 ? 0 : ybr > hSrc - 1 ? hSrc - 1 : ybr;

        // Swap pixels if in wrong order.
        xtlVrf = Math.Min(xtlVrf, xbrVrf);
        ytlVrf = Math.Min(ytlVrf, ybrVrf);
        xbrVrf = Math.Max(xtlVrf, xbrVrf);
        ybrVrf = Math.Max(ytlVrf, ybrVrf);

        int wTrg = 1 + xbrVrf - xtlVrf;
        int hTrg = 1 + ybrVrf - ytlVrf;
        if (wTrg > 0 && hTrg > 0)
        {
            int trgLen = wTrg * hTrg;
            T[] target = new T[trgLen];
            for (int i = 0; i < trgLen; ++i)
            {
                int xTrg = i % wTrg;
                int yTrg = i / wTrg;
                int xSrc = xtlVrf + xTrg;
                int ySrc = ytlVrf + yTrg;
                int j = xSrc + ySrc * wSrc;
                target[i] = source[j];
            }
            return (pixels: target, w: wTrg, h: hTrg, x: xtlVrf, y: ytlVrf);
        }

        return (pixels: new T[0], w: 0, h: 0, x: 0, y: 0);
    }

    /// <summary>
    /// Mirrors, or reflects, pixels from a source image horizontally
    /// across a pivot. The pivot is expected to be in [-1, width + 1].
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="pivot">pivot</param>
    /// <param name="flip">flip reflection flag</param>
    /// <returns>mirrored pixels</returns>
    public static T[] MirrorX<T>(
        in T[] source, in T[] target, in int w,
        in int pivot = -1, in bool flip = false)
    {
        int trgLen = target.Length;
        int flipSign = flip ? 1 : -1;

        for (int i = 0; i < trgLen; ++i)
        {
            int cross = i % w - pivot;
            if (flipSign * cross < 0)
            {
                target[i] = source[i];
            }
            else
            {
                int pxOpp = pivot - cross;
                if (pxOpp > -1 && pxOpp < w)
                {
                    target[i] = source[i / w * w + pxOpp];
                }
                else
                {
                    target[i] = default;
                }
            }
        }

        return target;
    }

    /// <summary>
    /// Mirrors, or reflects, pixels from a source image vertically
    /// across a pivot. The pivot is expected to be in [-1, height + 1].
    /// Positive y points down to the bottom of the image.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="pivot">pivot</param>
    /// <param name="flip">flip reflection flag</param>
    /// <param name="target">target pixels</param>
    /// <returns>mirrored pixels</returns>
    public static T[] MirrorY<T>(
        in T[] source, in T[] target, in int w, in int h,
        in int pivot = -1, in bool flip = false)
    {
        int trgLen = target.Length;
        int flipSign = flip ? 1 : -1;

        for (int i = 0; i < trgLen; ++i)
        {
            int cross = i / w - pivot;
            if (flipSign * cross < 0)
            {
                target[i] = source[i];
            }
            else
            {
                int pyOpp = pivot - cross;
                if (pyOpp > -1 && pyOpp < h)
                {
                    target[i] = source[pyOpp * w + i % w];
                }
                else
                {
                    target[i] = default;
                }
            }
        }
        return target;
    }

    /// <summary>
    /// Rotates the source pixel array 90 degrees counter-clockwise. The
    /// rotation is stored in the target pixel array. The rotated image's
    /// width and height are swapped.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>rotated pixels</returns>
    public static T[] Rotate90<T>(in T[] source, in T[] target, in int w, in int h)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            int srcLennh = srcLen - h;
            for (int i = 0; i < srcLen; ++i)
            {
                target[srcLennh + i / w - i % w * h] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Rotates the source pixel array 180 degrees counter-clockwise. The
    /// rotation is stored in the target pixel array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>rotated pixels</returns>
    public static T[] Rotate180<T>(in T[] source, in T[] target)
    {
        int srcLen = source.Length;
        if (source == target)
        {
            int srcHalfLen = srcLen / 2;
            for (int i = 0, j = srcLen - 1; i < srcHalfLen; ++i, --j)
            {
                (source[j], source[i]) = (source[i], source[j]);
            }
        }
        else if (srcLen == target.Length)
        {
            for (int i = 0, j = srcLen - 1; i < srcLen; ++i, --j)
            {
                target[j] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Rotates the source pixel array 270 degrees counter-clockwise. The
    /// rotation is stored in the target pixel array. The rotated image's
    /// width and height are swapped.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>rotated pixels</returns>
    public static T[] Rotate270<T>(in T[] source, in T[] target, in int w, in int h)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            int hn1 = h - 1;
            for (int i = 0; i < srcLen; ++i)
            {
                target[i % w * h + hn1 - i / w] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Sets a region of pixels on the target array sfrom the source.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">source pixels</param>
    /// <param name="wSrc">source image width</param>
    /// <param name="hSrc">source image height</param>
    /// <param name="wTrg">source image width</param>
    /// <param name="hTrg">source image height</param>
    public static void SetRegion<T>(
        in T[] source, in T[] target,
        in int wSrc, in int hSrc,
        in int wTrg, in int hTrg)
    {
        Pixels.SetRegion(
            source, target,
            wSrc, hSrc,
            wTrg, hTrg,
            0, 0, wTrg - 1, hTrg - 1);
    }

    /// <summary>
    /// Sets a region defined by a "write" top left and bottom right corners on
    /// the target array from the source.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">source pixels</param>
    /// <param name="wSrc">source image width</param>
    /// <param name="hSrc">source image height</param>
    /// <param name="wTrg">source image width</param>
    /// <param name="hTrg">source image height</param>
    /// <param name="xtlWrite">top left corner x</param>
    /// <param name="ytlWrite">top left corner y</param>
    /// <param name="xbrWrite">bottom right corner x</param>
    /// <param name="ybrWrite">bottom right corner y</param>
    public static void SetRegion<T>(
        in T[] source, in T[] target,
        in int wSrc, in int hSrc,
        in int wTrg, in int hTrg,
        in int xtlWrite, in int ytlWrite,
        in int xbrWrite, in int ybrWrite)
    {
        Pixels.SetRegion(
            source, target,
            wSrc, hSrc,
            wTrg, hTrg,
            xtlWrite, ytlWrite, xbrWrite, ybrWrite,
            0, 0, wSrc - 1, hSrc - 1);
    }

    /// <summary>
    /// Sets a region defined by a "write" top left and bottom right
    /// corners on the target array from the source. The "read" corners
    /// specify a region from the source image.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">source pixels</param>
    /// <param name="wSrc">source image width</param>
    /// <param name="hSrc">source image height</param>
    /// <param name="wTrg">source image width</param>
    /// <param name="hTrg">source image height</param>
    /// <param name="xtlWrite">top left corner x</param>
    /// <param name="ytlWrite">top left corner y</param>
    /// <param name="xbrWrite">bottom right corner x</param>
    /// <param name="ybrWrite">bottom right corner y</param>
    /// <param name="xtlRead">top left corner x</param>
    /// <param name="ytlRead">top left corner y</param>
    /// <param name="xbrRead">bottom right corner x</param>
    /// <param name="ybrRead">bottom right corner y</param>
    public static void SetRegion<T>(
        in T[] source, in T[] target,
        in int wSrc, in int hSrc,
        in int wTrg, in int hTrg,
        in int xtlWrite, in int ytlWrite,
        in int xbrWrite, in int ybrWrite,
        in int xtlRead, in int ytlRead,
        in int xbrRead, in int ybrRead)
    {
        // TODO: Test

        int xtlWrVrf = xtlWrite < 0 ? 0 : xtlWrite > wTrg - 1 ? wTrg - 1 : xtlWrite;
        int ytlWrVrf = ytlWrite < 0 ? 0 : ytlWrite > hTrg - 1 ? hTrg - 1 : ytlWrite;
        int xbrWrVrf = xbrWrite < 0 ? 0 : xbrWrite > wTrg - 1 ? wTrg - 1 : xbrWrite;
        int ybrWrVrf = ybrWrite < 0 ? 0 : ybrWrite > hTrg - 1 ? hTrg - 1 : ybrWrite;

        xtlWrVrf = Math.Min(xtlWrVrf, xbrWrVrf);
        ytlWrVrf = Math.Min(ytlWrVrf, ybrWrVrf);
        xbrWrVrf = Math.Max(xtlWrVrf, xbrWrVrf);
        ybrWrVrf = Math.Max(ytlWrVrf, ybrWrVrf);

        int xtlRdVrf = xtlRead < 0 ? 0 : xtlRead > wSrc - 1 ? wSrc - 1 : xtlRead;
        int ytlRdVrf = ytlRead < 0 ? 0 : ytlRead > hSrc - 1 ? hSrc - 1 : ytlRead;
        int xbrRdVrf = xbrRead < 0 ? 0 : xbrRead > wSrc - 1 ? wSrc - 1 : xbrRead;
        int ybrRdVrf = ybrRead < 0 ? 0 : ybrRead > hSrc - 1 ? hSrc - 1 : ybrRead;

        xtlRdVrf = Math.Min(xtlRdVrf, xbrRdVrf);
        ytlRdVrf = Math.Min(ytlRdVrf, ybrRdVrf);
        xbrRdVrf = Math.Max(xtlRdVrf, xbrRdVrf);
        ybrRdVrf = Math.Max(ytlRdVrf, ybrRdVrf);

        int wReg = Math.Min(1 + xbrWrVrf - xtlWrVrf,
                            1 + xbrRdVrf - xtlRdVrf);
        int hReg = Math.Min(1 + ybrWrVrf - ytlWrVrf,
                            1 + ybrRdVrf - ytlRdVrf);

        if (wReg > 0 && hReg > 0)
        {
            int trgLen = wReg * hReg;
            for (int i = 0; i < trgLen; ++i)
            {
                int xReg = i % wReg;
                int yReg = i / wReg;
                int xSrc = xtlWrVrf + xReg;
                int ySrc = ytlWrVrf + yReg;
                int j = xSrc + ySrc * wTrg;

                int xTrg = xtlRdVrf + xReg;
                int yTrg = ytlRdVrf + yReg;
                int k = xTrg + yTrg * wSrc;

                target[j] = source[k];
            }
        }
    }

    /// <summary>
    /// Skews the pixels of a source image horizontally by an integer rise. The
    /// run specifies the number of pixels to skip on the x axis for each rise.
    /// If the rise or run are zero, copies the source array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="wSrc">image width/param>
    /// <param name="hSrc">image height</param>
    /// <param name="rise">rise, or step/param>
    /// <param name="run">run, or skip/param>
    /// <returns>tuple</returns>
    public static (T[] pixels, int w, int h) SkewXInteger<T>(
        in T[] source,
        in int wSrc, in int hSrc,
        in int rise = 1, in int run = 1)
    {
        int srcLen = source.Length;
        if (rise == 0 || run == 0)
        {
            T[] t0 = new T[srcLen];
            System.Array.Copy(source, 0, t0, 0, srcLen);
            return (pixels: t0, w: wSrc, h: hSrc);
        }

        int absRun = run < 0 ? -run : run;
        int sgnRise = run < 0 ? -rise : rise;
        int absRise = sgnRise < 0 ? -sgnRise : sgnRise;

        int hn1Run = (hSrc - 1) / absRun;
        int offset = sgnRise < 0 ? 0 : hn1Run * sgnRise;
        int wTrg = wSrc + hn1Run * absRise;
        T[] target = new T[wTrg * hSrc];

        for (int i = 0; i < srcLen; ++i)
        {
            int xSrc = i % wSrc;
            int ySrc = i / wSrc;
            int shift = sgnRise * (ySrc / absRun);
            int xTrg = xSrc + offset - shift;
            target[xTrg + ySrc * wTrg] = source[i];
        }

        return (pixels: target, w: wTrg, h: hSrc);
    }

    /// <summary>
    /// Skews the pixels of a source image vertically by an integer rise. The
    /// run specifies the number of pixels to skip on the y axis for each rise.
    /// If the rise or run are zero, copies the source array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="wSrc">image width/param>
    /// <param name="hSrc">image height</param>
    /// <param name="rise">rise, or step/param>
    /// <param name="run">run, or skip/param>
    /// <returns>tuple</returns>
    public static (T[] pixels, int w, int h) SkewYInteger<T>(
        in T[] source,
        in int wSrc, in int hSrc,
        in int rise = 1, in int run = 1)
    {
        int srcLen = source.Length;
        if (rise == 0 || run == 0)
        {
            T[] t0 = new T[srcLen];
            System.Array.Copy(source, 0, t0, 0, srcLen);
            return (pixels: t0, w: wSrc, h: hSrc);
        }

        int absRun = run < 0 ? -run : run;
        int sgnRise = run < 0 ? -rise : rise;
        int absRise = sgnRise < 0 ? -sgnRise : sgnRise;

        int wn1Run = (wSrc - 1) / absRun;
        int offset = sgnRise < 0 ? 0 : wn1Run * sgnRise;
        int hTrg = hSrc + wn1Run * absRise;
        T[] target = new T[wSrc * hTrg];

        for (int i = 0; i < srcLen; ++i)
        {
            int xSrc = i % wSrc;
            int ySrc = i / wSrc;
            int shift = sgnRise * (xSrc / absRun);
            int yTrg = ySrc + offset - shift;
            target[xSrc + yTrg * wSrc] = source[i];
        }

        return (pixels: target, w: wSrc, h: hTrg);
    }

    /// <summary>
    /// Transposes the source array. The transposition is stored in the target
    /// pixel array. The transposed image's width and height are swapped.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <returns>transposed pixels</returns>
    public static T[] Transpose<T>(
        in T[] source, in T[] target,
        in int w, in int h)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i % w * h + i / w] = source[i];
            }
        }
        return target;
    }

    /// <summary>
    /// Blits a source image's pixels onto a target image's pixels, using
    /// integer floor modulo to wrap the source image. The source image can be
    /// offset horizontally and/or vertically, creating the illusion of
    /// parallax.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="wSrc">source image width</param>
    /// <param name="hSrc">source image height</param>
    /// <param name="wTrg">target image width</param>
    /// <param name="dx">horizontal pixel offset</param>
    /// <param name="dy">vertical pixel offset</param>
    /// <returns>wrapped pixels</returns>
    public static T[] Wrap<T>(
        in T[] source, in T[] target,
        in int wSrc, in int hSrc, in int wTrg,
        in int dx, in int dy)
    {
        int trgLen = target.Length;
        for (int i = 0; i < trgLen; ++i)
        {
            int yMod = (i / wTrg + dy) % hSrc;
            if ((yMod ^ hSrc) < 0 && yMod != 0) { yMod += hSrc; }

            int xMod = (i % wTrg - dx) % wSrc;
            if ((xMod ^ wSrc) < 0 && xMod != 0) { xMod += wSrc; }

            target[i] = source[xMod + wSrc * yMod];
        }
        return target;
    }

    #endregion

    #region Specific

    /// <summary>
    /// Multiplies a source pixels array's alpha by the input alpha. If the
    /// alpha is less than zero, sets the target array to clear black. If the
    /// alpha is approximately one, copies the source pixels to the target.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="alpha">alpha</param>
    /// <returns>adjusted pixels</returns>
    public static Rgb[] AdjustAlpha(in Rgb[] source, in Rgb[] target, in float alpha)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            if (alpha <= 0.0f)
            {
                Pixels.Fill(Rgb.ClearBlack, target);
                return target;
            }

            if (Utils.Approx(alpha, 1.0f))
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                target[i] = new Rgb(c.R, c.G, c.B, c.Alpha * alpha);
            }
        }
        return target;
    }

    /// <summary>
    /// Adjusts a source pixels array's colors in LCH. Assigns the results
    /// to a target array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="adjust">adjustment</param>
    /// <returns>adjusted pixels</returns>
    public static Rgb[] AdjustSrLch(in Rgb[] source, in Rgb[] target, in Lch adjust)
    {
        // TODO: Make AdjustSrLab2 method?
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            if (adjust.Alpha == 0.0f &&
                adjust.L == 0.0f &&
                adjust.C == 0.0f &&
                adjust.H == 0.0f)
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            Dictionary<Rgb, Rgb> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                if (Rgb.Any(c) && !dict.ContainsKey(c))
                {
                    Lch lch = Rgb.StandardToSrLch(c);
                    dict.Add(c, Rgb.SrLchToStandard(lch + adjust));
                }
            }

            if (dict.Count > 0)
            {
                for (int j = 0; j < srcLen; ++j)
                {
                    Rgb c = source[j];
                    target[j] = dict.GetValueOrDefault(c, Rgb.ClearBlack);
                }
            }
            else
            {
                Pixels.Fill(Rgb.ClearBlack, target);
            }
        }
        return target;
    }

    /// <summary>
    /// Adjusts the contrast of colors from a source pixels array by a factor.
    /// Uses the SR LAB 2 color space. The adjustment factor is expected to be
    /// in [-1.0, 1.0].
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="fac">factor</param>
    /// <returns>adjusted pixels</returns>
    public static Rgb[] AdjustContrast(in Rgb[] source, in Rgb[] target, in float fac)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            float valAdjust = 1.0f + Utils.Clamp(fac, -1.0f, 1.0f);
            if (Utils.Approx(valAdjust, 1.0f))
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            Dictionary<Rgb, Rgb> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                if (Rgb.Any(c))
                {
                    Rgb key = Rgb.Opaque(c);
                    if (!dict.ContainsKey(key))
                    {
                        Lab lab = Rgb.StandardToSrLab2(key);
                        Lab labAdj = new(
                            l: (lab.L - 50.0f) * valAdjust + 50.0f,
                            a: lab.A,
                            b: lab.B,
                            alpha: lab.Alpha);
                        dict.Add(key, Rgb.SrLab2ToStandard(labAdj));
                    }
                }
            }

            if (dict.Count > 0)
            {
                Rgb clearBlack = Rgb.ClearBlack;
                for (int j = 0; j < srcLen; ++j)
                {
                    Rgb c = source[j];
                    target[j] = Rgb.CopyAlpha(
                        dict.GetValueOrDefault(Rgb.Opaque(c),
                            clearBlack), c);
                }
            }
            else
            {
                Pixels.Fill(Rgb.ClearBlack, target);
            }
        }
        return target;
    }

    /// <summary>
    /// Blends backdrop and overlay pixels using a function provided by the
    /// caller. Forms a union of the bounding area of the two inputs.
    /// Returns a tuple containing the blended pixels, the image's width and
    /// height, and the top left corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <param name="func">blend function</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) Blend(
        in Rgb[] apx, in Rgb[] bpx,
        in int aw, in int ah, in int bw, in int bh,
        in Func<Rgb, Rgb, Rgb> func)
    {
        int wLrg = aw > bw ? aw : bw;
        int hLrg = ah > bh ? ah : bh;

        // The 0.5 is to bias the rounding.
        float cx = wLrg * 0.5f + 0.5f;
        float cy = hLrg * 0.5f + 0.5f;

        int ax = aw == wLrg ? 0 : (int)(cx - aw * 0.5f);
        int ay = ah == hLrg ? 0 : (int)(cy - ah * 0.5f);
        int bx = bw == wLrg ? 0 : (int)(cx - bw * 0.5f);
        int by = bh == hLrg ? 0 : (int)(cy - bh * 0.5f);

        return Pixels.Blend(apx, bpx,
            aw, ah, bw, bh,
            ax, ay, bx, by, func);
    }

    /// <summary>
    /// Blends backdrop and overlay pixels using a function provided by the
    /// caller. Forms a union of the bounding area of the two inputs.
    /// Returns a tuple containing the blended pixels, the image's width and
    /// height, and the top left corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <param name="ax">under top left corner x</param>
    /// <param name="ay">under top left corner y</param>
    /// <param name="bx">over top left corner x</param>
    /// <param name="by">over top left corner y</param>
    /// <param name="func">blend function</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) Blend(
        in Rgb[] apx, in Rgb[] bpx,
        in int aw, in int ah, in int bw, in int bh,
        in int ax, in int ay, in int bx, in int by,
        in Func<Rgb, Rgb, Rgb> func)
    {
        int cx = ax < bx ? ax : bx;
        int cy = ay < by ? ay : by;
        int cw = aw > bw ? aw : bw;
        int ch = ah > bh ? ah : bh;
        int cLen = cw * ch;

        int axud = ax - cx;
        int ayud = ay - cy;
        int bxud = bx - cx;
        int byud = by - cy;

        Rgb[] target = new Rgb[cLen];
        for (int i = 0; i < cLen; ++i)
        {
            int x = i % cw;
            int y = i / cw;

            Rgb aClr = Rgb.ClearBlack;
            int axs = x - axud;
            int ays = y - ayud;
            if (ays > -1 && ays < ah && axs > -1 && axs < aw)
            {
                aClr = apx[axs + ays * aw];
            }

            Rgb bClr = Rgb.ClearBlack;
            int bxs = x - bxud;
            int bys = y - byud;
            if (bys > -1 && bys < bh && bxs > -1 && bxs < bw)
            {
                bClr = bpx[bxs + bys * bw];
            }

            target[i] = func(aClr, bClr);
        }

        return (pixels: target, w: cw, h: ch, x: cx, y: cy);
    }

    /// <summary>
    /// Blends backdrop and overlay pixels in SR LAB 2. Forms a union of the
    /// bounding area of the two inputs. Returns a tuple containing the blended
    /// pixels, the image's width and height, and the top left corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) BlendSrLab2(
        in Rgb[] apx, in Rgb[] bpx,
        in int aw, in int ah, in int bw, in int bh)
    {
        int wLrg = aw > bw ? aw : bw;
        int hLrg = ah > bh ? ah : bh;

        // The 0.5 is to bias the rounding.
        float cx = wLrg * 0.5f + 0.5f;
        float cy = hLrg * 0.5f + 0.5f;

        int ax = aw == wLrg ? 0 : (int)(cx - aw * 0.5f);
        int ay = ah == hLrg ? 0 : (int)(cy - ah * 0.5f);
        int bx = bw == wLrg ? 0 : (int)(cx - bw * 0.5f);
        int by = bh == hLrg ? 0 : (int)(cy - bh * 0.5f);

        return Pixels.BlendSrLab2(apx, bpx,
            aw, ah, bw, bh,
            ax, ay, bx, by);
    }

    /// <summary>
    /// Blends backdrop and overlay pixels in SR LAB 2. Forms a union of the
    /// bounding area of the two inputs. Returns a tuple containing the blended
    /// pixels, the image's width and height, and the top left corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <param name="ax">under top left corner x</param>
    /// <param name="ay">under top left corner y</param>
    /// <param name="bx">over top left corner x</param>
    /// <param name="by">over top left corner y</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) BlendSrLab2(
        in Rgb[] apx, in Rgb[] bpx,
        in int aw, in int ah, in int bw, in int bh,
        in int ax, in int ay, in int bx, in int by)
    {
        // Find the bottom right corner for a and b.
        int abrx = ax + aw - 1;
        int abry = ay + ah - 1;
        int bbrx = bx + bw - 1;
        int bbry = by + bh - 1;

        Dictionary<Rgb, Lab> dict = new()
        {
            { Rgb.ClearBlack, Lab.ClearBlack }
        };

        // Blending only necessary at the intersection of a and b.
        int dx = ax > bx ? ax : bx;
        int dy = ay > by ? ay : by;
        int dbrx = abrx < bbrx ? abrx : bbrx;
        int dbry = abry < bbry ? abry : bbry;
        int dw = 1 + dbrx - dx;
        int dh = 1 + dbry - dy;
        if (dw > 0 && dh > 0)
        {
            // Find the difference between the intersection top left corner
            // and the top left corners of a and b.
            int axid = ax - dx;
            int ayid = ay - dy;
            int bxid = bx - dx;
            int byid = by - dy;

            int dLen = dw * dh;
            for (int h = 0; h < dLen; ++h)
            {
                int x = h % dw;
                int y = h / dw;

                int axs = x - axid;
                int ays = y - ayid;
                Rgb aClr = apx[axs + ays * aw];
                if (!dict.ContainsKey(aClr))
                {
                    dict.Add(aClr, Rgb.StandardToSrLab2(aClr));
                }

                int bxs = x - bxid;
                int bys = y - byid;
                Rgb bClr = bpx[bxs + bys * bw];
                if (!dict.ContainsKey(bClr))
                {
                    dict.Add(bClr, Rgb.StandardToSrLab2(bClr));
                }
            }
        }

        // The result dimensions are the union of a and b.
        int cx = ax < bx ? ax : bx;
        int cy = ay < by ? ay : by;
        int cbrx = abrx > bbrx ? abrx : bbrx;
        int cbry = abry > bbry ? abry : bbry;
        int cw = 1 + cbrx - cx;
        int ch = 1 + cbry - cy;
        int cLen = cw * ch;

        // Find the difference between the union top left corner and the
        // top left corners of a and b.
        int axud = ax - cx;
        int ayud = ay - cy;
        int bxud = bx - cx;
        int byud = by - cy;

        Lab clearBlack = Lab.ClearBlack;
        Rgb[] target = new Rgb[cLen];
        for (int i = 0; i < cLen; ++i)
        {
            int x = i % cw;
            int y = i / ch;

            Rgb aClr = Rgb.ClearBlack;
            int axs = x - axud;
            int ays = y - ayud;
            if (ays > -1 && ays < ah && axs > -1 && axs < aw)
            {
                aClr = apx[axs + ays * aw];
            }

            Rgb bClr = Rgb.ClearBlack;
            int bxs = x - bxud;
            int bys = y - byud;
            if (bys > -1 && bys < bh && bxs > -1 && bxs < bw)
            {
                bClr = bpx[bxs + bys * bw];
            }

            float t = bClr.Alpha;
            if (t >= 1.0f) { target[i] = bClr; }
            else if (t <= 0.0f) { target[i] = aClr; }
            else
            {
                float v = aClr.Alpha;
                if (v <= 0.0f) { target[i] = bClr; }
                else
                {
                    float u = 1.0f - t;
                    float uv = u * v;
                    float tuv = t + uv;

                    Lab orig = dict.GetValueOrDefault(aClr, clearBlack);
                    Lab dest = dict.GetValueOrDefault(bClr, clearBlack);
                    Lab cLab = new(
                        l: u * orig.L + t * dest.L,
                        a: u * orig.A + t * dest.A,
                        b: u * orig.B + t * dest.B,
                        alpha: tuv);
                    target[i] = Rgb.SrLab2ToStandard(cLab);
                }
            }
        }
        return (pixels: target, w: cw, h: ch, x: cx, y: cy);
    }

    /// <summary>
    /// Blurs an array of pixels by averaging each color with its
    /// neighbors in 8 directions. The step determines the size of the kernel,
    /// where the minimum step of 1 will make a 3x3, 9 pixel kernel. Averages
    /// the color's SR LAB 2 representation.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>  
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="step">kernel step</param>
    /// <returns>blurred pixels</returns>
    public static Rgb[] BlurBoxSrLab2(
        in Rgb[] source, in Rgb[] target,
        in int w, in int h, in int step = 1)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            Dictionary<Rgb, Lab> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb srgb = source[i];
                if (!dict.ContainsKey(srgb))
                {
                    dict.Add(srgb, Rgb.StandardToSrLab2(srgb));
                }
            }

            int stepVal = step < 1 ? 1 : step;
            int wKrn = 1 + stepVal * 2;
            int krnLen = wKrn * wKrn;
            float denom = 1.0f / krnLen;

            for (int i = 0; i < srcLen; ++i)
            {
                // Subtract step to center the kernel in the inner for loop.
                int xSrc = i % w - stepVal;
                int ySrc = i / w - stepVal;
                Rgb srgb = source[i];

                float lSum = 0.0f;
                float aSum = 0.0f;
                float bSum = 0.0f;
                float tSum = 0.0f;

                for (int j = 0; j < krnLen; ++j)
                {
                    int xComp = xSrc + j % wKrn;
                    int yComp = ySrc + j / wKrn;
                    if (yComp > -1 && yComp < h && xComp > -1 && xComp < w)
                    {
                        Lab labNgbr = dict[source[xComp + yComp * w]];
                        lSum += labNgbr.L;
                        aSum += labNgbr.A;
                        bSum += labNgbr.B;
                        tSum += labNgbr.Alpha;
                    }
                    else
                    {
                        // When the kernel is out of bounds, sample the
                        // central color but do not tally alpha.
                        Lab labCtr = dict[srgb];
                        lSum += labCtr.L;
                        aSum += labCtr.A;
                        bSum += labCtr.B;
                    }
                }

                Lab labAvg = new(
                    l: lSum * denom,
                    a: aSum * denom,
                    b: bSum * denom,
                    alpha: tSum * denom);
                target[i] = Rgb.SrLab2ToStandard(labAvg);
            }
        }
        return target;
    }

    /// <summary>
    /// Clamps a pixel array to a lower and upper bounds.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="lb">lower bound</param>
    /// <param name="ub">upper bound</param>
    /// <returns>clamped array</returns>
    public static Rgb[] Clamp(
        in Rgb[] source, in Rgb[] target,
        in float lb = 0.0f, in float ub = 1.0f)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = Rgb.Clamp(source[i], lb, ub);
            }
        }
        return target;
    }

    /// <summary>
    /// Filters an array of colors according to a filter function. The lower
    /// bound and upper bound are both inclusive. Returns a dictionary wherein
    /// the key is a color that passed the filter and the value is an array of
    /// indices that reference that color's occurrence in the source array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="lb">lower bound</param>
    /// <param name="ub">upper bound</param>
    /// <param name="f">filter function</param>
    /// <returns>dictionary</returns>
    public static Dictionary<Rgb, List<int>> Filter(
        in Rgb[] source,
        in float lb, in float ub,
        in Func<Rgb, float, float, bool> f)
    {
        int srcLen = source.Length;
        HashSet<Rgb> visited = new(256);
        Dictionary<Rgb, bool> included = new(256);
        Dictionary<Rgb, List<int>> dict = new(256);
        for (int i = 0; i < srcLen; ++i)
        {
            Rgb c = source[i];
            bool include;
            if (visited.Contains(c))
            {
                include = included[c];
            }
            else
            {
                include = f(c, lb, ub);
                if (include) { dict.Add(c, new(32)); }
                visited.Add(c);
                included.Add(c, include);
            }
            if (include) { dict[c].Add(i); }
        }
        return dict;
    }

    /// <summary>
    /// Generates a linear gradient from an origin point to a destination point.
    /// The origin and destination should be in the range [-1.0, 1.0]. The
    /// scalar projection is clamped to [0.0, 1.0].
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <returns>linear gradient</returns>
    public static Rgb[] GradientLinear(
            in Rgb[] target, in int w, in int h, in ClrGradient cg)
    {
        return Pixels.GradientLinear(target, w, h, cg,
            (x, y, z) => Lab.Mix(x, y, z));
    }

    /// <summary>
    /// Generates a linear gradient from an origin point to a destination point.
    /// The origin and destination should be in the range [-1.0, 1.0]. The
    /// scalar projection is clamped to [0.0, 1.0].
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <param name="orig">origin</param>
    /// <param name="dest">destination</param>
    /// <returns>linear gradient</returns>
    public static Rgb[] GradientLinear(
            in Rgb[] target, in int w, in int h,
            in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing,
            in Vec2 orig, in Vec2 dest)
    {
        return Pixels.GradientLinear(target, w, h, cg, easing,
            orig.X, orig.Y, dest.X, dest.Y);
    }

    /// <summary>
    /// Generates a linear gradient from an origin point to a destination point.
    /// The origin and destination should be in the range [-1.0, 1.0]. The
    /// scalar projection is clamped to [0.0, 1.0].
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <param name="xOrig">x origin</param>
    /// <param name="yOrig">y origin</param>
    /// <param name="xDest">x destination</param>
    /// <param name="yDest">y destination</param>
    /// <returns>linear gradient</returns>
    public static Rgb[] GradientLinear(
        in Rgb[] target, in int w, in int h,
        in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing,
        in float xOrig = -1.0f, in float yOrig = 0.0f,
        in float xDest = 1.0f, in float yDest = 0.0f)
    {
        float bx = xOrig - xDest;
        float by = yOrig - yDest;

        float bbInv = 1.0f / MathF.Max(Utils.Epsilon, bx * bx + by * by);

        float bxbbinv = bx * bbInv;
        float bybbinv = by * bbInv;

        float xobx = xOrig * bxbbinv;
        float yoby = yOrig * bybbinv;
        float bxwInv2 = 2.0f / (w - 1.0f) * bxbbinv;
        float byhInv2 = 2.0f / (h - 1.0f) * bybbinv;

        Dictionary<float, Rgb> visited = new();
        int trgLen = target.Length;
        for (int i = 0; i < trgLen; ++i)
        {
            float fac = Utils.Clamp(xobx + bxbbinv - bxwInv2 * (i % w)
                                 + (yoby + byhInv2 * (i / w) - bybbinv));
            if (visited.ContainsKey(fac))
            {
                target[i] = visited[fac];
            }
            else
            {
                Lab lab = ClrGradient.Eval(cg, fac, easing);
                Rgb rgb = Rgb.SrLab2ToStandard(lab);
                target[i] = rgb;
                visited[fac] = rgb;
            }
        }
        return target;
    }

    /// <summary>
    /// Maps the colors of a source pixels array to those of a gradient using
    /// the source's luminance. Retains the original color's transparency.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="cg">color gradient</param>
    /// <returns>gradient map</returns>
    public static Rgb[] GradientMap(
        in Rgb[] source, in Rgb[] target, in ClrGradient cg)
    {
        return Pixels.GradientMap(source, target, cg,
            (x, y, z) => Lab.Mix(x, y, z));
    }

    /// <summary>
    /// Maps the colors of a source pixels array to those of a gradient using
    /// the source's luminance. Retains the original color's transparency.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <returns>gradient map</returns>
    public static Rgb[] GradientMap(
        in Rgb[] source, in Rgb[] target,
        in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            Dictionary<Rgb, Rgb> dict = new(256);
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                if (Rgb.Any(c))
                {
                    Rgb key = Rgb.Opaque(c);
                    if (!dict.ContainsKey(key))
                    {
                        float lum = Rgb.StandardLuminance(key);
                        float fac = lum <= 0.0031308f ?
                            lum * 12.92f :
                            MathF.Pow(lum, 1.0f / 2.4f) * 1.055f - 0.055f;
                        Lab lab = ClrGradient.Eval(cg, fac, easing);
                        dict.Add(key, Rgb.SrLab2ToStandard(lab));
                    }
                }
            }

            if (dict.Count > 0)
            {
                Rgb clearBlack = Rgb.ClearBlack;
                for (int i = 0; i < srcLen; ++i)
                {
                    Rgb c = source[i];
                    target[i] = Rgb.CopyAlpha(
                        dict.GetValueOrDefault(Rgb.Opaque(c),
                            clearBlack), c);
                }
            }
            else
            {
                Pixels.Fill(Rgb.ClearBlack, target);
            }
        }
        return target;
    }

    /// <summary>
    /// Generates a radial gradient from an origin point. The origin should be
    /// in the range [-1.0, 1.0]. Does not account for aspect ratio, so an image
    /// that isn't 1:1 will result in an ellipsoid.
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <param name="xOrig">x origin</param>
    /// <param name="yOrig">y origin</param>
    /// <param name="radius">radius</param>
    /// <returns>linear gradient</returns>
    public static Rgb[] GradientRadial(
        in Rgb[] target, in int w, in int h,
        in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing,
        in float xOrig, in float yOrig, in float radius)
    {
        float hInv2 = 2.0f / (h - 1.0f);
        float wInv2 = 2.0f / (w - 1.0f);

        float r2 = radius + radius;
        float rsqInv = 1.0f / MathF.Max(Utils.Epsilon, r2 * r2);

        float yon1 = yOrig - 1.0f;
        float xop1 = xOrig + 1.0f;

        Dictionary<float, Rgb> visited = new();
        int trgLen = target.Length;
        for (int i = 0; i < trgLen; ++i)
        {
            float ay = yon1 + hInv2 * (i / w);
            float ax = xop1 - wInv2 * (i % w);
            float fac = 1.0f - (ax * ax + ay * ay) * rsqInv;
            if (visited.ContainsKey(fac))
            {
                target[i] = visited[fac];
            }
            else
            {
                Lab lab = ClrGradient.Eval(cg, fac, easing);
                Rgb rgb = Rgb.SrLab2ToStandard(lab);
                target[i] = rgb;
                visited[fac] = rgb;
            }
        }
        return target;
    }

    /// <summary>
    /// Generates a sweep gradient, where the factor rotates on the z axis
    /// around an origin point. Best used with square images; for other aspect
    /// ratios, the origin should be adjusted accordingly.
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <returns>sweep gradient</returns>
    public static Rgb[] GradientSweep(
        in Rgb[] target, in int w, in int h, in ClrGradient cg)
    {
        return Pixels.GradientSweep(target, w, h, cg,
            (x, y, z) => Lab.Mix(x, y, z));
    }

    /// <summary>
    /// Generates a sweep gradient, where the factor rotates on the z axis
    /// around an origin point. Best used with square images; for other aspect
    /// ratios, the origin should be adjusted accordingly.
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <param name="radians">radians</param>
    /// <param name="orig">origin</param>
    /// <returns>sweep gradient</returns>
    public static Rgb[] GradientSweep(
        in Rgb[] target, in int w, in int h,
        in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing,
        in Vec2 orig, in float radians = 0.0f)
    {
        return Pixels.GradientSweep(target, w, h, cg, easing,
            orig.X, orig.Y, radians);
    }

    /// <summary>
    /// Generates a sweep gradient, where the factor rotates on the z axis
    /// around an origin point. Best used with square images; for other aspect
    /// ratios, the origin should be adjusted accordingly.
    /// </summary>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cg">color gradient</param>
    /// <param name="easing">easing function</param>
    /// <param name="radians">radians</param>
    /// <param name="xOrig">x origin</param>
    /// <param name="yOrig">y origin</param>
    /// <returns>sweep gradient</returns>
    public static Rgb[] GradientSweep(
        in Rgb[] target, in int w, in int h,
        in ClrGradient cg, in Func<Lab, Lab, float, Lab> easing,
        in float xOrig = 0.0f, in float yOrig = 0.0f,
        in float radians = 0.0f)
    {
        float aspect = w / (float)h;
        float wInv = aspect / (w - 1.0f);
        float hInv = 1.0f / (h - 1.0f);
        float xo = (xOrig * 0.5f + 0.5f) * aspect * 2.0f - 1.0f;
        float yo = yOrig;
        float rd = radians;

        Dictionary<float, Rgb> visited = new();
        int trgLen = target.Length;
        for (int i = 0; i < trgLen; ++i)
        {
            float xn = wInv * (i % w);
            float yn = hInv * (i / w);
            float fac = Utils.OneTau * Utils.WrapRadians(
                MathF.Atan2(
                    1.0f - (yn + yn + yo),
                    xn + xn - xo - 1.0f) - rd);
            if (visited.ContainsKey(fac))
            {
                target[i] = visited[fac];
            }
            else
            {
                Lab lab = ClrGradient.Eval(cg, fac, easing);
                Rgb rgb = Rgb.SrLab2ToStandard(lab);
                target[i] = rgb;
                visited[fac] = rgb;
            }
        }
        return target;
    }

    /// <summary>
    /// Inverts colors from a source pixels array in SR LAB 2.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>inverted pixels</returns>
    public static Rgb[] InvertSrLab2(
        in Rgb[] source, in Rgb[] target)
    {
        return Pixels.InvertSrLab2(source, target,
            true, true, true, false);
    }

    /// <summary>
    /// Inverts colors from a source pixels array in SR LAB 2.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="inv">inversion</param>
    /// <returns>inverted pixels</returns>
    public static Rgb[] InvertSrLab2(
        in Rgb[] source, in Rgb[] target, in Vec4 inv)
    {
        return Pixels.InvertSrLab2(source, target,
            inv.Z != 0.0f, inv.X != 0.0f, inv.Y != 0.0f, inv.W != 0.0f);
    }

    /// <summary>
    /// Inverts colors from a source pixels array in SR LAB 2.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="l">invert lightness</param>
    /// <param name="a">invert a</param>
    /// <param name="b">invert b</param>
    /// <param name="alpha">invert alpha</param>
    /// <returns>inverted pixels</returns>
    public static Rgb[] InvertSrLab2(
        in Rgb[] source, in Rgb[] target,
        in bool l, in bool a, in bool b, in bool alpha)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            if (!alpha && !l && !a && !b)
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            float aSign = a ? -1.0f : 1.0f;
            float bSign = b ? -1.0f : 1.0f;

            Dictionary<Rgb, Rgb> dict = new(256);
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                if (!dict.ContainsKey(c))
                {
                    Lab lab = Rgb.StandardToSrLab2(c);
                    Lab inv = new(
                        l: l ? 100.0f - lab.L : lab.L,
                        a: lab.A * aSign,
                        b: lab.B * bSign,
                        alpha: alpha ? 1.0f - lab.Alpha : lab.Alpha);
                    dict.Add(c, Rgb.SrLab2ToStandard(inv));
                }
            }

            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = dict[source[i]];
            }
        }
        return target;
    }

    /// <summary>
    /// Masks backdrop according to overlay pixels. Returns a tuple containing
    /// the blended pixels, the image's width and height, and the top left
    /// corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) Mask(
            in Rgb[] apx, in Rgb[] bpx,
            in int aw, in int ah, in int bw, in int bh)
    {
        int wLrg = aw > bw ? aw : bw;
        int hLrg = ah > bh ? ah : bh;

        // The 0.5 is to bias the rounding.
        float cx = wLrg * 0.5f + 0.5f;
        float cy = hLrg * 0.5f + 0.5f;

        int ax = aw == wLrg ? 0 : (int)(cx - aw * 0.5f);
        int ay = ah == hLrg ? 0 : (int)(cy - ah * 0.5f);
        int bx = bw == wLrg ? 0 : (int)(cx - bw * 0.5f);
        int by = bh == hLrg ? 0 : (int)(cy - bh * 0.5f);

        return Pixels.Mask(apx, bpx,
            aw, ah, bw, bh,
            ax, ay, bx, by);
    }

    /// <summary>
    /// Masks backdrop according to overlay pixels. Returns a tuple containing
    /// the blended pixels, the image's width and height, and the top left
    /// corner x and y.
    /// </summary>
    /// <param name="apx">under image pixels</param>
    /// <param name="bpx">over image pixels</param>
    /// <param name="aw">under width</param>
    /// <param name="ah">under height</param>
    /// <param name="bw">over width</param>
    /// <param name="bh">over height</param>
    /// <param name="ax">under top left corner x</param>
    /// <param name="ay">under top left corner y</param>
    /// <param name="bx">over top left corner x</param>
    /// <param name="by">over top left corner y</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) Mask(
        in Rgb[] apx, in Rgb[] bpx,
        in int aw, in int ah, in int bw, in int bh,
        in int ax, in int ay, in int bx, in int by)
    {
        // Find the bottom right corner for a and b.
        int abrx = ax + aw - 1;
        int abry = ay + ah - 1;
        int bbrx = bx + bw - 1;
        int bbry = by + bh - 1;

        // The result dimensions are the intersection of a and b.
        int cx = ax > bx ? ax : bx;
        int cy = ay > by ? ay : by;
        int cbrx = abrx < bbrx ? abrx : bbrx;
        int cbry = abry < bbry ? abry : bbry;
        int cw = 1 + cbrx - cx;
        int ch = 1 + cbry - cy;
        int cLen = cw * ch;

        // Find the difference between the target top left and the inputs.
        int axd = ax - cx;
        int ayd = ay - cy;
        int bxd = bx - cx;
        int byd = by - cy;

        Rgb[] target = new Rgb[cLen];
        for (int i = 0; i < cLen; ++i)
        {
            int x = i % cw;
            int y = i / ch;
            int bxs = x - bxd;
            int bys = y - byd;
            if (bys > -1 && bys < bh && bxs > -1 && bxs < bw)
            {
                Rgb bClr = bpx[bxs + bys * bw];
                if (bClr.Alpha > 0.0f)
                {
                    int axs = x - axd;
                    int ays = y - ayd;
                    if (ays > -1 && ays < ah && axs > -1 && axs < aw)
                    {
                        Rgb aClr = apx[axs + ays * aw];
                        if (aClr.Alpha > 0.0f)
                        {
                            target[i] = new Rgb(aClr.R, aClr.G, aClr.B,
                               aClr.Alpha * bClr.Alpha);
                        }
                        else { target[i] = Rgb.ClearBlack; }
                    }
                    else { target[i] = Rgb.ClearBlack; }
                }
                else { target[i] = Rgb.ClearBlack; }
            }
            else { target[i] = Rgb.ClearBlack; }
        }
        return (pixels: target, w: cw, h: ch, x: cx, y: cy);
    }

    /// <summary>
    /// Mirrors, or reflects, pixels from a source image across the axis
    /// described by an origin and destination. Coordinates are expected to be
    /// in the range [-1.0, 1.0]. Out-of-bounds pixels are omitted from the
    /// mirror.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="orig">origin</param>
    /// <param name="dest">destination</param>
    /// <param name="flip">flip mirror</param>
    /// <returns>palette</returns>    
    public static Rgb[] MirrorBilinear(
        in Rgb[] source, in Rgb[] target,
        in int w, in int h,
        in Vec2 orig, in Vec2 dest,
        in bool flip = false)
    {
        return Pixels.MirrorBilinear(
            source, target, w, h,
            orig.X, orig.Y, dest.X, dest.Y,
            flip);
    }

    /// <summary>
    /// Mirrors, or reflects, pixels from a source image across the axis
    /// described by an origin and destination. Coordinates are expected to be
    /// in the range [-1.0, 1.0]. Out-of-bounds pixels are omitted from the
    /// mirror.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="xOrig">x origin</param>
    /// <param name="yOrig">y origin</param>
    /// <param name="xDest">x destination</param>
    /// <param name="yDest">y destination</param>
    /// <param name="flip">flip mirror</param>
    /// <returns>palette</returns>    
    public static Rgb[] MirrorBilinear(
        in Rgb[] source, in Rgb[] target,
        in int w, in int h,
        in float xOrig = -1.0f, in float yOrig = 0.0f,
        in float xDest = 1.0f, in float yDest = 0.0f,
        in bool flip = false)
    {
        float wfn1 = w - 1.0f;
        float hfn1 = h - 1.0f;
        float wfp1Half = (w + 1.0f) * 0.5f;
        float hfp1Half = (h + 1.0f) * 0.5f;

        float ax = (xOrig + 1.0f) * wfp1Half - 0.5f;
        float bx = (xDest + 1.0f) * wfp1Half - 0.5f;
        float ay = (yOrig - 1.0f) * -hfp1Half - 0.5f;
        float by = (yDest - 1.0f) * -hfp1Half - 0.5f;

        float dx = bx - ax;
        float dy = by - ay;
        bool dxZero = Utils.Approx(dx, 0.0f, 0.5f);
        bool dyZero = Utils.Approx(dy, 0.0f, 0.5f);

        if (dxZero && dyZero)
        {
            System.Array.Copy(source, 0, target, 0,
                source.Length < target.Length ?
                source.Length : target.Length);
            return target;
        }

        if (dxZero)
        {
            return Pixels.MirrorX(source, target, w,
                Utils.Round(bx), flip ? ay > by : by > ay);
        }

        if (dyZero)
        {
            return Pixels.MirrorY(source, target, w, h,
                Utils.Round(by), flip ? bx > ax : ax > bx);
        }

        float dMagSqInv = 1.0f / (dx * dx + dy * dy);
        float flipSign = flip ? -1.0f : 1.0f;

        int trgLen = target.Length;
        for (int k = 0; k < trgLen; ++k)
        {
            float cy = k / w;
            float ey = cy - ay;
            float cx = k % w;
            float ex = cx - ax;

            float cross = ex * dy - ey * dx;
            if (flipSign * cross < 0.0f)
            {
                target[k] = source[k];
            }
            else
            {
                float t = (ex * dx + ey * dy) * dMagSqInv;
                float u = 1.0f - t;

                float pyProj = u * ay + t * by;
                float pyOpp = pyProj + pyProj - cy;
                float pxProj = u * ax + t * bx;
                float pxOpp = pxProj + pxProj - cx;

                // Default to omitting pixels that are out-of-bounds, rather than
                // wrapping with floor modulo or clamping.
                if (pyOpp >= 0.0f && pyOpp <= hfn1 &&
                    pxOpp >= 0.0f && pxOpp <= wfn1)
                {
                    target[k] = Pixels.SampleBilinear(source, pxOpp, pyOpp, w, h);
                }
                else
                {
                    target[k] = Rgb.ClearBlack;
                }
            }
        }
        return target;
    }

    /// <summary>
    /// Sets the alpha channel of all colors from the source array to 1.0.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>opaque array</returns>
    public static Rgb[] Opaque(in Rgb[] source, in Rgb[] target)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i) { target[i] = Rgb.Opaque(source[i]); }
        }
        return target;
    }

    /// <summary>
    /// Extracts a palette from a source pixels array using an octree in SR
    /// LAB 2. The size of the palette depends on the capacity of each node in
    /// the octree. Doesn't retain alpha component of image pixels. The
    /// threshold describes the minimum number of unique colors in the image
    /// beneath which it is preferable to not engage the octree.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="capacity">node capacity</param>
    /// <param name="threshold">threshold</param>
    /// <returns>palette</returns>
    public static Lab[] PaletteExtract(
        in Rgb[] source,
        in int capacity,
        in int threshold)
    {
        SortedSet<Rgb> uniqueRgbs = new();
        int srcLen = source.Length;
        for (int h = 0; h < srcLen; ++h)
        {
            Rgb rgbSrc = source[h];
            if (Rgb.Any(rgbSrc))
            {
                uniqueRgbs.Add(Rgb.Opaque(rgbSrc));
            }
        }

        // If under threshold, do not engage octree.
        int valThresh = threshold < 3 ? 3 : threshold;
        int uniquesLen = uniqueRgbs.Count;
        if (uniquesLen < valThresh)
        {
            // Prepend alpha mask.
            Lab[] under = new Lab[1 + uniquesLen];
            under[0] = Lab.ClearBlack;
            int i = 0;
            foreach (Rgb unique in uniqueRgbs)
            {
                under[++i] = Rgb.StandardToSrLab2(unique);
            }
            return under;
        }

        Octree oct = new(Bounds3.Lab, capacity);
        foreach (Rgb unique in uniqueRgbs)
        {
            Lab lab = Rgb.StandardToSrLab2(unique);
            oct.Insert(new(x: lab.A, y: lab.B, z: lab.L));
        }
        oct.Cull();

        List<Vec3> centers = new();
        Octree.CentersMean(oct, centers, false);
        int centersLen = centers.Count;

        Lab[] over = new Lab[1 + centersLen];
        over[0] = Lab.ClearBlack;
        for (int j = 0; j < centersLen; ++j)
        {
            Vec3 center = centers[j];
            Lab lab = new(l: center.Z, a: center.X, b: center.Y, alpha: 1.0f);
            over[1 + j] = lab;
        }
        return over;
    }

    /// <summary>
    /// Applies a palette to an array of pixels using an Octree to find the
    /// nearest match in Euclidean space. Retains the original color's
    /// transparency.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="palette">color palette</param>
    /// <param name="radius">query radius</param>
    /// <param name="capacity">node capacity</param>
    /// <returns>palette map</returns>
    public static Rgb[] PaletteMap(
        in Rgb[] source, in Rgb[] target, in Lab[] palette,
        in float radius = 128.0f, in int capacity = 16)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            Octree oct = new(Bounds3.Lab, capacity);
            oct.Subdivide(1, capacity);

            Dictionary<Vec3, Rgb> lookup = new();

            int palLen = palette.Length;
            for (int h = 0; h < palLen; ++h)
            {
                Lab lab = palette[h];
                if (Lab.Any(lab))
                {
                    Vec3 point = new(x: lab.A, y: lab.B, z: lab.L);
                    oct.Insert(point);
                    lookup.Add(point, Rgb.SrLab2ToStandard(lab));
                }
            }
            oct.Cull();

            Dictionary<Rgb, Rgb> dict = new();
            SortedList<float, Vec3> found = new(32);
            static float distFunc(Vec3 o, Vec3 d)
            {
                float da = d.X - o.X;
                float db = d.Y - o.Y;
                return MathF.Abs(d.Z - o.Z)
                    + MathF.Sqrt(da * da + db * db);
            }

            for (int i = 0; i < srcLen; ++i)
            {
                Rgb srgb = source[i];
                if (Rgb.Any(srgb))
                {
                    Rgb opaque = Rgb.Opaque(srgb);
                    if (dict.ContainsKey(opaque))
                    {
                        target[i] = Rgb.CopyAlpha(dict[opaque], srgb);
                    }
                    else
                    {
                        found.Clear();
                        Lab lab = Rgb.StandardToSrLab2(opaque);
                        Vec3 point = new(x: lab.A, y: lab.B, z: lab.L);
                        Octree.Query(oct, point, radius, distFunc, found);
                        if (found.Count > 0)
                        {
                            Vec3 nearest = found.Values[0];
                            if (lookup.ContainsKey(nearest))
                            {
                                Rgb match = lookup[nearest];
                                target[i] = Rgb.CopyAlpha(match, srgb);
                                dict.Add(opaque, match);
                            }
                            else { target[i] = Rgb.ClearBlack; }
                        }
                        else { target[i] = Rgb.ClearBlack; }
                    }
                }
                else { target[i] = Rgb.ClearBlack; }
            }
        }
        return target;
    }

    /// <summary>
    /// Multiplies the red, green and blue color channels of a pixel array by
    /// the alpha channel. If alpha is less than or equal to zero, returns
    /// clear black.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>premultiplied array</returns>
    public static Rgb[] Premul(in Rgb[] source, in Rgb[] target)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = Rgb.Premul(source[i]);
            }
        }
        return target;
    }

    /// <summary>
    /// Replaces the colors in the source pixels that approximately match the
    /// fromColor within a tolerance. Two colors are converted to LAB, then the
    /// distance between them is compared against the tolerance.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="fromColor">from color</param>
    /// <param name="toColor">to color</param>
    /// <returns>replaced array</returns>
    public static Rgb[] ReplaceColorApprox(
        in Rgb[] source, in Rgb[] target,
        in Rgb fromColor, in Rgb toColor,
        in float tolerance = 0.0f)
    {
        if (tolerance <= 0.0f)
        {
            return Pixels.ReplaceColorExact(source, target, fromColor, toColor);
        }

        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            Lab fromLab = Rgb.StandardToSrLab2(fromColor);
            Dictionary<Rgb, Rgb> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb srcClr = source[i];
                if (dict.ContainsKey(srcClr))
                {
                    target[i] = dict[source[i]];
                }
                else
                {
                    Lab srcLab = Rgb.StandardToSrLab2(srcClr);
                    Rgb trgClr = Lab.Dist(srcLab, fromLab) <= tolerance ?
                        toColor : srcClr;
                    dict.Add(srcClr, trgClr);
                    target[i] = trgClr;
                }
            }
        }
        return target;
    }

    /// <summary>
    /// Replaces the colors in the source pixels that equal the fromColor with
    /// the toColor.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="fromColor">from color</param>
    /// <param name="toColor">to color</param>
    /// <returns>replaced array</returns>
    public static Rgb[] ReplaceColorExact(
        in Rgb[] source, in Rgb[] target,
        in Rgb fromColor, in Rgb toColor)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                if (Rgb.EqSatArith(source[i], fromColor))
                {
                    source[i] = toColor;
                }
                else
                {
                    source[i] = target[i];
                }
            }
        }
        return target;
    }

    /// <summary>
    /// Resizes an image to the target width and height.
    /// Uses a bilinear filter.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="wSrc">image width</param>
    /// <param name="hSrc">image height</param>
    /// <param name="wTrg">rescaled width</param>
    /// <param name="hTrg">rescaled height</param>
    /// <returns>resized array</returns>
    public static Rgb[] ResizeBilinear(
        in Rgb[] source,
        in int wSrc, in int hSrc,
        in int wTrg, in int hTrg)
    {
        int srcLen = source.Length;
        if (wSrc == wTrg && hSrc == hTrg)
        {
            Rgb[] sameSize = new Rgb[srcLen];
            System.Array.Copy(source, 0, sameSize, 0, srcLen);
            return sameSize;
        }

        float wDenom = wTrg - 1.0f;
        float hDenom = hTrg - 1.0f;
        float tx = wDenom != 0.0f ? (wSrc - 1.0f) / wDenom : 0.0f;
        float ty = hDenom != 0.0f ? (hSrc - 1.0f) / hDenom : 0.0f;
        float ox = wDenom != 0.0f ? 0.0f : 0.5f;
        float oy = hDenom != 0.0f ? 0.0f : 0.5f;

        int trgLen = wTrg * hTrg;
        Rgb[] target = new Rgb[trgLen];
        for (int i = 0; i < trgLen; ++i)
        {
            target[i] = Pixels.SampleBilinear(source,
                (i % wTrg) * tx + ox, (i / wTrg) * ty + oy,
                wSrc, hSrc);
        }
        return target;
    }

    /// <summary>
    /// Rotates the pixels of a source image around the image center by an angle
    /// in radians. Where the angle is approximately 0, 90, 180 and 270 degrees,
    /// resorts to faster methods.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="radians">angle in radians</param>
    /// <returns>rotated pixels</returns>
    public static (Rgb[] pixels, int w, int h) RotateBilinear(
        in Rgb[] source, in int w, in int h, in float radians)
    {
        int srcLen = source.Length;
        int deg = Utils.RemFloor(Utils.Round(radians * Utils.RadToDeg), 360);

        switch (deg)
        {
            case 0:
                {
                    Rgb[] t0 = new Rgb[srcLen];
                    System.Array.Copy(source, 0, t0, 0, srcLen);
                    return (pixels: t0, w, h);
                }
            case 90:
                {
                    Rgb[] t90 = Pixels.Rotate90(source, new Rgb[srcLen], w, h);
                    return (pixels: t90, w: h, h: w);
                }
            case 180:
                {
                    Rgb[] t180 = Pixels.Rotate180(source, new Rgb[srcLen]);
                    return (pixels: t180, w, h);
                }
            case 270:
                {
                    Rgb[] t270 = Pixels.Rotate270(source, new Rgb[srcLen], w, h);
                    return (pixels: t270, w: h, h: w);
                }
            default:
                {
                    return Pixels.RotateBilinear(source, w, h,
                        MathF.Cos(radians), MathF.Sin(radians));
                }
        }
    }

    /// <summary>
    /// Rotates the pixels of a source image around the image center by an angle
    /// in radians. Assumes that the sine and cosine of the angle have already
    /// been calculated and simple cases (0, 90, 180, 270 degrees) have been
    /// filtered out.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="cosa">cosine of the angle</param>
    /// <param name="sina">sine of the angle</param>
    /// <returns>rotated pixels</returns>
    public static (Rgb[] pixels, int w, int h) RotateBilinear(
        in Rgb[] source, in int w, in int h,
        in float cosa, in float sina)
    {
        // Do not generalize to pass in filter as a function,
        // as nearest neighbor may need to calculate the center
        // and round from decimals to integers differently.

        float wSrcf = w;
        float hSrcf = h;
        float absCosa = MathF.Abs(cosa);
        float absSina = MathF.Abs(sina);

        int wTrg = (int)(0.5f + hSrcf * absSina + wSrcf * absCosa);
        int hTrg = (int)(0.5f + hSrcf * absCosa + wSrcf * absSina);
        float wTrgf = wTrg;
        float hTrgf = hTrg;

        float xSrcCenter = wSrcf * 0.5f;
        float ySrcCenter = hSrcf * 0.5f;
        float xTrgCenter = wTrgf * 0.5f;
        float yTrgCenter = hTrgf * 0.5f;

        int trgLen = wTrg * hTrg;
        Rgb[] target = new Rgb[trgLen];
        for (int i = 0; i < trgLen; ++i)
        {
            float ySgn = i / wTrg - yTrgCenter;
            float xSgn = i % wTrg - xTrgCenter;
            target[i] = Pixels.SampleBilinear(
                source,
                xSrcCenter + cosa * xSgn - sina * ySgn,
                ySrcCenter + cosa * ySgn + sina * xSgn,
                w, h);
        }

        return (pixels: target, w: wTrg, h: hTrg);
    }

    /// <summary>
    /// Internal helper function to sample a source image with a bilinear color
    /// mix. Clamps the result to [0.0, 1.0].
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="xSrc">x source</param>
    /// <param name="ySrc">y source</param>
    /// <param name="wSrc">image width</param>
    /// <param name="hSrc">image height</param>
    /// <returns>mixed color</returns>
    static Rgb SampleBilinear(
        in Rgb[] source, in float xSrc, in float ySrc,
        in int wSrc, in int hSrc)
    {
        // Find the truncation, floor and ceiling.
        int yi = (int)ySrc;
        int yf = (ySrc > 0.0f) ? yi : (ySrc < -0.0f) ? yi - 1 : 0;
        int yc = (ySrc > 0.0f) ? yi + 1 : (ySrc < -0.0f) ? yi : 0;

        bool yfInBounds = yf >= 0 && yf < hSrc;
        bool ycInBounds = yc >= 0 && yc < hSrc;

        int xi = (int)xSrc;
        int xf = (xSrc > 0.0f) ? xi : (xSrc < -0.0f) ? xi - 1 : 0;
        int xc = (xSrc > 0.0f) ? xi + 1 : (xSrc < -0.0f) ? xi : 0;

        bool xfInBounds = xf >= 0 && xf < wSrc;
        bool xcInBounds = xc >= 0 && xc < wSrc;

        // Pixel corners colors.
        Rgb c00 = yfInBounds && xfInBounds ? source[yf * wSrc + xf] : Rgb.ClearBlack;
        Rgb c10 = yfInBounds && xcInBounds ? source[yf * wSrc + xc] : Rgb.ClearBlack;
        Rgb c11 = ycInBounds && xcInBounds ? source[yc * wSrc + xc] : Rgb.ClearBlack;
        Rgb c01 = ycInBounds && xfInBounds ? source[yc * wSrc + xf] : Rgb.ClearBlack;

        float xErr = xSrc - xf;

        float a0 = 0.0f;
        float r0 = 0.0f;
        float g0 = 0.0f;
        float b0 = 0.0f;

        float a00 = c00.Alpha;
        float a10 = c10.Alpha;
        if (a00 > 0.0f || a10 > 0.0f)
        {
            float u = 1.0f - xErr;
            a0 = u * a00 + xErr * a10;
            if (a0 > 0.0f)
            {
                r0 = u * c00.R + xErr * c10.R;
                g0 = u * c00.G + xErr * c10.G;
                b0 = u * c00.B + xErr * c10.B;
            }
        }

        float a1 = 0.0f;
        float r1 = 0.0f;
        float g1 = 0.0f;
        float b1 = 0.0f;

        float a01 = c01.Alpha;
        float a11 = c11.Alpha;
        if (a01 > 0.0f || a11 > 0.0f)
        {
            float u = 1.0f - xErr;
            a1 = u * a01 + xErr * a11;
            if (a1 > 0.0f)
            {
                r1 = u * c01.R + xErr * c11.R;
                g1 = u * c01.G + xErr * c11.G;
                b1 = u * c01.B + xErr * c11.B;
            }
        }

        if (a0 > 0.0f || a1 > 0.0f)
        {
            float yErr = ySrc - yf;
            float u = 1.0f - yErr;
            float a2 = u * a0 + yErr * a1;
            if (a2 > 0.0f)
            {
                return Rgb.Clamp(new Rgb(
                    u * r0 + yErr * r1,
                    u * g0 + yErr * g1,
                    u * b0 + yErr * b1, a2));
            }
        }

        return Rgb.ClearBlack;
    }

    /// <summary>
    /// Scales the source pixels. The scalar is expected to be a non-zero
    /// positive value. Returns a tuple containing the pixels, the new width
    /// and height.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width/param>
    /// <param name="h">image height</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h) ScaleBilinear(
    in Rgb[] source, in int w, in int h, in float scalar)
    {
        float sVrf = MathF.Abs(scalar);
        int wTrg = Utils.Max(1, (int)(0.5f + w * sVrf));
        int hTrg = Utils.Max(1, (int)(0.5f + h * sVrf));
        return (pixels: Pixels.ResizeBilinear(source, w, h,
            wTrg, hTrg), w: wTrg, h: hTrg);
    }

    /// <summary>
    /// Scales the source pixels. The scalar is expected to be a non-zero
    /// positive value. Returns a tuple containing the pixels, the new width
    /// and height.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width/param>
    /// <param name="h">image height</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h) ScaleBilinear(
    in Rgb[] source, in int w, in int h, in Vec2 scalar)
    {
        Vec2 sVrf = Vec2.Abs(scalar);
        int wTrg = Utils.Max(1, (int)(0.5f + w * sVrf.X));
        int hTrg = Utils.Max(1, (int)(0.5f + h * sVrf.Y));
        return (pixels: Pixels.ResizeBilinear(source, w, h,
            wTrg, hTrg), w: wTrg, h: hTrg);
    }

    /// <summary>
    /// Skews the pixels of a source image horizontally. If the angle is
    /// approximately 0 degrees, copies the source array. If the angle is
    /// approximately 90 degrees, returns an empty array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width/param>
    /// <param name="h">image height</param>
    /// <param name="radians">radians</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h) SkewXBilinear(
        in Rgb[] source, in int w, in int h, in float radians)
    {
        int srcLen = source.Length;
        int deg = Utils.Round(radians * Utils.RadToDeg);
        int deg180 = Utils.RemFloor(deg, 180);

        switch (deg180)
        {
            case 0:
                {
                    Rgb[] t0 = new Rgb[srcLen];
                    System.Array.Copy(source, 0, t0, 0, srcLen);
                    return (pixels: t0, w, h);
                }
            case 26:
            case 27:
                {
                    return Pixels.SkewXInteger(source, w, h, 1, 2);
                }
            case 45:
                {
                    return Pixels.SkewXInteger(source, w, h, 1, 1);
                }
            case 63:
            case 64:
                {
                    return Pixels.SkewXInteger(source, w, h, 2, 1);
                }
            case 88:
            case 89:
            case 90:
            case 91:
            case 92:
                {
                    Rgb[] t90 = Pixels.Fill(Rgb.ClearBlack, new Rgb[srcLen]);
                    return (pixels: t90, w, h);
                }
            case 116:
            case 117:
                {
                    return Pixels.SkewXInteger(source, w, h, -2, 1);
                }
            case 135:
                {
                    return Pixels.SkewXInteger(source, w, h, -1, 1);
                }
            case 153:
            case 154:
                {
                    return Pixels.SkewXInteger(source, w, h, -1, 2);
                }
            default:
                {
                    float wSrcf = w;
                    float hSrcf = h;

                    float tana = MathF.Tan(radians);
                    int wTrg = (int)(0.5f + wSrcf + MathF.Abs(tana) * hSrcf);
                    float wTrgf = wTrg;
                    float yCenter = hSrcf * 0.5f;
                    float xDiff = (wSrcf - wTrgf) * 0.5f;

                    int trgLen = wTrg * h;
                    Rgb[] target = new Rgb[trgLen];
                    for (int i = 0; i < trgLen; ++i)
                    {
                        float yTrg = i / wTrg;
                        target[i] = Pixels.SampleBilinear(
                            source, xDiff + i % wTrg + tana
                           * (yTrg - yCenter), yTrg, w, h);
                    }

                    return (pixels: target, w: wTrg, h);
                }
        }
    }

    /// <summary>
    /// Skews the pixels of a source image vertically. If the angle is
    /// approximately 0 degrees, copies the source array. If the angle is
    /// approximately 90 degrees, returns an empty array.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width/param>
    /// <param name="h">image height</param>
    /// <param name="radians">radians</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h) SkewYBilinear(
    in Rgb[] source, in int w, in int h, in float radians)
    {
        int srcLen = source.Length;
        int deg = Utils.Round(radians * Utils.RadToDeg);
        int deg180 = Utils.RemFloor(deg, 180);

        switch (deg180)
        {
            case 0:
                {
                    Rgb[] t0 = new Rgb[srcLen];
                    System.Array.Copy(source, 0, t0, 0, srcLen);
                    return (pixels: t0, w, h);
                }
            case 26:
            case 27:
                {
                    return Pixels.SkewYInteger(source, w, h, 1, 2);
                }
            case 45:
                {
                    return Pixels.SkewYInteger(source, w, h, 1, 1);
                }
            case 63:
            case 64:
                {
                    return Pixels.SkewYInteger(source, w, h, 2, 1);
                }
            case 88:
            case 89:
            case 90:
            case 91:
            case 92:
                {
                    Rgb[] t90 = Pixels.Fill(Rgb.ClearBlack, new Rgb[srcLen]);
                    return (pixels: t90, w, h);
                }
            case 116:
            case 117:
                {
                    return Pixels.SkewYInteger(source, w, h, -2, 1);
                }
            case 135:
                {
                    return Pixels.SkewYInteger(source, w, h, -1, 1);
                }
            case 153:
            case 154:
                {
                    return Pixels.SkewYInteger(source, w, h, -1, 2);
                }
            default:
                {
                    float wSrcf = w;
                    float hSrcf = h;

                    float tana = MathF.Tan(radians);
                    int hTrg = (int)(0.5f + hSrcf + MathF.Abs(tana) * wSrcf);
                    float hTrgf = hTrg;
                    float xCenter = wSrcf * 0.5f;
                    float yDiff = (hSrcf - hTrgf) * 0.5f;

                    int trgLen = w * hTrg;
                    Rgb[] target = new Rgb[trgLen];
                    for (int i = 0; i < trgLen; ++i)
                    {
                        float xTrg = i % w;
                        target[i] = Pixels.SampleBilinear(
                            source, xTrg, yDiff + i / w + tana
                           * (xTrg - xCenter), w, h);
                    }
                    return (pixels: target, w, h);
                }
        }
    }

    /// <summary>
    /// Finds the minimum, maximum and mean lightness in a source pixels array.
    /// If factor is positive, stretches color to maximum lightness range in
    /// [0.0, 100.0]. If factor is negative, compresses color to mean. Assigns
    /// result to target array. The factor is expected to be in [-1.0, 1.0].
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="fac">factor</param>
    /// <returns>adjusted pixels</returns>
    public static Rgb[] StretchContrast(
        in Rgb[] source, in Rgb[] target, in float fac)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            float valFac = Utils.Clamp(fac, -1.0f, 1.0f);
            if (Utils.Approx(valFac, 0.0f))
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            // Beware MinValue in C# differs from MIN_VALUE in Java.
            // https://learn.microsoft.com/en-us/dotnet/api/system.single.minvalue
            float lumMin = float.MaxValue;
            float lumMax = float.MinValue;
            float lumSum = 0.0f;

            Dictionary<Rgb, Lab> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb c = source[i];
                if (Rgb.Any(c))
                {
                    Rgb key = Rgb.Opaque(c);
                    if (!dict.ContainsKey(key))
                    {
                        Lab lab = Rgb.StandardToSrLab2(key);
                        dict.Add(key, lab);

                        float lum = lab.L;
                        if (lum < lumMin) { lumMin = lum; }
                        if (lum > lumMax) { lumMax = lum; }
                        lumSum += lum;
                    }
                }
            }

            int dictLen = dict.Count;
            if (dictLen > 0)
            {
                float diff = MathF.Abs(lumMax - lumMin);
                if (diff > Utils.Epsilon)
                {
                    float t = MathF.Abs(valFac);
                    float u = 1.0f - t;
                    bool gtZero = valFac > 0.0f;
                    bool ltZero = valFac < -0.0f;

                    float lumAvg = lumSum / dictLen;
                    float tDenom = t * (100.0f / diff);
                    float lumMintDenom = lumMin * tDenom;
                    float tLumAvg = t * lumAvg;

                    Dictionary<Rgb, Rgb> stretched = new();
                    foreach (KeyValuePair<Rgb, Lab> kv in dict)
                    {
                        Lab srcLab = kv.Value;
                        Lab stretchedLab = srcLab;
                        if (gtZero)
                        {
                            stretchedLab = new(
                                l: u * srcLab.L
                                    + tDenom * srcLab.L - lumMintDenom,
                                a: srcLab.A,
                                b: srcLab.B,
                                alpha: srcLab.Alpha);
                        }
                        else if (ltZero)
                        {
                            stretchedLab = new(
                                l: u * stretchedLab.L + tLumAvg,
                                a: srcLab.A,
                                b: srcLab.B,
                                alpha: srcLab.Alpha);
                        }

                        stretched.Add(kv.Key, Rgb.SrLab2ToStandard(stretchedLab));
                    }

                    Rgb clearBlack = Rgb.ClearBlack;
                    for (int j = 0; j < srcLen; ++j)
                    {
                        Rgb c = source[j];
                        target[j] = Rgb.CopyAlpha(
                            stretched.GetValueOrDefault(Rgb.Opaque(c),
                                clearBlack), c);
                    }
                }
                else
                {
                    System.Array.Copy(source, 0, target, 0, srcLen);
                    return target;
                }
            }
            else
            {
                Pixels.Fill(Rgb.ClearBlack, target);
            }
        }
        return target;
    }

    /// <summary>
    /// Mixes an image with a color in SR LAB 2 according to a factor.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="tint">tint color</param>
    /// <param name="fac">factor</param>
    /// <param name="preserveLight">preserve lightness</param>
    /// <returns>tinted pixels</returns>
    public static Rgb[] TintSrLab2(
        in Rgb[] source,
        in Rgb[] target,
        in Lab tint,
        in float fac = 0.5f,
        in bool preserveLight = false)
    {
        return Pixels.TintSrLab2(
            source, target,
            tint, fac, preserveLight,
            (x) => 1.0f);
    }

    /// <summary>
    /// Mixes an image with a color in SR LAB 2 according to a factor.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="tint">tint color</param>
    /// <param name="fac">factor</param>
    /// <param name="preserveLight">preserve lightness</param>
    /// <param name="tonePreset">tone preset</param>
    /// <returns>tinted pixels</returns>
    public static Rgb[] TintSrLab2(
        in Rgb[] source,
        in Rgb[] target,
        in Lab tint,
        in float fac = 0.5f,
        in bool preserveLight = false,
        Tone tonePreset = Tone.Highlight
                        | Tone.Midtone
                        | Tone.Shadow)
    {
        return Pixels.TintSrLab2(
            source, target,
            tint, fac, preserveLight,
            Pixels.ToneFuncFromPreset(tonePreset));
    }

    /// <summary>
    /// Mixes an image with a color in SR LAB 2 according to a factor.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <param name="tint">tint color</param>
    /// <param name="fac">factor</param>
    /// <param name="preserveLight">preserve lightness</param>
    /// <param name="toneResponse">tone response function</param>
    /// <returns>tinted pixels</returns>
    public static Rgb[] TintSrLab2(
        in Rgb[] source,
        in Rgb[] target,
        in Lab tint,
        in float fac,
        in bool preserveLight,
        in Func<float, float> toneResponse)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            if (fac <= 0.0f)
            {
                System.Array.Copy(source, 0, target, 0, srcLen);
                return target;
            }

            float facVrf = fac >= 1.0f ? 1.0f : fac;
            float lTint = tint.L;
            float aTint = tint.A;
            float bTint = tint.B;

            Dictionary<Rgb, Rgb> dict = new();
            for (int i = 0; i < srcLen; ++i)
            {
                Rgb rgbSrc = source[i];
                Rgb rgbTrg;
                if (dict.ContainsKey(rgbSrc))
                {
                    rgbTrg = dict[rgbSrc];
                }
                else
                {
                    if (rgbSrc.Alpha <= 0.0f)
                    {
                        rgbTrg = Rgb.ClearBlack;
                    }
                    else
                    {
                        Lab labSrc = Rgb.StandardToSrLab2(rgbSrc);
                        float lightSrc = labSrc.L;
                        float lightFac = toneResponse(lightSrc * 0.01f);
                        float t = lightFac * facVrf;
                        float u = 1.0f - t;

                        Lab labTrg = new(
                        preserveLight ? lightSrc : u * lightSrc + t * lTint,
                        u * labSrc.A + t * aTint,
                        u * labSrc.B + t * bTint,
                        labSrc.Alpha);
                        rgbTrg = Rgb.SrLab2ToStandard(labTrg);
                    }

                    dict.Add(rgbSrc, rgbTrg);
                }
                target[i] = rgbTrg;
            }
        }
        return target;
    }

    /// <summary>
    /// Generates a tone response curve from a preset.
    /// </summary>
    /// <param name="preset">tone preset</param>
    /// <returns>tone function</returns>
    public static Func<float, float> ToneFuncFromPreset(Tone preset)
    {
        switch (preset)
        {
            case Tone.Shadow:
                return x =>
                {
                    if (x <= 0.0f) { return 1.0f; }
                    if (x >= 0.5f) { return 0.0f; }
                    float y = 1.0f - 2.0f * x;
                    return y * y * (3.0f - 2.0f * y);
                };

            case Tone.Midtone:
                return x =>
                {
                    if (x <= 0.0f || x >= 1.0f) { return 0.0f; }
                    float y = MathF.Abs(2.0f * x - 1.0f);
                    return 1.0f - y * y * (3.0f - 2.0f * y);
                };

            case Tone.Highlight:
                return x =>
                {
                    if (x <= 0.5f) { return 0.0f; }
                    if (x >= 1.0f) { return 1.0f; }
                    float y = 2.0f * x - 1.0f;
                    return y * y * (3.0f - 2.0f * y);
                };

            case Tone.Shadow | Tone.Midtone: // 1 | 2 = 3
                return x =>
                {
                    if (x <= 0.0f) { return 1.0f; }
                    if (x >= 0.75f) { return 0.0f; }
                    float y = 1.0f - Utils.FourThirds * x;
                    return y * y * (3.0f - 2.0f * y);
                };

            case Tone.Shadow | Tone.Highlight: // 1 | 4 = 5
                return x =>
                {
                    if (x <= 0.0f || x >= 1.0f) { return 1.0f; }
                    float y = MathF.Abs(1.0f - 2.0f * x);
                    return y * y * (3.0f - 2.0f * y);
                };

            case Tone.Midtone | Tone.Highlight: // 2 | 4 = 6
                return x =>
                {
                    if (x <= 0.25f) { return 0.0f; }
                    if (x >= 1.0f) { return 1.0f; }
                    float y = Utils.FourThirds * x - Utils.OneThird;
                    return y * y * (3.0f - 2.0f * y);
                };

            case (Tone)~0: // Due to how Unity handles flag enums.
            case Tone.Shadow | Tone.Midtone | Tone.Highlight: // 1 | 2 | 4 = 7
                return x =>
                {
                    if (x <= 0.0f) { return 0.0f; }
                    if (x >= 1.0f) { return 1.0f; }
                    return x * x * (3.0f - 2.0f * x);
                };

            default: return x => { return 1.0f; };
        }
    }

    /// <summary>
    /// Tone maps all colors in an array according to ACES curve.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>tone mapped array</returns>
    public static Rgb[] ToneMapAces(in Rgb[] source, in Rgb[] target)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = Rgb.ToneMapAcesStandard(source[i]);
            }
        }
        return target;
    }

    /// <summary>
    /// Tone maps all colors in an array according to Hable.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>tone mapped array</returns>
    public static Rgb[] ToneMapHable(in Rgb[] source, in Rgb[] target)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = Rgb.ToneMapHableStandard(source[i]);
            }
        }
        return target;
    }

    /// <summary>
    /// Removes excess transparent pixels from an array of pixels. Adapted from
    /// the implementation by Oleg Mikhailov: https://stackoverflow.com/a/36938923 .
    /// Returns a tuple with the trimmed image's pixels, width, height and
    /// top left corner.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="w">image width</param>
    /// <param name="h">image height</param>
    /// <param name="tlx">top left corner x</param>
    /// <param name="tly">top left corner y</param>
    /// <returns>tuple</returns>
    public static (Rgb[] pixels, int w, int h, int x, int y) TrimAlpha(
        in Rgb[] source, in int w, in int h, in int tlx = 0, in int tly = 0)
    {
        int srcLen = source.Length;

        if (w < 2 && h < 2)
        {
            Rgb[] trg0 = new Rgb[srcLen];
            System.Array.Copy(source, 0, trg0, 0, srcLen);
            return (pixels: trg0, w, h, x: tlx, y: tly);
        }

        int wn1 = w > 1 ? w - 1 : 0;
        int hn1 = h > 1 ? h - 1 : 0;

        int minRight = wn1;
        int minBottom = hn1;

        // Top search. y is outer loop, x is inner loop.
        int top = -1;
        bool goTop = true;
        while (goTop && top < hn1)
        {
            ++top;
            int wtop = w * top;
            int x = -1;
            while (goTop && x < wn1)
            {
                ++x;
                if (Rgb.Any(source[wtop + x]))
                {
                    minRight = x;
                    minBottom = top;
                    goTop = false;
                }
            }
        }

        // Left search. x is outer loop, y is inner loop.
        int left = -1;
        bool goLeft = true;
        while (goLeft && left < minRight)
        {
            ++left;
            int y = h;
            while (goLeft && y > top)
            {
                --y;
                if (Rgb.Any(source[y * w + left]))
                {
                    minBottom = y;
                    goLeft = false;
                }
            }
        }

        // Bottom search. y is outer loop, x is inner loop.
        int bottom = h;
        bool goBottom = true;
        while (goBottom && bottom > minBottom)
        {
            --bottom;
            int wbottom = w * bottom;
            int x = w;
            while (goBottom && x > left)
            {
                --x;
                if (Rgb.Any(source[wbottom + x]))
                {
                    minRight = x;
                    goBottom = false;
                }
            }
        }

        // Right search. x is outer loop, y is inner loop.
        int right = w;
        bool goRight = true;
        while (goRight && right > minRight)
        {
            --right;
            int y = bottom + 1;
            while (goRight && y > top)
            {
                --y;
                if (Rgb.Any(source[y * w + right]))
                {
                    goRight = false;
                }
            }
        }

        int wTrg = 1 + right - left;
        int hTrg = 1 + bottom - top;
        if (wTrg < 1 || hTrg < 1)
        {
            Rgb[] trg1 = new Rgb[srcLen];
            System.Array.Copy(source, 0, trg1, 0, srcLen);
            return (pixels: trg1, w, h, x: tlx, y: tly);
        }

        int trgLen = wTrg * hTrg;
        Rgb[] target = new Rgb[trgLen];
        for (int i = 0; i < trgLen; ++i)
        {
            target[i] = source[w * (top + i / wTrg) + left + i % wTrg];
        }
        return (pixels: target, w: wTrg, h: hTrg, x: tlx + left, y: tly + top);
    }

    /// <summary>
    /// Divides the red, green and blue color channels of a pixel array by the
    /// alpha channel. Reverses pre-multiplication. If alpha is less than or
    /// equal to zero, returns clear black.
    /// </summary>
    /// <param name="source">source pixels</param>
    /// <param name="target">target pixels</param>
    /// <returns>premultiplied array</returns>
    public static Rgb[] Unpremul(in Rgb[] source, in Rgb[] target)
    {
        int srcLen = source.Length;
        if (srcLen == target.Length)
        {
            for (int i = 0; i < srcLen; ++i)
            {
                target[i] = Rgb.Unpremul(source[i]);
            }
        }
        return target;
    }

    #endregion
}