CN108135566A - 一种增强超声造影图像的方法、系统及超声造影成像设备 - Google Patents

一种增强超声造影图像的方法、系统及超声造影成像设备 Download PDF

Info

Publication number
CN108135566A
CN108135566A CN201680039006.7A CN201680039006A CN108135566A CN 108135566 A CN108135566 A CN 108135566A CN 201680039006 A CN201680039006 A CN 201680039006A CN 108135566 A CN108135566 A CN 108135566A
Authority
CN
China
Prior art keywords
information
threshold
weight coefficient
contrast
radiography
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.)
Granted
Application number
CN201680039006.7A
Other languages
English (en)
Other versions
CN108135566B (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.)
Chinese PLA General Hospital
Beijing Shen Mindray Medical Electronics Technology Research Institute Co Ltd
Original Assignee
Chinese PLA General Hospital
Beijing Shen Mindray Medical Electronics Technology Research Institute Co Ltd
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 Chinese PLA General Hospital, Beijing Shen Mindray Medical Electronics Technology Research Institute Co Ltd filed Critical Chinese PLA General Hospital
Publication of CN108135566A publication Critical patent/CN108135566A/zh
Application granted granted Critical
Publication of CN108135566B publication Critical patent/CN108135566B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Theoretical Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)

Abstract

一种增强超声造影图像的方法、系统和超声造影成像设备,根据从超声回波信号中获得的非线性的造影信息和线性的组织信息,计算非线性参量,得到非线性参量分布图,分别采用两极阈值对非线性参量分布图进行分割,每次划分都分出造影剂区域和组织残留区域,对造影剂区域进行增强处理,对组织残留区域抑制处理,从而可实现对绝对属于造影剂区域进行大程度的增强,对绝对属于组织残留的区域进行大程度的抑制,对介于造影剂和组织残留的区域增强和抑制的程度稍微弱些,最后再将两次分类和权重调整后的结果进行融合处理,使对造影剂的增强、对造影剂的抑制,对介于造影剂和组织残留的保守处理,都能达到比较理想的造影图像增强效果。

Description

