王宜賢，石德樂，楊　寧

（中國空間技術(shù)研究院513所，山東煙臺(tái)264003）

幀差法和Mean shift算法融合的高速無人機(jī)目標(biāo)跟蹤

王宜賢，石德樂，楊寧

（中國空間技術(shù)研究院513所，山東煙臺(tái)264003）

為實(shí)時(shí)跟蹤高速飛行無人機(jī)，圖像跟蹤算法必須滿足快速性和準(zhǔn)確性要求。文章給出一個(gè)融合算法，將幀差法和Mean shift算法的優(yōu)勢結(jié)合起來。2個(gè)算法平行運(yùn)行，差幀法實(shí)現(xiàn)快速跟蹤，Mean shift算法則用于對幀差法結(jié)果進(jìn)行準(zhǔn)確度修正。還利用Kalman濾波技術(shù)對計(jì)算周期內(nèi)的無人機(jī)運(yùn)動(dòng)位移進(jìn)行補(bǔ)償，進(jìn)一步提高實(shí)時(shí)跟蹤的準(zhǔn)確性，并給出Matlab仿真例子驗(yàn)證本文方法的有效性。

無人機(jī)；目標(biāo)跟蹤；幀差法；Mean shift算法；融合算法

1　問題的提出

對處于飛行狀態(tài)的無人機(jī)進(jìn)行無線充能需要2個(gè)環(huán)節(jié)：首先，利用成像器件捕捉無人機(jī)圖像序列，通過圖像處理技術(shù)檢測無人機(jī)目標(biāo)并獲得目標(biāo)坐標(biāo)；其次，將無人機(jī)坐標(biāo)信息傳給姿態(tài)定位及隨動(dòng)跟蹤系統(tǒng)實(shí)施充能動(dòng)作。如何利用圖像處理技術(shù)獲得準(zhǔn)確的無人機(jī)實(shí)時(shí)坐標(biāo)，這是無人機(jī)無線充能關(guān)鍵的一環(huán)。

在圖像處理領(lǐng)域，對動(dòng)態(tài)目標(biāo)的檢測、跟蹤一直是一個(gè)備受關(guān)注的問題［1-8］。基于不同的理論模型，不同的跟蹤算法被提出，主流的算法有幀差法、背景差分法、模板匹配法、光流法、Mean shift算法、Camshift算法等［9-15］。這些方法中，不同的算法有各自的優(yōu)勢，準(zhǔn)確度高的算法計(jì)算周期長，實(shí)時(shí)性較差，而快速的算法則往往跟蹤精度較差。沒有一種算法能兼?zhèn)淙克惴ǖ膬?yōu)點(diǎn)。

近年來，為了解決復(fù)雜的實(shí)際問題，人們開始將不同的圖像算法融合起來，通過不同算法的優(yōu)勢組合，獲得更好的結(jié)果［16-19］。文獻(xiàn)［16］將幀差法與Mean shift算法進(jìn)行融合，利用幀差法檢測運(yùn)動(dòng)目標(biāo)，確定初始窗口及設(shè)定跟蹤模板，然后用Mean shift算法進(jìn)行目標(biāo)跟蹤。文獻(xiàn)［17］將幀差法和光流法融合，使用幀差法對圖像進(jìn)行分割，然后利用光流場算法對運(yùn)動(dòng)目標(biāo)進(jìn)行檢測。文獻(xiàn)［18］將幀差法和背景差分法進(jìn)行融合用于運(yùn)動(dòng)目標(biāo)檢測，而文獻(xiàn)［19］將幀差法和分水嶺算法進(jìn)行融合。

與跟蹤快速無人機(jī)的要求相比，文獻(xiàn)［16-19］給出的融合算法不能滿足快速性的要求，比如，文獻(xiàn)［17］中的光流法運(yùn)算復(fù)雜，需要特定的硬件支持。文獻(xiàn)［18］給出的算法雖然運(yùn)算量小，可以實(shí)時(shí)運(yùn)算，但是準(zhǔn)確性不能滿足要求。

