張耀楠，金應(yīng)東,，邱維寶

1 西安思源學院電子信息工程學院，西安市，710038

2 東北大學中荷生物醫(yī)學與信息工程學院，沈陽市，110169

3 中國科學院深圳先進技術(shù)研究院生物醫(yī)學與健康工程研究所，深圳市，518055

超聲平面波兩種數(shù)據(jù)壓縮方法初步比較

【作 者】張耀楠1,2，金應(yīng)東2,3，邱維寶3

1 西安思源學院電子信息工程學院，西安市，710038

2 東北大學中荷生物醫(yī)學與信息工程學院，沈陽市，110169

3 中國科學院深圳先進技術(shù)研究院生物醫(yī)學與健康工程研究所，深圳市，518055

超聲平面波成像是超聲成像研究的熱點之一, 超聲平面波帶來超高幀頻的同時也帶來了超大數(shù)據(jù)量。如果將這些數(shù)據(jù)量直接通過硬件傳輸，勢必對硬件帶來極大的負擔，也對硬件處理數(shù)據(jù)的能力提出了極高的要求。為了減少大數(shù)據(jù)量高幀頻所帶來的問題，該文通過數(shù)據(jù)壓縮的方式來降低采集數(shù)據(jù)的傳輸量。分別從有損壓縮和無損壓縮方向進行了研究，選擇了LZW壓縮算法來嘗試無損壓縮的效果，而根據(jù)平面波特點，提出頻域切割壓縮算法作為有損壓縮方法。將LZW算法在FGPA進行了實施，通過對超聲仿真數(shù)據(jù)進行測試，可以達到一定的壓縮比。在已知分數(shù)帶寬的條件下，可以求出頻譜中需要保留的帶寬，而其余部分頻率信號對圖像影響的意義不大。通過對超聲仿真數(shù)據(jù)的測試，頻域切割壓縮算法可以達到較高的壓縮比，而比較原始圖像和壓縮重建圖像，其差異是微小的。所以從壓縮比角度，頻域切割壓縮算法是進行超聲數(shù)據(jù)壓縮的優(yōu)先選擇。

超聲；平面波；數(shù)據(jù)壓縮；LZW

0 引言

超聲成像具有安全、方法簡單、價格便宜、能夠區(qū)分不同軟組織等優(yōu)點，使得超聲成像在醫(yī)院中得到廣泛應(yīng)用。傳統(tǒng)的超聲成像通常采用電子掃描的方式，進行多次聚焦發(fā)射來獲得一幅完整的圖像，所以成像幀頻較低。近年來，超聲平面波成像是研究的熱點之一[1-3]。這種技術(shù)是通過向介質(zhì)中發(fā)射一個寬度范圍很大的波束，然后記錄下反射的回波，最后運用數(shù)字平行接收波束合成[4]技術(shù)或者時間逆轉(zhuǎn)聚焦技術(shù)[5-6]來最終形成需要的超聲圖像，很大程度上提高了成像的幀頻。

平面波帶來超高幀頻優(yōu)勢的同時也帶來了超大數(shù)據(jù)量。如果將這些數(shù)據(jù)量直接通過硬件傳輸，這勢必對硬件帶來極大的負擔，也對硬件處理數(shù)據(jù)的能力提出了極高的要求。假設(shè)傳輸?shù)臄?shù)據(jù)一幀為4 000 ×128像素，而一幀的數(shù)據(jù)量就是4 000 ×128×2 byte。在超快成像領(lǐng)域一般采集速度可以達到上千幀，這里以1 000幀為例，則超聲采集設(shè)備每秒鐘理論需要傳輸?shù)臄?shù)據(jù)量為：4 000×128×2×1 000=1 GB，而當前USB3.0的數(shù)據(jù)傳輸速度最大能夠達到的極限速率為300 MB/s，所以，對于海量的平面波數(shù)據(jù)來說，采集和傳輸?shù)乃俾什怀杀壤@就給超聲實時成像造成了極大的瓶頸。為了減少大數(shù)據(jù)量高幀頻所帶來的問題，本文通過數(shù)據(jù)壓縮的方式來降低采集數(shù)據(jù)的傳輸量。目前國內(nèi)外關(guān)于超聲平面波數(shù)據(jù)壓縮的研究還較少，本文描述超聲平面波數(shù)據(jù)壓縮的兩種方法并加以比較。本文選擇了LZW(Lempel Ziv Welch)壓縮算法來嘗試無損壓縮的效果，而根據(jù)平面波特點首先提出頻域切割壓縮算法作為有損壓縮方法。

