§3.2.3哈达玛变换 正反向哈达玛变换核(离散情况) 维正向哈达玛变换核的定义 g(x,u)=(1/N)(-1)2b(xbu;i=0,…,n-1,(N=2m 式中b(z)是z的二进制表达中的第k位(取0或1值); 维哈达玛变换为 H(u)=(1/N)∑f(x)(-1)2b(xb(;x=0,…,N-1 与沃尔什变换类似,哈达玛变换核组成的矩阵,对称且行列正交; 维哈达玛反向变换核的定义 h(x,u)=(-1)2b(xb0o;i=0,n-1 与正变换核相差常数1/N,算法可通用;
§3.2.3 哈达玛变换 • 一、正反向哈达玛变换核(离散情况) • 一维正向哈达玛变换核的定义: g(x,u)= (1/N)(-1) bi (x)bi(u);i=0,…,n-1,(N=2n) 式中bk (z)是z的二进制表达中的第k位(取0或1值); 一维哈达玛变换为 H(u)= (1/N) f(x)(-1) bi (x)bi(u) ; x=0,…,N-1 与沃尔什变换类似,哈达玛变换核组成的矩阵,对称且行列正交; 一维哈达玛反向变换核的定义: h(x,u)= (-1) bi (x)bi(u);i=0,…,n-1 与正变换核相差常数1/N,算法可通用;
§3.2.3哈达玛变换(续1) 维哈达玛反变换:(与正变换相差常数1/N) f(x)=∑H(u)(-1)2bi(xbu0;(u=0…,N-1) 二维哈达玛正变换核和反变换核 g(xyuV)=(1/N2)(-1)2lb(xbio)+by);(i=0,…,n-1) g1(, u) g2(y,V); h(x,y,u, v)=(1)2[bi(x)(u)+ bi (y)bi(v); (i=O,. n-1) h,(,u) h2(y,V) 维哈达玛正变换和反变换分别为: H(Uv) 1/N2)∑∑f(xy)(-1)2bi(xbu+bybi(y (Xy=0,…n-1) f(xy)=∑∑H(u,v)(-1)2(xb(u)+byb)1;(uV=0,…,n-1)
§3.2.3 哈达玛变换(续1) 一维哈达玛反变换:(与正变换相差常数1/N) f(x) = H(u)(-1) bi (x)bi(u) ; (u=0,…,N-1) 二、二维哈达玛正变换核和反变换核: g(x,y,u,v)= (1/N2)(-1) [bi (x)bi(u)+ bi (y)bi(v) ];(i=0,…,n-1) = g1(x,u)g2(y,v); h(x,y,u,v)= (-1) [bi (x)bi(u)+ bi (y)bi(v) ];(i=0,…,n-1) = h1(x,u)h2(y,v); 二维哈达玛正变换和反变换分别为: H(u,v)= ( 1/N2 ) f(x,y) (-1) [bi (x)bi(u)+ bi (y)bi(v) ]; (x,y=0,…,n-1) f(x,y) = H(u,v) (-1) [bi (x)bi(u)+ bi (y)bi(v) ]; (u,v=0,…,n-1)
832.3哈达玛变换(续2) N=8时,一维哈达玛变换核组成的矩阵如下(略去常数,仅用符 号表示+1/N、-1/N): +++++++ 符号变换次数07341625(竖直方向由+变-或由-变+) 可见,与沃尔什变换矩阵的区别行和列的次序不同; N=2时可混合使用
§3.2.3 哈达玛变换(续2) N=8时,一维哈达玛变换核组成的矩阵如下(略去常数,仅用符 号表示+1/N、-1/N): + + + + + + + + + - + - + - + - + + - - + + - - + - - + + - - + + + + + - - - - + - + - - + - + + + - - - - + + + - - + - + + - 符号变换次数0 7 3 4 1 6 2 5 (竖直方向由+变- 或 由-变+) 可见,与沃尔什变换矩阵的区别行和列的次序不同; N=2n时可混合使用
83.2.4离散余弦变换 DCT变换在图像处理中应用最广泛 维DCT的正向变换核和反向变换核 1D正变换核g(xu)=(2/N)12cos(2x+1)um/2N);x=0,…,N-1 =a(ucos(2X+ 1)uT/2N) 式中,a(u)=(1/N)12,当u=0,9(x0)=1/N12;a(u)=(2/N) 当u≠0; 1D反变换核h(Xu)=g(×,u) 由上述变换核有一维DCT的正反变换如下 C(u)=a(u)∑f(x)cos(2X+1)u/2N);X,u=0…N-1 f(x)=∑a(u)C(u)Cos(2X+1)uπ/2N);uX=0,…,N-1
§3.2.4 离散余弦变换 DCT变换在图像处理中应用最广泛; • 一、一维DCT的正向变换核和反向变换核 • 1D正变换核 g(x,u)=(2/N)1/2 cos((2x+1)u/2N);x=0,…,N-1 • =a(u)cos ((2x+1)u/2N); • 式中,a(u)= (1/N)1/2 ,当u=0,g(x,0)= 1/N1/2;a(u) = (2/N)1/2 , 当u0; • 1D反变换核 h(x,u) = g(x,u) • 由上述变换核有一维DCT的正反变换如下: • C(u)= a(u) ∑f(x)cos((2x+1)u/2N);x,u=0,…,N-1 • f(x) = ∑ a(u) C(u) cos((2x+1)u/2N);u,x =0,…,N-1
83.2.4离散余弦变换(续1) DCT计算 已知一维DCT C (u)=a(u Retz f(x) exp((2X+ 1)uT/2N)1 a(u)Re{exp(-juπ/2N)·∑f(x)exp(2Xuπ/2N)} 式中求和项就是2N个点的离散傅里叶变换 因此,DCT可以用FFT程序计算,反DCT也可以用反F计算 三、DCT的有关性质 DCT为实正交变换,即a=a*;a1=aT; DCT对高相关数据有能量集中的优势;能量集中在高频 正弦变换与余弦变换类似
§3.2.4 离散余弦变换(续1) • 二、DCT计算 • 已知一维DCT • C(u)= a(u)Re{ f(x)exp ((2x+1)u/2N) } • = a(u)Re{exp(-ju/2N)· f(x)exp(2xu/2N) } • 式中求和项就是2N个点的离散傅里叶变换, • 因此,DCT可以用FFT程序计算,反DCT也可以用反FFT计算; • 三、DCT的有关性质 • DCT为实正交变换,即a=a*;a -1 = a T ; • DCT对高相关数据有能量集中的优势;能量集中在高频。 • 正弦变换与余弦变换类似