李正周，侯倩，戴真，付紅霞，葛豐增，金鋼

（1.重慶大學(xué)通信工程學(xué)院，重慶400044；2.中國空氣動力研究與發(fā)展中心，四川綿陽621000）

基于空時稀疏表示的紅外小目標(biāo)檢測算法

李正周1，侯倩1，戴真1，付紅霞1，葛豐增1，金鋼2

（1.重慶大學(xué)通信工程學(xué)院，重慶400044；2.中國空氣動力研究與發(fā)展中心，四川綿陽621000）

提出了一種基于過完備空時字典及其稀疏表示的紅外小弱目標(biāo)運動檢測算法。采用K奇異值分解算法學(xué)習(xí)連續(xù)多幀圖像的運動信息和形態(tài)特征，構(gòu)建自適應(yīng)形態(tài)過完備空時字典；利用高斯運動模型檢驗自適應(yīng)形態(tài)過完備空時字典，將其劃分為能分別描述目標(biāo)與背景的目標(biāo)過完備空時字典和背景過完備空時字典；將連續(xù)多幀圖像分別在目標(biāo)過完備空時字典和背景過完備空時字典上稀疏分解，利用幾個最大稀疏系數(shù)及其空時原子重構(gòu)信號，增強二者殘差來檢測小目標(biāo)信號。實驗結(jié)果表明，該過完備空時字典不僅能同時描述目標(biāo)與背景的運動信息和形態(tài)特征，極大地提高信號表示的稀疏程度，而且能有效增強目標(biāo)與背景的特征差異，提高小運動目標(biāo)的探測能力。

信息處理技術(shù)；小弱目標(biāo)檢測；空時超完備字典；目標(biāo)空時字典；背景空時字典；信號稀疏重構(gòu)

0 引言

當(dāng)探測器與目標(biāo)之間的距離較遠時，目標(biāo)在成像上表現(xiàn)為只占若干個像素的小目標(biāo)，并且易于淹沒在各種背景雜波和強噪聲中，這給目標(biāo)檢測與跟蹤帶來了很大的難度［1-2］。小弱運動目標(biāo)檢測方法可分為基于圖像濾波和基于模式識別（內(nèi)容學(xué)習(xí)）的兩類檢測算法［3-4］?；趫D像濾波的檢測算法首先白化圖像信號［5］，如Top-Hat［6］、Two-dimensional Least Mean Square（TDLMS）［7］和小波變換［8-9］等，然后采取諸如恒虛警率方法判定圖像中是否存在目標(biāo)。跟蹤前檢測（TBD）算法對于淹沒在強噪聲和強雜波中的弱小運動目檢測往往性能不佳［10-11］。TBD算法避免對每幀圖像進行門限判決，而在時域中沿目標(biāo)信號的疑似軌跡上進行包括目標(biāo)信號幅度和運動狀態(tài)在內(nèi)的窮盡假設(shè)，然后以此假設(shè)條件尋找小弱目標(biāo)，代表方法主要有Hough變換［12］、最大似然估計［13］、動態(tài)規(guī)劃算法［14］、假設(shè)檢驗［15］以及三維匹配濾波器［16-17］等，共同之處在于TBD算法沿目標(biāo)軌跡積累和提高目標(biāo)能量和信噪比，將目標(biāo)檢測轉(zhuǎn)換為目標(biāo)軌跡檢測。然而，TBD算法搜尋的疑似軌跡數(shù)量眾多，勢必造成計算量過大和計算復(fù)雜度過高，應(yīng)用范圍受到限制。