说明书
发明名称:一种增强超声造影图像的方法、 系统及超声造影成像设 备
技术领域
[0001] 本发明涉及超声造影成像设备, 具体涉及增强超声造影图像的方法及系统。
[0002]
[0003] 背景技术
[0004] 超声造影成像与常规超声成像相比, 能增强超声回波信号强度, 显示常规超声 无法看见的细小血管, 随着临床对超声检査的要求越来越高, 超声造影成像越 来越多地被临床应用, 且用来检査器官也越来越广, 其在临床诊断中的地位越 来越高。
[0005] 超声造影成像是注入造影剂后再超声成像, 是一种增强超声成像, 造影剂是一 种能增强超声回波信号强度的增强剂。 超声造影剂的种类多样, 应用较多的是 一种大小和红细胞接近的空心球体样的微泡, 由于大小和红细胞接近所以微泡 能到达全身所有血管。 由于是空心球样所以微泡的声阻抗差大, 极大的增加反 射回波信号强度, 使微小血管或毛细血管能显像。 大多情况下造影剂微泡通过 静脉注入, 随血流分布到全身, 造影微泡与红细胞行为相似, 检测或者跟踪微 泡, 能反映出正常和异常组织的血流灌注情况, 超声造影成像技术也可以称之 为检测微泡信号技术。
[0006] 衡量造影图像的好坏, 一般会从以下几个指标考虑, SNR(signal noise
ratio) CTR(contrast to tissue ratio)、 吋间分辨力、 空间分辨力。 SNR也可以理解 为灵敏度或敏感性, 就是对微泡的检测程度; CTR是造影图像特有的特征, 是造 影图像中造影剂信号和组织信号强度之比。 理论上造影图像只是对微泡显影, 没有注入造影剂吋造影图像是没有内容的, 但是实际情况下由于系统的相消性 限制, 即使在没有注入造影剂吋造影图像也有少量显影, 该少量显影称为组织 残留, 这些组织残留在造影剂微泡灌注后也存在, 与微泡信号混在一起, 影响 用户对造影剂微泡的识别判断。 CTR越大, 组织残留越小, 微泡强度和组织残留 强度的差距越大, 纯造影剂信号的动态范围更大, 造影图像的层次更丰富, 微 泡的灌注和效果过程显示的更清楚, 图像的对比分辨力也更好。 吋间分辨力和 空间分辨力和常规超声图像的要求一样, 分辨力越高越好。
[0007] 随着超声造影的使用越来越广, 临床对造影图像的要求也越来越高, 不再是能 看见微泡显影, 而是希望能显示更多细节。 这就要求提高造影图像质量以能够 显示更多的细节。
[0008] 在超声造影图像中, 提高造影图像质量的关键是提高 SNR和 CTR。 SNR是由探 头灵敏度和系统平台决定的, SNR高, 意味着能检测出更微小的微泡, 而微泡越 小信号幅度越小, 有些小微泡的信号幅度和组织残留幅度相当甚至更小, 因此 , 即使 SNR很高, 但如果 CTR不够高, 这些小微泡反映出的信号将淹没在组织残 留信息中, 在造影图像中也观察不到。 因此提高造影图像质量最关键的是要尽 量提高 CTR。
[0009]
[0010] 发明内容
[0011] 本发明要解决的技术问题是提供另一种对造影信息进行增强的方案, 从而提高 造影图像的 CTR。
[0012] 根据第一方面, 一种实施例中提供一种增强超声造影图像的方法,包括:
[0013] 获取信息, 根据造影成像过程中的超声回波信号获得非线性的造影信息和线性 的组织信息;
[0014] 计算非线性参量, 比较所述造影信息和组织信息得到一帧造影信息的非线性参 量分布图, 所述非线性参量是衡量超声波在媒质中传播吋产生的非线性效应大 小的参量;
[0015] 第一阈值分割, 采用第一阈值对非线性参量分布图进行分割, 小于第一阈值的 区域为组织残留区域, 大于第一阈值的区域为造影剂区域, 所述第一阈值使得 非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧中对应的 点属于组织残留信息;
[0016] 形成第一图像, 对该造影信息帧中属于组织残留区域的信息进行抑制处理, 对 该造影信息帧中属于造影剂区域的信息进行增强处理, 并基于处理后的结果形 成第一图像数据;
[0017] 第二阈值分割, 采用第二阈值对非线性参量分布图进行分割, 小于第二阈值的 区域为组织残留区域, 大于第二阈值的区域为造影剂区域, 所述第二阈值使得 非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧中对应的 点属于造影剂信息, 所述第二阈值大于第一阈值;
[0018] 形成第二图像, 对该造影信息帧中属于组织残留区域的信息进行抑制处理, 对 该造影信息帧中属于造影剂区域的信息进行增强处理, 并基于处理后的结果形 成第二图像数据;
[0019] 图像融合, 将第一图像数据和第二图像数据进行融合, 形成增强后的造影图像 数据。
[0020] 根据第二方面, 一种实施例中提供一种超声造影图像增强系统,包括:
[0021] 信息获取单元, 用于根据造影成像过程中的超声回波信号获得非线性的造影信 息和线性的组织信息;
[0022] 非线性参量计算单元, 用于比较所述造影信息和组织信息得到一帧造影信息的 非线性参量分布图, 所述非线性参量是衡量超声波在媒质中传播吋产生的非线 性效应大小的参量;
[0023] 第一阈值分割单元, 用于采用第一阈值对非线性参量分布图进行分割, 小于第 一阈值的区域为组织残留区域, 大于第一阈值的区域为造影剂区域, 所述第一 阈值使得非线性参量分布图中小于第一阈值的区域中绝大多数在该造影信息帧 中对应的点属于组织残留信息;
[0024] 第一图像形成单元, 用于对该造影信息帧中属于组织残留区域的信息进行抑制 处理, 对该造影信息帧中属于造影剂区域的信息进行增强处理, 并基于处理后 的结果形成第一图像数据;
[0025] 第二阈值分割单元, 用于采用第二阈值对非线性参量分布图进行分割, 小于第 二阈值的区域为组织残留区域, 大于第二阈值的区域为造影剂区域, 所述第二 阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息帧 中对应的点属于造影剂信息, 所述第二阈值大于第一阈值;
[0026] 第二图像形成单元, 用于对该造影信息帧中属于组织残留区域的信息进行抑制 处理, 对该造影信息帧中属于造影剂区域的信息进行增强处理, 并基于处理后 的结果形成第二图像数据;
[0027] 图像融合单元, 用于将第一图像数据和第二图像数据进行融合, 形成增强后的
[0028] 根据第三方面, 一种实施例中提供一种超声造影成像设备, 包括:
[0029] 探头, 用于向组织发射超声波并接收组织反射回的超声回波;
[0030] 前信号处理模块, 对超声回波信号进行处理, 从超声回波信号中提取反映造影 剂微泡信息的非线性成分和反映组织解剖特征的线性成分;
[0031] 第一信号处理模块, 用于对线性成分进行处理;
[0032] 组织图像生成模块, 用于将第一信号处理模块的输出处理成组织图像数据; [0033] 第二信号处理模块, 用于对非线性成分进行处理;
[0034] 造影图像生成模块, 用于将第二信号处理模块的输出处理成造影图像数据; [0035] 用于对造影信息进行增强的超声造影图像增强系统, 超声造影图像增强系统如 上述所述, 以得到增强后的造影图像数据;
[0036] 显示模块, 用于对组织图像数据和增强后的造影图像数据进行可视化显示。
[0037] 本发明实施例中, 对造影信息帧的非线性参量分布图采用两极阈值分割, 每次 划分都分出造影剂区域和组织残留区域, 对造影剂区域进行增强处理, 对组织 残留区域抑制处理, 从而可实现对绝对属于造影剂区域进行大程度的增强, 对 绝对属于组织残留的区域进行大程度的抑制, 对介于造影剂和组织残留的区域 增强和抑制的程度稍微弱些, 最后再将两次分类和权重调整后的结果进行融合 处理, 使对造影剂的增强、 对造影剂的抑制, 对介于造影剂和组织残留的保守 处理, 都能达到比较理想的造影图像增强效果。
[0038]
[0039] 附图说明
[0040] 图 1为超声造影成像设备的第一种结构示意图;
[0041] 图 2为超声造影成像设备的第二种结构示意图;
[0042] 图 3为超声造影成像设备的第三种结构示意图;
[0043] 图 4为一种实施例中超声造影图像增强系统的结构示意图; [0044] 图 5为一种实施例中超声造影图像增强系统的流程图;
[0045] 图 6为一种实施例中采用直方图方式表示的非线性参量分布图;
[0046] 图 7为另一种实施例中采用二维方式表示的非线性参量分布图。
[0047]
[0048] 具体实施方式
[0049] 下面通过具体实施方式结合附图对本发明作进一步详细说明。
[0050] 如图 1所示, 超声造影成像设备包括探头 10、 前信号处理模块 13、 第一信号处 理模块 16、 第二信号处理模块 17、 组织图像生成模块 18、 造影图像生成模块 19
、 超声造影图像增强系统 20和显示模块 100。
[0051] 探头 10用于向人体组织 11发射超声波并接收组织 11反射回的超声回波, 超声回 波可存在信号接收模块 12中。
[0052] 前信号处理模块 13用于对超声回波信号进行处理, 提取有用的信息。 前信号处 理模块 13主要包含放大器、 波束合成器和滤波器, 提取更多有用的、 不同成分 的信号。 在超声造影成像设备中, 前信号处理模块 13需要提取两种不同的信号
, 一种是线性成分 14, 反映组织解剖特征; 一种是非线性成分 15, 反映造影剂 微泡信息。
[0053] 第一信号处理模块 16用于对线性成分进行处理, 第二信号处理模块 17用于对非 线性成分进行处理。 第一信号处理模块 16和第二信号处理模块 17的处理类似, 主要包括解调、 滤波、 降采样、 求包络、 求 Log、 动态范围变换等, 只是参数的 选择不同。
[0054] 组织图像生成模块 18用于将第一信号处理模块 16的输出处理成组织图像数据; 造影图像生成模块 19用于将第二信号处理模块 17的输出处理成造影图像数据。
[0055] 超声造影图像增强系统 20用于对造影信息进行增强, 以得到增强后的造影图像 数据; 显示模块 100用于对组织图像数据和增强后的造影图像数据进行可视化显 示。 大多数造影图像旁会同吋显示一幅组织图像, 用于定位病灶位置。 为了得 到微泡的非线性特征, 造影成像通常会在同一位置发射多次相位或者幅度不同 的脉冲波, 组织和造影图像都是用这多次发射波形中的部分进行成像, 所以组 织图像和造影图像显示的是同一吋刻同一位置的信息, 组织图像显示的是回波 的线性特性, 造影图像显示的是回波中的非线性特性。
[0056] 在采用造影增强技术对造影图像进行增强吋, 需要通过阈值将图像划分成组织 残留区域和造影剂区域, 然后对划归为造影剂区域的点进行增强, 但在对造影 增强的继续研究中, 发明人认识到组织残留区域和造影剂区域并没有明确的边 界, 任何阈值都无法正确的区分组织残留区域和造影剂区域。 因此, 在本发明 实施例中, 首先得到组织的线性信息和造影剂微泡的非线性信息, 再比较非线 性信息和线性信息得到非线性参量值, 接着根据非线性参量值对造影图像进行 两次分类, 每次都分成组织残留和造影剂两类, 对两类按照不同的权重增强或 抑制, 根据分类的不同, 权重系数也不同; 最后将两次分类和权重调整后的结 果进行融合处理。
[0057] 在具体实施例中, 超声造影图像增强系统 20可以对组织图像生成模块 18和造影 图像生成模块 19输出的图像域的造影信息进行增强处理, 如图 1所示, 超声造影 图像增强系统 18的输入端分别连接组织图像生成模块 18和造影图像生成模块 19 , 输出端连接到显示模块 100, 根据组织图像生成模块 18和造影图像生成模块 19 输出的组织图像数据和造影图像数据对造影图像数据进行增强处理, 并将增强 结果输出到显示模块 100。 超声造影图像增强系统 20还可以对前信号处理模块 13 输出的射频域的造影信息进行增强处理。 如图 2所示, 超声造影图像增强系统 20 连接在前信号处理模块 13和第二信号处理模块 17之间, 接收前信号处理模块 13 输出的非线性成分 15和线性成分 14, 根据非线性成分 15和线性成分 14对非线性 成分进行增强处理, 并将增强结果输出到第二信号处理模块 17。 超声造影图像 增强系统 18还可以对第二信号处理模块中任一环节的造影信息进行增强处理。 如图 3所示, 超声造影图像增强系统 20的输入端分别连接第一信号处理模块 16和 第二信号处理模块 17的任一环节, 例如解调后信号, 求包络后, Log变换后等, 超声造影图像增强系统 20的输出端连接造影图像生成模块 19。
[0058] 超声造影图像增强系统 20的结构图如图 4所示, 包括信息获取单元 201、 非线性 参量计算单元 202、 第一阈值分割单元 203、 第一图像形成单元 204、 第二阈值分 割单元 205、 第二图像形成单元 206和图像融合单元 207。 信息获取单元 201用于 根据造影成像过程中的超声回波信号获得非线性的造影信息和线性的组织信息 , 即造影信息和线性的组织信息是从同一超声回波信号中提取出的两个信息, 非线性信息为造影信息, 线性信息为组织残留信息; 非线性参量计算单元 202用 于比较同一吋刻的造影信息和组织信息得到一帧造影信息的非线性参量分布图 , 该非线性参量分布图称为该造影信息帧的非线性参量分布图, 所述非线性参 量是衡量超声波在媒质中传播吋产生的非线性效应大小的参量; 本实施例中, 非线性参量分布图根据同一吋刻的一帧非线性造影信息和一帧线性组织信息的 比值或差值计算得出。 第一阈值分割单元 203用于采用第一阈值对非线性参量分 布图进行分割, 小于第一阈值的区域为组织残留区域, 大于第一阈值的区域为 造影剂区域, 所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝 大多数在该造影信息帧中对应的点属于组织残留信息; 第一图像形成单元 204用 于对该造影信息帧中属于组织残留区域的信息进行抑制处理, 对该造影信息帧 中属于造影剂区域信息进行增强处理, 并基于处理后的结果形成第一图像数据 ; 第二阈值分割单元 205用于采用第二阈值对非线性参量分布图进行分割, 小于 第二阈值的区域为组织残留区域, 大于第二阈值的区域为造影剂区域, 所述第 二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多数在该造影信息 帧中对应的点属于造影剂信息, 所述第二阈值大于第一阈值; 第二图像形成单 元 206用于对该造影信息帧中属于组织残留区域的信息进行抑制处理, 对该造影 信息帧中属于造影剂区域的信息进行增强处理, 并基于处理后的结果形成第二 图像数据; 图像融合单元 207用于将第一图像数据和第二图像数据进行融合, 形 成增强后的造影图像数据。
[0059] 在另一实施例中, 超声造影图像增强系统 20还包括噪声识别单元 208, 噪声识 别单元 208用于在非线性参量计算单元计算非线性参量之前判断造影信息中的噪 声, 将一帧造影信息分为噪声信号和除去噪声的非噪声造影信号; 非线性参量 计算单元 202在计算非线性参量吋采用非噪声造影信号计算非线性参量。
[0060] 超声造影图像增强系统的具体处理流程如图 5所示, 包括以下步骤:
[0061] 一、 去噪
[0062] 去噪的目的是还原信号自身大小, 避免 TGC (信号增益补偿) 及噪声对后面处 理带来影响。 一般超声系统中用 ATGC (analog time gain compensation,前端模拟 信号增益补偿) 、 DTGC (digital time gain compensation,后端数字域增益补偿) 用来补偿信号传播过程中深度方向的衰减来补偿信号沿传播方向的信号衰减, 如果没有 ATGC和 DTGC的作用, 探头空打吋的噪声是近似为白噪声, 加上了 AT GC和 DTGC的影响后, 噪声均值曲线是随深度递增的曲线。 为了得到信号真实 的大小, 减去噪声均值曲线不仅仅是去噪, 也消除了 ATGC和 DTGC对信号的影 响。
[0063] 步骤 21到 24是实现去噪。 去噪首先要获得噪声信号, 本实施例中获得噪声信号 的方法是暂停发射, 在超声造影图像增强系统 20前的环节执行步骤 21和 22, 缓 存一帧非线性造影噪声信号和一帧线性组织噪声信号, 即控制探头不发射, 然 后根据回波得到一帧非线性造影噪声信号和一帧线性组织噪声信号。 处理的具 体步骤是计算一帧噪声信号每行的均值得到噪声纵向均值曲线, 再经过低通滤 波得到光滑的噪声均值曲线。
[0064] 在步骤 23中, 将超声造影图像增强系统 20输入的造影信息减去非线性造影噪声 信号的噪声均值曲线, 得到去噪后的造影信号。
[0065] 在步骤 24中, 将超声造影图像增强系统 20输入的组织信息减去线性组织噪声信 号的噪声均值曲线, 得到去噪后的组织信号。
[0066] 二、 检测噪声
[0067] 步骤 25, 检测去噪后造影图像中的噪声, 将一帧造影图像划分成噪声区域 26和 非噪声造影剂区域 27。 噪声检测可用阈值分割的方法, 小于阈值的就是噪声, 阈值是根据系统特性经验确定的。
[0068] 三、 计算非线性参量
[0069] 本实施例中, 非线性参量是衡量超声波在媒质中传播吋产生的非线性效应大小 的重要参量。
[0070] 步骤 28中, 比较同一吋刻的一帧造影信息和一帧组织信息, 得到该造影信息帧 的非线性参量分布图。 理论上, 非线性参量为媒质物态方程泰勒展幵式中多次 项系数与线性系数的比值或差值, 例如, 媒质物态方程泰勒展幵式中, 线性系 数为 A, 二次项系数为 B, 三次项系数为 C, 非线性参量可以是 B/A、 C/A或 (B+C) /A, 及 C其他组合与 A的比值。 本实施例中, 非线性参量通过计算同一吋刻 造影信号和组织信号的比值或差值得出。 对于一帧造影信息, 其由 MxN个像素 值组成, 每个像素值对应显示界面上的一个像素点, 该像素值可能属于造影剂 产生的信号, 也可能属于组织残留信号。 对于一帧造影信息而言, 每个像素点 都会有对应的一个非线性参量, 如果该像素值是组织残留信息, 则其对应的非 线性参量值较小, 如果该像素值是造影剂信息, 则其对应的非线性参量值较大 , 因此, 通过计算造影信息帧中各像素值的非线性参量, 即可将属于组织残留 信息的像素值和属于造影剂信息的像素值区分幵。 根据同一吋刻的一帧非线性 造影信息和一帧线性组织信息的非线性参量即可得出该造影信息帧中各像素点 的非线性参量, 对这些非线性参量进行统计, 则可得到该造影信息帧的非线性 参量分布图, 得出的非线性参量分布图可以采用直方图方式表示, 也可以采用 二维图方式表示。
[0071] 四、 两极阈值划分
[0072] 如图 6所示为采用直方图方式表示的非线性参量分布图, 包括一个大峰和一个 小峰, 直方图的横坐标就是计算出的非线性参量值, 单位是 dB, 纵坐标表示造 影信息帧中与该非线性参量值对应的点的数量。 经试验证明, 大峰区域中绝大 部分区域所对应的造影信息帧中点反映的是造影剂微泡信息, 小峰区域中绝大 部分区域所对应的造影信息帧中点反映的是组织残留信息。
[0073] 在确定两极阈值的另一实施例中, 非线性参量分布图为二维图, 如图 7所示, 其以色彩表示造影信息帧中各点非线性参量值, 深蓝色的表示噪声区域, 非线 性参量值由大到小采用的颜色表示为红、 黄、 绿、 蓝, 因组织残留的非线性参 量值较小, 造影剂的非线性声参量值较大, 故可以确定约小于 -40dB的一定是组 织残留, 约大于 -20dB的一定是造影剂, 这里 -40dB就是较小阈值, -20dB就是较 大阈值。 图 7只是一个具体的例子, 实际中, 需要根据系统平台, 通过较多非线 性参量二维图, 找出确定为组织残留的点的非线性参量值大小, 其最大的值就 是较小阈值, 找出确定为造影剂的点的非线性参量值的大小, 其最小值就是较 大阈值。
[0074] 本实施例中, 通过两极阈值来划分非线性参量分布图, 每次划分都分为造影剂 区域和组织残留区域, 对造影剂区域进行增强显示, 对组织残留区域抑制显示 。 两极阈值就是极小、 极大两个阈值, 非线性参量不同, 阈值的选择也不同, 阈值的大小和系统的相消性有关, 大阈值要大于系统的相消值, 小阈值要小于 系统的相消值, 相消值简单讲就是注入造影剂最强吋刻和没有注入造影剂吋两 个非线性成分大小的差值。 具体如下:
[0075] 在步骤 29中, 采用第一阈值对非线性参量分布图进行分割, 如图 6中, 左边的 直线阈值为第一阈值, 该第一阈值使得非线性参量分布图中小于第一阈值的区 域中绝大多数在该造影信息帧中对应的点属于组织残留信息, 因此第一阈值也 称为极小阈值。 划分后可认为: 在如图 6所示的非线性参量分布图中, 小于第一 阈值的左边的区域为组织残留区域 32, 组织残留区域 32中的点在造影信息帧中 的对应点被认为属于组织残留信息。 大于第一阈值的右边的区域为造影剂区域 3 1, 造影剂区域 31中的点在造影信息帧中的对应点被认为可能属于造影剂产生的 信息。 此吋的小于第一阈值的几乎都是组织残留, 而大于第一阈值的可能属于 造影剂。 再结合噪声区域 26, 这样就把造影图像分成噪声、 造影剂和组织残留 三个不同的区域。
[0076] 步骤 35, 为了达到增强造影剂、 抑制组织残留的目的, 对该造影信息帧中属于 组织残留区域的信息进行抑制处理, 对该造影信息帧中属于造影剂区域的信息 进行增强处理。 例如对上述三个区域分别乘以或加上不同的权重系数 Pl、 P2和 P 3, 如果是乘以权重系数, 则造影剂区域的系数 P2大于 1, 组织残留区域的系数 P 3应小于 1, 噪声区域系数 P1可以等于 1, 如果想抑制噪声, P1也可以小于 1。 同 吋, P2大于 1, 也增加了造影信号的动态范围。 如果是加上权重系数, 则 P3<0<P 2, 这样能增强造影信号强度, 减弱组织残留强度, P1≤0, P1<0也达到抑制噪声 的目的。 对上述三个区域分别乘以或加上不同的权重系数后再形成一个新的图 像, 称为第一图像数据。
[0077] 在步骤 30中, 采用第二阈值对非线性参量分布图进行分割, 如图 6中, 右边直 线阈值为第二阈值, 该第二阈值使得非线性参量分布图中大于第二阈值的区域 中绝大多数在该造影信息帧中对应的点属于造影剂信息, 因此, 第二阈值也称 为极大阈值, 且第二阈值大于第一阈值。 划分后可认为: 在如图 6所示的非线性 参量分布图中, 小于第二阈值的区域为组织残留区域 34, 组织残留区域 34中的 点在造影信息帧中的对应点被认为可能属于组织残留信息。 大于第二阈值的区 域为造影剂区域 33, 即造影剂区域 33中的点在造影信息帧中的对应点被认为属 于造影剂产生的信息。 此吋的大于第二阈值的几乎都是造影剂, 而小于第二阈 值的则可能属于组织残留。 再结合噪声区域 26, 也将造影图像分成噪声、 造影 剂和组织残留三个不同的区域。
[0078] 步骤 36, 为了达到增强造影剂、 抑制组织残留的目的, 对该造影信息帧中属于 组织残留区域的信息进行抑制处理, 对该造影信息帧中属于造影剂区域的信息 进行增强处理。 例如对步骤 30中的三个区域分别乘以或加上不同的权重系数 P1 、 P4和 P5。 同样如果乘以权重系数, 则造影剂区域的系数 P4大于 1, 组织残留区 域的系数 P5应小于 1, 即 P5<1<P4, 而噪声区域系数 PI等于或小于 1。 如果加上权 重系数, 则 P5<0<P4, P1≤0。 对上述三个区域分别乘以或加上不同的权重系数后 再形成一个新的图像, 称为第二图像数据。
[0079] 为了达到增强造影图像的效果, 对绝对属于造影剂区域的像素点增强的程度大 , 对绝对属于组织残留区域的像素点抑制的程度大, 对介于造影剂和组织残留 的区域增强和抑制的程度稍微弱些。 如果乘以权重系数, 增强和抑制的参数满 足 P4>P2>1, P3<P5<1 ; 如果加上权重系数, 增强和抑制的参数满足 P4>P2>0, P 3<P5<0。
[0080] 五、 图像融合
[0081] 因为任何阈值都无法完全正确的区分造影剂和组织残留, 所以第一图像数据和 第二图像数据都有些极端, 作为最终的增强图像都有些不正确。 因此在接下来 的步骤 37中将两种图像进行融合处理, 融合的方法有很多, 本实施例提出了一 种简单的融合方法, 将第一图像数据乘以权重系数 K, 第二图像数据乘以权重系 数 1-Κ, 0<Κ<1,再组合显示。 为了达到增强造影剂的目的, 1-K≥K, K≤0.5。
[0082] 本实施例最后通过对两个图像数据进行融合处理, 相互补充, 使对造影剂的增 强、 对造影剂的抑制, 对介于造影剂和组织残留的保守处理, 都能达到比较理 想的造影图像增强效果。
[0083]
[0084] 本领域技术人员可以理解, 上述实施方式中各种方法的全部或部分步骤可以通 过程序来指令相关硬件完成, 该程序可以存储于一计算机可读存储介质中, 存 储介质可以包括: 只读存储器、 随机存储器、 磁盘或光盘等。
[0085]
[0086] 以上应用了具体个例对本发明进行阐述, 只是用于帮助理解本发明, 并不用以 限制本发明。 对于本领域的一般技术人员, 依据本发明的思想, 可以对上述具 体实施方式进行变化。
技术问题
问题的解决方案
发明的有益效果

