王曉春 楊軍 計建軍 周盛

0 引言

眼科超聲成像技術(shù)目前是臨床上廣為應(yīng)用的一種診斷方式，相對于傳統(tǒng)應(yīng)用的10 MHz頻率，20 MHz及以上頻段超聲波對于顯示視網(wǎng)膜[1]、脈絡(luò)膜[2]、黃斑病變[3]等眼球壁細微結(jié)構(gòu)，以及白內(nèi)障晶體[4]成像具有更好圖像分辨效果。但由于其頻率較高，衰減速度快，對于深部眼眶病變組織，包括視神經(jīng)、肌肉、脂肪及相應(yīng)腫瘤的成像受到限制[5-6]。而專用眼科超聲影像設(shè)備對于時間平均聲強、脈沖平均聲強和機械指數(shù)等聲輸出參數(shù)的要求更為嚴格[7]，不能單純采用提高瞬間發(fā)射能量的方式。因此尋找一種安全、高效的方式，在保證較高分辨能力的前提下，提高20 MHz眼科超聲圖像的信噪比和探測深度，具有重要的臨床意義。

超聲編碼激勵技術(shù)作為提高聲波穿透力、獲取微弱信息、改善圖像質(zhì)量的關(guān)鍵手段，其安全有效性已獲得普遍認可，在醫(yī)學超聲成像與測量中具備很好的研究潛力與應(yīng)用前景。例如：在旋轉(zhuǎn)式血管內(nèi)超聲成像中[8]，可獲得滿意的圖像質(zhì)量和旁瓣噪聲抑制比，有效改善血管壁的層次信息,突出冠狀動脈血管層；內(nèi)窺鏡超聲成像中[9]，可在提高圖像分辨率和信噪比的同時，限制超聲波的平均功率和峰值功率，減少超聲熱效應(yīng)和空化效應(yīng)給患者帶來的傷害；超聲瞬時彈性成像技術(shù)中[10]，可提高剪切波的信噪比和探測深度，尤其適用于肥胖肝硬化患者的診斷；在評價骨質(zhì)疏松程度與骨折深度的研究中[11]，可增強背向散射信號的幅值，提高抗噪聲性能，節(jié)約硬件系統(tǒng)成本等。

本研究基于實驗室自行設(shè)計的數(shù)字化眼科超聲成像系統(tǒng)[12]。以本系統(tǒng)為核心，課題組與天津邁達醫(yī)學科技股份有限公司共同合作完成了MD-2300S眼科A/B型超聲診斷儀的設(shè)計與研發(fā)工作，該儀器已廣泛應(yīng)用于眼科臨床。為彌補眼科臨床應(yīng)用中20 MHz超聲波成像深度不足的缺陷，通過采用8位Golay互補序列編碼激發(fā)超聲換能器，回波信息經(jīng)高速采集和匹配濾波處理后，使用正、反編碼掃描線復用的方法，完成解碼壓縮運算。在確保顯示分辨力和幾何位置精度一致，且不降低每幅圖像的掃描線數(shù)和最終成像幀頻的前提下，提高圖像信噪比和探測深度，從而拓展20 MHz超聲波成像在眼科臨床中的應(yīng)用范圍。

1 方法

1.1 系統(tǒng)方案設(shè)計

系統(tǒng)硬件電路以大規(guī)模現(xiàn)場可編程門陣列(field programmable gate array，F(xiàn)PGA)及USB 3.0接口電路設(shè)計為核心，并實現(xiàn)與上位機的通訊。輸出編碼脈沖時序序列通過電平轉(zhuǎn)換后激勵換能器發(fā)射20 MHz高頻超聲波。輸出馬達驅(qū)動時序信號控制步進電機轉(zhuǎn)動，從而帶動換能器實現(xiàn)扇形掃描?；夭ㄐ畔⒔?jīng)前置放大與模數(shù)變換后完成數(shù)字化，其增益倍數(shù)和采樣參數(shù)可由上位機通過FPGA進行設(shè)定。硬件系統(tǒng)結(jié)構(gòu)如圖1所示。

圖1 硬件系統(tǒng)結(jié)構(gòu)圖Figure 1 Diagram of the hardware system

1.2 編、解碼掃描方式

Golay互補編碼序列是一種非常理想的編碼方式，其沒有碼長限制，最大優(yōu)勢在于可完全消除旁瓣達到最佳脈沖壓縮效果[13]。但由于Golay互補編碼序列雙次發(fā)射的機制，會使每幅圖形掃描線線數(shù)減少一半。如圖2所示，選用8位Golay互補序列，激勵超聲換能器產(chǎn)生超聲波。通過一對Golay正序列編碼X1、X2、X3、…、Xn和Golay反序列編碼Y1、Y2、Y3、…、Yn依次交替激發(fā)單陣元換能器，每一幀圖像依次發(fā)射2n次超聲波信號，形成2n條掃描線。如按照傳統(tǒng)的解碼方式處理，解碼后每幀圖像只剩余n條掃描線進行重建，這將嚴重影響最終的成像效果。因此本文設(shè)計中，在每次發(fā)射結(jié)束后，通過換能器接收每條掃描線上的回波，采用相鄰正、反編碼掃描線數(shù)據(jù)復用的方法交替相加，完成解碼運算，形成2n-1條回波掃描線Z1～Z2n-1，最終得到扇形的眼部超聲波成像，疊加示意圖見圖3。

圖2 掃描線示意圖Figure 2 Schematic diagram of scanning lines

圖3 相鄰掃描線回波疊加的示意圖Figure 3 Schematic diagram of echo superposition of adjacent scanning lines

1.3 FPGA的設(shè)計與實現(xiàn)

實驗中采用大規(guī)?，F(xiàn)場可編程門陣列FPGA，實現(xiàn)數(shù)據(jù)的采集、存儲與處理。圖4為FPGA內(nèi)部信號流程圖，由AD數(shù)字采樣后的回波信號先進入FPGA進行緩存，信號流隨后進入正反碼匹配濾波模塊進行濾波，由選擇器分別在奇數(shù)線和偶數(shù)線時選擇相應(yīng)的濾波結(jié)果。然后每條線的濾波結(jié)果由讀寫使能信號控制4個FIFO緩存的讀寫時序。從緩存讀出的結(jié)果兩兩相加，再由選擇器根據(jù)不同時序選擇輸出對應(yīng)的每條掃描線的結(jié)果。

圖4 基于FPGA 的回波信號處理流程圖Figure 4 Flowchart of echo signal processing based on FPGA

具體FIFO的讀寫時序如圖5所示。緩存1～4采用FIFO先進先出數(shù)據(jù)存儲器。具體存儲過程如圖5所示：第i線表示奇數(shù)線，i=0,1,2,…,n；第j線表示偶數(shù)線，j=0,1,2,…,n； 其中i=j。

圖5 FIFO讀寫時序設(shè)計Figure 5 Read and write timing design of FIFO

2 初步驗證及結(jié)果

為了驗證Golay互補編碼序列激勵在20 MHz頻段超聲成像的應(yīng)用效果，實驗中設(shè)計了編碼發(fā)射激勵實驗、分辨率檢測實驗和仿組織超聲體模實驗。

2.1 編碼發(fā)射激勵波形

本文基于Altera公司的Quartus II 13.0軟件開發(fā)環(huán)境和Verilog硬件描述語言實現(xiàn)了Golay互補序列的激勵。圖6為8位Golay正、反序列的FPGA激勵時序圖，其中ep是線同步信號，pa、pb、pc、pd和na、nb、nc、nd為正反序列的四路電平觸發(fā)信號。圖7為由Tektrnoix公司的MSO 4104示波器測得的激勵換能器的脈沖波形。其中，發(fā)射脈寬為50 ns，脈沖幅值為±90 V。實驗結(jié)果證明了由FPGA產(chǎn)生的時序信號以及激勵超聲換能器的觸發(fā)波形符合本實驗設(shè)計要求。

圖6 8位Golay正、反序列激勵時序圖Figure 6 Excitation sequence diagram of the 8-bit Golay positive and negative sequences

圖7 Golay正、反序列激勵脈沖Figure 7 Excitation pulses of Golay positive and negative sequences

2.2 分辨率測試實驗