基于學(xué)習(xí)的檢測算法則是將目標(biāo)檢測問題轉(zhuǎn)化為模式識別中的二類分類問題，它對目標(biāo)模型和背景模型進行訓(xùn)練，并根據(jù)規(guī)則判定圖像是否含有目標(biāo)，如主成分分析（PCA）法［17］。近年來，基于超完備稀疏表示理論的紅外小目標(biāo)檢測方法發(fā)展成為一種新的學(xué)習(xí)算法，它通過學(xué)習(xí)與構(gòu)建超完備字典，將信號的能量只集中于較少且能表示信號的主要特征和內(nèi)在結(jié)構(gòu)的原子上。目前，具有代表性的超完備字典有高斯字典、Gabor字典、Gabor多成分字典等。高斯字典［18］適合于高斯分布的小弱目標(biāo)，然而小弱目標(biāo)形態(tài)動態(tài)變化，高斯字典難以適應(yīng)非高斯分布等非結(jié)構(gòu)形態(tài)?；谛螒B(tài)成分分析的自適應(yīng)信號稀疏表示法根據(jù)圖像信號構(gòu)造反應(yīng)目標(biāo)和噪聲的自適應(yīng)超完備字典，增強了目標(biāo)和背景的稀疏表示系數(shù)差異，但自適應(yīng)超完備字典其表示原子同時包含了背景和目標(biāo)信息，稀疏表示系數(shù)較小，難以分辨出目標(biāo)與背景；同時，構(gòu)造自適應(yīng)超完備字典采用單幀靜態(tài)圖像訓(xùn)練樣本，其原子只能表示目標(biāo)形態(tài)特征，難以反映出目標(biāo)運動信息和軌跡信息，因此基于此類空域超完備字典的檢測算法探測能力有限。研究能同時體現(xiàn)目標(biāo)在時域和空域稀疏性的超完備字典，或者將表示目標(biāo)運動特征和形態(tài)特征的超完備字典與跟蹤前檢測策略相結(jié)合，有望進一步提高紅外小弱目標(biāo)檢測能力。

1 紅外小弱目標(biāo)圖像的稀疏模型

紅外小弱目標(biāo)圖像f由目標(biāo)、背景和噪聲組成，即

式中:ft、fb和n分別表示目標(biāo)信號、背景信號和噪聲信號。

根據(jù)信號表現(xiàn)出來的特征差異，可將圖像的每個像素點看作目標(biāo)信號或者背景信號:

TBD算法通過探測目標(biāo)軌跡來檢測出目標(biāo)，即處理連續(xù)M幀圖像，即

式中:Ft（x，y，k）和Fb（x，y，k）分別表示連續(xù)的目標(biāo)信號和背景信號。

信號稀疏表示模型假設(shè)每個信號都能由相同類型的原子及其稀疏系數(shù)進行重構(gòu)。對于紅外圖像背景信號Fb（x，y，k），它可以通過背景空時原子線性表示［19］為

式中:Db表示背景過完備空時字典，；表示背景空時原子，每個空時原子都包含了M個空域背景原子；α表示背景信號在背景過完備空時字典的稀疏表示系數(shù)。

類似的，對于紅外圖像目標(biāo)信號Ft（x，y，k），它可以通過目標(biāo)空時原子線性表示為

式中:Dt表示目標(biāo)過完備空時字典，；表示目標(biāo)空時原子，每個空時原子都包含了M個空域目標(biāo)原子；β表示目標(biāo)信號在目標(biāo)過完備空時字典的稀疏表示系數(shù)。

紅外圖像序列F（x，y，k）可以通過聯(lián)合目標(biāo)過完備空時字典和背景過完備空時字典的所有原子線性稀疏表示為

式中:D=［DbDt］表示包含Dt和Db的過完備空時字典；γ=［α′Tβ′T］T為該字典的稀疏表示系數(shù)，是一個（Nt+Nb）M維的向量。如果F是目標(biāo)信號，則它不能由背景過完備空時字典稀疏表示，即α應(yīng)為零向量而β是一個稀疏向量；類似的，如果F是背景信號，則它不能由目標(biāo)過完備空時字典稀疏表示，即α應(yīng)為稀疏向量而β是一個零向量。因此，紅外圖像序列F可以被目標(biāo)和背景聯(lián)合的過完備空時字典稀疏表示。

2 自適應(yīng)形態(tài)過完備空時字典

構(gòu)建能表征紅外目標(biāo)圖像稀疏性的自適應(yīng)過完備空時字典D，并對其進行分類，分為能表征目標(biāo)成分的目標(biāo)過完備字典Dt和能表征背景噪聲成分的背景過完備字典Db.能否從字典D有效區(qū)分目標(biāo)過完備字典Dt和背景過完備字典Db則是提高目標(biāo)稀疏表示性能和小弱目標(biāo)檢測能力的關(guān)鍵。

2.1 自適應(yīng)形態(tài)成分空時域字典

采用K均值奇異值分解（K-SVD）從大量訓(xùn)練樣本中學(xué)習(xí)自適應(yīng)空時過完備字典D.空時過完備字典D的訓(xùn)練模型如下:

式中:‖·‖0和‖·‖2分別表示范數(shù)和范數(shù)。在殘差能量‖Dγ-F‖限制在一定條件下，圖像序列信號F可被過完備空時字典D中少量的原子及其稀疏系數(shù)g重構(gòu)。構(gòu)建自適應(yīng)空時過完備字典是一迭代過程，每次迭代包含稀疏編碼和字典更新兩個步驟。

1）稀疏編碼:首先給定初始字典D和紅外序列圖像F，采用正交匹配追蹤算法計算序列圖像在字典中的稀疏系數(shù)g，即求解

式中:ε是限定的逼近誤差。

2）字典D更新:每次更新字典D的一列，即更新一個原子dk.原子dk逼近序列圖形F的誤差為

圖1為紅外圖像及訓(xùn)練所得的過完備空時字典。圖1（a）為云層背景下的紅外圖像，圖1中方框所在中心位置為小弱目標(biāo)，淹沒于云層之中；圖1（b）～圖1（d）為一個空時原子?？蓮膱D1中看出，空時原子不僅能描述小弱目標(biāo)和背景雜波的形態(tài)特性，還能描述弱目標(biāo)和背景雜波的運動信息。然而，圖1（b）和圖1（c）表示運動目標(biāo)，而圖1（d）則表示運動的云層背景。表示目標(biāo)與背景的原子混雜在空時字典中，為信號檢測帶來困難。

2.2 過完備空時分類字典

小弱目標(biāo)的形態(tài)通常為點狀或者斑狀，其在圖像序列上占有若干少量像素，且幅度較均勻，目標(biāo)信號的幅度可采用點擴展函數(shù)描述。

式中:［dx（k），dy（k）］表示目標(biāo)的水平和垂直寬度；amax為目標(biāo)信號幅度；（x0，y0）為目標(biāo)所在位置；vx和vy分別為目標(biāo)在水平和垂直方向的速度。小弱運動目標(biāo)形態(tài)通常表現(xiàn)為運動的點狀或者斑狀，因此，它可由此空時三維高斯函數(shù)表示為

式中:I（i，j，k）為在k時刻、位于坐標(biāo)（x，y）處的像素值；Imax為目標(biāo)中心像素值，（i0，j0）為模型的中心坐標(biāo)；sx和sy分別表示水平和垂直方向的散布參數(shù)；ux和uy分別控制目標(biāo)在水平和垂直方向的速度。以樣本圖像的左上角為原點，（x0，y0）為目標(biāo)圖像的中心坐標(biāo)，隨機生成像素值的峰值Imax，調(diào)整均方差兩個參數(shù)，得到一系列紅外小目標(biāo)樣本圖像；然后，移動其中心位置到下一個像素點，重新生成新的樣本圖像，得到的每一副樣本圖像即為高斯過完備空域字典中的原子；最后，隨機生成多組不同取值的ux和uy，將高斯過完備空域字典中的原子以這些不同的速度在相應(yīng)的方向上移動這些原子，構(gòu)造出多個高斯過完備空域字典，這些高斯空域字典一起構(gòu)成高斯過完備空時字典Dgaussian.

圖1 空時冗余字典Fig.1 Spatial-temporal redundant dictionary

由于目標(biāo)通常為運動的點狀或者斑狀信號，表示目標(biāo)信號的空時原子可通過少量高斯空時字典稀疏重構(gòu)，相反，背景空時原子則難以通過少量高斯空時字典進行稀疏重構(gòu)，或者重構(gòu)后背景殘余能量大?；谶@一思想，對過完備空時字典D中原子dk在高斯空時字典Dgaussian中進行稀疏分解，通過殘余能量來判斷是否為目標(biāo)原子?？諘r原子dk在高斯過完備空時字典Dgaussian中稀疏分解表示為

式中:n為常數(shù)，比如3.

將r（dk）與閾值進行比較，判斷dk為背景空時原子還是目標(biāo)空時原子。如殘余能量r（dk）大于閾值，則判斷原子dk為背景原子，相反，dk則為目標(biāo)原子。

將每個空時原子進行判斷分類，分為背景過完備空時字典Db和目標(biāo)過完備空時字典Dt.然而，當(dāng)很大時，背景原子也能夠被高斯過完備字典稀疏表示，所以，應(yīng)該被限制在合適的范圍。

3 空時稀疏表示小弱目標(biāo)檢測