Claims (18)

  1. 权利要求书
  2. [权利要求 1] 一种增强超声造影图像的方法,其特征在于包括:
    获取信息, 根据造影成像过程中的超声回波信号获得非线性的造影信 息和线性的组织信息;
    计算非线性参量, 比较所述造影信息和组织信息得到一帧造影信息的 非线性参量分布图, 所述非线性参量是衡量超声波在媒质中传播吋产 生的非线性效应大小的参量;
    第一阈值分割, 采用第一阈值对非线性参量分布图进行分割, 小于第 一阈值的区域为组织残留区域, 大于第一阈值的区域为造影剂区域, 所述第一阈值使得非线性参量分布图中小于第一阈值的区域中绝大多 数在该造影信息帧中对应的点属于组织残留信息; 形成第一图像, 对该造影信息帧中属于组织残留区域的信息进行抑制 处理, 对该造影信息帧中属于造影剂区域的信息进行增强处理, 并基 于处理后的结果形成第一图像数据;
    第二阈值分割, 采用第二阈值对非线性参量分布图进行分割, 小于第 二阈值的区域为组织残留区域, 大于第二阈值的区域为造影剂区域, 所述第二阈值使得非线性参量分布图中大于第二阈值的区域中绝大多 数在该造影信息帧中对应的点属于造影剂信息, 所述第二阈值大于第 一阈值;
    形成第二图像, 对该造影信息帧中属于组织残留区域的信息进行抑制 处理, 对该造影信息帧中属于造影剂区域的信息进行增强处理, 并基 于处理后的结果形成第二图像数据;
    图像融合, 将第一图像数据和第二图像数据进行融合, 形成增强后的 造影图像数据。
  3. [权利要求 2] 如权利要求 1所述的方法, 其特征在于, 所述增强处理包括乘以大于 1 的权重系数或加上大于 0的权重系数, 所述抑制处理包括乘以小于 1的 权重系数或加上小于 0的权重系数。
  4. [权利要求 3] 如权利要求 2所述的方法, 其特征在于, 形成第一图像数据过程中进 行抑制处理的权重系数小于形成第二图像数据过程中进行抑制处理的 权重系数, 形成第二图像数据过程中进行增强处理的权重系数大于形 成第一图像数据过程中进行增强处理的权重系数。
  5. [权利要求 4] 如权利要求 1所述的方法, 其特征在于, 在计算非线性参量之前还包 括: 判断造影信息中的噪声, 将一帧造影信息分为噪声信号和除去噪 声的非噪声造影信号; 在计算非线性参量吋采用非噪声造影信号计算 非线性参量。
  6. [权利要求 5] 如权利要求 1所述的方法, 其特征在于, 所述非线性参量分布图根据 同一吋刻的一帧非线性造影信息和一帧线性组织信息的比值或差值计 算得出。
  7. [权利要求 6] 如权利要求 1至 5中任一项所述的方法, 其特征在于, 所述图像融合包 括: 将第一图像数据和第二图像数据分别乘以第一权重系数和第二权 重系数, 所述第一权重系数与第二权重系数的和为 1且第二权重系数 大于或等于第一权重系数。
  8. [权利要求 7] —种超声造影图像增强系统,其特征在于包括:
    信息获取单元, 用于根据造影成像过程中的超声回波信号获得非线性 的造影信息和线性的组织信息;
    非线性参量计算单元, 用于比较所述造影信息和组织信息得到一帧造 影信息的非线性参量分布图, 所述非线性参量是衡量超声波在媒质中 传播吋产生的非线性效应大小的参量;
    第一阈值分割单元, 用于采用第一阈值对非线性参量分布图进行分割 , 小于第一阈值的区域为组织残留区域, 大于第一阈值的区域为造影 剂区域, 所述第一阈值使得非线性参量分布图中小于第一阈值的区域 中绝大多数在该造影信息帧中对应的点属于组织残留信息; 第一图像形成单元, 用于对该造影信息帧中属于组织残留区域的信息 进行抑制处理, 对该造影信息帧中属于造影剂区域的信息进行增强处 理, 并基于处理后的结果形成第一图像数据;
    第二阈值分割单元, 用于采用第二阈值对非线性参量分布图进行分割 , 小于第二阈值的区域为组织残留区域, 大于第二阈值的区域为造影 剂区域, 所述第二阈值使得非线性参量分布图中大于第二阈值的区域 中绝大多数在该造影信息帧中对应的点属于造影剂信息, 所述第二阈 值大于第一阈值;
    第二图像形成单元, 用于对该造影信息帧中属于组织残留区域的信息 进行抑制处理, 对该造影信息帧中属于造影剂区域的信息进行增强处 理, 并基于处理后的结果形成第二图像数据;
    图像融合单元, 用于将第一图像数据和第二图像数据进行融合, 形成 增强后的造影图像数据。
  9. [权利要求 8] 如权利要求 7所述的系统, 其特征在于, 所述增强处理包括乘以大于 1 的权重系数或加上大于 0的权重系数, 所述抑制处理包括乘以小于 1的 权重系数或加上小于 0的权重系数。
  10. [权利要求 9] 如权利要求 8所述的系统, 其特征在于, 第一图像形成单元进行抑制 处理的权重系数小于第二图像形成单元进行抑制处理的权重系数, 第 二图像形成单元进行增强处理的权重系数大于第一图像形成单元进行 增强处理的权重系数。
  11. [权利要求 10] 如权利要求 7所述的系统, 其特征在于, 还包括噪声识别单元, 所述 噪声识别单元用于在非线性参量计算单元计算非线性参量之前判断造 影信息中的噪声, 将一帧造影信息分为噪声造影和除去噪声的非噪声 造影信号; 非线性参量计算单元在计算非线性参量吋采用非噪声信号 计算非线性参量。
  12. [权利要求 11] 如权利要求 7所述的系统, 其特征在于, 所述非线性参量计算单元根 据同一吋刻的一帧非线性造影信息和一帧线性组织信息的比值或差值 计算该造影信息帧的非线性参量分布图。
  13. [权利要求 12] 如权利要求 7至 11中任一项所述的系统, 其特征在于, 所述图像融合 单元将第一图像数据和第二图像数据分别乘以第一权重系数和第二权 重系数后进行融合, 所述第一权重系数与第二权重系数的和为 1且第 二权重系数大于或等于第一权重系数。 [权利要求 13] 如权利要求 7至 11中任一项所述的系统, 其特征在于, 所述非线性参 量分布图为直方图或二维图。
  14. [权利要求 12] 如权利要求 7至 11中任一项所述的系统, 其特征在于, 所述图像融合 单元将第一图像数据和第二图像数据分别乘以第一权重系数和第二权 重系数后进行融合, 所述第一权重系数与第二权重系数的和为 1且第 二权重系数大于或等于第一权重系数。 [权利要求 13] 如权利要求 7至 11中任一项所述的系统, 其特征在于, 所述非线性参 量分布图为直方图或二维图。
  15. [权利要求 14] 一种超声造影成像设备, 其特征在于包括:
    探头, 用于向组织发射超声波并接收组织反射回的超声回波; 前信号处理模块, 对超声回波信号进行处理, 从超声回波信号中提取 反映造影剂微泡信息的非线性成分和反映组织解剖特征的线性成分; 第一信号处理模块, 用于对线性成分进行处理; 组织图像生成模块, 用于将第一信号处理模块的输出处理成组织图像 数据;
    第二信号处理模块, 用于对非线性成分进行处理; 造影图像生成模块, 用于将第二信号处理模块的输出处理成造影图像 数据;
    如权利要求 7至 13中任一项所述的超声造影图像增强系统, 用于对造 影信息进行增强, 以得到增强后的造影图像数据; 显示模块, 用于对组织图像数据和增强后的造影图像数据进行可视化 显示。
  16. [权利要求 15] 如权利要求 14所述的超声造影成像设备, 其特征在于, 所述超声造影 图像增强系统对前信号处理模块输出的射频域的造影信息进行增强处 理。
  17. [权利要求 16] 如权利要求 14所述的超声造影成像设备, 其特征在于, 所述超声造影 图像增强系统对组织图像生成模块和造影图像生成模块输出的图像域 的造影信息进行增强处理。
  18. [权利要求 17] 如权利要求 14所述的超声造影成像设备, 其特征在于, 所述超声造影 图像增强系统对第二信号处理模块中任一环节的造影信息进行增强处 理。
