第五章 基本图像变换ppt课件

1、长安大学地测学院长安大学地测学院1图像处理图像处理彭晓东图像工程第五章第五章 基本图像变换基本图像变换长安大学地测学院长安大学地测学院2第五章 基本图像变换v本章实质:频率域变换v重中之重:傅里叶变换长安大学地测学院长安大学地测学院3第5章 图象变换基础 为了有效和快速地对图象进行处理，常常为了有效和快速地对图象进行处理，常常 需要将原定义在图象空间的图象以某种形式转换需要将原定义在图象空间的图象以某种形式转换到另外一些空间，并利用在这些空间的特有性质到另外一些空间，并利用在这些空间的特有性质方便地进行一定的加工，最后再转换回图象空间方便地进行一定的加工，最后再转换回图象空间以得到所需的效果。

2、这些转换方法就是本章要着以得到所需的效果。这些转换方法就是本章要着重介绍和讨论的图象变换技术重介绍和讨论的图象变换技术 变换是双向的，或者说需要双向的变换。变换是双向的，或者说需要双向的变换。在图象处理中，一般将从图象空间向其他空间的在图象处理中，一般将从图象空间向其他空间的变换称为正变换，而将从其他空间向图象空间的变换称为正变换，而将从其他空间向图象空间的变换称为反变换或逆变换变换称为反变换或逆变换 长安大学地测学院长安大学地测学院4第5章 图象变换基础长安大学地测学院长安大学地测学院55.1 可分离和正交图象变换v分离变换目的:简化计算,长安大学地测学院长安大学地测学院65.1 可分离和正

