前言：尋找寫作靈感？中文期刊網用心挑選的Barker碼在眼科的應用研究，希望能為您的閱讀和創作帶來靈感，歡迎大家閱讀并分享。

作者：周盛 王曉春 計建軍 楊軍 王延群 單位：中國醫學科學院生物醫學工程研究所

1方法

1．1發射編碼

發射編碼按其內容分為相位編碼和頻率編碼，也可以按照發射次數分為單次發射碼和多次發射碼，如表1所示。Chirp碼脈沖壓縮效果受衰減和波束合成等非線性因素的影響相對較小，而且只需要一次發射就可以完成脈沖壓縮，受組織運動影響較小，因此是超聲成像系統編碼信號的理想選擇之一。但是由于Chirp碼需要多幅值發射電壓，對發射電壓和硬件的要求較高。因此在數字化高頻眼科超聲診斷設備中，基于發射編碼、系統復雜程度等考慮，一般采用二相編碼信號作為發射編碼。Golay互補序列雖然可以完全消除旁瓣達到最佳脈沖壓縮效果，但兩次發射的機制，容易受到組織運動的影響，引入誤差。雖然眼睛是基本靜止的生物組織，但由于我們所選用的是機械扇形掃描探頭，所以還是會在一定程度上影響成像效果。M序列在碼長較短時的脈沖壓縮效果較barker碼差;又根據文獻［10－11］給出的最優編碼序列的概念、判斷標準和計算方法，通過窮舉搜索可以得到不同碼長的最優編碼序列［8］，Barker碼除了碼長N=11時不是最優編碼序列，其余都是。而且文獻［12］中也對13位Barker碼應用于高頻超聲成像中做了一定的實驗室研究。綜合比較，最終采用13位Barker碼作為本項目數字超聲眼科診斷系統中編碼激勵的發射碼型。

1．2解碼壓縮方式

目前醫學超聲編碼激勵系統中常用的脈沖壓縮方法分為匹配濾波和非匹配濾波兩種。對于Golay碼和Chirp碼，匹配濾波都能取得很好的脈沖壓縮效果，并且設計及實現方法相對簡單。而對于Barker碼，采用匹配濾波脈沖壓縮方法具有最低的距離旁瓣，但主瓣幅度也只是旁瓣幅度的N倍，其中N是Barker碼的長度。對于最長的13位Barker碼，匹配脈沖壓縮得到的信噪比增益為11．1dB，而峰值旁瓣水平(peaksidelobelevel，PSL)也只有(22dB，對距離旁瓣的抑制不夠好［9］。為此設計出了逆濾波器［13－14］、維納濾波器［15］、尖峰濾波器［16］以及失配濾波器等非匹配濾波脈沖壓縮方法［12，17］，以損失很小的信噪比增益為代價來換取更好的距離旁瓣脈沖壓縮效果。對比各種解碼壓縮方式的優缺點，本研究對匹配濾波與非匹配濾波分別進行了仿真實驗，結果表明。失配濾波方法可以達到最好的脈沖壓縮效果。通過此方法，可以使平均距離旁瓣水平(integratedsidelobelevel，ISL)和PSL都盡可能的低。

1．3仿真實驗

二進制13位Barker碼表示為:［+1，+1，+1，+1，+1，－1，－1，+1，+1，－1，+1，－1，+1］，由于采樣頻率一般選用信號回波頻率的8倍，這里用以下序列來構造每周期8個采樣點的正弦波，分別用［0，－7，－10，－7，0，7，10，7］表示－1，用［0，7，10，7，0，－7，－10，－7］表示+1，圖1～圖3為Barker碼的Matlab仿真波形。其中，圖1為13位Barker碼波形。圖2為經過匹配濾波后得到的結果，主瓣峰峰值達到了±312，而旁瓣峰峰值為±24，主副比為13。圖3為13位Barker碼匹配濾波結果的峰值包絡波形，可見旁瓣衰減達到了－45dB。13位Barker碼匹配濾波的FPGA實現相對簡單，需要開辟一個深度13×8，14bit的存儲空間進行卷積解碼。將解碼濾波系數序列中為+1的系數對應的回波數據保持原數，將解碼濾波系數序列中為(1的系數對應的回波數據進行補碼取反，然后進行對應數據并行累加，即可得到匹配濾波的結果。圖4是用Modelsim仿真得到的結果圖，上端波形為13位Barker碼編碼激勵信號，下端波形為經過匹配濾波解碼壓縮后的信號?？梢?，Barker碼匹配濾波解碼壓縮后存在較大的旁瓣。圖5為13位Barker碼經過失配濾波后得到的結果，實驗中我們設計的失配濾波器的階數為P=3N=39。主瓣峰峰值達到了±127，而旁瓣峰峰值為±1，主副比為127，是匹配濾波主副比的9．7倍。由此說明失配濾波雖然降低了主瓣峰值，但有效地降低了旁瓣水平，提高了信噪比。圖6為13位Barker碼失配濾波結果的峰值包絡波形，可見旁瓣衰減達到了－60dB。在FPGA的設計中，具體的實現步驟如圖7所示。首先將失配濾波器系數進行量化，乘以27=28，并取整存入FPGA的RAM中，然后構建實時延時回波矩陣，將系數分別與回波矩陣各行相乘，隨后將矩陣39行相乘的結果累加，并除以量化倍數，即可得到失配濾波解碼壓縮的結果。圖8是用Modelsim仿真的結果圖，上端波形為13位Barker碼編碼激勵信號，下端波形為經過失配濾波解碼壓縮后的信號?？梢姡?3位Barker碼失配濾波解碼壓縮后有效地壓縮了旁瓣，在保證分辨率的前提下可以有效地提升信噪比。

1．4系統方案與設計

編碼激勵系統和傳統的脈沖回波成像系統不同之處在于:(1)發射電路采用編碼發射激勵，激勵脈沖的長度遠大于單脈沖;(2)接收電路需要對回波信號進行脈沖壓縮，脈沖壓縮算法的計算量很大。因此，編碼發射的實驗系統設計有以下特殊要求:(1)發射電路能夠產生各種不同的編碼激勵;(2)該系統成像幀頻是10幀/s，每幀由512條線組成，每條線包含768個點，要求FPGA能夠實時地對編碼激勵的回波進行脈沖壓縮并且成像。根據以上系統設計要求，圖9為系統的整體框架結構圖。主要包括上位機主機、硬件電路板、3個超聲探頭、LCD液晶顯示、腳閘、USB外設等。系統中的上位機選用嵌入式的工控主板，這樣可以實現儀器的小型化和便攜性。上位機主要用來實現對底層數據的接收和處理、以及軟件的實時控制。實驗系統采用單陣元機械扇形掃描探頭，3個超聲探頭分別為中心頻率為10MHz的B超探頭、中心頻率為20MHz的B超探頭和中心頻率為10MHz的A超探頭，其分別用來作為眼球全景成像、眼部組織高頻成像和眼球生物測量等模式的檢測探頭。LCD液晶顯示用來對圖像實時成像。雙鍵腳閘實現探頭的啟動、凍結。USB外設包括存儲器、打印機、鼠標和鍵盤。USB外設的使用，使得整機更加人性化和提高了可操作性。

硬件電路板是整個系統的核心部件，其結構如圖10所示。整個硬件電路以FPGA及其外圍電路為核心進行設計。FPGA采用Altera公司的CycloneIII系列EP3C55型產品。機械扇掃探頭的偏轉控制、圖像顯示的同步、USB接口邏輯等系統控制信號均由FPGA產生。FPGA產生的4組編碼脈沖序列通過MOSFET觸發電路將脈沖幅度由3.3V轉化為12V，以激勵發射電路，如圖11所示。發射電路由兩組雙極性脈沖發射芯片構成，并帶有歸零功能。編碼序列通過發射電路后，成為雙極性脈沖序列，激勵超聲換能器晶片?；夭ㄐ盘柦浉邏焊綦x電路后進入前置放大與可變增益控制電路。前置放大電路增益范圍為10～15dB，FPGA產生TGC曲線數據，控制可變增益放大器，可以使增益變化50dB左右。實驗系統時鐘為80MHz，高速A/D芯片按照不同模式工作在80MHz和120MHz，采樣精度為14bits，其采樣模式由FPGA通過SPI接口完成設定。FPGA完成的超聲回波數字信號處理包括:高頻數據采集與存儲、數字濾波、數字檢波、對數變換、二次采樣等，如圖12所示。數字信號處理之后的回波數據通過SRAM緩存后，經USB2．0接口電路上傳至計算機，并由人機界面完成整個系統的操作與控制。#p#分頁標題#e#

2結果

為了對比13位Barker碼脈沖壓縮中的匹配濾波和失配濾波在實際應用中的效果，設計了單反射面回波實驗。編碼激勵信號由FPGA產生，經發射電路激勵10MHz超聲換能器，換能器在編碼信號激勵下，發射超聲波，超聲信號經反射面反射驅動換能器產生回波信號?；夭ㄐ盘柦浄糯蠛笸ㄟ^A/D采樣進入FPGA進行數據處理，脈沖壓縮解碼，并最終通過D/A輸出到示波器上顯示。圖13顯示了發射電壓為±80V時，13位Barker碼的編碼激勵波形。圖14是編碼激勵的單反射面回波信號。圖15和圖16顯示了13位Barker碼的編碼激勵回波經過匹配濾波和失配濾波后的脈沖壓縮回波。

由于該系統中回波信號經A/D采樣后為14位輸入，經匹配濾波后輸出為18位，經失配濾波后輸出為22位，而系統所選用的D/A只有10位，圖15和圖16是示波器實時凍結的回波圖像，只能截取適當位數通過DA輸出至示波器顯示比對，無法如仿真軟件做到精確量化比對值。從圖15中可見，13位Barker碼激勵回波經過匹配濾波解碼后單反射面回波的峰峰值為0．6V，而且存在一些距離旁瓣，最大旁瓣峰峰值達到了0．15V，在超聲回波信號中，強回波旁瓣幅值過高會淹沒弱小回波信號，影響信號的動態范圍，因此抑制脈壓旁瓣有著重要的現實意義。而圖16中13位Barker碼激勵回波經過失配濾波解碼后單反射面回波就對旁瓣做了有效的抑制，回波峰峰值為0．72V，而最大距離旁瓣為0．06V，增大了主副比，并且距離旁瓣隨著離主瓣距離的增大而明顯減小。因此在最終的B超成像中，所選用的是失配濾波的脈沖壓縮算法。圖17為10MHz探頭眼科B超的眼球全景圖，其中(a)圖為單脈沖發射接收的眼球B超圖;(b)圖為13位Barker碼編碼發射且用失配濾波的脈沖壓縮方法實現的眼球B超圖。對比兩圖的二維眼球圖像以及下端的A超波形，且實驗中用分辨力測試線靶進行檢測，采用編碼發射和解碼壓縮之后的眼科B超圖在保持縱向分辨率為0．2mm和基底噪聲的同時，大大提升了信號的探測深度，提高了圖像的信噪比。

3討論和結論

本研究搭建了數字化眼科高頻超聲成像實驗平臺，基于Barker碼實現了編碼激勵和解碼壓縮技術，最終將實驗所得眼科超聲圖像數據，通過USB接口實時的在上位機上處理和顯示，并獲得了定量的實驗數據以及成像結果。實驗結果表明，編碼激勵的回波信噪比明顯優于單脈沖激勵。采用編碼激勵成像的圖像分辨率與單脈沖激勵系統相當，且圖像沒有出現明顯的距離旁瓣，而圖像的探測深度和噪聲水平明顯優于單脈沖激勵的圖像。編碼激勵系統還對微弱信號有明顯的提升作用，由此可以檢測出眼球玻璃體混濁等病變，在眼科診斷中具有重要意義。

本設計初步在系統平臺上實現了編碼激勵技術。但還有較多有待拓展、深入之處，主要包括:對于不同帶寬的超聲換能器，應對發射編碼進行相應的加窗處理，使得加窗后形成的激勵信號頻帶與換能器的頻帶接近，此時可充分發揮超聲換能器的轉換能力;在之后的研究中，將應用基礎序列調制技術來提高長編碼發射信號的質量;在脈沖壓縮仿真過程中，強發射界面還存在一定的旁瓣影響，在后續的研究中將借鑒雷達信號處理中的線性規劃法、最小二乘法以及迭代算法等旁瓣抑制方法，以達到更好的脈沖壓縮效果;為了進一步降低脈沖平均聲強(ISPPA．3)和機械指數(MI)等聲輸出參數，下一步將降低發射電壓，并觀察對成像質量的影響，以便選擇更合適的編碼及解碼方式。