CN201680039006.7A 2016-06-30 2016-06-30 一种增强超声造影图像的方法、系统及超声造影成像设备 Active CN108135566B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2016/087979 WO2018000359A1 (zh) 2016-06-30 2016-06-30 一种增强超声造影图像的方法、系统及超声造影成像设备

Publications (2)

Publication Number Publication Date
CN108135566A true CN108135566A (zh) 2018-06-08
CN108135566B CN108135566B (zh) 2020-11-24

Family

ID=60785845

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680039006.7A Active CN108135566B (zh) 2016-06-30 2016-06-30 一种增强超声造影图像的方法、系统及超声造影成像设备

Country Status (2)

Country Link
CN (1) CN108135566B (zh)
WO (1) WO2018000359A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110276717A (zh) * 2019-06-26 2019-09-24 纳米视觉(成都)科技有限公司 一种图像的拼接方法及终端
CN111200699A (zh) * 2018-11-19 2020-05-26 瑞昱半导体股份有限公司 影像调整方法
CN113164160A (zh) * 2020-05-18 2021-07-23 深圳迈瑞生物医疗电子股份有限公司 超声造影成像方法、装置和存储介质
CN113538299A (zh) * 2021-09-13 2021-10-22 深圳瀚维智能医疗科技有限公司 超声图像去噪方法、装置、设备及计算机可读存储介质
CN113768547A (zh) * 2021-09-14 2021-12-10 南京超维景生物科技有限公司 冠状动脉成像方法及装置、存储介质及电子设备

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112837254A (zh) * 2021-02-25 2021-05-25 普联技术有限公司 一种图像融合方法、装置、终端设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010039381A1 (en) * 2000-02-16 2001-11-08 Burns Peter N. Ultrasonic harmonic image segmentation
US6371914B1 (en) * 2000-04-13 2002-04-16 Bracco Research S.A. Single-shot phase cancellation ultrasound contrast imaging
US20030097069A1 (en) * 2001-11-21 2003-05-22 Ge Medical Systems Global Technology Company, Llc Computationally efficient noise reduction filter for enhancement of ultrasound images
CN104720850A (zh) * 2013-12-23 2015-06-24 深圳迈瑞生物医疗电子股份有限公司 一种超声造影成像方法及造影图像的区域检测、显像方法
CN105631867A (zh) * 2015-12-25 2016-06-01 中国科学院深圳先进技术研究院 一种全自动超声造影图像分割方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4662380A (en) * 1985-10-31 1987-05-05 General Electric Company Adaptive time gain compensation system for ultrasound imaging
US6679844B2 (en) * 2002-06-20 2004-01-20 Acuson Corporation Automatic gain compensation for multiple mode or contrast agent imaging
CN101987023B (zh) * 2009-07-31 2014-05-07 深圳迈瑞生物医疗电子股份有限公司 一种超声成像增益补偿及图像优化方法及其装置和系统
CN103845077B (zh) * 2012-12-05 2016-01-20 深圳迈瑞生物医疗电子股份有限公司 超声图像增益优化方法及超声成像增益自动优化装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010039381A1 (en) * 2000-02-16 2001-11-08 Burns Peter N. Ultrasonic harmonic image segmentation
US6371914B1 (en) * 2000-04-13 2002-04-16 Bracco Research S.A. Single-shot phase cancellation ultrasound contrast imaging
US20030097069A1 (en) * 2001-11-21 2003-05-22 Ge Medical Systems Global Technology Company, Llc Computationally efficient noise reduction filter for enhancement of ultrasound images
CN104720850A (zh) * 2013-12-23 2015-06-24 深圳迈瑞生物医疗电子股份有限公司 一种超声造影成像方法及造影图像的区域检测、显像方法
CN105631867A (zh) * 2015-12-25 2016-06-01 中国科学院深圳先进技术研究院 一种全自动超声造影图像分割方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111200699A (zh) * 2018-11-19 2020-05-26 瑞昱半导体股份有限公司 影像调整方法
CN111200699B (zh) * 2018-11-19 2022-04-26 瑞昱半导体股份有限公司 影像调整方法
CN110276717A (zh) * 2019-06-26 2019-09-24 纳米视觉(成都)科技有限公司 一种图像的拼接方法及终端
CN110276717B (zh) * 2019-06-26 2023-05-05 图码思(成都)科技有限公司 一种图像的拼接方法及终端
CN113164160A (zh) * 2020-05-18 2021-07-23 深圳迈瑞生物医疗电子股份有限公司 超声造影成像方法、装置和存储介质
WO2021232192A1 (zh) * 2020-05-18 2021-11-25 深圳迈瑞生物医疗电子股份有限公司 超声造影成像方法、装置和存储介质
CN113164160B (zh) * 2020-05-18 2022-09-16 深圳迈瑞生物医疗电子股份有限公司 超声造影成像方法、装置和存储介质
CN113538299A (zh) * 2021-09-13 2021-10-22 深圳瀚维智能医疗科技有限公司 超声图像去噪方法、装置、设备及计算机可读存储介质
CN113538299B (zh) * 2021-09-13 2022-01-11 深圳瀚维智能医疗科技有限公司 超声图像去噪方法、装置、设备及计算机可读存储介质
CN113768547A (zh) * 2021-09-14 2021-12-10 南京超维景生物科技有限公司 冠状动脉成像方法及装置、存储介质及电子设备
CN113768547B (zh) * 2021-09-14 2024-03-22 南京超维景生物科技有限公司 冠状动脉成像方法及装置、存储介质及电子设备

Also Published As

Publication number Publication date
WO2018000359A1 (zh) 2018-01-04
CN108135566B (zh) 2020-11-24

Similar Documents

Publication Publication Date Title
CN108135566A (zh) 一种增强超声造影图像的方法、系统及超声造影成像设备
Karaman et al. An adaptive speckle suppression filter for medical ultrasonic imaging
Abd-Elmoniem et al. Real-time speckle reduction and coherence enhancement in ultrasound imaging via nonlinear anisotropic diffusion
US8202220B2 (en) Ultrasound diagnostic apparatus
CN108836392B (zh) 基于超声rf信号的超声成像方法、装置、设备及存储介质
WO2012051216A1 (en) Direct echo particle image velocimetry flow vector mapping on ultrasound dicom images
JP7041147B2 (ja) 造影剤流の肝灌流を特徴付けるシステム及び方法
US20090024033A1 (en) Ultrasound diagnostic apparatus
Hu et al. Acoustic shadow detection: study and statistics of B-mode and radiofrequency data
Tasnim et al. Study of speckle noise reduction from ultrasound B-mode images using different filtering techniques
US8047991B2 (en) Automatic identification of orientation in medical diagnostic ultrasound
CN104720850B (zh) 一种超声造影成像方法及造影图像的区域检测、显像方法
Abd-Elmoniem et al. Real time adaptive ultrasound speckle reduction and coherence enhancement
CN107169978B (zh) 超声图像边缘检测方法及系统
JPH08511184A (ja) 超音波走査のための方法及び装置
WO2016161574A1 (zh) 一种超声造影成像方法及装置
JPH0611568A (ja) 画像処理システム及び超音波走査イメージングシステム
Shen et al. Deep learning-based wavelet threshold function optimization on noise reduction in ultrasound images
Tsui et al. An adaptive threshold filter for ultrasound signal rejection
CN111476790A (zh) 一种超声穿刺中穿刺针增强显示方法及装置
Jayanthi Sree et al. De-speckling of ultrasound images using local statistics-based trilateral filter
Hourani et al. Block-wise ultrasound image deconvolution from fundamental and harmonic images
Wang et al. A novel edge enhancement method for ultrasound imaging
CN112263278B (zh) 基于随机孔径的心脏超声杂波抑制方法及系统
Gupta et al. Linear despeckle approach for ultrasound carotid artery images

Legal Events

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