数字图像变换技术及MATLAB实现2

1、3.3.4 沃尔什函数的性质沃尔什函数的性质 沃尔什函数有如下一些主要性质：沃尔什函数有如下一些主要性质： （）（） 在区间在区间0,1内有下式成立内有下式成立 1), 0(10dttwal(4119)2 , 1 0),(10idttiwal(4120)3 ,2, 1 ,0 1),(2itiwal(4121)walt( , )0（）（） 在区间在区间0,1的第一小段时间内（通常称的第一小段时间内（通常称为时隙）沃尔什函数总是取。为时隙）沃尔什函数总是取。（）（） 沃尔什函数有如下乘法定理沃尔什函数有如下乘法定理 wal i twalj twal ij t( , )( , )(, )(4122)

2、 并且，该定理服从结合律并且，该定理服从结合律( , )( , )( , )( , ) ( , )( , ), , ,()wal i twal j twal k twal i twal j twal k ti j kp 0 1 221(4123) hKhKjgigpkjgpkigpktkRtkRtkRtjwaltiwal)()(10)(10)(10), 1(), 1(), 1(),(),( 证明：由定义式证明：由定义式 但是但是10100101011211),1(),1(),1(),1(1),1(),1(1),1(),1(0,1)( ,)(tkRtkRtkRtkRtkRtkRtkRtkRjgi

3、gkk),(), 1(), 1(),(),()()(10)()(10tjiwaltkRtkRtjwaltiwalkkkkjgigpkjgigpk因此因此以上便是乘法定理的证明。以上便是乘法定理的证明。（）（） 沃尔什函数有归一化正交性沃尔什函数有归一化正交性 i 1 i 0),(),(10jjdttjwaltiwal(4124) 证明：由乘法定理有证明：由乘法定理有dttlwaldttjiwaldttjwaltiwal),(),(),(),(101010其中 ijl 由于由于 ,3 ,2, 1 0),(1),0(1010idttiwaldttwal所以当所以当l=0，即，即i=j时，则时，则1

4、),(),(10dttjwaltiwal而当而当 ，即，即 时，则时，则l 0ij0),(),(10dttjwaltiwal正交性得证。正交性得证。 3.3.5 沃尔什变换沃尔什变换 离散沃尔什变换可由下二式表达离散沃尔什变换可由下二式表达10),()(1)(NttiwaltfNiW (4127)10( )( )( , )Nif tW iwal i t (4128)离散沃尔什变换解析式写成矩阵式可得到沃尔什变换离散沃尔什变换解析式写成矩阵式可得到沃尔什变换矩阵式矩阵式)1()1 ()0()()1()1 ()0(NWWWNwalNfff WWW NNwal NfffN( )( )()()( )(

5、 )()0111011 (4130)式中)(NWal代表N阶沃尔什矩阵。另外，沃尔什函数可写成如下形式1011)()1(),(phkkkpiittiwal式中 因此，可得到指数形式的沃尔什变换式二进)(01221tttttttkpp二进)(01221tttiiiikppNp 2 (4133)以上是离散沃尔什变换的三种定义，其中矩阵式最以上是离散沃尔什变换的三种定义，其中矩阵式最为简洁。为简洁。11101()0( )( )( 1)ppkkkkNittif tW i 10)(1011) 1()(1)(NtiitpkkkhptfNiW (4132) 3.3.6 离散沃尔什离散沃尔什-哈达玛变换哈达玛

6、变换 由沃尔什函数的定义可知，按哈达玛排列的沃尔什由沃尔什函数的定义可知，按哈达玛排列的沃尔什函数与按沃尔什排列的沃尔什函数相比较只是排列函数与按沃尔什排列的沃尔什函数相比较只是排列顺序不同，其本质并没有什么不同。但是哈达玛矩顺序不同，其本质并没有什么不同。但是哈达玛矩阵具有简单的递推关系，也就是高阶矩阵可用低阶阵具有简单的递推关系，也就是高阶矩阵可用低阶矩阵的直积得到，这就使得沃尔什一哈达玛变换有矩阵的直积得到，这就使得沃尔什一哈达玛变换有许多方便之处。因此，用得较多的是沃尔什哈达许多方便之处。因此，用得较多的是沃尔什哈达玛变换。玛变换。 离散沃尔什哈达玛变换的定义可直接由沃尔什离散沃尔什哈

7、达玛变换的定义可直接由沃尔什变换得到，只要用按哈达玛排列的沃尔什函数去变换得到，只要用按哈达玛排列的沃尔什函数去代替沃尔什排列的沃尔什函数，就可以得其矩阵代替沃尔什排列的沃尔什函数，就可以得其矩阵式如下：式如下： (4-134)WWWNNHNfffN( )( )()()( )( )()0111011 式中，式中， 是沃尔什哈是沃尔什哈达玛变换系数序列达玛变换系数序列，是时间序列，是时间序列， ，p为正整数。式为正整数。式(3134)的逆的逆变换式如下变换式如下 ) 1(),2(),1 (),0(NWwWW )1(),2().1 (),0(NffffNP 2fff NH NWWW N( )( )

8、()()( )( )()011011 (4135) 例：将时间序列1 ,1 ,0,0,1 ,1 ,0,0做沃尔什哈达玛变换及反变换。 0000021-021 11001100 1-1 1 1-1 1-1-11 1 1-1-1-1-1 11 1-1 1-1-1 1-11-1-1-1-1 1 1 11 1-1-1 1 1-1-11-1-1 1 1-1-1 11-1 1-1 1-1 1-11 1 1 1 1 1 1 181)7()6()5()4()3(2)1()0(WWWWWWWWffffffff( )( )( )( )( )( )( )01345671 11 11 1(2) 1 1 1 1 1 1

9、- 1-1 1-1 1-1 -1-1 1 1-1-11-1-1 1 1-1-1 11 1 1 1-1-1-1-11-1 1-1-1 1-1 11 1-1-1-1-1 1 11-1-1 1-1 1 1-1 -1 1 2020000000110011反变换为3.3.7 离散沃尔什变换的性质离散沃尔什变换的性质 离散沃尔什变换有许多性质。下面把主要性质列举离散沃尔什变换有许多性质。下面把主要性质列举于下。为叙述方便起见，用于下。为叙述方便起见，用 表示时间序表示时间序列，用列，用 表示变换系数序列，以表示变换系数序列，以 表示沃尔什变换对应关系。表示沃尔什变换对应关系。 )(tf)(nW)()(nW

10、tf (4136) ( )( )f tW n11( )( )ftW n22af tafta W na W n11221122 ( )( )( )( )（）线性（）线性 如果如果 则则 其中其中为常数为常数。21aa ，（）（） 模模2移位性质移位性质 将时间序列将时间序列 作作L位模位模2移位所得到的序列，移位所得到的序列，我们称为模我们称为模2移位序列。移位序列。)(tf模模2移位是这样实现的：移位是这样实现的： 设：设： 是周期长度为是周期长度为N的序列。作一个新的序列的序列。作一个新的序列)1(),2(),1 (),0()(Nfffftf () ( ), ( ), ( )()z mzzz

11、z Nl0121 (4137) 其中其中 。此时，称。此时，称 是序列是序列f(t)的位模的位模2移位序列。移位序列。)()(ltfmzlmz)( )7 (),6 (),5 (),4 (),3 (),2 (),1 (),0 ()(fffffffftf 例：例：8NL=2， 则：)27 (),26 (),25 (),24(),23 () 22(),21 (),20 ()(2ffffffffmz 由于B0) 1 (0=0) 1 (00) 0 0(20=(2)D B1) 1 (0=0) 1 (01) 0 0(21=(3)D B0) 0 (0=0) 1 (00) 1 0(22=(0)D B1) 0 (

12、0=0) 1 (01) 1 0(23=(1)D B0) 1 (1=0) 1 (00) 0 1 (24=(6)D B1) 1 (1=0) 1 (01) 0 1 (25=(7)D B0) 0 (1=0) 1 (00) 1 1 (26=(4)D B1) 0 (1=0) 1 (01) 1 1 (27=(5)D所以 ( ) ( ), ( ), ( ) ( ), ( ), ( ), ( ), ( )z mffffffff223016745同理 ( ) ( ), ( ), ( ) ( ), ( ), ( ), ( ), ( )z mffffffff332107654用矩阵表示为 zMf11式中 ( ) (

13、)ff tz m z11M10 1 0 0 0 0 0 01 0 0 0 0 0 0 00 0 0 1 0 0 0 00 0 1 0 0 0 0 00 0 0 0 0 1 0 00 0 0 0 1 0 0 00 0 0 0 0 0 0 10 0 0 0 0 0 1 0(4139) (4140) IMMT zMf22按照模按照模2和的性质，可知和的性质，可知这里这里I是么阵。是么阵。(4141) 0 0 1 0 0 0 0 00 0 0 1 0 0 0 01 0 0 0 0 0 0 00 1 0 0 0 0 0 00 0 0 0 0 0 1 00 0 0 0 0 0 0 10 0 0 0 1 0

14、 0 00 0 0 0 0 1 0 02M模模2移位性质是指下面的关系移位性质是指下面的关系：如果如果 ，并且，并且 是是 的模的模2移位序列，移位序列，)()(nWtf1)( tz)(tf则则 ( )( )z tW nz1ltz )(式中：式中： ， 是矩阵是矩阵 中中的第的第n行第行第l列的元素；列的元素； n=0，1，2，(N-1)；t=0，1，2，(N-1)； ，p是正整数。是正整数。)(),()(nWlnwalnWz),( lnWalPwal2Np 2此定理可证明如下：此定理可证明如下： 令 为 的元素， 是 的模2移位序列，则)(tzltz )(ltz)()( tf1010),()