為實(shí)現(xiàn)高速無人機(jī)實(shí)時(shí)跟蹤的目的，本文提出一個(gè)將幀差法和Mean shift算法融合的算法。其思路是并行運(yùn)算幀差法和Mean shift算法，利用幀差法對目標(biāo)進(jìn)行快速檢測和定位，然后用Mean shift算法對幀差法的結(jié)果進(jìn)行精度修正。為進(jìn)一步補(bǔ)償無人機(jī)在計(jì)算周期內(nèi)飛行所導(dǎo)致的跟蹤計(jì)算結(jié)果的偏差，將修正的結(jié)果看作無人機(jī)飛行過程的觀測量，建立Kalman濾波方程對無人機(jī)位置進(jìn)行實(shí)時(shí)預(yù)報(bào)，從而實(shí)現(xiàn)對高速無人機(jī)的實(shí)時(shí)跟蹤。

2　幀差法和Mean shift算法的融合

2.1幀差法

幀差法又叫閾值幀間差分法（TIFD），是運(yùn)動(dòng)目標(biāo)檢測中使用最多的算法。它的突出特點(diǎn)是算法簡單、計(jì)算量小、實(shí)時(shí)性高，在大多數(shù)情況下檢測效果較好。幀差法基本原理是設(shè)定一個(gè)閾值，然后將前后兩幀圖像對應(yīng)像素點(diǎn)的灰度相減。對應(yīng)灰度差小于閾值的像素點(diǎn)標(biāo)記為靜止的，對應(yīng)灰度值差大于等于閾值的像素點(diǎn)標(biāo)記為運(yùn)動(dòng)的。所有標(biāo)記為運(yùn)動(dòng)的像素點(diǎn)可以形成運(yùn)動(dòng)目標(biāo)的輪廓。

以T為閾值，It-1（x）、It（x）分別為相鄰2個(gè)時(shí)刻采集的圖像在像素點(diǎn)x處的灰度值，做差分并得到差分圖像：

利用矩法計(jì)算目標(biāo)輪廓的質(zhì)心［20］。令：

質(zhì)心的計(jì)算公式為：

式中，yf被看作是無人機(jī)的位置坐標(biāo)。

2.2Mean shift算法

Mean shift算法用概率密度來描述目標(biāo)和候選目標(biāo)的特征，定義兩者間的相似系數(shù)，度量它們間的匹配程度，算法的目標(biāo)是尋找候選目標(biāo)，使其與目標(biāo)有最大的相似系數(shù)。Mean shift算法具有很高的穩(wěn)定性，它能夠適應(yīng)目標(biāo)的形狀、大小的連續(xù)變化，而且抗干擾能力強(qiáng)。為獲得準(zhǔn)確度高的結(jié)果，Mean shift算法需要迭代多次，這使得計(jì)算周期變長，影響實(shí)時(shí)性。

Mean shift算法跟蹤算法主要包括3個(gè)步驟：目標(biāo)模型和候選模型的建立、相似性度量、迭代計(jì)算。

選定目標(biāo)區(qū)域?yàn)榘瑹o人機(jī)的一個(gè)矩形區(qū)域，區(qū)域中心選為無人機(jī)圖像的質(zhì)心。設(shè)該區(qū)域內(nèi)有n個(gè)像素點(diǎn)，將圖像的灰度空間劃分m個(gè)不同的區(qū)間，分別以1，2，…，m索引，則目標(biāo)模型被表示為灰度索引u的直方分布密度為：

式（4）中，qu被定義為：

式（5）中：C是歸一化參數(shù)；k（·）是核函數(shù)的輪廓函數(shù)；h為帶寬；b（xi）表示像素點(diǎn)xi所屬區(qū)間的索引；x0是待跟蹤目標(biāo)的中心，選為無人機(jī)的圖像質(zhì)心；δ是Kronecker delta函數(shù)。

候選模型的特征：候選模型由像素點(diǎn)y點(diǎn)確定，它仍被描述為灰度索引u的直方分布密度為：

其中，pu（y）被定義為：

式（7）中：y表示候選目標(biāo)的位置；其他參數(shù)與目標(biāo)模型相同。

這樣，搜尋相似系數(shù)的極值問題就是尋求ρ（y）的優(yōu)化問題。這個(gè)搜索y的過程可以通過迭代進(jìn)行，從上一步y(tǒng)0到下一步y(tǒng)1的迭代公式為：

為保證盡可能獲得相似系數(shù)的優(yōu)化值，一般需迭代30次，迭代結(jié)束時(shí)輸出的結(jié)果ym即被看作無人機(jī)的坐標(biāo)位置。

