背景技术
医疗超声系统通过从由临床医生放置在受试对象上的换能器发送超声波束,扫描受试对象的内部解剖结构。超声波在具有不同声阻抗的内部组织的界面上反射,产生回声。换能器接收回声并将它们转换为电超声信号。超声系统将一序列的处理步骤应用到超声信号以产生在控制台上显示的图像或一系列图像,用以由临床医生分析。基于接收的回声的强度形成的图像称为B-模式图像。此外,系统可以测量超声信号的多普勒频移(Doppler shifts)以产生指示比如血液之类的流体的流动的彩色图像,以及执行对诊断有用的额外的分析。
常规的医疗超声换能器包括在由电信号驱动时发送超声波,接收返回的回声并且将接收的回声转换为多个模拟信号的压电元件阵列。多个模拟到数字转换器(analog to digital converter,简称ADC)对模拟信号采样,每个都产生数字信号样本流。信号样本的典型的数字信号处理包括波束形成、下变频转换、B-模式(亮度)处理和/或多普勒处理、扫描转换以及用于显示的图像处理。波束形成延迟并组合信号样本流以形成对应于视场中的特定方向的波束形成样本阵列。波束形成器能够产生对应于视场中的数个方向的数个波束形成样本阵列。取决于期望的诊断信息类型,随后在波束形成样本上执行B-模式处理和/或多普勒处理以形成B-模式检测样本和/或多普勒检测样本。检测样本的空间坐标仍然对应于波束形成样本的波束几何形状。扫描转换器执行检测样本的坐标转换以产生具有适合用于显示的光栅格式的数据帧。将额外的图像处理应用到样本帧以允许它们作为二维(2-D)或三维(3-D)图像显示。
用以改进医疗超声系统的当前努力针对于提高控制台/车系统的诊断能力以及开发具有提高的图像质量的较小的便携式设备。对于高端控制台或车系统,增加换能器元件的数量以产生更高分辨率和/或3-D图像以扩展诊断能力是期望的。增加换能器元件的数量会增加从换能器头传达到控制台处理器的数据量,其可能需要更高带宽的通信信道以及更大的电缆连接。换能器头的数据采集能力由操纵需要和外形要素所限制。手提或手持式超声设备对于在小诊所、移动治疗单元以及家庭中的使用是经济的和期望的。对于这些设备,电池寿命也是一种限制。在超声系统中更加高效的处理、传输和存储超声信号数据可以节省功率、数据传输带宽和存储器容量。
超声信号数据的压缩能够对控制台/车系统和便携式系统两者提供益处。益处包括降低系统的数据传输带宽、存储器容量和功率要求。对于便携式或者手提式超声系统,这些益处会降低重量并且增加电池寿命。对于控制台系统,压缩会减轻增加由换能器头获取的数据量并传送数据到超声信号处理器的影响。在计算上高效的压缩在对系统复杂度具有较低影响或者没有影响的情况下引入压缩的益处。
本描述使用术语“压缩(compression)”来指代超声信号样本的数据压缩,其中代表信号样本的比特数目被减少并且信号样本随后在用于显示的处理之前解压缩。超声成像系统的某些描述使用术语“压缩”来表示“脉冲压缩”而不是数据压缩。脉冲压缩是指过滤和/或调制发送的超声脉冲,以及逆向过滤和/或解调接收的超声脉冲。(例如,见作者为V.Behar与D.Adam,于2004年在Ultrasonics的卷号42,页码1101-1109中发表的“Parameter optimization of pulsecompression in ultrasound imaging system with coded excitation”)。超声成像系统的某些描述使用术语“压缩”来表示“对数压缩”而不是数据压缩。在该上下文中,对数压缩是指在显示之前计算处理的超声数据,典型地为检测的数据的幅度的对数。(例如,见作者为A.Murtaza等,于2008年11月在Texas Instruments SPRAB 12的页码1-26中发表的“Signal Processing Overview of Ultrasound”)。脉冲压缩和对数压缩两者都有意地改变发送的或接收的超声信号在时域和频域中的特征。在随后跟随有解压缩的接收的超声信号样本的数据压缩是在时域和频域中保存信号特征的处理。本描述指的是超声信号样本的无损及有损压缩。在无损压缩中,解压缩的样本具有与原始样本相一致的值。在有损压缩中,解压缩的样本类似,但不同于原始样本。本描述使用术语“帧”以指代最终将被处理以形成用于显示的超声图像的原始的或处理的超声数据阵列。在本领域中超声成像系统的某些描述还使用术语“屏”以指代超声数据帧。在本描述中,“实时”表示至少与数字信号的采样速率一样快的速率。可以使用术语“实时”来描述数字信号的处理、传输以及存储的速率。采样速率是ADC在模拟信号的转换期间形成数字信号样本的速率。在本领域中超声成像系统的某些描述使用“实时”以指代显示超声图像的帧速率。本描述将实时与采样速率,而不是帧速率解释相关联。
在超声系统中的数据压缩的从前的申请包括了在用于图像形成的扫描转换之前和之后的数据压缩的替代方式。在公告于2001年11月13日,名称为“Ultrasonic Diagnostic Device”的美国专利6,315,722中,Yaegashi描述了用于存储输出自ADC单元的超声信号样本的时间轴扩展单元。时间轴扩展单元以输出自ADC单元的速率写数据并以较低的速率读数据。时间轴扩展单元存储用于一屏或一帧的信号样本,并且可以使用先入先出(first-in first-out,简称FIFO)存储器实施。数据压缩单元压缩从时间轴扩展单元读取的信号样本。Yaegashi描述图像压缩技术的应用,这样的方法基于用于利用一个数据帧之内的空间相关性的离散余弦变换(discrete cosine transform,简称DCT)或者用于多数据帧的MPEG压缩方法。(MPEG是指由Moving Picture Experts Group开发的视频数据压缩标准)。压缩的样本存储于大容量存储设备,比如硬盘中。数据压缩减少在大容量存储设备中所需的存储容量。为了产生图像,数据扩展单元解压缩取回自大容量存储设备的压缩的样本。将包括过滤、对数转换、检测以及数字扫描转换在内的常规操作应用到解压缩的样本,用于图像形成和显示。Yaegashi没有公开在处理序列中的波束形成。
在公开号为2008/0114246,名称为“Transducer Array ImagingSystem”的美国专利公开中,Randall等人描述了在波束形成之前和/或之后使用映射、重采样和/或数据开窗而压缩超声数字数据。映射可以包括再量化或者裁剪信号样本。例如,所需比特的数量随深度单调地减少从而基于深度更少的比特可以分配给每个样本。在某些实施方式中,来自延伸出发送与接收孔径的接收通道的信号样本可以被截去。为了对感兴趣区域(region of interest,简称ROI)成像,信号采集时间可以与深度范围成正比,从而使在最小采样时间之前和/或最大采样时间之后获取的数据能够被截去--如果它们不贡献于图像像素的形成。在某些实施方式中,如果显示分辨率低于全分辨率成像的需要,那么数据可以被重采样至更少的样本,从而减少传输的样本的数量。
在公告于2000年3月28日,名称为“Medical DiagnosticUltrasound System and Method for Transform Ultrasound Processing”的美国专利6,042,545中,Hossack等人描述用于在波束形成之后的超声数据的变换压缩技术。波束形成的替代方式包括ADC之前的模拟波束形成或者ADC之后的数字波束形成。波束形成器产生同相位和正交的(in-phase和quadrature,简称I和Q)样本,或者替代地,射频(RF)样本。对应于二维(2-D)帧的波束形成样本被过滤并变换以产生变换域表示。变换域样本被量化和/或编码用以压缩。压缩可以是无损的或有损的。可以应用任何变换,比如DCT或离散小波变换(Discrete Wavelet Transform,简称DWT)、量化功能,以及编码功能,用以压缩数据帧。例如,JPEG压缩包括将数据帧分割成2-D数据块、在每个块上使用2-D DCT进行变换、量化变换域样本、差分编码在块之间DC(零频率)变换样本,以及熵编码量化的变换域样本的2-D块(例如,Huffman编码)。JPEG压缩算法可以配置为有损的或无损的。(JEPG压缩是指由Joint Photographic ExpertsGroup开发的标准图像压缩方法)。在变换域中对于比如过滤之类的各种图像处理功能的额外的操作在变换域中比在空间域中在计算上更加高效。例如,空间域中的2-D过滤使用2-D卷积运算。在变换域中2-D过滤通过变换域过滤系数使用更加高效的乘法。压缩的变换域数据可被存储用于之后的图像形成。为了解压缩,在显示处理之前应用逆向编码与变换功能。
在公告于2005年2月15日,名称为“Diagnostic InformationGeneration Apparatus and Ultrasonic Diagnostic System”的美国专利6,855,113中,Amemiya等人描述了在从超声波单元到信息单元的无线传输之前的超声数据帧的压缩。超声波单元包括换能器和处理器(用于随后的波束形成、B-模式成像以及多普勒成像)。将通用数据压缩标准应用到B-模式成像数据或多普勒成像数据,例如用于单帧的JPEG压缩或者用于多帧的MPEG压缩。压缩的数据使用标准无线通信方式发送到信息单元。信息单元包括根据压缩标准解压缩接收的数据的中央处理单元(central processing unit,简称CPU)。CPU此外处理解压缩的B-模式成像数据和解压缩的多普勒成像数据用于显示。
在公开于1997年3月20日,名称为“Ultrasonic DiagnosticApparatus for Compressing and Storing Data in CINE Memory”的PCT公开申请,国际公开号为WO 97/09930中,Lee描述了在CINE存储器中存储之前对超声数据进行压缩,并且对从CINE存储器中取回的数据进行解压缩。CINE存储器包括几个由时间组织的库。在此系统中,超声探测器在ADC之前执行波束形成,因此ADC输出数据代表波束形成样本。对数据帧应用压缩并且可以在扫描转换之前或者之后应用。应用Lempel-Ziv-Welch(LZW)算法用以压缩和解压缩。LZW算法基于检测数据中比特的重复模式并且为重复模式分配代码。取回自CINE存储器的帧的压缩的数据被解压缩并且进一步处理用于显示。
在公开于2005年3月31日,名称为“Ultrasonograph andUltrasonic Data Compression Method”,公开号为2005-081082的日本专利申请中,Akihiro描述了三种用于在模拟波束形成之后压缩超声数据的实施方式。在第一实施方式中,ADC生成模拟波束形成器输出信号的I和Q样本。压缩器计算相邻波束的I、Q样本的差异并在随后通过对差异进行游程长度编码而形成压缩的数据。压缩的数据存储于存储器中。取回自存储器的压缩的数据被解压缩和处理用于图像显示。在第二实施方式中,ADC生成模拟波束形成器输出样本的RF样本。压缩器计算相邻波束的RF样本之间的差异并在随后进行游程长度编码。压缩的样本被存储于存储器中、取回、解压缩并处理用于图像显示。在第三实施方式中,波束形成器输出被进一步处理以在压缩之前生成B-模式图象帧和多普勒图像帧。压缩器计算帧到帧的差异以产生压缩的数据帧。压缩的数据帧被存储于存储器中、取回、解压缩并进行进一步处理用于显示。
在公告于1988年6月21日,名称为“Ultrasonic BloodstreamDiagnostic Apparatus with Dual Displays of Velocity Profiles andAverage Flow Velocity”的美国专利4,751,929中,Hayakawa等人描述了对多普勒频率检测数据的压缩。压缩器在计算频谱样本的实部和虚部的幅度平方的平方和加法电路的输出上操作。压缩器重新编码输出自加法器的每个样本的比特,以减少在表示中的比特数量。压缩器在加法器的输出样本上操作,以编码尾数中的最重要比特的位置,保留固定数量的最重要比特并且清除剩余的最不重要的比特。由此产生的每个样本的压缩的字都包括固定数量的最重要比特以及表示从原始样本中清除的最不重要比特的数量的编码。将可变数量的最不重要比特从每个样本中去除,所以压缩是有损的。
在标题为“A Novel B-Mode Ultrasound Image CompressionMethod Based on Beam Forming Data”,1998Proc.Intl.Conf.IEEEEngineering in Medicine and Biology Society,Vol.20 No.3,pp.1274-76中,Li等人描述了对用于在远程超声系统中发送的波束形成样本的压缩。DWT应用在128x512的波束形成样本帧之中。在竖直方向中的子图像系数使用算术编码而被量化和编码。在解压缩之后,将扫描转换应用于128x512的解压缩的样本帧以形成512x512的样本帧用于显示。
数篇论文描述了用于在用以图像形成的扫描转换之后压缩超声图像的不同方法。一些示例包括以下。在标题为“Comparative Surveyof Ultrasound Images compression Methods Dedicated to aTele-Echography Robotic System”,2001Proc.23rd Annual IEEEEngineering in Medicine and Biology Society Intl.Conf.,pp.2461-64的论文中,Delgorge等人描述了对超声图像应用不同的压缩方法。方法包括傅立叶变换(Fourier transform)、DCT、四叉树分解、DWT、分形、直方图阈值以及游程长度编码。方法在扫描转换之后应用在512x512的超声图像。在标题为“Despeckling of Medical UltrasoundImages Using Data and Rate Adaptive Lossy Compression”,IEEE Trans.Medical Imaging,vol.24,No.6,2005年6月,pp.743-54的论文中,Gupta等人描述了将压缩与算法相结合以从超声图像中去除斑点。在斑点去除算法、量化以及熵编码后进行DWT。在标题为“ATele-Operated Mobile Ultrasound Scanner Using a Light-WeightRobot”,IEEE Trans.Information Technology in Biomedicine,Vol.9,No.1,2005年3月,pp.50-58的论文中,Delgorge等人描述了对超声图像应用多种无损和有损压缩方法。无损方法包括Huffman、算术编码、Lempel-Ziv、游程长度编码以及Fano编码。有损方法包括各种JPEG版本,包括JPEG、JPEG-LS以及JPEG2000。在标题为“MaximumLikelihood Motion Estimation in Ultrasound Image Sequences”,IEEESignal Processing Letters,Vol.4,No.6,1997年6月,pp.156-7的论文中,Strintzis等人描述了对超声图像序列应用MPEG压缩。方法包括检测在图像序列中连续帧之间的8x8的像素块的运动矢量。运动矢量被编码用于帧到帧的MPEG压缩。
公告日期为2006年3月7日的名称为“Adaptive Compression andDecompression of Bandlimited Signals”的共同所有的美国专利7,009,533(‘533专利)描述用于特定带限信号的压缩与解压缩的算法。公告日期为2007年8月8日的名称为“Enhanced Data ConvertersUsing Compression and Decompression”的共同所有的美国专利7,088,276(‘276专利)描述对输出自实施于单一集成电路之中的ADC的信号样本应用无损或有损压缩。提交于2008年5月15日的名称为“Digital Interface for Data Converters”的共同所有及共同未决的序列号码为12/120988(‘988申请)的美国专利申请描述复用并行输出自多个ADC的数据,以减少在数字接口的活动数据端口的数量。
存在对超声成像系统的部件之间的超声信号数据的高效数据传输和存储的需求。存在对超声信号数据的计算上高效的数据压缩以在对系统复杂度产生最小影响的情况下提高数据传输和存储容量的需求。
具体实施方式
图1为依据现有技术的典型医疗超声系统的示例的框图。发送波束形成器104为现有技术中的已知构造,比如数字或模拟波束形成器。发送波束形成器104响应于系统控制器102生成一个或多个激励信号。激励信号具有相关联的中心频率,典型地在1-20MHz范围。来自发送波束形成器104的激励信号经由发送/接收切换器112提供至超声换能器110。超声换能器110包含换能器元件110i的阵列。超声换能器110是允许将超声波耦合到被检查的受试对象的现有技术中的已知构造。换能器元件110i发射并且接收超声波。发送/接收切换器112包括用于发送和接收模式的交换电路。对于发送模式,发送/接收切换器112将来自发送波束形成器104的激励信号耦合至换能器110。对于接收模式,发送/接收切换器112将接收的超声信号从换能器110耦合至模拟前端(AFE)114。为了传输,换能器元件110i转换激励信号以产生超声波声波波形。特别是,换能器110响应于发送波束形成器104,将激励信号转换成在受试对象内的方向上前行的超声波波形。具有拥有不同声阻抗的界面的散射位置反射超声波波形,产生会返回换能器110的回声。多个换能器元件110i接收回声并将它们转换成多个模拟超声信号。发送/接收切换器112在采样窗口期间将来自换能器110的多个模拟超声信号耦合到AFE114。采样窗口对应于时间间隔,在其中接收的回声代表来自在受试对象中期望的深度范围内的散射位置的反射。控制器102依据用户输入或者扫描协议设置采样窗口并将定时控制信息提供给发送/接收切换器112。发送/接收切换器112在采样窗口期间输出并行的多个模拟超声信号。AFE114将多个模拟超声信号进行放大并过滤以对模拟到数字转换进行准备。对于每个模拟信号通道113i,AFE114可以包括低噪声放大器(low noise amplifier,简称LNA)、可变增益放大器(variable gain amplifier,简称VGA)以及低通滤波器。因为接收的信号强度随时间减弱,所以VGA应用作为时间的函数增加增益的增益分布。信号强度随时间减弱是由超声波在穿过更多组织的更长路径上的衰减造成的。ADC库120包括多个ADC以将在采样窗口期间接收的多个模拟超声信号转换为多个并行的超声信号样本序列。在每个ADC输入通道115i的模拟超声信号被转换为在相应的ADC输出通道121i的超声信号样本流。超声信号样本具有非零中心频率,通常对应于与换能器的压电材料的本质、谐振频率相关的接收的超声信号的射频(RF)。
超声信号处理器130执行从超声信号样本中提取诊断信息的操作,包括关于图2进一步描述的波束形成、下变频转换、B-模式处理和/或多普勒处理。超声信号处理器130可以实施为一个或多个处理器,所述处理器比如数字信号处理器(digital signal processor,简称DSP)、现场可编程门阵列(field-programmable gate array,简称FPGA)、微处理器或者计算机的CPU。扫描转换器140执行处理的样本帧的坐标变换以产生具有光栅(正交直线)格式的超声图像样本帧。图像处理器150能够对超声图像样本应用额外的图像增强操作。显示器160提供二维或三维图像供用户分析。
图2为依据现有技术的超声信号处理器130的框图。医疗超声系统能够在输出自ADC库120的RF超声信号样本上执行数字波束形成操作。接收波束形成器160对超声信号样本应用延迟、变迹(衰减)以及加法操作以形成对应于视场中的特定方向的波束形成样本,或者波束的1-D阵列。接收波束形成器160产生对应于视场中的多个方向的多个波束形成样本的1-D阵列。接收波束形成器160包括波束形成器(BF)存储器162以存储用于波束形成(BF)计算器164的延迟和加法操作的超声信号样本。针对同一接收的脉冲,BF计算器164能够使用取回自BF存储器162的超声信号样本形成多个波束。BF计算器164能够在相邻超声信号样本之间应用插值算法,以提高计算的波束的相位分辨率。接收波束形成器160还能够在加法操作之前对样本应用加权函数,以实施空间窗口功能或者变迹。为每个角度计算的波束形成样本被提供给相应的波束形成器输出通道161i。接收波束形成器160通常具有少于输入通道的,包含ADC输出通道121i的输出通道161i。在这种配置中,波束形成样本具有RF中心频率。数字下变频转换器(digital down converter,简称DDC)170将波束形成样本解调至基带,以生成每个波束的复基带I和Q样本。作为DDC170的替代或补充,可以对在以期望的频率为中心的频带处的波束形成样本应用带通滤波器,或者DDC170可以将波束形成样本解调为中间频率(intermediate frequency,简称IF)而不是基带。在本领域中的替代结构包括在模拟到数字转换之前的模拟波束形成以及在波束形成之前的超声信号样本的数字下变频转换。
诊断信息处理器180为期望的类型的超声图像在I、Q样本上执行合适的操作。B-模式处理生成代表回声信号强度的信息。I、Q样本的幅度可被计算以形成用于B-模式成像的检测样本。多普勒处理从I、Q样本中估算速度、速度变化以及能量以形成多普勒检测样本。B-模式检测样本和多普勒检测样本的空间坐标对应于波束形成样本的几何形状。扫描转换器140执行检测样本的坐标变换以产生具有适合于显示的光栅格式的数据帧。图像处理器150在将样本帧作为二维或三维图象显示之前执行样本帧的额外的图像处理。
图3为在波束形成之前的超声信号样本的绘图。绘图显示由ADC阵列中的一个ADC所采样的四个脉冲回声的同相位样本。对于这个示例,数字下变频转换在波束形成之前应用到输出自ADC的超声信号样本,以形成I、Q样本。图4为对应于一个波束形成器输出通道的波束形成超声信号的同相位样本的绘图。对于这个示例,波束形成器通过对I、Q样本执行延迟与加权函数而结合输出自多个ADC的多个I、Q样本序列。
图5为依据优选实施方式的,包括压缩输出自ADC库120的超声信号样本的超声系统的框图。压缩器210包括多个压缩单元210i。多个压缩单元210i压缩代表输出自相应的换能器元件(例如图1的110i)的模拟超声信号的多个信号样本序列,以形成在压缩器输出端211i的多个压缩的样本序列。压缩单元210i接收来自相应的ADC120i的ADC输出的121i的,代表在采样窗口期间输出自相应的换能器元件110i的模拟超声信号的信号样本序列,其中采样窗口为在其中模拟超声信号代表从受试对象中的相应深度范围反射的回声的时间间隔。压缩单元210i独立于代表在相同的采样窗口期间输出自其他换能器元件的模拟超声信号的信号样本,对在它的输入端的特定信号样本序列应用压缩操作,以在压缩器输出端211i产生相应的压缩的样本序列。或者,压缩单元210i可以压缩来自一个以上ADC的样本,其中来自特定ADC 120i的信号样本独立于代表在相同的采样窗口期间输出自其他换能器元件的模拟超声信号的信号样本而被压缩。当前描述假设每个ADC 120i产生实数值样本序列。或者,如果ADC 120i执行正交采样以在每个ADC输出端121i生成I和Q样本序列,那么每个I、Q样本序列被独立地压缩。压缩在波束形成之前应用到超声信号样本。压缩的样本跨越数字接口220传输到在其中它们被解压缩用于超声处理操作的超声信号处理器130。数字接口220用于传输压缩的样本所需的数据传输带宽相较于用于传输未压缩的超声信号样本流的带宽是减小的。
图6为依据替代配置的,包括压缩超声信号样本并且在捕获存储器中存储压缩的样本的超声系统的框图。在跨越数字接口220传输之后,压缩的样本存储于捕获存储器230之中。足够用于存储压缩的样本的捕获存储器230的容量减少,从而节省系统资源。解压缩器240解压缩取回自捕获存储器的压缩的样本。接收波束形成器160对解压缩的样本应用波束形成操作以在每个波束形成器输出通道161i形成波束形成样本序列。如关于图2描述的,可以在波束形成操作之前对解压缩的样本应用数字下变频转换和/或过滤。
图7为依据替代配置的,其中压缩的样本存储于波束形成器存储器之中的超声系统的框图。接收波束形成器160将压缩的样本存储于波束形成器(BF)存储器162之中并且当波束形成操作需要时取回压缩的样本。解压缩器240处理取回自BF存储器的压缩的样本并且将解压缩的样本提供给BF计算器164。BF计算器164对解压缩的样本应用波束形成操作。BF计算器164还能够在解压缩的样本上执行关于图2描述的,用于空间过滤等的额外的操作。
由于压缩的样本相对于原始超声信号样本具有较少的每样本比特,压缩的样本能够使用比未压缩的样本的较少的数据端口跨越数字接口220传输。图8为包括复用压缩的样本用以跨越较少数据端口传输的超声系统的框图。对于这个示例,ADC库120包括N个ADC 120i,ADC1到ADCN,在N个ADC输出通道120i产生信号样本流X1到XN。压缩单元210i在N压缩器输出端211i产生相应的压缩的样本流Y1到YN。具有降低的比特率的压缩的样本被提供给数据端口270i。优选地,数据端口270i提供低电压差分信令(low voltagedifferential signaling,简称LVDS)数据传输。或者,数据端口270i可以使用串行器-解串行器(serializer-deserializer,简称SerDes)接口用于数据传输。由National Semiconductor于2008年1月出版的,标题为“LVDS Owner’s Manual Including High-Speed CML and SignalConditioning”的文件第四版中描述了LVDS设备与结构。LVDS数据传输具有期望的特性,包括3.125Gbps的最高数据传输速率、低噪声以及低功耗。差分信令需要在每个通道输出端具有两个I/O针脚,一个用于正差分输出而另一个用于负差分输出,被称为LVDS对。当在压缩器输出端Yi的压缩的样本的比特率低于端口的最高数据传输速率时,数据端口270i具有过剩带宽。通过组合通过给定数据端口270i传输的多个压缩样本流,可以利用过剩带宽。复用器250将N个压缩的样本的序列进行组合以形成M个复用的压缩的样本的序列,其中M<N,用以在M个数据端口270i上传输。可以组合的压缩的样本流的数量受限于数据传输端口270i的带宽。例如,对于N=16ADC,其中每个ADC以每秒50兆样本(megasamples per sec.,简称Msps)的采样速率产生具有每样本12比特的超声信号样本,在每个ADC输出通道121i的比特率为600Mbps。假设数据端口270i具有高达800Mbps的数据传输速率并且压缩器210产生3∶1的压缩比。在压缩器输出端121i的压缩的样本序列Yi具有200Mbps的比特率。在这种情况下,一个数据端口270i具有足够的带宽以800Mbps的数据传输速率传输来自四个压缩器输出端的压缩的样本。对于这个示例,复用器250组合来自四个压缩的序列Yi、Yi+1、Yi+2和Yi+3的组中的压缩的样本,以在一个复用器输出端251i形成相应的复用的压缩的样本序列Zj。在这个示例中,复用的压缩的样本序列与相应的活动数据端口的数量为M=N/4=4。串行器260向相应的数据端口270i提供比特序列。
复用压缩的样本以使用M个数据端口的优点包括使用较少的物理数据端口,这反过来减少数据端口的连接并且降低其功耗。在其中压缩的样本的比特率被固定的应用中,复用器250具有对于固定数量的物理数据端口270i的固定数量的复用器输出端251i。或者,灵活架构能够依据压缩的样本的比特率提供可变数量的活动数据端口。对于灵活架构,压缩器210依据用户可选的压缩比参数,在不同的比特率提供压缩的样本。压缩控制器(未在图8中示出)为压缩单元210i提供压缩控制参数,从而操作产生具有对应于期望的压缩比的比特率的压缩的样本。压缩控制器为复用器250提供复用控制参数以指示需要在复用器输出端251i产生的复用的压缩的样本序列的数量。压缩控制器还能够响应于用户输入断电不活动数据端口,从而进一步节能。
数据端口270i仍可以对于特定的比特率和每样本比特的值具有过剩带宽,即便不压缩信号样本的情况下也是。当采样速率与每样本比特数量的乘积的结果为在ADC输出端121i的比特率低于数据端口270i的最高数据传输速率的比特率时,数据端口270i具有过剩带宽。例如,假设ADC时钟频率为50兆赫兹(MHz),对应于50Msps,并且样本宽度为每样本12比特从而每个ADC 120i每秒产生600兆比特(Mbps)。对于每个数据端口270i,数据传输速率将会是600Mbps。对于具有800Mbps的最高数据传输速率的LVDS端口,每个端口有200Mbps的未使用带宽。数据端口270i在600Mbps上传输数据时将会消耗与当它们在800Mbps的满容量上传输数据时相同的功率,因为LVDS接口是恒定电流(从而是恒定功率)连接。端口集中通过组合来自多个ADC输出端的比特并且在较少的端口上以更快的速率传输组合的数据而利用过剩数据传输容量。每个数据端口在高于ADC比特率并且直到每个数据端口270i的最高数据传输速率的速率上传输组合的数据。
图9为对ADC输出端121i应用端口集中的超声系统的框图。在N个ADC输出端X1到XN的N个超声信号样本序列被输入到复用器252。复用器252组合来自一组ADC输出端组的样本以形成M条复用的信号样本通道R1到RM,其中M<N。每条通道Ri具有高于在ADC输出端121i的ADC比特率,但是低于或者等于数据端口270i的以每秒比特计的最高数据传输速率的比特率。复用器250通过将每个采样周期的N个信号样本的比特映射到M个比特子集而组合样本。串行器262串行化每个比特子集用以由相应的数据端口270i传输。考虑上面的例子,其中额外地,ADC的数量N为16。以600Mbps的ADC比特率输出自16个ADC 120i的信号样本可以组合成12个比特流,每个具有800Mbps的比特率。12个活动数据端口270i将会传输这12个每一个在800Mbps(800Mbpsx12=600Mbpsx16)的数据传输速率上的比特流。取决于采样速率、每样本比特数量以及最高数据传输速率,可以有几种配置。
图10A为包括用于传输来自十六个ADC的数据的十二个数据端口的端口集中配置的框图。配置能够容纳对应于高达50MHz的ADC时钟频率的高达50Msps的采样速率。ADC输入时钟118带有频率对应于ADC采样速率的时钟信号。ADC 120响应于ADC输入时钟118,以采样速率对输入模拟信号进行采样。复用器252接收输出自ADC库120的信号样本。四个ADC 120h、120i、120j和120k为一个复用器252i提供信号样本。复用器252i包括三个输出端253i、253j和253k。复用器252i重新排序它在每个采样周期期间接收的样本的48个比特并且分配重新排序的比特到三个串行器262i、262j和262k。每个串行器262i、262j和262k接收48个重新排序的比特子集,在此例中为16个,并且串行化它们用以由相应的端口270i、270j和270k传输。锁相环(phase lock loop,简称PLL)280产生具有16倍于ADC时钟频率的频率的数据时钟281。PLL280在ADC输入时钟上操作以产生用于串行器262和数据端口270的数据时钟281。数据时钟频率对应于输出自每个串行器262i的比特的数据传输速率。数据时钟频率为ADC时钟频率乘以每采样周期中复用的比特的数量,在此例中为16。12个数据端口270以16倍于采样速率的数据传输速率传输被串行化的比特。
图10B示出针对图10A的配置的一个采样周期期间提供至数据端口270的比特排布。表290描绘在一个采样周期期间输出自16个ADC的样本比特的比特映射。对于具有四个复用器252和三个串行器262的配置,复用器252将每个信号样本的比特分组入4比特的块。复用器252i将第一比特块引导向串行器262i,将第二比特块引导向串行器262j并将第三比特块引导向串行器262k。线290i、290j和290k相应地代表对串行器262i、262j和262k的比特映射。以每个串行器262在一个采样周期期间输出16个映射的比特为条件,可以使用产生不同比特顺序的其他比特映射。在跨越数字接口220传输之后,接收的复用的比特依据逆向比特映射而重新排序以恢复原始信号样本序列。
图11为重新排序接收到的由端口集中所产生的复用的比特流,以恢复原始信号样本序列的框图。在端口集中模式中,接收器需要较少的输入数据端口310,从而节省电能并且降低复杂度。输入数据端口310为LVDS接收器。每个LVDS接收器接收由LVDS传输所产生的差分信号并且生成再现相应的接收的比特的逻辑状态。解串行器315将接收的串行比特转换为M个并行的复用的信号样本通道R1到RM。解复用器320将复用的信号样本的比特重新排列为信号样本的原始序列X1到XN的顺序。超声信号处理器130随后对再现的信号样本序列应用操作,比如波束形成或数字下变频转换。对于端口集中模式的其他备选方案在‘988申请中描述。
由压缩单元210i应用的压缩方法包括块浮点编码以及由块浮点编码所跟随的对信号样本的一阶或高阶导数的计算。Huffman或者其他类型的编码可以是块浮点编码的替代方法。
压缩单元210i的优选实施方式对来自ADC输出端121i的连续的信号样本组应用块浮点编码,每组具有N_GROUP个样本。N_GROUP个样本的最高指数被编码,并且N_GROUP个样本依照以下步骤编码。
对于第一组N_GROUP个样本:
1)比如通过计算每组N_GROUP个样本中的最高幅度的log2,为具有最高幅度的样本确定指数(以2为底)。这指示每个被编码的样本中的比特数量或者n_exp(0)。
2)使用S比特对第一组的指数n_exp(0)绝对编码,其中S为每样本比特的原始数量。
3)使用每样本n_exp(0)比特对N_GROUP个样本编码。
对于第i组N_GROUP个样本(i>0):
4)为具有最高幅度的样本确定第i(以2为底)指数,其指示在第i组中的每个被编码样本的比特数量,或者n_exp(i);
5)通过从n_exp(i-1)中减去n_exp(i)以确定第i差值,从而差分编码第i指数。使用相应的令牌编码第i差值,其中较短的令牌代表更常见的差值而较长的令牌代表较不常见的差值。
6)使用每样本n_exp(i)比特编码第i组N_GROUP个样本。
对于第一组样本,指数n_exp(0)被直接编码。例如,指数n_exp(0)可以按如下编码,其中S为每样本比特的原始数量:
a.0:n_exp(0)=0(全部4个样本值为零)
b.1:n_exp(0)=2(每样本2比特)
c.2:n_exp(0)=3(每样本3比特)
d.其他直到S-1:n_exp(0)=S(每样本S比特)
对于第i组,指数n_exp(i)使用前缀码差分编码,其中没有码字为另一个码字的前缀。优选的差分编码如下:
1.计算差分:e_diff=n_exp(i) -n_exp(i-1)
2.编码e_diff如下:
a.0:e_diff=e(i)-e(i-1)
b.101:e_diff=+1
c.110:e_diff=-1
d.1001:e_diff=+2
e.1110:e_diff=-2
f.其他.
或者,指数n_exp(i)可以被Huffman编码而不是被差分编码。
图12为在其中N_GROUP=4的块浮点编码器的框图。如在步骤1和步骤4中,指数计算器402为N_GROUP个样本确定比特最大指数n_exp。如在步骤2和步骤5中,指数令牌生成器404编码n_exp值。如在步骤3和步骤6中,尾数打包器406为N_GROUP个样本编码尾数。图13示出了选择用于编码的信号样本的n比特的示例。输入信号样本420由K比特代表。样本420的n_exp较低比特被选择用于编码。样本的符号比特被附加到选择的比特,而产生的比特序列代表尾数。返回图12,复用器408打包由N_GROUP个尾数所跟随的编码的指数令牌411以形成代表N_GROUP个压缩的样本的压缩组410。对于这个示例,压缩组410包括由四个打包的尾数序列412-0、412-1、412-2和412-3所跟随的指数令牌411。压缩单元210i连接连续的压缩组以在压缩器输出端211i形成压缩包的数据部分。N_GROUP的优选的大小为每组三个或四个样本。但是,也可以使用可变的组大小。
分别编码尾数和指数能够提供额外的压缩并且减轻压缩误差。连续指数的差值被计算并编码。指数缓慢的变化,所以有一些由零值串所分离的非零值。可以通过仅代表非零差值以及它们相应的位置而高效的编码指数差值。位置可以由相应的索引值或者相对于最后的非零差值的位置代表。指数差值的编码是无损的,这防止相对较大的误差。为了解码指数,指数值通过积分指数差值并解码相应的位置而重建。为了解码尾数,每个重建的尾数值被限制使得它不改变被解码的样本的相应的指数的值。对于n_exp的解码的指数,重建的尾数可以具有2n_exp-1的最大值。这防止了在尾数中的压缩误差对指数值的改变。
替代的块浮点编码方法包括降低代表尾数的比特数量以及如上所述的指数的差分编码。代表N_GROUP个样本的尾数的比特数量通过从每个尾数选择性地去除的一定数量的最不重要比特(LSB)(取决于组的n_exp值)而减少。图14示出了选择代表减少的尾数的比特的示例。n_exp如同在以上步骤1和步骤4中那样确定。比特的一定数量的n_LSB从最不重要的比特开始去除,而不是对尾数的全部n_exp较低比特编码。剩余的m_exp比特(m_exp=n_exp-n_LSB)被选择用以编码。n_LSB的值依据公式或者列表,取决于n_exp的值。图15为n_exp、n_LSB以及m_exp的示例值的列表。对于较大的n_exp值,更多LSB通过截尾或者舍入而去除,以形成具有m_exp比特的减少的尾数。例如,如果n_exp为12,那么去除3个LSB以便保留9个尾数比特用以打包N_GROUP减少的尾数。压缩器210能够存储n_exp、n_LSB和m_exp的值的查询表。或者,压缩器210可以将n_LSB和m_exp表示为n_exp的函数并且在需要时计算它们的值。图16为使用减少的尾数的块浮点编码的框图。对于每一组N_GROUP个样本,指数计算器402如以上所述确定最大指数n_exp。尾数比特计算器414使用查询表或公式确定在减少的尾数中的比特的数量m_exp。减少的尾数打包器416为N_GROUP个样本中的每一个选择m_exp比特。复用器408随后打包由减少的尾数419-0、419-1、419-2和419-3所跟随的指数令牌411以形成压缩组418。对于某些条件,不从N_GROUP个样本组中去除LSB。例如,当N_GROUP个样本中的一个或多个的样本的幅度小于可接受的最低值时,包括原始LSB的N_GROUP个尾数将会被打包。压缩的样本序列可以包括具有或不具有减少的尾数的压缩组。
压缩控制器为压缩单元210i提供压缩控制参数用于块浮点编码。可以有多种用于n_LSB、m_exp和n_exp的替代的查询表或公式。压缩控制参数包括N_GROUP以及用于n_LSB、m_exp和n_exp的替代的查询表或公式的选择参数。压缩控制参数对于所有的压缩单元210i可以是一致的。或者,压缩控制参数对于不同的压缩单元210i可以具有不同的值。压缩控制器能够响应于用户输入以选择压缩控制参数。
可以将压缩的样本插入压缩包的数据部分用以在数字接口220上传输。对应于由ADC 120i所数字化的接收的脉冲的压缩的样本序列可以安排在一个或多个压缩包中。或者,如关于图9所述,可以将来自多个ADC 120i的压缩的样本序列相组合以形成压缩包用以在给定的数据端口270i上传输。压缩包的头部部分包含包的标识信息。头部还可以包含代表在包中压缩的样本的压缩控制参数的控制数据。解压缩器240可以使用关于压缩控制参数的信息以配置解压缩操作。
为了解压缩,解压缩器240对压缩的样本序列应用块浮点解码。对于每组N_GROUP个压缩的样本,解压缩器240解码指数令牌以确定n_exp的值。被差分编码的指数被积分以形成n_exp的值。N_GROUP个尾数随后通过为来自压缩组410或418的每个尾数的比特解包以及将比特映射到解压缩的样本而重建。解压缩的样本可以由原始数量的每样本比特或者不同数量的每样本比特代表,取决于超声信号处理器130的下游处理要求。对于使用减少的尾数的块浮点编码器,解压缩器240还包括用于基于n_exp的解码值确定n_LSB的值的查询表或公式。减少的尾数的解包的比特由n_LSB比特附加,其可为零或抖动值,以近似原始样本值。
‘553专利描述用于某些带限信号的压缩和解压缩的算法。以下描述的某些替代的压缩方法为用于超声信号样本的‘533专利的算法的改进。
用于压缩超声信号样本的另一种替代方式为计算由编码所跟随的差分。计算超声信号样本的一阶或更高阶差分可以产生具有小于原始信号样本的幅度的差分样本。编码差分样本可以造成比编码样本本身更高的压缩。如上所述,计算ADC输出端121i的连续的样本的差分能够由差分样本的块浮点编码所跟随。在这种情况下,输入块浮点编码器的是差分样本,而不是ADC输出121i的信号样本。或者,可以对差分样本应用Huffman编码或者其他编码。
图17示出包括差分操作的压缩单元210i的框图。压缩单元210i接收来自ADC 120i的超声信号样本。压缩控制器340为每个压缩单元210i的差分运算器330i和编码器332i提供压缩控制参数。用于差分运算器330i的压缩控制参数可以选择一阶、二阶或者更高阶差分。差分运算器330i应用选择的差分阶以产生差分样本。压缩控制参数还可以选择绕过差分操作以便编码器332i编码信号样本而不是差分样本。用于编码器332i的压缩控制参数能够如以上所述,为块浮点编码器,或者为Huffman编码器或者另一个编码器指示参数。压缩控制参数对于不同的压缩单元210i可以是相同的或者不同的。
图18为解压缩器240的框图。根据系统架构,解压缩器接收来自数字接口220、捕获存储器230或者波束形成器存储器162的压缩的样本。解码器352逆转编码器332i的操作以形成解码的样本。例如,解码器352执行块浮点解码、Huffman解码或者其他解码。积分运算器354添加解码的差分样本,以逆转为压缩所执行的一阶或更高阶差分。如果未执行差分用于压缩,积分运算器354将被绕过。解压缩控制器350将控制参数提供给解码器352和积分运算器354。解压缩控制器350能够从压缩的数据包的头部提取控制数据用以为解压缩操作确定控制参数。
用于压缩的另一种替代方式基于中心频率和采样速率对超声信号样本应用算术运算。图19给出了示出在用于压缩具有不同的中心频率的信号样本的替代方式之下的原理的示例。以对应于图19中的标记为“带1”的行的基带信号的示例开始,中心频率接近于DC(0Hz)并且连续的样本之间的相位增加小于10度。第一相量图710示出,因为连续的样本之间的相位变化很小,所以连续的样本的差分的幅度相较于样本本身的幅度相对较小。第一示例序列712对应于带1基带信号的样本。由于连续的样本之间的差分相对于样本幅度较小,计算一阶或更高阶差分,或者差分编码产生具有小于原始样本的数据宽度的差分样本。对于基带(带1)示例,使用关于图17描述的差分编码压缩是高效的。
图19还给出其中中心频率在DC之上,但是低于奈奎斯特(Nyquist)频率的fs/2的采样信号的示例。对于带2,中心频率为接近fs/6而连续样本之间的相位增加为大约60度。第二相量图720示出,由180度,或者三个采样间隔所分隔的样本对,如样本对(720-0,720-3)、(720-1,720-4)和(720-2,720-5)所示,具有相似的幅度以及相反的极性。反相对中的一个样本[或者乘以(-1)]提供对中的另一个样本的接近估计。第二示例序列722还示出由三个采样间隔所分隔的样本具有相似的幅度以及相反的正负号。例如,样本722-0的值为32767而样本722-3的值为-32756。对于带2,在由三个采样间隔所分隔的样本对上的操作产生具有较小数据宽度的变更的样本。将对中的样本相加在一起的操作产生能够被更高效编码的具有较小数据宽度的变更的样本。
对于图19中的带3的示例,中心频率为接近fs/4而连续样本之间的相位增加为大约90度。第三相量图730示出,由180度,或者两个采样间隔所分隔的样本,具有相似的幅度和相反的极性。第三示例序列732还示出,每隔一个样本都具有相似的幅度和相反的极性。对于带3,将每隔一个样本相加将会产生具有能够以比原始样本更高的效率编码的较小数据宽度的变更的样本。
对于图19中的带4的示例,中心频率为接近fs/3而连续样本之间的相位增加为大约120度。第四相量图740示出,由360度,或者三个采样间隔所分隔的样本,将具有相似的幅度。第四示例序列742示出,每三个样本具有相似的幅度。在这种情况中,在由3个采样间隔所分隔的样本之间形成差分将会给出能够以比原始样本更高的效率编码的具有较小数据宽度的变更的样本。
对于图19中的带5的示例,中心频率为fs/2而连续样本之间的相位增加为大约180度。第五相量图750示出,由180度,或者一个采样间隔所分隔的样本,将具有相似的幅度以及相反的极性。第五示例序列752示出连续样本具有相似的幅度和相反的极性。在这种情况中,将两个连续样本相加将形成能够以比原始样本更高的效率编码的具有较小数据宽度的变更的样本。
上述对于图19描述的示例示出,幅度降低能够通过根据采样速率与中心频率的比率,在由1、2或3个采样间隔所分隔的信号样本上,执行比如加(或者反相随后减)或减(或者反相随后加)的操作而实现。产生的变更的样本随后被编码以形成压缩的样本。类似的操作可以根据中心频率与采样速率的比率,应用到由四个或更多个采样间隔所分隔的样本,以产生具有比原始样本较小的数据宽度的变更的样本。
图20为基于信号样本的中心频率的压缩算法的框图。ADC 120i将超声信号样本提供给重新排序解复用器810。重新排序解复用器810选择信号样本,从而使选择的样本依据压缩控制参数852由合适数量的采样间隔所分隔,以形成解复用器输出812。算术运算器830依据压缩控制参数856在解复用器输出样本对812上执行加或减操作,以形成变更的样本832。算术运算器830还可以被配置以在解复用器输出样本812上执行更高阶差分。编码器840对变更的样本832进行编码,以形成压缩的信号样本。编码器840应用块浮点编码、Huffman编码或者其他编码以形成压缩的样本。针对块浮点编码,变更的样本832,而不是ADC输出端121i的信号样本,被输入块浮点编码器。
压缩控制器860基于采样速率与信号样本的中心频率的比率,将控制参数提供给压缩器元件。重新排序解复用器810和算术运算器830相应地响应于压缩控制参数852和856,以执行合适的操作。图21示出基于中心频率产生变更的样本832的操作。第一列871给出对于这个示例的可能的中心频率。第二列872给出对每个中心频率的相应的频带指示符。指示符可作为压缩控制852和856的参数使用。第三列873给出可依据压缩控制参数852产生的在重新排序解复用器输出812的样本x(i)和x(j)的不同分隔。第四列874示出依据压缩控制参数856选择加或减的算术运算的结果。当逆转器在延迟的样本“上”时,x(i-j)被减去。第五列875示出产生变更的样本832,或y(i)的算术运算器830的数学结果。压缩控制器860还提供编码器840的控制。压缩控制参数858能够指示块浮点编码或者其他编码技术的参数。
图22给出如关于图15和16描述的对于不同中心频率计算出的图19的示例的信号样本x(i)与x(i-j)的和或者差。信号样本的示例序列与图19中的那些相同。在示例序列912和942中的DIFF行以及在示例序列922、932和952中的SUM行中的样本具有基本上低于相应的信号样本,或x(i)的幅度。DIFF样本和SUM样本为输入至图20中的编码器840的变更的样本932的示例。
图23为对于关于图20描述的压缩器210的解压缩器240所执行的操作的框图。压缩的样本由解码器910根据系统架构接收自数字接口220、捕获存储器230或者波束形成器存储器162。解码器910在压缩的数据上解包并且执行解码操作,譬如块浮点解码,以形成解码的变更的样本。逆转算术运算器920执行算术运算器830的逆转运算以从解码的变更的样本重建信号样本。复用器930将原始次序恢复到解压缩的信号样本以重建超声信号样本序列。解压缩控制器940将控制参数提供给解码器910、逆转运算器920以及复用器930。解压缩控制器940能够从压缩的数据包的头部中提取控制数据以确定用于解压缩操作的控制参数。
压缩器210的实施方式应用能够实时的,或者以至少与采样速率一样快的速率压缩输出自ADC的样本的简单操作。差分运算器330i(图17)包括一个或多个减法器。块浮点编码(图13和16)使用比较器、减法器以及查询表。或者,Huffman编码使用查询表为值分配编码。关于图20描述的压缩操作包括解复用、加以及减。用于解压缩器240的实施方式应用简单的操作以解压缩被压缩的样本。解压缩器240包括用于块浮点解码的查询表和加法器。积分运算器354(图18)包括一个或多个用于积分解码的样本的加法器。图23中的解压缩器240的操作包括加、减和复用。
在超声系统中本发明的优选实施方式将压缩器210和ADC库120集成在单独的专用集成电路(application specific integratedcircuit,简称ASIC)设备中。参考图5的框图,ADC库120与压缩器210集成在具有模拟输入端和数字输出端的混合信号集成电路设备中。N(多个)独立的ADC 120i并行地将N个输入模拟超声信号转换为N个数字超声信号。ADC 120i可以由包含多个快闪式转换器级的流水线数据转换器或者其他ADC架构实施。用于ADC的知识产权(intellectual property,简称IP)核商业上可得到用于ASIC实施。ADC输出通道121i与实施压缩器210的数字逻辑耦合。压缩器210的优选实施方式包括多个并行的压缩核,其中每个压缩核与ADC120i中的一个耦合并且实施一个压缩单元210i的压缩操作。或者,一个压缩核可以实施多个压缩单元210i以压缩来自一个以上ADC120i的信号样本。对于这种替代方式,压缩核包括缓冲器以存储来自不同ADC 120i的信号样本,直到它们被处理。压缩操作还可以在现场可编程门阵列(FPGA)中实施。如图8、9和10中所示,压缩的样本可以通过LVDS端口270i输出到数字接口220。LVDS接口的IP核商业上可得到用于ASIC和FPGA实施。替代实施包括在耦合到并行ADC设备的输出端的独立的设备中的压缩操作。压缩操作可以在ASIC、FPGA或者比如数字信号处理器(DSP)、微处理器、微控制器、多核CPU(比如IBM Cell),或者图形处理单元(GPU;比如Nvidia GeForce)等可编程处理器中实施。
根据超声系统架构,解压缩器240可被并入超声信号处理器130的同一设备或者作为与其不同的设备。解压缩操作可以在ASIC或者FPGA中实施。或者,解压缩操作可以在可由比如DSP、微处理器、微控制器、CPU或者GPU等可编程处理器执行的软件或固件程序中实施。解压缩器240的优选实施为具有可由GPU执行的解压缩操作的指令的软件程序。GPU也可以被编程以在解压缩的样本上实施超声信号处理器130的至少一部分操作,比如波束形成操作。或者,解压缩的样本可以传输到另一个可编程处理器,比如CPU,用于超声信号处理操作。
图24为GPU中的解压缩器240以及其他用以生成超声图像的操作的实施的框图。当前的GPU架构包括多个被优化用于并行计算的处理核。例如,Nvidia GeForce GTS 150GPU包括128个处理核。Nvidia的“CUDA”(Compute Unified Device Architecture,计算统一设备架构)是包括用于在GPU的处理核上实施并行算法的C语言的扩展的应用编程接口(application programming interface,简称API),其在标题为“Getting Started with CUDA”的文件(Ruetsch和Oster,Nvidia Corp.,2008)中描述。以下所述的比如OpenCL和Larrabee之类的替代的GPU和编程方法可以提供实施平台。对于图24中描绘的实施,GPU设备1000可被编程以执行解压缩器240、超声信号处理器130(波束形成、B-模式处理以及多普勒处理)、扫描转换器140以及图像处理器150的操作。GPU设备1000可以包括可由并行处理单元访问的动态随机访问存储器(dynamic random accessmemory,简称DRAM)1002。DRAM1002可以存储由GPU的其他处理操作所产生的数据和压缩的和/或解压缩的样本。系统控制器1010提供用于从接收的数据产生超声图像的任务的协调并响应于用户指令。CPU1012可以实施支持解压缩的操作,比如从压缩包的头部解码压缩控制参数并将它们提供给GPU设备1000用于解压缩操作的配置。DRAM1016可以存储接收自数字接口220的压缩的样本以及其他CPU操作所需的数据。通信控制器1014将接收自数字接口220的压缩包引导至DRAM1002或DRAM1016并且管理系统控制器1010和GPU设备1000之间的数据交换。
在优选的系统架构中,系统控制器1010可以在具有用于显示的屏幕160的计算机的主板中实施。GPU设备1000可以在包括DRAM1002,与系统控制器1010通过PCIe(Peripheral ComponentInterconnect Express,外围组件互连高速)背板链路而通信的图形卡中实施。或者,GPU设备1000可以在安装在主板上的IC中实施。在其中ADC库120和压缩器210安装在数据采集卡中的系统架构中,数字接口220可以由PCIe背板链路所实施。
应该认识到,依照摩尔定律(Moore’s Law),更高水平的集成将提供更加紧凑的设备,从而系统控制器1010和GPU设备1000可以在单一的IC中实施。例如,Intel Corporation正在开发包括多个以向量处理单元增强的x86 CPU核的实例的多核IC架构。此架构被称为Larrabee,在标题为“Larrabee:A Many-Core x86 Architecture forVisual Computing”的文件(Seiler等人,ACM Transactions on Graphics,Vol.27,No.3,Article 18,August 2008)中描述。Larrabee架构支持需要包括图形处理的并行处理的应用。“Larrabee Native”编程模型包括用于并行编程和向量化的C/C++和API。
图25为CPU/GPU设备中的解压缩器240以及其他用以产生超声图像的操作的实施的框图。CPU/GPU设备1020可以实施解压缩器240、超声信号处理器130、扫描转换器140和图像处理器150的操作。CPU/GPU控制器1022协调在接收自数字接口220的压缩的样本上的处理操作并响应用户输入。CPU/GPU设备1020可以由Larrabee平台或者其他具有集成的CPU和GPU功能的可编程设备所实施。
在其中ADC库120和压缩器210被安置于换能器头中的系统架构中,数字接口220可以为有线或无线通信链路。对于有线通信链路,数字接口可以由PCIe电缆链路或者光纤链路所实施。对于无线通信链路,数字接口可以提供数字调制以及经由射频信道的压缩包的传输以及接收的压缩包的数字解调。无线链路可以遵守无线通信协议,比如WiFi(IEEE 802.11)或者UWB(ultra-wideband,超宽带)格式。
虽然示出和描述了本发明的优选实施方式,应该清楚,本发明不仅仅限定在这些实施方式。众多修改,变化,变异,替代和等效将会在不背离在权利要求中所描述的本发明的精髓和范围的前提下对本领域中的技术人员来说是明显的。