3、交图象变换1 , , 1 , 0),()()(10NuuxhxfuTNx正向变换核1 , , 1 , 0),()()(10NxuxkuTxfNu反向变换核长安大学地测学院长安大学地测学院75.1 可分离和正交图象变换 1010),(),( ),(NxNyvuyxhyxfvuT 1010),(),( ),(NuNvvuyxkvuTyxf反向变换核正向变换核变换核与原始函数及变换后函数无关长安大学地测学院长安大学地测学院8),(),(),(21vyhuxhvuyxh),(),(),(11vyhuxhvuyxh102),(),(),(NyvyhyxfvxT101),(),(),(NxuxhvxTvu

4、T5.1 可分离和正交图象变换长安大学地测学院长安大学地测学院9AFAT BAFABBTB 图象矩阵对称变换矩阵反变换矩阵变换结果BTBF BAFABF 1 AB1 AB5.1 可分离和正交图象变换反变换反变换完全恢复不完全恢复长安大学地测学院长安大学地测学院10BTBF 1 AB5.1 可分离和正交图象变换*1AAT1AA长安大学地测学院长安大学地测学院115.2傅里叶变换5.2.1 2-D傅里叶变换 5.2.2 傅里叶变换定理 5.2.3 快速傅里叶变换长安大学地测学院长安大学地测学院125.2傅里叶变换周期为 的周期函数用一系列三角函数 的和来表示为:2)sin(0tnAn10)sin(

5、)(nnntnAAtf这种展开称为谐波分析。其中， 为直流分量， 为一次谐波又做基波），而 依次称为二次谐波，三次谐波等等。0A)sin(11tA)2sin(22tA)3sin(33tA傅里叶级数三角级数）傅里叶级数三角级数）长安大学地测学院长安大学地测学院135.2 傅里叶变换若函数若函数 ( )f x以以2l为周期的光滑或分段光滑函数，即为为周期的光滑或分段光滑函数，即为(2 )( )f xlf x( )f x将将展开为傅里叶级数展开为傅里叶级数01( )(cossin)2kkkak xk xf xabll1( )cos()lnln xaf xdxll1( )sin()lnln xbf x

6、dxll0,1,2n 1,2,3n 长安大学地测学院长安大学地测学院145.2傅里叶变换欧拉公式欧拉公式:cossinixexix长安大学地测学院长安大学地测学院155.2 傅里叶变换01( )(cossin)kkkk xk xf xaablliiii1cos()21sin()2ixxxxxeexeei( )k xlkkf xC eii*11( ) d( )d22k xk xllllkllCf x exf x exll+复数形式的傅里叶级数长安大学地测学院长安大学地测学院165.2 傅里叶变换傅里叶展开的复数形式上式的物理意义为上式的物理意义为: 一个周期为一个周期为2l 的函数的函数 ( )

7、f x可以分解可以分解为频率为为频率为nl，复振幅为，复振幅为 nc的复简谐波的叠加的复简谐波的叠加 nl称为谱点，称为谱点，所有谱点的集合称为谱对于周期函数所有谱点的集合称为谱对于周期函数 ( )f x而言，谱是离散的而言，谱是离散的i( )k xlkkf xC e1( )2n xillklCf x edxl长安大学地测学院长安大学地测学院175.2 傅里叶变换傅里叶变换对(三种):ii1111( )( )d , ( )( )d22xxFf x exf xFeii221( )( )d , ( )( )d2xxFf x exf xFei2i233( )( )d , ( )( )dxxFf t

8、exf xFe(1)(1)(2)(2)(3)(3)12311( )( )()222FFF三者之间的关系:长安大学地测学院长安大学地测学院185.2傅里叶变换傅里叶谱频谱，幅度）傅里叶相位角傅里叶功率谱22 1/2|(*)| FRIarctanIR222(*) |(*)|PFRI长安大学地测学院长安大学地测学院195.2 傅里叶变换v傅立正反变换对是怎么确立的？假设非周期函数假设非周期函数 ( )f x是一个周期函数是一个周期函数 ( )g x的周期的周期 2l 时的极限情况。由此，时的极限情况。由此， ( )g x的傅里叶级数展开式的傅里叶级数展开式01( )(cossin)kkkk xk x

9、g xaabll在在 l 时的极限形式就是所要寻找的非周期函数时的极限形式就是所要寻找的非周期函数 ( )f x的傅里叶展开的傅里叶展开长安大学地测学院长安大学地测学院205.2傅里叶变换 设不连续的参量设不连续的参量 1 (0,1,2,), kkkkkkll则傅立叶级数写为：则傅立叶级数写为： 01( )(cossin)kkkkkg xaaxbx傅里叶系数为傅里叶系数为1( )cosd1( )sind lkklklkklaf xxxlbf xxxl2 (0)1 (0)kkk长安大学地测学院长安大学地测学院215.2傅里叶变换 对于系数对于系数 0a，假设，假设 lim( )dlllf xx有

10、限，那么有限，那么 01limlim( )d02llllaf xxl对于余弦部分：对于余弦部分：当当 ,0kll ，不连续参变量，不连续参变量 k变为变为 连续参量，以符号连续参量，以符号 替代对替代对 k的求和变为对连续参量的求和变为对连续参量 的积分。的积分。长安大学地测学院长安大学地测学院225.2傅里叶变换11001(cos)( )cosdcos1( )cosd cosd1( )cosd cosdlkklkkk xk xaf xxxlllk xf xxxkllf xxxx同理可得正弦部分同理可得正弦部分01( )sind sindf xxxx长安大学地测学院长安大学地测学院235.2傅

11、里叶变换若令若令1( )( )cosd1( )( )sindAf xxxBf xxx00( )( )cosd( )sindf xAxBx则傅里叶变换可表示为：长安大学地测学院长安大学地测学院245.2 傅里叶变换iiii1cos(), sin()22xxxxixeexee 由欧拉公式得：00( )( )cosd( )sindf xAxBxii0011( ) ( )i ( )d ( )i ( )d22xxf xABeABe0ii011( ) ( )i ( )d (|)i (|)d22xxf xABeABei( )( )dxf xFe反变换反变换长安大学地测学院长安大学地测学院255.2.1 2-

12、D傅里叶变换傅里叶变换 1-D正变换对1个连续函数 f (x) 等间隔采样20( )f x13101525450575685790 x10/j2exp)(1)()(NxNuxxfNuFxfF长安大学地测学院长安大学地测学院265.2.1 2-D傅里叶变换傅里叶变换 1-D反变换 变换表达频谱幅度）相位角101/j2exp)()()(NuNuxuFxfuFF)( j exp )( )(j)()(uuFuIuRuF21 22)()( )( uIuRuF)()(arctan)( uRuIu 长安大学地测学院长安大学地测学院275.2.1 2-D傅里叶变换傅里叶变换 二维傅里叶变换频谱幅度）相位角功率

13、谱 1010/ )(j2exp),( 1),(NuNvNvyuxvuFNyxf21 22 ),(),( ),( vuIvuRvuF),(),(arctan),(vuRvuIvu 1010/ )(j2exp),( 1),(NxNyNvyuxyxfNvuF),(),( ),( ),(222 vuIvuRvuFvuP长安大学地测学院长安大学地测学院285.2.1 2-D傅里叶变换傅里叶变换傅里叶频谱图，就是图像梯度的分布图。（实际上图像上某一点与邻域点差异的强弱，即梯度的大小，也即该点的频率的大小 长安大学地测学院长安大学地测学院295.2.1 2-D傅里叶变换傅里叶变换v 频率域v 由傅立叶变换和

14、频率变量( u, v)定义的空间v 基本性质v （1变化最慢的频率成分(u=0, v=0)对应一幅图像的平均灰度v （2低频原点附近对应图像灰度变化慢的像素v （3高频远离原点对应图像灰度变化快的像素长安大学地测学院长安大学地测学院305.2.1 2-D傅里叶变换傅里叶变换长安大学地测学院长安大学地测学院315.2.2 傅里叶变换定理傅里叶变换定理 0、分离性质 1次2-D 2次1-D O(N 4)减为O(N 2) (0,0)XV(0,0)XY(0,0)UV列变换乘以行变换f (x, y)F(x, v)F(u, v)N(N 1)(N 1)-(N 1)-(N 1)-(N 1)-(N 1)-10/

15、xj2exp),(1),(NxNuvxFNvuF/exp),(=),(=101NyNvyj2yxfNNvxF长安大学地测学院长安大学地测学院321、平移定理 ),(/ )(j2exp),(dvcuFNdycxyxf5.2.2 傅里叶变换定理傅里叶变换定理 )( ),(/ )()(j2exp),( 1/ )(j2exp/ )(j2exp),( 1/ )(j2exp),(10101010dvcuFNydvxcuyxfNNvyuxNdycxyxfNNdycxyxfNxNyNxNy F),(),(vuFyxf长安大学地测学院长安大学地测学院335.2.2 傅里叶变换定理傅里叶变换定理2、旋转定理、旋转

16、定理借助极坐标：cossinxryr00( ,)( ,)f rF 说明：对 旋转 对应于将其傅里叶变换 也旋转。反之亦然。( , )f x y00( , )F x y长安大学地测学院长安大学地测学院345.2.2 傅里叶变换定理傅里叶变换定理v旋转定理示例旋转定理示例例5.2.3长安大学地测学院长安大学地测学院355.2.2 傅里叶变换定理傅里叶变换定理3、尺度定理(相似定理)( , )( , )af x yaF u1(,)(,)|uf ax byFaba bu原函数在幅度方面变化导致对其傅里叶变换在幅度方面产生对应的尺度变化。u原函数在空间尺度方面的缩放导致对其傅立叶变换在频域方面的相反缩放

17、。长安大学地测学院长安大学地测学院365.2.2 傅里叶变换定理傅里叶变换定理v尺度定理例5.2.4长安大学地测学院长安大学地测学院374、剪切定理（水平方向纯剪切 （垂直方向纯剪切 5.2.2 傅里叶变换定理傅里叶变换定理 byxxyy ),(),(buvuFybyxfxx ydxy),(),(vduuFydxxf图5.2.5 (a)与(b)长安大学地测学院长安大学地测学院385、组合剪切定理平移旋转尺度 水平剪切 垂直剪切 5.2.2 傅里叶变换定理傅里叶变换定理 bdvbubddvuFbdydxbyxf1,111),(xx11db101b101d图5.2.5长安大学地测学院长安大学地测学

18、院396、仿射定理u = (eu dv)/D和v = ( bu + av)/D 5.2.2 傅里叶变换定理傅里叶变换定理 bdaeedbavuebdavu长安大学地测学院长安大学地测学院407、卷积定理 2-D ( , )( , )( , ) ( , )f x yg x yF u v G u v( , ) ( , )( , )( , )f x y g x yF u vG u v5.2.2 傅里叶变换定理傅里叶变换定理 ( )( )( ) ( )f xg xF u G u( ) ( )( )( )f x g xF uG u-( , )( , ) ( , ) (,)d df x yg x yf p

19、 q g xp yqp q u两函数在空间的卷积与其傅立叶变换在频率域的乘积构成一对变换u两函数在空间的乘积与其傅里叶变换在频率域的卷积构成一对变换长安大学地测学院长安大学地测学院418、相关定理互相关：f (x) 与g(x)不是同一个函数 自相关：f (x) = g(x)2-D5.2.2 傅里叶变换定理傅里叶变换定理 *( )( )( ) ()df xg xfz g xzz),(),(*),(),(vuGvuFyxgyxf),(),(),(),(*vuGvuFyxgyxf qpqypxgqpfyxgyxfdd),(),(),(),(长安大学地测学院长安大学地测学院425.2.3 快速傅里叶变

20、换快速傅里叶变换直接进行一个N N的2-D傅里叶变换需要N4次复数乘法运算和N2(N2 1) 次复数加法运算1-D：复数乘法和加法的次数都正比于N2 快速傅里叶变换FFT）： 将复数乘法和加法的次数减少为正比于N log2N 逐次加倍法：复数乘法次数由N2减少为(N log2 N)/2 复数加法次数由N2减少为N log2 N 1 , , 1 , 0/j2exp)(1)()(10NuNuxxfNuFxfNxF长安大学地测学院长安大学地测学院435.2.3 快速傅里叶变换快速傅里叶变换 考察一维有限长序列x(n)(0=n=N-1)的傅立叶变换： 1, 1 ,0e1210NmnxNmXNmnjNn

21、其中：10*110*010212)()1()()0()()(e,eNnnNnnNnmnNjNjWnxXWnxXWnxmXWW即有：则令或记为或记为Wn，Wn -1长安大学地测学院长安大学地测学院445.2.3 快速傅里叶变换快速傅里叶变换) 1(.) 1 () 0 (.) 1(.) 1 () 0 () 1() 1(1) 1(0) 1() 1(11101) 1(01000NxxxWWWWWWWWWNXXXNNNNNN计算复杂性：一个频率分量需计算复杂性：一个频率分量需N N次乘法，次乘法，N-1N-1次加法次加法 整个变换需整个变换需N2N2次乘法，次乘法，N(N-1)N(N-1)次加法次加法矩

22、阵表示：矩阵表示：10)()(NnmnWnxmX长安大学地测学院长安大学地测学院455.2.3 快速傅里叶变换快速傅里叶变换为为周周期期是是以以，知知由由 2NWeWnmNj 结论：系数多数相同，且具有对称性结论：系数多数相同，且具有对称性mnNmnmnNNNjNWWWWeW 222221 9630642032100000WWWWWWWWWWWWWWWW例：例：N=40213192604;WWWWWWWWWW ；对对称称性性：；周周期期性性： 1010000010100000WWWWWWWWWWWWWWWW最小无重最小无重复运算复运算0W2W3W1WN=4长安大学地测学院长安大学地测学院465

23、.2.3 快速傅里叶变换快速傅里叶变换12/0212/011021)()()()()12, 1 ,0()12()()2()(:)(2NnmnNNnmnNNnmnNWnxWnxWnxmXNnnxnxnxnxnxN划分为为整数），现将序列（设)()()12()2()12()2(2112/02/12/02/12/0)12(12/0)2(mXWmXWWnxWnxWnxWnxmNmNNnmnNNnmnNNnnmNNnnmNmnNmNmnNjmNjmmnNjnmNmnNmnNjmnNjnmNWWeeeWWeeW2/2/22)2(2)12(2/2/222)2(N点的点的DFT转化转化为两个求为两个求N/2点

24、的点的DFT 1965 1965年年，库利，库利- -图图基提出把基提出把原始的原始的N N点点序列依次序列依次分解成一分解成一系列短序系列短序列，减少列，减少乘法运算乘法运算长安大学地测学院长安大学地测学院47)3,2,1 ,0()()4();()4(4)()()7()7()7()6()6()6()5()5()5()4()4()4()3()3()3()2()2()2()1()1()1()0()0()0(8848822112127812681258124812381228121812081mWWWmXmXmXmXDFTmXmXXWXXXWXXXWXXXWXXXWXXXWXXXWXXXWXXNm

25、mm的对称性有：的是周期为和又例子：)()()(21mXWmXmXmN5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院48)3()3()7()2()2()6() 1 () 1 ()5()0()0()4()3()3()3()2()2()2() 1 () 1 () 1 ()0()0()0(23812281218120812381228121812081XWXXXWXXXWXXXWXXXWXXXWXXXWXXXWXX则则有有： X1(0) X(0) 08W 08W X2(0) X(4) 蝶蝶形形运运算算单单元元 5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学

26、地测学院49X2(0)X2(1)X2(2)X2(3)08W18W28W38W-1-1-1-11111X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)N/2点点DFTN/2点点DFTX1(0)X1(1)X1(3)X1(2) 3 () 3 () 7 () 2 () 2 () 6 () 1 () 1 () 5 () 0 () 0 () 4 () 3 () 3 () 3 () 2 () 2 () 2 () 1 () 1 () 1 () 0 () 0 () 0 (2381228121812081238122812181208

27、1XWXXXWXXXWXXXWXXXWXXXWXXXWXXXWXX5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院50DFTMXmXmXmXmXWmXmXmXWmXmXNnnxnxnxnxNnnxnxnxnxDFTmHmGmNmN点的点的都是都是、而而，同样奇偶分组，同样奇偶分组点的点的都是都是和和又又2)()()()()()()()()()() 14/,.,1 , 0() 12()()2()() 14/,.,1 , 0() 12()()2()(4)()(654362524231262514135.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院51

28、 X5(0) X2 (0) X5 (1) X2 (1) X6 (0) X2 (2) X6 (1) X2 (3) X3(0) X1 (0) X3 (1) X1 (1) X4 (0) X1 (2) X4 (1) X1 (3) ) 1 () 1 () 3 () 0() 0() 2() 1 () 1 () 1 () 0() 0() 0() 1 () 1 () 3 () 0() 0() 2() 1 () 1 () 1 () 0() 0() 0(2)()()()(62852608526285260852428314083142831408316543XWXXXWXXXWXXXWXXXWXXXWXXXWXX

29、XWXXDFTMXmXmXmX点点的的都都是是、而而5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院52)6()2()1 ()6()2()0()4()0()1 ()4()0()0(284284083083xWxXxWxXxWxXxWxX)7()3()1 ()7()3()0()5()1 ()1 ()5()1 ()0(386386185185xWxXxWxXxWxXxWxX完整的蝶形图如下页：完整的蝶形图如下页：5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院53)7()3()5() 1 ()6()2()4()0(xxxxxxxx ) 1 ()0()