根據(jù)國家醫(yī)療器械行業(yè)標準《YY 0773—2010眼科B型超聲診斷儀通用技術(shù)條件》的要求，使用分辨力測試線靶模塊對常規(guī)單脈沖激勵與8位Golay互補編碼序列激勵所獲得圖像的軸向分辨力和側(cè)向分辨力進行對比。將測試線靶置于底部敷有吸聲材料的水槽中，探頭置于測試靶上方，平行靶線排列方向分別與探頭端面保持垂直或平行。探頭掃查平面保持在與靶線垂直的平面上進行掃描。增益、亮度、對比度等調(diào)節(jié)機構(gòu)置于最佳位置，調(diào)節(jié)探頭與靶線的距離至最佳值，微動探頭角度，從顯示器上觀察圖像，直至平行靶線能被區(qū)分開時，則判定分辨力達到了靶線間距離值。如圖8所示，靶線采用10 μm±1 μm鎢絲，使用光學顯微鏡制作線靶模塊并驗證靶線間隔。目前實驗室自行設(shè)計的20 MHz單脈沖數(shù)字化眼科超聲成像系統(tǒng)的軸向分辨率和側(cè)向分辨率分別為80 μm和150 μm。因此在這個編碼激勵的實驗中，測量軸向分辨力時，間隔設(shè)置為80 μm；測量側(cè)向分辨力時，間隔設(shè)置為150 μm。圖像掃描結(jié)果如圖9所示，Golay互補編碼模式可達到與單脈沖一致的分辨能力。

圖8 分辨率測試裝置示意圖Figure 8 Schematic diagrams of resolution testing device

圖9 分辨率測試結(jié)果Figure 9 Resolution test results

2.3 仿組織超聲體模實驗

為驗證編碼激勵的方式對探測深度的提升，本實驗中分別采用單脈沖激勵和8位Golay編碼激勵的方式對仿組織超聲體模(中國科學院聲學研究所制造的KS107BG)進行圖像采集。該體模適用于工作頻率≥5 MHz超聲設(shè)備，其填充的特殊材料具有與人體軟組織相似的超聲特性，縱向靶群相鄰兩線中心距離均為10 mm。

對比結(jié)果如圖10所示，當系統(tǒng)增益倍數(shù)設(shè)置為相同時，采用編碼激勵后的圖像探測深度有顯著提升，圖10(a)中的3 cm靶點回波信號只是隱約可見，小信號探測深度也只能達到2.5 cm，而圖10(b)中3 cm靶點已清晰可見，小信號探測深度也接近3 cm。整體圖像的信噪比也有了顯著改善。

圖10 體模成像對比圖Figure 10 Phantom imaging results

3 討論與結(jié)論

為平衡醫(yī)學高頻超聲成像中縱向分辨率和探測深度這兩項指標之間的矛盾，本研究以數(shù)字化眼科高頻超聲成像系統(tǒng)為平臺，設(shè)計采用Golay互補序列編碼激勵，并通過正、反編碼掃描線復用的方法，完成解碼壓縮運算。在保證每幅圖像掃描線數(shù)的前提下，使得成像幀頻與單脈沖激勵成像方式保持一致。實驗中由FPGA實現(xiàn)了該方法，并同步傳輸至上位機進行實時成像。通過掃查鎢絲靶線模塊和仿組織超聲體模實驗，與傳統(tǒng)的單脈沖激勵方式相比，證明了編碼激勵的方法可在保持成像分辨力不變的前提下，顯著提高信噪比和探測深度，從而可拓展20 MHz眼科超聲成像在臨床中的應(yīng)用范圍。

編碼激勵方法很多，本實驗中選用Golay編碼序列也是考慮了軟、硬件系統(tǒng)復雜度、實時處理能力與眼睛屬于基本靜止器官的緣由。但其最大的缺陷在于需要發(fā)射2次才可以完成脈沖壓縮，雖然在本實驗方案中采用了正、反掃描線復用的方法，但是當掃查線較為稀疏時，也會由于聲束位置偏移而產(chǎn)生誤差。因此，在下一步的深入研究中，將采用更好的碼型，如M序列[14]、MAC碼[15]、chirp碼[16]等其他相位編碼模式。同時，可以引入新的脈沖壓縮方法，如失配濾波器[17]、維納濾波[18]等，并進一步研究和改進實現(xiàn)方案，以提高信噪比，提升成像質(zhì)量。