一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法及系統(tǒng)的制作方法
【專利摘要】本發(fā)明提供了一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法及系統(tǒng),用于處理干涉合成孔徑聲納相位噪聲嚴(yán)重的問(wèn)題,所述方法包含:步驟101)經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)主接收陣和輔接收陣圖像之間的相位差別,校正主輔圖像之間的相位偏移,完成主接收陣圖像和輔接收陣圖像對(duì)的配準(zhǔn);步驟102)對(duì)完成配準(zhǔn)的主接收陣圖像和輔接收陣圖像進(jìn)行相干處理,得到干涉相位圖;步驟103)將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,再依據(jù)圓周期中值濾波對(duì)干涉相位圖進(jìn)行濾波;步驟104)基于濾波后的干涉相位值,使用最小二乘解纏,得到干涉相位,該相位取值位于(-∞,+∞)之間;步驟105)根據(jù)纏繞相位和InSAS模型進(jìn)行高程重建,完成對(duì)干涉合成孔徑聲納信號(hào)的相位處理。
【專利說(shuō)明】一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法及系統(tǒng)
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及合成孔徑聲納信號(hào)處理領(lǐng)域,特別涉及一種干涉合成孔徑聲納相位處理方法及系統(tǒng)。
【背景技術(shù)】
[0002]合成孔徑聲納(Synthetic Aperture Sonar,SAS)是一種水下高分辨率成像技術(shù)。干涉合成孔徑聲納(Interferometric Synthetic Aperture Sonar, InSAS)在合成孔徑聲納的基礎(chǔ)上發(fā)展而來(lái),給水下三維成像技術(shù)提供了一種新的選擇。干涉合成孔徑聲納使用兩副接收陣,每副陣分別進(jìn)行合成孔徑聲納成像,并使用比相測(cè)高的方法得到目標(biāo)場(chǎng)景的數(shù)字高程模型。聲波傳輸速度較慢,聲納平臺(tái)運(yùn)動(dòng)穩(wěn)定性較差,水體存在多徑等現(xiàn)象。這些惡劣的條件增加了 InSAS相位處理方法的復(fù)雜性。
[0003]干涉合成孔徑技術(shù)在雷達(dá)上已有一定應(yīng)用,但在聲納中發(fā)展相對(duì)滯后。各國(guó)的InSAS技術(shù)還存在試驗(yàn)階段。本文在研究InSAS現(xiàn)有算法的基礎(chǔ)上,給出了一種簡(jiǎn)單易行的InSAS相位處理方法。
【發(fā)明內(nèi)容】
[0004]本發(fā)明的目的在于,為克服現(xiàn)有技術(shù)的不足,從而提供一種干涉合成孔徑聲納相位處理方法及系統(tǒng)。
[0005]為了克服上述問(wèn)題,本發(fā)明提供一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,用于處理干涉合成孔徑聲納相位噪聲嚴(yán)重的問(wèn)題,所述方法包含:
[0006]步驟101)經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)主接收陣和輔接收陣圖像之間的相位差別,校正主輔圖像之間的相位偏移,完成主接收陣圖像和輔接收陣圖像對(duì)的配準(zhǔn);
[0007]步驟102)對(duì)完成配準(zhǔn)的主接收陣圖像和輔接收陣圖像進(jìn)行相干處理,得到干涉相位圖;
[0008]步驟103)將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,再依據(jù)圓周期中值濾波對(duì)干涉相位圖進(jìn)行濾波;
[0009]步驟104)基于濾波后的干涉相位值,使用最小二乘解纏,得到干涉相位,該相位取值位于(-OO,+ OO )之間;
[0010]步驟105)根據(jù)纏繞相位和InSAS模型進(jìn)行高程重建,完成對(duì)干涉合成孔徑聲納信號(hào)的相位處理。
[0011]上述技術(shù)方案中,所述步驟101)之前還包含獲取主接收陣圖像和輔接收陣圖像的步驟,該步驟具體為:
[0012]使用一副發(fā)射陣列、兩副平行放置的SAS陣列;所述發(fā)射陣以一定的脈沖重復(fù)周期發(fā)射信號(hào),到達(dá)目標(biāo)場(chǎng)景返回,兩副接收陣以相近的視角接收回波信號(hào)。分別經(jīng)過(guò)頻域成像算法,得到目標(biāo)場(chǎng)景的主輔圖像對(duì)。
[0013]上述技術(shù)方案中,所述步驟103)和步驟104)之間還包含:使用分支割線方法補(bǔ)償由于相位的周期特性弓I起的干涉相位模糊的步驟。
[0014]上述技術(shù)方案中,所述步驟102)進(jìn)一步包含:
[0015]步驟102-1)基于配準(zhǔn)后的主接收陣圖像和輔接收陣圖像得到復(fù)干涉圖;
[0016]步驟102-2)對(duì)復(fù)干涉圖求取相位得到干涉相位圖。
[0017]上述技術(shù)方案中,所述步驟104進(jìn)一步包含:
[0018]首先,掃描干涉圖,判斷每一個(gè)像素點(diǎn)是否為殘余點(diǎn);搜索干涉圖中未被支切線連接的殘余點(diǎn),找到符合條件的點(diǎn)后,在以該點(diǎn)為中心大小為3 X 3的搜索窗口內(nèi),搜索其他殘余點(diǎn);若未搜索到,則設(shè)置支切線過(guò)程結(jié)束;若在窗口內(nèi)搜索到殘余點(diǎn),則不管該殘余點(diǎn)的殘數(shù)是否和搜索中心點(diǎn)的殘數(shù)相反,將兩殘余點(diǎn)用支切線相連,如果搜索到的殘余點(diǎn)未被其他支切線連接,則求支切線的殘數(shù),即支切線上各殘余點(diǎn)的殘數(shù)和;若殘數(shù)為0,整個(gè)支切線上正負(fù)殘余點(diǎn)相互抵消,則重新搜索干涉圖中未被支切線連接的殘余點(diǎn)。若在窗口內(nèi)搜索到干涉圖邊界,則用支切線連接搜索中心點(diǎn)與干涉圖邊界,將該支切線殘數(shù)置零,重新搜索干涉圖中未被支切線連接的殘余點(diǎn);若未搜索到殘余點(diǎn)或邊界,或搜索到殘余點(diǎn)后支切線上的殘數(shù)不為零,則以支切線上其余殘余點(diǎn)為中心搜索,直至搜索完整個(gè)支切線上所有殘余點(diǎn),則將搜索窗口大小增加2,繼續(xù)搜索。
[0019]基于上述方法本發(fā)明提供了一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理系統(tǒng),用于處理干涉合成孔徑聲納相位噪聲嚴(yán)重的問(wèn)題,所述系統(tǒng)包含:
[0020]配準(zhǔn)模塊,用于經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)主接收陣和輔接收陣圖像之間的相位差別,校正主輔圖像之間的相位偏移,完成主接收陣圖像和輔接收陣圖像對(duì)的配準(zhǔn);
[0021]干涉相位圖獲取模塊,用于對(duì)完成配準(zhǔn)的主接收陣圖像和輔接收陣圖像對(duì)進(jìn)行相干處理,得到干涉相位圖;
[0022]濾波模塊,用于將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,再依據(jù)圓周期中值濾波對(duì)干涉相位圖進(jìn)行濾波;
[0023]纏繞相位獲取模塊,用于基于濾波后的干涉相位值,使用最小二乘解纏,得到干涉相位,該相位取值位于(-OO,+ OO )之間;
[0024]建模模塊,用于根據(jù)纏繞相位和InSAS模型進(jìn)行高程重建,完成對(duì)干涉合成孔徑聲納信號(hào)的相位處理。
[0025]上述技術(shù)方案中,所述系統(tǒng)還包含:
[0026]主輔圖像獲取模塊,用于使用一副發(fā)射陣列、兩副平行放置的SAS陣列;所述發(fā)射陣以一定的脈沖重復(fù)周期發(fā)射信號(hào),到達(dá)目標(biāo)場(chǎng)景返回,兩副接收陣以相近的視角接收回波信號(hào)。分別經(jīng)過(guò)頻域成像算法,得到目標(biāo)場(chǎng)景的主輔圖像對(duì)。
[0027]上述技術(shù)方案中,所述系統(tǒng)還包含:
[0028]補(bǔ)償模塊,用于使用分支割線方法補(bǔ)償由于相位的周期特性引起的干涉相位模糊的步驟。
[0029]上述技術(shù)方案中,所述干涉相位圖獲取模塊進(jìn)一步包含:
[0030]復(fù)干涉圖獲取模塊,用于基于配準(zhǔn)后的主接收陣圖像和輔接收陣圖像得到復(fù)干涉圖;
[0031]干涉相位圖獲取模塊,用于對(duì)復(fù)干涉圖求取相位得到干涉相位圖。[0032]與現(xiàn)有技術(shù)相比,本發(fā)明的技術(shù)優(yōu)勢(shì)在于:
[0033]本發(fā)明建立了一套簡(jiǎn)單有效的干涉合成孔徑聲納信號(hào)的相位處理方法,該處理方法該適合干涉合成孔徑聲納相位噪聲嚴(yán)重的特點(diǎn),切實(shí)可行。
【專利附圖】
【附圖說(shuō)明】
[0034]圖1本發(fā)明提供的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法的流程圖;
[0035]圖2本發(fā)明的干涉合成孔徑聲納干涉模型;
[0036]圖3是本發(fā)明技術(shù)方案采用的解纏繞方法的流程圖。
【具體實(shí)施方式】
[0037]下面結(jié)合附圖對(duì)本發(fā)明的技術(shù)方案進(jìn)行詳細(xì)的闡述和說(shuō)明。
[0038]為了在水下二維高精度成像的基礎(chǔ)上獲得高程信息,本文給出了一種InSAS處理方法。InSAS系統(tǒng)使用一副發(fā)射陣列、兩副平行放置的SAS陣列。發(fā)射陣以一定的脈沖重復(fù)周期發(fā)射信號(hào),到達(dá)目標(biāo)場(chǎng)景返回,兩副接收陣以相近的視角接收回波信號(hào)。分別經(jīng)過(guò)頻域成像算法,得到目標(biāo)場(chǎng)景的主輔圖像對(duì)。經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)兩幅圖像的相位差別,校正主輔圖像之間的偏移。主輔圖像對(duì)為復(fù)圖像,經(jīng)過(guò)相干處理,可以得到干涉圖。干涉圖中包含大量由于運(yùn)動(dòng)誤差等引起的相位噪聲,使用圓周期均值方法,對(duì)干涉圖濾波,可以提高干涉圖質(zhì)量。使用最小二乘解纏方法,可以恢復(fù)纏繞相位,將相位由(_π, Ji]恢復(fù)到(-OO,+ OO )之間。
[0039]本發(fā)明目的在于建立一套切實(shí)可行的InSAS處理方法,由InSAS主輔圖像得到高
程信息。
`[0040]本發(fā)明提出的一種InSAS相位處理方法,主要步驟包括:
[0041]如圖1所示,該圖給出了本發(fā)明提供的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法的流程圖,所述方法包含如下步驟:
[0042]1.配準(zhǔn)
[0043]本發(fā)明首先使用雙線性插值方法對(duì)主輔圖像進(jìn)行10倍插值。在主輔圖像中每隔100行100列選取一個(gè)點(diǎn)作為控制點(diǎn)。以相鄰像素間的平均相位差為評(píng)價(jià)標(biāo)準(zhǔn)進(jìn)行干涉圖配準(zhǔn)
[0044]/ = Σ Σ (\ρ O +./) — P 仏./)1 + \ρ./.+ l)-P (i,./)|) / 2
* 1(I)
[0045]式中,i為目標(biāo)點(diǎn)在干涉圖中的方位向坐標(biāo),j為距離向坐標(biāo)。p(i,j)為圖像中坐標(biāo)為(i,j)點(diǎn)處的干涉相位值。|p(i,j) I表示求取相位值P (i,j)的絕對(duì)值。當(dāng)f取值最小時(shí),就得到了控制點(diǎn)的配準(zhǔn)偏移量。根據(jù)配準(zhǔn)偏移量使用二次多項(xiàng)式求取主輔圖像之間的變換模型
[0046]IX = ^ + % + 叫,+ V + 明,+(2)
[7=4,+ f\x + b2y + byx_ + b4xy + b5jr
[0047]根據(jù)變換模型,對(duì)主輔圖像雙線性插值得到配準(zhǔn)干涉圖相對(duì)。
[0048]其中,X和y是變換之前的距離向和方位向坐標(biāo)。X和Y是變換之后的距離向和方位向坐標(biāo)。a。、a” a2、a3、a4、a5>b0> b2、b3、b4和b5是二次多項(xiàng)式變換模型的參數(shù)。[0049]2.生成干涉圖
[0050]經(jīng)過(guò)配準(zhǔn)步驟后,通過(guò)將主圖像和中各像素乘以輔圖像中對(duì)應(yīng)像素的共軛,就可以得到復(fù)干涉圖。對(duì)復(fù)干涉圖求取相位,得到干涉相位圖。
[0051]3.濾波
[0052]干涉圖具有圓周期性,可以將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,進(jìn)行圓周期中值濾波:
[0053].? I xii^n2(exp[/?^(/?,?)]] re ?)
爐(jr,j)=M/^ —— ΣΣ ai'g1-jiTi— r+ar8.);) r
I Μ.Ν χ (Λ/ I?>2 Kr (..V 1) 2 I /(**,>') JJ (3)
[0054]其中
[0055]
-V!{M I) 2ι-!(Λr 1)-.2
/(X,J7)= XΣ argfexp[>^(OT,?)]}(4)
m χ (Χ? I).2 ? r-(.V-1),2
為原始的丨:涉_相位值,為濾波后的P涉相位值,M和況為濾波窗口距離向和方位向尺寸,j為虛數(shù)單位,exp{}為指數(shù)運(yùn)算,arg{}為角度運(yùn)算,w{}為相位纏繞算子。圓周期均值濾波的效果,取決于濾波窗口,濾波窗口越大,濾波結(jié)果越平滑,但分辨率越低。
[0056]4.解纏繞
[0057]如圖3所示,干涉相位的取值在(_π,π]之間,這里使用分支割線方法,恢復(fù)干涉相位測(cè)量值與真實(shí)值存在的2 π整數(shù)倍的差別。首先,掃描干涉圖,判斷每一個(gè)像素點(diǎn)是否為殘余點(diǎn)。搜索干涉圖中未被支切線連接的殘余點(diǎn)。找到符合條件的點(diǎn)后,在以該點(diǎn)為中心大小為3X3的搜索窗口內(nèi),搜索其他殘余點(diǎn)。若未搜索到,則設(shè)置支切線過(guò)程結(jié)束。若在窗口內(nèi)搜索到殘余點(diǎn),則不管該殘余點(diǎn)的殘數(shù)是否和搜索中心點(diǎn)的殘數(shù)相反,將兩殘余點(diǎn)用支切線相連,如果搜索到的殘余點(diǎn)未被其他支切線連接,則求支切線的殘數(shù),即支切線上各殘余點(diǎn)的殘數(shù)和。若殘數(shù)為0,整個(gè)支切線上正負(fù)殘余點(diǎn)相互抵消,則重新搜索干涉圖中未被支切線連接的殘余點(diǎn)。若在窗口內(nèi)搜索到干涉圖邊界,則用支切線連接搜索中心點(diǎn)與干涉圖邊界,將該支切線殘數(shù)置零,重新搜索干涉圖中未被支切線連接的殘余點(diǎn)。若未搜索到殘余點(diǎn)或邊界,或搜索到殘余點(diǎn)后支切線上的殘數(shù)不為零,則以支切線上其余殘余點(diǎn)為中心搜索,直至搜索完整個(gè)支切線上所有殘余點(diǎn),則將搜索窗口大小增加2,繼續(xù)搜索。
[0058]5.生成數(shù)字高程模型
[0059]得到解纏相位Φ后,根據(jù)InSAS干涉模型,可以進(jìn)行高程重建。如圖3,ΤΧ為發(fā)射陣,RX1和RX2S兩副接收陣。RX1和RX2在垂直航跡平面內(nèi)分開(kāi),距離即為基線B。X為地距向,Y為方位向,Z為高度向。目標(biāo)點(diǎn)P的高度為h。R0為從發(fā)射陣TX到目標(biāo)點(diǎn)P的斜距,R1和R2分別為從目標(biāo)點(diǎn)P到接收陣RX1和RX2的斜距。R1與豎直方向的夾角為θ,基線B與垂直方向的夾角為4。則
I
[0060]AR=^w
2 TT(4)
【權(quán)利要求】
1.一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,用于處理干涉合成孔徑聲納相位噪聲嚴(yán)重的問(wèn)題,所述方法包含: 步驟101)經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)主接收陣和輔接收陣圖像之間的相位差另IJ,校正主輔圖像之間的相位偏移,完成主接收陣圖像和輔接收陣圖像對(duì)的配準(zhǔn); 步驟102)對(duì)完成配準(zhǔn)的主接收陣圖像和輔接收陣圖像進(jìn)行相干處理,得到干涉相位圖; 步驟103)將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,再依據(jù)圓周期中值濾波對(duì)干涉相位圖進(jìn)行濾波; 步驟104)基于濾波后的干涉相位值,使用最小二乘解纏,得到干涉相位,該相位取值位于(-OO,+ OO )之間; 步驟105)根據(jù)纏繞相位和InSAS模型進(jìn)行高程重建,完成對(duì)干涉合成孔徑聲納信號(hào)的相位處理。
2.根據(jù)權(quán)利要求1所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,其特征在于,所述步驟101)之前還包含獲取主接收陣圖像和輔接收陣圖像的步驟,該步驟具體為: 使用一副發(fā)射陣列、兩副平行放置的SAS陣列;所述發(fā)射陣以一定的脈沖重復(fù)周期發(fā)射信號(hào),到達(dá)目標(biāo)場(chǎng)景返回,兩副接收陣以相近的視角接收回波信號(hào)。分別經(jīng)過(guò)頻域成像算法,得到目標(biāo)場(chǎng)景的主輔圖像對(duì)。
3.根據(jù)權(quán)利要求1所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,其特征在于,所述步驟103)和步驟104)之間還包含:使用分支割線方法補(bǔ)償由于相位的周期特性引起的干涉相位模糊的步驟。
4.根據(jù)權(quán)利要求1所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,其特征在于,所述步驟102)進(jìn)一步包含: 步驟102-1)基于配準(zhǔn)后的主接收陣圖像和輔接收陣圖像得到復(fù)干涉圖; 步驟102-2)對(duì)復(fù)干涉圖求取相位得到干涉相位圖。
5.根據(jù)權(quán)利要求1所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理方法,其特征在于,所述步驟104進(jìn)一步包含: 首先,掃描干涉圖,判斷每一個(gè)像素點(diǎn)是否為殘余點(diǎn);搜索干涉圖中未被支切線連接的殘余點(diǎn),找到符合條件的點(diǎn)后,在以該點(diǎn)為中心大小為3X3的搜索窗口內(nèi),搜索其他殘余點(diǎn);若未搜索到,則設(shè)置支切線過(guò)程結(jié)束;若在窗口內(nèi)搜索到殘余點(diǎn),則不管該殘余點(diǎn)的殘數(shù)是否和搜索中心點(diǎn)的殘數(shù)相反,將兩殘余點(diǎn)用支切線相連,如果搜索到的殘余點(diǎn)未被其他支切線連接,則求支切線的殘數(shù),即支切線上各殘余點(diǎn)的殘數(shù)和;若殘數(shù)為0,整個(gè)支切線上正負(fù)殘余點(diǎn)相互抵消,則重新搜索干涉圖中未被支切線連接的殘余點(diǎn)。若在窗口內(nèi)搜索到干涉圖邊界,則用支切線連接搜索中心點(diǎn)與干涉圖邊界,將該支切線殘數(shù)置零,重新搜索干涉圖中未被支切線連接的殘余點(diǎn);若未搜索到殘余點(diǎn)或邊界,或搜索到殘余點(diǎn)后支切線上的殘數(shù)不為零,則以支切線上其余殘余點(diǎn)為中心搜索,直至搜索完整個(gè)支切線上所有殘余點(diǎn),則將搜索窗口大小增加2,繼續(xù)搜索。
6.一種針對(duì)干涉合成孔徑聲納信號(hào)的相位處理系統(tǒng),用于處理干涉合成孔徑聲納相位噪聲嚴(yán)重的問(wèn)題,所述系統(tǒng)包含: 配準(zhǔn)模塊,用于經(jīng)過(guò)相位差波動(dòng)函數(shù)配準(zhǔn)方法,根據(jù)主接收陣和輔接收陣圖像之間的相位差別,校正主輔圖像之間的相位偏移,完成主接收陣圖像和輔接收陣圖像對(duì)的配準(zhǔn);干涉相位圖獲取模塊,用于對(duì)完成配準(zhǔn)的主接收陣圖像和輔接收陣圖像對(duì)進(jìn)行相干處理,得到干涉相位圖; 濾波模塊,用于將干涉相位轉(zhuǎn)變?yōu)閺?fù)單位矢量,再依據(jù)圓周期中值濾波對(duì)干涉相位圖進(jìn)行濾波; 纏繞相位獲取模塊,用于基于濾波后的干涉相位值,使用最小二乘解纏,得到干涉相位,該相位取值位于(-00,+ 00 )之間; 建模模塊,用于根據(jù)纏繞相位和InSAS模型進(jìn)行高程重建,完成對(duì)干涉合成孔徑聲納信號(hào)的相位處理。
7.根據(jù)權(quán)利要求6所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理系統(tǒng),其特征在于,所述系統(tǒng)還包含: 主輔圖像獲取模塊,用于使用一副發(fā)射陣列、兩副平行放置的SAS陣列;所述發(fā)射陣以一定的脈沖重復(fù)周期發(fā)射信號(hào),到達(dá)目標(biāo)場(chǎng)景返回,兩副接收陣以相近的視角接收回波信號(hào)。分別經(jīng)過(guò)頻域成像算法,得到目標(biāo)場(chǎng)景的主輔圖像對(duì)。
8.根據(jù)權(quán)利要求6所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理系統(tǒng),其特征在于,所述系統(tǒng)還包含: 補(bǔ)償模塊,用于使用分支割線方法補(bǔ)`償由于相位的周期特性引起的干涉相位模糊的步驟。
9.根據(jù)權(quán)利要求6所述的針對(duì)干涉合成孔徑聲納信號(hào)的相位處理系統(tǒng),其特征在于,所述干涉相位圖獲取模塊進(jìn)一步包含: 復(fù)干涉圖獲取模塊,用于基于配準(zhǔn)后的主接收陣圖像和輔接收陣圖像得到復(fù)干涉圖; 干涉相位圖獲取模塊,用于對(duì)復(fù)干涉圖求取相位得到干涉相位圖。
【文檔編號(hào)】G01S7/52GK103487792SQ201210261997
【公開(kāi)日】2014年1月1日 申請(qǐng)日期:2012年7月26日 優(yōu)先權(quán)日:2012年6月12日
【發(fā)明者】秦曉東, 劉紀(jì)元 申請(qǐng)人:中國(guó)科學(xué)院聲學(xué)研究所