15、(1 ),()(1)(NtNtztnwalltfNtnwaltzNnW令令 ，则有，则有 ，并且当，并且当 t 取值由取值由0到到N-1时，时，r 也取同样的值，只不过取值的顺序不也取同样的值，只不过取值的顺序不同而已。于是可写成如下形式：同而已。于是可写成如下形式：rtl trl )(),(),()(1),(),()(1),(),(),()(1 ),()(1)(10101010nWlnwaltnwaltfNlnwalrnwalrfNlnwallnwalrnwalrfNlrnwalrfNnWNtNrNrNrz所以，证明所以，证明 又因为又因为 ，这说，这说明明 与与l无关。也就是说，模无关。也

16、就是说，模2移位后的序列，作移位后的序列，作沃尔什变换后，所得到的第沃尔什变换后，所得到的第n个系数的平方个系数的平方 与与模模2移位的移位位数无关。移位的移位位数无关。 仍然等于仍然等于 。)()(nWtzzl2222z)n(W)n(W)l ,n(wal)n(W2)(nWz2)(nWz2)(nWz2)(nW 因此，模2移位定理（或称为并元移位定理）又可表达为输入序列 模2移位后的功率谱是不变的。)(tf 例如：设输入序列 ，对此序列作l=3的模2移位，得1 1, 0, 0, 1, 1, 0, , 0)(tf ( )() , , , , , , , z tf t3311 0 0 11 0 0作

17、沃尔什变换得0 0, 0, 0, ,0, -0, 21,21)(nW根据)(),()(nWlnwalnWz 可得0 0, 0, 0, ,0, 0, 21,21)()3,()(),()(nWnwalnWlnwalnWz41)21()0(22W41)21()3(22W41)21()0(22zW可见n相同时，功率也相同，也就是说功率列率谱是不变的。41)21()3(22zW（） 模2移位卷积定理（时间） 在讨论下面的定理之前，首先说明一下模2移位卷积与模2移位相关的概念。 令 和 是两个长度相同的周期性序列。用下式来定义两个序列的模2移位卷积：)(1tf)(2tf)()(1)()(1)(210121

18、0112ltflfNltflfNtCNlNl (4143) )()(1)(210112ltflfNtKNl(4144) )(12tC式中 为模2卷积， 为模2减运算符，它的运算结果与模2加一样。 模2移位相关的定义式如式(3144)所示其中其中 表示模表示模2移位相关，移位相关， 是是 的模的模2移移位序列。位序列。 )(12tK)(2ltf)(2tf由式由式(4143)和和(4144)可见，模可见，模2移位卷积和模移位卷积和模2移位移位相关具有相同的结果相关具有相同的结果,即：即： )()(1)()(21011212ltflfNtCtKNl如果用如果用W 代表作沃尔什变换，则：代表作沃尔什变

19、换，则：( )( )( )CtW nW n1212则(4145) ( )( )( )( )ftWnftWn1122 下面讨论模下面讨论模2移位卷积定理移位卷积定理 如果如果)()(),()()(),()()(),()()(),()()(),()()(nWnWlnwallfNnWlnwalnWlfNtnwalltfNlfNtnwalltflfNNtnwaltCNtCNlNlNlNtNtNlNt12101210211010211010211012121 1 11 11 1W所以证明( )( )( )CtW nW n1212 （）（） 模模2移位列率卷积定理移位列率卷积定理模2移位列率卷积由下式来表

20、示1121201120( )( )( )()( )()NrNrW nW nW rW r nW rW rn(4146) 依照模2时间卷积定理，模2移位列率卷积定理如下 如果 ( )( )( )( )ftWnftWn1122( )( )( )( )W nW nftft1212(4147) 仿照模2移位时间卷积定理的证明方法可得到证明。则（）（） 模模2移位自相关定理移位自相关定理 从模从模2 移位时间卷积移位时间卷积(相关相关)定理可以得到模定理可以得到模2移位移位自相关定理。只要把定理中的自相关定理。只要把定理中的 和和 换成换成 和和 便立即可以得到模便立即可以得到模2移位移位自相关定理。自相

21、关定理。)(2tf)(2nW)(1tf)(1nW() () KtWn1 112(4148) 其证明方法也与模其证明方法也与模2移位时间卷积定理的证明一样。移位时间卷积定理的证明一样。 从式从式(4148)可以建立一个重要概念：可以建立一个重要概念：模模2移位自相关序列的沃尔什变换等于序列的功移位自相关序列的沃尔什变换等于序列的功率谱。率谱。也就是说，模也就是说，模2移位下的自相关序列的移位下的自相关序列的沃尔什变换正好与序列的功率谱相符合。沃尔什变换正好与序列的功率谱相符合。与傅与傅里叶变换相比较，模里叶变换相比较，模2 2移位下的自相关与沃尔移位下的自相关与沃尔什谱的关系相当于线性移位下的自

22、相关序列的什谱的关系相当于线性移位下的自相关序列的离散傅里叶变换与其功率谱的关系。离散傅里叶变换与其功率谱的关系。 （）帕斯维尔定理（）帕斯维尔定理 如果 )()(11nWtf则 10211021)()(1NnNtnWtfN(4149) 证明：设 )()(11nWtK)()(2111nWtK因为因为是自相关函数，所以是自相关函数，所以)(11tK10211011),()(),()()(NnNntnwalnWtnwalnWtK 则 又由于又由于 )()(1)(111011ltflfNtKNl所以所以),()()()(110211110tnwalnWltflfNNlNl 如果令如果令t=0，则11

23、2011201NflWnlNnN( )()10211021)()(1NnNtnWtfN由于由于 l 仅是求和运算的变量，因此将仅是求和运算的变量，因此将 l 换成换成 t ，即可得：即可得：（）循环移位定理（）循环移位定理 把序列把序列 循环地向左移若干位，例如移循环地向左移若干位，例如移l位，位，l=1，2，N-1，这样得到的序列叫循环移位，这样得到的序列叫循环移位序列。如果用序列。如果用 来表示循环移位序列，来表示循环移位序列，)( tf)(ltz1 , 2, , 1 )1(),2(),1( ),()(Nllflflflftzl(4150) ( ) ( ), ( ), ( ),() ( )

24、, ( ), ( ), ( )z tffffffff334567012 ( ) ( ), ( ), ( ),() ( ), ( ), ( ), ( )z tffffffff556701234 ( ) ( ), ( ), ( ),() ( ), ( ), ( ), ( )f tffffffff01234567 例如：有一个N=8的序列， 当l=5，l=3 的循环移位序列分别为循环移位定理的内容如下：循环移位定理的内容如下： 如果如果 和它的循环移位序列和它的循环移位序列 的沃尔的沃尔什哈达玛变换分别是什哈达玛变换分别是 和和 ，则，则)(tf)(ltz)(nWf)(nWz1222122222)1

