CN101933814B - 超声系统中的波束形成后压缩 - Google Patents
超声系统中的波束形成后压缩 Download PDFInfo
- Publication number
- CN101933814B CN101933814B CN2010101854538A CN201010185453A CN101933814B CN 101933814 B CN101933814 B CN 101933814B CN 2010101854538 A CN2010101854538 A CN 2010101854538A CN 201010185453 A CN201010185453 A CN 201010185453A CN 101933814 B CN101933814 B CN 101933814B
- Authority
- CN
- China
- Prior art keywords
- wave beam
- sample
- compression
- forms
- group
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52034—Data rate converters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
- G10K11/341—Circuits therefor
- G10K11/346—Circuits therefor using phase variation
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
- H03M7/40—Conversion to or from variable length codes, e.g. Shannon-Fano code, Huffman code, Morse code
Abstract
在对接收到的超声信号样本应用波束形成器以形成由波束形成样本的阵列代表的一个或多个波束的超声成像系统中,一种方法和装置独立于其他阵列而压缩波束形成样本的每个阵列以形成压缩的波束。多个模拟到数字转换器采样由换能器阵列所产生的多个模拟超声信号以将多个超声信号样本流提供给波束形成器。压缩的波束经由数字接口传输至信号处理器。在信号处理器处,对压缩的波束进行解压缩以形成解压缩的波束。信号处理器进一步处理解压缩的波束用于供诊断成像(比如B-模式和多普勒成像)以及扫描转换,以准备生成的用以显示的超声图像。此摘要不限定如权利要求中所描述的本发明的范围。
Description
技术领域
本发明涉及在超声成像系统中由应用在接收的超声信号样本的接收波束形成器所产生的波束形成样本的压缩,特别是在用于图像形成的处理之前独立地压缩每个波束的波束形成样本以及解压缩。
背景技术
医疗超声系统通过从由临床医生放置在受试对象上的换能器发送超声波束,扫描受试对象的内部解剖结构。超声波在具有不同声阻抗的内部组织的界面上反射,产生回声。换能器接收回声并将它们转换为电超声信号。超声系统将一序列的处理步骤应用到超声信号以产生在控制台上显示的图像或一系列图像,用以由临床医生分析。基于接收的回声的强度形成的图像称为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 PhotographicExperts Group开发的标准图像压缩方法)。在变换域中对于比如过滤之类的各种图像处理功能的额外的操作在变换域中比在空间域中在计算上更加高效。例如,空间域中的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 AnnualIEEEEngineering 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专利)描述用于特定带限信号的压缩与解压缩的算法。提交于2000年6月2日的名称为“Ultrasound Signal Compression”的共同拥有和共同未决的美国专利申请序列号12/477062(‘062申请)描述了在波束形成操作之前对输出自模拟到数字转换器的超声信号样本的压缩和解压缩。本申请针对于对产生自应用到超声信号样本的波束形成操作的超声波束形成样本的压缩。
存在对超声成像系统的部件之间的超声波束形成样本的高效数据传输的需求。存在对超声波束形成样本在计算上高效的数据压缩以在对系统复杂度产生最小影响的情况下改进数据传输。
发明内容
本发明的实施方式考虑到上述的传统问题而做出。本发明提供了用于压缩代表一个或多个由超声成像系统中的接收波束形成器所产生的波束形成样本的装置和方法。将接收波束形成器应用到多个在采样窗口期间接收的超声信号样本序列以形成一个或多个波束,每个波束由波束形成样本阵列代表。超声信号样本序列通过由换能器元件在采样窗口期间输出的模拟超声信号的模拟到数字转换而产生。本发明的方法和装置提供波束的压缩以形成压缩的波束,包括压缩代表特定波束的阵列中的波束形成样本以形成相应的压缩的波束的压缩的波束形成样本。阵列中特定波束形成样本的压缩部分地取决于阵列中至少一个其他波束形成样本的特征。每个波束独立于对应于采样窗口的另一波束而压缩。压缩的波束跨越数字接口传输到信号处理器。在信号处理器处,压缩的波束在用于图像形成的处理之前被解压缩。
在本发明的装置的实施方式中,压缩器包括一个或多个压缩单元。相应的压缩单元独立于对应于采样窗口的其他波束,对输出自接收波束形成器的相应的波束应用压缩操作。相应的压缩单元配置以部分基于相应的波束的至少一个其他波束形成样本的特征来压缩特定的波束形成样本。
在本发明的另一个方面,压缩操作包括对特定波束的波束形成样本组应用块浮点编码以形成压缩的波束形成样本组。组复用器将对应于两个或更多个压缩的波束的组相交错以形成复用的序列用以跨越数字接口传输。对接收到的复用的序列的压缩组进行解压缩以形成解压缩的波束形成样本组。组解复用器重新布置解压缩的组以恢复相应的解压缩的波束的原始组次序。在本发明的另一种实施方式中,压缩操作包括计算相应波束的波束形成样本之间的差分,并且跟随以差分样本的块浮点编码,以形成压缩的波束。
波束内的波束形成样本之间的相关性通过部分根据同一波束中的另一波束形成样本的特征来压缩特定的波束形成样本而利用。对于包括块浮点编码的实施方式,特征为波束内的波束形成样本组中具有最大幅度的波束形成样本的指数值。块浮点编码基于指数值以及波束形成样本本身的值来代表组中的每个波束形成样本。对于包括计算波束中的波束形成样本对之间的差分的实施方式,特征为该对中的每个波束形成样本的值。独立于采样窗口中的其他波束来压缩每个波束允许对于产生压缩的波束的高效计算和低延迟。
压缩超声成像系统中波束形成样本的优点包括降低了用于在数字接口上向信号处理器传输压缩的波束形成样本所需的带宽。对于其中将接收波束形成器和压缩器安置于换能器头之中的系统架构,数字接口可以是有线或无线的通信链路。对于有线通信链路,数字接口可以通过较低成本的电缆部件,比如PCIe(Peripheral ComponentInterconnect Express,外围组件互连高速)电缆链路或者光纤链路而实施。对于无线通信链路,降低的所需带宽可以降低实施无线链路的发射器和接收器的成本和复杂度。对于其中数字接口的带宽固定的系统架构,波束形成样本的压缩允许压缩形式中的更多波束传输至信号处理器用于图像形成。信号处理器可以使用更多数量的波束形成样本以产生具有提高的分辨率或更大视场的超声图像。
本发明的另一个优点包括压缩的波束形成样本的高效存储。对于其中波束形成样本存储于存储器或其他存储介质中的系统架构,用于存储压缩的波束形成样本所需的存储容量相较于未压缩的波束形成样本是减少的。
本发明的另一个优点包括在计算上高效的压缩和解压缩。压缩和解压缩的实现较不复杂,降低了系统资源上的负担并且降低成本。本发明的其他方面和优点可以通过参阅以下附图、详细描述以及权利要求而看出。
附图说明
图1为依据现有技术的典型医疗超声系统的示例的框图。
图2为依据现有技术的超声信号处理器130的框图。
图3为依据现有技术的在波束形成之前的超声信号样本的绘图。
图4为依据现有技术的对应于一个波束形成器输出通道的波束形成超声信号的同相位样本的绘图。
图5为包括输出自接收波束形成器的波束形成样本的压缩的超声成像系统的框图。
图6为包括复用压缩的波束形成样本用以跨越较少数据端口传输的超声成像系统的框图。
图7为在其中N_GROUP=4的块浮点编码器的框图。
图8示出了选择用于编码的样本的n比特的示例。
图9示出了选择用以代表减少的尾数的比特的示例。
图10为n_exp、n_LSB以及m_exp的示例值的表。
图11为使用减少的尾数的块浮点编码的框图。
图12示出了复用压缩的波束的集合的压缩组的示例。
图13示出了对应于图12的复用示例的解复用示例。
图14示出了包括差分操作的压缩单元的框图。
图15为与关于图14描述的压缩操作相对应的解压缩操作的框图。
图16为解复用解码的差分样本组的框图。
图17给出了示出在用于压缩具有不同的中心频率的波束形成样本的替代方式之下的原理的示例。
图18为基于波束形成样本的中心频率的压缩算法的框图。
图19示出了基于中心频率产生修改的样本832的操作。
图20给出图17的示例的样本x(i)与x(i-j)的和或者差。
图21为对于关于图18描述的压缩方法的由解压缩器240所执行的操作的框图。
图22为解复用解码的修改的样本组的框图。
图23为GPU中的解压缩器以及其他用以产生超声图像的操作的实施的框图。
图24为CPU/GPU设备中的解压缩器以及其他用以产生超声图像的操作的实施的框图。
具体实施方式
图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的多个模拟超声信号耦合到AFE 114。采样窗口对应于时间间隔,在其中接收的回声代表来自在受试对象中期望的深度范围内的散射位置的反射。控制器102依据用户输入或者扫描协议设置采样窗口并将定时控制信息提供给发送/接收切换器112。发送/接收切换器112在采样窗口期间输出并行的多个模拟超声信号。AFE 114将多个模拟超声信号进行放大并过滤以对模拟到数字转换进行准备。对于每个模拟信号通道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对超声信号样本应用延迟模式(取决于所计算的波束的方向)。接收波束形成器160产生对应于视场中的多个方向的多个波束形成样本的1-D阵列。接收波束形成器160包括波束形成器(BF)存储器162以存储用于波束形成(BF)计算器164的延迟和加法操作的超声信号样本。针对同一接收的脉冲,BF计算器164能够使用取回自BF存储器162的超声信号样本形成多个波束。BF计算器164还可以对超声信号样本应用差值操作,以提高计算的波束的相位分辨率。接收波束形成器160还可以在加法操作之前对信号样本应用加权函数,以实施空间窗口功能或者变迹。将为每个角度计算的波束形成样本提供给相应的波束形成器输出通道161i。接收波束形成器160通常具有少于输入通道的,包含ADC输出通道121i的输出通道161i。在这种配置中,波束形成样本具有RF中心频率。数字下变频转换器(digitaldown converter,简称DDC)170将波束形成样本解调至基带,以生成每个波束的复基带I和Q样本。作为DDC 170的替代或补充,可以对在以期望的频率为中心的频带处的波束形成样本应用带通滤波器,或者DDC 170可以将波束形成样本解调为中间频率(intermediatefrequency,简称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为依据优选实施方式的包括压缩由接收波束形成器160所产生的波束形成样本的超声成像系统的框图。压缩器210包括多个压缩单元210i,每个压缩单元对相应的波束形成器输出通道161i应用压缩操作。压缩单元210i对代表特定波束的波束形成样本的1-D阵列应用压缩操作以在压缩器输出211i处产生相应的压缩的波束形成样本序列。特定波束的波束形成样本独立于代表其他波束的波束形成样本而压缩。备选地,压缩单元210i可以压缩代表接收自多个波束形成器输出通道161i的多个波束的波束形成样本。在这种备选方式中,压缩单元210i独立地对代表不同波束的波束形成样本应用压缩操作。压缩的波束形成样本跨越数字接口220传输至诊断信息处理器180。解压缩器240解压缩接收的压缩的数据以重建波束形成样本的1-D阵列用于由诊断信息处理器180进行进一步信号处理操作。用于压缩的波束形成样本的传输所需的数字接口220的数据传输带宽相较于用于未压缩的波束形成样本的传输所需的带宽是减小的。
在处理序列中的某个点处,将处理的超声信号被下变频转换到IF或基带。下变频转换可以在波束形成之后和压缩之前,或者在解压缩之后,由ADC库120的模拟到数字转换之前或之后应用。当在压缩之前应用下变频转换时,输入至压缩器210的波束形成样本将具有在中间频率(IF)或在零频率(基带)处的中心频率。当在解压缩之后应用下变频转换时,输入至压缩器210的波束形成样本将具有RF中心频率。当前描述假设波束形成样本为实值样本。备选地,如果ADC 120i执行正交采样或者如果在压缩前应用正交下变频转换,那么波束形成样本将具有I和Q分量。在这种情况中,压缩器210将独立地处理波束形成样本的I和Q分量。
由于压缩的波束形成样本具有更少的每样本比特,用于连接到数字接口220的数据端口可以具有过剩带宽。可以复用压缩的波束形成样本以形成较少的用于通过数字接口220传输的数据通道。图6为包括复用压缩的波束形成样本用以跨越较少数据端口传输的超声系统的框图。接收波束形成器160从相应ADC输出121i接收N个超声信号样本流X1到XN,并且在波束形成器输出通道161i处产生J个波束B1到BJ。压缩单元210i产生J个相应的压缩的波束形成样本流C1到CJ。具有降低的比特率的压缩的波束形成样本提供给数据端口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对。当在压缩器输出Ci处的压缩的波束的比特率低于端口的最高数据传输速率时,数据端口270i具有过剩带宽。过剩带宽可以通过组合多个压缩的波束用以通过给定的数据端口270i传输而利用。复用器250将J个压缩的波束形成样本序列相组合以形成M个复用的序列Dj,其中M<J,用以在M个数据端口270i上传输。可以结合的压缩的波束的数量受限于数据传输端口270i的带宽。串行器260将复用的比特序列提供给相应的数据端口270i。
复用压缩的波束形成样本以使用M个数据端口的优点包括使用较少的物理数据端口,这继而减少了数据端口的功耗和连接。在其中压缩的波束形成样本的比特率被固定的应用中,复用器250具有对于固定数量的物理数据端口270i的固定数量的复用器输出251i。备选地,灵活架构可以依据压缩的波束形成样本的比特率来支持固定数量的物理数据端口270i之中的可变数量的活动数据端口。对于灵活架构,压缩器210依据用户可选的压缩比参数,以不同的比特率提供压缩的波束形成样本。压缩控制器(未在图5中示出)为压缩单元210i提供压缩控制参数,以便操作产生具有对应于期望的压缩比的比特率的压缩的波束形成样本。压缩控制器为复用器250提供复用控制参数以指示需要产生的复用的序列的数量M。压缩控制器还可以响应于用户输入而断电不活动数据端口,从而进一步节能。
可被结合用于一个数据端口的压缩的波束的数量由数据端口270i的带宽所限制。例如,对于其中每个波束具有每波束形成样本16比特以及每秒50兆样本的采样速率的J=8波束Bi,每个波束形成器输出通道161i处的比特率为800Mbps。假设数据端口270i具有高达800Mbps的数据传输速率而压缩器210产生2∶1的压缩比。压缩器输出211i处的压缩的波束Ci具有400Mbps的比特率。在这种情况中,一个数据端口270i具有足够的带宽以800Mbps的数据传输速率传输2个压缩的波束Ci。复用器250结合来自2个压缩的波束集合Ci和Ci+1的压缩的波束形成样本,以在一个复用器输出251i形成相应的复用的压缩的样本序列Dj。在这个示例中,提供给相应的活动数据端口的复用的序列数量为M=J/2=4。串行器260向相应的数据端口270i提供比特序列。在数字接口220上传输之后,解复用器可以接收M个复用的序列并将它们解复用以恢复提供给解压缩器240的J个压缩的波束。备选地,复用器250可以交错具有固定每组样本数量的压缩的波束形成样本组。对于这种备选方式,如以下关于图12和图13所描述的,解压缩器240可以在解复用之前解压缩压缩的波束形成样本组。
由压缩单元210i应用的压缩方法独立于其他波束而对每个波束形成器输出通道161i的波束形成样本进行操作。应用到波束形成样本的压缩方法取决于波束中其他样本的值。压缩方法包括块浮点编码以及由块浮点编码所跟随的计算波束形成样本的一阶或更高阶导数。Huffman或者其他编码类型可以是块浮点编码的备选方式。
压缩单元210i的优选实施方式对来自波束形成器输出通道161i的连续的波束形成样本组应用块浮点编码,其中每组具有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编码而不是差分编码。
图7为其中N GROUP=4的块浮点编码器的框图。如在步骤1和步骤4中,指数计算器402为N_GROUP个样本以比特、n_exp的形式确定最大指数。如在步骤2和步骤5中,指数令牌生成器404编码n_exp值。如在步骤3和步骤6中,尾数打包器406编码N_GROUP个样本尾数。图8示出了选择用于编码的样本的n比特的示例。输入样本420由K比特代表。样本420的n_exp较低比特被选择用于编码。样本的符号比特被附加到选择的比特,而产生的比特序列代表尾数。返回图7,复用器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值)而减少。图9示出了选择代表减少的尾数的比特的示例。n_exp如同在以上步骤1和步骤4中那样确定。比特的一定数量的n_LSB从最不重要的比特开始去除,而不是对尾数的全部n_exp较低比特编码。将剩余的m_exp比特(m_exp=n_exp-n_LSB)选择用以编码。n_LSB的值依据公式或者表,取决于n_exp的值。图10为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的函数,并且在需要时计算它们的值。图11为使用减少的尾数进行块浮点编码的框图。对于每一组N_GROUP个样本,指数计算器402如以上所述确定最大指数n_exp。尾数比特计算器414使用查询表或公式确定减少的尾数中的比特的数量m_exp。减少尾数打包器416为N_GROUP个样本中的每一个选择m_exp比特。复用器408随后打包由减少的尾数419-0、419-1、419-2和419-3所跟随的指数令牌以形成压缩组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上传输。对应于采样窗口的压缩的波束的压缩的波束形成样本可以布置在一个或多个压缩包中。备选地,如关于图6所述,可以将对应于多个压缩的波束的压缩的波束形成样本序列相结合以形成压缩包用以在给定的数据端口270i上传输。压缩包的报头部分包含包的识别信息。报头还可以包含代表包中的压缩的波束形成样本的压缩控制参数的控制数据。解压缩器240可以使用关于压缩控制参数的信息以配置解压缩操作。
图6的复用器250的优选实施方式复用N_GROUP个压缩的波束形成样本的组。复用器250通过交错J个压缩的波束的压缩组而产生M个复用的序列。将J个压缩的波束分为压缩的波束的M个集合。对于每个集合,复用器250交错集合中的压缩的波束的压缩组以形成相应的复用的序列Dm。复用的序列Dm是来自由复用器250布置在组次序中的相应集合的压缩组序列。复用器250可以作为M个并行的复用器实施,每一个复用器接收相应集合的压缩的波束。
图12示出了复用压缩的波束C1和C2的集合的压缩组的示例。组复用器252交错来自压缩的波束C1的压缩组G1j和来自压缩的波束C2的压缩组G2j以形成复用的序列D1。在这个示例中,复用的序列D1的组的次序在压缩的波束之间交替以形成序列G1j、G2j、G1(j+1)、G2(j+1),等。压缩的组Gij可以代表对应于图7的块浮点压缩组410、图11的压缩组418或者由比如Huffman编码的另一种编码方法所产生的压缩组的打包的比特。压缩组具有每组相同的压缩的波束形成样本数量,即,N_GROUP个样本。然而,由于每个压缩的样本的比特数量可能变动,所以每个压缩组的比特数量也可能变动。对于上述的块浮点编码方法,每个压缩组的比特数量取决于n_exp的值。
为了解压缩,解压缩器240解码并且解包对应于每个压缩的波束的压缩组的压缩的波束形成样本。对于每组N_GROUP个压缩的波束形成样本,解压缩器240解码指数令牌以确定n_exp的值。对差分编码的指数进行积分以确定n_exp的值。N_GROUP个尾数随后通过将来自压缩组410或418的每个尾数的比特解包以及将比特映射到相应的解压缩的波束形成样本以形成解压缩组而重建。取决于下游处理要求,解压缩的波束形成样本可以由每波束形成样本的原始比特数量或者每波束形成样本的不同比特数量代表。对于使用减少的尾数的块浮点编码器,解压缩器240还包括用于基于n_exp的解码值确定n_LSB的值的查询表或公式。减少的尾数的解包的比特通过n_LSB比特附加,其可为零或抖动值,以接近原始波束形成样本值。根据特定的压缩的波束计算的解压缩组序列形成相应的解压缩的波束。
为了复用的压缩组序列的解压缩,解压缩器240可以在解复用之前应用到复用的序列。参考图13,复用的序列D1的每个压缩组Gij包括提供用以解压缩那个组的N_GROUP个压缩的波束形成样本的必要信息的指数令牌。对于差分编码的指数令牌,解压缩器使用复用的序列D1的组次序以选择对应于用于积分的相同的压缩的波束的指数令牌以确定n_exp的相应的值。组次序可以是对于复用和解复用的固定参数,在这种情况下其由解压缩器240和组解复用器254所知。备选地,组次序可以是可调参数,在这种情况下它可以随同其他压缩参数包括在压缩包的报头部分。如前所述,压缩组可以具有不同的比特数量。然而,一旦组被解压缩,每一组的N_GROUP个解压缩的波束形成样本将具有相同的每样本比特数量。解压缩组的组次序与压缩组的相同。组解复用器254将解压缩组重新排序到它们的原始次序以形成解压缩的波束的集合,每个具有按照原始的连续次序的解压缩的波束形成样本。
图13示出了对应于图12的复用示例的解复用示例。解压缩器240对复用的序列D1应用块浮点解码以形成解压缩组的序列241。每组Rij中的解压缩的波束形成样本具有相同的每样本比特数量。组解复用器254恢复解压缩的波束形成样本的原始次序以形成解压缩的波束253和255的集合。组解复用器254使用组次序,即,R1j、R2j R1(j+1)、R2(j+2)等,以将解压缩的波束形成样本组分离到它们相应的解压缩的波束253和255。由于每组的样本数量为N_GROUP,所以组解复用器254将每组Rij的N_GROUP个解压缩的波束形成样本附加到相应的解压缩的波束。
用于波束形成样本的压缩和解压缩的备选方法包括在产生具有较低幅度的修改的样本之前,对每个波束内的波束形成样本执行计算,并因此需要较少的比特来代表它们。‘553专利描述了用于某些带限信号的压缩和解压缩的算法。以下描述的某些备选的压缩方法为用于波束形成样本的‘533专利的算法的修改。
用于波束形成样本的压缩的备选方法包括跟随以编码而独立于其他波束计算每个波束中的波束形成样本之间的差分。计算波束形成样本的一阶或更高阶差分可产生具有小于原始波束形成样本的幅度的差分样本。编码差分样本可产生比编码波束形成样本本身更大的压缩。
图14示出了包括差分操作的压缩单元210i的框图。压缩单元210i接收来自波束形成器输出通道161i的波束形成样本。压缩控制器340为每个压缩单元210i的差分运算器330i和编码器332i提供压缩控制参数。用于差分运算器330i的压缩控制参数可以选择一阶、二阶或者更高阶差分。差分运算器330i应用选择的差分阶以产生差分样本。压缩控制参数还可以选择绕过差分操作以便编码器332i编码波束形成样本而不是差分样本。如上所述,编码器332i可以应用差分样本的块浮点编码。在这种情况下,将差分样本提供给图7或图11的块浮点编码器的输入401。备选地,编码器332i可以对差分样本应用Huffman编码或者其他编码。用于编码器332i的压缩控制参数可以如以上所述指示块浮点编码器的参数,或者指示Huffman编码器或者另一编码器的参数。压缩控制参数对于不同的压缩单元210i可以是相同的或者不同的。
如以上关于图6和图12所述,对应于不同的压缩的波束的压缩器输出221i可以在数字接口220上传输之前复用至较少的输出通道。编码器332i可以对差分样本应用块浮点编码以产生压缩组Gij,其在图12中代表包括压缩的波束C1和C2的N_GROUP个编码的差分样本的组。备选地,编码器332i可以对差分样本应用Huffman编码从而使Gij代表具有每个压缩组N_GROUP个压缩的波束形成样本的Huffman编码的差分样本组。组复用器252结合压缩的波束C1和C2的组以形成复用的序列D1。
图15为对应于关于图14所述的压缩操作的解压缩操作的框图。解压缩器240接收来自数字接口220的压缩的波束。解码器352逆转编码器332i的操作以形成解码的差分样本。例如,解码器352执行块浮点解码、Huffman解码或者其他解码。积分运算器354添加解码的差分样本以逆转由压缩所执行的一阶或更高阶差分,以形成提供给诊断信息处理器180的解压缩的波束。如果未执行差分用于压缩,积分运算器354将被绕过。解压缩控制器350将控制参数提供给解码器352和积分运算器354。解压缩控制器350可以从压缩的数据包的报头提取控制数据用以为解压缩操作确定控制参数。
如关于图13所述,图15的解压缩器240可以在解复用之前应用到不同采用的压缩组的复用的序列。在这种情况中,解码器352解码复用的序列中的压缩组Gij以形成在组次序中的解码的差分样本组。对于块浮点解码,如果指数令牌进行差分编码,那么解码器352使用组次序以积分对应于相同波束的指数令牌,以确定相应组的n_exp值。积分运算器354使用组次序和n_exp来确定对应于相同波束的解码差分样本组。积分运算器354积分相应的组解码差分样本以形成241的序列的解压缩的组Rij。组解复用器254随后将解压缩的组Rij分为相应的解压缩的波束253和255。
备选地,如图16中所示,解码的差分样本组可以在积分之前解复用。解码器352接收比如由图12中的D1所代表的复用的序列。解码器352逆转编码器332i的操作以解包并且解码压缩组以形成在组次序中的解码的差分样本组。对于其中指数令牌被差分编码的块浮点解码,解码器352使用组次序以积分对应于相同波束的指数令牌。解码的差分样本组具有相同的每样本比特数量并且在对应于压缩组的复用的序列的组次序中。组解复用器254分离解码的差分样本组以形成对应于每个波束的解码的差分样本的阵列。每个阵列中的解码的差分样本的次序对应于输出自差分运算器330i的差分样本的次序。积分运算器354i和354j执行解码的差分样本的相应阵列的一阶或更高阶积分以形成第i和第j解压缩的波束。
压缩的另一备选方式基于中心频率和采样速率对波束形成样本应用算术操作用以产生具有较低幅度的样本。图17给出了示出在用于压缩具有不同的中心频率的波束形成样本的备选方式之下的原理的示例。以对应于图17中的标记为“带1(band1)”的行的基带信号的示例开始,中心频率接近于DC(0Hz)并且连续的样本之间的相位增加小于10度。第一相量图710示出,因为连续的样本之间的相位变化很小,所以连续的样本的差分的幅度相较于样本本身的幅度相对较小。第一示例序列712对应于带1基带信号的样本。由于连续的样本之间的差分相对于样本幅度较小,所以计算一阶或更高阶差分,或者差分编码产生具有小于原始样本的数据宽度的差分样本。对于基带(带1)示例,使用关于图14描述的差分编码压缩是高效的。在超声信号样本或波束形成样本下变频转换至基带的超声系统中,压缩单元210i可以应用差分编码。
图17还给出了中心频率在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,对由三个采样间隔所分隔的样本的操作产生具有较小数据宽度的修改的样本。将对中的样本相加在一起的操作产生可以进行更高效编码的具有较小数据宽度的修改的样本。
对于图17中的带3的示例,中心频率接近fs/4而连续样本之间的相位增加为大约90度。第三相量图730示出,由180度,或者两个采样间隔所分隔的样本,具有相似的幅度和相反的极性。第三示例序列732还示出,每隔一个样本具有相似的幅度和相反的极性。对于带3,将每隔一个样本相加将会产生具有可以进行比原始样本更高效率编码的具有较小数据宽度的修改的样本。
对于图17中的带4的示例,中心频率接近fs/3而连续样本之间的相位增加为大约120度。第四相量图740示出,由360度,或者三个采样间隔所分隔的样本,将具有相似的幅度。第四示例序列742示出,每三个样本具有相似的幅度。在这种情况中,在由3个采样间隔所分隔的样本之间形成差分将会给出可以进行比原始样本更高效率编码的具有较小数据宽度的修改的样本。
对于图17中的带5的示例,中心频率接近fs/2而连续样本之间的相位增加为大约180度。第五相量图750示出,由180度,或者一个采样间隔所分隔的样本,将具有相似的幅度以及相反的极性。第五示例序列752示出连续样本具有相似的幅度和相反的极性。在这种情况中,将两个连续样本相加将形成可以进行比原始样本更高效率编码的具有较小数据宽度的变更的样本。
上述对于图17描述的示例示出,根据采样速率与中心频率的比率,幅度降低能够通过在由1、2或3个采样间隔所分隔的波束形成样本上,执行比如加(或者反相随后减)或减(或者反相随后加)的操作而实现。产生的修改的样本随后被编码以形成压缩的样本。可以根据中心频率与采样速率的比率,将类似的操作应用到由四个或更多个采样间隔所分隔的样本,以产生具有比原始样本更小的数据宽度的修改的样本。
图18为基于波束形成样本的中心频率的压缩算法的框图。对应于特定波束的波束形成器输出通道161i输入到重新排序解复用器810。重新排序解复用器810选择波束形成样本,从而使选择的样本依据压缩控制参数852由合适数量的采样间隔所分隔,以形成解复用器输出812。算术运算器830依据压缩控制参数856在解复用器输出样本对812上执行加或减运算,以形成修改的样本832。算术运算器830还可以被配置以在解复用器输出样本812上执行更高阶差分。编码器840编码修改的样本832以形成压缩的波束形成样本。编码器840应用块浮点编码、Huffman编码或者其他编码以形成压缩的样本。针对关于图7和11描述的块浮点编码,将修改的样本832提供给块浮点编码器的输入401。
压缩控制器860基于采样速率与波束形成样本的中心频率的比率,将控制参数提供给压缩器元件。重新排序解复用器810和算术运算器830相应地响应于压缩控制参数852和856,以执行合适的操作。图19示出了基于中心频率产生修改的样本的操作。第一列871给出了对于这个示例的可能的中心频率。第二列872给出了每个中心频率的相应的频带指示符。指示符可以作为压缩控制参数852和856使用。第三列873给出了依据压缩控制参数852在重新排序解复用器输出812处产生的样本x(i)和x(j)的不同分隔。第四列874示出了依据压缩控制参数856选择加或减的算术运算的结果。当反相器在延迟的样本“上”时,x(i-j)被减去。第五列875示出了产生修改的样本832或y(i)的算术运算器830的数学结果。压缩控制器860还提供编码器840的控制。压缩控制参数858可以指示块浮点编码、Huffman编码或者其他编码技术的参数。
图20给出了如关于图18和19描述的对于不同中心频率计算出的图17的示例的样本x(i)与x(i-j)的和或者差。样本的示例序列与图17中的相同。示例序列912和942中的DIFF行以及示例序列922、932和952中的SUM行中的样本具有基本上低于相应的样本或x(i)的幅度。DIFF样本和SUM样本为输入至图18中的编码器840的修改的样本932的示例。
如以上关于图6和图12所描述,在图18中压缩器输出211i处由编码器840产生的压缩的波束在数字接口220上传输之前可以与其他压缩的波束复用以形成较少的输出通道。编码器840可以对修改的样本应用块浮点编码从而使图12中的压缩组Gij代表N_GROUP个编码的修改的样本组以形成压缩的波束C1和C2。备选地,编码器840可以对修改的样本应用Huffman编码从而使Gij代表具有每组N_GROUP个Huffman编码的修改的样本的压缩组。组复用器252结合压缩的波束C1和C2的压缩组以形成复用的序列D1。
图21为对于关于图18描述的压缩方法的由解压缩器240所执行的操作的框图。解码器910对压缩的数据解包并且执行解码(例如块浮点解码)操作,以形成解码的修改的样本。逆算术运算器920执行对算术运算器830的逆运算以从解码的修改的样本重建波束形成样本。复用器930将原始样本次序恢复到解压缩的波束形成样本以形成相应的解压缩的波束。解压缩控制器940将控制参数提供给解码器910、逆运算器920以及恢复次序复用器930。解压缩控制器940可以从压缩的数据包的报头中提取控制数据以确定用于解压缩操作的控制参数。
图21的解压缩器240可于解复用组之前应用到修改的样本的压缩组的复用的序列。解码器910逆转编码器840的操作以解码压缩组,按照组次序形成解码的修改的样本组。对于其中指数令牌被差分编码的块浮点解码,解码器910使用组次序来积分指数令牌以确定对应于相同波束的压缩组的n_exp值。逆运算器920以及恢复次序复用器930也使用组次序和n_exp来将它们的相应操作应用到对应于相同波束的解码的修改的样本组。恢复次序复用器930将每组内的解压缩的波束形成样本重新排序到原始样本次序以形成解压缩组;然而,解压缩组仍然按照组次序。参考图13,恢复次序复用器930的输出对应于解压缩的序列241的解压缩组Rij。组解复用器254将解压缩组Rij分为相应的解压缩的波束253和255。
备选地,如图22中所示,解码的修改的样本组可以在逆运算器920之前解复用。解码器910逆转编码器840的操作来解包和解码压缩组以形成按照组次序的相应的解码的差分样本组。对于其中指数令牌被差分编码的块浮点解码,解码器910使用组次序来积分指数令牌以确定对应于相同波束的压缩组的n_exp值。解码的修改的样本组具有相同的每个修改的样本的比特数量。组解复用器254分离解码的修改的样本组以形成对应于每个波束的解码的修改的样本的阵列。逆运算器920i与920j应用到相应的解码的修改的样本的阵列并且恢复次序复用器930i与930j恢复相应的原始样本次序以分别形成第i与第j解压缩的波束。
压缩器210的实施方式对波束形成样本应用简单操作。差分运算器330i(图14)包括一个或多个减法器。块浮点编码(图7和11)使用比较器、减法器以及查询表。备选地,Huffman编码使用查询表为值分配代码。关于图18描述的压缩操作包括解复用、加以及减。用于解压缩器240的实施方式应用简单操作以对压缩的样本进行解压缩。解压缩器240包括用于块浮点解码的查询表和加法器。积分运算器354(图15)包括一个或多个用于积分解码的样本的加法器。图21中的解压缩器240的操作包括加、减和复用。
本发明的优选实施方式在超声系统中将压缩器210和接收波束形成器160集成在具有耦合以接收来自ADC输出通道121i的超声信号样本的输入数据通道的单个专用集成电路(application specificintegrated circuit,简称ASIC)设备中。压缩器210的优选实施包括多个并行的压缩核,其中每个压缩核耦合到波束形成器输出通道161i中的一个并且对相应的波束实施一个压缩单元210i的压缩操作。备选地,一个压缩核可以实施多个压缩单元210i以压缩多个波束。对于这种备选方式,压缩核包括缓冲器用以存储对应于不同波束的波束形成样本,直到它们被处理。波束形成与压缩操作还可以在现场可编程门阵列(field programmable gate array,简称FPGA)中实施。压缩的波束形成样本可以在LVDS端口270i上输出至数字接口220。用于LVDS接口的IP核可商业得到用于ASIC和FPGA实施。备选架构包括在与接收波束形成器160分离的设备中实施压缩器210。压缩器210可以在ASIC、FPGA或者可编程处理器,比如数字信号处理器(digital signal processor,简称DSP)、微处理器、微控制器、多核CPU(比如IBM Cell),或者图形处理单元(graphicsprocessing unit,简称GPU;比如Nvidia GeForce)中实施。
根据超声系统架构,可以将解压缩器240并入诊断信息处理器180的同一设备或者作为与其不同的设备。解压缩操作可以在ASIC或者FPGA中实施。备选地,解压缩操作可以在可由比如DSP、微处理器、微控制器、CPU或者GPU等可编程处理器执行的软件或固件程序中实施。解压缩器240的优选实施为具有用于可由GPU执行的解压缩操作的指令的软件程序。GPU也可以被编程以实施诊断信息处理器180和扫描转换器140的至少一部分操作。备选地,解压缩的波束形成样本可以传输到另一可编程处理器,比如CPU,用于额外的信号处理操作。
图23为GPU中的解压缩器以及其他用以生成超声图像的操作的实施的框图。当前的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和编程方法可以提供实施平台。对于图23中描绘的实施,GPU设备1000可被编程以执行解压缩器240、诊断信息处理器180(B-模式处理以及多普勒处理)、扫描转换器140以及图像处理器150的操作。GPU设备1000可以包括可由并行处理单元访问的动态随机访问存储器(dynamic random access memory,简称DRAM)1002。DRAM 1002可以存储由GPU的其他处理操作所产生的数据和压缩的和/或解压缩的波束形成样本。系统控制器1010提供从接收的数据产生超声图像的任务的协调并响应于用户指令。CPU1012可以实施支持解压缩的操作,比如从压缩包的报头解码压缩控制参数并将它们提供给GPU设备1000用于配置解压缩操作。DRAM1016可以存储接收自数字接口220的压缩的波束形成样本以及其他CPU操作所需的数据。通信控制器1014将接收自数字接口220的压缩包引导至DRAM 1002或DRAM 1016并且管理系统控制器1010和GPU设备1000之间的数据交换。
在优选的系统架构中,系统控制器1010可以在具有用于显示的屏幕160的计算机的主板中实施。GPU设备1000可以在包括DRAM1002,与系统控制器1010通过PCIe(Peripheral ComponentInterconnect Express,外围组件互连高速)背板链路而通信的图形卡中实施。备选地,GPU设备1000可以在安装在主板上的IC中实施。在ADC库120、接收波束形成器160以及压缩器210安装在数据采集卡中的系统架构中,数字接口220可以由PCIe背板链路来实现。
应该认识到,依照摩尔定律(Moore’s Law),更高水平的集成将提供更加紧凑的设备,从而使系统控制器1010和GPU设备可以在单个IC中实施。例如,Intel公司正在开发包括以向量处理单元增强的x86CPU核的多个实例的多核IC架构。此架构称为Larrabee,在标题为“Larrabee:A Many-Core x86CPU Architecture for VisualComputing”的文件(Seiler等人,ACM Transactions on Graphics,Vol.27,No.3,Article 18,August 2008)中描述。Larrabee架构支持需要包括图形处理的并行处理的应用。“Larrabee Native”编程模型包括用于并行编程和向量化的C/C++和API。
图24为CPU/GPU设备中的解压缩器240以及其他用以产生超声图像的操作的实施的框图。CPU/GPU设备1020可以实施解压缩器240、诊断信息处理器130、扫描转换器140和图像处理器150的操作。CPU/GPU控制器1022协调对接收自数字接口220的压缩的波束形成样本的处理操作并响应于用户输入。CPU/GPU设备1020可以由Larrabee平台或者其他具有集成的CPU和GPU功能的可编程设备所实施。
在ADC库120、接收波束形成器160以及压缩器210安置于换能器头中的系统架构中,数字接口220可以是有线或无线通信链路。对于有线通信链路,数字接口可以由PCIe电缆链路或者光纤链路来实施。对于无线通信链路,数字接口可以提供数字调制以及经由射频信道的压缩包的传输以及接收的压缩包的数字解调。无线链路可以遵守无线通信协议,比如WiFi(IEEE 802.11)或者UWB(ultra-wideband,超宽带)格式。
虽然示出和描述了本发明的优选实施方式,但是应该清楚,本发明不仅仅限定在这些实施方式。众多修改,变化,变体,替代和等同将会在不脱离权利要求中所描述的本发明的精神和范围的前提下对本领域技术人员变得明显。
Claims (39)
1.一种在超声成像系统中处理波束的方法,所述超声成像系统包括应用到在采样窗口期间接收的多个超声信号样本序列以形成一个或多个波束的接收波束形成器,其中每个波束由波束形成样本的阵列代表,其中在采样窗口期间由多个超声换能器元件所输出的多个模拟超声信号的模拟到数字转换产生多个超声信号样本序列,所述方法包括:
压缩通过对在所述采样窗口期间接收到的超声信号样本应用所述接收波束形成器而产生的波束以形成压缩的波束,包括压缩代表特定波束的阵列中的波束形成样本以形成相应的压缩波束的压缩的波束形成样本,其中压缩所述阵列中的特定波束形成样本部分取决于所述阵列中的至少一个其他波束形成样本的特征,其中每个波束独立于对应于所述采样窗口的另一波束而进行压缩;以及
跨越数字接口向信号处理器传输压缩的波束。
2.根据权利要求1的方法,其中在压缩步骤之前,将所述波束形成样本的阵列下变频转换到中间频率或者基带。
3.根据权利要求1的方法,进一步包括:
解压缩接收自所述数字接口的相应的压缩的波束的压缩的波束形成样本的至少一部分,以形成相应的解压缩的波束的解压缩的波束形成样本,其中独立于另一压缩的波束将解压缩应用到所述相应的压缩的波束,其中所述信号处理器进一步处理所述解压缩的波束形成样本。
4.根据权利要求1的方法,其中压缩步骤对代表特定波束的波束形成样本阵列应用块浮点编码,以形成所述相应的压缩的波束。
5.根据权利要求4的方法,进一步包括:
通过对所述相应的压缩的波束的相应部分应用块浮点解码来解压缩接收自所述数字接口的所述压缩的波束形成样本的至少一部分,以形成解压缩的波束的解压缩的波束形成样本,其中所述信号处理器进一步处理所述解压缩的波束形成样本。
6.根据权利要求1的方法,其中压缩步骤进一步包括:
定义对应于所述特定波束的波束形成样本的阵列中连续的波束形成样本组,每组具有预设定数量的波束形成样本;
确定在组中具有最大幅度的波束形成样本的指数值;
对所述组的所述指数值进行编码以形成指数令牌;
形成具有减少的比特数量的尾数来代表所述组中的每个波束形成样本,其中比特的减少的数量基于所述指数值,以及
使用所述指数令牌和预设定数量的尾数来代表所述组的压缩的波束形成样本,以形成相应的压缩组,其中多个压缩组形成所述相应的压缩的波束。
7.根据权利要求6的方法,其中形成尾数的步骤进一步包括:
通过截断或舍入所述波束形成样本而从所述组中的每个波束形成样本中去除一定数量的最不重要比特LSB,以形成减少的尾数,其中去除的LSB的数量基于所述指数值,其中代表所述压缩的波束形成样本的步骤使用所述指数令牌和所述预设定数量的减少的尾数来形成所述相应的压缩组。
8.根据权利要求6的方法,其中所述接收波束形成器产生J个波束,其中J至少是2,其中相应波束的波束形成样本的组按照原始组次序,所述方法进一步包括:
复用对应于所述J个波束的压缩组,以形成M个压缩组的复用序列,其中M小于J,包括交错对应于至少两个压缩的波束的集合的压缩组以形成相应的复用的序列,其中所述相应的复用的序列中的压缩组按照组次序布置,其中传输步骤传输所述M个复用的序列。
9.根据权利要求8的方法,进一步包括:
解压缩接收自所述数字接口的所述相应的复用的序列的压缩组,以形成包含解压缩的波束形成样本的解压缩组,其中所述解压缩组按照相应的复用的序列的组次序;以及
解复用相应的复用的序列的解压缩组,以形成对应于所述压缩的波束的集合的解压缩的波束的集合,包括根据所述组次序重新布置所述解压缩组以将原始组次序恢复至所述解压缩的波束,其中将所述解压缩的波束提供用于进一步的信号处理。
10.根据权利要求6的方法,进一步包括:
通过选择用于解压缩的相应的压缩组来解压缩接收自所述数字接口的所述压缩的波束形成样本的至少一部分,以形成提供用于进一步信号处理的相应的解压缩的波束形成样本组。
11.一种在超声成像系统中处理波束的方法,所述超声成像系统包括应用到在采样窗口期间接收到的多个超声信号样本序列以形成一个或多个波束的接收波束形成器,其中每个波束由波束形成样本的阵列代表,其中在所述采样窗口期间由多个超声换能器元件输出的多个模拟超声信号的模拟到数字转换产生多个超声信号样本序列,所述方法包括:
压缩通过对在采样窗口期间接收到的超声信号样本应用接收波束形成器而形成的波束,以形成压缩的波束,所述压缩包括:
计算代表特定波束的阵列中的连续的波束形成样本之间的差分,以形成相应的差分样本序列,其中所述差分为一阶或更高阶差分;以及
编码相应的差分样本序列,以形成相应的压缩的波束的压缩的波束形成样本;以及
跨越数字接口向信号处理器传输所述压缩的波束。
12.根据权利要求11的方法,其中在压缩步骤之前,将波束形成样本的阵列下变频转换至基带或者中间频率。
13.根据权利要求11的方法,其中计算差分的步骤进一步包括:
计算在由预设定数量的采样间隔所分隔的阵列中的波束形成样本对中的波束形成样本之间的差分,以形成相应的差分样本序列;以及
编码所述相应的差分样本序列,以形成所述相应的压缩的波束的压缩的波束形成样本。
14.根据权利要求13的方法,进一步包括:
在计算差分的步骤之前,依据逆转控制参数,选择性地逆转由预设定数量的采样间隔所分隔的每一对波束形成样本中的一个波束形成样本。
15.根据权利要求11的方法,进一步包括解压缩接收自所述数字接口的所述压缩的波速形成样本的至少一部分,所述解压缩进一步包括:
解码所述相应的压缩波束的所述压缩的波束形成样本,以形成解码的差分样本;以及
积分所述解码的差分样本,以形成相应的解压缩的波束的解压缩的波束形成样本,其中所述积分计算一阶或更高阶积分,以逆转计算差分步骤的差分操作。
16.根据权利要求11的方法,其中编码步骤对相应的差分样本序列应用块浮点编码,以形成所述相应的压缩的波束。
17.根据权利要求16的方法,进一步包括解压缩接收自所述数字接口的所述压缩的波束形成样本的至少一部分,所述解压缩进一步包括:
通过对所述相应的压缩的波束的所述压缩的波束形成样本应用块浮点解码来解码所述压缩的波束形成样本,以形成解码的差分样本;以及
积分所述解码的差分样本,以形成相应的解压缩的波束的解压缩的波束形成样本,其中所述积分计算一阶或更高阶积分,以逆转计算差分步骤的差分操作。
18.根据权利要求11的方法,其中编码步骤进一步包括:
定义在所述相应的差分样本序列中连续的差分样本组,每组具有预设定数量的差分样本;
确定所述组中具有最大幅度的差分样本的指数值;
编码所述组的所述指数值,以形成指数令牌;
形成具有减少的数量的比特的尾数来代表所述组中的每个差分样本,其中比特的减少的数量基于所述指数值;以及
使用所述指数令牌和预设定数量的尾数来代表所述组中的所述差分样本,以形成相应的压缩组,其中多个压缩组形成所述相应的压缩的波束。
19.根据权利要求18的方法,其中形成尾数的步骤进一步包括:
通过截断或舍入所述差分样本而从所述组中的每个差分样本中去除一定数量的最不重要比特LSB,以形成减少的尾数,其中去除的LSB的数量基于所述指数值,其中代表所述差分样本的步骤使用所述指数令牌和预设定数量的减少的尾数来形成所述相应的压缩组。
20.根据权利要求18的方法,其中所述接收波束形成器产生N个波束,其中N至少是2,其中相应波束的差分样本组按照原始组次序,所述方法进一步包括:
复用对应于N个压缩的波束的压缩组来形成M个复用的压缩组的序列,其中M小于N,包括交错对应于至少两个压缩的波束的集合的压缩组以形成相应的复用的序列,其中所述相应的复用的序列中的压缩组按照组次序布置,其中传输步骤传输M个复用的序列。
21.根据权利要求20的方法,进一步包括:
解码接收自所述数字接口的所述相应的复用的序列的压缩组,以形成解码的差分样本的解码组,其中所述解码组按照组次序;
根据所述组次序来积分对应于相同压缩的波束的所述解码的差分样本,以形成包含解压缩的波束形成样本的解压缩组,其中所述解压缩组按照所述相应的复用的序列的组次序;以及
解复用所述相应的复用的序列的解压缩组,以形成对应于所述压缩的波束的集合的解压缩的波束的集合,包括根据所述组次序重新布置解压缩组以将原始组次序恢复到所述解压缩的波束,其中将所述解压缩的波束提供用于进一步信号处理。
22.根据权利要求20的方法,进一步包括:
解码接收自所述数字接口的所述相应的复用的序列的压缩组,以形成包含解码的差分样本的解码组,其中所述解码组按照所述相应的复用的序列的组次序;
解复用所述解码组,以形成对应于所述压缩的波束的集合的解码的差分样本的集合,包括按照所述组次序重新布置所述解码组以将原始组次序恢复到所述解码的差分样本序列;以及
积分每个所述解码的差分样本序列,以形成相应的解压缩的波束的解压缩的波束形成样本,其中所述积分计算一阶或更高阶积分以逆转计算差分步骤的差分操作。
23.根据权利要求18的方法,进一步包括:
通过选择用于解压缩的相应的压缩组而解压缩接收自所述数字接口的所述压缩的波束形成样本的至少一部分,以形成提供用于进一步信号处理的解压缩的波束形成样本的相应组。
24.一种用于在超声成像系统中处理波束的装置,所述超声成像系统包括应用到在采样窗口期间接收到的多个超声信号样本序列以形成一个或多个波束的接收波束形成器,其中每个波束由波束形成样本的阵列代表,其中多个模拟到数字转换器采样在采样窗口期间由多个超声换能器元件输出的多个模拟超声信号以产生提供给所述接收波束形成器的多个超声信号样本序列,所述装置包括:
压缩器,耦合至所述接收波束形成器,以接收通过对采样窗口期间接收到的所述超声信号样本应用所述接收波束形成器而形成的一个或多个波束形成样本阵列,所述压缩器包括一个或多个压缩单元,其中相应的压缩单元独立于对应于所述采样窗口的另一波束而压缩相应波束的波束形成样本,以形成相应的压缩的波束的压缩的波束形成样本,其中所述相应的压缩单元配置为部分地基于相应波束的至少一个其他的波束形成样本的特征来压缩特定的波束形成样本,其中所述压缩器将所述压缩的波束提供给数字接口用以传输至信号处理器。
25.根据权利要求24的装置,进一步包括:
耦合至所述接收波束形成器与所述压缩器之间的一个或多个下变频转换器,每个下变频转换器接收来自所述接收波束形成器的相应波束并且将所述相应波束的波束形成样本转换至中间频率或基带,以形成提供给相应的压缩单元的下变频转换波束。
26.根据权利要求24的装置,进一步包括:
解压缩器,耦合以接收来自所述数字接口的压缩的波束并且将解压缩的波束提供给所述信号处理器,其中所述解压缩器独立于对应于采样窗口的另一压缩的波束而解压缩特定的压缩的波束的压缩的波束形成样本,以形成相应的解压缩的波束的解压缩的波束形成样本。
27.根据权利要求24的装置,其中所述压缩单元进一步包括:
块浮点编码器,其应用到波束形成样本的阵列以产生相应的压缩的波束。
28.根据权利要求27的装置,进一步包括:
解压缩器接收来自所述数字接口的压缩的波束并且将解压缩的波束提供给所述信号处理器,所述解压缩器包括应用到相应的压缩的波束的压缩的波束形成样本的至少一部分以形成相应的解压缩的波束的解压缩的波束形成样本的块浮点解码器。
29.根据权利要求24的装置,其中所述压缩单元进一步包括:
差分运算器,用以计算相应的波束的连续的波束形成样本之间的一阶或更高阶差分,以产生差分样本序列;以及
编码器,用以编码所述差分样本序列,以形成相应的压缩的波束的压缩的波束形成样本。
30.根据权利要求29的装置,其中所述差分运算器计算相应的波束的波束形成样本对中的波束形成样本之间的差分,其中每一对的波束形成样本由预设定数量的采样间隔所分隔以形成所述差分样本序列。
31.根据权利要求29的装置,进一步包括解压缩器,其接收来自所述数字接口的压缩的波束并且将解压缩的波束提供给所述信号处理器,所述解压缩器包括:
解码器,其应用到相应的压缩的波束的压缩的波束形成样本的至少一部分以产生解码的差分样本序列;以及
积分器,用于计算解码的差分样本序列的一阶或更高阶积分以形成相应的解压缩的波束的解压缩的波束形成样本。
32.根据权利要求30的装置,其中所述压缩单元进一步包括:
逆转器,依据逆转控制参数,选择性地逆转每一对由预设定数量的采样间隔所分隔的波束形成样本中的一个波束形成样本,以产生用于所述差分运算器的波束形成样本对。
33.根据权利要求24的装置,其中所述接收波束形成器将J个波束提供给J个压缩单元以形成J个压缩的波束,其中相应的压缩单元产生用于相应的压缩的波束的压缩的波束形成样本的多个压缩组,其中压缩组按照原始组次序,每个压缩组具有预设定数量的压缩的波束形成样本,所述装置进一步包括:
组复用器,其接收J个压缩的波束并且将M个复用的压缩组的序列提供给所述数字接口,其中M小于J,所述组复用器交错对应于至少两个压缩的波束的集合的压缩的样本组以形成相应的复用的序列,其中相应的复用的序列中的压缩组按照组次序布置。
34.根据权利要求33的装置,进一步包括:
解压缩器,其接收来自所述数字接口的M个复用的序列并且形成M个解压缩组序列,其中所述解压缩器解压缩相应的复用的序列的压缩组以按照所述组次序形成解压缩组;以及
组解复用器,其耦合以接收所述M个解压缩组序列并且将N个解压缩的波束提供给所述信号处理器,其中相应的复用的序列的解压缩组从组次序重新布置为原始次序,以形成对应于压缩的波束的集合的解压缩的波束的集合。
35.根据权利要求24的装置,其中压缩器实施于专用集成电路ASIC或者现场可编程门阵列FPGA之中。
36.根据权利要求26的装置,其中所述解压缩器实施于FPGA中。
37.根据权利要求26的装置,其中所述信号处理器实施于FPGA中,所述FPGA进一步包括:
实施所述解压缩器的解压缩核;
数据通道,用于接收来自所述数字接口的压缩的波束以将压缩的波束形成样本提供给所述解压缩核;以及
数据通道,用于将解压缩的波束提供给所述信号处理器。
38.根据权利要求26的装置,其中所述解压缩器至少部分地实施于图形处理单元GPU中。
39.根据权利要求38的装置,其中所述信号处理器至少部分地实施于GPU中。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/494,184 | 2009-06-29 | ||
US12/494,184 US8317706B2 (en) | 2009-06-29 | 2009-06-29 | Post-beamforming compression in ultrasound systems |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101933814A CN101933814A (zh) | 2011-01-05 |
CN101933814B true CN101933814B (zh) | 2012-07-25 |
Family
ID=43381492
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101854538A Active CN101933814B (zh) | 2009-06-29 | 2010-05-20 | 超声系统中的波束形成后压缩 |
CN2010202049382U Expired - Fee Related CN201855272U (zh) | 2009-06-29 | 2010-05-20 | 超声成像系统中的装置 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010202049382U Expired - Fee Related CN201855272U (zh) | 2009-06-29 | 2010-05-20 | 超声成像系统中的装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US8317706B2 (zh) |
EP (1) | EP2448496A2 (zh) |
JP (1) | JP5866282B2 (zh) |
KR (1) | KR101629423B1 (zh) |
CN (2) | CN101933814B (zh) |
WO (1) | WO2011008408A2 (zh) |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013506917A (ja) * | 2009-09-30 | 2013-02-28 | サンプリファイ システムズ インコーポレイテッド | 圧縮及び復元を用いたマルチ・プロセッサの波形データ交換の改善 |
US8301803B2 (en) * | 2009-10-23 | 2012-10-30 | Samplify Systems, Inc. | Block floating point compression of signal data |
KR101335724B1 (ko) * | 2010-08-09 | 2013-12-04 | 삼성전자주식회사 | 초음파 진단기 및 그 제어방법 |
EP2661329B1 (en) * | 2011-01-06 | 2014-09-10 | The Lubrizol Corporation | Improvements in or relating to ultrasound generating apparatus, and methods for generating ultrasound |
TW201232476A (en) * | 2011-01-27 | 2012-08-01 | Univ Nat Taiwan | Detection system and signal processing method thereof |
US8451158B2 (en) * | 2011-06-30 | 2013-05-28 | Lsi Corporation | Analog to digital converter with generalized beamformer |
US9030908B2 (en) * | 2011-09-23 | 2015-05-12 | Texas Instruments Incorporated | Programmable wavelet tree |
JP5769678B2 (ja) * | 2011-11-10 | 2015-08-26 | 富士フイルム株式会社 | 超音波診断装置および超音波画像生成方法 |
CN103284753B (zh) * | 2012-02-22 | 2015-12-09 | 香港理工大学 | 超声波成像系统及成像方法 |
KR101242368B1 (ko) * | 2012-04-23 | 2013-03-15 | 주식회사 힐세리온 | 2차원 배열 데이터를 이용한 모바일 초음파 진단 프로브 장치 및 이를 이용한 모바일 초음파 진단 시스템 |
US11126394B2 (en) | 2012-05-01 | 2021-09-21 | Lisnr, Inc. | Systems and methods for content delivery and management |
US11452153B2 (en) | 2012-05-01 | 2022-09-20 | Lisnr, Inc. | Pairing and gateway connection using sonic tones |
KR101319033B1 (ko) * | 2012-06-25 | 2013-10-15 | 주식회사 힐세리온 | 2차원 배열 데이터를 이용한 모바일 초음파 진단 시스템, 이를 위한 모바일 초음파 진단 프로브 장치, 및 초음파 진단 장치 |
WO2014125371A1 (en) * | 2013-02-12 | 2014-08-21 | Urs-Us Medical Technology Inc. | Analog store digital read ultrasound beamforming system and method |
US9042462B2 (en) | 2013-04-24 | 2015-05-26 | Commscope Technologies Llc | Differential signal transmission |
KR101978178B1 (ko) * | 2013-05-24 | 2019-05-15 | 삼성전자주식회사 | 초음파 데이터를 처리하는 데이터 처리 장치 및 방법 |
KR101496168B1 (ko) * | 2014-04-28 | 2015-02-26 | 주식회사 힐세리온 | 다양한 해상도를 갖는 무선통신단말기에 적용이 가능한 초음파 진단시스템 및 진단방법 |
US9059970B1 (en) * | 2014-04-28 | 2015-06-16 | Innowireless Co., Ltd. | Method of compressing and decompressing I/Q sample data in cloud radio access network |
EP3207639A4 (en) * | 2014-10-15 | 2018-04-11 | Lisnr, Inc. | Inaudible signaling tone |
WO2016159413A1 (ko) * | 2015-04-02 | 2016-10-06 | 알피니언메디칼시스템 주식회사 | 보간으로 인한 흐름 현상을 완화한 스캔 변환 방법 및 장치 |
US11684346B2 (en) | 2015-05-29 | 2023-06-27 | Siemens Medical Solutions Usa, Inc. | Ultrasound beamformer-based channel data compression |
KR102560033B1 (ko) * | 2015-11-09 | 2023-07-26 | 삼성전자주식회사 | 프로브 장치 및 그 제어 방법 |
US10677901B2 (en) | 2016-04-08 | 2020-06-09 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the ultrasound diagnosis apparatus |
EP4145176A1 (en) | 2016-04-29 | 2023-03-08 | R2Sonic, LLC | Sonar data compression |
US10132924B2 (en) * | 2016-04-29 | 2018-11-20 | R2Sonic, Llc | Multimission and multispectral sonar |
WO2017207815A1 (en) * | 2016-06-02 | 2017-12-07 | Koninklijke Philips N.V. | Ultrasound systems with time compression and time multiplexing of received ultrasound signals |
US10154253B2 (en) * | 2016-08-29 | 2018-12-11 | Disney Enterprises, Inc. | Multi-view displays using images encoded with orbital angular momentum (OAM) on a pixel or image basis |
CA3034741A1 (en) * | 2016-09-28 | 2018-04-05 | Covidien Lp | System and method for parallelization of cpu and gpu processing for ultrasound imaging devices |
GB2560167B (en) * | 2017-02-28 | 2021-12-22 | Imv Imaging Uk Ltd | Ultrasound Imaging apparatus and methods |
WO2019010022A1 (en) * | 2017-07-03 | 2019-01-10 | R2Sonic, Llc | SYSTEM AND METHOD FOR MULTI-PERSPECTIVE INSONIFICATION |
US10067228B1 (en) * | 2017-09-11 | 2018-09-04 | R2Sonic, Llc | Hyperspectral sonar |
US11189295B2 (en) | 2017-09-28 | 2021-11-30 | Lisnr, Inc. | High bandwidth sonic tone generation |
CN111479514A (zh) | 2017-11-15 | 2020-07-31 | 蝴蝶网络有限公司 | 超声设备和用于制造超声装置的方法 |
US10816652B2 (en) * | 2018-02-28 | 2020-10-27 | Codaoctopus Group | Method of compressing sonar data |
US10211906B1 (en) * | 2018-03-22 | 2019-02-19 | Intel Corporation | Two-stage spatial compression method |
US10718865B2 (en) * | 2018-05-14 | 2020-07-21 | Coda Octopus Group | Method of compressing beamformed sonar data |
US20190361102A1 (en) * | 2018-05-22 | 2019-11-28 | Analog Devices, Inc. | Delay and apodization control interface for ultrasound beamformer |
TWI676043B (zh) * | 2018-11-08 | 2019-11-01 | 立積電子股份有限公司 | 超寬頻雷達收發器及其運作方法 |
CN111914045B (zh) * | 2020-07-09 | 2023-06-30 | 珠海云洲智能科技股份有限公司 | 一种数据压缩方法、装置、终端设备及存储介质 |
US11882201B2 (en) * | 2022-03-30 | 2024-01-23 | Itron, Inc. | Data compression techniques for efficient network management |
US20240012140A1 (en) * | 2022-07-08 | 2024-01-11 | Cirrus Logic International Semiconductor Ltd. | Systems and methods for generating a signal |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61170442A (ja) | 1985-01-23 | 1986-08-01 | 松下電器産業株式会社 | 超音波ドツプラ血流装置 |
JPH07181996A (ja) * | 1993-12-22 | 1995-07-21 | Sony Corp | 情報処理方法、情報処理装置、及びメディア |
KR0180058B1 (ko) | 1995-09-13 | 1999-04-01 | 이민화 | Cine메모리에 데이타를 압축 저장하는 초음파진단장치 |
JP3300287B2 (ja) | 1998-06-10 | 2002-07-08 | 松下電器産業株式会社 | 超音波診断装置 |
WO2000010638A2 (en) | 1998-08-24 | 2000-03-02 | Baskent University | An asynchronous oversampling beamformer |
US6042545A (en) | 1998-11-25 | 2000-03-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for transform ultrasound processing |
AU1726900A (en) | 1998-11-24 | 2000-06-13 | Acuson Corporation | Processing ultrasound data during or after compression |
JP4210386B2 (ja) | 1999-03-30 | 2009-01-14 | テルモ株式会社 | 体腔内超音波診断装置 |
US7115093B2 (en) | 2001-11-21 | 2006-10-03 | Ge Medical Systems Global Technology Company, Llc | Method and system for PDA-based ultrasound system |
JP3910860B2 (ja) | 2002-02-05 | 2007-04-25 | 株式会社日立メディコ | 超音波撮像装置 |
JP2003265468A (ja) | 2002-03-19 | 2003-09-24 | Ge Medical Systems Global Technology Co Llc | 診断情報生成装置および超音波診断装置 |
JP4049791B2 (ja) * | 2003-04-28 | 2008-02-20 | 日本電信電話株式会社 | 浮動小数点形式ディジタル信号可逆符号化方法、及び復号化方法と、その各装置、その各プログラム |
JP4469583B2 (ja) | 2003-09-11 | 2010-05-26 | 株式会社東芝 | 超音波診断装置 |
US7998072B2 (en) | 2003-12-19 | 2011-08-16 | Siemens Medical Solutions Usa, Inc. | Probe based digitizing or compression system and method for medical ultrasound |
US7009533B1 (en) | 2004-02-13 | 2006-03-07 | Samplify Systems Llc | Adaptive compression and decompression of bandlimited signals |
US20060074320A1 (en) | 2004-08-27 | 2006-04-06 | Yoo Yang M | Home ultrasound system |
US20070239019A1 (en) | 2006-02-13 | 2007-10-11 | Richard William D | Portable ultrasonic imaging probe than connects directly to a host computer |
JP4713382B2 (ja) * | 2006-03-28 | 2011-06-29 | 富士フイルム株式会社 | 超音波診断装置及びデータ解析計測装置 |
JP4908928B2 (ja) * | 2006-05-30 | 2012-04-04 | 日立アロカメディカル株式会社 | ワイヤレス超音波診断装置 |
US20080114246A1 (en) | 2006-11-10 | 2008-05-15 | Penrith Corporation | Transducer array imaging system |
US7844097B2 (en) * | 2007-12-03 | 2010-11-30 | Samplify Systems, Inc. | Compression and decompression of computed tomography data |
-
2009
- 2009-06-29 US US12/494,184 patent/US8317706B2/en not_active Expired - Fee Related
-
2010
- 2010-05-20 CN CN2010101854538A patent/CN101933814B/zh active Active
- 2010-05-20 CN CN2010202049382U patent/CN201855272U/zh not_active Expired - Fee Related
- 2010-06-17 EP EP10800225A patent/EP2448496A2/en not_active Withdrawn
- 2010-06-17 JP JP2012518545A patent/JP5866282B2/ja not_active Expired - Fee Related
- 2010-06-17 KR KR1020127002321A patent/KR101629423B1/ko active IP Right Grant
- 2010-06-17 WO PCT/US2010/039081 patent/WO2011008408A2/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
KR20120108959A (ko) | 2012-10-05 |
WO2011008408A3 (en) | 2011-04-07 |
JP5866282B2 (ja) | 2016-02-17 |
JP2012531966A (ja) | 2012-12-13 |
US20100331689A1 (en) | 2010-12-30 |
WO2011008408A2 (en) | 2011-01-20 |
US8317706B2 (en) | 2012-11-27 |
CN101933814A (zh) | 2011-01-05 |
EP2448496A2 (en) | 2012-05-09 |
KR101629423B1 (ko) | 2016-06-10 |
CN201855272U (zh) | 2011-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101933814B (zh) | 超声系统中的波束形成后压缩 | |
CN101904751B (zh) | 超声信号压缩 | |
AU2018222929B2 (en) | Ultrasound signal processing circuitry and related apparatus and methods | |
CN106456115B (zh) | 超声成像压缩方法及设备 | |
CN104254773B (zh) | 使用二维阵列数据的移动超声诊断探测设备以及使用该设备的移动超声诊断系统 | |
US6760486B1 (en) | Flash artifact suppression in two-dimensional ultrasound imaging | |
CN106210719B (zh) | 用于超声成像的系统 | |
US20060074320A1 (en) | Home ultrasound system | |
CN101569539A (zh) | 远程图像传输方法、远程超声波诊断系统及装置 | |
JP2000149015A (ja) | 画像のエッジ強調方法及びイメ―ジング装置 | |
US20230240664A1 (en) | Method and system for data transfer reduction in ultrasound imaging | |
Ibrahim et al. | Single-FPGA complete 3D and 2D medical ultrasound imager | |
Levesque et al. | Novel low-power ultrasound digital preprocessing architecture for wireless display | |
Dallai et al. | Real-time bilateral filtering of ultrasound images through highly optimized DSP implementation | |
CN105559819A (zh) | 基于精简硬件的超声成像系统 | |
EP2510381B1 (en) | Method of re-sampling ultrasound data | |
Tae-Sung et al. | A Comparative Study of Compression Methods and the Development of CODEC Program of Biological Signal for Emergency Telemedicine Service | |
JP2000316851A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
ASS | Succession or assignment of patent right |
Owner name: ALTERA CORP. Free format text: FORMER OWNER: SAMPLIFY SYSTEMS INC. Effective date: 20140707 |
|
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20140707 Address after: American California Patentee after: Altera Corp. Address before: American California Patentee before: Samplify Systems Inc. |