dct.php

<?php

function dct1D($in){
    $results = array();
    $N = count($in);
    for($k = 0; $k < $N; $k++){
        $sum = 0;
        for($n = 0; $n < $N; $n++){
             $sum += $in[$n] * cos($k * pi() * ($n + 0.5) / ($N));
        }
        $sum *= sqrt(2 / $N);
        if($k == 0){
            $sum *= 1 / sqrt(2);
        }
        $results[$k] = $sum;
    }
    return $results;
}

function optimizedImgDTC($img){
    $results = array();

    $N1 = imagesx($img);
    $N2 = imagesy($img);

    $rows = array();
    $row = array();
    for($j = 0; $j < $N2; $j++){
        for($i = 0; $i < $N1; $i++)
            $row[$i] = imagecolorat($img, $i, $j);
        $rows[$j] = dct1D($row);
    }

    for($i = 0; $i < $N1; $i++){
        for($j = 0; $j < $N2; $j++)
            $col[$j] = $rows[$j][$i];
        $results[$i] = dct1D($col);
    }
    return $results;
}