25、()1()()()0()0(rrrrnfnzfznWnWWW (4151) 式中 1,log, 2, 12 N,1,2,l Np p r 这个定理把序列的沃尔什哈达玛变换系数与循环这个定理把序列的沃尔什哈达玛变换系数与循环移位序列的沃尔什哈达玛变换系数联系了起来。即移位序列的沃尔什哈达玛变换系数联系了起来。即某些之某些之 和和 与与 之和是相等的。所以这之和是相等的。所以这个定理又称为个定理又称为沃尔什哈达玛变换的循环移位不变沃尔什哈达玛变换的循环移位不变性性。下面用一个例子来说明本定理的意义。下面用一个例子来说明本定理的意义。)(2nWf)(2nWz 例如设 ，经沃尔什哈达玛变换后的系数序列

26、为 1 1, 0, 0, 1, 1, 0, , 0)(tf( ) ,W nf12 0, -12, 0, 0, 0, 0, 0现将 做l=3的循环移位，则 )(tf ( ) ,z t31 0, 0, 1, 1, 0, 0, 1此序列经沃尔什哈达玛变换后的系数序列为0 0, 0, 0, ,0,21 0, ,21)(nWz0)1()(0)1()(2212222122)1()1(ffnzznWnWWnWrrrr当r=1时，则： WWzf220014( )( )(nWf)(nWz21) 0 (,21) 0 (fzWW 所以 0) 1 () 1 (22fzWW当2r时，则： 41)3()2()(41)3(

27、)2()(222122222122)1()1(fffnzzznWWnWWWnWrrrr所以 41)3()2()3()2(2222ffzzWWWW 0)7()6()5()4()(222122221fffnffWWWWnWrr0)7()6()5()4()(222122221zzznzzWWWWnWrr当r=3时，则：所以：所以： )7() 6() 5 () 4()7() 6() 5 () 4(22222222ffffzzzzWWWWWWWW显然，这些关系符合循环移位定理。显然，这些关系符合循环移位定理。 需要特别指出的是这个定理只适用于沃尔什哈达需要特别指出的是这个定理只适用于沃尔什哈达玛变换。此

28、定理的更加一般性的证明，请参阅有关书玛变换。此定理的更加一般性的证明，请参阅有关书籍。籍。4.3.8 快速沃尔什变换快速沃尔什变换 离散付里哀变换有快速算法。同样，离散沃尔什变离散付里哀变换有快速算法。同样，离散沃尔什变换也有快速算法。利用快速算法，完成一次变换只换也有快速算法。利用快速算法，完成一次变换只须须 次加减法，运算速度可大大提高。次加减法，运算速度可大大提高。当然快速算法只是一种运算方法，就变换本身来说当然快速算法只是一种运算方法，就变换本身来说快速变换与非快速变换是没有区别的。快速变换与非快速变换是没有区别的。NN2log 由于沃尔什哈达玛变换有清晰的分解过程，而由于沃尔什哈达玛

29、变换有清晰的分解过程，而且快速沃尔什变换可由沃尔什哈达玛变换修改且快速沃尔什变换可由沃尔什哈达玛变换修改得到，所以下面着重讨论沃尔什哈达玛快速变得到，所以下面着重讨论沃尔什哈达玛快速变换。换。式中式中 为正整数。为正整数。pNp,2W nNHf tH( ) ( )1 (4152) 由离散沃尔什哈达玛变换的定义可知由离散沃尔什哈达玛变换的定义可知 这里以这里以8阶沃尔什哈达玛变换为例，讨论其分解阶沃尔什哈达玛变换为例，讨论其分解过程及快速算法。由克罗内克积可知过程及快速算法。由克罗内克积可知： 210444422222222222244442222222244444444442480000000

30、000000000000000000000GGGIIIIII II IIIIH H H HIIIIHH HH HHHHIIIIHHHHHHHHH - - - - - - - 0 0 - (3153) 其中GHHHH02222000 0 0 0 0 0 0 0 0 0 (4154) 222222221- 0 0 0 0 0 0 - 00 IIIIIIIIG(4155) GIIII24444 -(4156) 其中 均为么阵42 , II由上面的分解有WnGGGf tH( )( )18012 (4157) 令 ftGftftGftftGft12211302( )( )( )( )( )( ) 则 W

31、 nf tH( )( )183 下面是具体计算 的公式及流程图。)( , )( , )(321tftftfftGft12( )( )fffffffffffff111111110123456710123( )( )( )( )( )( )( )( )( )( )( )( )( 0 0 0 1 0 0 00 1 0 0 0 1 0 00 0 1 0 0 0 1 00 0 0 1 0 0 0 11 0 0 0 -1 0 0 00 1 0 0 0 -1 0 00 0 1 0 0 0 -1 00 0 0 1 0 0 0 -145670415263704152637)( )( )( )( )( )( )(

32、 )( )( )( )( )( )( )( )( )( )( )( )( )fffffffffffffffffff(4158) )7()5()6()4()7()5()6()4()3() 1 ()2()0()3() 1 ()2()0()7()6()5()4()3()2() 1 ()0( 1-0 1 0 0 0 0 00 1-0 1 0 0 0 01 0 1 0 0 0 0 00 1 0 1 0 0 0 00 0 0 0 1-0 1 00 0 0 0 0 1-0 10 0 0 0 1 0 1 00 0 0 0 0 1 0 1)7()6()5()4()3()2() 1 ()0(11111111111

33、111111111111122222222ffffffffffffffffffffffffffffffffftGft211( )( )ftGft302( )( )fffffffffff33333333222012345671012( )( )( )( )( )( )( )( )( )( )( ) 1 0 0 0 0 0 01-1 0 0 0 0 0 00 0 1 1 0 0 0 00 0 1-1 0 0 0 00 0 0 0 1 1 0 00 0 0 0 1-1 0 00 0 0 0 0 0 1 10 0 0 0 0 0 1-1 fffffffffffffffffffff22222222222

34、2222222222345670101232345456767( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )( )(4160) 图图 414 快速沃尔什快速沃尔什-哈达玛变换信号流图（哈达玛变换信号流图（N=8）离散沃尔什哈达玛变换的特点：离散沃尔什哈达玛变换的特点： 1）、）、WHT只有加减运算，没有乘除运算，运只有加减运算，没有乘除运算，运算速度快；算速度快； 2）、）、H是对称矩阵，是对称矩阵，H=H ,所以，正所以，正反变换均用一样的公式，一样的运算程序，甚反变换均用一样的公式，一样的运算程序，甚至用一样的硬件，

35、给工程带来极大方便。至用一样的硬件，给工程带来极大方便。4.3.9 多维变换多维变换图像处理多用二维变换，二维变换的定义图像处理多用二维变换，二维变换的定义：H)v,u(WH)y .x(f H)y,x(fHN1)v,u(wH2H 二维沃尔什哈达玛变换可用一维沃尔什哈达玛变换来计二维沃尔什哈达玛变换可用一维沃尔什哈达玛变换来计算，其步骤如下：算，其步骤如下： (1)、以、以 ，对，对 中中 个列中的每一列做变换，个列中的每一列做变换，得到得到 ； (2)、以、以 对对 中中 行的每一行作变换，行的每一行作变换，即可得到二维变换系数即可得到二维变换系数 。根据这一步骤，便可以。根据这一步骤，便可以

36、利用一维快速沃尔什哈达玛变换来完成二维沃尔什哈利用一维快速沃尔什哈达玛变换来完成二维沃尔什哈达玛变换的计算。达玛变换的计算。xNN),(yxfyN),(yuWxyNN),(yuWxxN),(vuWxy另外一种计算方法是将二维沃尔什哈达玛变换当做另外一种计算方法是将二维沃尔什哈达玛变换当做一维来计算。这种方法是将数据矩阵的各列依次顺序一维来计算。这种方法是将数据矩阵的各列依次顺序排列，这样就形成由排列，这样就形成由 个元素的列矩阵。然后再个元素的列矩阵。然后再按照一维沃尔什哈达玛变换方法来计算。下面用实按照一维沃尔什哈达玛变换方法来计算。下面用实例说明一下两种计算方法。例说明一下两种计算方法。y

37、xNN 例：设数据矩阵如下例：设数据矩阵如下 (,)fxy 1 3 12 1 2 21),(yxf),(yxf求求 的二维沃尔什哈达玛变换的二维沃尔什哈达玛变换。 W uHx( , ) 01212121121232 11-1 -1首先对首先对 的每一列作变换：的每一列作变换： 第一列 Wux( ,) 11-1 2121321251第三列第三列 W ux( , ) 1121111220 11-1 第二列第二列 第四列第四列 Wux( ,) 11-1 -1312112123( , )W u yx123 2 5 3-1 0 1-1 所以所以 Wux( , )014121325318133 1 1 1

38、1-1 1-11 1-1-11-1-1 1 -3-1第一行第一行 对对 每一行作变换每一行作变换 ),(yuWx第二行第二行 W ux( , ) 11412118 1 1 11-1 1-11 1-1-11-1-1 1 -1 0 1-1-1 1-1-3最后得到二维变换系数矩阵最后得到二维变换系数矩阵 3- 1- 1 1-1- 3- 3 1381),(vuWxy以上是采用第一种算法得到的结果。以上是采用第一种算法得到的结果。 第二种算法如下：第二种算法如下： 将将 改写成列矩阵改写成列矩阵Y，即，即 ),(yxf21231121 Y对对 Y做一维变换做一维变换 WWWWWWWWYyyyyyyy(

39、)( )()( )()( )( )()01234567181 1 1 1 1 1 1 11 - 1 1 - 1 1 - 1 1 - 11 1 - 1 - 1 1 1 - 1 - 11 - 1 - 1 1 1 - 1 - 1 11 1 1 1 - 1 - 1 - 1 - 11 - 1 1 - 1 - 1 1 - 1 11 1 - 1 - 1 - 1 - 1 1 11 - 1 - 1 1 - 1 1 1 - 1 -1 3 1-3-1-1-3121132121813然后重排一下然后重排一下3- 1- 1 1-1- 3- 3 1381),(vuWxy 显然，与第一种算法得到的结果一致。显然，与第一种算

40、法得到的结果一致。 哈达玛变换哈达玛变换函数函数哈达玛变换矩阵函数哈达玛变换矩阵函数 H=hadamard(n)返回返回nxn矩阵矩阵hadamard变换矩阵，变换矩阵， H*H=n*I Hadmard变换的应用变换的应用 sig=imread(lena.bmp); sig=double(sig)/255; subplot(121);imshow(sig) m_sig,n_sig=size(sig); sizi=8; snum=64; T=hadamard(sizi); hdcoe=blkproc(sig,sizi sizi,P1*x*P2,T,T); coe=im2col(hdcoe,siz

41、i sizi,distinct); coe_temp=coe; Y,Ind=sort(coe); m,n=size(coe); snum=m-snum; for i=1:n coe_temp(Ind(1:snum),i)=0; end re_hdcoe=col2im(coe_temp,sizi sizi,m_sig,n_sig,distinct); re_sig=blkproc(re_hdcoe,sizi sizi,P1*x*P2,T,T); subplot(122);imshow(re_sig,)( , )f x y dxy在 上的投影；Radon变换变换Radon变换是计算图象在某一指定角度射线上投影变换是计算图象在某一指定角度射线上投影的变换的变换f(x,y)在某方向上的投影在该方向上的线积分：在某方向上的投影在该方向上的线积分：( , )f x y dyx在 上的投影；cossinsincosxxyy Radon变换变换f(x,y)在某方向上的投影在该方向上的线积分：在某方向上的