1 LZW壓縮算法及其硬件實現(xiàn)

1.1 LZW算法

LZW[7]算法有三個重要的對象，輸入數(shù)據(jù)流，輸出數(shù)據(jù)流和一張用于編碼的字符串表。輸入數(shù)據(jù)流指被壓縮的數(shù)據(jù)；輸出數(shù)據(jù)流指壓縮后輸出的代碼流；字符串表存儲的是數(shù)據(jù)的索引號，重復(fù)的數(shù)據(jù)流就是由同一索引號存儲，從而實現(xiàn)了數(shù)據(jù)的壓縮。

1.2 LZW算法硬件實現(xiàn)流程

LZW算法在FPGa[8]上具體的實現(xiàn)方式為：首先將FPGa初始化，設(shè)置兩個大小為4 096的字符串表，他們輪流工作，即當其中一個字符串表滿時就轉(zhuǎn)到另一個字符串表工作，同時對剛才的字符串表進行清空操作，為下一次填表做準備。當FPGa準備就緒后，在第一個時鐘脈沖上升沿到來時，向FIFo 1發(fā)出第一個讀信號，則這個通道輸送第一個待壓縮的數(shù)據(jù)到FPGa的查字符串模塊。經(jīng)過數(shù)據(jù)比較單元比較之后，判斷接收到的數(shù)據(jù)是否位于存儲器中，如果在，則把這個數(shù)據(jù)左移八位放到寄存器作為數(shù)據(jù)前綴，同時向FIFo 1再發(fā)一個讀信號，F(xiàn)IFo 1向FPGa輸送下一個字節(jié)的數(shù)據(jù)，F(xiàn)PGa將移位寄存器中的前綴與這個字節(jié)的數(shù)據(jù)相加后形成新的字符串，再判斷這個字符串是否位于字符串表中。如果不在字符串表中，則將這個新字符串添加到字符串表中，同時字符串表指針加一，然后將新的字符串表輸出到FIFo 2存儲器；如果在字符串表中，則將字符串表中的索引號作為前綴，繼續(xù)接受下一個字節(jié)。如此循環(huán)，直至第一個字符串表滿。當?shù)谝黄址頋M時，第一片的片選端置零，第二片選通，第二片字符串表開始工作。同時對第一片字符串表進行清空操作，為下一次工作做準備。當?shù)诙址頋M時，再選用第一片存儲器，如此循環(huán)進行操作。當輸入信號為空，即FIFo 1向FPGa發(fā)出空信號，表示數(shù)據(jù)傳送完畢，壓縮完成。

1.3 LZW算法用于超聲平面波數(shù)據(jù)壓縮實例

下面介紹LZW算法在超聲平面波數(shù)據(jù)壓縮中的應(yīng)用，通過實際血管的仿真數(shù)據(jù)進行實驗說明。假設(shè)從探頭中得到的數(shù)據(jù)集位寬是16位的，也就是說數(shù)據(jù)數(shù)值分布在-1215～215-1。而LZW算法的一個要求是輸入數(shù)據(jù)應(yīng)該是大于1的整數(shù)，需要將原始數(shù)據(jù)轉(zhuǎn)換成整數(shù)。

首先我們將數(shù)據(jù)分為4 096個數(shù)據(jù)塊，比如(1～8)，(9～16)，······，然后將原始數(shù)據(jù)映射到每個數(shù)據(jù)塊上。具體映射法如原始數(shù)據(jù)中某一位數(shù)值為2，可以發(fā)現(xiàn)它位于第一個數(shù)據(jù)塊中，我們用4或者5也就是數(shù)據(jù)塊的中間元素來替代它。對于原始數(shù)據(jù)的每一個元素都進行這樣的操作，最終得到的新數(shù)據(jù)的冗余度將大大增加。最后使用LZW算法對數(shù)據(jù)進行壓縮處理。處理結(jié)果如表1所示。

表1 LZW算法結(jié)果Tab.1 LZW compress result

LZW無損壓縮算法的壓縮比相較于頻域切割壓縮算法（將在后面介紹）并不理想，但是也能達到一定的壓縮比。產(chǎn)生此差別主要有三方面的原因，第一方面是無損壓縮本身壓縮比沒有有損壓縮算法高；第二方面是由于超聲波信號的特有性質(zhì)決定了信號數(shù)據(jù)的冗余度并不高，使得LZW壓縮算法沒有發(fā)揮其特點，在GIF動態(tài)圖壓縮方面，LZW壓縮算法就能將自己的優(yōu)勢發(fā)揮最大；第三方面是由于LZW算法除了輸出壓縮數(shù)據(jù)外，還需要傳輸解壓縮的字典串表，這使得輸出數(shù)據(jù)量增加了一倍，這也是壓縮比不高的主要原因之一。

2 頻域切割壓縮算法

LZW算法是針對于信號的時域部分來進行壓縮的一種算法。LZW算法不僅需要傳輸壓縮后的數(shù)據(jù)，還需要將解壓縮的字典串表一并發(fā)送給解壓縮端，這增加了數(shù)據(jù)的傳輸量也就縮小了壓縮比。而頻域切割算法則是針對超聲數(shù)據(jù)的頻域特征來進行數(shù)據(jù)壓縮的一種方法，而與LZW壓縮算法不同的是，頻域壓縮算法屬于有損壓縮。

2.1 算法描述

在敘述算法之前，首先說明一個分數(shù)帶寬的概念(Fractional Bandwidth, FBW)。分數(shù)帶寬指標主要描述一個超聲信號重要的主要頻率所占總體頻率的一個比例，其定義如下：

帶寬BW通常指信號所占據(jù)的頻帶寬度，帶寬有時候也稱為頻寬，以Hz為單位，而fc為中心頻率。分數(shù)帶寬更多是描述主頻部分信號占據(jù)帶寬的比例。如果分數(shù)帶寬確定，中心頻率越大，帶寬也就越大。所以，當一個系統(tǒng)確定時，它的中心頻率是確定的，可以通過公式求解出有效的帶寬。在超聲系統(tǒng)中，一般認為下降到-40 dB時的分數(shù)帶寬就是主頻信息基本保存完整的特征參數(shù)。所以在已知分數(shù)帶寬的條件下，我們可以通過此公式求出頻譜中需要保留的帶寬，而其余部分頻率信號對圖像影響的意義不大。算法實質(zhì)就是通過保留這段主要的頻帶信息通過后期數(shù)據(jù)的添加合成重建圖像。每個超聲探頭都有其自己的特征參數(shù)，如圖1示例。

圖1 探頭的頻域參數(shù)，中心頻率fc=9.2 MHz，及在-6 dB和-20 dB處的值Fig.1 The spectrum parameter of an aperture with fc=9.2 MHz, and value at -6 dB and -20 dB

如圖1所示，這個探頭發(fā)射超聲波的中心頻率為9.2 MHz，在-6 dB和-20 dB處的FBW值分別為88.4%和137.2%。如果想保證超聲成像不失真且主頻信號完整，本文規(guī)定將參數(shù)控制在157%左右，也就是-40 dB處。結(jié)合公式(1)，在知道探頭中心頻率的基礎(chǔ)上，可以計算出可通過帶寬，然后結(jié)合帶寬計算出相應(yīng)的起始點和結(jié)束點數(shù)。

頻域切割算法的主要流程如下：

1) 將原始一個FIFo中的RF-data（Radio Frequency data）進行FFT傅里葉變換；

2) 由探頭的中心頻率和FBW計算出BW，由頻譜圖計算起始和結(jié)束點；

3) 在頻域上將結(jié)束點以后的所有點都置為0，相當于截斷無用頻率信號點；

4) 將置零的個數(shù)存儲，只將主頻信號的數(shù)據(jù)和這個標記位上傳；

5) 將主頻數(shù)據(jù)進行IFFT傅里葉反變換，重建時域信號數(shù)據(jù)。

2.2 實驗結(jié)果

針對上述算法，本文使用Field II（一種超聲成像仿真軟件）產(chǎn)生的囊腫仿真數(shù)據(jù)進行測試，將128振元的數(shù)據(jù)先提取一條掃描線的數(shù)據(jù)來進行試驗，其他的數(shù)據(jù)壓縮效果基本一致。Field II相關(guān)仿真參數(shù)為：fc=3.5 MHz, fs=50 MHz, elements numbers=128, kerf=0.05/1000。

由于超聲圖像每條掃描線中所包含的頻率信息相似，所以，沒有必要將所有的掃描線都進行算法分析，但是評判后期圖像重建效果，每條掃描線進行算法實現(xiàn)是必要的，但是這并不影響壓縮比的計算。通過對一條掃描線數(shù)據(jù)的FFT變化，經(jīng)過對數(shù)壓縮操作后得到了探頭特征參數(shù)指標圖，如圖2所示。

圖2 囊腫仿真超聲數(shù)據(jù)頻域參數(shù)Fig.2 The spectrum parameter of the aperture with cyst simulated ultrasound data

本文設(shè)置-40 dB時的FBW為主頻通過的特征參數(shù)。下面通過圖3來具體闡述所進行頻域切割壓縮算法的步驟。

圖3中有四個單獨的信號圖，分別為a、b 、c 、d。圖3a表示一條掃描線的原始時域采集信號，我們通過在128中隨機選取一條掃描線的數(shù)據(jù)來進行試驗測試，發(fā)現(xiàn)圖像縱向深度大于5 000個點。然后對于圖3a的時域信號進行FFT傅里葉變換，得到頻譜圖。在頻譜圖中，可以看出信號的頻率以3.5 MHz為中心。然后通過計算BW的結(jié)束點差不多位于552個點，通過將其后的數(shù)據(jù)點置零得到壓縮后的頻域數(shù)據(jù)。

這里需要特別注意的一點是，圖3b中顯示的只是經(jīng)過FFT后對稱的一半信號，所以，置零的過程是將正負頻率部分的點都換成0。也就是說，不僅僅是500～5 000個點置為0，5 000以后與前端信號對稱的相應(yīng)4 500個點都得置為0。這是由于FFT變換會產(chǎn)生信號的正負頻率兩部分，而這兩部分的復(fù)數(shù)實部完全相同，但是虛部互為相反數(shù)。為了能夠正確地將切斷后的信號通過傅里葉反變換回去，就必須對應(yīng)置零，否則信號將無法正確反變換。

圖3 頻域切割壓縮過程Fig.3 The process of the frequency domain compression

圖3c顯示的就是由圖3b頻率信號經(jīng)552點以后切斷的頻譜圖。對稱的另一半頻率圖缺省顯示。圖3d就是由圖3c所截斷置零以后的頻域信號經(jīng)過IFFT傅里葉反變換還原的時域信號。

以此類推，另外的128個通道的數(shù)據(jù)以同樣的方式進行FFT變化，通過計算FBW和相應(yīng)的帶通點數(shù)，將結(jié)束點后的點進行舍棄操作，并將數(shù)據(jù)量用標志位進行記錄。最后將壓縮后的數(shù)據(jù)上傳電腦后在電腦中實現(xiàn)數(shù)據(jù)的重建，最后成像。圖4為原始圖像與重建圖像的對比。

從圖4中，如果僅僅用人眼來區(qū)分，很難發(fā)現(xiàn)兩幅圖像存在的細微差別。這也正好說明了頻域壓縮算法的可行性。但是為了更具說服力，本文通過計算兩幅圖像的后期評價指標來進行量化比較。

圖5顯示了不同深度下的縱向分辨率。它們大多集中在246.4 μm，原始圖像的分辨率稍好于重建以后的圖像。圖中小正方形的數(shù)值代表具體的分辨率大小。從五個不同深度的分辨率可以得出結(jié)論，頻域切割壓縮算法對于圖像分辨率的影響微乎其微。

圖4 原始圖像與重建圖像的對比Fig.4 comparison of the original image and the reconstructed image

圖5 不同深度的縱向分辨率Fig.5 axial resolutions at different depths

表2顯示了原始圖像與重建圖像在中間深度對比度與cNR (contrast to Noise Ratio)的值，從表2中可知，對比度比值為64:61，cNR則只相差0.04 dB。在這兩項指標上，原始圖像和壓縮重建以后的圖像差別十分微小。所以不管是從對比度，分辨率還是cNR上，我們可以得出結(jié)論：重建后的圖像與原始圖像相似度極高，頻域切割壓縮算法的可行性和正確性得到了驗證。

表2 對比度與cNR在中間深度的比較Tab.2 The contrast and cNR in the mid-depth comparison of the two images

通過這樣的壓縮算法到底為我們節(jié)省了多少資源？壓縮比的大小見表3。

表3 頻域切割算法結(jié)果壓縮結(jié)果Tab.3 one line data compression result

一條掃描線的總數(shù)據(jù)量為5 010個點，頻域切割點從552個點開始。由于記錄的點數(shù)是5 010個點，硬件上實現(xiàn)時只需要記錄下來沒有截取的數(shù)據(jù)和后續(xù)需要截取多少個0的標志位，而標志位在壓縮比計算式忽略不計。所以壓縮比的大小僅是數(shù)據(jù)總體長度5 010個點和截取剩下的552個點的比值，達到了9.076。而根據(jù)這個主頻數(shù)據(jù)點和置零的標志位我們可以在Pc機上順利地進行傅里葉反變換，重建出原始數(shù)據(jù)的5 010個點。對于128振元的通道依次處理將大大壓縮平面波所帶來的大數(shù)據(jù)量。

3 結(jié)束語

本文主要從時域和頻域兩個方面提出了兩種對超聲平面波大數(shù)據(jù)量進行壓縮的方法，主要工作是在科學院深圳先進技術(shù)研究院生物醫(yī)學與健康工程研究所完成的。兩種方法各有它們的特點，主要針對當前硬件處理芯片無法及時傳輸大批量數(shù)據(jù)所提出。當使用頻域切割壓縮算法時，必須充分了解自己所要處理的超聲信號特點，這樣才能夠給出正確的傅里葉變換和傅里葉反變換的方法。這種處理方式不管是在硬件還是軟件上都是易于實現(xiàn)的。而LZW算法不需要將數(shù)據(jù)轉(zhuǎn)換到頻域進行計算，在硬件上實現(xiàn)操作也不是很難，但是LZW算法的壓縮比由于需要字符串表的傳輸而大大減少。

從壓縮比角度，頻域切割壓縮算法是進行超聲數(shù)據(jù)壓縮的算法優(yōu)先選擇。但要注意的是，頻域切割壓縮算法是一種有損壓縮方法，雖然本文根據(jù)后期評價指標比較原始圖像和壓縮重建圖像的差異來看，結(jié)果是微小的，但這種比較還是初步的。最終的比較還要根據(jù)最終的應(yīng)用來判斷頻域切割壓縮算法是否有影響。

[1] Tanter M, Fink M. Ultrafast imaging in biomedical ultrasound[J]. Ieee Trans Ultrason Ferroelectr Freq control, 2014, 61(1):102-119.

[2] couture o, Fink M, Tanter M. Ultrasound contrast plane wave imaging[J]. Ieee Trans Ultrason Ferroelectr Freq control, 2012, 59(12): 2676-2683.

[3] Garcia D, Le Tarnec L, Muth S,et al. Stolt's f-k migration for plane wave ultrasound imaging[J]. Ieee Trans Ultrason Ferroelectr Freq control, 2013, 60(9): 1853-1867.

[4] owen K, Fuller MI, Hossack Ja. application of x-y separable 2-D array beamforming for increased frame rate and energy effciency in handheld device[J]. Ieee Trans Ultrason Ferroelectr Freq control, 2012, 59(7): 1332-1343.

[5] asl BM, Mahloojifar a. eigenspace-based minimum variance beamforming applied to medical ultrasound imaging[J]. Ieee Trans Ultrason Ferroelectr Freq control, 2010, 57(11): 2381-2390.

[6] Gallot T, Perge c, Grenard V, et al. Ultrafast ultrasonic imaging coupled to rheometry: Principle and illustration[J]. Rev Sci Instrum, 2013, 84(4): 045107.

[7] 彭勇. 一種改進的LZW編碼數(shù)據(jù)壓縮算法的設(shè)計與仿真[J]. 信息與電腦, 2012, (2): 55-57.

[8] 劉荔. 全數(shù)字B超診斷儀中數(shù)字信號處理的FPGa設(shè)計與實現(xiàn)[D]. 杭州: 浙江大學, 2006.

Preliminary Comparison of Two Compression Methods for Ultrasonic Plane Wave Data

【 Writers 】ZHANG Yaonan1,2, JIN Yingdong2,3, QIU Weibao3
1 College of Electronics and Information Engineering, Xi’an Siyuan University, Xi’an, 710038
2 Sino-Dutch Biomedical and Information Engineering School, Northeastern University, Shenyang, 110169
3 Institute of Biomedical and Health Engineering, Shenzhen Institutes of Advanced Technology, Chinese Academy of Science, Shenzhen, 518055

Ultrasound, plane wave, data compression, LZW

R 445.1

A

10.3969/j.issn.1671-7104.2016.03.003

1671-7104(2016)03-0164-05

2016-01-29

遼寧省自然科學基金資助項目(201202071)

張耀楠，教授，e-mail: johanzyn@qq.com

【 Abstract 】Ultrasonic plane wave imaging is one of the hot topics in the study of ultrasound imaging, but the plane wave ultrahigh frequency ultrasound has brought a large amount of data. If these data are transmitted directly through the hardware, it will not only bring a great burden on the hardware, but also put a requirement for the ability of the hardware processing data. In order to reduce the problem caused by the high frame frequency and large amount of data, the data compression methods are used to reduce the amount of data acquisition and transmission. In this paper both lossy compression and lossless compression are studied, where LZW compression algorithm is chosen to test lossless compression effect, and based on the plane wave characteristics, the frequency domain cutting compression algorithm is developed as a lossy compression method. In this paper, the LZW algorithm is implemented in FGPA, and it can achieve a certain compression ratio after testing the ultrasonic simulation data. Under the condition of known fractional bandwidth, it is found that the main bandwidth can be preserved in the spectrum, while the rest of the frequency signal has little effect on the image. By testing the ultrasonic simulation data, the frequency domain cutting algorithm can achieve high compression ratio, and the difference between the original image and the reconstructed image is relatively small. So from the point of the compression ratio, the frequency domain cutting algorithm is the primary choice of the ultrasonic data compression algorithm.