一種帶運動補償?shù)暮铣煽讖匠暢上穹椒?br>【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及合成孔徑成像領(lǐng)域，特別涉及一種帶運動補償?shù)暮铣煽讖匠暢上穹?法。
【背景技術(shù)】
[0002] 合成發(fā)射孔徑通過單陣元發(fā)射，全孔徑接收，經(jīng)波束形成得到低分辨率圖像 (LRI)，再通過將所有發(fā)射陣元得到的LRI相干相加得到高質(zhì)量圖像（HRI)。由于HRI對發(fā) 射和接收都實現(xiàn)了動態(tài)聚焦，所以其較傳統(tǒng)的超聲成像方法就有更高的分辨率。但是，這種 合成發(fā)射孔徑成像方法面臨三個主要問題：1)單陣元發(fā)射造成的穿透性差和回波信號信 噪比低；2)低分辨率圖像間的運動會降低其相干性，從而影響成像質(zhì)量；3)系統(tǒng)實現(xiàn)過程 中需要存儲和傳輸大量的射頻回波（RF)線，系統(tǒng)復(fù)雜度很高。
[0003] 對于合成發(fā)射孔徑成像方法的上述三個問題，目前有多種改進方法。合成孔徑序 列波束方法就是一種比較容易實現(xiàn)的方法。合成孔徑序列波束方法由于不需要存儲和傳輸 大量的單陣元接收的RF線數(shù)據(jù)，其實現(xiàn)復(fù)雜度較合成發(fā)射孔徑方法大大降低。此外，由于 是多陣元發(fā)射，所以也不存在合成發(fā)射孔徑方法單陣元發(fā)射造成的穿透性差和回波信號信 噪比低的問題。但是，合成孔徑序列波束方法作為一種合成孔徑方法，其中仍然存在運動影 響成像質(zhì)量的問題。
[0004] 傳統(tǒng)的運動補償算法需要增加多次發(fā)射次數(shù)，然后利用增加的接收回波進行運動 估計，再進行補償。這樣的方法的缺點是降低了系統(tǒng)成像的掃描幀頻，會造成圖像不連貫， 高速運動物體成像不佳等問題。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的在于解決現(xiàn)有技術(shù)在存在軸向位移時成像質(zhì)量較差的問題，從而提 供一種無需增加發(fā)射次數(shù)卻能得到高分比率圖像的方法.
[0006] 為了實現(xiàn)上述目的，本發(fā)明提供了一種帶運動補償?shù)暮铣煽讖匠暢上穹椒?，?yīng) 用于帶有N個陣元的合成孔徑超聲換能器陣列中，各個陣元之間等間距，且各個陣元均能 發(fā)射信號與接收信號；該方法包括：
[0007] 步驟1)、采用發(fā)射、接收都是單焦點，且發(fā)射焦點與接收焦點相同的線性掃描方式 得到第一低分辨率掃描線；采用偏移陣元接收的方式得到運動估計掃描線；
[0008] 步驟2)、利用步驟1)得到的運動估計掃描線進行估計，所得到的估計值用于對步 驟1)所得到的第一低分辨率掃描線做補償，從而得到第二低分辨率掃描線；
[0009] 步驟3)、將虛擬源開角內(nèi)其它第二低分辨率掃描線上的相應(yīng)采樣點乘以加權(quán)值后 與虛擬源對應(yīng)第二低分辨率掃描線上的相應(yīng)采樣點相加，得到高分辨率掃描線，進而得到 高分辨率圖像；其中，所述虛擬源為合成孔徑超聲換能器陣列中各個陣元的發(fā)射焦點。
[0010] 上述技術(shù)方案中，所述步驟1)包括：
[0011] 采用合成孔徑超聲換能器陣列中的第1到第Μ個陣元發(fā)射信號，然后采用第i到 第Μ+i-l個陣元發(fā)射信號，以此類推，直到第N-M+1到第N個陣元發(fā)射信號，共N-M+1次發(fā) 射信號，其中，將第i次發(fā)射的階段記為!\，1 <i<N-M+1 ;在接收信號時，則會依據(jù)陣列中 發(fā)射信號的陣元的位置來選取合適的接收陣元；具體包括：
[0012] 當(dāng)i= 1，即第1\次發(fā)射時，用所述陣列中的第1到第Μ個陣元發(fā)射信號，發(fā)射焦 點為VSi;由所述陣列中的第1到第Μ+1個陣元接收信號，每個陣元的采樣點數(shù)為Ρ，陣列中 的陣元在接收信號時需經(jīng)歷正常接收與右移接收兩個階段；
[0013] 正常接收：對第1到第Μ陣元所接收的信號進行波束形成，接收焦點也為VSi，得到 第1條第一低分辨率掃描線Li，該掃描線有P個采樣點；
[0014] 右移接收：對第2到第M+1陣元所接收的信號進行波束形成，接收焦點也為VRi， VRi = VS^ △，得到第1條右移焦點掃描線LR1，該掃描線有P個采樣點；
[0015] 當(dāng)l〈i〈N_M+l時，即第?\次發(fā)射時，用所述陣列中的第i到第Μ+i-l個陣元發(fā)射 信號，發(fā)射焦點為VSi，由所述陣列中的第i-Ι到第M+i個陣元接收信號，每個陣元的采樣點 數(shù)為P，陣列中的陣元在接收信號時需經(jīng)歷左移接收、正常接收與右移接收三個階段：
[0016] 左移接收：對第i-Ι到第M+i-2陣元所接收的信號進行波束形成，接收焦點也為 VQ，VQ=VS「Λ，得到第i-Ι條左移焦點掃描線Lu，該掃描線有P個采樣點；
[0017] 正常接收：對第i到第Μ+i-l陣元所接收的信號進行波束形成，接收焦點也為VSp 得到第i條第一低分辨率掃描線Q，該掃描線有P個采樣點；
[0018] 右移接收：對第i+Ι到第M+i陣元所接收的信號進行波束形成，接收焦點也為VRp VRi=VSJ△，得到第i條右移焦點掃描線LRl，該掃描線有P個采樣點；
[0019] 當(dāng)i=N-M+1時，即第TNM+1次發(fā)射時，用所述陣列中的第N-M+1到第N個陣元發(fā) 射信號，發(fā)射焦點為VSNM+1，由所述陣列中的第N-M到第N個陣元接收信號，每個陣元的采樣 點數(shù)為P，陣列中的陣元在接收信號時需經(jīng)歷左移接收與正常接收兩個階段；
[0020] 左移接收：對第N-Μ到第N-ι陣元所接收的信號進行波束形成，接收焦點也為 VLNM+1，VLNM+1 =VSNΜ+「Λ，得到第N-Μ條左移焦點掃描線LUNM+1)，該掃描線有P個采樣點；
[0021] 正常接收：對第N-M+1到第N陣元所接收的信號進行波束形成，接收焦點也為 VSNM+1，得到第N-M+1條第一低分辨率掃描線LNM+1，該掃描線有P個采樣點；
[0022] 依次重復(fù)上述步驟，直至得到N-M+1條第一低分辨率掃描線；其中，第一低分辨率 掃描線LpL2,……，LNM+1構(gòu)成一幅低分辨率圖像LRI;所得到的左移焦點聚焦掃描線U2， LL3,……，LUnM+1)用于與右移焦點聚焦掃描線LR1，LR2，……，L_M)做運動估計處理。
[0023] 上述技術(shù)方案中，所述步驟2)進一步包括：
[0024] 步驟2-1)、計算第一條左移焦點聚焦掃描線與第一條右移焦點聚焦掃描線LRi 之間的時延，得到第一時延值Di;
[0025] 步驟2-2)、計算第i條左移焦點聚焦掃描線LUl+1)與第i條右移焦點聚焦掃描線 LRW之間的時延，得到第i個時延值Di;
[0026] 步驟2-3)、第一條第一低分辨率掃描線Q不做補償，直接作為第一條第二低分辨 率掃描線，記作CQ;
[0027] 步驟2-4)、對于第二條第一低分辨率掃描線L2，利用第一時延值01進行延時處理， 得到補償后的第二條第二低分辨率掃描線CL2;
[0028] 步驟2-5)、利用第1至第i時延值之和D1+D2+~+D1對第i+1條第一低分辨率掃 描線L1+1進行延時處理，得到補償后的第i+1條第二低分辨率掃描線CLn+1，直至i的值達(dá)到N-M，得到補償后的第N-M+1條第二低分辨率掃描線CLNM+1 ;
[0029] 其中，CL。CL2……CLNM+1稱為第二低分辨率掃描線。
[0030] 上述技術(shù)方案中，所述步驟3)進一步包括：
[0031]步驟3-1)、對第一條第二低分辨率掃描線進行第二次波束形成，即將VSi開角內(nèi)其 它第二低分辨率掃描線上的相應(yīng)采樣點乘以加權(quán)值后與第一條第二低分辨率掃描線上相 應(yīng)采樣點相加，得到第一條高分辨率掃描線氏；
[0032]步驟3-2)、對第二條第二低分辨率掃描線進行第二次波束形成，即將￥&開角內(nèi)其 它第二低分辨率掃描線上的相應(yīng)采樣點乘以加權(quán)值后與第二條第二低分辨率掃描線上相 應(yīng)采樣點相加，得到第二條高分辨率掃描線H2;
[0033]步驟3-3)、重復(fù)上述處理，直到對第N-M+1條第二低分辨率掃描線進行第二次波 束形成，即將VSNM+1開角內(nèi)其它第二低分辨率掃描線上的相應(yīng)采樣點乘以加權(quán)值后與第 N-M+1條第二低分辨率掃描線上相應(yīng)采樣點相加，得到第N-M+1條高分辨率掃描線HNM+1。
[0034] 上述技術(shù)方案中，將虛擬源開角內(nèi)其它第二低分辨率掃描線上的相應(yīng)采樣點乘以 加權(quán)值后與虛擬源對應(yīng)第二低分辨率掃描線上相應(yīng)采樣點相加，得到高分辨率掃描線的具 體計算公式為 ：
[0036] 其中，CLXfc(td(fc))表示第二低分辨率掃描線CLXfc上延時td(k)對應(yīng)的值，W為加 權(quán)函數(shù)，其值隨深度Z和掃描線的水平位置xk而改變；K是深度z的函數(shù)，K(z)表示虛擬源 開角內(nèi)的低分辨率掃描線數(shù)量；
[0037]Κ(ζ)的計算公式如下：
[0039] 其