2.3算法融合

幀差法雖然計(jì)算速度快，但相對于真實(shí)目標(biāo)，這一方法檢測到的目標(biāo)輪廓往往范圍較大，算出的質(zhì)心與真實(shí)目標(biāo)位置會(huì)有較大的偏差。Mean shift算法雖然跟蹤精度比較好，但需要做大量迭代運(yùn)算，這樣就難以保證實(shí)時(shí)性。為將幀差法的快速性和Mean shift算法的精度優(yōu)勢結(jié)合起來，對2個(gè)算法進(jìn)行融合。方法是同時(shí)并行運(yùn)行2個(gè)算法，用Mean shift算法的結(jié)果修正幀差法的結(jié)果，使修正后的結(jié)果盡可能接近Mean shift算法的結(jié)果。由于在一個(gè)Mean shift算法計(jì)算周期內(nèi)，幀差法會(huì)執(zhí)行多個(gè)計(jì)算周期，用上一Mean shift算法的結(jié)果修正當(dāng)前Mean shift計(jì)算周期內(nèi)幀差法的所有結(jié)果。

假設(shè)成像部件采集圖像的周期為T，在kT（k=0，1，2，…）時(shí)刻采集的無人機(jī)圖像組成一個(gè)圖像序列，記為G0，G1，G2，…，Gk，…；幀差法記為Pr，計(jì)算周期為 Tr＜T；Mean shift算法記為 Ph，計(jì)算周期為Th＜NT。

幀差法Pr對kT時(shí)刻采集的圖像Gk進(jìn)行處理，在kT+Tr時(shí)刻得到無人機(jī)在kT時(shí)刻的位置坐標(biāo)，結(jié)果記為這一過程表示為：

Mean shift算法Ph的計(jì)算周期看作NT，它只能對kNT時(shí)刻采集的圖像GkN處理，在kNT+Th時(shí)刻計(jì)算出無人機(jī)在kNT時(shí)刻的位置坐標(biāo)，得到結(jié)果記為，這一過程表示為：

在kNT時(shí)刻，幀差法和Mean shift算法處理同一幅圖像GkN，結(jié)果分別為和，其差即為Mean shift算法相對幀差法的改進(jìn)值。這個(gè)差反映2個(gè)算法間算法質(zhì)量的差別，在背景近似不變的條件下，差也近似不變。

因Mean shift算法在（k+1）NT時(shí)刻完成，這個(gè)差在（k+1）NT 得到，所以以修 正時(shí)段內(nèi)幀差法得到的N個(gè)結(jié)果，其中Ak+1為修正系數(shù)矩陣，即：

Ak+1的選擇應(yīng)使修正后的結(jié)果盡可能接近Mean shift算法的結(jié)果。在Mean shift算法計(jì)算周期內(nèi)，沒有新的Mean shift算法的結(jié)果用于接近程度的評(píng)價(jià)，這種評(píng)價(jià)只能建立在已有的計(jì)算結(jié)果上。按如下方式選擇 Ak+1：假設(shè)在（k-1）NT時(shí)刻獲得后，用修正［kNT，（k+1）NT-T］時(shí)間段內(nèi)幀差法的計(jì)算結(jié)果，那么Mean shift算法在（k+1）NT時(shí)刻得出結(jié)果后，“理想的修正值”應(yīng)滿足，即A應(yīng)滿足：

將上一時(shí)段的“理想的系數(shù)矩陣A”取做當(dāng)前時(shí)間段的Ak+1，它可表示為：

3　Kalman濾波

在執(zhí)行圖像處理期間，處于飛行狀態(tài)的無人機(jī)會(huì)飛行一段距離，這樣計(jì)算出的坐標(biāo)不是無人機(jī)的實(shí)時(shí)位置。下面引入Kalman濾波方法消除這種計(jì)算時(shí)延導(dǎo)致的跟蹤誤差。首先為高速無人機(jī)系統(tǒng)建立Kalman濾波所需的狀態(tài)方程和觀測方程。

在無人機(jī)目標(biāo)跟蹤過程中，視頻圖像采集時(shí)間間隔較短，無人機(jī)在一個(gè)圖像采集周期內(nèi)的運(yùn)動(dòng)可看作是勻速地。

