李瑞東,孫協(xié)昌,李 勐
?
空間目標(biāo)紅外特征提取與識別技術(shù)
李瑞東,孫協(xié)昌,李 勐
(北京衛(wèi)星信息工程研究所,北京 100086)
針對空間目標(biāo)紅外識別中典型的類內(nèi)變化大、類間變化小等問題,提出了一種用支持向量機(jī)(Support Vector Machine,SVM)分類器結(jié)合局部二值模式(Local Binary Pattern,LBP)直方圖特征和灰度直方圖特征的空間目標(biāo)紅外特征提取與識別方法。以國內(nèi)某衛(wèi)星和國外某衛(wèi)星為研究對象,提取它們紅外圖像的LBP直方圖特征以及灰度直方圖特征;使用紅外仿真軟件生成兩個目標(biāo)在不同姿態(tài)、不同分辨率下的樣本圖,并分成兩部分,分別作為SVM分類器的訓(xùn)練集和測試集。實(shí)驗(yàn)結(jié)果表明,LBP直方圖特征和灰度直方圖特征均能夠以較高的準(zhǔn)確率對空間目標(biāo)進(jìn)行識別,且其識別效能與目標(biāo)紅外圖像的分辨率以及SVM核函數(shù)有關(guān)。
紅外識別;特征提取;空間目標(biāo)
自從1957年發(fā)射第一顆人造地球衛(wèi)星以來,人類在航天領(lǐng)域的活動也逐漸頻繁起來,隨之而來的則是不斷增加的在軌航天器以及空間碎片,造成空間環(huán)境日益惡化。惡劣的空間環(huán)境不僅對在軌航天器構(gòu)成直接威脅[1],甚至對地面上人類的活動也形成了潛在的威脅。據(jù)不完全統(tǒng)計(jì),目前在環(huán)繞地球的軌道上總共有超過6000噸的人造空間物體,其中能夠工作的航天器僅占7%,其余的均為空間碎片,且碎片的數(shù)量仍然逐年增加。
為了保障衛(wèi)星等航天器的安全,急需開展空間目標(biāo)監(jiān)視識別技術(shù)研究。當(dāng)前已有許多國家取得較好的成效,美國、俄羅斯、英國、法國、德國和日本等國家都已經(jīng)具有探測空間目標(biāo)的能力,特別是美國己經(jīng)建立了較完善的空間目標(biāo)監(jiān)視識別系統(tǒng)。
紅外探測經(jīng)過幾十年的發(fā)展,已經(jīng)成為一種非常重要的光學(xué)探測手段。熱輻射是紅外輻射的來源主要,其本質(zhì)是物體內(nèi)部分子和原子的運(yùn)動所產(chǎn)生的,而且溫度越高,這種微觀運(yùn)動越劇烈,所產(chǎn)生的紅外輻射也越強(qiáng),因此紅外探測對于目標(biāo)與背景之間的溫差具有較高的敏感性[2];其次,物體的紅外輻射能量不完全來自于環(huán)境的入射能量,其散發(fā)出的能量也不局限于反射能量,因此無論在白天還是夜間都可以進(jìn)行目標(biāo)紅外探測[3],可全天時工作;與有源的電磁探測手段相比[4],紅外探測手段本身還具有隱蔽性的特點(diǎn)。因此,從目標(biāo)的紅外特征入手,研究空間目標(biāo)紅外特征提取識別技術(shù),對空間目標(biāo)監(jiān)視識別系統(tǒng)的構(gòu)建具有重要意義。
目標(biāo)特征提取是目標(biāo)識別的前提和保障,其基本任務(wù)是從大量的觀測數(shù)據(jù)中提取最有效的特征,為空間目標(biāo)識別、特性編目、參考圖制備和交會分析提供信息基礎(chǔ)。
目前計(jì)算機(jī)視覺領(lǐng)域發(fā)展了一些用于目標(biāo)檢測的、性能優(yōu)良的特征描述子,比如Haar特征、LBP特征、SIFT描述子等。
1)Haar特征
Haar特征是一種簡單而有效的特征,相對比較簡單,計(jì)算速度較快,對于有固定結(jié)構(gòu)的目標(biāo)檢測有非常好的效果。Haar特征模板由白色和黑色兩種矩形組成,其特征值為白色矩形在圖像子窗口對應(yīng)區(qū)域的灰度值之和減去黑色矩形在圖像子窗口對應(yīng)區(qū)域的灰度值之和,反映了圖像局部的灰度變化情況。圖像的結(jié)構(gòu)位置以及輪廓線或關(guān)鍵點(diǎn)都可以通過這種對簡單矩形特征的計(jì)算進(jìn)行提取。
常用的Haar特征包括4類特征,分別為邊緣特征、線型特征、中心特征和對角特征。由于矩形模板可以在圖像中任意位置、任意尺度計(jì)算,因此可以認(rèn)為Haar特征主要由3個因素決定,分別是矩形模板類別、起始位置和矩形大小。也正是由于Haar特征這種非常靈活的特性,因而可以在圖像中提取大量的Haar特征。此外,還可以通過積分圖實(shí)現(xiàn)Haar特征的快速計(jì)算。
2)LBP算子
LBP特征不僅能描述一些局部微小特征,包括亮點(diǎn)、暗點(diǎn)、穩(wěn)定區(qū)域及紋理細(xì)節(jié)等,而且能反映出這些特征的分布情況。另外,LBP具有任意單調(diào)變換不變性和圖像旋轉(zhuǎn)不變性。
3)SIFT特征
SIFT算法由D.G.Lowe 1999年提出[7],并在2004年完善總結(jié)。SIFT算法是一種基于計(jì)算機(jī)視覺的算法,主要用來對影像中的局部性特征進(jìn)行偵測與描述,并在空間尺度中尋找極值點(diǎn),同時將其位置、尺度和旋轉(zhuǎn)不變量提取出來。
SIFT算法[8]信息量豐富,適用于對大量特征進(jìn)行快速、準(zhǔn)確的匹配;在目標(biāo)物體很少的條件下,該算法也能夠產(chǎn)生大量的SIFT特征向量;不僅如此,SIFT算法經(jīng)過優(yōu)化之后可能達(dá)到實(shí)時運(yùn)算的效果,且能夠與多種其它形式的特征向量聯(lián)合使用,其拓展性也十分良好。
從理論上說,SIFT特征是一種相似不變量,即對圖像尺度變化和旋轉(zhuǎn)是不變量。SIFT特征對圖像的視點(diǎn)變化、仿射變換、噪聲和光照變化具有了較強(qiáng)的適應(yīng)性,定位精度比較高。
針對空間目標(biāo)紅外特征提取識別問題,以國內(nèi)某衛(wèi)星(以下簡稱國內(nèi)衛(wèi)星)和國外某衛(wèi)星(以下簡稱國外衛(wèi)星)為研究對象,并對其進(jìn)行特征提取和識別。初始考慮如下3種特征提取方案:
1)提取紅外圖像的Harris角點(diǎn),對每個角點(diǎn)生成描述符,然后利用詞袋(BOW)模型進(jìn)行編碼形成每幅圖像的描述子;
2)提取紅外圖像的LBP直方圖特征;
3)提取紅外圖像的灰度直方圖特征。
經(jīng)過實(shí)驗(yàn)采用后兩種特征。對于方案1),實(shí)驗(yàn)發(fā)現(xiàn)在圖像分辨率較高的情況下,確實(shí)能夠提取到很多角點(diǎn),但是當(dāng)圖像分辨率降低到[20 20]以下之后,角點(diǎn)檢測算法基本失效,在圖像上檢測到的角點(diǎn)個數(shù)非常少、甚至檢測不到,因此后續(xù)無法利用BOW模型對圖像進(jìn)行編碼,圖1給出了一組角點(diǎn)檢測失敗的例子。
對于方案2),鑒于最小的圖像分辨率為[10 10],在提取LBP直方圖特征時設(shè)定cell大小為[5 5],并最終對整幅圖像所有cell上的LBP直方圖取平均值,然后做L1歸一化,作為該圖像的LBP直方圖特征。圖2和圖3給出了國內(nèi)衛(wèi)星和國外衛(wèi)星部分圖像及其LBP直方圖特征。從圖中可以看出,衛(wèi)星在不同的姿態(tài)下其LBP直方圖特征也不相同;而且LBP特征直方圖的bin數(shù)量較多,圖像局部特征差異提取較為充分,更有利于微小差異的檢測。
對于方案3),直接提取灰度直方圖特征,直方圖的bin數(shù)設(shè)為16,然后做L1歸一化。圖4和圖5給出了國內(nèi)衛(wèi)星和國外衛(wèi)星部分紅外圖像及其灰度直方圖特征。相比于LBP直方圖特征,灰度直方圖bin的數(shù)量大大減少,但是直方圖特征之間的差異也更為明顯,對于目標(biāo)整體之間的差異檢測非常有效。
圖1 在分辨率為[50 50]下角點(diǎn)檢測失敗的圖像對
圖2 國外衛(wèi)星紅外圖像及其LBP直方圖特征
Fig.2 Infrared images and LBP histogram characteristics of the foreign satellite
圖3 國內(nèi)衛(wèi)星紅外圖像及其LBP直方圖特征
圖4 國外衛(wèi)星紅外圖像及其灰度直方圖特征
圖5 國內(nèi)衛(wèi)星紅外圖像及其灰度直方圖特征
目標(biāo)識別的基本原理[9]是利用檢測手段獲取被識別對象的各種特征參數(shù),并在此基礎(chǔ)上建立參考樣本,然后將被識別對象與參考樣本進(jìn)行比較辨認(rèn)。因此,參考樣本的準(zhǔn)確性和唯一性對于目標(biāo)的準(zhǔn)確識別具有重大影響。
對于空間目標(biāo)紅外識別而言,難點(diǎn)主要在于觀測獲取的目標(biāo)圖像無約束姿態(tài)變化和分辨率低。同一目標(biāo)在不同姿態(tài)下得到的觀測圖像往往差異很大,而不同衛(wèi)星在某些姿態(tài)下的觀測圖像又很相似,即目標(biāo)紅外識別是典型的類內(nèi)變化大,類間變化小的目標(biāo)識別問題。
利用紅外仿真軟件分別生成國內(nèi)衛(wèi)星和國外衛(wèi)星在不同姿態(tài)、不同分辨率下的圖像作為訓(xùn)練樣本和測試樣本。針對每個目標(biāo),分別生成[10 10], [20 20], [30 30], …, [100 100]等10種分辨率下的全姿態(tài)空間樣本圖像,在特定的分辨率下,分別在,,三個軸上每隔45°采樣一個樣本,總共生成8×8×8=512幅圖像,并針對每種分辨率按照統(tǒng)一的方法分別訓(xùn)練分類器,從而給出分辨率對識別結(jié)果的影響。圖6為兩個目標(biāo)在不同像素大小下的樣本示例。
圖6 不同分辨率下國外衛(wèi)星和國內(nèi)衛(wèi)星樣本示例圖
目前常用的分類器有歐氏距離最近鄰法和最近中心鄰法,方法簡單直觀且計(jì)算量小,但是對于復(fù)雜目標(biāo)的分類能力有限[10];神經(jīng)網(wǎng)絡(luò)分類器具有自適應(yīng)性、自組織和自學(xué)習(xí)等優(yōu)點(diǎn),但是其收斂速度慢,且對樣本依賴性較強(qiáng);而SVM分類器在解決小樣本和高維模式識別問題上表現(xiàn)優(yōu)異[11],且分類速度較快。對于空間目標(biāo)這種難以長時間持續(xù)跟蹤獲得大量樣本的識別問題而言,SVM分類器將具有更好的適應(yīng)性。
SVM分類器通過訓(xùn)練學(xué)習(xí)得到樣例的輸入輸出對應(yīng)關(guān)系[12],然后根據(jù)新輸入樣本預(yù)測輸出結(jié)果。其基本思想可以概括為[13]通過內(nèi)積函數(shù)定義的非線性變化將輸入空間變換到一個高維空間,然后在這個新空間中求解最優(yōu)分類面。
1)SVM基本原理
SVM是一個兩類分類器,對于給定標(biāo)記好的訓(xùn)練樣本,SVM輸出一個最優(yōu)化的分隔超平面,最優(yōu)化度量可以定義為超平面距離所有訓(xùn)練樣本的最小距離,或稱為間隔(margin),最優(yōu)分割超平面能夠最大化訓(xùn)練數(shù)據(jù)的間隔。
下式定義經(jīng)典超平面表達(dá)式:
|+T|=1 (2)
式中:稱為權(quán)重向量;稱為偏置;式中表示離超平面最近的點(diǎn)。點(diǎn)到超平面(,)的距離為:
對于經(jīng)典超平面,表達(dá)式中的分子為1,因此支持向量到經(jīng)典超平面的距離為:
間隔(margin)表示為,取值是最近距離的兩倍:
最大化轉(zhuǎn)化為附加限制條件下最小化函數(shù),限制條件為超平面將所有訓(xùn)練樣本x正確分類:
即
式中:y表示樣本的類別標(biāo)記。上式可以通過拉格朗日乘數(shù)法得到最優(yōu)超平面的權(quán)重向量和偏置。
2)軟間隔SVM分類器
上述SVM學(xué)習(xí)方法是針對線性可分問題的,并不適用于線性不可分訓(xùn)練數(shù)據(jù),因?yàn)榇藭r上述不等式約束并不能都成立。因此,需要修改間隔最大化,使其成為軟間隔最大化。假設(shè)給定一個特征空間上的訓(xùn)練數(shù)據(jù)集:
={(1,1), (2,2), …, (x,y)} (8)
式中:x為第個特征向量,y為x的類標(biāo)記。由于訓(xùn)練數(shù)據(jù)集是線性不可分的,說明訓(xùn)練數(shù)據(jù)中存在一些特異點(diǎn),將這些特異點(diǎn)除去后,其余樣本組成的集合將仍然是線性可分的。為此可以對每個樣本點(diǎn)(x,y)引進(jìn)一個松弛變量≥0,使函數(shù)間隔加上松弛變量大于等于1。約束條件變?yōu)椋?/p>
即
y(Tx+)≥1-,≥0,=1, 2, …,(12)
3)核函數(shù)與SVM分類器
針對線性問題,線性SVM是有效辦法;而針對非線性問題,則需要將原空間的數(shù)據(jù)映射到新空間,在新空間里使用線性分類學(xué)習(xí)方法。此時,涉及度量樣本間相似性的徑向基函數(shù)(RBF):
式中:、為兩個樣本;為經(jīng)驗(yàn)系數(shù)。
識別過程如下:
1)在每種分辨率條件下,國內(nèi)衛(wèi)星和國外衛(wèi)星的紅外圖像被隨機(jī)劃分為兩份,其中512幅圖像(國內(nèi)衛(wèi)星和國外衛(wèi)星各256幅)作為訓(xùn)練樣本,剩余512幅作為測試樣本。
2)對訓(xùn)練樣本分別進(jìn)行LBP直方圖特征和灰度直方圖特征的提取。
3)利用提取的兩類特征對線性核SVM分類器和RBF核SVM分類器進(jìn)行訓(xùn)練,分別確定其最優(yōu)超平面的權(quán)重向量和偏置。
4)將訓(xùn)練好的兩種SVM分類器對測試樣本進(jìn)行分類識別,并統(tǒng)計(jì)兩種SVM分類器分別用LBP直方圖特征和灰度直方圖特征進(jìn)行目標(biāo)識別的準(zhǔn)確率、目標(biāo)樣本的召回率以及虛警率。
假設(shè)測試樣本中共有個國外衛(wèi)星樣本,個國內(nèi)衛(wèi)星樣本,經(jīng)過分類之后,共有P個國外衛(wèi)星樣本被正確識別,N的國外衛(wèi)星樣本被錯誤識別,有:
=P+N(14)
共有N個國內(nèi)衛(wèi)星樣本被正確識別,P個國內(nèi)衛(wèi)星樣本被錯誤識別,有:
=N+P(15)
識別準(zhǔn)確率定義為:
Acc=(P+N)/(+) (16)
召回率定義為:
Rec=P/(P+N) (17)
虛警率定義為:
FalA=P/(P+P) (18)
5)將上述過程通過隨機(jī)劃分樣本的方式對每種分辨率下的紅外圖像進(jìn)行100次實(shí)驗(yàn),統(tǒng)計(jì)不同核函數(shù)的SVM分類器和不同特征的識別結(jié)果。
圖7為使用線性核SVM和LBP直方圖特征,針對不同分辨率圖像,在訓(xùn)練集和測試集上的識別準(zhǔn)確率、召回率和虛警率曲線,表1給出相應(yīng)的性能指標(biāo)數(shù)據(jù)。
圖7 線性核SVM+LBP直方圖特征的識別結(jié)果曲線
表1 線性核SVM+LBP直方圖特征的識別結(jié)果
圖8為使用RBF核SVM和LBP直方圖特征,針對不同分辨率圖像,在訓(xùn)練集和測試集上的識別準(zhǔn)確率、召回率和虛警率曲線,表2給出相應(yīng)的性能指標(biāo)數(shù)據(jù)。
圖8 RBF核SVM+LBP直方圖特征的識別結(jié)果曲線
表2 RBF核SVM+LBP直方圖特征的識別結(jié)果
圖9為使用線性核SVM和灰度直方圖特征,針對不同分辨率圖像,在訓(xùn)練集和測試集上的識別準(zhǔn)確率、召回率和虛警率曲線,表3給出相應(yīng)的性能指標(biāo)數(shù)據(jù)。
圖10為使用RBF核SVM和灰度直方圖特征,針對不同分辨率圖像,在訓(xùn)練集和測試集上的識別準(zhǔn)確率、召回率和虛警率曲線,表4給出相應(yīng)的性能指標(biāo)數(shù)據(jù)。
圖9 線性核SVM+灰度直方圖特征的識別結(jié)果曲線
識別過程中,不同方案的時間復(fù)雜度其區(qū)別主要體現(xiàn)在兩個方面,一是提取特征的復(fù)雜度,二是SVM分類器核函數(shù)的復(fù)雜度。本次實(shí)驗(yàn)選取了目標(biāo)紅外圖像的LBP直方圖特征和灰度直方圖特征,由于LBP直方圖特征需要先提取目標(biāo)的LBP特征,然后提取其直方圖特征,對于×大小的圖像,其時間復(fù)雜度為(92),而灰度直方圖特征的時間復(fù)雜度為(2),因此相比于灰度直方圖特征,LBP直方圖特征的時間復(fù)雜度更高;其次,在核函數(shù)的選擇上,由于RBF核函數(shù)涉及到非線性運(yùn)算,因此其時間復(fù)雜度也要高于線性核函數(shù)。
表3 線性核SVM+灰度直方圖特征的識別結(jié)果
圖10 RBF核SVM+灰度直方圖特征的識別結(jié)果曲線
表4 RBF核SVM+灰度直方圖特征的識別結(jié)果
通過上述實(shí)驗(yàn)結(jié)果可知:
1)LBP直方圖特征和灰度直方圖特征均具有目標(biāo)識別的能力。從圖7和圖9中可以看出,在使用線性核SVM分類器的條件下,LBP直方圖特征在紅外圖像分辨率達(dá)到[10 10]的時候,識別率已經(jīng)大于85%;灰度直方圖特征在分辨率為[30 30]的時候也已經(jīng)大于85%的識別率;而從圖8和圖10中可以看出,在使用RBF核SVM分類器的條件下,在紅外圖像分辨率為[10 10]的時候,LBP直方圖特征和灰度直方圖特征均達(dá)到90%以上;
2)無論是LBP直方圖特征還是灰度直方圖特征,使用RBF核SVM分類器的識別效能要優(yōu)于線性核SVM分類器。空間目標(biāo)識別并非線性可分問題,通過RBF核函數(shù)將非線性問題映射到線性空間,再進(jìn)行分類識別,結(jié)果大大提高了識別概率,也說明了SVM分類器處理線性不可分問題時的良好性能,當(dāng)然其時間復(fù)雜度也更高;此外,隨著紅外圖像分辨率的增加,識別率整體呈上升趨勢。這是由于隨著分辨率的增加,紅外圖像包含的信息越多,與直觀預(yù)期相符合;
3)在使用線性核SVM分類器的時候,分辨率較低時,LBP直方圖特征的識別率要優(yōu)于灰度直方圖特征;而分辨率較高的時候,灰度直方圖特征優(yōu)于LBP直方圖特征。
[1] 王藝. 基于激光技術(shù)空間碎片的主動移除[D]. 南昌: 南昌航空大學(xué), 2015.
WANG Yi. Active Space Debris Removal Based on Laser Technology [D]. Nanchang: Nanchang Hangkong University, 2015.
[2] 蘇剛. 云背景空中目標(biāo)紅外仿真研究[D]. 武漢: 華中科技大學(xué), 2014.
SU Gang. The Study on the Infrared Simulation of Aerial Target under Cloud Background [D]. Wuhan: Huazhong University of Science and Technology, 2014.
[3] 季圓圓. 典型地物紅外特征仿真關(guān)鍵技術(shù)研究[D]. 大連: 大連海事大學(xué), 2014.
JI YuanYuan. Study on Key Techniques of Simulation for Infrared Characteristics of Typical Ground Objects [D]. Dalian: Dalian Maritime University, 2014.
[4] 馬偉. 衛(wèi)星熱輻射特性及其空間輻照環(huán)境效應(yīng)研究[D]. 南京: 南京理工大學(xué), 2011.
MA Wei. Research on Thermal Radiation Characteristics and Effects of Space Exposure [D]. Nanjing: Nanjing University of Science and Technology, 2011.
[5] 曹凱. 基于特征描述子的人臉識別方法研究[D]. 西安: 西安電子科技大學(xué), 2014.
CAO Kai. A Study of Feature Descriptors For Face Recognition [D]. Xian: Xidian University, 2014.
[6] Y. Rodriguez, S. Marcel. Face Authentication Using Adapted Local Binary Pattern Histograms[C]//,, 2006, 3954(4): 321-332.
[7] D. G. Lowe. Distinctive image features from scale invariant key points[J]., 2004, 60(2): 91-110.
[8] 劉朵. 基于特征描述子的目標(biāo)跟蹤研究[D]. 西安: 西安電子科技大學(xué), 2013.
LIU Duo. Research on Feature Descriptor-based Object Tracking [D]. Xi’an: Xidian University, 2013.
[9] 王彥. 基于圖像的空間目標(biāo)識別關(guān)鍵技術(shù)研究[D]. 西安: 西安電子科技大學(xué), 2010.
WANG Yan. Key Technology of Space Object Identification Based on Image[D]. Xi’an: Xidian University, 2010.
[10] 譚敏潔. 基于壓縮感知的雷達(dá)一維距離像目標(biāo)識別[D]. 成都: 電子科技大學(xué), 2015.
TAN MinJie. High Range Resolution Profile Target Recognition Based on Compressed Sensing[D]. Chengdu: University of Electronic Science and Technology of China, 2015.
[11] 李瑞帥. 基于Adaboost和SVM的快速行人檢測方法研究[D]. 桂林: 桂林電子科技大學(xué), 2014.
LI RuiShuai. Research on Fast Pedestrian Detection Method Based on Adaboost and SVM[D]. Guilin: Guilin University of Electronic Technology, 2014.
[12] 壽娜. 基于圖像的行人檢測系統(tǒng)研究[D]. 浙江: 浙江工業(yè)大學(xué), 2013.
SHOU Na. The Research of Pedestrian Detection System For Images [D]. Zhejiang: Zhejiang University of Technology, 2013.
[13] 張迪飛, 張金鎖, 姚克明,等. 基于SVM分類的紅外艦船目標(biāo)識別[J].紅外與激光工程, 2016, 45(1): 167-172.
ZHANG DiFei, ZHANG JinSuo, YAO KeMing, et al. Infrared Ship-Target Recognition Based on SVM Classification[J]., 2016, 45(1): 167-172.
Infrared Feature Extraction and Recognition Technology of Space Target
LI Ruidong,SUN Xiechang,LI Meng
(,100086,)
Dissimilarity is usually large among the same space target but small among different ones for infrared recognition. Aiming at these problems, a method of space-target infrared feature extraction and recognition with the use of SVM classifier combining with LBP histogram feature and gray histogram feature is proposed. A domestic satellite and a foreign satellite are taken as the objects of study. LBP histogram feature and gray histogram feature of their infrared images are extracted. Sample images under different attitudes and resolutions are generated by the simulation software, and they are divided into two parts, one for SVM classifier training set and the other for test set. The test results indicate that not only can both LBP histogram feature and gray histogram feature recognize the space targets with a high accuracy, but also their recognition performances relate to the resolution of infrared images and the kernel functions of SVM.
infrared recognition,feature extraction,space target
TN219
A
1001-8891(2017)05-0427-09
2017-01-04;
2017-04-18。
李瑞東(1992-),男,湖北武穴人,碩士研究生,主要從事空間遙感信息處理方向的工作。
國家863計(jì)劃課題(2014AA7032031E);火箭軍裝備預(yù)研項(xiàng)目。