求解紅外圖像信號F在空時字典D=［DbDt］中的表示系數(shù)，即分別在背景空時字典Db和目標(biāo)空時字典Dt中的稀疏表示系數(shù)α和β.通過迭代的貪婪算法（匹配追蹤算法）求解它在一定容許誤差σ內(nèi)的逼近解:

如前所述，目標(biāo)信號由目標(biāo)空時域字典稀疏分解，但不能由背景空時域字典稀疏分解；因此，目標(biāo)在目標(biāo)字典中稀疏分解后，目標(biāo)信號就可由最大的幾個稀疏表示系數(shù)所對應(yīng)的目標(biāo)空時原子重構(gòu)，其殘余能量很??；但是，如果目標(biāo)信號由相同數(shù)目的背景空時原子重構(gòu)，其殘余能量卻很大。目標(biāo)空時原子重構(gòu)殘差為

式中:βi向量表示序列圖像信號F在目標(biāo)空時字典中稀疏表示系數(shù)；m為常數(shù)，比如5.

類似地，背景信號可由背景空時字典稀疏表示，但不能由目標(biāo)空時字典稀疏表示。背景信號可以由最大的幾個稀疏表示系數(shù)所對應(yīng)的空時背景原子重構(gòu)，其殘余能量很?。坏?，如果背景信號由相同數(shù)目的目標(biāo)空時原子重構(gòu)，其殘余能量卻很大。背景空時原子重構(gòu)殘差為

式中:αi向量代表序列圖像信號在背景空時字典的表示系數(shù)。比較圖像信號F在背景空時字典和目標(biāo)空時字典上稀疏分解后的殘余能量差

如果殘余能量差大于閾值η，則判斷F為目標(biāo)信號，相反，F(xiàn)為背景信號。

4 實驗結(jié)果與分析

圖2為目標(biāo)圖像序列在空時字典分解后的稀疏系數(shù)。圖2（a）表示的背景稀疏系數(shù)存在若干非零的稀疏系數(shù)；圖2（b）表示的目標(biāo)稀疏系數(shù)只存在幾個較大的分解系數(shù)，而其他原子分解系數(shù)基本為0.這表明目標(biāo)圖像序列在目標(biāo)過完備空時字典中的稀疏度較在背景空時字典中更稀疏。目標(biāo)圖像分別通過3個最大非零系數(shù)所對應(yīng)的目標(biāo)空時原子和背景空時域原子重構(gòu)，重構(gòu)后的歸一化殘差rt（F）和rb（F）分別為0.016和0.830，可看出，目標(biāo)空時域字典對該圖像是否含有目標(biāo)具有很強的區(qū)分度。

圖3為背景圖像序列及其在空時字典分解后稀疏表示系數(shù)。圖3（a）表示的背景稀疏系數(shù)只有3個非零稀疏系數(shù)；圖3（b）表示的目標(biāo)稀疏系數(shù)非零值很多。這表明背景圖像信號在目標(biāo)空時域字典中不稀疏，而在背景空時字典中則很稀疏。背景圖像分別通過圖3（a）和圖3（b）中最大3個非零系數(shù)所對應(yīng)的目標(biāo)空時原子和背景空時原子重構(gòu)，重構(gòu)后的歸一化殘差rt（F）和rb（F）分別為0.780和0.058，由此可以檢測出圖像是否含有背景雜波。

圖2 目標(biāo)在空時字典的稀疏系數(shù)Fig.2 Sparse coefficients of target image on spatial-temporal dictionary

對紅外圖像分別采用空域字典、空時字典和空時分類字典對目標(biāo)進行檢測，統(tǒng)計不同虛警概率下目標(biāo)檢測概率，即從接收機特征曲線，如圖4所示。從曲線可以看出，在虛警概率為0.1時，空時分類字典檢測概率最高大約為0.93，空時字典的檢測概率大約為0.80，空域字典的檢測概率大約為0.60.空域字典缺點在于目標(biāo)信號被分解為含有目標(biāo)和背景雜波的若干原子中，能量分散太廣泛。小弱目標(biāo)可在目標(biāo)空時字典中稀疏表示，但不能由背景字典稀疏表示；背景信號由背景字典稀疏表示，但不能由目標(biāo)字典稀疏表示，且目標(biāo)稀疏分解系數(shù)與背景稀疏分解系數(shù)的差異很大，比空時字典更易從雜波信號中檢測出目標(biāo)。

