Skip to content

ISCC - Discrete Cosine Transform#

alg_dct(v) #

Discrete cosine transform.

See: nayuki.io.

Parameters:

Name Type Description Default
v Sequence[float]

Input vector for DCT calculation.

 required

Returns:

Type Description
List

DCT Transformed vector.
Source code in iscc_core\dct.py 
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
def alg_dct(v):
    # type: (Sequence[float]) -> List
    """
    Discrete cosine transform.

    See: [nayuki.io](https://www.nayuki.io/page/fast-discrete-cosine-transform-algorithms).

    :param Sequence[float] v: Input vector for DCT calculation.
    :return: DCT Transformed vector.
    :rtype: List
    """

    n = len(v)
    if n == 1:
        return list(v)
    elif n == 0 or n % 2 != 0:
        raise ValueError()
    else:
        half = n // 2
        alpha = [(v[i] + v[-(i + 1)]) for i in range(half)]
        beta = [
            (v[i] - v[-(i + 1)]) / (math.cos((i + 0.5) * math.pi / n) * 2.0) for i in range(half)
        ]
        alpha = alg_dct(alpha)
        beta = alg_dct(beta)
        result = []
        for i in range(half - 1):
            result.append(alpha[i])
            result.append(beta[i] + beta[i + 1])
        result.append(alpha[-1])
        result.append(beta[-1])
        return result