第13章离散图像处理.ppt
DFT仅仅是数字图像处理中的一种变换，其实仅仅是数字图像处理中的一种变换，其实还有很多种变换。还有很多种变换。定义定义x是是N1的向量，的向量，T是是NN的矩阵，则：的矩阵，则：或定义了向量或定义了向量x的一的一个线性变换。个线性变换。1,0,0,1Nii jjjyt xiN其中y=Tx核矩阵例：二维坐标系统中的一个向量旋转例：二维坐标系统中的一个向量旋转1122cossinsincosyxyxT是非奇异的，则原向量是非奇异的，则原向量。对上例来说，对上例来说，相当于该向量反向旋转。相当于该向量反向旋转。若若T是酉矩阵，则是酉矩阵，则对对T的每个元素取共轭复数的每个元素取共轭复数转置转置当当T的所有元素都是实数时，的所有元素都是实数时，TTt的第的第(i,j)元素是元素是T的第的第i行与行与Tt的第的第j列（也就是列（也就是T的的第第j行）的内积，行）的内积，i=j时为时为1，否则为，否则为0。因此，。因此，例：一维例：一维DFT就是酉变换就是酉变换1201iNjkNkiif eNFF=Wf酉阵酉阵2,1ijkNi kweN线性酉变换产生一个有线性酉变换产生一个有N个变换系数的向量个变换系数的向量y，每个每个变换系数都是输入向量变换系数都是输入向量x和变换矩阵和变换矩阵T的某一行的内积。的某一行的内积。反变换也类似。反变换也类似。正变换可看作是一个分解过程：将信号向量分解成正变换可看作是一个分解过程：将信号向量分解成它的各个基元分量，这些基元分量自然以基向量的形式它的各个基元分量，这些基元分量自然以基向量的形式表示，变换系数规定了在原信号中各分量所占的量。表示，变换系数规定了在原信号中各分量所占的量。反变换可看作是一个合成过程：通过将各分量相加反变换可看作是一个合成过程：通过将各分量相加来合成原始向量。来合成原始向量。上述过程的上述过程的：任一个向量都能唯一地分解：任一个向量都能唯一地分解为分别具有为分别具有“合适合适”幅度的一组基向量，然后通过将这幅度的一组基向量，然后通过将这些分量相加可以重构原向量。变换系数的个数与向量的些分量相加可以重构原向量。变换系数的个数与向量的元素个数是相同的。元素个数是相同的。变换后的向量是原始向量的一种表示，可由它完整变换后的向量是原始向量的一种表示，可由它完整地恢复出原始向量。因此它是原始向量的另一种形式。地恢复出原始向量。因此它是原始向量的另一种形式。将一个将一个NN的矩阵的矩阵F变换成另一个变换成另一个NN阵阵G。11,00(,),0,1NNm ni kikGFi k m ni k m nNT变换的变换的，是，是N2N2的块矩阵，每行的块矩阵，每行N块，共块，共N行，行，m,n用于寻块，用于寻块，i,k用于块内寻元素用于块内寻元素(,)(,)(,)rci k m nT i m T k nT若：若：11,00(,)(,)NNm ni kcrikGF T k nT i m 则：则：m=1m=2m=Nn=1n=2n=N例：二维例：二维DFT，是可分离的、对称的酉阵。是可分离的、对称的酉阵。正变换：正变换：GWFW，反变换：反变换：FW*tGW*t与与FT不同，许多变换在其核矩阵不同，许多变换在其核矩阵T中只有实元素，中只有实元素，而实数酉阵是正交的，因此，而实数酉阵是正交的，因此，FTtGTt。若若T是对称阵，正反变换相同，则：是对称阵，正反变换相同，则：GTFT，FTGT再进一步，如果两个分量相同，则变换是对称的：再进一步，如果两个分量相同，则变换是对称的：(,)(,)(,)i k m nT i m T k nT11,00(,)(,)NNm ni kikGT i mF T k n则：则：记为记为反变换：反变换：核矩阵的各行构成了核矩阵的各行构成了N维向量空间的一组基向维向量空间的一组基向量，这些行是正交的，即：量，这些行是正交的，即：TT*tI或：或：1*,0Nj ii kj kiT T其中其中 j,k是是Kronecker函数：函数：当当j=k时时 j,k=1，而而当当j k时时 j,k=0。任一组正交向量集都可用于一个线性变换，但任一组正交向量集都可用于一个线性变换，但通常通常。如。如FT用复用复指数作基函数。指数作基函数。二维反变换可以看作是通过将一组被适当地加二维反变换可以看作是通过将一组被适当地加权的基图像求和而重构原图像。变换矩阵权的基图像求和而重构原图像。变换矩阵G中的每中的每个元素就是其对应的基本图像在求和时所乘的倍个元素就是其对应的基本图像在求和时所乘的倍（系）数（即权值）。（系）数（即权值）。一幅基图像可通过对只含有一个非零元素（令一幅基图像可通过对只含有一个非零元素（令其值为其值为1）的系数矩阵进行反变换而产生，）的系数矩阵进行反变换而产生，N2个这个这样的矩阵产生样的矩阵产生N2幅基本图像。设其中一个系数矩阵幅基本图像。设其中一个系数矩阵为：为：其中其中i,j分别为行和列的下标，分别为行和列的下标，p,q是标明非零元素位置的整数。是标明非零元素位置的整数。,p qi p j qG反变换：反变换：11,00(,)(,)(,)(,)NNm nip k qikFT i mT k nT p m T q n这样，对于一个可分离的酉变换，每幅基本图这样，对于一个可分离的酉变换，每幅基本图像就是变换矩阵某两行的外积像就是变换矩阵某两行的外积。基图像可看作是分解原图像所得的单位集分量，基图像可看作是分解原图像所得的单位集分量，同时也是组成原图像的基本结构单元。同时也是组成原图像的基本结构单元。正变换通过确定系数来实现分解，反变换通过正变换通过确定系数来实现分解，反变换通过将基图像加权求和来实现重构。将基图像加权求和来实现重构。由于存在着无限多组基图像集，从而也就存在由于存在着无限多组基图像集，从而也就存在着无限多的变换。而某一组特定的基图像集仅对相着无限多的变换。而某一组特定的基图像集仅对相应的变换有重要的意义。应的变换有重要的意义。DFT的核矩阵：的核矩阵：0,00,11,01,1NNNNwwwwW2,1ikjNi kweN虚指数具有周期性，因此虚指数具有周期性，因此W是酉矩阵。是酉矩阵。一维一维DFT：，N1的的向量向量N1的的向量向量上图是当上图是当f是实向量时，谱向量是实向量时，谱向量F中各频率分量所中各频率分量所处的位置。零频和最高频率仅出现一次，其他分量以处的位置。零频和最高频率仅出现一次，其他分量以共轭复数的形式出现两次。如果共轭复数的形式出现两次。如果Ft被看作是一个行向量，被看作是一个行向量，则前面的则前面的N/2+1个元素是谱的右半边，后个元素是谱的右半边，后N/2-1个元素个元素在左半边。当在左半边。当fN是奈奎斯特折叠频率（采样频率的一半）是奈奎斯特折叠频率（采样频率的一半）时，对应于时，对应于F的第的第i个元素的频率是：个元素的频率是：20/22()/2 11NiNifiNNsNifNiNN 0 1 i N/2 N-1 N/2 N/2 N/2 N/2 N/2 fN fN fN fN fN0 i N/2 -(N-i)-如果如果f的后的后N/2个元素是第一个元素到第个元素是第一个元素到第N/2-1个元个元素的镜像，则素的镜像，则F为实。为了生成一个适于画出频谱的向为实。为了生成一个适于画出频谱的向量，可对量，可对F循环右移（或左移）循环右移（或左移）N/2个元素，这样零频个元素，这样零频率元素就会位于率元素就会位于N/2。而它两边的频率分别向两个方向而它两边的频率分别向两个方向递增。奈奎斯特频率元素仅在递增。奈奎斯特频率元素仅在f0出现。出现。也可以利用傅立叶变换的平移定理，即：也可以利用傅立叶变换的平移定理，即：020()()()()()(1)()ujxj xxNF uf xF uuef xef xf x 平移量平移量u0N/2时时上式意味着，在执行上式意味着，在执行DFT之前，改变之前，改变f(x)的奇号元的奇号元素的符号，可使谱移到适于绘图的位置。素的符号，可使谱移到适于绘图的位置。GWFWFW*tGW*t矩阵矩阵矩阵矩阵矩阵矩阵11223344F0,0对对4个像限重新排列，使显示更为方便。此时，零个像限重新排列，使显示更为方便。此时，零频落在矩阵中心，并沿径向增长。频落在矩阵中心，并沿径向增长。(,)(,)(/2,/2)(1)(,)x yF u vf x yF uNvNf x y 通过改变图像矩阵通过改变图像矩阵F中一半元素的符号可中一半元素的符号可以得到所要的平移。以得到所要的平移。1100(21)(21)(,)()()(,)coscos22NNcikimknG m na m b ng i kNN1100(21)(21)(,)()()(,)coscos22NNcmnimkng i ka m b n G m nNN12(0),(),1aa mmNNNGcCgC矩阵矩阵矩阵矩阵矩阵矩阵,(21)()cos2i mimCa mN11002(1)(1)(1)(1)(,)(,)sinsin111NNsikimknG m ng i kNNN11002(1)(1)(1)(1)(,)(,)sinsin111NNsmnimkng i kG m nNNN,2(1)(1)sin11i kikTNN核矩阵的元素：核矩阵的元素：11,0012cas()NNm ni kikGgimknNN基函数：基函数：11,0012cas()NNi km nmngGimknNN形式相同形式相同cascossin2cos(/4)核矩阵：核矩阵：,1cas 2i kikTNN哈特利变换是相应傅立叶变换的实部减去虚部。哈特利变换是相应傅立叶变换的实部减去虚部。傅立叶变换是哈特利变换的偶部减去傅立叶变换是哈特利变换的偶部减去j乘以奇部。乘以奇部。()()()()()()()()eog xf xh xG vF v H vFv H v对称、可分离的酉变换，元素均为对称、可分离的酉变换，元素均为 1，且，且N=2n，（，（n整数）整数）211111122H/2/2/2/211NNNNNNNHHHHH8111111111111111111111111111111111111111112 2111111111111111111111111H07341625相应于矩阵相应于矩阵行的符号变行的符号变化次数。每化次数。每一行的这个一行的这个数都不同。数都不同。称为列率称为列率（sequency）81 11 111111 11 111111 11111111 11111111111 111112 2111 1111111 11111111 111111H01234567按列率重新按列率重新安排各行，安排各行，得到得到，也称作沃尔也称作沃尔什（什（Walsh）核矩阵核矩阵哈达玛变换的基函数就是沃尔什函数，因此，哈哈达玛变换的基函数就是沃尔什函数，因此，哈达玛变换也叫作沃尔什变换。达玛变换也叫作沃尔什变换。有的书中仅将有有的书中仅将有序哈达玛变换称为沃尔什变换。序哈达玛变换称为沃尔什变换。斜变换的酉矩阵从斜变换的酉矩阵从22的哈尔或哈达玛阵开始的哈尔或哈达玛阵开始2111112S/2/210100000010010120000NNNNNNNNNNNababSSSbabaIIII222222314141NNNNabNNN=8时的斜变换基函数时的斜变换基函数傅立叶变换的基函数间仅是傅立叶变换的基函数间仅是不同。而哈尔函数不同。而哈尔函数在在和和上都是不同的。（上都是不同的。（）基函数索引：基函数索引：和和的的。令整数令整数0 k N-1由其它两个整数由其它两个整数p和和q唯一决定，即：唯一决定，即：21pkqk和和p、q互为函数。对任意互为函数。对任意k0，2p是使是使2p k的的2的最大幂，而的最大幂，而q1是余数。是余数。定义哈尔函数：定义哈尔函数：01()h xN22112222112()2220ppppkppqqxqqh xxNelsek p q0 0 01 0 12 1 13 1 24 2 15 2
