CN104688271B - 合成聚焦超声成像方法和装置 - Google Patents

合成聚焦超声成像方法和装置 Download PDF

Info

Publication number
CN104688271B
CN104688271B CN201510142749.4A CN201510142749A CN104688271B CN 104688271 B CN104688271 B CN 104688271B CN 201510142749 A CN201510142749 A CN 201510142749A CN 104688271 B CN104688271 B CN 104688271B
Authority
CN
China
Prior art keywords
channel data
sparse
data collection
signal
synthetic focusing
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
Application number
CN201510142749.4A
Other languages
English (en)
Other versions
CN104688271A (zh
Inventor
罗建文
刘敬
何琼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201510142749.4A priority Critical patent/CN104688271B/zh
Publication of CN104688271A publication Critical patent/CN104688271A/zh
Application granted granted Critical
Publication of CN104688271B publication Critical patent/CN104688271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • A61B8/145Echo-tomography characterised by scanning multiple planes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8997Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography

Abstract

本发明提供一种合成聚焦超声成像方法和装置。该方法包括:激励超声探头的多个阵元多次发射平面波,其中所述多个阵元每次发射平面波时的发射变迹与其元素成随机分布的测量矩阵中的相应行相对应;每发射一次平面波后,激励所述超声探头的所有阵元接收回波信号,以获得通道数据;根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢复合成聚焦通道数据集;以及对所述合成聚焦通道数据集进行波束合成,以生成超声波图像。利用上述方法和装置所获得的超声波图像不仅具有很高分辨率;而且因为超声波的发射次数较少,每次发射能量高,因此图像还具有高帧频、高对比度的优点。

Description

合成聚焦超声成像方法和装置
技术领域
[0001] 本发明涉及超声波图像成像领域,具体地,涉及一种合成聚焦超声成像方法和装 置。
背景技术
[0002] 超声成像因为具有实时,廉价,非侵入性和非电离辐射等优点而广泛地用于临床 诊断。空间分辨率,时间分辨率和对比度是三个评价超声图像的指标,高时空分辨率和高对 比度的超声图像能更好的辅助临床诊断。但可惜的是这三者并不能同时达到最优。
[0003] 传统的超声成像中,通过发射聚焦波而得到一条图像扫描线。这种模式所得到的 图像具有较高的空间分辨率和对比度,特别是在发射聚焦区域。但是因为一幅图像通常需 要上百条扫描线从而需要上百次的发射次数,所以这种模式的帧频,即时间分辨率,不够 高。虽然减小每次发射的聚焦区域的个数可以有效地提高帧频,但这样会降低图像的空间 分辨率和对比度。
[0004] 合成聚焦法作为一种合成孔径方法,早在1992年就被提出用于提高超声图像的空 间分辨率。在合成聚焦法中,对于具有128个阵元的超声探头,依次激励每个阵元发射超声 波,每个阵元发射超声波后,所有阵元一起接收回波信号。图1示出了现有的合成聚焦超声 成像以及其中原始信号 x的提取的示意图。斜直线表示超声探头。斜直线上的黑色矩形框表 示合成聚焦模式下每次发射时的激活阵元,而探头上的其他阵元未被激活。向上实线箭头 表示回波信号的接收过程,图1示出共进行了 η次发射和接收,向右虚线箭头表示原始信号X 的提取过程。对于128个阵元的超声探头,η取128。每一次发射和接收形成一幅低质量的超 声图像。根据η幅低质量的图像即可得到一幅高质量的图像。具体地,用矩阵Xi表示所有阵 元第i次接收到的合成聚焦通道数据。该矩阵Xi有2*d*fs/c行,128列,其中,d为采样深度, fs为采样频率,c为声速。矩阵Xi中的每列对应相应阵元接收到的回波信号。经η次发射、接 收后,数据集X: {X1;X2,.. .XJ即为合成聚焦通道数据集,其为合成聚焦通道数据&的集合。
[0005] 因为合成聚焦法能够实现发射和接收的双向动态聚焦,所以该方法获得的超声图 像具有很高的空间分辨率。其分辨率等价于传统模式下有无限个发射聚焦区域的情况。但 因为合成聚焦法中每次只有一个阵元被激活,能量非常低,所以获得的超声图像的对比度 不高。基于这种情况,合成发射孔径(synthetic transmit aperture)法被提出用于提高超 声图像的对比度。基于合成发射孔径法的超声成像中,将阵列划分为多个子孔径,每个子孔 径由多个连续的阵元构成。通过依次激活每一个子孔径发射超声波,所有阵元接收回波信 号来得到超声图像。这种方法和合成聚焦法相比,因为每次激活的阵元数多而提高了图像 对比度,而且其发射次数也比合成聚焦少,所以能得到更高的帧率。但是这种模式却牺牲了 空间分辨率。
发明内容
[0006] 为了至少部分地解决上述技术问题,根据本发明的一个方面,提供一种合成聚焦 超声成像方法,包括:
[0007] 激励超声探头的多个阵元多次发射平面波,其中所述多个阵元每次发射平面波时 的发射变迹(transmit apodization)与其元素成随机分布的测量矩阵中的相应行相对应;
[0008] 每发射一次平面波后,激励所述超声探头的所有阵元接收回波信号,以获得通道 数据;
[0009] 根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢复合成聚焦通道数据 集;以及
[0010] 对所述合成聚焦通道数据集进行波束合成,以生成超声波图像。
[0011] 可选地,所述随机分布是0至1的随机分布。
[0012] 可选地,所述恢复合成聚焦通道数据集的步骤进一步包括:
[0013] 根据所述通道数据集、所述测量矩阵和稀疏基,利用所述压缩感知重建算法计算 待恢复信号的稀疏表达;以及
[0014] 利用所述稀疏基对所述待恢复信号的稀疏表达进行稀疏反变换,从而恢复出所述 合成聚焦通道数据集。
[0015] 可选地,所述稀疏基是小波基。
[0016] 可选地,所述计算待恢复信号的稀疏表达是利用1-范数最小求解方法。
[0017] 可选地,所述计算待恢复信号的稀疏表达之前,所述通道数据被归一化。
[0018] 可选地,对于通道数据集中每个矩阵的第m行、第η列的元素全部为0的情况,将所 述合成聚焦通道数据集的每个矩阵的第m行、第η列的元素全部直接置为0。
[0019] 根据本发明另一方面,还提供了 一种合成聚焦超声成像装置,包括:
[0020] 激励发射模块,用于激励超声探头的多个阵元多次发射平面波,其中所述多个阵 兀每次发射平面波时的发射变迹与其兀素成随机分布的测量矩阵中的相应行相对应;
[0021] 激励采集模块,每发射一次平面波后,激励所述超声探头的所有阵元接收回波信 号,以获得通道数据;
[0022] 数据恢复模块,用于根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢 复合成聚焦通道数据集;以及
[0023] 图像重建模块,用于对所述合成聚焦通道数据集进行波束合成,以生成超声波图 像。
[0024] 可选地,所述随机分布是0至1的随机分布。
[0025] 可选地,所述数据恢复模块进一步包括:
[0026] 第一计算模块,用于根据所述通道数据集、所述测量矩阵和稀疏基,利用所述压缩 感知重建算法计算待恢复信号的稀疏表达;以及
[0027] 第二计算模块,用于利用所述稀疏基对所述待恢复信号的稀疏表达进行稀疏反变 换,从而恢复出所述合成聚焦通道数据集。
[0028]利用上述方法和装置所获得的超声波图像不仅具有很高分辨率;而且因为超声波 的发射次数较少,每次发射能量高,因此图像还具有高帧频、高对比度的优点。
[0029] 在发明内容中引入了一系列简化形式的概念,这将在具体实施方式部分中进一步 详细说明。本发明内容部分并不意味着要试图限定出所要求保护的技术方案的关键特征和 必要技术特征,更不意味着试图确定所要求保护的技术方案的保护范围。
[0030] 以下结合附图,详细说明本发明的优点和特征。
附图说明
[0031] 本发明的下列附图在此作为本发明的一部分用于理解本发明。附图中示出了本发 明的实施方式及其描述,用来解释本发明的原理。在附图中,
[0032] 图1示出了现有的合成聚焦超声成像以及其中原始信号的提取的示意图;
[0033]图2示出了根据本发明一个具体实施例的合成聚焦超声成像以及其中测量信号的 提取的不意图;
[0034] 图3是根据本发明一个具体实施例的合成聚焦超声成像方法的流程图;
[0035] 图4示出了根据本发明一个具体实施例的测量信号的示意图;
[0036] 图5示出了根据本发明一个具体实施例的所恢复信号的示意图;
[0037] 图6示出了根据本发明一个具体实施例的待恢复信号的稀疏表达的示意图;以及
[0038] 图7示出了根据本发明一个具体实施例的合成聚焦超声成像装置的示意性框图。
具体实施方式
[0039] 在下文的描述中,提供了大量的细节以便能够彻底地理解本发明。然而,本领域技 术人员可以了解,如下描述仅涉及本发明的较佳实施例,本发明可以无需一个或多个这样 的细节而得以实施。此外,为了避免与本发明发生混淆,对于本领域公知的一些技术特征未 进行描述。
[0040] 压缩感知理论表明,如果信号xeIR"是稀疏的或者在某个变换域是可压缩的,那么 其可以从远低于奈奎斯特采样频率的信号中被高概率的恢复。因为这个性质,压缩感知可 以应用于数据压缩,通道编码,逆问题求解和数据采集。目前压缩感知理论在超声领域的应 用主要集中在减少数据量的问题上。Liebgott和Eldar小组在这方面做了很多的研究,通过 将压缩采样应用在快速时间域(fast time domain),通道数据量能大大减少。其中,快速时 间域是沿着超声回波信号采集的方向。数据量的减少对便携式超声系统非常有帮助,但并 不能提高帧频。帧频可以按c/2dn计算,其中c是声速,d是成像深度,η是发射次数。一旦成像 目标和成像深度确定,而在快速时间域进行的压缩采样也不能改变发射次数η,即帧频不会 发生变化。所以期望提高帧频,可以通过降低发射次数的方式,即将压缩采样应用在慢速时 间域(slow time domain)。慢速时间域是沿着发射-接收重复方向。
[0041] 本发明将压缩感知这种新理论应用到超声成像领域,以同时得到高帧频,高分辨 率和高对比度的超声图像。图2示出了根据本发明一个具体实施例的合成聚焦超声成像以 及其中测量信号y的提取的示意图。本发明提供的合成聚焦超声成像方法将压缩采样应用 到了图2所示的慢速时间域。图3是根据本发明一个具体实施例的合成聚焦超声成像方法 300的流程图。
[0042] 如图3所示,合成聚焦超声成像方法300包括以下步骤。
[0043] S310,激励超声探头的多个阵元发射平面波。其中该多个阵元发射平面波时的发 射变迹与其元素成随机分布的测量矩阵Φ中的相应行相对应。
[0044]测量矩阵Φ的列数等于探头阵元个数。在本发明的示例中,以一个具有128个阵元 的探头为例,假设其全部阵元都参与发射平面波,则测量矩阵Φ的列数为128。测量矩阵Φ 中的元素服从随机分布。优选地,测量矩阵Φ中的元素服从0至1的随机分布。通过0至1的随 机分布可以显著减少最终获得的超声波图像的伪影,使图像更准确。更优选地,测量矩阵Φ 中的元素服从0至1的均匀随机分布。0至1的均匀随机分布可以保证成像过程中超声波的能 量较高,最终获得的超声波图像对比度高,便于医生等用户观察。可选地,测量矩阵Φ可以 是随机高斯矩阵。
[0045] 如图2所示,矩形框表示阵元,矩形框的不同灰度表示发射变迹大小。对比上面描 述的合成聚焦模式,其需要进行例如128次发射,每次发射只有一个阵元被完全激活。根据 线性声场理论,一次平面波发射得到的通道数据是128个单阵元中的每个发射一次超声波 得到的通道数据集的线性组合,其中的线性系数即为探头发射平面波时的发射变迹。每次 发射平面波时的发射变迹对应于测量矩阵Φ的相应行。即,测量矩阵Φ中的第1行是第一次 发射平面波时的发射变迹,第2行是第二次发射平面波时的发射变迹,依次类推。因此,根据 本发明的实施例发射一次有发射变迹的平面波就得到了所有128次单阵元发射的信息,也 即在慢时间域完成了压缩采样。
[0046] S320,每发射一次平面波后,激励超声探头的所有阵元接收回波信号,以获得通道 数据Yi。
[0047] 重复上述步骤S310和步骤S320,直至完成总发射次数,例如64。由此,获得通道数 据的集合:通道数据集Y: {Y^Ys,. . .Y64},也就是结合压缩感知重建算法的合成聚焦通道数 据集。其中Yi为矩阵。通道数据Yi和合成聚焦通道数据Xi的维度一样,行数等于2*d*fs/c。其 中,d为米样深度,fs为米样频率,c为声速。列数和探头阵兀数目相等,在此实施例中为128。 每列表示平面波发射后,每个阵元接收到的回波信号。
[0048] 根据步骤S310的描述,本领域普通技术人员可以理解,测量矩阵Φ的行数等于总 发射次数。
[0049] 本领域普通技术人员可以理解,总发射次数可以是小于阵元个数的任意数。优选 地,总发射次数取32至112次。成像目标为简单的点目标或重建算法是基追踪算法(Basis P UrsUit,BP)时,总发射次数优选为32至64次。成像目标为复杂的生物组织或重建算法是正 交匹配追踪算法(Orthogonal Matching Pursuit,OMP)时,总发射次数优选为64至112次。
[0050] 本领域普通技术人员可以理解,如果在上述步骤S310中,不是所有阵元都参加平 面波的发射,那么其余未参加平面波发射的阵元仍可以依次各自发射超声波,并且在其中 任何一个阵元发射了超声波之后,超声探头的所有阵元都接收回波信号,以获得对应的通 道数据供成像计算,如现有的合成聚焦超声成像方法中一样。
[0051] S330,根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢复合成聚焦通 道数据集。
[0052] 合成聚焦通道数据集X: {X^Xs,. . .X128}是利用现有的合成聚焦超声成像方法经 128次超声波发射和接收得到的通道数据集。利用测量矩阵Φ,应用压缩感知重建算法即可 从本实施例的通道数据集Y: {Yi,Y2, . . .Υ64}中恢复出合成聚焦通道数据集。具体地,遍历所 有接收通道,所有采样深度,对通道数据集Y沿慢时间域提取测量信号y。根据测量信号y和 测量矩阵Φ,利用压缩感知重建算法即可恢复待恢复信号X。将这个恢复过程遍历所有接收 通道,所有采样深度即可得到完整的合成聚焦通道数据集。
[0053] 采用上面提到的数学符号,步骤S330即从通道数据集Y: m,Y2, . . .Y64}恢复合成 聚焦通道数据集X: {&,心,...Xn}的过程,其中η等于探头阵元数,例如是128。
[0054] 对Y^Ys,.. .Υ64每个矩阵的第一行,第一列,沿慢时间域提取测量信号,得第一行、 第一列的测量信号y= [yi,y2,.. .y64](即第一个采样深度,第一个接收通道的测量信号), yi,y2, . . .y64的提取如图2所示。利用压缩感知重建算法从测量信号y= [yi,y2,. . .y64]恢复 出其对应的待恢复信号x= [X1,X2, . . .,χ128],并将Xl写入X中每个矩阵第一行和第一 列;对所有行和所有列重复上述过程,则合成聚焦通道数据集X被恢复出来。如此,则利用压 缩感知实现了从64次平面波发射得到的通道数据集Y: {Y^Ys,.. .Y64}恢复出128次传统合 成聚焦发射得到的合成聚焦通道数据集X: {Xi,X2,. . .Χη}。图4和图5分别示出了根据本发明 一个具体实施例的测量信号和所恢复信号的示意图。图4和图5中的纵坐标均为信号的幅 度。图4的横坐标为平面波发射索引。图5的横坐标表示合成聚焦发射通道。
[0055] 压缩感知理论表示,若信号X是稀疏的或在某一个变换域是可压缩的,即χ=Ψν (公式1)是稀疏基,ν只有少数几个是非零的或其元素的顺序排列服从一个指数 衰减,那么该信号X可以从其少数几个线性测量信号y= Φχ (公式2)中恢复出来,Φ是测量 矩阵。在本实施例中,根据线性声场理论,测量矩阵Φ的每行对应发射变迹。上述稀疏基Ψ 可以为小波基。例如sys8小波基。小波基对在慢时间域提取的超声信号有较好的稀疏表达, 信号恢复效果更好,而且成像的帧率更高。图6示出了根据本发明一个具体实施例的待恢复 信号的稀疏表达v的示意图,其中纵坐标表示信号的幅度,横坐标表示排序后变换域的索 引。
[0056] S340,对合成聚焦通道数据集进行波束合成,以生成超声波图像。
[0057] 该步骤S340采用现有的波束合成方法即可,本领域普通技术人员可以理解该步骤 的实现过程,为了简洁,在此不再赘述。
[0058] 本发明的实施例中,利用一次平面波发射是完整合成聚焦的一次压缩采样的原 理,相比现有的合成聚焦超声成像方法,可以发射较少次数的超声波,成功地在慢时间域实 现了压缩采样,提高了图像帧频。同时因为一次平面波发射所有阵元都被激活,能量较高, 所以也提高了图像对比度。因为上述方法得到是合成聚焦图像,所以保留了其高分辨率的 特点。总之,利用压缩感知的技术,得到了高帧频,高分辨率和高对比度的超声图像。
[0059] 此外,上述方法对超声成像领域具有普适性,特别是一些对成像帧频和成像质量 同时有较高要求的应用,如实时三维超声成像,心血管成像等。除了传统的结构成像,该方 法还能帮助提供高质量的超声功能成像,如心血管弹性成像等。
[0060] 根据本发明一个具体实施例,上述步骤S330恢复合成聚焦通道数据集的操作可以 进一步包括以下子步骤。
[0061] S331,根据通道数据集Y、测量矩阵Φ和稀疏基Ψ,计算待恢复信号X的稀疏表达V。 [0062] 从通道数据集Y中沿慢时间域提取测量信号y。将上述公式1代入公式2,得测量信 号γ=φψν= 0V (公式3)。根据公式3即可计算获得待恢复信号X的稀疏表达V。
[0063] 龙魅ν的讨趕可以利用1-范数最小求解方法。
Figure CN104688271BD00071
其中 I ly-Θ ν| |2<ε
[0065] 其中ε为容许误差。1-范数最小的求解方法有很多。优选地,采用基追踪算法求解。 该方法可以使得测量信号y的计算结果更准确,从而获得更准确的重建图像。可选地,采用 正交匹配追踪求解。1-范数最小的优化是对0-范数最小问题(稀疏问题)的良好近似,利用 1-范数最小的优化方法可以高概率的解出高难度的0-范数最小问题。由此,可以确保所生 成的超声图像的分辨率更高。
[0066] S332,利用稀疏基ψ对待恢复信号X的稀疏表达ν进行稀疏反变换,从而恢复出合 成聚焦通道数据集X。
[0067] 根据待恢复信号X的稀疏表达ν和稀疏基Ψ,代入公式(1) χ= Ψν即得所恢复信号 X。所有所恢复信号X即组成了合成聚焦通道数据集X。
[0068] 可选地,在计算待恢复信号的稀疏表达之前,从通道数据集Υ沿慢时间域提取的测 量信号y被归一化,并且,在稀疏反变换之后,对稀疏反变换的结果进行量值恢复。可选地, 先将测量信号y归一化,并保存最大值y max,设稀疏基为sys8小波基,ε为le-3,将上述稀疏反 变换获得的结果乘上ymax的实部即为所恢复信号X。如此操作,即使信号强度非常小,也能够 保证正确求解。对于信号强度非常大的情况,能够将容许误差ε合理化。归一化的操作更好 的保证了信号的准确重建,使成像效果更稳定。
[0069] 可选地,对于通道数据集Υ中每个矩阵的第m行、第η列的元素全部为0的情况,将合 成聚焦通道数据集X的每个矩阵的第m行、第η列的元素全部直接置为0。这样可以简化计算, 提尚效率。
[0070] 根据本发明另一方面,还提供了 一种合成聚焦超声成像装置。图7示出了根据本发 明一个具体实施例的合成聚焦超声成像装置700的框图。如图7所示,该装置700包括激励发 射模块710、激励采集模块720、数据恢复模块730和图像重建模块740。
[0071] 其中,激励发射模块710用于激励超声探头的多个阵元多次发射平面波,其中该多 个阵兀每次发射平面波时的发射变迹与其兀素成随机分布的测量矩阵中的相应行相对应。 激励采集模块720用于每发射一次平面波后,激励超声探头的所有阵元接收回波信号,以获 得通道数据。数据恢复模块730用于根据通道数据集和测量矩阵,利用压缩感知重建算法恢 复合成聚焦通道数据集。图像重建模块740用于对合成聚焦通道数据集进行波束合成,以生 成超声波图像。
[0072] 可选地,上述随机分布是0至1的随机分布。
[0073]可选地,上述数据恢复模块730进一步包括第一计算模块和第二计算模块。第一计 算模块用于根据通道数据集、测量矩阵和稀疏基,利用压缩感知重建算法计算待恢复信号 的稀疏表达。第二计算模块用于利用稀疏基对待恢复信号的稀疏表达进行稀疏反变换,从 而恢复出合成聚焦通道数据集。
[0074] 本领域普通技术人员通过阅读上文关于合成聚焦超声成像方法的描述能够理解 该合成聚焦超声成像装置的构成、实现及其技术效果,因此为了简洁,这里不再赘述。
[0075] 本发明已经通过上述实施例进行了说明,但应当理解的是,上述实施例只是用于 举例和说明的目的,而非意在将本发明限制于所描述的实施例范围内。此外本领域技术人 员可以理解的是,本发明并不局限于上述实施例,根据本发明的教导还可以做出更多种的 变型和修改,这些变型和修改均落在本发明所要求保护的范围以内。本发明的保护范围由 附属的权利要求书及其等效范围所界定。

Claims (10)

1. 一种合成聚焦超声成像方法,包括: 激励超声探头的多个阵元多次发射平面波,其中所述多个阵元每次发射平面波时的发 射变迹与其元素成随机分布的测量矩阵中的相应行相对应; 每发射一次平面波后,激励所述超声探头的所有阵元接收回波信号,以获得通道数据; 根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢复合成聚焦通道数据集; 以及 对所述合成聚焦通道数据集进行波束合成,以生成超声波图像。
2. 如权利要求1所述的方法,其特征在于,所述随机分布是0至1的随机分布。
3. 如权利要求1或2所述的方法,其特征在于,所述恢复合成聚焦通道数据集的步骤进 一步包括: 根据所述通道数据集、所述测量矩阵和稀疏基,利用所述压缩感知重建算法计算待恢 复信号的稀疏表达;以及 利用所述稀疏基对所述待恢复信号的稀疏表达进行稀疏反变换,从而恢复出所述合成 聚焦通道数据集。
4. 如权利要求3所述的方法,其特征在于,所述稀疏基是小波基。
5. 如权利要求3所述的方法,其特征在于,所述计算待恢复信号的稀疏表达是利用1-范 数最小求解方法。
6. 如权利要求3所述的方法,其特征在于,所述计算待恢复信号的稀疏表达之前,从所 述通道数据集沿慢时间域提取的测量信号被归一化;并且,在所述稀疏反变换之后,对所述 稀疏反变换的结果进行量值恢复。
7. 如权利要求1或2所述的方法,其特征在于,对于通道数据集中每个矩阵的第m行、第η 列的元素全部为0的情况,将所述合成聚焦通道数据集的每个矩阵的第m行、第η列的元素全 部直接置为0。
8. -种合成聚焦超声成像装置,包括: 激励发射模块,用于激励超声探头的多个阵元多次发射平面波,其中所述多个阵元每 次发射平面波时的发射变迹与其兀素成随机分布的测量矩阵中的相应行相对应; 激励采集模块,用于每发射一次平面波后,激励所述超声探头的所有阵元接收回波信 号,以获得通道数据; 数据恢复模块,用于根据通道数据集和所述测量矩阵,利用压缩感知重建算法恢复合 成聚焦通道数据集;以及 图像重建模块,用于对所述合成聚焦通道数据集进行波束合成,以生成超声波图像。
9. 如权利要求8所述的装置,其特征在于,所述随机分布是0至1的随机分布。
10. 如权利要求8所述的装置,其特征在于,所述数据恢复模块进一步包括: 第一计算模块,用于根据所述通道数据集、所述测量矩阵和稀疏基,利用所述压缩感知 重建算法计算待恢复信号的稀疏表达;以及 第二计算模块,用于利用所述稀疏基对所述待恢复信号的稀疏表达进行稀疏反变换, 从而恢复出所述合成聚焦通道数据集。
CN201510142749.4A 2015-03-27 2015-03-27 合成聚焦超声成像方法和装置 Active CN104688271B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510142749.4A CN104688271B (zh) 2015-03-27 2015-03-27 合成聚焦超声成像方法和装置

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN201510142749.4A CN104688271B (zh) 2015-03-27 2015-03-27 合成聚焦超声成像方法和装置
US15/526,305 US20170336500A1 (en) 2015-03-27 2015-09-06 Method and device for ultrasonic imaging by synthetic focusing
PCT/CN2015/089005 WO2016155239A1 (zh) 2015-03-27 2015-09-06 合成聚焦超声成像方法和装置
EP15887196.2A EP3187113B1 (en) 2015-03-27 2015-09-06 Method and device for ultrasonic imaging by synthetic focusing

Publications (2)

Publication Number Publication Date
CN104688271A CN104688271A (zh) 2015-06-10
CN104688271B true CN104688271B (zh) 2017-04-26

Family

ID=53336135

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510142749.4A Active CN104688271B (zh) 2015-03-27 2015-03-27 合成聚焦超声成像方法和装置

Country Status (4)

Country Link
US (1) US20170336500A1 (zh)
EP (1) EP3187113B1 (zh)
CN (1) CN104688271B (zh)
WO (1) WO2016155239A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111134719A (zh) * 2019-12-19 2020-05-12 西安交通大学 一种聚焦超声辐照相变纳米液滴的主被动超声复合成像方法及系统

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104688271B (zh) * 2015-03-27 2017-04-26 清华大学 合成聚焦超声成像方法和装置
CN111278363A (zh) * 2017-10-16 2020-06-12 北京深迈瑞医疗电子技术研究院有限公司 超声成像设备、系统及其超声造影成像的图像增强方法
US11194014B1 (en) 2018-02-22 2021-12-07 United States Of America As Represented By The Secretary Of The Air Force System, method and apparatus for recovering polarization radar data
CN108693252B (zh) * 2018-03-31 2021-04-16 大连交通大学 三角矩阵聚焦成像的工件探伤检测方法
CN109363714B (zh) * 2018-09-17 2020-12-15 深圳迈瑞生物医疗电子股份有限公司 一种超声成像设备及其超声成像方法
CN109350112B (zh) * 2018-11-13 2020-06-12 飞依诺科技(苏州)有限公司 超声图像的成像方法、装置及医疗设备
CN109766646B (zh) * 2019-01-16 2021-06-04 北京大学 一种基于稀疏通道回波数据重建的超声成像方法及装置
CN113509208B (zh) * 2021-09-14 2021-11-30 西南石油大学 一种基于相位约束的超高速超声成像的重建方法

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6135961A (en) * 1996-06-28 2000-10-24 Sonosite, Inc. Ultrasonic signal processor for a hand held ultrasonic diagnostic instrument
US6674687B2 (en) * 2002-01-25 2004-01-06 Navcom Technology, Inc. System and method for navigation using two-way ultrasonic positioning
US6679847B1 (en) * 2002-04-30 2004-01-20 Koninklijke Philips Electronics N.V. Synthetically focused ultrasonic diagnostic imaging system for tissue and flow imaging
EP1788940A4 (en) * 2004-09-10 2009-09-09 Ind Res Ltd SYNTHETIC FOCUSING PROCESS
US7173555B2 (en) * 2004-11-16 2007-02-06 Gmr Research & Technology, Inc. Method and system for nonlinear and affine signal processing
US7271747B2 (en) * 2005-05-10 2007-09-18 Rice University Method and apparatus for distributed compressed sensing
WO2009058915A1 (en) * 2007-10-29 2009-05-07 The Trustees Of The University Of Pennsylvania Computer assisted diagnosis (cad) of cancer using multi-functional, multi-modal in-vivo magnetic resonance spectroscopy (mrs) and imaging (mri)
CN105182263A (zh) * 2008-04-28 2015-12-23 康奈尔大学 分子mri中的磁敏度精确量化
US8693597B2 (en) * 2008-12-12 2014-04-08 Blackberry Limited Sensor-based wireless communication systems using compressive sampling
US7834795B1 (en) * 2009-05-28 2010-11-16 Bae Systems Information And Electronic Systems Integration Inc. Compressive sensor array system and method
US9638798B2 (en) * 2010-09-20 2017-05-02 B-K Medical Aps Imaging transducer array
CN102028499B (zh) * 2010-12-24 2012-10-31 飞依诺科技(苏州)有限公司 超声探头自适应成像系统
US9172476B2 (en) * 2012-03-09 2015-10-27 The United States Of America As Represented By The Secretary Of The Army Method and system for removal of noise in signal
US8861588B2 (en) * 2011-04-04 2014-10-14 The United States Of America As Represented By The Secretary Of The Army Apparatus and method for sampling and reconstruction of wide bandwidth signals below Nyquist rate
US8824544B2 (en) * 2012-03-09 2014-09-02 The United States Of America As Represented By The Secretary Of The Army Method and system for recovery of missing spectral information in wideband signal
CN102288968B (zh) * 2011-05-17 2014-03-12 西安电子科技大学 小孔径高分辨相控阵超声探测成像方法
CN102279394B (zh) * 2011-05-17 2013-06-26 西安电子科技大学 低频窄带高分辨超声探测成像方法
CN102435992B (zh) * 2011-09-26 2013-10-23 重庆博恩富克医疗设备有限公司 基于广义相干系数的合成聚焦成像方法
EP2785253A4 (en) * 2011-12-01 2015-07-15 Maui Imaging Inc MOTION DETECTION USING AN OPENING DOPPLER BASED ON ACOUSTIC PULSE EMITTED AND MULTIPLE OPENINGS
US9857463B2 (en) * 2012-11-12 2018-01-02 Sony Corporation Radar apparatus and method
CN102920478B (zh) * 2012-11-22 2014-04-16 山东大学 一种合成聚焦的便携式b型超声成像方法
CN103679762A (zh) * 2013-12-13 2014-03-26 江苏大学 一种基于稀疏数据的超声信号重建方法
US11125866B2 (en) * 2015-06-04 2021-09-21 Chikayoshi Sumi Measurement and imaging instruments and beamforming method
US10349917B2 (en) * 2014-06-11 2019-07-16 The Johns Hopkins University Synthetic aperture ultrasound system
CN104318619B (zh) * 2014-10-20 2017-07-18 西北工业大学 面向无损检测的自适应压缩感知的重建方法
CN104306022B (zh) * 2014-10-24 2016-05-25 西安电子科技大学 用gpu实现压缩感知超声成像的方法
CN104306023B (zh) * 2014-10-24 2016-05-25 西安电子科技大学 基于压缩感知的超声成像快速实现方法
CN104688271B (zh) * 2015-03-27 2017-04-26 清华大学 合成聚焦超声成像方法和装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111134719A (zh) * 2019-12-19 2020-05-12 西安交通大学 一种聚焦超声辐照相变纳米液滴的主被动超声复合成像方法及系统
CN111134719B (zh) * 2019-12-19 2021-01-19 西安交通大学 一种聚焦超声辐照相变纳米液滴的主被动超声复合成像方法及系统

Also Published As

Publication number Publication date
EP3187113B1 (en) 2020-12-16
EP3187113A1 (en) 2017-07-05
US20170336500A1 (en) 2017-11-23
EP3187113A4 (en) 2017-11-22
CN104688271A (zh) 2015-06-10
WO2016155239A1 (zh) 2016-10-06

Similar Documents

Publication Publication Date Title
CN104688271B (zh) 合成聚焦超声成像方法和装置
CN102920478B (zh) 一种合成聚焦的便携式b型超声成像方法
CN102551791B (zh) 一种超声成像方法和装置
CN101190134A (zh) 超声波诊断系统中的多波束发射和接收方法及其装置
Burshtein et al. Sub-Nyquist sampling and Fourier domain beamforming in volumetric ultrasound imaging
CN104777484B (zh) 压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统
CN104739448A (zh) 一种超声成像方法及装置
Wagner et al. Xampling in ultrasound imaging
US20210333379A1 (en) Ultrasound system for high-speed and high resolution imaging applications
Besson et al. Ultrafast ultrasound imaging as an inverse problem: Matrix-free sparse image reconstruction
CN106680825B (zh) 一种声学阵列成像系统与方法
Madore et al. Accelerated focused ultrasound imaging
US9384530B2 (en) Enhanced ultrasound image formation using qualified regions of overlapping transmit beams
Liu et al. Compressed sensing based synthetic transmit aperture imaging: Validation in a convex array configuration
US20160061950A1 (en) Synthetic Aperture Imaging Methods And Systems
Hemmsen et al. Tissue harmonic synthetic aperture ultrasound imaging
US20100191115A1 (en) Ultrasound imaging system and method
CN111134719B (zh) 一种聚焦超声辐照相变纳米液滴的主被动超声复合成像方法及系统
Wang et al. An easily-achieved time-domain beamformer for ultrafast ultrasound imaging based on compressive sensing
Ramkumar et al. Compressed Sensing Approach for Reducing the Number of Receive Elements in Synthetic Transmit Aperture Imaging
US8235906B2 (en) System and method for accelerated focused ultrasound imaging
Liu et al. Compressed sensing for high frame rate, high resolution and high contrast ultrasound imaging
CN103565469B (zh) 一种提高二维阵列实时超声体成像质量的方法
CN107003403A (zh) 声成像的方法和设备
CN106037805A (zh) 超声成像的方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant