專利名稱：：基于隱馬爾科夫樹模型的sar圖像去噪方法
技術(shù)領(lǐng)域：
：本發(fā)明屬于圖像處理
技術(shù)領(lǐng)域：
，特別是一種涉及SAR圖像去噪的.方法，可用于對'SAR圖像、自然圖像及醫(yī)學(xué)圖像的去噪。
背景技術(shù)：
：合成孔徑雷達SAR作為主動雷達，具有不受光照、天氣條件影響的特點，可以全天候、全天時對地觀測，還可以透過地表和植被獲取信息，在農(nóng)業(yè)、林業(yè)、地質(zhì)、環(huán)境、水文、海洋、災(zāi)害、測繪與軍事領(lǐng)域得到了廣泛的應(yīng)用。SAR圖像中由于成像散射體散射回波的相干作用造成的斑點噪聲的存在不利于圖像中場景的自動分析和SAR圖像的理解，使得圖像的解譯工作變得困難，尤其是對SAR圖像中的點目標和邊緣的影響非常明顯。因此斑點噪聲的去除對SAR圖像的后續(xù)處理例如邊緣檢測是非常重要的。近年來，SAR圖像相干斑噪聲抑制技術(shù)飛速發(fā)展，可分為成像前多視平滑技術(shù)和成像后濾波技術(shù)兩大類，后處理技術(shù)中空域濾波和變換域方法得到了廣泛地研究?？沼驗V波如Lee，Gamma-MAP等對SAR圖像進行去噪實現(xiàn)過程簡單易于理解，一定程度上可以有效的減弱噪聲的影響，但此類濾波器都不同程度的對圖像產(chǎn)生了過平滑作用，使圖像變模糊，細節(jié)信息丟失嚴重。近年來，變換域的方法也得到了發(fā)展，其中基于小波變換的方法成為研究的焦點，如Crouse等人將小波變換與隱馬爾科夫模型聯(lián)系起來，提出了小波域隱馬爾科夫模型，開辟了多尺度變換域統(tǒng)計信號處理這一新的研究領(lǐng)域。但對于具有線狀奇異的目標函數(shù)，例如圖像中的邊緣，小波系數(shù)不再稀疏，因此小波在SAR圖像去噪中并不能夠很好地保持圖像中的細節(jié)信息。為了解決小波基函數(shù)各向同性帶來的問題，Donoho等人提出了多尺度幾何分析，其能夠有效表示和處理高維奇異函數(shù)，并在圖像處理得到了廣泛的應(yīng)用，Contouriet便是其中一種分析工具。隨著多尺度幾何分析的發(fā)展，眾多學(xué)者提出了多尺度幾何變換域下的HMT模型和各種改善的HMT模型。Contourlet變換作為新的多尺度信號表示方法，除了具有小波變換的優(yōu)點外，還具有多方向和各向異性的特點。以Contourlet系數(shù)的統(tǒng)計特征為基礎(chǔ)，Po和Do提出了Contourlet域用于圖像去噪的HMT模型，Contourlet域HMT模型和小波域HMT模型相比，不但可以描述尺度間的相關(guān)性，而且可以描述不同方向間系數(shù)的相關(guān)性。但這種模型只捕獲了微小的尺度內(nèi)的依賴性，用于SAR圖像去噪時，一些細節(jié)、紋理信息被平滑，同質(zhì)區(qū)域內(nèi)的噪聲去除不夠徹底導(dǎo)致同質(zhì)區(qū)域等效視數(shù)低，去噪效果并不理想。
發(fā)明內(nèi)容本發(fā)明的目的在于克服上述已有方法存在的細節(jié)、紋理信息被平滑及同質(zhì)區(qū)域等效視數(shù)低的缺點，提出一種基于隱馬爾可夫樹模型的SAR圖像去噪方法，以達到更徹底的去噪效果，并很好的保持圖像的細節(jié)和紋理信息。本發(fā)明的技術(shù)方案是以Contourlet系數(shù)鄰域間和尺度間都具有很大的相關(guān)性為基礎(chǔ)，利用隱馬爾可夫樹模型HMT模型來捕獲系數(shù)尺度間的依賴性，利用背景隱馬爾科夫模型CHMM來捕獲系數(shù)鄰域內(nèi)的依賴性，把HMT和CHMM結(jié)合起來建立Contourlet域改進的HMT統(tǒng)計模型，模型結(jié)構(gòu)如圖3所示，用該統(tǒng)計模型對SAR圖像進行去噪處理。具體實現(xiàn)步驟如下-(1)對輸入SAR圖像依次進行對數(shù)變換和Contourlet變換，獲得不同方向的Contourlet變換系數(shù)；(2)采用單方向傳遞HMT模型和多方向傳遞HMT模型分別對不同方向的Contourlet系數(shù)進行尺度間HMT建模；(3)用EM算法對建立的HMT模型進行訓(xùn)練，得到均值、方差、狀態(tài)轉(zhuǎn)移概率和狀態(tài)概率最優(yōu)估計參數(shù)，并利用該最優(yōu)估計參數(shù)按照貝葉斯最小均方誤差準則對Contouriet系數(shù)進行萎縮，得到初始無斑系數(shù)；(4)對初始無斑系數(shù)建立尺度內(nèi)的背景隱馬爾科夫模型CHMM，利用EM算法對尺度內(nèi)的CHMM進行訓(xùn)練，得到尺度內(nèi)的估計參數(shù)，利用尺度內(nèi)的估計參數(shù)對初始無斑系數(shù)進行萎縮，得到最終無斑系數(shù)；(5)對最終無斑系數(shù)依次進行Contourlet逆變換和指數(shù)變換，得到初步的去噪圖像/A(6)采用各向異性擴散方法對原始SAR圖像和去噪圖像/7的差值圖像進行去噪處理，得到二次去噪圖像/2;'(7)用模極大值方法對初步去噪圖像和二次去噪圖像/2進行融合-(8)對融合后的圖像用CycleSpi皿ing方法進行旋轉(zhuǎn)平移處理，獲得最終去噪圖像，并輸出。本發(fā)明與現(xiàn)有技術(shù)相比具有如下優(yōu)點1)本發(fā)明充分利用了Contourlet系數(shù)之間的相關(guān)性傳統(tǒng)的HMT去噪方法，只考慮了Contourlet系數(shù)微小的尺度內(nèi)相關(guān)性，沒有充分利用系數(shù)之間的相關(guān)性，導(dǎo)致去噪圖像的一些細節(jié)、紋理信息被平滑，圖像噪聲去除的不徹底，等效視數(shù)較低，本發(fā)明提出的改進HMT算法克服了系數(shù)相關(guān)性沒有充分利用的缺點，利用了Contourlet系數(shù)尺度間、尺度內(nèi)和方向內(nèi)的相關(guān)性，得到了較好的去噪結(jié)果。2、本發(fā)明提供了更好的相干斑抑制結(jié)果傳統(tǒng)的空域濾波不同程度的對圖像產(chǎn)生了過平滑作用，細節(jié)信息丟失嚴重，傳統(tǒng)的Contourlet域去噪算法，由于Contourlet變換非平移不變性，會出現(xiàn)劃痕即偽吉布斯現(xiàn)象；本發(fā)明采用各向異性擴散方法對差值圖像進行去噪，補充了去噪中丟失的圖像細節(jié)信息，使得融合后圖像的細節(jié)特征也得到了較好的保持；引入了CycleSpinning，克服Contourlet變換非平移不變性帶來的偽吉布斯現(xiàn)象，使得同質(zhì)區(qū)域內(nèi)噪聲去除的徹底，用于評價去噪后同質(zhì)區(qū)域內(nèi)噪聲去除多少的等效視數(shù)得到提高。仿真實驗的結(jié)果表明，本發(fā)明方法更好的保持了圖像的細節(jié)和紋理信息，同質(zhì)區(qū)域內(nèi)噪聲去除的較徹底。圖1是本發(fā)明的主要操作過程示意圖；圖2是本發(fā)明Contourlet變換域尺度間HMT模型圖；圖3是本發(fā)明改進的Contourlet變換域HMT模型圖；圖4是用本發(fā)明與已有方法對SARI圖像去噪效果對比圖；圖5是用本發(fā)明與已有方法對SAR2圖像去噪效果對比圖。具體實施例方式參照圖l,本發(fā)明的具體實現(xiàn)步驟如下步驟l:對SAR圖像作對數(shù)變換，將乘性噪聲轉(zhuǎn)換成為加性高斯白噪聲來處理logj=logz+logx,其中；;表示輸入sar圖像，z表示噪聲圖像，^表示不含噪聲的圖像，采用方向不相同的兩組方向濾波器，該兩組方向濾波器的方向分別為4，4,4和4，8，8，對經(jīng)對數(shù)變換后的數(shù)據(jù)進行Contourlet分解，得到的Contourlet變換系數(shù)分別為乂和J^。步驟2:對變換系數(shù)Yi建立單方向傳遞HMT模型，如圖2(a)所示，黑色方形為父節(jié)點，四個空方形是其子節(jié)點，該模型的參數(shù)集為1;對變換系數(shù)^2建立多方向傳遞HMT模型，如圖2(b)所示，黑色方形為父節(jié)點，四個空方形是其子節(jié)點，每一個父節(jié)點對應(yīng)的四個子節(jié)點分布在兩個不同的方向子帶內(nèi)，該模型的參數(shù)集為2。步驟3:對包括均值、方差、狀態(tài)轉(zhuǎn)移概率及根節(jié)點的初始狀態(tài)概率的兩個參數(shù)集0i和2進行初始化，再用EM算法對單方向傳遞HMT模型和多方向傳遞HMT模型進行訓(xùn)練，得到A和2的最優(yōu)估計參數(shù);利用該最優(yōu)估計參數(shù)按照貝葉斯最小均方誤差準則對Contourlet系數(shù)進行萎縮，得到初始無斑系數(shù)，萎縮公式如下其中m表示節(jié)點的狀態(tài)，取值為1或2，p(5^,^m)是J尺度，A:方向，Z位置的節(jié)點狀態(tài)概率，《".,M是狀態(tài)為m，/尺度，*方向，H立置的節(jié)點的方差，a^^^是尺度/的噪聲方差，".是/尺度，*方向，/位置的Contourlet系數(shù)，"w是萎縮后的初始無斑系數(shù)。步驟4:對初始無斑系數(shù)建立尺度內(nèi)的背景隱馬爾科夫模型CHMM，如圖3所示，黑色的節(jié)點是Contourlet變換系數(shù)，白色的圓形節(jié)點是Contourlet的系數(shù)的隱狀態(tài)變量，狀態(tài)數(shù)為2，菱形節(jié)點為Contouriet的系數(shù)的背景變量，一個父節(jié)點的狀態(tài)對應(yīng)4個子節(jié)點的狀態(tài)，y尺度，/t方向，/位置的節(jié)點的背景,是由鄰域內(nèi)八個節(jié)點的局部能量平均值決定的，如下式其中，A^,為系數(shù)c,t,的八領(lǐng)域內(nèi)Contourlet系數(shù)的局部能量平均值，《是尺度y內(nèi)系數(shù)的平均能量；利用EM算法對尺度內(nèi)的CHMM進行訓(xùn)練，先將步驟(3)得到的均值、方差和狀態(tài)概率最優(yōu)估計參數(shù)作為CHMM訓(xùn)練的初始參數(shù)，然后進行訓(xùn)練，得到尺度內(nèi)的估計參數(shù),利用尺度內(nèi)的估計參數(shù)對初始無斑系數(shù)進行萎縮，得到最終無斑系數(shù)。步驟5:對第(4)步得到的最終無斑系數(shù)依次進行Contourlet逆變換和指數(shù)變換，得到初步的去噪圖像//。步驟6:用各向異性擴散方法對原始SAR圖像和去噪圖像//的差值圖像進行去噪處理，得到二次去噪圖像/2。對差值圖像利用梯度算子來區(qū)分邊緣和噪聲，然后對差值圖像進行擴散去噪，擴散去噪方程如下:其中《是進行各向異性擴散之前的差值圖像像素值，《1是進行=《"Zr，s=—1(r，"-(O,O)2,2r+s10次各向異性擴散后的差值圖像像素值，r是尺度參數(shù)，取值范圍是0"<1/6,r,s是目標點鄰域系數(shù)坐標相對于目標點的位置，g是擴散率函數(shù)，在上式中使用的擴散率函數(shù)是Perona-Malik擴散率gW=1/(1+jc2/—)，〃是擴散率中的對比參數(shù)。對于AT的選取，根據(jù)多次實驗，通過比較發(fā)現(xiàn)，廳取50時效果最好。步驟7:用模極大值方法對初步去噪圖像/7和二次去噪圖像/2進行融合，具體步驟如下-(7a)對初步去噪圖像/7和二次去噪圖像/2分別進行Contourlet分解，選擇的濾波器方向均為4，4，4,得到初步去噪圖像//和二次去噪圖像/2的Contourlet變換系數(shù)；(7b)對兩個去噪圖像的Contourlet變換系數(shù)的細節(jié)項進行模極大值融合得到融合圖像的Contourlet系數(shù)的細節(jié)項，兩個去噪圖像的Contourlet變換系數(shù)的模糊項直接相加得到融合圖像Contourlet系數(shù)的模糊項；(7c)對融合圖像的Contourlet系數(shù)進行Contourlet反變換得到融合后的圖像。步驟8:對融合后的圖像用CycleSpi皿ing方法進行旋轉(zhuǎn)平移處理，具體步驟如下(8a)通過行和列方向上的每組平移量得到一個不同的去噪結(jié)果其中，S是循環(huán)平移算子，T是Contourlet變換算子，T-1是Contourlet逆變換算子，A為Contourlet域HMT去噪算子，下標-/，-_/，/和_/分別為行和列方向上的平移量；(8b)對所有去噪結(jié)果進行線性平均,得到抑制偽吉布斯現(xiàn)象的去噪結(jié)果-《1《2i=0，_/=0，其中，《，i^分別表示行和列方向上的最大平移量，均為4。步驟9:把步驟(8)中得到的去噪結(jié)果&作為最終去噪結(jié)果，并輸出。以下給出仿真實驗結(jié)果分析，以進一步說明本發(fā)明的效果1、仿真條件本發(fā)明中的兩幅測試SAR圖像均是從美國SandiaLabs網(wǎng)站上獲得的HorsetrackSAR圖像上截取的，大小均為256x256，分別命名為SARI和SAR2。采用Gamma—MAP濾波方法、Lee濾波方法、小波域HMT去噪方法、Contourlet域HMT去噪方法和本發(fā)明方法1改進單方向傳遞HMT模型及本發(fā)明方法2改進多方向傳遞HMT模型對SARI和SAR2圖像進行去噪處理。2、仿真結(jié)果分析仿真的結(jié)果如圖4和圖5，其中圖4(a)為SAR1原圖像，圖4(b)為SARI圖像經(jīng)Gamma一MAP濾波后的去噪圖像，圖4(c)為SARI圖像經(jīng)Lee濾波后的去噪圖像，圖4(d)為SARI圖像經(jīng)小波域HMT去噪后的去噪圖像圖4(e)為SARI圖像經(jīng)Contourlet域HMT去噪后的去噪圖像，圖4(f)為SARI圖像經(jīng)本發(fā)明改進單方向傳遞HMT模型方法去噪后的去噪圖像，圖4(g)為SARI圖像經(jīng)本發(fā)明改進多方向傳遞HMT模型方法去噪后的去噪圖像。圖5(a)為SAR2原圖像，圖5(b)為SAR2圖像經(jīng)Gamma—MAP濾波后的去噪圖像，圖5(c)為SAR2圖像經(jīng)Lee濾波后的去噪圖像，圖5(d)為SAR2圖像經(jīng)小波域HMT去噪后的去噪圖像，圖5(e)為SAR2圖像經(jīng)Contourlet域HMT去噪后的去噪圖像，圖5(f)為SAR2圖像經(jīng)本發(fā)明改進單方向傳遞HMT模型方法去噪后的去噪圖像，圖5(g)為SAR2圖像經(jīng)本發(fā)明改進多方向傳遞HMT模型方法去噪后的去噪圖像。從圖4(b)和圖4(c)可以看出，Gamma_MAP濾波和Lee濾波對SARI圖像的斑點噪聲去除的較好，但SAR1圖像的邊緣被模糊，圖像細節(jié)信息丟失較嚴重；由圖4(d)和圖4(e)，從同質(zhì)區(qū)域內(nèi)噪聲的去除可以看出Contourlet域HMT方法要優(yōu)于小波域HMT方法，但這兩種方法對SARI圖像去噪后的同質(zhì)區(qū)域內(nèi)噪聲去除的都不夠徹底。從圖4(f)和圖4(g)，可以看出SAR1去噪后同質(zhì)區(qū)域內(nèi)噪聲去除的很徹底，細節(jié)信息、邊緣信息及紋理信息得到了很好的保持。從圖5可以看出，本發(fā)明在徹底去除噪聲的同時更多的保持了SAR2圖像的邊緣和紋理細節(jié)；可見，本發(fā)明在細節(jié)，紋理保持方面有相當(dāng)大的優(yōu)勢。對SAR圖像去噪性能的評價指標包括等效視數(shù)KV丄、均值及均值比，這些評價標準的定義如下(i)等效視數(shù)(五M:):五^￡=」~其中"和cr分別為相干斑抑制后圖像中某同質(zhì)區(qū)域的均值和標準差，對于強度圖像^=1，對于幅度圖像^=4"-1。它是衡量圖像中相干斑相對強度的一種指標，相干斑抑制程度越深，等效視數(shù)越大。對于SAR1圖像選擇尺寸為40x40，50x50，40x40的同質(zhì)區(qū)域作為計算五7仏的測試數(shù)據(jù)，如圖4(a)中的方形框標示區(qū)域；對于SAR2圖像選擇尺寸為30x30,30x30，30x30的同質(zhì)區(qū)域作為計算^V丄的測試數(shù)據(jù)，如圖5(a)中的方形框標示區(qū)域。(2)圖像的均值均值反映圖像的平均亮度，相干斑抑制前后的圖像均值要基本保持一致。(3)均值比原圖像/與相干斑抑制后圖像"均值之比，理想情況下比率y為純斑點噪聲的統(tǒng)計幅值1。實際的r值與l相差越大，表明輻射度失真越厲害。SAR1和SAR2圖像去噪后的評價指標如表1和表2所示表1對SAR1圖像楚e行去噪后去噪圖像的評'價指標比較SAR1圖像區(qū)域l區(qū)域2冊丄區(qū)域3五m:均值均值比Gamma—MAP217.1888122.2114107.955983.12470.9913Lee濾波217.1888122.2114107.955983.66540.9978小波HMT78.818558.168436.542980.31430.9577ContourletHMT120.216080.063952.256079.94570.9533本發(fā)明方法1327.6217158.3174140.774283,78110.9992本發(fā)明方法2338.5326153.0471143.338783,77410.999114表2對SAR2圖像進行去噪后去噪圖像的i5<table>tableseeoriginaldocumentpage15</column></row><table>價指標比較從表1及表2可以看出(a)本發(fā)明的等效視數(shù)值均比其它幾種經(jīng)典方法的等效視數(shù)值大，表明本發(fā)明對SAR圖像相干斑噪聲抑制程度較深，斑點噪聲去除的較徹底，去噪效果較好；(b)本發(fā)明的均值和相應(yīng)原圖像的均值基本相同，表明本發(fā)明去噪后圖像的平均亮度基本保持不變；(c)本發(fā)明的均值比與1的差距比其它經(jīng)典方法小，表明本發(fā)明去噪過程中輻射度失真很小，細節(jié)信息和紋理信息保持的較好。綜上，本發(fā)明方法既較好地保持了SAR圖像中的點目標及其邊緣等細節(jié)特征和紋理信息，也大幅度的減少了SAR圖像同質(zhì)區(qū)域內(nèi)的斑點噪聲。權(quán)利要求1、一種基于隱馬爾科夫樹模型HMT的SAR圖像去噪方法，包括如下步驟(1)對輸入SAR圖像依次進行對數(shù)變換和Contourlet變換，獲得不同方向的Contourlet變換系數(shù)；(2)采用單方向傳遞模型和多方向傳遞模型分別對不同方向的Contourlet系數(shù)進行尺度間HMT建模；(3)用EM算法對建立的HMT模型進行訓(xùn)練，得到均值、方差、狀態(tài)轉(zhuǎn)移概率和狀態(tài)概率最優(yōu)估計參數(shù)，并利用該最優(yōu)估計參數(shù)按照貝葉斯最小均方誤差準則對Contourlet系數(shù)進行萎縮，得到初始無斑系數(shù)；(4)對初始無斑系數(shù)建立尺度內(nèi)的背景隱馬爾科夫模型CHMM，利用EM算法對尺度內(nèi)的CHMM進行訓(xùn)練，得到尺度內(nèi)的估計參數(shù)，利用尺度內(nèi)的估計參數(shù)對初始無斑系數(shù)進行萎縮，得到最終無斑系數(shù)；(5)對最終無斑系數(shù)依次進行Contourlet逆變換和指數(shù)變換，得到初步的去噪圖像I1；(6)采用各向異性擴散方法對原始SAR圖像和去噪圖像I1的差值圖像進行去噪處理，得到二次去噪圖像I2；(7)用模極大值方法對初步去噪圖像I1和二次去噪圖像I2進行融合；(8)對融合后的圖像用CycleSpinning方法進行旋轉(zhuǎn)平移處理，獲得最終去噪圖像，并輸出。2、根據(jù)權(quán)利要求i所述的的SAR圖像去噪方法，萁中步驟(1)中所述的Contourlet變換，采用方向不相同的兩組方向濾波器，該兩組方向濾波器的方向分別為4，4，4禾Q4，8，8，得到的Contourlet變換系數(shù)分別為兀和h。3、根據(jù)權(quán)利要求1所述的SAR圖像去噪方法，其中步驟(2)所述的采用單方向傳遞模型和多方向傳遞模型分別對不同方向的Contourlet系數(shù)進行尺度間HMT建模，是對變換系數(shù)力建立單方向傳遞模型，該模型的參數(shù)集為1;對變換系數(shù)^建立多方向傳遞模型，該模型的參數(shù)集為2。4、根據(jù)權(quán)利要求1所述的SAR圖像去噪方法，其中步驟(3)中所述的用EM算法對建立的HMT模型進行訓(xùn)練，是先對包括均值和方差，狀態(tài)轉(zhuǎn)移概率，根節(jié)點的初始狀態(tài)概率的兩個參數(shù)集^和02進行初始化，再用EM算法對HMT模型進行訓(xùn)練，得到A和2的最優(yōu)估計參數(shù)。5、根據(jù)權(quán)利要求1所述的SAR圖像去噪方法，其中步驟(4)中所述的利用EM算法對尺度內(nèi)的CHMM進行訓(xùn)練，是將步驟(3)得到的均值、方差、狀態(tài)轉(zhuǎn)移概率和狀態(tài)概率最優(yōu)估計參數(shù)作為CHMM訓(xùn)練的初始參數(shù)。6、根據(jù)權(quán)利要求1所述的SAR圖像去噪方法，其中步驟(8)中所述的對融合后的圖像用CycleSpinning方法進行旋轉(zhuǎn)平移處理，按如下步驟進行(6a)通過行和列方向上的每組平移量得到一個不同的去噪結(jié)果、其中，S是循環(huán)平移算子，T是Contourlet變換算子，T-1是Contourlet逆變換算子，A為Contourlet域HMT去噪算子，下標-/,，/和_/分別為行和列方向上的平移量；(6b)對所有去噪結(jié)果進行線性平均,得到抑制偽吉布斯現(xiàn)象的去噪其中，《，《2分別表示行和列方向上的最大平移量。全文摘要本發(fā)明公開一種基于隱馬爾科夫樹模型的SAR圖像去噪方法，它涉及圖像處理領(lǐng)域，主要解決圖像細節(jié)、紋理信息被平滑及同質(zhì)區(qū)域等效視數(shù)低的問題。其步驟為(1)對數(shù)變換和Contourlet分解；(2)對Contourlet系數(shù)進行HMT建模，并訓(xùn)練；(3)利用估計參數(shù)對Contourlet系數(shù)進行校正；(4)建立背景隱馬爾科夫訓(xùn)練模型，用該估計參數(shù)重新校正系數(shù)，并進行Contourlet逆變換和指數(shù)變換，得到一次去噪圖像；(5)對差值圖像進行去噪，得到二次去噪圖像；(6)對兩次去噪圖像進行融合，并對融合后的圖像進行旋轉(zhuǎn)平移，獲得最終去噪圖像并輸出。本發(fā)明較好地保持了SAR圖像的細節(jié)和紋理信息，減少SAR圖像同質(zhì)區(qū)域內(nèi)斑點噪聲，可用于對SAR圖像進行去噪處理。文檔編號G06T5/00GK101685158SQ20091002317公開日2010年3月31日申請日期2009年7月2日優(yōu)先權(quán)日2009年7月2日發(fā)明者彪侯,張向榮,焦李成,爽王,田福苓,馬文萍申請人:西安電子科技大學(xué)