專利名稱：一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，屬于遙感數(shù)據(jù)處理方法與應(yīng)用技術(shù)領(lǐng)域，適用于高光譜傳感器輻射定標(biāo)的理論方法和應(yīng)用技術(shù)研究。
背景技術(shù)：
遙感器輻射定標(biāo)是遙感數(shù)據(jù)定量化應(yīng)用的前提，該類技術(shù)的發(fā)展能夠大大推動(dòng)遙感數(shù)據(jù)的應(yīng)用深度和拓展遙感數(shù)據(jù)的應(yīng)用廣度。輻射定標(biāo)作為一種衛(wèi)星遙感數(shù)據(jù)預(yù)處理技術(shù)的一個(gè)重要環(huán)節(jié)，其目的在于建立傳感器遙感器入瞳處輻亮度與其輸出的數(shù)字計(jì)數(shù)值 (DN值)之間的關(guān)系，以確定遙感器的輸出信號(hào)對(duì)入瞳輻亮度的響應(yīng)。在輻射定標(biāo)中，設(shè)DN 是遙感影像的數(shù)字計(jì)數(shù)值，L是待定標(biāo)傳感器的入瞳處輻亮度。則輻射定標(biāo)即是通過(guò)如下關(guān)系，確定傳感器探測(cè)單元的定標(biāo)系數(shù)G和B :L = G*DN+B其中，G和B分別是定標(biāo)系數(shù)中的增益(Gain)和偏移(Bias)。輻射定標(biāo)方法包括衛(wèi)星發(fā)射前實(shí)驗(yàn)室定標(biāo)、在軌星上定標(biāo)、場(chǎng)地定標(biāo)及交叉定標(biāo)等。不同的定標(biāo)方法分別在衛(wèi)星的不同階段發(fā)揮著一定的作用，但每一種定標(biāo)方法都有其特點(diǎn)與局限適用性。衛(wèi)星發(fā)射前的實(shí)驗(yàn)室定標(biāo)，關(guān)鍵是評(píng)估工程設(shè)計(jì)指標(biāo)和建立滿足不同應(yīng)用要求的定標(biāo)系統(tǒng)；在軌星上定標(biāo)基于星載定標(biāo)器的飛行中定標(biāo)采用星上人工輻射源定標(biāo)、天體或外太空，基于太陽(yáng)和漫反射板定標(biāo)，利用星上黑體和冷空間作為輻射參考輻射源定標(biāo)等手段進(jìn)行飛行中定標(biāo)；對(duì)于沒(méi)有星上定標(biāo)系統(tǒng)的衛(wèi)星，為了動(dòng)態(tài)跟蹤傳感器輻射性能的變化，美國(guó)、法國(guó)在這一方面遙遙領(lǐng)先。在軌運(yùn)行期間可以采用基于陸地(或海面)特性的“替代定標(biāo)”，即場(chǎng)地定標(biāo)。利用一個(gè)經(jīng)過(guò)嚴(yán)格鑒定的地表面試驗(yàn)場(chǎng)輻射校正場(chǎng)的地面同步實(shí)驗(yàn)和大氣輻射傳輸模式模型，可以實(shí)現(xiàn)確定傳感器入瞳輻亮度的輻射定標(biāo)。但由于場(chǎng)地定標(biāo)需要耗費(fèi)大量的人力、物力與財(cái)力，故通常定標(biāo)間隔周期較長(zhǎng)。因此需要一種更加經(jīng)濟(jì)與便捷的方法。利用高輻射精度的衛(wèi)星傳感器對(duì)精度相對(duì)較低的衛(wèi)星傳感器進(jìn)行交叉輻射定標(biāo)方法(簡(jiǎn)稱“交叉定標(biāo)”)應(yīng)運(yùn)而生。交叉輻射定標(biāo)是新近發(fā)展起來(lái)的一種無(wú)場(chǎng)地定標(biāo)技術(shù)，得到了國(guó)際遙感界的廣泛重視，被國(guó)際地球觀測(cè)系統(tǒng)委員會(huì)(CE0Q定標(biāo)檢驗(yàn)組 (WGCV)所推薦。交叉定標(biāo)是利用一臺(tái)標(biāo)定好的傳感器與待標(biāo)定傳感器在同一時(shí)間和相似的幾何條件進(jìn)行觀測(cè)，然后進(jìn)行相對(duì)比較，實(shí)現(xiàn)定標(biāo)。交叉輻射定標(biāo)方法無(wú)需建立地面校正場(chǎng)，該方法最大的優(yōu)點(diǎn)是定標(biāo)成本低，可以實(shí)現(xiàn)高頻次、多傳感器間的輻射定標(biāo)，并且可以實(shí)現(xiàn)對(duì)歷史數(shù)據(jù)的輻射定標(biāo)。但是該方法， 對(duì)參考傳感器和被定標(biāo)傳感器之間的通道設(shè)置、通道光譜響應(yīng)函數(shù)、空間分辨率、過(guò)境時(shí)間、回訪周期、幾何配準(zhǔn)精度等方面有嚴(yán)格要求，并且定標(biāo)精度極大的依賴于參考傳感器自身的絕對(duì)輻射定標(biāo)精度。對(duì)于缺乏星上定標(biāo)系統(tǒng)的衛(wèi)星傳感器，只能借助場(chǎng)地定標(biāo)或交叉定標(biāo)方法。特別是對(duì)于技術(shù)水平相對(duì)落后的發(fā)展中國(guó)家的衛(wèi)星遙感具有重要的作用。國(guó)內(nèi)外關(guān)于星-星交叉定標(biāo)的開(kāi)展了較多的研究，并取得了較好的效果。但是現(xiàn)有的研究大都是基于多光譜傳感器之間的輻射定標(biāo)，而目前在軌的星載高光譜傳感器極少，僅有美國(guó)的EOl/Hyperion、歐空局PR0BA/CHRIS與中國(guó)的HJ/HIS等少數(shù)幾個(gè)。受高光譜數(shù)據(jù)的可獲得性、訪問(wèn)時(shí)間、空間覆蓋范圍、天氣條件、空間分辨率的差異等諸多因素影響，通常很難找到匹配的兩套高光譜影像。因此對(duì)于高光譜的在軌輻射定標(biāo)需要尋找新的技術(shù)方法。本發(fā)明就是為解決高光譜在軌輻射定標(biāo)的一種新型技術(shù)。

發(fā)明內(nèi)容
本發(fā)明的目的是提出一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，它克服了傳統(tǒng)高光譜數(shù)據(jù)由于不能找到適合的相匹配的參考高光譜影像數(shù)據(jù)而不能進(jìn)行交叉定標(biāo)的局限，有效地提高了高光譜遙感數(shù)據(jù)的可用性，是一種穩(wěn)定性好、可靠性高、精確度高的高光譜傳感器交叉輻射定標(biāo)方法。本發(fā)明的技術(shù)解決方案是一種在不能獲得輻射定標(biāo)精度較高的高光譜衛(wèi)星影像的情況下，利用輻射定標(biāo)精度較高的多光譜影像實(shí)現(xiàn)對(duì)高光譜數(shù)據(jù)輻射定標(biāo)的方法。該方法主要是基于優(yōu)化算法重構(gòu)高光譜傳感器入瞳處輻亮度曲線實(shí)現(xiàn)輻射定標(biāo)，利用6S大氣輻射傳輸模型對(duì)兩傳感器進(jìn)行入瞳處輻亮度和氣體特征吸收模擬，并通過(guò)計(jì)算高光譜有效波段的光譜匹配因子，利用光譜匹配因子實(shí)現(xiàn)有效波段的輻射定標(biāo)，最后利用優(yōu)化算法由部分有效波段的輻亮度值獲得高光譜傳感器的入瞳處輻亮度曲線，從而完成輻射定標(biāo)。本發(fā)明一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，如

圖1所示，其步驟如下(1)、選擇無(wú)云的高光譜影像數(shù)據(jù)，影像中包含幾種大塊均勻地物，并對(duì)影像進(jìn)行必要的預(yù)處理(如壞波段剔除)后，作為待定標(biāo)的影像；O)、選擇與高光譜傳感器成像時(shí)間接近且覆蓋同一區(qū)域的輻射定標(biāo)精度高的多光譜遙感影像作為參考影像；(3)、裁剪出兩幅影像即待定標(biāo)影像和參考影像的共同覆蓋區(qū)域；選擇同名特征點(diǎn)，對(duì)兩幅影像進(jìn)行幾何精校正處理；在影像上選擇均勻地物作為交叉定標(biāo)的感興趣區(qū) (ROI)；(4)、根據(jù)兩幅影像的成像條件，利用6S大氣輻射傳輸模型計(jì)算兩傳感器各波段的入瞳輻亮度和氣體吸收特征曲線，并分別對(duì)兩傳感器的入瞳輻亮度光譜曲線去除氣體吸收特征；(5)、根據(jù)步驟(4)獲得的去除氣體吸收特征后的光譜曲線，使用優(yōu)化算法，以重構(gòu)曲線與模擬曲線最接近為判斷依據(jù)，篩選出最優(yōu)波段與最優(yōu)插值方案；(6)、根據(jù)中心波長(zhǎng)最接近的規(guī)則用多光譜傳感器作為參考，對(duì)高光譜傳感器最優(yōu)波段進(jìn)行交叉輻射定標(biāo)；(7)、采用步驟( 得到的最優(yōu)插值方案，對(duì)定標(biāo)后的高光譜波段輻亮度進(jìn)行插值，重構(gòu)高光譜曲線，然后對(duì)重構(gòu)曲線加載氣體吸收特征，作為參考傳感器的輻亮度；(8)、在待定標(biāo)的高光譜影像的感興趣區(qū)內(nèi)，選擇若干個(gè)像元各波段DN值和步驟 (7)求得的定標(biāo)參考傳感器輻亮度曲線進(jìn)行擬合，得到高光譜傳感器各波段定標(biāo)系數(shù)(增益和偏移量)，從而實(shí)現(xiàn)多光譜傳感器對(duì)高光譜傳感器的交叉輻射定標(biāo)。其中，步驟(1)中所述的高光譜影像數(shù)據(jù)選擇，要求盡可能選擇大氣條件與地面光譜特性穩(wěn)定的輻射定標(biāo)場(chǎng)，以減少兩幅影像因成像時(shí)間差異以及傳感器光譜響應(yīng)差異所帶來(lái)的誤差，有利于衛(wèi)星傳感器交叉輻射定標(biāo)。其中，步驟O)中所述的多光譜遙感影像的選擇，要求盡可能選擇與高光譜成像時(shí)間一致，以減少因成像時(shí)間不同，導(dǎo)致大氣條件的差異所帶來(lái)的誤差。其中，步驟(3)中所述的幾何精校正，分別在參考影像和待校正的影像上選擇幾何特征明顯的地物作為控制點(diǎn)，對(duì)控制點(diǎn)的空間分布應(yīng)該盡量均勻覆蓋整幅影像，控制點(diǎn)的數(shù)據(jù)量通常應(yīng)該大于20個(gè)，幾何校正的整體誤差應(yīng)該控制在半個(gè)像素以內(nèi)。在待定標(biāo)的高光譜影像和已定標(biāo)的參考多光譜遙感影像上，盡量在星下點(diǎn)附近選擇均勻地物作為感興趣區(qū)(ROI)，減少地表因非朗伯特性和大視場(chǎng)觀測(cè)所帶來(lái)的不確定性；另外ROI盡量覆蓋遙感器的整個(gè)動(dòng)態(tài)響應(yīng)范圍；其中，步驟中所述的“利用6S大氣輻射傳輸模型計(jì)算兩傳感器各波段的入瞳輻亮度和氣體吸收特征曲線，并分別對(duì)其去除氣體吸收特征”，為了在后續(xù)的優(yōu)化插值環(huán)節(jié)上盡可能減小誤差，并考慮到輻亮度光譜曲線上所有氣體特征吸收谷是由大氣分子的物理特性決定的，僅依靠單純的數(shù)學(xué)插值算法不能將其形狀復(fù)原，本發(fā)明采用大氣輻射傳輸模型中計(jì)算得出的氣體吸收函數(shù)f(0 < fi < 1)對(duì)輻亮度曲線L進(jìn)行變換，公式如下L' = L/f上式中L是大氣輻射傳輸模型模擬得到的某波段入瞳處輻亮度，f是大氣輻射傳輸模型模擬得到的對(duì)應(yīng)波段的氣體吸收函數(shù)，L'是某波段去除氣體吸收特征后的輻亮度。其中，其步驟(5)中所述的幾“根據(jù)步驟(4)獲得的去氣體吸收特征后的光譜曲線，使用優(yōu)化算法，以重構(gòu)曲線與模擬曲線最接近為判斷依據(jù)，篩選出最優(yōu)波段與最優(yōu)插值方案”，其含義說(shuō)明如下使用優(yōu)化算法，從多光譜傳感器波段可完全覆蓋的高光譜傳感器對(duì)應(yīng)波段(即有效波段)中選擇若干波段，通過(guò)設(shè)置不同的最優(yōu)波段數(shù)量和不同插值算法， 最后獲得與去除氣體吸收特征曲線誤差最小的插值曲線，被選出的對(duì)應(yīng)波段為最優(yōu)波段， 對(duì)應(yīng)的插值方法為最優(yōu)插值方法。其中，步驟(6)中所述的“根據(jù)中心波長(zhǎng)最接近的規(guī)則用多光譜傳感器作為參考對(duì)高光譜傳感器最優(yōu)波段進(jìn)行交叉輻射定標(biāo)”，其含義說(shuō)明如下兩傳感器對(duì)應(yīng)通道的表觀輻亮度在6S大氣輻射傳輸模型中存在著一定的差異和特定的關(guān)系，這種關(guān)系在對(duì)應(yīng)的遙感影像中應(yīng)同樣成立。故Lh6s = K^Ltm6sLh map = K*LTM mapLh 6s和Ltm 6s分別表示利用6S大氣輻射傳輸模型模擬的高光譜和多光譜傳感器入瞳輻亮度，同理，Lh map和LTM—map分別表示影像DN值對(duì)應(yīng)的高光譜和多光譜傳感器入瞳輻亮度L，其中K表示光譜匹配因子。其中，步驟⑶中所述的“針對(duì)不同波段，將所有ROI在待定標(biāo)的高光譜影像的DN 值和步驟(7)求得的定標(biāo)輻亮度曲線進(jìn)行線性擬合，可求得高光譜傳感器各波段的定標(biāo)系數(shù)(增益和偏移量)，從而實(shí)現(xiàn)多光譜傳感器對(duì)高光譜傳感器的交叉輻射定標(biāo)?！?，為了有效降低定標(biāo)系數(shù)的不確定性，本發(fā)明采用基于區(qū)域平均值的定標(biāo)系數(shù)生成法，本發(fā)明采取與原定標(biāo)系數(shù)相同的處理流程獲取定標(biāo)系數(shù)如下根據(jù)灰度值的不同，將感興趣區(qū)進(jìn)行分區(qū)， 各波段分別計(jì)算各分區(qū)的平均灰度值M,基于M和步驟(7)求得的各分區(qū)的定標(biāo)輻亮度L進(jìn)行線性擬合，可得第i波段定標(biāo)系數(shù)Ai
權(quán)利要求
1.一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于包含以下步驟(1)、選擇無(wú)云的高光譜影像數(shù)據(jù)，影像中包含幾種大塊均勻地物，并對(duì)影像進(jìn)行預(yù)處理后，作為待定標(biāo)的影像；O)、選擇與高光譜傳感器成像時(shí)間接近且覆蓋同一區(qū)域的輻射定標(biāo)精度高的多光譜遙感影像作為參考影像；(3)、裁剪出兩幅影像的共同覆蓋區(qū)域；選擇同名特征點(diǎn)，對(duì)兩幅影像進(jìn)行幾何精校正處理；在影像上選擇均勻地物作為交叉定標(biāo)的感興趣區(qū)，即ROI ；(4)、根據(jù)兩幅影像的成像條件，利用6S大氣輻射傳輸模型計(jì)算兩傳感器各波段的入瞳輻亮度和氣體吸收特征曲線，并分別對(duì)兩傳感器的入瞳輻亮度光譜曲線去除氣體吸收特征；(5)、根據(jù)步驟(4)獲得的去除氣體吸收特征后的光譜曲線，使用優(yōu)化算法，以重構(gòu)曲線與模擬曲線最接近為判斷依據(jù)，篩選出最優(yōu)波段與最優(yōu)插值方案；(6)、根據(jù)中心波長(zhǎng)最接近的規(guī)則用多光譜傳感器作為參考對(duì)高光譜傳感器最優(yōu)波段進(jìn)行交叉輻射定標(biāo)；(7)、采用步驟( 得到的最優(yōu)插值方案，對(duì)定標(biāo)后的高光譜波段輻亮度進(jìn)行插值，重構(gòu)高光譜曲線，然后對(duì)重構(gòu)曲線加載氣體吸收特征，作為參考傳感器的輻亮度；(8)、在待定標(biāo)的高光譜影像的感興趣區(qū)內(nèi)，選擇若干個(gè)像元各波段DN值和步驟(7)求得的定標(biāo)參考傳感器輻亮度曲線進(jìn)行擬合，得到高光譜傳感器各波段定標(biāo)系數(shù)，從而實(shí)現(xiàn)多光譜傳感器對(duì)高光譜傳感器的交叉輻射定標(biāo)。
2.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟(1)中所述的選擇無(wú)云的高光譜影像數(shù)據(jù)，要求選擇大氣條件與地面光譜特性穩(wěn)定的輻射定標(biāo)場(chǎng)，以減少兩幅影像因成像時(shí)間差異以及傳感器光譜響應(yīng)差異所帶來(lái)的誤差，有利于衛(wèi)星傳感器交叉輻射定標(biāo)。
3.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟O)中要求選擇與高光譜成像時(shí)間一致的多光譜遙感影像，以減少因成像時(shí)間不同，導(dǎo)致大氣條件的差異所帶來(lái)的誤差。
4.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟(3)中所述的幾何精校正，分別在參考影像和待校正的影像上選擇幾何特征明顯的地物作為控制點(diǎn)；對(duì)控制點(diǎn)的空間分布應(yīng)該均勻覆蓋整幅影像，控制點(diǎn)的數(shù)據(jù)量應(yīng)該大于20個(gè)，幾何校正的整體誤差應(yīng)該控制在半個(gè)像素以內(nèi)；在待定標(biāo)的高光譜影像和已定標(biāo)的參考多光譜遙感影像上，在星下點(diǎn)附近選擇均勻地物作為R0I，減少地表因非朗伯特性和大視場(chǎng)觀測(cè)所帶來(lái)的不確定性；另外ROI覆蓋遙感器的整個(gè)動(dòng)態(tài)響應(yīng)范圍。
5.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟中所述的“利用6S大氣輻射傳輸模型計(jì)算兩傳感器各波段的入瞳輻亮度和氣體吸收特征曲線，并分別對(duì)其去除氣體吸收特征”，具體如下為了在后續(xù)的優(yōu)化插值環(huán)節(jié)上盡可能減小誤差，并考慮到輻亮度光譜曲線上所有氣體特征吸收谷是由大氣分子的物理特性決定的，僅依靠單純的數(shù)學(xué)插值算法不能將其復(fù)原的形狀，故采用大氣輻射傳輸模型中計(jì)算得出的氣體吸收函數(shù)f(0 < < 1)對(duì)輻亮度曲線L進(jìn)行變換，公式如下L' = L/f上式中L是大氣輻射傳輸模型模擬得到的某波段入瞳處輻亮度，f是大氣輻射傳輸模型模擬得到的對(duì)應(yīng)波段的氣體吸收函數(shù)，L'是某波段去除氣體吸收特征后的輻亮度。
6.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于其步驟(5)中所述的幾“根據(jù)步驟(4)獲得的去氣體吸收特征后的光譜曲線，使用優(yōu)化算法，以重構(gòu)曲線與模擬曲線最接近為判斷依據(jù)，篩選出最優(yōu)波段與最優(yōu)插值方案”， 具體如下使用優(yōu)化算法，從多光譜傳感器波段可完全覆蓋的高光譜傳感器對(duì)應(yīng)波段中選擇若干波段，通過(guò)設(shè)置不同的最優(yōu)波段數(shù)量和不同插值算法，最后獲得與去除氣體吸收特征曲線誤差最小的插值曲線，被選出的對(duì)應(yīng)波段為最優(yōu)波段，對(duì)應(yīng)的插值方法為最優(yōu)插值方法。
7.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟(6)中所述的“根據(jù)中心波長(zhǎng)最接近的規(guī)則用多光譜傳感器作為參考對(duì)高光譜傳感器最優(yōu)波段進(jìn)行交叉輻射定標(biāo)”，具體如下兩傳感器對(duì)應(yīng)通道的表觀輻亮度在6S 大氣輻射傳輸模型中存在著一定的差異和特定的關(guān)系，這種關(guān)系在對(duì)應(yīng)的遙感影像中應(yīng)同樣成立；故
8.根據(jù)權(quán)利要求1所述的一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，其特征在于步驟(8)中所述的“針對(duì)不同波段，將所有ROI在待定標(biāo)的高光譜影像的DN值和步驟(7)求得的定標(biāo)輻亮度曲線進(jìn)行線性擬合，可求得高光譜傳感器各波段的定標(biāo)系數(shù)， 從而實(shí)現(xiàn)多光譜傳感器對(duì)高光譜傳感器的交叉輻射定標(biāo)”，為了有效降低定標(biāo)系數(shù)的不確定性，采用基于區(qū)域平均值的定標(biāo)系數(shù)生成法，采取與原定標(biāo)系數(shù)相同的處理流程獲取定標(biāo)系數(shù)如下根據(jù)灰度值的不同，將感興趣區(qū)進(jìn)行分區(qū)，各波段分別計(jì)算各分區(qū)的平均灰度值■，基于M和步驟(7)求得的各分區(qū)的定標(biāo)輻亮度L進(jìn)行線性擬合，可得第i波段定標(biāo)系數(shù)Ai:Aj=DN/L +Br,其中，Bi表示偏移量。
全文摘要
本發(fā)明一種基于優(yōu)化算法的高光譜傳感器交叉輻射定標(biāo)方法，(1)、選擇待定標(biāo)的影像；(2)、選擇參考影像；(3)、裁剪出兩幅影像共同覆蓋區(qū)域，選擇均勻地物作為交叉定標(biāo)的感興趣區(qū)；(4)、計(jì)算兩傳感器各波段的入瞳輻亮度和氣體吸收特征曲線，去除氣體吸收特征；(5)、篩選出最優(yōu)波段與最優(yōu)插值方案；(6)、對(duì)高光譜傳感器最優(yōu)波段進(jìn)行交叉輻射定標(biāo)；(7)、對(duì)定標(biāo)后的高光譜波段輻亮度進(jìn)行插值，重構(gòu)高光譜曲線，加載氣體吸收特征，作為參考傳感器的輻亮度；(8)、實(shí)現(xiàn)多光譜傳感器對(duì)高光譜傳感器的交叉輻射定標(biāo)。本發(fā)明方法能在缺乏參考高光譜影像的條件下，能實(shí)現(xiàn)對(duì)星載成像光譜儀的在軌交叉輻射定標(biāo)，具有業(yè)務(wù)化推廣前景。
文檔編號(hào)G01S7/497GK102540166SQ201110398868
公開(kāi)日2012年7月4日 申請(qǐng)日期2011年12月5日 優(yōu)先權(quán)日2011年12月5日
發(fā)明者周冠華, 姜禾, 徐武健, 賈斐, 趙慧潔 申請(qǐng)人:北京航空航天大學(xué)