定義nT時(shí)刻卡爾曼濾波的系統(tǒng)狀態(tài)為Xn=（xn-1，yn-1，xn，yn）T，（xn-1，yn-1）和（xn，yn）分別為無人機(jī)在（n-1）T時(shí)刻和nT時(shí)刻所在的位置。設(shè)為 nT時(shí)刻無人機(jī)在x軸和y軸方向的運(yùn)動(dòng)速度，則有：

則對于n+1時(shí)刻的無人機(jī)的坐標(biāo)有：

于是，可以建立無人機(jī)的狀態(tài)轉(zhuǎn)移方程：

式中，Wn為狀態(tài)轉(zhuǎn)移過程中存在的干擾白噪聲。

式中，Vn為觀測過程中存在的干擾白噪聲。

關(guān)于轉(zhuǎn)移方程和觀測方程中的干擾白噪聲，假定它們滿足下面統(tǒng)計(jì)特性：

式（19）中：Qk是白噪聲Wk的4階協(xié)方差矩陣；Rk是白噪聲Vk的2階協(xié)方差陣；δkj是Kronecker函數(shù)。

根據(jù)建立的狀態(tài)方程和觀測方程，可以執(zhí)行標(biāo)準(zhǔn)的Kalman濾波步驟。由關(guān)于（n-1）T時(shí)刻觀測值，可以得到當(dāng)前nT時(shí)刻的預(yù)報(bào)值，這個(gè)預(yù)報(bào)值可以作為無人機(jī)的當(dāng)前時(shí)刻的位置，傳遞給隨動(dòng)系統(tǒng)。

4　仿真實(shí)驗(yàn)

根據(jù)問題實(shí)際背景的要求，成像設(shè)備采集的無人機(jī)飛行圖像的分辨率設(shè)為480像素×600像素，圖像采集頻率為48Hz，周期T≈20.8 ms，無人機(jī)坐標(biāo)定位的誤差小于1個(gè)像素。為保證對圖像跟蹤的實(shí)時(shí)性，要求跟蹤算法的處理周期在20 ms以內(nèi)。

幀差法的平均計(jì)算周期約為Tr=13.2 ms＜T，小于圖像采集周期，Mean shift算法的平均運(yùn)算周期約為Th=126.1 ms＜7T。

利用成像設(shè)備在真實(shí)場景下采集無人機(jī)飛行的圖像序列，在Matlab環(huán)境下對融合算法進(jìn)行仿真驗(yàn)證，得到的結(jié)果在圖1給出。圖1中的方框表示跟蹤波門，它以Kalman預(yù)報(bào)值為中心。

圖1　無人機(jī)融合算法跟蹤仿真Fig.1 Simulation of UAV tracking with fusion algorithm

5　結(jié)語

在高速無人機(jī)圖像跟蹤問題中，跟蹤算法需要同時(shí)滿足快速性和準(zhǔn)確性的要求，這2個(gè)目標(biāo)難以用一個(gè)現(xiàn)有常規(guī)的算法單獨(dú)實(shí)現(xiàn)。

本文給出一個(gè)融合算法，即將幀差法和Mean shift算法的優(yōu)勢結(jié)合起來。2個(gè)算法平行運(yùn)行，差幀法實(shí)現(xiàn)快速跟蹤，Mean shift算法用于對幀差法結(jié)果的進(jìn)行精度修正。同時(shí)，本文還利用Kalman濾波技術(shù)對計(jì)算周期內(nèi)的無人機(jī)運(yùn)動(dòng)位移進(jìn)行補(bǔ)償，進(jìn)一步提高實(shí)時(shí)跟蹤的準(zhǔn)確性。雖然本文提出的融合算法能提高幀差法的準(zhǔn)確性，實(shí)現(xiàn)對無人機(jī)快速和準(zhǔn)確跟蹤，但也存在不足。當(dāng)背景出現(xiàn)較強(qiáng)干擾時(shí)，幀差法會(huì)將干擾源處理成無人機(jī)輪廓使計(jì)算結(jié)果出現(xiàn)偏差，這種計(jì)算偏差無法用本文的融合算法修正。

［1］崔家瑞，陳鵬，胡廣大.空間站脈沖編碼被動(dòng)式交會(huì)對接合作目標(biāo)三維跟蹤［J］.宇航學(xué)報(bào)，2016，37（1）：102-111. CUI JIARUI，CHEN PENG，HU GUANGDA.Cooperative target 3D tracking for space station pulse coded passive rendezvous and docking［J］.Journal of Astronautics，2016，37（1）：102-111.（in Chinese）

