[求助]有关图象处理的问题

snowfox

dct变换的步骤：1。用i=imread读取图象的值
2。电位平移动。（将步骤一的各数据均减去128）
i=double(i)-128;
3.将系数分块为若干8*8的系数矩阵，（多出来的部分做补零处理）
4。对每个块做二维fdct，得到dct系数。用函数 dct2(i)
5。对每个8*8矩阵用公式c(k,l)=INT((X(k,l)/Q(k,l)+1/2)对矩阵进行量化。得到量化系数矩阵。
其中int为取整，k,l=0,1……7，X(k,l)为步骤4中的一个8*8矩阵，Q(k,l)为特定矩阵，其值为
Q =[ 16 11 10 16 24 40 51 61
12 12 14 19 26 58 60 55
14 13 16 24 40 57 69 56
14 17 22 29 51 87 80 62
18 22 37 56 68 109 103 77
24 35 55 6 81 104 113 92
49 64 78 87 103 121 120 101
72 92 95 98 112 100 103 99 ]
6.量化恢复矩阵，用公式X'(k,l)=c(k.l)*Q(k,l)
7.重建图象数据，（把若干个X'(k,l)组成一个新的数组i'，其i'的大小与i一样。
8.把新数据还原成图象。 (包括电位平移和精度转换）
注意点：不是直接调用函数blkproc，而是用for语句实现分块功能。
图象的重组也用for语句实现。
先本人的问题是如果实现图象的分块处理与合成，希望有高手能提供下代码！谢谢。