专利名称：：基于四元数分解和相关函数判断彩色图像相关性的方法
技术领域：
：本发明涉及一种判断彩色图像相关性的方法，主要利用四元数分解和相关函数的计算，属于图像处理领域。
背景技术：
：四元数，也称Hamilton数，是1843年由英国数学家哈密顿(W.R.Hamilton)首先发现的，建立四元数理论最初的目的是为研究空间矢量找到类似解决平面问题中使用的复数方法。设g-a+^'+^+汲，a，6,c,dei是实数，/，/,fceC代表虚数单位，若有/2=_/2=A2=_1，=A:,y.A:=/力'=_/称g为Hamilton四元数，a为四元数g的实部，W+qZ+汲为g的虚部。特别地当c"-0时，g就是复数，当"c-"0时，g就是实数，故四元数是实数和复数的扩充。从Hamilton四元数的定义上可以看出，如果改变/JA之间的运算定义规则，就可以出现其它不同于Hamilton四元数的定义。如定义/2=/=42=-1，!y=》=jt,=^=汰=-_/，贝提有人称之为"超复数"的四元数。在相当长的一段时间里，四元数没有为人们所重视，更没有得到实际的应用。近20年来它逐渐成为国内外相关领域研究的热点，并且已经在刚体动力学、惯性导航、航天器姿态控制、机器人技术、计算机图形学、信号处理、图像处理等领域得到广泛应用。由于^/,乘法交换律对Hamilton四元数不再适用，直接计算一个四元数乘法需要计算16个实数乘法运算。针对四元数特有的性质，研究定义在四元数域上的计算方法及其快速算法具有重要的理论和实际意义。在彩色图像处理中，可以用纯四元数形式表示彩色图像，纯四元数g-W+q/'+汲中的6，c^分别表示彩色图像的红(Red)、绿(Green)、蓝(Blue)分量，这样可以在多维空间上把彩色图像的色彩分量作为矢量整体进行处理，可以充分考虑它们之间的色彩关联。近年来，四元数广泛应用于彩色图像处理，如彩色图像边缘检测、彩色图像滤波、彩色图像配准等。-在图像配准、图像压縮和图像融合等图像处理领域，图像的相关性有着重要的应用，如可以通过寻找相关性的最大值进行图像配准，可以通过相关性衡量两帧图像间的冗余度，进而为压縮和融合提供依据。如何快速准确地度量图像的相关性是图像处理领域一个重要课题。过去，传统的相关性技术只能用于灰度图像。当处理彩色图像时，必须将彩色图像按其颜色色彩空间(对应R、G、B)的每一维分别看成是一个灰度图像，并分别进行处理，这显然忽略了彩色图像各色彩之间的内在关系。1999年Sangwine等人提出了基于四元数的彩色图像相关性度量，在多维空间上进行处理，比传统互相关性更能体现图像的色彩关联，它表达了两帧图像色彩之间的映射、旋转等，同时可以减少光照变化等的影响。对于大小为MxiV的两幅彩色图像"m，n)和g(m,w)，其相关函数可表示成M-lW-l_ct(附,")=SZi《"，-附，"w)(1)其中《(w，w)=r(w，w)/+s(m，w)y'+"m，w)A,g(/w,=w)/+c(;w，w)/+d(附,《)A:。这里(附，")为时域坐标，表示图像所处的位置，6,W分别表示彩色图像的红、绿、蓝分量，^^表示g(m，w)的共轭复数。对于式(1)这样的值域为复数域的函数，可以利用巻积定理即时域上的巻积等价于频域上相乘的性质计算其相关函数，但由于四元数的乘法不满足交换律，使得无法用类似手段计算四元数的相关函数。2000年ELL和Sangwine提出把四元数Fourier变换(傅立叶变换)分解成平行和垂直的分量来计算两个四元数的相关函数。即对于两个纯四元数《(m，m)=《)/+s(/w,m)_/+"/w，禾口g(/w,w)=6(/w，w)f+c(m,w)_/+，《(/w，w)为单位四元数(艮卩—+^+^=1)，可)l每g分解为g二g〃+g丄，其中.2(2)而且满足<formula>formulaseeoriginaldocumentpage7</formula>公式(3)可以理解为四元数的"乘法交换律"，分解公式(2)需要21次实数乘法。对于上述给定的两幅彩色图像"/n，w)和g(m，w),ELL等人给出了相关函数cr(附，w)的计算公式为<formula>formulaseeoriginaldocumentpage7</formula>或cr<formula>formulaseeoriginaldocumentpage7</formula>其中有关变量定义如下对于二维实值函数变量x(附,w)，假如其大小为Mx7V，其右Fourier变换、右Fourier逆变换、左Fourier变换、左Fourier逆变换分别定义如下<formula>formulaseeoriginaldocumentpage7</formula>这里(附,")表示时域坐标；(",v)表示频域坐标；//为复数虚部单位，这里定义为特定单位四元数向量，>"=1±^。'对于四元数而言，由于分解后的水平和垂直分量可以满足一定意义下的"乘法交换律"，使得利用巻积定理来计算相关函数的思想得以推广到四元数空间，采用Sangwine提出的方法分解四元数的水平和垂直分量需要21次实数乘法，利用Sangwine等提出的基于平行垂直分解的方法计算两幅彩色图像的相关函数，最好的结果计算量为層(6.51og2層+105)，计算速度低，仍然有待进一步提高。
发明内容本发明的目的是解决采用ELL提出的方法分解四元数的水平和垂直分量，以及基于平行垂直分解的方法计算两幅彩色图像的相关函数计算速度低的问题，提供了一种基于四元数分解和相关函数判断彩色图像相关性的方法。本发明方法具体包括以下步骤步骤一、大小均为MxiV像素的两幅彩色图像的纯四元数表达式分别为二夂(w,")/+a(/w,w)y'+《(w,)女禾卩")=62(/w，")/+c2(附,w)/+<i2(7W，，其中，W附，")、c，(w，")4(附，")ei是实数，分别代表第一幅彩色图像p(w，")的红、绿、蓝分量；W附，")A(m,"X(附，")ei是实数，分别代表第二幅彩色图像《(/n,")的红、绿、蓝分量；C是虚数单位，步骤二、将第二幅彩色图像沿//=^1方向分解为g(/w，w)=《〃(w,w)+《±(/，w)，其中，第二幅彩色图像《(w，w)的水平分量^(m，")////，第二幅彩色图像《(m，")的垂直分量&(m，")丄//，;/为复数虚部单位，这里定义为特定单位四元数向量，第二幅彩色图像《(m,")的水平分量^(m，w)和垂直分量^(m，")按下述公式计《丄(附，w)，其中，=[62O，")+c2(w，")+d2(w，")]/J，步骤三、按下述公式获取两幅彩色图像的相关函数cr(m,n):ct(w，")=—厕F[F/("，vOF:"(",v)+(w,v)F^(w，v)]其中，(w,")表示时域坐标，(M,v)表示频域坐标，《("，v)为第一幅彩色图像函数的右傅立叶变换，F/("，力为第一幅彩色图像函数的右傅立叶逆变换，i^("，v)为第二幅彩色图像《(w，")的水平分量《,(m,")函数的右傅立叶逆变换，F^O/，v)为第二幅彩色图像《(w,")的垂直分量&(/n,")函数的右傅立叶逆变换，所述相关函数CT(m，W)的值表示两幅彩色图像的相关性。.在实际应用时，将本发明的方法获得的相关函数cr(;n,")的值与设定的阈值进行对比，当所述相关函数cr(附，")的值大于设定的阈值时，说明两幅彩色图像的相关性较强；否则，说明两幅彩色图像的相关性较弱。本发明的优点1、利用本发明方法平行垂直分解四元数的速度明显优于ELL等提出的方法，ELL理论中每个像素的分解需要3x3+3x4-21次实数乘法运算，而本发明方法仅需要3+1=4次实数乘法；2、利用本发明方法判断两幅彩色图像的相关性，计算相关函数的速度明显优于ELL理论中相关函数的计算速度，ELL等提出的理论最好的结果计算量为MV(6.51og2層+105)，本发明方法的计算量为層(5.51og2層+43)。图1是本发明方法的流程图。具体实施例方式具体实施方式一下面结合图1说明本实施方式，本实施方式方法包括以下步骤步骤一、大小均为MxiV像素的两幅彩色图像的纯四元数表达式分别为JE(7W，W)=Z\(W，《)/+q(W，/2)7.+《(W，W)A禾口《(W，W)=办2(W，+C2(777，W)/+"2(W，W)允,其中，W附,")、q(m,")4(m,")ei是实数，分别代表第一幅舉色图像p(m，")的红、绿、蓝分量；62(附，")《2(/，")、《(附，")€及是实数，分别代表第二幅彩色图像《(w,")的红、绿、蓝分量；eC是虚数单位，步骤二、将第二幅彩色图像《(m，")沿//=^4^方向分解为g(w，w)-《〃(w，w)+《丄(/w，w)，其中，第二幅彩色图像《(w，")的水平分量^(/n，")〃//，第二幅彩色图像《(m，")的垂直分量^(m，")丄^，//为复数虚部单位，这里定义为特定单位四元数向量，第二幅彩色图像《(w，")的水平分量^(/n,")和垂直分量^—,")按下述公式计算j一,("1《丄(附,")=《-"叫//其中，cr叫=[62(附,m)+c2(附,w)+<(附,")]/V^(8)步骤三、按下述公式获取两幅彩色图像的相关函数cr(w,"):cK附,")=—[F/(",v)F^(w，v)+(m，v);丄"(w，v)](9)其中，(w，")表示时域坐标，(m，v)表示频域坐标，"(",v)为第一幅彩色图像函数;^,")的右傅立叶变换，F/("，v)为第一幅彩色图像函数K附，")的右傅立叶逆变换，'F^("，v)为第二幅彩色图像《(m，")的水平分量^(m,")函数的右傅立叶逆变换，.《(M，v)为第二幅彩色图像《(w,")的垂直分量&(m,")函数的右傅立叶逆变换，所述相关函数cr(附，")的值表示两幅彩色图像的相关性。下面分两部分介绍本发明的方法一、本发明提出的对第二幅彩色图像g(m，")沿;/方向分解的公式(7)与ELL等人提出的公式(2)是等价的，具体推导过程如下ELL提出的公式(2)为'根据z2=/=a2一i,y=、#=a/=，可得10《g《=+s2+f)6-2r(r6++,+j'((一+s2+f2)c-2s(r6+sc+,+他2++一)c-2/(r6+sc+,由于《为单位四元数，有2"2=1记~=r6+sc+^/则利用公式(3)的乘法交换律，将&《g-2J代入公式(2)，可得g〃="J(10)公式(10)表明，纯四元数g沿单位四元数《分解的快速算法与ELL提出的公式(2)是等价的，用公式(10)分解的速度要明显优于ELL提出的四元数分解方法的速度。在本发明中用四元数《(附^)=62(附，一+(:2(/，")_/+《(附，")*代替ELL理论中的gOw,")^^,")"^^,")^^^)*,二者都是纯四元数，涉及各参量代表的意义相同；用更为特殊的特定单位四元数向量//=(/+_/+^)/^代替￡1^理论中的单位四元数"m，")-Kw,")Z+s(m，"X/'+"/w，")it，很显然，根据推导出的公式(10)可得公式(7)成立，艮P:上述关于本发明对彩色图像的分解方法与ELL提出的分解方法等价的论述，以及本发明提出的两幅彩色图像的相关函数与ELL提出的两幅彩色图像的相关函数等价的论述，说明本发明方法的结论与ELL提出的理论是等价的，但本发明方法的计算量比现有方法要少很多，进而提高了运算速度，节省了计算所用的资源、提高了工作效率。ELL理论中分解水平垂直分量所涉及公式为式(2),每个像素的分解需要3x3+3x4=21次实数乘法运算。利用本发明提出的分解公式(7)，应用于彩色图像时，得到同样的分解结果，仅需要3+1=4次实数乘法；应用于实部不为零的四元数分解是也只需要3+3=6次实数乘法。从结果上分析，本发明分解的速度明显优于ELL提出的理论。作为上述公式(10)的扩展，在上面分解(10)和(9)中，如果g的实部不为零，公式(10)的分解依然有效，只是可将g的实部保留在水平分量中，相应的(3)式变为一,,g/"W=&,g丄f(11)其中5表示g的共轭复数，即1=-汲。由于彩色图像处理中，g为纯四元数，因此下面只讨论g为纯四元数的情况。二、基于时域的四元数相关函数cr(w，n)计算方法针对ELL等人提出的公式(4)或(5)，本发明提出了用于彩色图像处理的基于时域分解的四元数相关函数计算公式，'使相关函数的计算更加方便。对于上述给定的两幅大小为MxiV的彩色图像P(/n，")和《(m，")，其相关函数计算公式(9)为—,")=-VM尸[F/(M,《("，v)+("，《(",v)]公式(9)与ELL提出的公式(4)是等价的，利用第一部分所述四元数分解方法，并考虑p(m,")、g(m,")为纯四元数的情况，推导过程力口下大小为MxW的第二幅彩色图像)沿//方向分解成g(m，打)=g〃(am,m)+g±(m，m)其中g〃(m，w)〃；>(w，w)，g丄(m，w)丄；(m，w)由公式(6)的定义可知:M—1W-1附-0tt=0(12)其中//为复数虚部单位，这里定义为特定单位四元数向量，//=/+y+A:将公式(1)中g换成《，《换成p，并带入公式(12)得:<formula>formulaseeoriginaldocumentpage13</formula>-S2>""SIX","),Mw(18)jc=0"0m=0根据公式(6)，将公式(18)整理为.;c=0^=0x=0少-O=-V^尸/(w,v)F,0,v)_VMF/(m,v)&:"(m,v)=-(",v)i^("，v)+(w,(M,v))于是获得本发明方法所述的公式(9):cr(/n，")=—v]i^F-k(",v)i^("，v)+尸/(",y)/^("，v)]。第二幅彩色图像《(附，")的水平分量^(w，")函数的右傅立叶逆变换F广(",v)采用实数快速傅立叶变换方法获得，由于^(m，")////，A=，，可以把/z看成常规复数虚部单位，按常规实数右Fourier逆变换处理，其计算量为顯0.51og2(磨)+l)。本发明计算两幅彩色图像相关函数cr(附，")的公式(9)时涉及第一幅彩色图像函数P(附，")的右傅立叶变换F/(",力、第一幅彩色图像函数;;(m，")的右傅立叶逆变换F/(",v)、第二幅彩色图像《(m，")的水平分量^(m，")函数的右傅立叶逆变换i^(",力和第二幅彩色图像g(m,/i)的垂直分量l(rn,")函数的右傅立叶逆变换F^(",v)的计算量分析如下1、第一幅彩色图像函数的右傅立叶变换F/("，v)和第一幅彩色图像函数p(m，")的右傅立叶逆变换F/(",v)的计算按如下方法进行/(w,w)=61(m,")/+c1(/w,")/+《(m，")A;，实数&,(w,")，c,(m，")，《(m，")的Fourier变换分别记为尸m(w，v)，0，v)，i^(w，v)，即FM(M，V)=FXm,"))=~^ZJ，")e17"(19)1m-lw-l一"2ff，+!)i^(w，v卜FK(c,(m,"))-"^tJ^(m,")eM"(20)VMATw=o"=o>/MVm=o"=o则(w,v)=F61(w，v>'+FclO,v)_/+Fdl(w，v)A(22)("，v)=i^("，V>.+Fcl(",v)_/+Fdl(",v》(23)在利用式(22)和式(23)计算相关函数时，可利用上述性质减少计算量，计算量为層(l.5log2(磨)+3)。式(22)和式(23)的推导过程如下."("，vh"^(《S6(附，")eM,m=0/=0m=0fi=0=FA(m,v)/+Fc(w，v)+Fd(m,v)A,力""、—_f_〃"VVW附M、y、A/w'、",v卜"^('工》(附，")em=0w=0附=0=《+Fc(m,v》.+(m,v)A:2、第二幅彩色图像《(m，")的水平分量、〃(m，")函数的右傅立叶逆变换i^(",v)采用实数快速傅立叶变换方法获得，由于^(m，")////，A=，，可以把/z看成常规复数虚部单位，按常规实数右Fourier逆变换处理，其计算量为層(0.51og2(層)+l)。3、第二幅彩色图像《(m，")的垂直分量&(m，")函数的右傅立叶逆变换获取的方法为第二幅彩色图像《(附,")的垂直分量&(附，")函数的纯四元数表达式为<formula>formulaseeoriginaldocumentpage16</formula>e及是实数，/,,A:eC是虚数单位，第二幅彩色图像《(m,n)的垂直分量i(w，n)函数的右傅立叶逆变换/^(",v)按如下公式获得'<formula>formulaseeoriginaldocumentpage16</formula>其中，F/表示实数"(m,w)的右傅立叶逆变换；F/表示实数风m,w)的右傅立叶逆变换；F〃表示实数/(w，")的右傅立叶逆变换；w"/0代表取其复数的实部；^MO代表取其复数的虚部。，计算量为MV(1.51og2(層)+3)。综上，两幅彩图像相关函数cr(附，")中涉及的,/("，v)、i^(",v)vP/("，v)和K/(",v)]按如上所述的方法进行计算，则cr(附，")的计算最后涉及的F，……]按如下方法计算设<formula>formulaseeoriginaldocumentpage16</formula>则<formula>formulaseeoriginaldocumentpage16</formula>两幅彩图像相关函数CT(加，/0按本发明方法进行计算，最终的计算量为層(5.510&層+43)，而利用ELL等提出的基于平行垂直分解的方法，两幅彩图像相关函数cr(附，")最好的结果计算量为廢(6.51og2MV+105)，江淑红提出的不基于平行垂直分解的方法(参见江淑红，郝明非，张建秋，胡波.快速超复数傅氏变换和超复数互相关的新算法及应用。电子学报，2008:36(1),100-105.)，计算量为MV(71o^MV+36)。从结果分析，本发明方法明显优于它们。表1、表2是本发明的方法与上述两种方法计算彩色图像相关函数计算量比较。表1计算彩色图像相关函数的计算量(所需实数乘法次数)8x8320x240512x5121200x1920本发明的方去48641.0157x1073.7224x1073.6690x101ELL等的方法92161.6165x1075.8160x1075.5845x10:江淑红等的方法49921.1489x1074.2467x1074.2382x10'表2本发明的方法与已有<table>tableseeoriginaldocumentpage17</column></row><table>权利要求1、基于四元数分解和相关函数判断彩色图像相关性的方法，其特征在于，实现该方法包括以下步骤步骤一、大小均为M×N像素的两幅彩色图像的纯四元数表达式分别为p(m，n)＝b1(m，n)i+c1(m，n)j+d1(m，n)k和q(m，n)＝b2(m，n)i+c2(m，n)j+d2(m，n)k，其中，b1(m，n)、c1(m，n)、d1(m，n)∈R是实数，分别代表第一幅彩色图像p(m，n)的红、绿、蓝分量；b2(m，n)、c2(m，n)、d2(m，n)∈R是实数，分别代表第二幅彩色图像q(m，n)的红、绿、蓝分量；i，j，k∈C是虚数单位，步骤二、将第二幅彩色图像q(m，n)沿μ方向分解为q(m，n)＝q//(m，n)+q⊥(m，n)，其中，第二幅彩色图像q(m，n)的水平分量q//(m，n)//μ，第二幅彩色图像q(m，n)的垂直分量q⊥(m，n)⊥μ，μ为复数虚部单位<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>μ</mi><mo>=</mo><mfrac><mrow><mi>i</mi><mo>+</mo><mi>j</mi><mo>+</mo><mi>k</mi></mrow><msqrt><mn>3</mn></msqrt></mfrac><mo>,</mo></mrow>]]></math>id="icf0001"file="A2009100723190002C1.tif"wi="23"he="10"top="137"left="73"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>第二幅彩色图像q(m，n)的水平分量q//(m，n)和垂直分量q⊥(m，n)按下述公式计算<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>q</mi><mrow><mo>/</mo><mo>/</mo></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>α</mi><mi>qμ</mi></msub><mi>μ</mi></mtd></mtr><mtr><mtd><msub><mi>q</mi><mo>⊥</mo></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>q</mi><mo>-</mo><msub><mi>α</mi><mi>qμ</mi></msub><mi>μ</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math></maths>其中，<mathsid="math0003"num="0003"><math><![CDATA[<mrow><msub><mi>α</mi><mi>qμ</mi></msub><mo>=</mo><mo>[</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>/</mo><msqrt><mn>3</mn></msqrt><mo>,</mo></mrow>]]></math>id="icf0003"file="A2009100723190002C3.tif"wi="73"he="6"top="191"left="44"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>步骤三、按下述公式获取两幅彩色图像的相关函数cr(m，n)<mathsid="math0004"num="0004"><math><![CDATA[<mrow><mi>cr</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msqrt><mi>MN</mi></msqrt><msup><mi>F</mi><mrow><mo>-</mo><mi>R</mi></mrow></msup><mo>[</mo><msubsup><mi>F</mi><mi>p</mi><mi>R</mi></msubsup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msubsup><mi>F</mi><msub><mi>q</mi><mrow><mo>/</mo><mo>/</mo></mrow></msub><mrow><mo>-</mo><mi>R</mi></mrow></msubsup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>F</mi><mi>p</mi><mrow><mo>-</mo><mi>R</mi></mrow></msubsup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msubsup><mi>F</mi><msub><mi>q</mi><mo>⊥</mo></msub><mrow><mo>-</mo><mi>R</mi></mrow></msubsup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>]</mo></mrow>]]></math></maths>其中，(m，n)表示时域坐标，(u，v)表示频域坐标，FpR(u，v)为第一幅彩色图像函数p(m，n)的右傅立叶变换，Fp-R(u，v)为第一幅彩色图像函数p(m，n)的右傅立叶逆变换，id="icf0005"file="A2009100723190002C5.tif"wi="7"he="6"top="246"left="29"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>(u，v)为第二幅彩色图像q(m，n)的水平分量q//(m，n)函数的右傅立叶逆变换，id="icf0006"file="A2009100723190003C1.tif"wi="15"he="5"top="29"left="31"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>为第二幅彩色图像q(m，n)的垂直分量q⊥(m，n)函数的右傅立叶逆变换，所述相关函数cr(m，n)的值表示两幅彩色图像的相关性。2、根据权利要求1所述的基于四元数分解和相关函数判断彩色图像相关性的方法，其特征在于，第一幅彩色图像函数P(m，")的右傅立叶变换g("，v)按如下公式获得-F/("，v)-尸w(M,v》+i^("，v)7'+F力(w,v)A:其中<formula>formulaseeoriginaldocumentpage3</formula>3、根据权利要求1所述的基于四元数分解和相关函数判断彩色图像相关性的方法，其特征在于，第一幅彩色图像函数^(m,")的右傅立叶逆变换Fp-、,v)按如下公式获得i^(",v)=FM(M,v>'+0，v)/+Frfl("，v)A:其中，《,(M,v)与Fw(",v)共轭，Fd(",v)与Fd("，v)共轭，i^(M，v)与F^"，v)共轭。4、根据权利要求1所述的基于四元数分解和相关函数判断彩色图像相关性的方法，其特征在于，第二幅彩色图像《(/n,w)的水平分量^(m，w)函数的右傅立叶逆变换F^(",v)采用实数快速傅立叶变换方法获得。5、根据权利要求1所述的基于四元数分解和相关函数判断彩色图像相关性的方法，其特征在于，获得第二幅彩色图像"加，w)的垂直分量^(加，")函数的右傅立叶逆变换尸广(",力的方法为第二幅彩色图像《(m，")的垂直分量^(w，w).函数的纯四元数表达式为《丄(附,w)=or(/,w)/+"(w，w)_/+;K附，or(附,")，风w,")，y(w,")e及是实数，！》UeC是虚数单位，第二幅彩色图像《(附，")的垂直分量^(附^)函数的右傅立叶逆变换々/(",v)按如下公式获得F;丄及("，v)=—/(F,)+—g(F/》+(F/)+—g(尸/))其中，F^表示实数"(m,")的右傅立叶逆变换；F/表示实数/(w，")的右傅立叶逆变换；F"表示实数W附，n)的右傅立叶逆变换；"/0代表取其复数的实部；//^g()代表取其复数的虚部。全文摘要基于四元数分解和相关函数判断彩色图像相关性的方法，属于图像处理领域，本发明的目的是解决采用ELL提出的方法分解四元数的水平和垂直分量，以及基于平行垂直分解的方法计算两幅彩色图像的相关函数计算速度低的问题。本发明的方法包括步骤一、大小均为M×N像素的两幅彩色图像的纯四元数表达式分别为p(m，n)＝b<sub>1</sub>(m，n)i+c<sub>1</sub>(m，n)j+d<sub>1</sub>(m，n)k和q(m，n)＝b<sub>2</sub>(m，n)i+c<sub>2</sub>(m，n)j+d<sub>2</sub>(m，n)k；步骤二、将第二幅彩色图像q(m，n)沿μ＝(i+j+k)/√3方向分解为q(m，n)＝q<sub>∥</sub>(m，n)+q<sub>⊥</sub>(m，n)；步骤三、按公式cr(m，n)＝-√MNF<sup>-R</sup>[F<sup>R</sup><sub>p</sub>(u，v)F<sup>-R</sup><sub>q∥</sub>(u，v)+F<sup>-R</sup><sub>p</sub>(u，v)F<sup>-R</sup><sub>q⊥</sub>(u，v)]获取两幅彩色图像的相关函数cr(m，n)，所述相关函数cr(m，n)的值表示两幅彩色图像的相关性。本发明方法用于判断两幅彩色图像的相关性强弱。文档编号G06F17/14GK101582167SQ200910072319公开日2009年11月18日申请日期2009年6月18日优先权日2009年6月18日发明者朱春辉,毅沈,晶金申请人:哈尔滨工业大学