采用檢測概率和虛警概率兩個評價指標(biāo)對幀差法、背景減法、光流法和本文算法的評估結(jié)果如表1所示?？梢钥闯霰疚乃惴軌蚋吒怕实臋z測出目標(biāo)，同時保持較低的虛警率，性能明顯優(yōu)于其他3種檢測算法，但是本文算法的檢測時間相對較長，可以考慮采用更加快速的算法進行求解。

圖3 背景圖像在空時字典的稀疏系數(shù)Fig.3 Sparse coefficients of background on spatial-temporal redundant dictionary

圖4 稀疏表示檢測目標(biāo)概率Fig.4 Target detection probability

表1 4種算法性能對比Tab.1 Performance comparison of four algorithms

5 結(jié)論

本文提出了一種基于過完備空時稀疏表示的小弱目標(biāo)檢測方法，能描述出小弱目標(biāo)和背景雜波的運動特征和形態(tài)信息，減小檢測的不確定性。通過對空時字典進行分類，小弱運動目標(biāo)和背景雜波可分別在目標(biāo)空時字典和背景空時字典中稀疏表示，卻難以在相反的空時字典中稀疏表示，增強了目標(biāo)與背景的稀疏特征差異，提高小運動目標(biāo)的探測能力。

（

）

［1］ Bar-Shalom Y，Kirubarajan T，Lin X.Probabilistic data association techniques for target tracking with applications to sonar，radar and EO sensors［J］.Aerospace&Electronic Systems Magazine IEEE，2005，20（8）:37-56.

［2］ Liou R J，Azimi-Sadjadi M R.Dim target detection using high order correlation method［J］.IEEE Transactions on Aerospace and Electronic Systems，1993，29（3）:841-856.

［3］ Orlando D，Venturino L，Lops M，et al.Track-before-detectstrategies for STAP radars［J］.IEEE Transactions on Signal Processing，2010，58（2）:933-938.

［4］ Li Z，Qi L，Li W，et al.Track initiation for dim small moving infrared target based on spatial-temporal hypothesis testing［J］.Journal of Infrared，Millimeter and Terahertz Waves，2009，30（5）:513-525.

［5］ Zeng M，Li J，Peng Z.The design of top-hat morphological filter and application to infrared target detection［J］.Infrared Physics& Technology，2006，48（1）:67-76.

［6］ Cao Y，Liu R M，Yang J.Small target detection using two-dimensional least mean square（TDLMS）filter based on neighborhood analysis［J］.International Journal of Infrared and Millimeter Waves，2008，29（2）:188-200.

［7］ Ting W，Yang S.Weak and small infrared target automatic detection based on wavelet transform［C］//International Symposium on Intelligent Information Technology Application.Piscataway，NJ，US:IEEE，2008:697-701.

［8］ Davidson G，Griffiths H D.Wavelet detection scheme for small targets in sea clutter［J］.Electronics letters，2002，38（19）: 1128-1130.

［9］ 宋相法，焦李成.基于稀疏表示及光譜信息的高光譜遙感圖像分類［J］.電子與信息學(xué)報，2012，34（2）:268-272. SONG Xiang-fa，JIAO Li-cheng.Classification of hyperspectral remote sensing image based on sparse representation and spectral information［J］.Journal of Electronics&Information Technology，2012，34（2）:268-272.（in Chinese）

［10］ Buzzi S，Lops M，Venturino L.Track-before-detect procedures for early detection of moving target from airborne radars［J］. IEEE Transaction on Aerospace and Electronic Systems，2005，41（3）:937-954.

［11］ 劉建軍，吳澤彬，韋志輝，等.基于空間相關(guān)性約束稀疏表示的高光譜圖像分類［J］.電子與信息學(xué)報，2012，34（11）: 2666-2671. LIU Jian-jun，WU Ze-bin，WEI Zhi-hui，et al.Spatial correlation constrained sparse representation for hyperspectral image classification［J］.Journal of Electronics&Information Technology，2012，34（11）:2666-2671.（in Chinese）

［12］ Carlson B D，Evans E D，Wilson S L.Search radar detection and track with the Hough transform.I.system concept［J］.IEEE Transactions on Aerospace&Electronic Systems，1994，30（1）: 102-108.