30、 1 ()0() 1 ()0() 1 ()0(66554433XXXXXXXX )3()2() 1 ()0()3()2() 1 ()0(22221111XXXXXXXX )7()6()5()4()3()2() 1 ()0(XXXXXXXX 5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院54时间复杂性：Nlog2N注意：输出X(m)是以m从小到大排列，而输入序列x(n)不是以n从小到大排列。采用码位倒序法排序：将自然顺序数转换成二进制数，然后首尾位倒序，再转换成十进制数，那么十进制数就是输入序列中的位置。如：N=8X(3)的位置是：011110(6)5.2.3 快速傅里

31、叶变换快速傅里叶变换长安大学地测学院长安大学地测学院55二维FFT 101021010210/210/2e ),(*1),(*e ),(1),(e )(*1)(*1e )()(NuNvNvyuxjNuNvNvyuxjNuNuxjNuNuxjvuFNyxfvuFNyxfuFNxfNuFxf 逆逆FFTFFT112 ()001122001( , )( , )e1( , )e( , )eux vyNNjNxyuxvyNNjjNNxyF u vf x yNF u vf x yNFFT根据可分离性有：二维变换用两次一维实现5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院56v傅立

32、叶变换注意的问题v 两个缺点：v (1)要进行复数运算，计算比较费时v 实用中还采用如沃尔什(Walsh)变换等。v (2)很多图像的高频项衰减的很快，在频域不清楚。v 解决方法：vuFvuD,1lg,5.2.3 快速傅里叶变换快速傅里叶变换长安大学地测学院长安大学地测学院575.3沃尔什/哈达玛变换5.3.1 沃尔什变换 5.3.2 哈达玛变换 5.3.3 关于两种变换的讨论沃尔什和哈达码变换都是可分离和正交变换 长安大学地测学院长安大学地测学院585.3.1 沃尔什变换沃尔什变换10)( )(1) 1( 1),(niubxbiniNuxh长安大学地测学院长安大学地测学院595.3.1 沃尔

33、什变换沃尔什变换N=4N=4和和8 8时的沃尔什变换核时的沃尔什变换核长安大学地测学院长安大学地测学院605.3.1 沃尔什变换沃尔什变换1010)( )(1) 1( )(1)(NxniubxbinixfNuW10)( )(1) 1(),(niubxbiniuxk1010)( )(1) 1( )()(NuniubxbiniuWxf长安大学地测学院长安大学地测学院6110)( )()( )(11) 1(1),(nivbybubxbiniiniNvuyxh 101010)( )()( )(11) 1( ),( 1),(NxNynivbybubxbiniiniyxfNvuW 101010)( )()