［2］DOU JIANFANG，LI JIANXUN.Robust visual tracking base on adaptively multi-feature fusion and particle filter ［J］.Optik，2014，125：1680-1686.

［3］鄭遠(yuǎn)力，胡志坤.基于滑動(dòng)掃描框的高速物體的圖像實(shí)時(shí)跟蹤算法［J］.計(jì)算機(jī)科學(xué)，2015，42（10）：287-291. ZHENG YUANLI，HU ZHIKUN.Real-time tracking algorithm for fast target based on dynamical scanning boxes ［J］.Computer Science，2015，42（10）：287-291.（in Chinese）

［4］GAO MINGLIANG，SHEN JIN，YIN LIJU.A novel visual tracking method using bat algorithm［J］.Neurocomputing，2016，177：612-619.

［5］胡晉銘，張魯濱，任瑞，等.飛行試驗(yàn)復(fù)雜背景下的圖像跟蹤技術(shù)［J］.海軍航空工程學(xué)院學(xué)報(bào)，2010，25（1）：11-14. HU JINMING，ZHANG LUBIN，REN RUI，et al.Video tracking technology for the flight test under complex background［J］.Journal of Naval Aeronautical and Astronautical University，2010，25（1）：11-14.（in Chinese）

［6］GIUSEPPE GUIDO，ALESSANDRO VITALE，F(xiàn)RANK F SACCOMANNO.Vehicle tracking system based on videotaping data［J］.Procedia-Social and Behavioral Sciences，2014，111：1123-1132.

［7］史延科，高曉光.圖像目標(biāo)跟蹤及其在UCAV中的應(yīng)用［J］.飛行力學(xué)，2005，23（1）：79-81. SHI YANKE，GAO XIAOGUANG.The image target tracking and its applications in UCAV［J］.Flight Dynamics，2005，23（1）：79-81.（in Chinese）

［8］江山，張銳，韓廣良，等.復(fù)雜背景灰度圖像下的多特征融合運(yùn)動(dòng)目標(biāo)跟蹤［J］.中國光學(xué)，2016，9（3）：320-328. JIANG SHAN，ZHANG RUI，HAN GUANGLIANG，et al.Moving object tracking based on multi-feature fusion in the complex background gray image［J］.Chinese Optics，2016，9（3）：320-328.（in Chinese）

［9］李超，于飛，康曉軍.一種模板匹配的快速實(shí)現(xiàn)方法［J］.航天返回與遙感，2016，37（1）：63-70. LI CHAO，YU FEI，KANG XIAOJUN.A fast method of template matching［J］.Spacecraft Recovery&Remote Sensing，2016，37（1）：63-70.（in Chinese）

［10］崔陸月，寧媛，袁浩.運(yùn)動(dòng)物體的跟蹤檢測算法的研究［J］.工業(yè)控制計(jì)算機(jī)，2015，28（8）：88-90. CUI LUYUE，NING YUAN，YUAN HAO.Study of detection algorithm of moving objects［J］.Industrial Control Computer，2015，28（8）：88-90.（in Chinese）

［11］顧幸方，茅耀斌，李秋潔.基于Mean shift的視覺目標(biāo)跟蹤算法綜述［J］.計(jì)算機(jī)科學(xué)，2012，39（12）：16-24. GU XINGFANG，MAO YAOBIN，LI QIUJIE.Survey on visual tracking algorithms based on mean shift［J］.Computer Science，2012，39（12）：16-24.（in Chinese）

［12］MUSAB CO?KUN，SENCER üNAL.Implementation of tracking of a moving object based on camshift approach with a UAV［J］.Procedia Technology，2016，22：556-561.

［13］邱道尹，張文靜，顧波，等.幀差法在運(yùn)動(dòng)目標(biāo)實(shí)時(shí)跟蹤中的應(yīng)用［J］.華北水利水電學(xué)院學(xué)報(bào)，2009，30（3）：45-46. QIU DAOYIN，ZHANG WENJING，GU BO，et al.Appli-cation of frame difference methods in real-time moving target tracking［J］.Journal of North China Institute of Water Conservancy and Hydroelectric Power，2009，30（3）：45-46.（in Chinese）

［14］劉洪彬，常發(fā)亮.權(quán)重系數(shù)自適應(yīng)光流法運(yùn)動(dòng)目標(biāo)檢測［J］.光學(xué)精密工程，2016，24（2）：460-468. LIU HONGBIN，CHANG FALIANG.Moving object detection based on adaptive by optical flow method weight coefficient［J］.Optics and Precision Engineering，2016，24（2）：460-468.（in Chinese）

［15］李培華.一種改進(jìn)的Mean Shift跟蹤算法［J］.自動(dòng)化學(xué)報(bào)，2007，33（4）：347-354. LI PEIHUA.An improved mean shift algorithm for object tracking［J］.Acta Automatica Sinica，2007，33（4）：347-354.（in Chinese）

［16］鄧文浩，唐立才，張南鋒，等.基于改進(jìn)的幀差法和Mean shift結(jié)合的運(yùn)動(dòng)目標(biāo)自動(dòng)檢測與跟蹤［J］.現(xiàn)代電子技術(shù)，2016，39（4）：108-111. DENG WENHAO，TANG LICAI，ZHANG NANFENG，et al.Moving target automatic detection and tracking based on integration of modified frame difference method and Mean shift algorithm［J］.Modern Electronics Technique，2016，39（4）：108-111.（in Chinese）

［17］李建，藍(lán)金輝，李杰.一種新型快速運(yùn)動(dòng)目標(biāo)檢測算法［J］.中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，44（3）：978-984. LI JIAN，LAN JINHUI，LI JIE.A novel fast moving target detection method［J］.Journal of Central South University：Science and Technology，2013，44（3）：978-984.（in Chinese）

［18］屈晶晶，辛云宏.連續(xù)幀間差分與背景差分相融合的運(yùn)動(dòng)目標(biāo)檢測方法［J］.光子學(xué)報(bào)，2014，43（7）：213-220. QU JINGJING，XIN YUNHONG.Combined continuous frame difference with background difference method for moving object detection［J］.Acta Photonica Sinica，2014，43（7）：213-220.（in Chinese）

［19］黃應(yīng)清，趙哲，閆興鵬，等.改進(jìn)幀差與分水嶺算法融合的運(yùn)動(dòng)目標(biāo)檢測法［J］.激光與紅外，2014，44（12）：1369-1373. HUANG YINGQING，ZHAO ZHE，YAN XINGPENG，et al.Moving target detection based on improved frame difference and watershed algorithm［J］.Laser and Infrared，2014，44（12）：1369-1373.（in Chinese）

［20］宋克歐，黃鳳崗，蘭曉亭.二值圖像目標(biāo)質(zhì)心快速搜索跟蹤算法［J］.模式識(shí)別與人工智能，1998，11（2）：161-168. SONG KEOU，HUANG FENGGAN，LAN XIAOTING. A fast algorithm for searching and tracking object centroids in binary images［J］.Pattern Recognition and Artificial Intelligence，1998，11（2）：161-168.（in Chinese）

Fusion Algorithm of Frame Difference and Mean Shift for High-Speed UAVs’Tracking

WANG Yixian，SHI Dele，YANG Ning
（Shandong Institute of Aerospace Electronic Technology，Yantai Shandong 264000，China）

In order to track high-speed UAVs in real time，the visual tracking algorithms need to satisfy the requirements of accuracy and rapidity.In this paper，a fusion algorithm was presented that combined the advantages of two typical algorithms，frame difference and Mean shift.The two algorithms were run simultaneously the frame difference was used for rapid tracking，while Mean shift was used for improving the precision of the tracking of frame difference.To get more precise tracking，the Kalman filter was used to compensate the displacement of UAVs moving during the implement of the tracking algorithm.The effectiveness of fusion algorithm was illustrated by a simulation example in Matlab environment.

UAV；object tracking；frame difference；Mean shift algorithm；fusion algorithm

TN949.2；V249.1

A

1673-1522（2016）04-0437-05

10.7682/j.issn.1673-1522.2016.04.006

2016-05-06；

2016-06-10

國家高技術(shù)研究發(fā)展計(jì)劃（863計(jì)劃）基金資助項(xiàng)目（2012AA120605）

王宜賢（1982-），男，工程師，大學(xué)。