［13］ Chen H，Bar-shalom Y，Pattipati K R，et al，MDL approach for multiple low observable targettrack initiation［J］.IEEE Transaction on Aerospace and Electronic Systems，2003，39（3）:862-882.

［14］ Johnston L A，Krishnamurthy V.Performance analysis of a dynamic programming track before detect algorithm［J］.IEEE Transactions onAerospace&Electronic Systems，2002，38（1）: 228-242.

［15］ Blackman S S.Multiple hypothesis tracking for multiple target tracking［J］.Aerospace&Electronic Systems Magazine IEEE，2004，19（1）:5-18.

［16］ Reed I S，Gagliardi R M，Stotts L B.Optical moving target detection with 3-D matched filtering［J］.IEEE Transactions onAerospace&Electronic Systems，1988，24（4）:327-336.

［17］ Zhang J D，Zhu D Y，Zhang G.Adaptive compressed sensing radar oriented toward cognitive detection in dynamic sparse target scene［J］.IEEE Transaction on Signal Processing，2012，60（4）: 1718-1729.

［18］ 趙佳佳，唐崢遠，楊杰，等.基于圖像稀疏表示的紅外小目標(biāo)檢測算法［J］.紅外與毫米波學(xué)報，2011，30（2）:156-161. ZHAO Jia-jia，TANG Zheng-yuan，YANG Jie，et al.Infrared small target detection based on image sparse representation［J］. Journal of Infrared and Millimeter Waves，2011，30（2）:156-161.（in Chinese）

［19］ Chen Y，Nasrabadi N M，Tran T D，et al.Sparse representation for target detection in hyperspectral imagery［J］.IEEE Journal of Selected Topics in Signal Processing，2011，5（3）:629-640.

Dim Moving Target Detection Algorithm Based on Spatial-temporal Sparse Representation

LI Zheng-zhou1，HOU Qian1，DAI Zhen1，F(xiàn)U Hong-xia1，GE Feng-zeng1，JIN Gang2
（1.College of Communication Engineering，Chongqing University，Chongqing 400044，China；2.China Aerodynamics Research and Development Center，Mianyang 621000，Sichuan，China）

A dim moving target detection algorithm based on over-complete spatial-temporal dictionary and sparse representation is proposed.A spatial-temporal adaptive morphological over-complete dictionary is trained and constructed according to infrared image sequence.It can represent the motion information and morphological feature of target and background clutter.The spatial-temporal morphological over-complete dictionary is subdivided into two categories:target over-complete spatial-temporal dictionary for describing moving target，and background over-complete spatial-temporal dictionary for embedding background.The criteria adopted to distinguish the target spatial-temporal redundant dictionary from the background spatial-temporal redundant dictionary is that the atom in target over-complete spatial-temporal dictionary could be decomposed more sparsely over Gaussian over-complete spatial-temporal dictionary.Subsequently，the image sequence is decomposed on the target and background over-complete spatial-temporal dictionaries，respectively.The dim moving target and background clutter can be sparsely decomposed on their corresponding over-complete spatial-temporal dictionary，yet it couldn't be sparsely decomposedon their background over-complete spatial-temporal dictionary.Therefore，the target and background clutter would be reconstructed effectively by prescribed number of atoms with maximum sparse coefficients in their corresponding over-complete spatial-temporal dictionary，and their residuals would differ so visibly to distinguish target from background clutter.The results show that the proposed approach not only could improve the sparsity more efficiently for dim target image sequence，but also could improve the performance of small target detection.

information processing technology；dim target detection；spatial-temporal redundant dictionary；target spatial-temporal redundant dictionary；background spatial-temporal redundant dictionary；signal sparse reconstruction

TP391.4

A

1000-1093（2015）07-1273-07

10.3969/j.issn.1000-1093.2015.07.016

2014-11-14

國家自然科學(xué)基金項目（61071191）；中國科學(xué)院光束控制重點實驗室基金項目（2014LBC005）；中國博士后基金項目（2014M550455）；重慶博士后科研項目特別基金項目（XM201489）；中央高?；究蒲袠I(yè)務(wù)費專項資金項目（106112013CDJZR160007、106112014CDJZR165502）；2013年重慶高校創(chuàng)新團隊建設(shè)計劃項目（KJTD201331）

李正周（1974—），男，教授，博士生導(dǎo)師。E-mail:lizhengzhou@cqu.edu.cn