34、( )(11) 1( ),( 1),(NuNvnivbybubxbiniinivuWNyxf5.3.1 沃尔什变换沃尔什变换10)( )()( )(11) 1(1),(nivbybubxbiniiniNvuyxk长安大学地测学院长安大学地测学院6210)( )(10)( )(111111) 1(1) 1(1),(),(),(),(),(nivbybniubxbiniiniNNvykuxkvyhuxhvuyxh5.3.1 沃尔什变换沃尔什变换111( )( )( )( )01( , , , )( 1)iniininb x buby bvik x y u vN 正变换核反变换核长安大学地测学院长安

35、大学地测学院635.3.1 沃尔什变换沃尔什变换v沃尔什变换中,正向变换核与反向变换核只依赖于x,y,u,v，而与f(x,y)或W(u,v)的值无关。这些核可看一组基本函数，一旦图像尺寸确定，这些函数也就完全确定。N=4时的二维沃尔什变换图5.3.1长安大学地测学院长安大学地测学院645.3.1 沃尔什变换沃尔什变换长安大学地测学院长安大学地测学院655.3.1 沃尔什变换沃尔什变换v离散傅里叶变换是浮点数的运算，因此计算量会比较大，而且浮点数运算产生的误差会比较大；v沃尔什变换矩阵的系数是1或是1，只需要使用加法即可实现，计算复杂度小；v离散傅立叶转换相当于把信号拆解成在不同频率的正弦函数与

36、余弦函数的分量，而使用沃尔什转换相当于把信号拆解成在许多不同震荡频率的方波上，因而，除非所要分析的信号拥有类似方波组合的特性，否则沃尔什转换作频谱分析的效果会比使用离散傅立叶转换分析的效果要差，这是降低运算复杂度所要付出的代价。长安大学地测学院长安大学地测学院665.3.1 沃尔什变换沃尔什变换沃尔什变换具有某种能量集中。而且沃尔什变换具有某种能量集中。而且原始数据中数字越是均匀分布，经变换后原始数据中数字越是均匀分布，经变换后的数据越集中于矩阵的边角上。因此沃尔的数据越集中于矩阵的边角上。因此沃尔什变换可以压缩图像信息。且变换比傅立什变换可以压缩图像信息。且变换比傅立叶变换快。叶变换快。长安

37、大学地测学院长安大学地测学院6710)( )( ) 1(1),(niiiubxbNuxh5.3.2 哈达玛变换哈达玛变换10)()( 10) 1)(1)(NxubxbniiixfNuHN=8时，哈达玛变换核表5.3.2长安大学地测学院长安大学地测学院685.3.2 哈达玛变换哈达玛变换10)( )( ) 1(1),(niiiubxbNuxh10)( )( ) 1(),(niiiubxbuxk10)( )( 10) 1)()(NuubxbniiiuHxf长安大学地测学院长安大学地测学院695.3.2 哈达玛变换哈达玛变换10 )( )()( )( ) 1(1),(niiiiivbybubxbNv

38、uyxh10 )( )()( )( ) 1(1),(niiiiivbybubxbNvuyxk 1010 )( )()( )( 10) 1)(,( 1),(NxNyvbybubxbniiiiiyxfNvuH 1010 )( )()( )( 10) 1)(,( 1),(NuNvvbybubxbniiiiivuHNyxf长安大学地测学院长安大学地测学院705.3.2 哈达玛变换哈达玛变换v二维哈达玛变换正变换和反变换都可分成两个步骤计算，每个步骤用一个1维变换实现。长安大学地测学院长安大学地测学院715.3.3 关于两种变换的讨论关于两种变换的讨论v两种变换核里的数值都是1和-1，但是在行列的秩序上

39、两者有所不同。（下页说明）v在绝大多数图像变换应用中，常混合使用沃尔什变换和哈达玛变换，所以被统称为“沃尔什哈达玛变换，通常用来指两者中的任意一个。长安大学地测学院长安大学地测学院7210)( )( ) 1(1),(niiiupxbNuxg5.3.3 关于两种变换的讨论关于两种变换的讨论 )()()()()()()()(01121110ububupububupubupnnnn对比表5.3.2与表5.3.3长安大学地测学院长安大学地测学院73ux012345670+1+2+3+4+5+6+7+5.3.3 关于两种变换的讨论关于两种变换的讨论长安大学地测学院长安大学地测学院7411112HNNNNNHHHHH2NNHA1AFAT 11111111111111114H5.3.3 关于两种变换的讨论关于两种变换的讨论同样适合于哈达玛反变换长安